Replication study of multiple sclerosis (MS) susceptibility alleles and correlation of DNA-variants with disease features in a cohort of Austrian MS patients

Association of CD58 rs12044852 and rs2300747 polymorphisms with the risk of multiple sclerosis: A systematic review and meta-analysis

OBJECTIVE

Multiple sclerosis is a serious neurodegenerative disorder that causes disability in young adults. Genetic predisposition of multiple sclerosis is well documented and several single nucleotide polymorphisms (SNPs) of the CD58 were found to be associated with this disease. This systematic review and meta-analysis were done with the aim of finding the association between CD58 gene SNPs (rs12044852 and rs2300747) and the risk of multiple sclerosis (MS).

METHOD

A comprehensive search was done in PubMed, Google Scholar, Embase, and MSGene.org to find the relevant data. Our search yielded 13 relevant publications which were included for meta-analysis consisting of 5194 cases and 5766 controls. All the statistical analysis was conducted using meta and metafor packages in R studio. The odds ratio (OR) along with 95 % confidence intervals and p values were determined using the fixed effects and random effects model. The I2 test was done to measure heterogeneity. Subgroup analysis was performed along with analysis for publication bias.

RESULTS

We found significant association for both rs12044852 (allelic, dominant, over-dominant, heterozygous, and homozygous models) and rs2300747 (allelic, dominant, over-dominant, heterozygous models) with multiple sclerosis. Both the SNPs provided a protective effect for multiple sclerosis. Subgroup analysis indicated that rs12044852 polymorphism provided a protective effect in both Asians and Caucasians. However, for rs2300747, the Asian population showed no statistically significant association with the risk of MS.

CONCLUSION

Polymorphism of rs12044852 and rs2300747 of the CD58 gene provided a protective effect for multiple sclerosis. The protective effect is more prominent in Caucasian populations compared to Asians.

mRNA expression profile reveals differentially expressed genes in splenocytes of experimental autoimmune encephalomyelitis model

Experimental autoimmune encephalomyelitis (EAE) is a mouse model that can be used to investigate aetiology, pathogenesis, and treatment approaches for multiple sclerosis (MS). A novel integrated bioinformatics approach was used to understand the involvement of differentially expressed genes (DEGs) in the spleen of EAE mice through data mining of existing microarray and RNA-seq datasets. We screened differentially expressed mRNAs using mRNA expression profile data of EAE spleens taken from Gene Expression Omnibus (GEO). Functional and pathway enrichment analyses of DEGs were performed by Database for Annotation, Visualization, and Integrated Discovery (DAVID). Subsequently, the DEGs-encoded protein-protein interaction (PPI) network was constructed. The 784 DEGs in GSE99300 A.SW PP-EAE mice spleen mRNA profiles, 859 DEGs in GSE151701 EAE mice spleen mRNA profiles, and 646 DEGs in GSE99300 SJL/J PP-EAE mice spleen mRNA profiles were explored. Functional enrichment of 55 common DEGs among 3 sub-datasets revealed several immune-related terms, such as neutrophil extravasation, leucocyte migration, antimicrobial humoral immune response mediated by an antimicrobial peptide, toll-like receptor 4 bindings, IL-17 signalling pathway, and TGF-beta signalling pathway. In the screening of 10 hub genes, including MPO, ELANE, CTSG, LTF, LCN2, SELP, CAMP, S100A9, ITGA2B, and PRTN3, and in choosing and validating the 5 DEGs, including ANK1, MBOAT2, SLC25A21, SLC43A1, and SOX6, the results showed that SLC43A1 and SOX6 were significantly decreased in EAE mice spleen. Thus this study offers a list of genes expressed in the spleen that might play a key role in the pathogenesis of EAE.

Circadian Genes MBOAT2/CDA/LPCAT2/B4GALT5 in the Metabolic Pathway Serve as New Biomarkers of PACA Prognosis and Immune Infiltration

