Abstract
Hirschsprung’s disease (HSCR) is a neural crest disease that results from the failure of enteric neural crest cells (ENCCs) to migrate to the corresponding intestinal segment. The RET gene, which regulates enteric neural crest cell proliferation and migration, is considered one of the main risk factors for HSCR and is commonly used to construct HSCR mouse models. The epigenetic mechanism of m6A modification is involved in HSCR. In this study, we analyzed the GEO database (GSE103070) for differentially expressed genes (DEGs) and focused on m6A–related genes. Comparing the RNA-seq data of Wide Type and RET Null, a total of 326 DEGs were identified, of which 245 genes were associated with m6A. According to the CIBERSORT analysis, the proportion of Memory B-cell in RET Null was significantly higher than that of Wide Type. Venn diagram analysis was used to identify key genes in the selected memory B-cell modules and DEGs associated with m6A. Enrichment analysis showed that seven genes were mainly involved in focal adhesion, HIV infection, actin cytoskeleton organization and regulation of binding. These findings could provide a theoretical basis for molecular mechanism studies of HSCR.
1 Introduction
HSCR is a common digestive malformation in children, with a prevalence of 1/5000 to 1/2000, ranking second in digestive malformations (). It is more prevalent in males, with an average male-to-female ratio of 4:1 (). The main clinical manifestations of HSCR include absent feces or delayed fetal excretion, bilious vomiting, feeding difficulties, and intractable constipation. Hirschsprung-associated enterocolitis (HAEC) occurs in some children, with a prevalence of 18.5% preoperatively and 18.2% postoperatively (). The symptoms include abdominal distention, fever, diarrhea, vomiting, bloody stools, and constipation (). However, some children exhibit atypical results in preoperative barium enema and rectal-anal manometry tests. Additionally, intestinal tissue damage with bleeding and perforation may occur after biopsy. The underlying cause of HSCR is the absence of ganglion cells in the intermuscular plexus of the intestinal wall due to arrested development of ENCCs migration. Mutations in genes, such as RET, EDNRB, and SOX10 have been associated with the pathogenesis of HSCR(). Among them, RET is the major risk factor for HSCR, which is involved in the proliferation, migration, and differentiation of intestinal neural crest cells (). The tyrosine kinase receptor encoded by the RET proto-oncogene plays a critical role in the development of the enteric nervous system. Aberrant expression of the RET gene can result in abnormal colonization of ENCCs within the intestine and even promote ENCCs apoptosis, ultimately leading to the occurrence of HSCR. RET knockout is a conventional method to build the model of HSCR in mice ().
N6-methyladenosine (m6A) methylation is a common epigenetic modification in eukaryotic, accounting for 80% of RNA methylation modifications, and it exerts biological functions by affecting RNA metabolism (). The formation of m6A is catalyzed by a methyltransferase complex (MTC), which comprises genes such as METTL3, METTL14, METTL16, and WTAP. The m6A demethylase, also known as the “eraser,” is involved in the removal of m6A, while RNA reader proteins recognize and bind to m6A, thus achieving corresponding functions. Related genes involved in this process include YTHSCRC2, YTHSCRF1, LRPPRC, HNRNPA2B1, and IGFBP1(). Methyltransferases can catalyze the m6A modification of adenosine on mRNA, while the role of demethylases is to demethylate the bases that have undergone m6A modification. The function of reader proteins is to recognize the bases that have undergone m6A modification, thus activating downstream regulatory pathways, such as RNA degradation and miRNA processing. The m6A demethylase, ALKBH5, has been found to inhibit the proliferation and migration of enteric neural crest cells by upregulating TAGLN and is involved in the pathogenesis of HSCR().
B lymphocytes (B-cell) are a complex cell type that plays a central role in humoral immunity against infection, autoimmunity and transplantation. B-cell activation leads to the formation of memory B-cell populations (). Memory B-cell can be classified based on the immunoglobulin markers they express, including IgM+, IgG+, IgA+, and IgE+ memory B-cell. Tissue-resident memory B-cell have been reported to be predominantly present in the intestine and characterized by the expression of CD45RB and CD69 (), providing the first line of defense against reinfection. IgM+ memory B-cell that are distributed throughout the intestine are closely related to IgM+ plasma cells in the gut and play a unique role in mucosal protection. B-cell can impede the interaction between epithelial and stromal cells and play a negative role in repairing intestinal damage ().
In this study, we obtained RNA-seq data of Wide Type and RET Null mouse intestinal samples from the GEO database (GSE103070), and analyzed the DEGs associated with m6A modification. We used weighted gene co-expression analysis and the CIBERSORT algorithm to screen for differential genes associated with m6A methylation modifications in memory B-cell of RET Null mouse samples. Enrichment analysis was performed to develop a theoretical framework for the diagnosis and study of HSCR.
2 Materials and methods
2.1 Data collection and processing
The Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/) provides publicly available data. We identified the GSE103070 dataset by searching the GEO database using keywords such as RET and Gut. The contributors of the GSE103070 dataset have confirmed significant changes in gene expression during late-stage development of the intestine (). By analyzing data from three time points, E10.5, E12.5, and E14.5, we confirmed that there were substantial differences in gene expression at E14.5. Therefore, we focused our analysis and discussion on the data from E14.5.
2.2 Differential expression gene analysis
In R (4.2.0) software, data from the above RNA-seq were analyzed using the R package “limma” to identify differentially expressed genes (DEGs). Significantly DEGs were detected according to the following criteria: (1) |log2FC| > 0.5, (2) False discovery rate (FDR) adjusted p-value <0.05. The “ggplots” and “pheatmap” packages in R were utilized to create volcano plots and heatmaps depicting the DEGs.
2.3 Identification of m6A-related genes
Pearson correlation analysis was conducted to evaluate the association between m6A modifications and DEGs in Wide Type and RET Null mouse intestinal tissues, with a threshold of |Pearson R|>0.5 and p < 0.05.
2.4 Immune cell abundance analysis
CIBERSORT is a general deconvolution algorithm based on gene expression, which can estimate the relative proportion of 22 types of immune cells (), such as M1 macrophages, M2 macrophages, memory B-cell and mast cells from gene expression profiles. The relative abundance of the 22 immune cell types in the Wide Type and RET Null samples was assessed by CIBERSORT.
2.5 WGCNA analysis
To explore the set of modular genes associated with the immune cells in RET Null sample, genes with mean FPKM>5 were selected for sample clustering. The soft threshold for subsequent network construction was calculated by pickSoftThreshold in the WGCNA package based on the correlation coefficient between the β value and the mean gene connectivity. The topological overlap matrix (TOM) was constructed using BlockwiseModules, and clustering was performed using the dynamic tree shearing algorithm. Each co-expression module was assigned a unique color to distinguish genes with similar expression patterns.
2.6 Venn diagram and enrichment analysis
The overlap between immune cell-related genes and m6A-related DEGs was examined using a Venn diagram generated with the online tool (https://bioinfogp.cnb.csic.es/tools/venny/). For enrichment analysis of genes in the identified modules and to explore their biological functions, Metascape (https://metascape.org/gp/index.html) was employed, which offers extensive gene function annotation and pathway analysis capabilities.
3 Results
3.1 Screening for differentially expressed genes
We identified 326 genes (Supplementary Material S1) as differentially expressed genes by screening the RNA-seq data of E14.5 in the GSE103070 dataset with |log2FC| >0.5 and adjusted p < 0.05, including 139 upregulated genes and 187 downregulated genes. A volcano map was created using DEG analysis results, with upregulated genes shown in red and downregulated genes shown in green (Figure 1A). The heat map was also plotted, with each column of the heat map representing a sample and each row representing a gene (Figure 1B). Similar genes and samples were clustered in the horizontal and vertical coordinates. We were able to observe that the expression patterns of functionally related genes were also similar.
FIGURE 1
3.2 Analysis of differentially expressed genes associated with m6A
Using Pearson correlation analysis, we identified 245 genes associated with m6A regulators by applying the screening conditions of |Pearson R|>0.5 and p < 0.05. M6A-modified genes were depicted with red dots in the figure, and differentially expressed genes were shown with blue dots (Figure 2). In summary, 245 genes were aberrantly expressed in HSCR and may be related to m6A methylation modification.
FIGURE 2
3.3 Assessment of immune cell abundance
In our study, we analyzed the abundance of immune cells in both Wide Type and RET Null samples using the CIBERSORT algorithm. R language was used to visualize the distribution of 22 immune cell types, with red representing Wide Type and green representing RET Null (Figure 3). Our analysis revealed a significant difference in the proportion of memory B-cell between RET Null and control samples, suggesting that memory B-cell may play an important role in the development of HSCR.
FIGURE 3
3.4 Construction and screening of WGCNA modules
We constructed a gene co-expression network using WGCNA to identify modules associated with HSCR samples with the highest correlation with immune cells. The scale-free topological fit index R2 was 0.9 and βwas 12 (Figure 4A). After identifying the gene modules using the dynamic tree-cutting method, we calculated the feature vector values for each module. Modules with similar distances were merged (Figure 4B). Different types of immune cells were represented by the horizontal coordinates, and different module genes were represented by the vertical coordinates. Darker colors represented higher correlations. Red represents positive correlations and green represents negative correlations. Since there is a significant difference in the proportion of memory B-cell between RET Null mouse intestinal samples and Wide Type samples, subsequent studies will focus on the brown module of memory B-cell genes (Figure 5A).
FIGURE 4
FIGURE 5
3.5 Identification of key genes
To identify differentially expressed genes associated with memory B-cell and related with m6A modifications, we mapped 245 DEGs associated with m6A and 641 genes associated with the most relevant modules of memory B-cell in a Venn diagram (Figure 6A), resulting in seven (COL23A1, FBLIM1, GNG10, MYOCD, NIN, NRARP, and ZKSCAN1) differentially expressed genes. COL23A1, FBLIM1, GNG10, MYOCD, NIN and ZKSCAN1 were highly expressed in the RET Null group, while NRARP was low expressed (Figure 6B). Enrichment analysis of the brown module genes (Figure 5B) revealed that these seven genes were mainly involved in biological processes, such as local adhesion, HIV infection, actin cytoskeleton organization, binding regulation, and microtubule morphogenesis.
FIGURE 6
4 Discussion
HSCR is a common congenital gastrointestinal disorder in children, characterized by a congenital absence of intrinsic ganglion cells in the muscular and submucosa plexuses of the gastrointestinal tract (). The mechanism of HSCR is the disturbance of the proliferation and migration of ENCCs during embryonic development or alteration of the intestinal microenvironment (). Genetic studies have found that genes, such as RET and EDNRB are associated with HSCR(). RET is a major risk factor for HSCR that encodes tyrosine kinase and is closely associated with the proliferation, migration, and differentiation of ENCCs. M6A modification is involved in the occurrence and progression of HSCR, and this study focuses on the possible role of m6A modification-related genes in HSCR after RET knockdown.
In this study, we analyzed the GSE103070 dataset from the GEO database using the CIBERSORT algorithm. Our results showed that the proportion of memory B-cell was significantly increased in intestinal samples from RET Null compared to the Wide Type. Using WGCNA, we identified gene modules with the highest correlation with immune cells and intersected them with m6A-related differentially expressed genes to identify seven key genes (COL23A1, FBLIM1, GNG10, MYOCD, NIN, NRARP and ZKSCAN1) significantly associated with m6A and memory B-cell.
Studies have reported that the expression of memory B-cell is significantly different in patients with ulcerative colitis compared to normal subjects (). In hepatoblastoma, the relative abundance of memory B-cell is positively correlated with disease severity (). In HSCR, abnormalities in gut microbial regulation usually cause HAEC, resulting in delayed thymic development. This study found that memory B-cell were more predominant in the RET Null HSCR model, suggesting their close association with HSCR. Therefore, we performed a follow-up analysis of memory B-cell.
Enrichment analysis of memory B cell-associated genes revealed that the seven candidate genes are mainly involved in several biological processes, including focal adhesion, HIV infection, actin cytoskeleton organization, regulation of binding, tube morphogenesis, and the VEGFA-VEGFR2 signaling pathway. By checking through the RMDisease2 website, we found that six of the key genes can be modified by m6A. In multiple sclerosis (MS), CD38+ memory B-cell in the blood migrate toward the central nervous system by expressing the adhesion molecule ALCAM(). The downregulation of p130Cas in ENCCs in the presence of collagen ColVI affects FN-induced actin polymerization and focal adhesion formation (). Notch-regulated ankyrin repeat-containing protein (NRARP) is a small protein with two anchor protein repeat sequences that can regulate neural crest cell differentiation by regulating LEF1 stability (). NRARP protein is required for the anterior-posterior somite pattern in mice, and the reduced expression of constitutive NRARP blocks early thymic cell maturation progression (). In our findings, NRARP expression was significantly reduced in intestinal samples from RET Null mice, suggesting that the neural crest cell differentiation involving NRARP is closely related to HSCR. Serine-binding LIM protein 1 (FBLIM1) is involved in cell attachment and its encoded Migfilin is essential for cellular processes, such as cell morphology and migration. FBLIM1 can co-localize with kindlin-1 and kindlin-2 on the adhesion spots of normal human keratin-forming cells (). The regulation of Migfilin expression is essential for neuronal development (). G protein-coupled membrane receptors regulate intercellular transmitters and hormone secretion in neuroendocrine cells. GNG10 encodes a subunit of G proteins involved in cell cycle regulation and is closely associated with the intestinal microflora. Overexpression of GNG10 promotes colorectal cancer progression (), and we found that GNG10 expression was higher in the RET Null, suggesting an association between GNG10 and HSCR. ZKSCAN1 acts as a cell cycle regulator to regulate cell progression (Zhu et al., 2019). MYOCD can inhibit the developmental program through the NOTCH signaling pathway, transforming vascular smooth muscle cells (VSMC) into a dedifferentiated state, leading to complete transformation into macrophage-like cells (Zhang et al., 2022). The COL23A1 gene can be used as a surrogate marker for monitoring disease activity in ulcerative colitis and Crohn’s disease (), as well as a biomarker for papillary thyroid cancer () and clear cell renal cell carcinoma ().
Our findings suggest that the key genes associated with memory B-cell and m6A modification may influence the development of HSCR by participating in the inflammatory response and affecting neural crest cell proliferation, migration, and differentiation. However, due to limited data availability on human samples in HSCR, further exploration is required. We will conduct clinical validation of the genes identified in this study to determine their potential as new therapeutic targets, providing a theoretical basis for clinical diagnosis and a deeper understanding of HSCR mechanisms.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
Conceptualization, JH; Data curation, JH; Formal analysis, JH, TC, JW, and ZW; Funding acquisition, SH; Investigation, JH and TC; Methodology, JW; Software, ZW and JW; Validation, JW and JH; Writing–original draft, JH and TC; Writing–review and editing, JH, TC, and SH.
Funding
The work was supported by Suzhou Key Discipline of Medicine (SZXK202105) and Suzhou Municipal Health and Family Planning Commission’s Key Diseases Diagnosis and Treatment Program (LCZX202107).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1183467/full#supplementary-material
References
1
BachettiT.CeccheriniI. (2020). Causative and common PHOX2B variants define a broad phenotypic spectrum. Clin. Genet.97 (1), 103–113. 10.1111/cge.13633
2
BertoutJ. A.PatelS. A.FryerB. H.DurhamA. C.CovelloK. L.OliveK. P.et al (2009). Heterozygosity for hypoxia inducible factor 1alpha decreases the incidence of thymic lymphomas in a p53 mutant mouse model. Cancer Res.69 (7), 3213–3220. 10.1158/0008-5472.CAN-08-4223
3
BrahmeN. N.HarburgerD. S.Kemp-O'BrienK.StewartR.RaghavanS.ParsonsM.et al (2013). Kindlin binds migfilin tandem LIM domains and regulates migfilin focal adhesion localization and recruitment dynamics. J. Biol. Chem.288 (49), 35604–35616. 10.1074/jbc.M113.483016
4
ChatterjeeS.NandakumarP.AuerD. R.GabrielS. B.ChakravartiA. (2019). Gene- and tissue-level interactions in normal gastrointestinal development and Hirschsprung disease. Proc. Natl. Acad. Sci. U. S. A.116 (52), 26697–26708. 10.1073/pnas.1908756116
5
FrascaD.DiazA.RomeroM.GarciaD.BlombergB. B. (2020). B cell immunosenescence. Annu. Rev. Cell Dev. Biol.36, 551–574. 10.1146/annurev-cellbio-011620-034148
6
FredeA.CzarnewskiP.MonasterioG.TripathiK. P.BejaranoD. A.Ramirez FloresR. O.et al (2022). B cell expansion hinders the stroma-epithelium regenerative cross talk during mucosal healing. Immunity55 (12), 2336–2351.e12. 10.1016/j.immuni.2022.11.002
7
GaoT.XuW.ShengQ.XuT.WuW.LvZ. (2022). Clinical outcomes and risk factors for postoperative complications in children with Hirschsprung's disease. Am. J. Transl. Res.14 (7), 4830–4837.
8
GosainA.FrykmanP. K.CowlesR. A.HortonJ.LevittM.RothsteinD. H.et al (2017). Guidelines for the diagnosis and management of Hirschsprung-associated enterocolitis. Pediatr. Surg. Int.33 (5), 517–521. 10.1007/s00383-017-4065-8
9
HagensJ.ReinshagenK.TomuschatC. (2022). Prevalence of Hirschsprung-associated enterocolitis in patients with Hirschsprung disease. Pediatr. Surg. Int.38 (1), 3–24. 10.1007/s00383-021-05020-y
10
HeuckerothR. O. (2018). Hirschsprung disease - integrating basic science and clinical medicine to improve outcomes. Nat. Rev. Gastroenterol. Hepatol.15 (3), 152–167. 10.1038/nrgastro.2017.149
11
IshitaniT.MatsumotoK.ChitnisA. B.ItohM. (2005). Nrarp functions to modulate neural-crest-cell differentiation by regulating LEF1 protein stability. Nat. Cell Biol.7 (11), 1106–1112. 10.1038/ncb1311
12
IshizukaK.TabataH.ItoH.KushimaI.NodaM.YoshimiA.et al (2018). Possible involvement of a cell adhesion molecule, Migfilin, in brain development and pathogenesis of autism spectrum disorders. J. Neurosci. Res.96 (5), 789–802. 10.1002/jnr.24194
13
JiangX.LiuB.NieZ.DuanL.XiongQ.JinZ.et al (2021). The role of m6A modification in the biological functions and diseases. Signal Transduct. Target Ther.6 (1), 74. 10.1038/s41392-020-00450-x
14
LeT.AronowR. A.KirshteinA.ShahriyariL. (2021). A review of digital cytometry methods: Estimating the relative abundance of cell types in a bulk of cells. Brief. Bioinform22 (4), bbaa219. 10.1093/bib/bbaa219
15
LiuS.ZhengQ.ZhangR.LiT.ZhanJ. (2022). Construction of a combined random forest and artificial neural network diagnosis model to screening potential biomarker for hepatoblastoma. Pediatr. Surg. Int.38 (12), 2023–2034. 10.1007/s00383-022-05255-3
16
Manon-JensenT.SunS.LindholmM.DomislovicV.GiuffridaP.BrinarM.et al (2021). Elevated ectodomain of type 23 collagen is a novel biomarker of the intestinal epithelium to monitor disease activity in ulcerative colitis and Crohn's disease. United Eur. Gastroenterol. J.9 (2), 268–278. 10.1177/2050640620977371
17
McWilliamO.SellebjergF.MarquartH. V.von EssenM. R. (2018). B cells from patients with multiple sclerosis have a pathogenic phenotype and increased LTα and TGFβ1 response. J. Neuroimmunol.324, 157–164. 10.1016/j.jneuroim.2018.09.001
18
NagyN.KovacsT.StavelyR.HalasyV.SoosA.SzocsE.et al (2021). Avian ceca are indispensable for hindgut enteric nervous system development. Development148 (22), dev199825. 10.1242/dev.199825
19
NishidaS.YoshizakiH.YasuiY.KuwaharaT.KiyokawaE.KohnoM. (2018). Collagen VI suppresses fibronectin-induced enteric neural crest cell migration by downregulation of focal adhesion proteins. Biochem. Biophys. Res. Commun.495 (1), 1461–1467. 10.1016/j.bbrc.2017.11.184
20
OkamotoM.UesakaT.ItoK.EnomotoH. (2021). Increased RET activity coupled with a reduction in the RET gene dosage causes intestinal aganglionosis in mice. eNeuro8 (3), ENEURO.0534–20.2021. 10.1523/ENEURO.0534-20.2021
21
PakarinenM. (2018). Perioperative complications of transanal pull-through surgery for hirschsprung's disease. Eur. J. Pediatr. Surg.28 (2), 152–155. 10.1055/s-0038-1632393
22
PanY.WuL.HeS.WuJ.WangT.ZangH. (2021). Identification of hub genes in thyroid carcinoma to predict prognosis by integrated bioinformatics analysis. Bioengineered12 (1), 2928–2940. 10.1080/21655979.2021.1940615
23
PararasaC.ZhangN.TullT. J.ChongM. H. A.SiuJ. H. Y.GuesdonW.et al (2019). Reduced CD27(-)IgD(-) B cells in blood and raised CD27(-)IgD(-) B cells in gut-associated lymphoid tissue in inflammatory bowel disease. Front. Immunol.10, 361. 10.3389/fimmu.2019.00361
24
SergiC. M.CaluseriuO.McCollH.EisenstatD. D. (2017). Hirschsprung's disease: Clinical dysmorphology, genes, micro-RNAs, and future perspectives. Pediatr. Res.81 (1-2), 177–191. 10.1038/pr.2016.202
25
SoretR.SchneiderS.BernasG.ChristophersB.SouchkovaO.CharrierB.et al (2020). Glial cell-derived neurotrophic factor induces enteric neurogenesis and improves colon structure and function in mouse models of hirschsprung disease. Gastroenterology159 (5), 1824–1838.e17. 10.1053/j.gastro.2020.07.018
26
WangB.FangX.SunX.DuC.ZhouL.LvX.et al (2021a). m(6)A demethylase ALKBH5 suppresses proliferation and migration of enteric neural crest cells by regulating TAGLN in Hirschsprung's disease. Life Sci.278, 119577. 10.1016/j.lfs.2021.119577
27
WangN.LiJ.HeJ.JingY. G.ZhaoW. D.YuW. J.et al (2021b). Knockdown of lncRNA CCAT1 inhibits the progression of colorectal cancer via hsa-miR-4679 mediating the downregulation of GNG10. J. Immunol. Res.2021, 8930813. 10.1155/2021/8930813
28
WeiselN. M.WeiselF. J.FarberD. L.BorghesiL. A.ShenY.MaW.et al (2020). Comprehensive analyses of B-cell compartments across the human body reveal novel subsets and a gut-resident memory phenotype. Blood136 (24), 2774–2785. 10.1182/blood.2019002782
29
XuF.ChangK.MaJ.QuY.XieH.DaiB.et al (2017). The oncogenic role of COL23A1 in clear cell renal cell carcinoma. Sci. Rep.7 (1), 9846. 10.1038/s41598-017-10134-2
30
ZaccaraS.RiesR. J.JaffreyS. R. (2019). Reading, writing and erasing mRNA methylation. Nat. Rev. Mol. Cell Biol.20 (10), 608–624. 10.1038/s41580-019-0168-5
31
ZhangZ.HuangJ.WangY.ShenW. (2022). Transcriptome analysis revealed a two-step transformation of vascular smooth muscle cells to macrophage-like cells. Atherosclerosis346, 26–35. 10.1016/j.atherosclerosis.2022.02.021
32
ZhuY. J.ZhengB.LuoG. J.MaX. K.LuX. Y.LinX. M.et al (2019). Circular RNAs negatively regulate cancer stem cells by physically binding FMRP against CCAR1 complex in hepatocellular carcinoma. Theranostics9 (12), 3526–3540. 10.7150/thno.32796
Summary
Keywords
hirschsprung’s disease, M6A, WGCNA, CIBERSORT, enrichment analysis
Citation
Huang J, Chen T, Wang J, Wang Z and Huang S (2023) Weighted gene co-expression network analysis and CIBERSORT screening of key genes related to m6A methylation in Hirschsprung’s disease. Front. Genet. 14:1183467. doi: 10.3389/fgene.2023.1183467
Received
10 March 2023
Accepted
06 April 2023
Published
18 April 2023
Volume
14 - 2023
Edited by
Yi Zheng, Northwest A&F University, China
Reviewed by
Xinxin Li, Northwestern Polytechnical University, China
Pengfei Zhang, Guangzhou National Laboratory, China
Updates
Copyright
© 2023 Huang, Chen, Wang, Wang and Huang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shungen Huang, drhuang01@163.com
† These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.