Integrative analysis revealed potential causal genetic and epigenetic factors for multiple sclerosis

Biological sex affects functional variation across the human genome

Humans display sexual dimorphism across many traits, but little is known about underlying genetic mechanisms and impacts on disease. We utilized single-cell RNA-seq of 480 lymphoblastoid cell lines to reveal that the vast majority (79%) of sex-biased genes are targets of transcription factors that display sex-biased expression. Further, we developed a two-step regression method that identified sex-biased expression quantitative trait loci (sb-eQTL) across the genome. In contrast to previous work, these sb-eQTL are abundant (n=10,754; FDR 5%) and reproducible (replication up to π1=0.56). These sb-eQTL are enriched in over 600 GWAS phenotypes, including 120 sb-eQTL associated with the female-biased autoimmune disease multiple sclerosis. Our results demonstrate widespread genetic impacts on sexual dimorphism and identify possible mechanisms and clinical targets for sex differences in diverse diseases.

The Roles of RNA N6-methyladenosine Modifications in Systemic Lupus Erythematosus

N6-methyladenosine (m6A) modification is the most widespread RNA internal modification involved in RNA metabolism. M6A regulators consist of writers, erasers and readers. They exert their function by methylation, demethylation and recognization respectively, participating in cell biology and immune responses. Previously, the focus of m6A modification is its effect on tumor progress. Currently, extensive m6A-related studies have been performed in autoimmune diseases, such as RA, IBD and SLE, revealing that the unique influence of m6A modification in autoimmunity is undeniable. In this review, we summarize the function of m6A regulators, analyze their roles in pathogenic immune cells, summarize the m6A modification in SLE, and provide the potential m6A-targeting therapies for autoimmune diseases.

Identification of key genes with abnormal RNA methylation modification and selected m6A regulators in ankylosing spondylitis

BACKGROUND

N6-methyladenosine (m6A) has been identified as the most abundant modification of RNA molecules and the aberrant m6A modifications have been associated with the development of autoimmune diseases. However, the role of m6A modification in ankylosing spondylitis (AS) has not been adequately investigated. Therefore, we aimed to explore the significance of m6A regulator-mediated RNA methylation in AS.

METHODS

The methylated RNA immunoprecipitation sequencing (meRIP-seq) and digital RNA sequencing (Digital RNA-seq) were conducted using the peripheral blood mononuclear cells from three AS cases and three healthy controls, to identify genes affected by abnormal RNA methylation. The genes associated with different peaks were cross-referenced with AS-related genes obtained from the GeneCards Suite. Subsequently, the expression levels of shared differentially expressed genes (DEGs) and key m6A regulators in AS were evaluated using data from 68 AS cases and 36 healthy controls from two data sets (GSE25101 and GSE73754). In addition, the results were validated through quantitative polymerase chain reaction (qPCR).

RESULTS

The meRIP-seq and Digital RNA-seq analyses identified 28 genes with upregulated m6A peaks but with downregulated expression, and 52 genes with downregulated m6A peaks but with upregulated expression. By intersecting the genes associated with different peaks with 2184 AS-related genes from the GeneCards Suite, we identified a total of five shared DEGs: BCL11B, KAT6B, IL1R1, TRIB1, and ALDH2. Through analysis of the data sets and qPCR, we found that BCL11B and IL1R1 were differentially expressed in AS. Moreover, two key m6A regulators, WTAP and heterogeneous nuclear ribonucleoprotein C, were identified.

CONCLUSIONS

In conclusion, the current study revealed that m6A modification plays a crucial role in AS and might hence provide a new treatment strategy for AS disease.

Exploring the impact of m6A modification on immune diseases: mechanisms and therapeutic implication

N6-methyladenosine (m6A) is a chemical modification of RNA and has become a widely discussed topic among scientific researchers in recent years. It is distributed in various organisms, including eukaryotes and bacteria. It has been found that m6A is composed of writers, erasers and readers and is involved in biological functions such as splicing, transport and translation of RNA. The balance of the human immune microenvironment is important for human health abnormalities. Increasing studies have found that m6A affects the development of immune diseases such as inflammatory enteritis and systemic lupus erythematosus (SLE) by participating in the homeostatic regulation of the immune microenvironment in vivo. In this manuscript, we introduce the composition, biological function, regulation of m6A in the immune microenvironment and its progression in various immune diseases, providing new targets and directions for the treatment of immune diseases in clinical practice.

Potential biomarkers for multiple sclerosis stage from targeted proteomics and microRNA sequencing

Multiple sclerosis is a chronic demyelinating disease of the central nervous system. There is a need for new circulating biomarkers for multiple sclerosis, in particular, markers that differentiate multiple sclerosis subtypes (relapsing-remitting, secondary progressive and primary progressive multiple sclerosis), as this can help in making treatment decisions. In this study, we explore two classes of potential multiple sclerosis biomarkers-proteins and microRNAs-circulating in the cerebrospinal fluid and serum. Targeted medium-throughput proteomics (92 proteins) and microRNA sequencing were performed on serum samples collected in a cross-sectional case-control cohort (cohort I, controls n = 30, multiple sclerosis n = 75) and a prospective multiple sclerosis cohort (cohort II, n = 93). For cohort I, we also made these measurements in paired cerebrospinal fluid samples. In the cohort I cerebrospinal fluid, we observed differences between multiple sclerosis and controls for 13 proteins, including some previously described to be markers for multiple sclerosis [e.g. CD27, C-X-C motif chemokine 13 (CXCL13) and interleukin-7 (IL7)]. No microRNAs were significantly differentially expressed between multiple sclerosis and controls in the cerebrospinal fluid. In serum, 10 proteins, including angiopoietin-1 receptor (TIE2), and 16 microRNAs were significantly different between relapsing-remitting multiple sclerosis and secondary progressive multiple sclerosis after performing a meta-analysis combining both cohorts. In the prospective part of the study, participants with relapsing-remitting multiple sclerosis were followed for around 3 years, during which time 12 participants converted to secondary progressive multiple sclerosis. In these longitudinally collected serum samples, we observed a peak in granzyme B, A and H proteins around the time of conversion. Single-sample enrichment analysis of serum microRNA profiles revealed that the peak in granzyme B levels around conversion coincides with enrichment for microRNAs that are enriched in CD4+, CD8+ and natural killer cells (e.g. miRNA-150). We identified several proteins and microRNAs in serum that represent potential biomarkers for relapsing-remitting and secondary progressive multiple sclerosis. Conversion to secondary progressive disease is marked by a peak in granzyme B levels and enrichment for immune-related microRNAs. This indicates that specific immune cell-driven processes may contribute to the conversion of relapsing-remitting multiple sclerosis to secondary progressive multiple sclerosis.