Pancreatic cancer (PACA) is a highly malignant tumor with a poor prognosis. Recent studies have discovered substantial differences in the expression levels of several circadian genes in PACA samples compared to normal samples. The goal of this research was to find differentially expressed rhythm genes (DERGs) in PACA samples and determine their role in the development of PACA. A total of 299 DERGs were identified in PACA, including 134 downregulated genes and 165 upregulated genes. DERGs were significantly abundant in the metabolic pathway and immune response pathways, according to GO and KEGG analyses. Survival analyses showed that PACA patients who had higher expression levels of MBOAT2/CDA/LPCAT2/B4GALT5 had shorter overall survival times. Using cell assay verification, the mRNA levels of MBOAT2/CDA/LPCAT2/B4GALT5 in Patu-8988 and PNAC-1 cells were found to be significantly higher than those in HPDE6-C7 cells, which was in line with previous studies on PACA patient data. Through conducting univariate Cox analysis, it was determined that MBOAT2/CDA/LPCAT2/B4GALT5 expression, age and grade were all high-risk factors. The MBOAT2/CDA/LPCAT2/B4GALT5 genes were independently correlated with overall survival, according to the multivariate Cox analysis. The proportion of immune cells in PACA and normal samples significantly changed, according to the immune infiltration analysis. Furthermore, MBOAT2/CDA/LPCAT2/B4GALT5 expression levels were significantly related to the level of immune cell infiltration. The protein-protein interaction network of the MBOAT2/CDA/LPCAT2/B4GALT5 genes included 54 biological nodes and 368 interacting genes. In conclusion, the finding of these DERGs adds to the investigation of the molecular processes underlying the onset and progression of PACA. In the future, DERGs may serve as prognostic and diagnostic biomarkers as well as drug targets for chronotherapy in PACA patients.

Analysis of herpesvirus infection and genome single nucleotide polymorphism risk factors in multiple sclerosis, Volga federal district, Russia

Multiple sclerosis (MS) is a heterogeneous disease where herpesvirus infection and genetic predisposition are identified as the most consistent risk factors. Serum and blood samples were collected from 151 MS and 70 controls and used to analyze circulating antibodies for, and DNA of, Epstein Barr virus (EBV), human cytomegalovirus (HCMV), human herpes virus 6 (HHV6), and varicella zoster virus (VZV). The frequency of selected single nucleotide polymorphisms (SNPs) in MS and controls were studied. Herpesvirus DNA in blood samples were analyzed using qPCR. Anti-herpesvirus antibodies were detected by ELISA. SNPs were analyzed by the allele-specific PCR. For statistical analysis, Fisher exact test, odds ratio and Kruskall-Wallis test were used; p<0.05 values were considered as significant. We have found an association between circulating anti-HHV6 antibodies and MS diagnosis. We also confirmed higher frequency of A and C alleles in rs2300747 and rs12044852 of CD58 gene and G allele in rs929230 of CD6 gene in MS as compared to controls. Fatigue symptom was linked to AC and AA genotype in rs12044852 of CD58 gene. An interesting observation was finding higher frequency of GG genotype in rs12722489 of IL2RA and T allele in rs1535045 of CD40 genes in patient having anti-HHV6 antibodies. A link was found between having anti-VZV antibodies in MS and CC genotype in rs1883832 of CD40 gene.

Identification of MBOAT2 as an Unfavorable Biomarker Correlated with KRAS Activation and Reduced CD8+ T-Cell Infiltration in Pancreatic Cancer

Objectives

Limited research on the role of membrane-bound O-acyltransferase domain–containing 2 (MBOAT2) in cancer biology exists. In particular, the underlying role of MBOAT2 and its potential mechanisms in pancreatic cancer have not yet been explored. Further study of MBOAT2 could provide new ideas about the carcinogenesis and treatment of pancreatic cancer (PC).

Methods

In the current study, the potential biological and clinical significances of MBOAT2 were explored by bioinformatics analysis. Real-time quantitative polymerase chain reaction and western blot analysis were performed to determine the level of MBOAT2 in pancreatic ductal adenocarcinoma (PDAC) cell lines. MTT, colony formation, and Transwell assays and flow cytometry of cell cycle were performed to analyze PDAC cell proliferation, migration, and cycle progression. The potential relationship between MBOAT2 level and tumor immunity was analyzed using the ESTIMATE algorithm, CIBERSORT algorithm, and single-sample gene set enrichment analysis.

Results

