AUTHOR=Yang Jing , Xu Jing , Zhang Luting , Li Yingting , Chen Min 

TITLE=Identifying key m6A-methylated lncRNAs and genes associated with neural tube defects via integrative MeRIP and RNA sequencing analyses

JOURNAL=Frontiers in Genetics

VOLUME=13

YEAR=2022

URL=https://www.frontiersin.org/journals/genetics/articles/10.3389/fgene.2022.974357

DOI=10.3389/fgene.2022.974357

ISSN=1664-8021

ABSTRACT=<p><bold>Objective:</bold> N<sup>6</sup>-methyladenosine (m<sup>6</sup>A) is a common post-transcriptional modification of messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs). However, m<sup>6</sup>A-modified lncRNAs are still largely unexplored. This study aimed to investigate differentially m<sup>6</sup>A-modified lncRNAs and genes involved in neural tube defect (NTD) development.</p><p><bold>Methods:</bold> Pregnant Kunming mice (9–10 weeks of age) were treated with retinoic acid to construct NTD models. m<sup>6</sup>A levels and methyltransferase-like 3 (<italic>METTL3</italic>) expression were evaluated in brain tissues of the NTD models. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were performed on the NovaSeq platform and Illumina HiSeq 2,500 platform, respectively. Differentially m<sup>6</sup>A-methylated differentially expressed lncRNAs (DElncRNAs) and differentially expressed genes (DEGs) were identified, followed by GO biological process and KEGG pathway functional enrichment analyses. Expression levels of several DElncRNAs and DEGs were evaluated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) for validation.</p><p><bold>Results:</bold> m<sup>6</sup>A levels and <italic>METTL3</italic> expression levels were significantly lower in the brain tissues of the NTD mouse model than in controls. By integrating MeRIP-seq and RNA-seq data, 13 differentially m<sup>6</sup>A-methylated DElncRNAs and 170 differentially m<sup>6</sup>A-methylated DEGs were identified. They were significantly enriched in the Hippo signaling pathway and mannose-type <italic>O</italic>-glycan biosynthesis. The qRT-PCR results confirmed the decreased expression levels of lncRNAs, such as Mir100hg, Gm19265, Gm10544, and Malat1, and genes, such as <italic>Zfp236</italic>, <italic>Erc2</italic>, and <italic>Hmg20a</italic>, in the NTD group.</p><p><bold>Conclusion:</bold><italic>METTL3</italic>-mediated m<sup>6</sup>A modifications may be involved in NTD development. In particular, decreased expression levels of Mir100hg, Gm19265, Gm10544, Malat1, <italic>Zfp236</italic>, <italic>Erc2</italic>, and <italic>Hmg20a</italic> may contribute to the development of NTD.</p>