DDR1-targeted therapies: current limitations and future potential

Discoidin domain receptor (DDR)-1 has a crucial role in regulating vital processes, including cell differentiation, proliferation, adhesion, migration, invasion, and matrix remodeling. Overexpression or activation of DDR1 in various pathological scenarios makes it a potential therapeutic target for the treatment of cancer, fibrosis, atherosclerosis, and neuropsychiatric, psychiatric, and neurodegenerative disorders. In this review, we summarize current therapeutic approaches targeting DDR1 from a medicinal chemistry perspective. Furthermore, we analyze factors other than issues of low selectivity and risk of resistance, contributing to the infrequent success of DDR1 inhibitors. The complex interplay between DDR1 and the extracellular matrix (ECM) necessitates additional validation, given that DDR1 might exhibit complex and synergistic interactions with other signaling molecules during ECM regulation. The mechanisms involved in DDR1 regulation in cancer and inflammation-related diseases also remain unknown.

The role of m6A RNA methylation in autoimmune diseases: Novel therapeutic opportunities

N6-methyladenosine (m6A) modifications, as one of the most common forms of internal RNA chemical modifications in eukaryotic cells, have gained increasing attention in recent years. The m6A RNA modifications exert various crucial roles in various biological processes, such as embryonic development, neurogenesis, circadian rhythms, and tumorigenesis. Recent advances have highlighted that m6A RNA modification plays an important role in immune response, especially in the initiation and progression of autoimmune diseases. In this review, we summarized the regulatory mechanisms of m6A methylation and its biological functions in the immune system and mainly focused on recent progress in research on the potential role of m6A RNA methylation in the pathogenesis of autoimmune diseases, thus providing possible biomarkers and potential targets for the prevention and treatment of autoimmune diseases.

M6A methylation modification in autoimmune diseases, a promising treatment strategy based on epigenetics

BACKGROUND

N6-methyladenosine (m6A) methylation modification is involved in the regulation of various biological processes, including inflammation, antitumor, and antiviral immunity. However, the role of m6A modification in the pathogenesis of autoimmune diseases has been rarely reported.

METHODS

Based on a description of m6A modification and the corresponding research methods, this review systematically summarizes current insights into the mechanism of m6A methylation modification in autoimmune diseases, especially its contribution to rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE).

RESULTS

By regulating different biological processes, m6A methylation is involved in the pathogenesis of autoimmune diseases and provides a promising biomarker for the diagnosis and treatment of such diseases. Notably, m6A methylation modification is involved in regulating a variety of immune cells and mitochondrial energy metabolism. In addition, m6A methylation modification plays a role in the pathological processes of RA, and m6A methylation-related genes can be used as potential targets in RA therapy.

CONCLUSIONS

M6A methylation modification plays an important role in autoimmune pathological processes such as RA and SLE and represents a promising new target for clinical diagnosis and treatment, providing new ideas for the treatment of autoimmune diseases by targeting m6A modification-related pathways.

In Search of a Function for the N6-Methyladenosine in Epitranscriptome, Autophagy and Neurodegenerative Diseases

Changes in epitranscriptome with N6-methyladenine (m6A) modification could be involved in the development of multiple diseases, which might be a prevalent modification of messenger RNAs (mRNAs) in eukaryotes. The m6A modification might be performed through the action of methyltransferases, demethylases, and methylation-binding proteins. Importantly, the m6A methylation may be associated with various neurological disorders including Alzheimer's disease (AD), Parkinson's disease (PD), depression, aging-related diseases, and/or aging itself. In addition, the m6A methylation might functionally regulate the eukaryotic transcriptome by influencing the splicing, export, subcellular localization, translation, stability, and decay of mRNAs. Neurodegenerative diseases may possess a wide variety of phenotypes, depending on the neurons that degenerate on occasion. Interestingly, an increasing amount of evidence has indicated that m6A modification could modulate the expression of autophagy-related genes and promote autophagy in neuronal cells. Oxidative stresses such as reactive oxygen species (ROS) could stimulate the m6A RNA methylation, which may also be related to the regulation of autophagy and/or the development of neurodegenerative diseases. Both m6A modification and autophagy could also play critical roles in regulating the health condition of neurons. Therefore, a comprehensive understanding of the m6A and autophagy relationship in human diseases may benefit in developing therapeutic strategies in the future. This paper reviews advances in the understanding of the regulatory mechanisms of m6A modification in the occurrence and development of neurodegenerative diseases and/or aging, discussing the possible therapeutic procedures related to mechanisms of m6A RNA methylation and autophagy.

METTL1 promotes tumorigenesis through tRNA-derived fragment biogenesis in prostate cancer

Newly growing evidence highlights the essential role that epitranscriptomic marks play in the development of many cancers; however, little is known about the role and implications of altered epitranscriptome deposition in prostate cancer. Here, we show that the transfer RNA N7-methylguanosine (m7G) transferase METTL1 is highly expressed in primary and advanced prostate tumours. Mechanistically, we find that METTL1 depletion causes the loss of m7G tRNA methylation and promotes the biogenesis of a novel class of small non-coding RNAs derived from 5'tRNA fragments. 5'tRNA-derived small RNAs steer translation control to favour the synthesis of key regulators of tumour growth suppression, interferon pathway, and immune effectors. Knockdown of Mettl1 in prostate cancer preclinical models increases intratumoural infiltration of pro-inflammatory immune cells and enhances responses to immunotherapy. Collectively, our findings reveal a therapeutically actionable role of METTL1-directed m7G tRNA methylation in cancer cell translation control and tumour biology.