The level of MBOAT2 was remarkably upregulated in most tumors, especially pancreatic tumors, and was positively correlated with a greater rate of tumor recurrence, higher histologic grade, and worse overall survival. MBOAT2 overexpression was also closely correlated with the mutation status and expression level of driver genes, especially KRAS. Meanwhile, functional enrichment analysis demonstrated that MBOAT2 might be involved in cell–cell communication; cell cycling; the Ras signaling pathway; and immune-related biological functions such as the leukocyte activation involved in T-cell–receptor signaling pathway, the inflammatory response, and antigen processing and presentation. Furthermore, in vitro experiments demonstrated that MBOAT2 overexpression accelerated PC cell proliferation and migration. MBOAT2 overexpression also enhanced CDK2 and CCNA2 expression, leading to cell cycle progression from the G1 phase to the G2 phase. Lastly, MBOAT2 overexpression reduced the infiltration level of CD8⁺ T-cells, plasmacytoid dendritic cells, and activated dendritic cells but triggered a high type-2 T helper/type-1 T helper cell ration (Th2/Th1 ration) in PC.

Conclusion

Our findings suggest that MBOAT2 is a potential protooncogene in PDAC that predicts a poor prognosis and is related to KRAS activation and inferior infiltration of CD8⁺ T-cells in PC.

Novel Insights Into Gene Signatures and Their Correlation With Immune Infiltration of Peripheral Blood Mononuclear Cells in Behcet's Disease

Background

Behcet’s disease (BD) is a chronic inflammatory disease that involves systemic vasculitis and mainly manifests as oral and genital ulcers, uveitis, and skin damage as the first clinical symptoms, leading to gastrointestinal, aortic, or even neural deterioration. There is an urgent need for effective gene signatures for BD’s early diagnosis and elucidation of its underlying etiology.

Methods

We identified 82 differentially expressed genes (DEGs) in BD cases compared with healthy controls (HC) after combining two Gene Expression Omnibus datasets. We performed pathway analyses on these DEGs and constructed a gene co-expression network and its correlation with clinical traits. Hub genes were identified using a protein–protein interaction network. We manually selected CCL4 as a central hub gene, and gene-set enrichment and immune cell subset analyses were applied on patients in high- and low-CCL4 expression groups. Meanwhile, we validated the diagnostic value of hub genes in differentiating BD patients from HC in peripheral blood mononuclear cells using real-time PCR.

Results

Twelve hub genes were identified, and we validated the upregulation of CCL4 and the downregulation of NPY2R mRNA expression. Higher expression of CCL4 was accompanied by larger fractions of CD8 + T cells, natural killer cells, M1 macrophages, and activated mast cells. Receiver operator characteristic curves showed good discrimination between cases and controls based on the expression of these genes.

Conclusion

CCL4 and NPY2R could be diagnostic biomarkers for BD that reveal inflammatory status and predict vascular involvement in BD, respectively.

The S100A7 nuclear interactors in autoimmune diseases: a coevolutionary study in mammals

S100A7, a member of the S100A family of Ca²⁺-binding proteins, is considered a key effector in immune response. In particular, S100A7 dysregulation has been associated with several diseases, including autoimmune disorders. At the nuclear level, S100A7 interacts with several protein-binding partners which are involved in transcriptional regulation and DNA repair. By using the BioGRID and GAAD databases, S100A7 nuclear interactors with a putative involvement in autoimmune diseases were retrieved. We selected fatty acid-binding protein 5 (FABP5), autoimmune regulator (AIRE), cystic fibrosis transmembrane conductance regulator (CFTR), chromodomain helicase DNA-binding protein 4 (CHD4), epidermal growth factor receptor (EGFR), estrogen receptor 1 (ESR1), histone deacetylase 2 (HDAC2), v-myc avian myelocytomatosis viral oncogene homolog (MYC), protection of telomeres protein 1 (POT1), telomeric repeat-binding factor (NIMA-interacting) 1 (TERF1), telomeric repeat-binding factor 2 (TERF2), and Zic family member 1 (ZIC1). Linear correlation coefficients between interprotein distances were calculated with MirrorTree. Coevolution clusters were also identified with the use of a recent version of the Blocks in Sequences (BIS2) algorithm implemented in the BIS2Analyzer web server. Analysis of pair positions identified interprotein coevolving clusters between S100A7 and the binding partners CFTR and TERF1. Such findings could guide further analysis to better elucidate the function of S100A7 and its binding partners and to design drugs targeting for these molecules in autoimmune diseases.

