Identifying Key MA-methylated LncRNAs and Genes Associated with Neural Tube Defects Integrative MeRIP and RNA Sequencing Analyses

N-methyladenosine (mA) is a common post-transcriptional modification of messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs). However, mA-modified lncRNAs are still largely unexplored. This study aimed to investigate differentially mA-modified lncRNAs and genes involved in neural tube defect (NTD) development. Pregnant Kunming mice (9-10 weeks of age) were treated with retinoic acid to construct NTD models. mA levels and methyltransferase-like 3 () 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 mA-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. mA levels and 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 mA-methylated DElncRNAs and 170 differentially mA-methylated DEGs were identified. They were significantly enriched in the Hippo signaling pathway and mannose-type -glycan biosynthesis. The qRT-PCR results confirmed the decreased expression levels of lncRNAs, such as Mir100hg, Gm19265, Gm10544, and Malat1, and genes, such as , , and , in the NTD group. -mediated mA modifications may be involved in NTD development. In particular, decreased expression levels of Mir100hg, Gm19265, Gm10544, Malat1, , , and may contribute to the development of NTD.