N6-Methyladenosine Directly Regulates CD40L Expression in CD4+ T Lymphocytes

T cell activation is a highly regulated process, modulated via the expression of various immune regulatory proteins including cytokines, surface receptors and co-stimulatory proteins. N⁶-methyladenosine (m⁶A) is an RNA modification that can directly regulate RNA expression levels and it is associated with various biological processes. However, the function of m⁶A in T cell activation remains incompletely understood. We identify m⁶A as a novel regulator of the expression of the CD40 ligand (CD40L) in human CD4⁺ lymphocytes. Manipulation of the m⁶A 'eraser' fat mass and obesity-associated protein (FTO) and m⁶A 'writer' protein methyltransferase-like 3 (METTL3) directly affects the expression of CD40L. The m⁶A 'reader' protein YT521-B homology domain family-2 (YTHDF2) is hypothesized to be able to recognize and bind m⁶A specific sequences on the CD40L mRNA and promotes its degradation. This study demonstrates that CD40L expression in human primary CD4⁺ T lymphocytes is regulated via m⁶A modifications, elucidating a new regulatory mechanism in CD4⁺ T cell activation that could possibly be leveraged in the future to modulate T cell responses in patients with immune-related diseases.

Landscape analysis of m6A modification regulators related biological functions and immune characteristics in myasthenia gravis

BACKGROUND

N6-methyladenosine (m6A) modification has been recognized to play fundamental roles in the development of autoimmune diseases. However, the implication of m6A modification in myasthenia gravis (MG) remains largely unknown. Thus, we aimed to systematically explore the potential functions and related immune characteristics of m6A regulators in MG.

METHODS

The GSE85452 dataset with MG and healthy samples was downloaded from Gene Expression Omnibus (GEO) database. m6A modification regulators were manually curated. The targets of m6A regulators were obtained from m6A2Target database. The differential expressed m6A regulators in GSE85452 dataset were identified by "limma" package and were validated by RT-PCR. Function enrichment analysis of dysregulated m6A regulators was performed using "clusterProfiler" package. Correlation analysis was applied for analyzing the relationships between m6A regulators and immune characteristics. Unsupervised clustering analysis was used to identify distinct m6A modification subtypes. The differences between subtypes were analyzed, including the expression level of all genes and the enrichment degree of immune characteristics. Weighted gene co-expression network analysis (WGCNA) was conducted to obtain modules associated with m6A modification subtypes.

RESULTS

We found that CBLL1, RBM15 and YTHDF1 were upregulated in MG samples of GSE85452 dataset, and the results were verified by RT-PCR in blood samples from19 MG patients and 19 controls. The targeted genes common modified by CBLL1, RBM15, and YTHDF1 were mainly enriched in histone modification and Wnt signaling pathway. Correlation analysis showed that three dysregulated m6A regulators were closely associated with immune characteristics. Among them, RBM15 possessed the strongest correlation with immune characteristics, including CD56dim natural killer cell (r = 0.77, P = 0.0023), T follicular helper cell (r = - 0.86, P = 0.0002), Interferon Receptor (r = 0.78, P = 0.0017), and HLA-DOA (r = 0.64, P = 0.0200). Further two distinct m6A modification patterns mediated by three dysregulated m6A regulators was identified. Bioinformatics analysis found that there were 3029 differentially expressed genes and different immune characteristics between two m6A modification patterns. Finally, WGCNA analysis obtained a total of 12 modules and yellow module was the most positively correlated to subtype-2.

CONCLUSION

Our findings suggested that m6A RNA modification had an important effect on immunity molecular mechanism of MG and provided a new perspective into understanding the pathogenesis of MG.

Cross-sectional and longitudinal assessment of the association between DDR1 variants and processing speed in patients with early psychosis and healthy controls

Recent evidence indicates that DDR1 participates in myelination and that variants of DDR1 are associated with decreased cognitive processing speed (PS) in schizophrenia (SZ). Here, we explored whether DDR1 variants were associated with PS in subjects diagnosed with an early psychosis (EP), a condition often preceding SZ. Data from two Spanish independent samples (from Reus and Santander) including patients with EP (n = 75 and n = 312, respectively) and healthy controls (HCs; n = 57 and n = 160) were analyzed. The Trail Making Test part A was used to evaluate PS. Participants underwent genotyping to identify DDR1 variants rs1264323 and rs2267641. Cross-sectional data were analyzed with general linear models and longitudinal data were analyzed using mixed models. We examined the combined rs1264323AA-rs2267641AC/CC genotypes (an SZ-risk combination) on PS. The SZ-risk combined genotypes were associated with increased PS in EP patients but not in HCs in the cross-sectional analysis. In the longitudinal analysis, the SZ-risk combined genotypes were significantly associated with increased PS in both HCs and EP patients throughout the 10-year follow-up but no genotype × time interaction was observed. These results provide further evidence that DDR1 is involved in cognition and should be replicated with other samples.

Differential transcriptomic changes in the central nervous system and urinary bladders of mice infected with a coronavirus

Multiple sclerosis (MS) often leads to the development of neurogenic lower urinary tract symptoms (LUTS). We previously characterized neurogenic bladder dysfunction in a mouse model of MS induced by a coronavirus, mouse hepatitis virus (MHV). The aim of the study was to identify genes and pathways linking neuroinflammation in the central nervous system with urinary bladder (UB) dysfunction to enhance our understanding of the mechanisms underlying LUTS in demyelinating diseases. Adult C57BL/6 male mice (N = 12) received either an intracranial injection of MHV (coronavirus-induced encephalomyelitis, CIE group), or sterile saline (control group). Spinal cord (SC) and urinary bladders (UB) were collected from CIE mice at 1 wk and 4 wks, followed by RNA isolation and NanoString nCounter Neuroinflammation assay. Transcriptome analysis of SC identified a significantly changed expression of >150 genes in CIE mice known to regulate astrocyte, microglia and oligodendrocyte functions, neuroinflammation and immune responses. Two genes were significantly upregulated (Ttr and Ms4a4a), and two were downregulated (Asb2 and Myct1) only in the UB of CIE mice. Siglec1 and Zbp1 were the only genes significantly upregulated in both tissues, suggesting a common transcriptomic link between neuroinflammation in the CNS and neurogenic changes in the UB of CIE mice.