Predicting signaling pathways regulating demyelination in a rat model of lithium-pilocarpine-induced acute epilepsy: A proteomics study

Demyelination is observed in animal models of intractable epilepsy (IE). Epileptogenesis damages the myelin sheath and dysregulates oligodendrocyte precursor cell (OPC) development. However, the molecular pathways regulating demyelination in epilepsy are unclear. Here, we predicted the molecular mechanisms regulating demyelination in a rat model of lithium-pilocarpine hydrochloride-induced epilepsy. We identified DGKA/Mboat2/Inpp5j and NOS/Keratin 28 as the main target molecules that regulate demyelination via glycerolipid and glycerophospholipid metabolism, phosphatidylinositol signaling, and estrogen signaling in demyelinated forebrain slice cultures (FSCs). In seizure-like FCSs, the actin cytoskeleton was regulated by Cnp and MBP via Pak4/Tmsb4x (also known as Tβ4) and Kif5c/Kntc1. Tβ4 possibly prevented OPC differentiation and maturation and inhibited MBP phosphorylation via the p38MAPK/ERK1/JNK1 pathway. The MAPK signaling pathway was more likely activated in seizure-like FCSs than in demyelinated FCSs. pMBP expression was decreased in the hippocampus of lithium-pilocarpine hydrochloride-induced acute epilepsy rats. The expression of remyelination-related factors was suppressed in the hippocampus and corpus callosum in lithium-pilocarpine hydrochloride-induced epilepsy rats. These findings suggest that the actin cytoskeleton, Tβ4, and MAPK signaling pathways regulate the decrease in pMBP in the hippocampus in a rat model of epilepsy. Our results indicate that regulating the actin cytoskeleton, Tβ4, and MAPK signaling pathways may facilitate the prevention of demyelination in IE.

Adaptive Fisher method detects dense and sparse signals in association analysis of SNV sets

BACKGROUND

With the development of next generation sequencing (NGS) technology and genotype imputation methods, statistical methods have been proposed to test a set of genomic variants together to detect if any of them is associated with the phenotype or disease. In practice, within the set, there is an unknown proportion of variants truly causal or associated with the disease. There is a demand for statistical methods with high power in both dense and sparse scenarios, where the proportion of causal or associated variants is large or small respectively.

RESULTS

We propose a new association test - weighted Adaptive Fisher (wAF) that can adapt to both dense and sparse scenarios by adding weights to the Adaptive Fisher (AF) method we developed before. Using simulation, we show that wAF enjoys comparable or better power to popular methods such as sequence kernel association tests (SKAT and SKAT-O) and adaptive SPU (aSPU) test. We apply wAF to a publicly available schizophrenia dataset, and successfully detect thirteen genes. Among them, three genes are supported by existing literature; six are plausible as they either relate to other neurological diseases or have relevant biological functions.

CONCLUSIONS

The proposed wAF method is a powerful disease-variants association test in both dense and sparse scenarios. Both simulation studies and real data analysis indicate the potential of wAF for new biological findings.