The epitranscriptome in ageing and stress resistance: A systematic review

Modifications of RNA, collectively called the "epitranscriptome", might provide novel biomarkers and innovative targets for interventions in geroscience but are just beginning to be studied in the context of ageing and stress resistance. RNA modifications modulate gene expression by affecting translation initiation and speed, miRNA binding, RNA stability, and RNA degradation. Nonetheless, the precise underlying molecular mechanisms and physiological consequences of most alterations of the epitranscriptome are still only poorly understood. We here systematically review different types of modifications of rRNA, tRNA and mRNA, the methodology to analyze them, current challenges in the field, and human disease associations. Furthermore, we compiled evidence for a connection between individual enzymes, which install RNA modifications, and lifespan in yeast, worm and fly. We also included resistance to different stressors and competitive fitness as search criteria for genes potentially relevant to ageing. Promising candidates identified by this approach include RCM1/NSUN5, RRP8, and F33A8.4/ZCCHC4 that introduce base methylations in rRNA, the methyltransferases DNMT2 and TRM9/ALKBH8, as well as factors involved in the thiolation or A to I editing in tRNA, and finally the m⁶A machinery for mRNA.

Alterations of m6A RNA methylation regulators contribute to autophagy and immune infiltration in primary Sjögren's syndrome

N6-methyladenosine (m6A) RNA modification is a new epigenetic regulation mechanism on eukaryotic mRNA. Few autoimmune diseases focused on the role of m6A in their pathogenies, and m6A modulation in the pathological process of primary Sjögren's syndrome (pSS) is still unknown. In this work, three microarray datasets of pSS patients were downloaded from the GEO database: datasets #1 and #2 from the whole peripheral blood (PB) samples, dataset #3 from the labial salivary gland tissue samples, as well as a PB cohort collected from our hospital. Six differentially expressed m6A regulators were identified by comparing the PB dataset #1 of pSS and healthy controls using the Wilcox test and logistic regression analysis. Among them, four (ALKBH5, RBMX, RBM15B, and YTHDF1) were confirmed as down-regulated in PB dataset #2 and in our PB cohort by RT-PCR, and four (ALKBH5, METTL3, RBM15B, and YTHDF1) were confirmed as down-regulated in the dataset #3 of the labial gland tissue. In addition, discrepantly expressed m6A regulators accompanied by diverse immunocytes, including dendritic cells (DCs), T cells, and CD56dim natural killer cells, and among the regulators, ALKBH5 and METTL3 were comprehensively linked with the infiltrated immune cells. Notably, the most enriched autophagy mechanism mediated by m6A was observed in pSS using functional annotation analysis. Ten hub genes were identified using a protein-protein interaction network, and their expression in PB dataset #2 and the expression of three genes (PIK3CA, STAT1, and MAPK3) in the labial gland tissue dataset #3 were confirmed. Our study provides evidence that m6A methylation is widely involved in the immune infiltration and autophagy of pSS, thus contributing to the pathogenesis of this disease and potentially representing a novel therapeutic target.

The Role of Oxidative Stress in Epigenetic Changes Underlying Autoimmunity

Significance: Epigenetic dysregulation plays an important role in the pathogenesis and development of autoimmune diseases. Oxidative stress is associated with autoimmunity and is also known to alter epigenetic mechanisms. Understanding the interplay between oxidative stress and epigenetics will provide insights into the role of environmental triggers in the development of autoimmunity in genetically susceptible individuals. Recent Advances: Abnormal DNA and histone methylation patterns in genes and pathways involved in interferon and tumor necrosis factor signaling, cellular survival, proliferation, metabolism, organ development, and autoantibody production have been described in autoimmunity. Inhibitors of DNA and histone methyltransferases showed potential therapeutic effects in animal models of autoimmune diseases. Oxidative stress can regulate epigenetic mechanisms via effects on DNA damage repair mechanisms, cellular metabolism and the local redox environment, and redox-sensitive transcription factors and pathways. Critical Issues: Studies looking into oxidative stress and epigenetics in autoimmunity are relatively limited. The number of available longitudinal studies to explore the role of DNA methylation in the development of autoimmune diseases is small. Future Directions: Exploring the relationship between oxidative stress and epigenetics in autoimmunity will provide clues for potential preventative measures and treatment strategies. Inception cohorts with longitudinal follow-up would help to evaluate epigenetic marks as potential biomarkers for disease development, progression, and treatment response in autoimmunity. Antioxid. Redox Signal. 36, 423-440.

A Final Frontier in Environment-Genome Interactions? Integrated, Multi-Omic Approaches to Predictions of Non-Communicable Disease Risk

Epidemiological and associative research from humans and animals identifies correlations between the environment and health impacts. The environment-health inter-relationship is effected through an individual's underlying genetic variation and mediated by mechanisms that include the changes to gene regulation that are associated with the diversity of phenotypes we exhibit. However, the causal relationships have yet to be established, in part because the associations are reduced to individual interactions and the combinatorial effects are rarely studied. This problem is exacerbated by the fact that our genomes are highly dynamic; they integrate information across multiple levels (from linear sequence, to structural organisation, to temporal variation) each of which is open to and responds to environmental influence. To unravel the complexities of the genomic basis of human disease, and in particular non-communicable diseases that are also influenced by the environment (e.g., obesity, type II diabetes, cancer, multiple sclerosis, some neurodegenerative diseases, inflammatory bowel disease, rheumatoid arthritis) it is imperative that we fully integrate multiple layers of genomic data. Here we review current progress in integrated genomic data analysis, and discuss cases where data integration would lead to significant advances in our ability to predict how the environment may impact on our health. We also outline limitations which should form the basis of future research questions. In so doing, this review will lay the foundations for future research into the impact of the environment on our health.

N6-methyladenosine and Neurological Diseases

N6-methyladenosine (m6A) is a dynamic reversible methylation modification of the adenosine N6 position and is the most common chemical epigenetic modification among mRNA post-transcriptional modifications, including methylation, demethylation, and recognition. Post-transcriptional modification involves multiple protein molecules, including METTL3, METTL14, WTAP, KIAA1429, ALKBH5, YTHDF1/2/3, and YTHDC1/2. m6A-related proteins are expressed in almost all cells. However, the abnormal expression of m6A-related proteins may occur in the nervous system, thereby affecting neuritogenesis, brain volume, learning and memory, memory formation and consolidation, etc., and is implicated in the development of diseases, such as Parkinson's disease, Alzheimer's disease, multiple sclerosis, depression, epilepsy, and brain tumors. This review focuses on the functions of m6A in the development of central nervous system diseases, thus contributing to a deeper understanding of disease pathogenesis and providing potential clinical therapeutic targets for neurological diseases.

Genetics and functional genomics of multiple sclerosis

Multiple sclerosis (MS) is an inflammatory neurodegenerative disease with genetic predisposition. Over the last decade, genome-wide association studies with increasing sample size led to the discovery of robustly associated genetic variants at an exponential rate. More than 200 genetic loci have been associated with MS susceptibility and almost half of its heritability can be accounted for. However, many challenges and unknowns remain. Definitive studies of disease progression and endophenotypes are yet to be performed, whereas the majority of the identified MS variants are not yet functionally characterized. Despite these shortcomings, the unraveling of MS genetics has opened up a new chapter on our understanding MS causal mechanisms.

Genome-Wide Identification of N6-Methyladenosine Associated SNPs as Potential Functional Variants for Type 1 Diabetes

OBJECTIVES

N6-methyladenosine (m6A) is essential in the regulation of the immune system, but the role that its single nucleotide polymorphisms (SNPs) play in the pathogenesis of type 1 diabetes (T1D) remains unknown. This study demonstrated the association between genetic variants in m6A regulators and T1D risk based on a case-control study in a Chinese population.

METHODS

The tagging SNPs in m6A regulators were genotyped in 1005 autoantibody-positive patients with T1D and 1257 controls using the Illumina Human OmniZhongHua-8 platform. Islet-specific autoantibodies were examined by radioimmunoprecipitation in all the patients. The mixed-meal glucose tolerance test was performed on 355 newly diagnosed patients to evaluate their residual islet function. The functional annotations for the identified SNPs were performed in silico. Using 102 samples from a whole-genome expression microarray, key signaling pathways associated with m6A regulators in T1D were comprehendingly evaluated.

RESULTS

Under the additive model, we observed three tag SNPs in the noncoding region of the PRRC2A (rs2260051, rs3130623) and YTHDC2 (rs1862315) gene are associated with T1D risk. Although no association was found between these SNPs and islet function, patients carrying risk variants had a higher positive rate for ZnT8A, GADA, and IA-2A. Further analyses showed that rs2260051[T] was associated with increased expression of PRRC2A mRNA (P = 7.0E-13), and PRRC2A mRNA was significantly higher in peripheral blood mononuclear cell samples from patients with T1D compared to normal samples (P = 0.022). Enrichment analyses indicated that increased PRRC2A expression engages in the most significant hallmarks of cytokine-cytokine receptor interaction, cell adhesion and chemotaxis, and neurotransmitter regulation pathways. The potential role of increased PRRC2A in disrupting immune homeostasis is through the PI3K/AKT pathway and neuro-immune interactions.

CONCLUSION

This study found intronic variants in PRRC2A and YTHDC2 associated with T1D risk in a Chinese Han population. PRRC2A rs2260051[T] may be implicated in unbalanced immune homeostasis by affecting the expression of PRRC2A mRNA. These findings enriched our understanding of m6A regulators and their intronic SNPs that underlie the pathogenesis of T1D.

The Emerging Role of m6A Modification in Regulating the Immune System and Autoimmune Diseases

Over the past several decades, RNA modifications have rapidly emerged as an indispensable topic in epitranscriptomics. N6-methyladenosine (m6A), namely, methylation at the sixth position of an adenine base in an RNA molecule, is the most prevalent RNA modification in both coding and noncoding RNAs. m6A has emerged as a crucial posttranscriptional regulator involved in both physiological and pathological processes. Based on accumulating evidence, m6A participates in the pathogenesis of immune-related diseases by regulating both innate and adaptive immune cells through various mechanisms. Autoimmune diseases are caused by a self-destructive immune response in the setting of genetic and environmental factors, and recent studies have discovered that m6A may play an essential role in the development of autoimmune diseases. In this review, we focus on the important role of m6A modification in biological functions and highlight its contributions to immune cells and the development of autoimmune diseases, thereby providing promising epitranscriptomic targets for preventing and treating autoimmune disorders.

Exploring the Epigenetic Regulatory Role of m6A-Associated SNPs in Type 2 Diabetes Pathogenesis

Purpose

Genetic factors in type 2 diabetes (T2D) pathogenesis have been widely explored by the genome-wide association studies (GWAS), identifying a great amount of susceptibility loci. With the development of high-resolution sequencing, the N(6)-methyladenosine (m6A) RNA modification has been proved to be affected by genetic variation. In this study, we identified the T2D-associated m6A-SNPs from T2D GWAS data and explored the underlying mechanism of the pathogenesis of T2D.

Methods

We examined the association of m6A-SNPs with T2D among large-scale T2D GWAS summary statistics and further performed multi-omics integrated analysis to explore the potential role of the identified m6A-SNPs in T2D pathogenesis.

Results

Among the 15,124 T2D-associated m6A-SNPs, 71 of them reach the genome-wide significant threshold (5.0e-05). The leading SNP rs4993986 (C>G), which is located near the m6A modification site at the 3' end of the HLA-DQB1 transcript, is expected to participate in the pathogenesis of T2D by influencing m6A modification to regulate the HLA-DQB1 expression.

Conclusion

The current study has suggested a potential correlation between m6A-SNPs and T2D pathogenesis and also provided new insights into the pathogenic mechanism of the T2D susceptibility loci identified by GWAS.