[Association analysis of cytokine receptors' genes polymorphisms with clinical features of multiple sclerosis]

AIM

To study the association of polymorphisms in the IL2RA and TNFRSF1A genes with severity and early clinical manifestations of remitted multiple sclerosis (MS).

MATERIAL AND METHODS

Five hundred and eight patients of Russian ethnicity with bout-onset MS were genotyped for IL7RA (rs6897932), IL2RA (rs2104286) and TNFRSF1A (rs1800693) polymorphisms. Association analysis of the gene variants with disease severity, variants of MS manifestation, and first remission duration was performed.

RESULTS AND CONCLUSION

Dividing the MS patients by disease severity, estimated with the MSSS, we found a significant increase in the TNFRSF1A*T/T genotype carriage in patients with milder MS course (MSSS≤3), and, respectively, in the TNFRSF1A*C allele carriage in patients with moderate to severe MS (MSSS> 3). Dividing the MS patients into two groups according to their MS manifestation variants, we revealed a significant increase in the TNFRSF1A*T allele carriage in patients with favorable variants of MS manifestation (optic neuritis or sensory disturbances), and of the TNFRSF1A*C/C genotype in patients with unfavorable variants (motor disorders, brain stem disorders, impaired coordination, pelvic disorders, mental disorders or polysymptomatic onset). No associations with first remission duration were observed. Multi-locus analysis to search for allelic combinations associated with the studied clinical features of MS was applied. In this analysis, a polymorphic variant of CTLA4 gene (rs231775), for which we have previously reported the association of the CTLA4*G allele with short first remission (less than 1 year), was also included. The carriage of biallelic combination (CTLA4*G + TNFRSF1A*C) was associated with short first remission more significantly than the carriage of CTLA4*G by itself. One more biallelic combination associated with short first remission (CTLA4*G/G + IL7RA*T), was identified. No other biallelic combinations significantly associated with the clinical features studied were observed.

This Is the End: Regulation of Rab7 Nucleotide Binding in Endolysosomal Trafficking and Autophagy

Rab7 – or in yeast, Ypt7p – governs membrane trafficking in the late endocytic and autophagic pathways. Rab7 also regulates mitochondrion-lysosome contacts, the sites of mitochondrial fission. Like all Rab GTPases, Rab7 cycles between an “active” GTP-bound form that binds downstream effectors – e.g., the HOPS and retromer complexes and the dynactin-binding Rab-interacting lysosomal protein (RILP) – and an “inactive” GDP-bound form that cannot bind effectors. Accessory proteins regulate the nucleotide binding state of Rab7: guanine nucleotide exchange factors (GEFs) stimulate exchange of bound GDP for GTP, resulting in Rab7 activation, whereas GTPase activating proteins (GAPs) boost Rab7’s GTP hydrolysis activity, thereby inactivating Rab7. This review will discuss the GEF and GAPs that control Rab7 nucleotide binding, and thus regulate Rab7’s activity in endolysosomal trafficking and autophagy. It will also consider how bacterial pathogens manipulate Rab7 nucleotide binding to support intracellular invasion and immune evasion.

The Impact of Five VDR Polymorphisms on Multiple Sclerosis Risk and Progression: a Case-Control and Genotype-Phenotype Study

Vitamin D receptor polymorphisms have been the target of many studies focusing on multiple sclerosis. However, previously reported results have been inconclusive. The objective of this study was to investigate the association between five vitamin D receptor polymorphisms (EcoRV, FokI, ApaI, TaqI, and BsmI) and multiple sclerosis susceptibility and its course. The study was carried out as a case-control and genotype-phenotype study, consisted of 296 Czech multiple sclerosis patients and 135 healthy controls. Genotyping was carried out using polymerase chain reaction and restriction analysis. In multiple sclerosis men, allele and/or genotype distributions differed in EcoRV, TaqI, BsmI, and ApaI polymorphisms as compared to controls (EcoRV, p_a = 0.02; Taq, p_g = 0.02, p_a = 0.02; BsmI, p_g = 0.02, p_a = 0.04; ApaI, p_g = 0.008, p_a = 0.005). In multiple sclerosis women, differences in the frequency of alleles and genotypes were found to be significant in ApaI (controls vs multiple sclerosis women: p_g = 0.01, p_a = 0.05). Conclusive results were observed between multiple sclerosis women in the case of EcoRV [differences in Expanded Disability Status Scale (p = 0.05); CT genotype was found to increase the risk of primary progressive multiple sclerosis 5.5 times (CT vs CC+TT p_corr = 0.01, sensitivity 0.833, specificity 0.525, power test 0.823)] and FokI [borderline difference in Multiple Sclerosis Severity Score (p = 0.05)]. Our results indicate that the distribution of investigated vitamin D receptor polymorphisms is a risk factor for multiple sclerosis susceptibility and progression in the Czech population. The association between disease risk and polymorphisms was found to be stronger in men. The association of disease progression with polymorphisms was observed only in women.

[Clinical and molecular genetic analysis of a case of familial multiple sclerosis in the Republic of Bashkortostan]

AIM

To investigate clinical manifestations of multiple sclerosis (MS) and the genetic makeup of six affected members of one family.

MATERIAL AND METHODS

Six members of the family of Russian ethnic origin were examined. Pedigree analysis and genotyping of polymorphic markers of candidate genes for multiple sclerosis were performed.

RESULTS AND CONCLUSION

The accumulation of alleles that were associated with autoimmune diseases according to the results of genome-wide association studies (rs1109670*C, rs3129934*T, rs9523762*G, rs1570538*T) was found in the family. The results confirm the contribution of several genetic variants to familial forms of MS.

A candidate gene approach of the calcineurin pathway to identify variants associated with clinical outcomes in renal transplantation

AIM

To investigate the potential influence of variants in genes involved in the calcineurin pathway on the efficacy and toxicity of calcineurin inhibitors in renal transplantation.

MATERIALS & METHODS

Twenty-three polymorphisms in thirteen genes were tested in 381 renal transplant recipients receiving ciclosporin (n = 221) or tacrolimus (n = 160) and mycophenolate mofetil. Data were collected prospectively over the first year post-transplantation.

RESULTS

Multivariate survival analyses revealed no genetic associations with biopsy proven acute graft rejection and serious infections. Donor-recipient Cytomegalovirus mismatch was the only variable associated with serious infection.

CONCLUSION

This large exploratory study casts doubts on the potential interest of genetic biomarkers related to CNI pharmacodynamics but associations with other phenotypes in transplantation deserve further studies.

Genome-wide association study of posttraumatic stress disorder in a cohort of Iraq-Afghanistan era veterans

BACKGROUND

Posttraumatic stress disorder (PTSD) is a psychiatric disorder that can develop after experiencing traumatic events. A genome-wide association study (GWAS) design was used to identify genetic risk factors for PTSD within a multi-racial sample primarily composed of U.S. veterans.

METHODS

Participants were recruited at multiple medical centers, and structured interviews were used to establish diagnoses. Genotypes were generated using three Illumina platforms and imputed with global reference data to create a common set of SNPs. SNPs that increased risk for PTSD were identified with logistic regression, while controlling for gender, trauma severity, and population substructure. Analyses were run separately in non-Hispanic black (NHB; n = 949) and non-Hispanic white (NHW; n = 759) participants. Meta-analysis was used to combine results from the two subsets.

RESULTS

SNPs within several interesting candidate genes were nominally significant. Within the NHB subset, the most significant genes were UNC13C and DSCAM. Within the NHW subset, the most significant genes were TBC1D2, SDC2 and PCDH7. In addition, PRKG1 and DDX60L were identified through meta-analysis. The top genes for the three analyses have been previously implicated in neurologic processes consistent with a role in PTSD. Pathway analysis of the top genes identified alternative splicing as the top GO term in all three analyses (FDR q < 3.5 × 10(-5)).

LIMITATIONS

No individual SNPs met genome-wide significance in the analyses.

CONCLUSIONS

This multi-racial PTSD GWAS identified biologically plausible candidate genes and suggests that post-transcriptional regulation may be important to the pathology of PTSD; however, replication of these findings is needed.

Variants of MicroRNA Genes: Gender-Specific Associations with Multiple Sclerosis Risk and Severity

Multiple sclerosis (MS) is an autoimmune neuro-inflammatory disease arising from complex interactions of genetic, epigenetic, and environmental factors. Variations in genes of some microRNAs--key post-transcriptional regulators of many genes--can influence microRNAs expression/function and contribute to MS via expression changes of protein-coding target mRNA genes. We performed an association study of polymorphous variants of MIR146A rs2910164, MIR196A2 rs11614913, MIR499A rs3746444 MIR223 rs1044165 and their combinations with MS risk and severity. 561 unrelated patients with bout-onset MS and 441 healthy volunteers were enrolled in the study. We observed associations of MS risk with allele MIR223*T and combination (MIR223*T + MIR146A*G/G) carriage in the entire groups and in women at Bonferroni-corrected significance level (pcorr < 0.05). Besides, MIR146A*G/G association with MS was observed in women with nominal significance (pf = 0.025). No MS associations were found in men. A more severe MS course (MSSS value > 3.5) was associated with the carriage of MIR499A*C/T and, less reliably, of MIR499A*C (pcorr = 0.006 and pcorr = 0.024, respectively) and with the carriage of combinations (MIR499A*C/T + MIR196A2*C) and (MIR499A*C + MIR196A2*C) (pcorr = 0.00078 and pcorr = 0.0059, respectively). These associations also showed gender specificity, as they were not significant in men and substantially reinforced in women. The strongest association with MS severity was observed in women for combination (MIR499A*C/T + MIR196A2*C): pcorr = 4.43 × 10(-6) and OR = 3.23 (CI: 1.99-5.26).