A spotter's guide to SNPtic exons: The common splice variants underlying some SNP-phenotype correlations

BACKGROUND

Cryptic exons are typically characterised as deleterious splicing aberrations caused by deep intronic mutations. However, low-level splicing of cryptic exons is sometimes observed in the absence of any pathogenic mutation. Five recent reports have described how low-level splicing of cryptic exons can be modulated by common single-nucleotide polymorphisms (SNPs), resulting in phenotypic differences amongst different genotypes.

METHODS

We sought to investigate whether additional 'SNPtic' exons may exist, and whether these could provide an explanatory mechanism for some of the genotype-phenotype correlations revealed by genome-wide association studies. We thoroughly searched the literature for reported cryptic exons, cross-referenced their genomic coordinates against the dbSNP database of common SNPs, then screened out SNPs with no reported phenotype associations.

RESULTS

This method discovered five probable SNPtic exons in the genes APC, FGB, GHRL, MYPBC3 and OTC. For four of these five exons, we observed that the phenotype associated with the SNP was compatible with the predicted splicing effect of the nucleotide change, whilst the fifth (in GHRL) likely had a more complex splice-switching effect.

CONCLUSION

Application of our search methods could augment the knowledge value of future cryptic exon reports and aid in generating better hypotheses for genome-wide association studies.

Coexpression of the discoidin domain receptor 1 gene with oligodendrocyte-related and schizophrenia risk genes in the developing and adult human brain

BACKGROUND

Discoidin domain receptor tyrosine kinase 1 (DDR1) is present in multiple types of epithelial cells and is highly expressed in the nervous system. Previous studies have revealed that DDR1 is involved in schizophrenia (SCZ). Although the expression of DDR1 in oligodendrocytes has been described, its role in brain myelination is not well understood. In this study, we aimed to explore the coexpression network of DDR1 in the human brain and to compare the list of DDR1 coexpressing genes with the list of genes containing single nucleotide polymorphisms (SNPs) that are associated with SCZ.

MATERIALS AND METHODS

We used a weighted gene coexpression network analysis (WGCNA) of a dataset from four brain areas (the dorsolateral prefrontal cortex, primary motor cortex, hippocampus, and striatum) and from four different intervals (I) of life (I-1 = 10-38 weeks postconception, I-2 ≥0 to < 6 years, I-3 ≥ 6 to < 40 years, and I-4 ≥ 40 years of age). We compared the list of genes that are associated with SCZ in the GWAS Catalog with the list of genes coexpressing with DDR1 in each interval.

RESULTS

Our study revealed that DDR1 was coexpressed with oligodendrocyte-related genes mainly in I-2 (adjP = 5.66e-24) and I-3 (adjP = 2.8e-114), which coincided with the coexpression of DDR1 with myelination-related genes (adjP = 9.04e-03 and 2.51e-08, respectively). DDR1 was also coexpressed with astrocyte-related genes in I-1 (adjP = 1.11e-71), I-2 (adjP = 2.12e-20) and I-4 (adjP = 9.93e-52) and with type 2 microglia-related genes in I-1 (adjP = 2.84e-08), I-2 (adjP = 5.68e-16) and I-4 (adjP = 3.66e-10). Moreover, we observed significant enrichment of SCZ susceptibility genes within the coexpression modules containing DDR1 in I-1 and I-4 (P = 1e-04 and 0.0037, respectively), during which the DDR1 module showed the highest association with the astrocytes.

CONCLUSIONS

Our study confirmed that DDR1 is coexpressed with oligodendrocyte- and myelin-related genes in the human brain but suggests that DDR1 may contribute mainly to SCZ risk through its role in other glial cell types, such as astrocytes.

N6-Methyladenosine RNA modification in cerebrospinal fluid as a novel potential diagnostic biomarker for progressive multiple sclerosis

Background

Progressive multiple sclerosis (PMS) is an uncommon and severe subtype of MS that worsens gradually and leads to irreversible disabilities in young adults. Currently, there are no applicable or reliable biomarkers to distinguish PMS from relapsing–remitting multiple sclerosis (RRMS). Previous studies have demonstrated that dysfunction of N6-methyladenosine (m6A) RNA modification is relevant to many neurological disorders. Thus, the aim of this study was to explore the diagnostic biomarkers for PMS based on m6A regulatory genes in the cerebrospinal fluid (CSF).

Methods

Gene expression matrices were downloaded from the ArrayExpress database. Then, we identified differentially expressed m6A regulatory genes between MS and non-MS patients. MS clusters were identified by consensus clustering analysis. Next, we analyzed the correlation between clusters and clinical characteristics. The random forest (RF) algorithm was applied to select key m6A-related genes. The support vector machine (SVM) was then used to construct a diagnostic gene signature. Receiver operating characteristic (ROC) curves were plotted to evaluate the accuracy of the diagnostic model. In addition, CSF samples from MS and non-MS patients were collected and used for external validation, as evaluated by an m6A RNA Methylation Quantification Kit and by real-time quantitative polymerase chain reaction.

Results

The 13 central m6A RNA methylation regulators were all upregulated in MS patients when compared with non-MS patients. Consensus clustering analysis identified two clusters, both of which were significantly associated with MS subtypes. Next, we divided 61 MS patients into a training set (n = 41) and a test set (n = 20). The RF algorithm identified eight feature genes, and the SVM method was successfully applied to construct a diagnostic model. ROC curves revealed good performance. Finally, the analysis of 11 CSF samples demonstrated that RRMS samples exhibited significantly higher levels of m6A RNA methylation and higher gene expression levels of m6A-related genes than PMS samples.

Conclusions

The dynamic modification of m6A RNA methylation is involved in the progression of MS and could potentially represent a novel CSF biomarker for diagnosing MS and distinguishing PMS from RRMS in the early stages of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02981-5.

The m6A epitranscriptome on neural development and degeneration

N⁶-methyladenosine (m⁶A) is the most prevalent, conserved, and abundant RNA modification of the mRNAs of most eukaryotes, including mammals. Similar to epigenetic DNA modifications, m⁶A has been proposed to function as a critical regulator for gene expression. This modification is installed by m⁶A methylation "writers" (Mettl3/Mettl14 methyltransferase complex), and it can be reversed by demethylase "erasers" (Fto and Alkbh5). Furthermore, m⁶A can be recognized by "readers" (Ythdf and Ythdc families), which may be interpreted to affect mRNA splicing, stability, translation or localization. Levels of m⁶A methylation appear to be highest in the brain, where it plays important functions during embryonic stem cell differentiation, brain development, and neurodevelopmental disorders. Depletion of the m⁶A methylation writer Mettl14 from mouse embryonic nervous systems prolongs cell cycle progression of radial glia and extends cortical neurogenesis into postnatal stages. Recent studies further imply that dysregulated m⁶A methylation may be significantly correlated with neurodegenerative diseases. In this review, we give an overview of m⁶A modifications during neural development and associated disorders, and provide perspectives for studying m⁶A methylation.

Emerging Role of m6 A Methylome in Brain Development: Implications for Neurological Disorders and Potential Treatment

Dynamic modification of RNA affords proximal regulation of gene expression triggered by non-genomic or environmental changes. One such epitranscriptomic alteration in RNA metabolism is the installation of a methyl group on adenosine [N⁶-methyladenosine (m⁶A)] known to be the most prevalent modified state of messenger RNA (mRNA) in the mammalian cell. The methylation machinery responsible for the dynamic deposition and recognition of m⁶A on mRNA is composed of subunits that play specific roles, including reading, writing, and erasing of m⁶A marks on mRNA to influence gene expression. As a result, peculiar cellular perturbations have been linked to dysregulation of components of the mRNA methylation machinery or its cofactors. It is increasingly clear that neural tissues/cells, especially in the brain, make the most of m⁶A modification in maintaining normal morphology and function. Neurons in particular display dynamic distribution of m⁶A marks during development and in adulthood. Interestingly, such dynamic m⁶A patterns are responsive to external cues and experience. Specific disturbances in the neural m⁶A landscape lead to anomalous phenotypes, including aberrant stem/progenitor cell proliferation and differentiation, defective cell fate choices, and abnormal synaptogenesis. Such m⁶A-linked neural perturbations may singularly or together have implications for syndromic or non-syndromic neurological diseases, given that most RNAs in the brain are enriched with m⁶A tags. Here, we review the current perspectives on the m⁶A machinery and function, its role in brain development and possible association with brain disorders, and the prospects of applying the clustered regularly interspaced short palindromic repeats (CRISPR)–dCas13b system to obviate m⁶A-related neurological anomalies.

DNA Methylation and Expression Profiles of Whole Blood in Parkinson's Disease

Parkinson’s disease (PD) is the second most common age-related neurodegenerative disease. It is presently only accurately diagnosed at an advanced stage by a series of motor deficits, which are predated by a litany of non-motor symptoms manifesting over years or decades. Aberrant epigenetic modifications exist across a range of diseases and are non-invasively detectable in blood as potential markers of disease. We performed comparative analyses of the methylome and transcriptome in blood from PD patients and matched controls. Our aim was to characterize DNA methylation and gene expression patterns in whole blood from PD patients as a foundational step toward the future goal of identifying molecular markers that could predict, accurately diagnose, or track the progression of PD. We found that differentially expressed genes (DEGs) were involved in the processes of transcription and mitochondrial function and that PD methylation profiles were readily distinguishable from healthy controls, even in whole-blood DNA samples. Differentially methylated regions (DMRs) were functionally varied, including near transcription factor nuclear transcription factor Y subunit alpha (NFYA), receptor tyrosine kinase DDR1, RING finger ubiquitin ligase (RNF5), acetyltransferase AGPAT1, and vault RNA VTRNA2-1. Expression quantitative trait methylation sites were found at long non-coding RNA PAX8-AS1 and transcription regulator ZFP57 among others. Functional epigenetic modules were highlighted by IL18R1, PTPRC, and ITGB2. We identified patterns of altered disease-specific DNA methylation and associated gene expression in whole blood. Our combined analyses extended what we learned from the DEG or DMR results alone. These studies provide a foundation to support the characterization of larger sample cohorts, with the goal of building a thorough, accurate, and non-invasive molecular PD biomarker.

m6A RNA Methylation in Systemic Autoimmune Diseases-A New Target for Epigenetic-Based Therapy?

The general background of autoimmune diseases is a combination of genetic, epigenetic and environmental factors, that lead to defective immune reactions. This erroneous immune cell activation results in an excessive production of autoantibodies and prolonged inflammation. During recent years epigenetic mechanisms have been extensively studied as potential culprits of autoreactivity. Alike DNA and proteins, also RNA molecules are subjected to an extensive repertoire of chemical modifications. N6-methyladenosine is the most prevalent form of internal mRNA modification in eukaryotic cells and attracts increasing attention due to its contribution to human health and disease. Even though m6A is confirmed as an essential player in immune response, little is known about its role in autoimmunity. Only few data have been published up to date in the field of RNA methylome. Moreover, only selected autoimmune diseases have been studied in respect of m6A role in their pathogenesis. In this review, I attempt to present all available research data regarding m6A alterations in autoimmune disorders and appraise its role as a potential target for epigenetic-based therapies.

Emerging Perspectives of RNA N 6-methyladenosine (m6A) Modification on Immunity and Autoimmune Diseases

N ⁶-methyladenosine (m⁶A) modification, the addition of a methylation decoration at the position of N6 of adenosine, is one of the most prevalent modifications among the over 100 known chemical modifications of RNA. Numerous studies have recently characterized that RNA m⁶A modification functions as a critical post-transcriptional regulator of gene expression through modulating various aspects of RNA metabolism. In this review, we will illustrate the current perspectives on the biological process of m⁶A methylation. Then we will further summarize the vital modulatory effects of m⁶A modification on immunity, viral infection, and autoinflammatory disorders. Recent studies suggest that m⁶A decoration plays an important role in immunity, viral infection, and autoimmune diseases, thereby providing promising biomarkers and therapeutic targets for viral infection and autoimmune disorders.

Genome-Wide Detection of m6A-Associated Genetic Polymorphisms Associated with Ischemic Stroke

N⁶-Methyladenosine (m⁶A) methylation is the most abundant post-transcription modification in eukaryotes and plays a vital role in many pathological conditions including cerebral ischemia-reperfusion injury and vascular inflammation. Moreover, recent studies have reported that single-nucleotide polymorphisms (SNPs) can affect the m⁶A modification. Therefore, we investigated the relationship between m⁶A-SNPs and ischemic stroke (IS) risk through integrative analysis of an IS genome-wide association study and m⁶A-SNP list from the m6AVar database. Next, we performed eQTL and differential expression analysis to support these IS-associated m⁶A-SNPs. Finally, using the identified polymorphisms, a PPI network was constructed using the STRING database, and GO and pathway enrichment analyses were performed using the DAVID online tool. Accordingly, we identified 305 IS-associated SNPs that could affect m⁶A methylation. Next, 158 of these SNPs were determined to have eQTL signals on local genes. We further identified 84 local genes (containing a total of 87 SNPs) that were differentially expressed in IS patients. Finally, we identified several biological processes and pathways related to IS pathogenesis, such as "leukocyte migration" and "focal adhesion." In summary, our study detected dozens of m⁶A-SNPs as critical functional polymorphisms and novel genetic biomarkers for IS susceptibility and provided a new means of elucidating the biological mechanism underlying IS development.

Epitranscriptomic regulation by m6A RNA methylation in brain development and diseases

Cellular RNAs are pervasively tagged with diverse chemical moieties, collectively called epitranscriptomic modifications. The methylation of adenosine at N6 position generates N6-methyladenosine (m6A), which is the most abundant and reversible epitranscriptomic modification in mammals. The m6A signaling is mediated by a dedicated set of proteins comprised of writers, erasers, and readers. Contrary to the activation-repression binary view of gene regulation, emerging evidence suggests that the m6A methylation controls multiple aspects of mRNA metabolism, such as splicing, export, stability, translation, and degradation, culminating in the fine-tuning of gene expression. Brain shows the highest abundance of m6A methylation in the body, which is developmentally altered. Within the brain, m6A methylation is biased toward neuronal transcripts and sensitive to neuronal activity. In a healthy brain, m6A maintains several developmental and physiological processes such as neurogenesis, axonal growth, synaptic plasticity, circadian rhythm, cognitive function, and stress response. The m6A imbalance contributes to the pathogenesis of acute and chronic CNS insults, brain cancer, and neuropsychiatric disorders. This review discussed the molecular mechanisms of m6A regulation and its implication in the developmental, physiological, and pathological processes of the brain.

Summary-data-based Mendelian randomization prioritizes potential druggable targets for multiple sclerosis

Multiple sclerosis is a complex autoimmune disease caused by a combination of genetic and environmental factors. Translation of Genome-Wide Association Study findings into therapeutics and effective preventive strategies has been limited to date. We used summary-data-based Mendelian randomization to synthesize findings from public expression quantitative trait locus, methylation quantitative trait locus and Multiple Sclerosis Genome-Wide Association Study datasets. By correlating the effects of methylation on multiple sclerosis, methylation on expression and expression on multiple sclerosis susceptibility, we prioritize genetic loci with evidence of influencing multiple sclerosis susceptibility. We overlay these findings onto a list of 'druggable' genes, i.e. genes which are currently, or could theoretically, be targeted by therapeutic compounds. We use GeNets and search tool for the retrieval of interacting genes/proteins to identify protein-protein interactions and druggable pathways enriched in our results. We extend these findings to a model of Epstein-Barr virus-infected B cells, lymphoblastoid cell lines. We conducted a systematic review of prioritized genes using the Open Targets platform to identify completed and planned trials targeting prioritized genes in multiple sclerosis and related disease areas. Expression of 45 genes in peripheral blood was strongly associated with multiple sclerosis susceptibility (False discovery rate 0.05). Of these 45 genes, 20 encode a protein which is currently targeted by an existing therapeutic compound. These genes were enriched for Gene Ontology terms pertaining to immune system function and leucocyte signalling. We refined this prioritized gene list by restricting to loci where CpG site methylation was associated with multiple sclerosis susceptibility, with gene expression and where expression was associated with multiple sclerosis susceptibility. This approach yielded a list of 15 prioritized druggable target genes for which there was evidence of a pathway linking methylation, expression and multiple sclerosis. Five of these 15 genes are targeted by existing drugs and three were replicated in a smaller expression Quantitative Trait Loci dataset (CD40, MERTK and PARP1). In lymphoblastoid cell lines, this approach prioritized 7 druggable gene targets, of which only one was prioritized by the multi-omic approach in peripheral blood (FCRL3). Systematic review of Open Targets revealed multiple early-phase trials targeting 13/20 prioritized genes in disorders related to multiple sclerosis. We use public datasets and summary-data-based Mendelian randomization to identify a list of prioritized druggable genetic targets in multiple sclerosis. We hope our findings could be translated into a platform for developing targeted preventive therapies.

Genome-wide identification of m6A-associated single-nucleotide polymorphisms in Parkinson's disease

N⁶-methyladenosine (m⁶A)-associated single nucleotide polymorphisms (SNPs) play a vital role in several neurological diseases. However, little is known about the relationship between m⁶A modification and Parkinson's disease (PD). We investigated potential functional variants of m⁶A-SNPs from large-scale genome-wide association studies (GWAS) in PD patients. The candidate m⁶A-SNPs were further assessed by expression quantitative trait loci (eQTL) analysis and differential gene expression analysis. We identified 12 m⁶A-SNPs that were significantly associated with PD risk. Further, eQTL and expression analyses identified five of these m⁶A-SNPs (rs75072999 of GAK, rs1378602, rs4924839 and rs8071834 of ALKBH5, and rs1033500 of C6orf10) that were associated with altered gene expression in PD. Our results suggest that m⁶A-SNPs could play a role in PD risk. Future studies are needed to confirm these PD-associated m⁶A-SNPs and elucidate their mechanisms.