Next Article in Journal
Treatment-Emergent Cilgavimab Resistance Was Uncommon in Vaccinated Omicron BA.4/5 Outpatients
Previous Article in Journal
Pharmacotherapeutic Approaches to Treatment of Muscular Dystrophies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 13
Issue 10
10.3390/biom13101537
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Transcriptomic Profiling of m6A Modification in Age-Related Hearing Loss

by
Menglong Feng
1,2,3,4,
Xiaoqing Zhou
4,
Yaqin Hu
4,
Juhong Zhang
4,
Ting Yang
1,2,3,4,
Zhiji Chen
4 and
Wei Yuan
1,2,3,4,*
1
Chongqing Medical University, Chongqing 400016, China
2
Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400714, China
3
Chongqing School, University of Chinese Academy of Sciences, Chongqing 400714, China
4
Department of Otolaryngology & Head and Neck, Chongqing General Hospital, Chongqing 401147, China
*
Author to whom correspondence should be addressed.
Biomolecules 2023, 13(10), 1537; https://doi.org/10.3390/biom13101537
Submission received: 9 August 2023 / Revised: 27 September 2023 / Accepted: 7 October 2023 / Published: 18 October 2023
Browse Figures
Review Reports Versions Notes

Abstract

:
Age-related hearing loss (ARHL), also known as presbycusis, is one of the most common neurodegenerative disorders in elderly individuals and has a prevalence of approximately 70–80% among individuals aged 65 and older. As ARHL is an intricate and multifactorial disease, the exact pathogenesis of ARHL is not fully understood. There is evidence that transcriptional dysregulation mediated by epigenetic modifications is widespread in ARHL. However, the potential role of N6-methyladenosine (m6A) modification, as a crucial component of epigenetics, in ARHL progression remains unclear. In this study, we confirmed that the downregulation of m6A modification in cochlear tissues is related to ARHL and found that the expression of the m6A methylation regulators Wilms tumour suppressor-1-associated protein (WTAP), methyltransferase-like 3 (METTL3), ALKB homologous protein 5 (ALKBH5) and fat mass and obesity-associated protein (FTO) is decreased significantly at the mRNA and protein levels in ARHL mice. Then, we used methylated RNA immunoprecipitation sequencing (MeRIP-Seq) and RNA sequencing (RNA-Seq) to identify the differentially m6A-methylated genes in the cochlear tissues of ARHL mice. A total of 3438 genes with differential m6A methylation were identified, of which 1332 genes were m6A-hypermethylated and 2106 genes were m6A-hypomethylated in the ARHL group compared to the control group according to MeRIP-seq. Further joint analysis of RNA-Seq and MeRIP-Seq data showed that 262 genes had significant differences in both mRNA expression and m6A methylation. GO and KEGG analyses indicated that 262 unique genes were enriched mainly in the PI3K-AKT signalling pathway. In conclusion, the results of this study reveal differential m6A methylation patterns in the cochlear tissues of ARHL mice, providing a theoretical basis for further study of the pathogenesis of ARHL and potential therapeutic strategies.

1. Introduction

Age-related hearing loss (ARHL), one of the three most common chronic diseases among older people, has a prevalence of approximately 70–80% among individuals aged 65 and older and has gradually become a major global public health issue [1]. ARHL is caused by age-related degeneration of the cochlear and the central auditory system, and its main clinical manifestations are high-frequency hearing loss and decreased speech comprehension (especially in noisy environments) [2,3]. As the main cause of sensory disorders, ARHL limits the ability of elderly individuals to communicate with their families and society, leading to social isolation, anxiety, depression, and cognitive decline [4]. Although hearing aids and cochlear implant surgery have been shown to significantly improve hearing loss in patients with ARHL, the majority of patients are still not effectively treated due to the inconvenience of treatment or high medical expenses [5]. Hair cells and spiral ganglion neurons are the two main cell types responsible for hearing in the cochlea [6]. Current research has found that a reduction in cell numbers due to irreversible degeneration and death of inner hair cells and spiral ganglion neurons in the cochlea is the primary pathological change in ARHL [7]. In addition, oxidative stress, mitochondrial DNA deletions or mutations, reduced autophagy, disorders of ion transport across the membrane and damage to vascular endothelial cells have been found to be closely involved in the above pathological changes [8]. Nevertheless, the precise molecular mechanisms and signalling pathways underlying ARHL have not yet been fully elucidated.
RNA modifications, as the crucial components of epigenetics, have been found to play important regulatory roles in various physiological processes and disease progression [9]. N6-methyladenosine (m6A) methylation is widely acknowledged as the most common modification observed in eukaryotic mRNAs and exerts crucial influences on various aspects of mRNA biology, including RNA splicing, transport, translation, and decay [10]. The process of m6A methylation is characterised by a dynamic nature and reversibility and is regulated by three types of enzymes known as “writers”, “erasers”, and “readers” [11]. Writers are m6A methyltransferases, mainly including Wilms tumour suppressor-1-associated protein (WTAP), methyltransferase-like 14 (METTL14), and methyltransferase-like 3 (METTL3). Erasers, known as m6A demethylases, include only two enzymes: ALKB homologous protein 5 (ALKBH5) and fat mass and obesity-associated protein (FTO). Additionally, readers play an essential role in the biological function of m6A methylation by recognizing and binding m6A methylation sites to regulate mRNA metabolism. The main readers are insulin-like growth factor 2 mRNA binding proteins (IGF2BPs) and the YT521-B homology (YTH) domain family (YTHDC1/2 and YTHDF1–3) [12].
Increasing evidence indicates that m6A methylation and its regulators are related to multiple ageing-associated diseases. Guo et al. showed that METTL3 positively regulates JAK2 mRNA expression levels to promote angiogenesis and atherosclerosis by increasing m6A methylation of JAK2 mRNA, which is facilitated by the recognition and binding of IGF2BP1 to the m6A site in JAK2 mRNA [13]. Chen et al. discovered a decrease in the global m6A modification of mRNAs in PC12 cells induced by 6-OHDA and in the striatum of the Parkinson’s disease rat brain. Further research confirmed that the decrease in m6A level promoted the expression of N-methyl-d-aspartate receptor 1 and alleviated Ca2+ influx and oxidative stress, leading to dopaminergic neuron apoptosis [14]. Shafik et al. revealed that METTL3 is downregulated in the hippocampus in Alzheimer’s disease (AD) patients and that the levels of METTL3 and m6A modification of AD-related mRNAs are also reduced in 5XFAD mice [15], while upregulation of METTL3 can promote STUB1 gene expression through an m6A-IGF2BP1-dependent mechanism involving stabilisation of STUB1 mRNA and promotion of autophagic p-Tau clearance in Aβ1-42-treated cells, thus improving AD [16]. Moreover, a recent study found that there are decreased global m6A modification and METTL3 levels in the blood of ARHL patients and the cochleae of mice. Further studies have revealed that METTL3 reduces apoptosis in ARHL by regulating the m6A level of SIRT1 mRNA in an IGF2BP3-dependent manner [17]. However, m6A methylome profiles in ARHL have never been reported.
In this study, we mapped the trend in m6A level in the mouse cochlea with age and quantified the expression of METTL3, WTAP, FTO, and ALKBH5 genes in the cochleae of ARHL mice compared with control mice. In addition, we acquired, for the first time, the transcriptome-wide m6A profile of ARHL mice by using a combination of methylated RNA immunoprecipitation sequencing (MeRIP-Seq) and high-throughput RNA sequencing (RNA-Seq) and found considerable numbers of significantly differentially expressed m6A-modified transcripts and related pathways. Our study reveals the important roles of m6A-modified transcripts in the molecular mechanisms of ARHL and identifies potential therapeutic targets of ARHL.

2. Materials and Methods

2.1. Animals

All experiments were conducted in C57BL/6J male mice, and a total of 96 male C57BL/6J mice were included in our experiments. The mice were divided into five groups: six-week-old (6 w) (n = 39), 3-month-old (3 m) (n = 6), 6-month-old (6 m) (n = 6), 9-month-old (9 m) (n = 6) and 12-month-old (12 m) (n = 39) mice. Of these, 6 each of the 6 w, 3 m, 6 m, 9 m and 12 m mice were used for m6A quantification; 9 each of the 6 w and 12 m mice were used for quantitative real-time PCR (qRT–PCR) and western blotting; 15 each of the 6 w and 12 m mice were used for MeRIP-Seq and RNA-Seq; and 9 each of the 6 w and 12 m mice were used for MeRIP-qPCR. All procedures involving animal use and care were approved by the Chongqing Medical University Animal Welfare Committee.

2.2. Auditory Brainstem Response (ABR) Analysis

The mice were anaesthetised using 2% isoflurane for induction and 1.5% isoflurane for maintenance before each measurement. The mouse body temperature was maintained at 37–38 °C by keeping the anaesthetised mice on a heating pad. All tests were performed in a soundproof chamber. The acoustic signals of click and 4-, 8-, 16-, 24- or 32 kHz pure tone bursts were channelled into the mouse ear canals. The amplified brainstem responses were recorded from the scalps of the mice by Intelligent Hearing Systems ( Smart-EP software v3.30, Miami, FL, USA) through subdermal electrode needles inserted on the overlying and ventral vertex regions of the left and right bullae. The stimulus intensities were recorded in 5-dB SPL intervals down from a maximum stimulus intensity to identify the minimum stimulation level at which no waveform could be visualised, and the minimum stimulation level was determined as the ABR threshold.

2.3. RNA Extraction and qRT–PCR

Total RNA was extracted from mouse cochlear tissue using TRIzol Reagent (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. DNase I was used to digest DNA after RNA extraction. A260/A280 ratios were calculated to detect the quality of RNA using a NanoDrop® ND-1000. PrimeScript™ RT Master Mix (RR036A, Takara Biotechnology Co., Ltd., Dalian, China) was used to convert total RNA to cDNA following the manufacturer’s protocol. qRT–PCR was performed using a TB Green® Premix Ex Taq™ II (RR820A, Takara Biotechnology Co., Ltd., Dalian, China) and LightCycler 480 system (LightCycler480 software 1.5, Roche, Basel, Switzerland). The conditions of the qRT–PCR cycle included a predenaturation step at 95 °C for 30 s followed by 40 cycles of 95 °C for 5 s and 60 °C for 20 s. The mRNA primers used in our study are listed in Table A1.

2.4. m6A Quantification

The alterations in global m6A levels in total mRNA were detected with an EpiQuik m6A RNA Methylation Quantification Kit (Colorimetric) (Epigentek, Farmingdale, NY, USA) following the manufacturer’s protocol. Two hundred nanograms of poly-A-purified RNA was used for sample analysis. Signals of m6A were detected at 450 nm.

2.5. Western Blot Analysis

Total protein was extracted using RIPA lysis buffer containing 1% phosphatase inhibitors and 1% PMSF and quantified using a BCA protein assay kit (Beyotime, Shanghai, China). Protein samples (25 μg) were separated with 4–12% SDS–PAGE and transferred from the gel to PVDF membranes (Merck Millipore, Kenilworth, NJ, USA). Then, the PVDF membranes were blocked with QuickBlock™ western solution (Beyotime, Shanghai, China) for 10 min and incubated with METTL3, WTAP, ALKBH5, FTO, and GAPDH primary antibodies at 4 °C overnight (dilution 1:2000; Sanying, Wuhan, China). The membranes were subsequently incubated with the secondary antibody, goat anti-rabbit IgG (dilution 1:4000; Sanying, Wuhan, China), for 2 h at room temperature. SuperSignal™ West Femto Maximum Sensitivity Substrate (Thermo Fisher Scientific, Waltham, MA, USA) was used to develop the bands. ImageJ software (ImageJ software v1.6.0, NIH, Bethesda, MD, USA) was used to quantify western blot results.

2.6. MeRIP-Seq and Data Analysis

MeRIP-Seq was completed by Shanghai Cloud-Seq Biotech, and m6A immunoprecipitation was performed by using a m6A RNA Methylation Kit (GenSeq Inc., Shanghai, China). Fifteen mice were randomly selected from each of the 6 w and 12 m groups, and each group of mice was further randomly divided into three mixed biological samples (five in each sample). TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA) was used to extract total RNA from each sample according to the manufacturer’s instructions. A NanoDrop® ND-1000 and agarose gel electrophoresis were used to detect the quantity, integrity, and quality of total RNA from the above six samples. The RNA was fragmented into approximately 200-nt fragments by RNA Fragmentation Reagents (GenSeq Inc., Shanghai, China) as instructed by the manufacturer and reserved 1 μg of fragmented RNA as the input control. Briefly, the tubes with 1 μg/μL RNA and fragmentation buffer (metal-ion) were incubated at 70 °C for 6 min in a preheated thermal cycler block. Subsequently, 30 μL of 3M sodium acetate (NaOAc),1 μL of glycoBlue, and 750 μL volumes of 100% ethanol were added, mixed, and kept at −80 °C overnight. After centrifugation and pellet resuspension. Then, fragmented RNA was incubated with 5 μg of anti-m6A polyclonal antibody (Cat. No. 202003, Synaptic Systems, Goettingen, Germany) in IPP buffer for 2 h at 4 °C. The mixture was then immunoprecipitated by incubation with protein-A beads (Thermo Fisher Scientific, Waltham, MA, USA) at 4 °C for an additional 2 h. The RNA was eluted from the RNA-binding protein complexes with m6A in IPP buffer, and then the RNA was extracted with TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer’s instructions. Purified RNA was used for RNA-seq library generation with the GenSeq® Low Input Whole RNA Library Prep Kit (GenSeq Inc., Shanghai, China) and an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA) was used to assess the quality of the libraries. Both the input sample without immunoprecipitation and the m6A IP samples were subjected to 150 bp paired-end sequencing on an Illumina HiSeq 6000 sequencer (Illumina, San Diego, CA, USA).
After sequencing, the Q30 of the raw data was determined for assessment of sequencing quality. Cutadapt (v1.9.3) (https://journal.embnet.org/index.php/embnetjournal/article/view/200, accessed on 6 October 2023) software was used to remove adaptor sequences and low-quality reads. The clean reads of all samples were matched to the mouse reference genome (mm10) using HISAT2 (v2.0.4) software. MACS (v1.4.2) (http://liulab.dfci.harvard.edu/MACS/ accessed on 6 October 2023) software was used to identify methylated genes in each sample under the parameters: “(-g): 1.27 × 109, p-value < 0.00001”. Next, DREME (v5.3.0) (DREME Tutorial - MEME Suite (http://memesuite.org/, accessed on 6 October 2023) software was utilised to analyse the motifs of the m6A peaks with the criteria of E-value < 1.0 × 10−10. The diffReps (v1.55.6) (https://code.google.com/p/diffreps/under, accessed on 6 October 2023) software identified differentially methylated genes between 6 w mice and 12 m mice with the criteria of a p-value < 0.00001 and a fold change ≥ 2. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were performed on differentially methylated genes by using the KEGG (www.genome.jp/kegg, accessed on 6 October 2023) and GO (www.geneontology.org, accessed on 6 October 2023) databases.

2.7. RNA-Seq and Data Analysis

The same total RNA samples employed for MeRIP-Seq were used for RNA-Seq. RNA-Seq was also completed by Shanghai Cloud-Seq Biotech. A GenSeq® rRNA Removal Kit (GenSeq Inc., Shanghai, China) was used to remove the ribosomal RNA from RNA samples. After the removal of rRNA, the above samples were applied to construct an RNA sequencing library following the manufacturer’s instructions with a GenSeq® Low Input RNA Library Prep Kit (GenSeq Inc., Shanghai, China). Quality control and quantitation of the constructed RNA sequencing library were carried out with a Bioanalyzer 2100 system (Agilent Technologies, Inc., Santa Clara, CA, USA), and the libraries were subsequently sequenced (150 bp, paired-end reads) on an Illumina NovaSeq 6000 instrument (Illumina, San Diego, CA, USA).
After sequencing, the Q30 of the raw data was determined to assess sequencing quality. Cutadapt (v1.9.3) (https://journal.embnet.org/index.php/embnetjournal/article/view/200, accessed on 6 October 2023) software was used to remove adaptor sequences and low-quality reads. The clean reads of all samples were matched to the reference genome using HISAT2 (v2.0.4) (https://daehwankimlab.github.io/hisat2/, accessed on 6 October 2023) software. HTSeq (v0.9.1) (High-throughput sequence analysis in Python-HTSeq 2.0.4 documentation) software was used to calculate the raw counts. Then, edgeR was used for standardisation, and the fold change and p-value between the two groups of samples were calculated to screen for differentially expressed mRNAs with the criteria of a p-value < 0.05 and a fold change ≥ 2. KEGG and GO analyses of differentially expressed mRNAs were performed as described above.

2.8. MeRIP-qPCR

MeRIP-qPCR was performed by Shanghai Cloud-Seq Biotech. This experiment consisted of two parts: MeRIP and qRT–PCR. The experimental flow of MeRIP and qRT–PCR was performed as described previously. Briefly, the fragmented RNA was subjected to immunoprecipitation with 5 μg of anti-m6A polyclonal antibody (Cat. No. 202003, Synaptic Systems, Goettingen, Germany) coated magnetic beads. Methylated RNA bound to the magnetic beads was eluted and purified. Finally, the m6A modification was verified by qPCR analysis. The sequences of the primers used in MeRIP-qPCR are listed in Table A2. Percentage (IP/Input) was used to calculate the modification difference between the two groups.

3. Results

3.1. Increased Hearing Thresholds and Decreased Cochlear m6A Modification in Ageing Mice

The C57BL/6J mouse strain is known for early-onset hearing loss and is one of the most extensively used experimental models for studying ARHL. To explore the change in hearing with ageing, the hearing thresholds of click, 4 kHz, 8 kHz, 16 kHz and 32 kHz in 6-week-old (6 w), 3-month-old (3 m), 6-month-old (6 m), 9-month-old (9 m) and 12-month-old (12 m) C57BL/6J mice were measured by auditory brainstem response (ABR) assessment. We found that the 6 w mice showed normal mean hearing thresholds (<30 dB SPL), whereas the 12 m mice displayed severe hearing loss (>80 dB SPL) (Figure 1a–e). In addition, the hearing thresholds of mice showed a significantly upwards shift with ageing at all frequencies, and the hearing thresholds of 12 m mice were highest at all frequencies (Figure 1a–e). Previous reports have shown that C57BL/6J mice exhibit the classic pattern of ARHL by 12 to 15 months of age [18], which is similar to our above results. Therefore, 12-month-old C57BL/6J mice were selected as the experimental model of ARHL, and 6-week-old C57BL/6J mice were selected as the normal controls (NCs) in this study.
Moreover, to investigate whether ARHL is related to the alteration of m6A modification, the m6A levels in the cochleae of 6 w, 3 m, 6 m, 9 m, and 12 m mice were compared by colorimetric m6A quantification assay. As indicated in Figure 1f, the m6A levels in the cochleae of mice showed overall downward trends with ageing, and the total m6A levels in 9 m and 12 m mice were significantly lower than those in 6 w mice. These results suggest that the development of ARHL is accompanied by changes in the level of m6A in the cochlea.

3.2. Changes in m6A-Related Modification Enzymes in the Cochleae of 6 w and 12 m Mice

To explore the changes in m6A-related modification enzymes in ARHL, we used qRT–PCR and western blotting to measure the expression of genes associated with methylases and demethylases in the cochleae of 6 w and 12 m mice, including WTAP, METTL3, FTO, and ALKBH5 genes. As shown in Figure 2a–d, we found that the expression of WTAP, METTL3, FTO, and ALKBH5 genes was significantly downregulated in the 12 m mice compared with the 6 w mice at the mRNA level. In addition, we further measured the protein expression of the above four enzymes (Figure 3a–d). The levels of WTAP, METTL3, FTO, and ALKBH5 were also significantly downregulated in 12 m mice compared to 6 w mice, which was generally consistent with the pattern of changes at the mRNA level. Therefore, the combined effect of the downregulation of WTAP, METTL3, FTO, and ALKBH5 might be responsible for the decrease in total m6A in the cochlear RNA of ARHL mice.

3.3. General Characteristics of m6A Methylation Modification in ARHL

MeRIP-Seq was conducted to elucidate the characteristics of m6A methylation in 12 m mouse and 6 w mouse cochlear tissues. The bar graph shows the levels of m6A methylation in different groups. We found an average of 30,957 m6A peaks in 6-week-old mice and 28,147 m6A peaks in 12-month-old mice (Figure 4a). Subsequently, we mapped the above m6A peaks to 11,633 annotated genes in 6 w mouse cochlea tissues and 10,963 annotated genes in 12 m mouse cochlea tissues, of which 10,566 annotated genes were common between the two groups (Figure 4b). All m6A-methylated genes were split into groups based on the number of peaks in each gene (Figure 4c). Of note, most genes had two or more m6A peaks in the 6 w mice and 12 m mice. However, the number of peaks per gene in the 12 m mice was lower than that in the 6 w mice, which also demonstrated that the m6A modification level decreased after ARHL. In addition, our results revealed that the distribution of m6A peaks was different between the 6 w mice and 12 m mice, but both were enriched mainly in coding sequences (CDSs), stop codons (stopC), and start codons (startC) of mRNA (Figure 4d,e). More specifically, there were relative decreases in CDSs (6 w: 52.9%, 12 m: 47.5%) and stopC (6 w: 23.8%, 12 m: 16.5%), and relative increases in the numbers of m6A peaks in startC (6 w: 15.0%, 12 m: 16.7%), 3′ untranslated regions (3′UTR) (6 w: 1.9%, 12 m: 6.9%), and 5′ untranslated regions (5′UTR) (6 w: 6.4%, 12 m: 12.4%) (Figure 4e).
Considering that the binding of various binding proteins to specific motifs is the premise for the onset of RNA methylation and demethylation, we conducted motif enrichment analysis on the m6A peaks. As shown in Figure 4f, the top-ranking motif structure was characterised as DGAAGH (D = A/G/U; H = A/C/U) in both groups. Analysis of the m6A methylation enrichment in different chromosome loci showed that the m6A methylation peaks of chromosomes were increased on the X and Y chromosomes of the 12 m mice compared with the control mice but decreased in the remaining chromosomes (Figure 4g,h). Furthermore, in the 6 w and 12 m mice, the chromosomes with the highest m6A methylation frequency were chromosome 2, chromosome 11, and chromosome 1, with 3378 and 3268, 3173 and 3014, and 2997 and 2826 m6A methylation peaks, respectively (Figure 4h). However, further comparison showed that there was no significant difference in the enrichment number of m6A methylation peaks on chromosomes between the 12 m and 6 w groups.

3.4. Functional Enrichment and Pathway Analysis of Differentially m6A-Methylated Genes

We generated volcano plots and heatmaps to visually display the significant differentially methylated genes (12 m vs. 6 w). As shown in Figure 5a,b, a total of 3438 genes with differential m6A methylation were identified, of which 1332 genes were m6A-hypermethylated and 2106 genes were m6A-hypomethylated in the ARHL group compared to the control group (p-value < 0.00001 and |fold change| ≥ 2).
Furthermore, GO analysis was performed to gain further insights into the potential functions of the differentially methylated genes in 12 m mice. For the biological process (BP) category (Figure 6a,e), the genes with hypermethylated or hypomethylated m6A-modified peaks were all significantly enriched with metabolism-related terms, such as metabolic process, cellular metabolic process, and primary metabolic process. For the cellular component (CC) category (Figure 6b,f), the genes with hypermethylated or hypomethylated m6A-modified peaks were all significantly enriched in intracellular anatomical structure, organelle, intracellular organelle, etc. For the molecular function (MF) category (Figure 6c,g), the differentially methylated genes were associated with binding, protein binding, ion binding, etc.
KEGG analysis was further performed to explore the potential pathways associated with the differentially methylated genes. The genes with hypermethylated m6A-modified peaks were significantly enriched in several key pathways, including RNA degradation and the thyroid hormone signalling pathway (Figure 6d). However, the genes with hypomethylated m6A-modified peaks were significantly associated with several key pathways, including the MAPK signalling pathway, Ras signalling pathway, and Rap1 signalling pathway (Figure 6h).

3.5. Combined Analysis of MeRIP-Seq and RNA-Seq Data

The transcriptional profiles of the 6 w and 12 m mice were obtained from RNA-Seq. We found 1241 differentially expressed genes in the ARHL group compared with the control group, 290 upregulated and 951 downregulated (p-value < 0.05 and |fold change| ≥ 2) (Figure 7a,b). Further combined analysis of differentially expressed genes and differentially methylated m6A peaks identified a total of 262 genes that were significantly differentially expressed (Figure 7c). Of those, compared with the control group, 27 hypermethylated m6A peaks were distributed in 19 upregulated genes, 182 hypermethylated m6A peaks in 67 downregulated genes, 42 hypomethylated m6A peaks in 25 upregulated genes, and 276 hypomethylated m6A peaks in 137 downregulated genes (Figure 7c,d). Furthermore, there were both hypomethylated m6A peaks and hypomethylated m6A peaks in 1 upregulated gene and 13 downregulated genes (Figure 7c). To illustrate the role of m6A modification in ARHL, 262 genes with both differentially methylated m6A peaks and differential expression were selected for further GO and KEGG analyses. The top highly enriched GO terms are shown in Figure 7e–g. These included cellular metabolic process, metabolic process, regulation of cellular metabolic process (BP category), binding, protein binding, enzyme binding (MF category), and intracellular anatomical structure, organelle, and intracellular organelle (CC category). In addition, as shown in Figure 7h, the PI3K-AKT signalling pathway was highly enriched among the KEGG pathways of those genes.

3.6. Validation of Differentially m6A-Modified Genes by MeRIP-qPCR and qRT–PCR

To verify the reliability of the MeRIP-Seq and RNA-Seq results, MeRIP-qPCR and qRT–PCR were used to verify the modification levels of modification sites on the target genes and the gene expression levels, respectively. We selected the top two hypermethylated genes and the top two hypomethylated genes with the most significant differences in methylation levels as target genes from among 262 genes with significant differences in both methylated m6A peaks and mRNA expression levels. The MeRIP-qPCR and qRT–PCR results of RAPGEF6, BIRC6, RPS6KA3, and SH2D1B1 genes were in line with the sequencing results, and there were significant differences (Figure 8a–h). These results suggested that our sequencing results were accurate and credible. The information of the top 10 hypermethylated genes and top 10 hypomethylated genes among the 262 genes with differentially expressed mRNA and differentially methylated m6A peaks were listed in Table 1.

4. Discussion

ARHL, one of the most common sensory disorders, has become the main cause of communication deficits in elderly individuals and seriously affects their quality of life [19]. A growing number of studies have revealed that epigenetic mechanisms, including DNA methylation, histone modification, and noncoding RNA activities, are involved in the development of ARHL [20,21]. For example, Xu et al. have reported that DNA hypermethylation of the SLC26A4 gene results in an increased risk for ARHL in men [22], while DNA hypermethylation of the KCNQ5, ERBB3, and SOCS3 genes results in an increased risk for ARHL in women [23]. Wang et al. found that acetylated histone H3 can be detected in spiral ganglion cells and the organ of Corti in young mouse cochleae, while it is absent in aged mouse cochleae [20]. Conversely, dimethylated histone H3 can be detected in aged mice but not in young mice. In addition, Chen et al. indicated that upregulation of miRNA expression in pathways promoting apoptosis and inhibiting autophagy and downregulation of miRNA expression in pathways promoting proliferation and differentiation are involved in the progression of age-related hearing loss [24]. In recent years, one of the crucial components of epigenetics, m6A methylation, has been reported to be associated with a variety of ageing-associated degenerative diseases, such as chronic obstructive pulmonary disease, atherosclerosis, and age-related macular degeneration [13,25,26]. However, there are limited studies regarding m6A methylation in the pathological progression of ARHL.
To investigate whether m6A modification plays a role in the development of ARHL, we first investigated the association between m6A levels and ARHL. We found that ABR thresholds of experimental model mice increased with increasing age, but total RNA m6A modification in the cochleae of experimental model mice declined gradually. The level of m6A modification in the cochleae of 12 m mice was the lowest, and the ABR thresholds were greatest in these mice. In addition, Wang et al. reported similar findings that the levels of m6A modification in the total RNA of blood cells in ARHL normal patients were significantly lower than those in elderly patients without hearing loss [17]. These results suggest that the downregulation of m6A modification in the cochlea may be correlated with the development of ARHL. Considering that m6A modification is directly regulated by methyltransferases and demethylases [11], we measured the expression of genes encoding four common regulators in the cochleae of 12 m and 6 w mice, including WTAP, METTL3, FTO, and ALKBH5. Interestingly, at both the mRNA and protein levels, the expression of the methyltransferases WTAP and METTL3 genes was decreased in 12 m mice, and the expression of the demethylases FTO and ALKBH5 genes was also decreased significantly. This suggests that there may be a simultaneous occurrence of hypermethylation and hypomethylation of numerous genes in the cochlear tissues of ARHL mice and that the overall decrease in m6A level may be a result of coordinated regulation by multiple m6A enzymes.
To better clarify the changes in methylation modification and gene expression in cochlear tissues during ARHL, we constructed MeRIP-Seq and RNA-Seq libraries for high-throughput sequencing and performed a series of bioinformatics analyses on the results. We found that the transcripts of 11,633 genes had m6A modifications in 6 w mice and the transcripts of 10,963 genes in 12 m mice. In further analysis, 3438 differentially methylated genes between 6 w and 12 m mice, of which 1332 m6A hypermethylated genes and 2106 m6A hypomethylated genes were identified in 12 m mice. These results were consistent with our previous finding of a significant decrease in overall methylation level in 12 m mice, suggesting that our MeRIP-Seq results were reliable. Subsequent analysis revealed that although the degree of m6A methylation of genes was different between the two groups, the distribution and pattern of m6A methylation were similar, with major enrichment observed in the CDS region. Nevertheless, m6A deposition in the CDS region has been reported to tend to stabilize mRNA and be dynamic and reversible [27,28], which suggests that m6A modification may have a major role in regulating the expression of genes related to ARHL. Previous studies have revealed that the RRACH motif (R = A/G; H = A/C/U), the core sequence of m6A modification sites, is critical for m6A methylation [29]. The common m6A regulatory molecules WTAP, METTL14, and METTL3 have been shown to interact with the RNA motifs GACU, GGAC, and GGAC, respectively [30]. However, Gilbert et al. showed that not all RRACH sites in the body undergo m6A modification [31], which is consistent with our findings that the top-ranking motif structure in both groups was DGAAGH (D = A/G/U; H = A/C/U). This suggests that m6A modification in ARHL may be regulated by other molecular mechanisms and that further studies are needed. In addition, we also found many m6A motifs in 6 w and 12 m mouse tissues with some differences in sequence, further suggesting the possible presence of specific m6A methylation sites in ARHL.
In our study, GO and KEGG analyses of differentially m6A-modified genes were performed to explore their functions. Interestingly, in the GO analysis, we found that genes with significant hypomethylation or hypermethylation were mostly enriched with the same BP, MF, and CC terms, which suggests that these differentially methylated genes may have the same subcellular localisation and molecular biological functions. Of those, genes with significant hypomethylation or hypermethylation were mostly enriched in metabolic biological processes. Coling et al. reported that SOD1 gene deficiency increases the vulnerability of the cochlea to damage associated with ageing and noise through metabolic pathways involving the superoxide radical [32]. In addition, the overactivation of AMP-activated protein kinase leads to accelerated ARHL in Tg-mtTFB1 mice through redox imbalance and dysregulation of the apoptotic pathway, whereas a reduction in AMP-activated protein kinase can suppress cochlear cell apoptosis through the ROS-AMPK-BCL2 pathway [33]. These results are consistent with those of our GO analysis showing that metabolic enzymes and related regulatory factors may be involved in the occurrence of ARHL. Unlike the results obtained for GO analysis, the results of KEGG pathway analysis of hypomethylated genes were significantly different from those of hypermethylated genes. Among them, the genes with hypomethylated m6A-modified peaks were mainly associated with the MAPK signalling pathway. Previous evidence has shown that the MAPK signalling pathway plays an important role in regulating cell proliferation, apoptosis, autophagy, and differentiation [34], and the activated MAPK signalling pathway can promote the senescence of various cells and tissues [35,36]. However, inhibition of the MAPK signalling pathway has been confirmed to alleviate noise or drug-induced hair cell apoptosis and hearing loss [37,38], suggesting that ARHL might also be related to alteration of the MAPK signalling pathway. In addition, our results revealed that the genes with hypermethylated m6A-modified peaks were mainly enriched in RNA degradation. It is well known that RNA degradation is a major component of overall RNA metabolism and plays an important role in regulating levels of RNA expression [39]. However, RNA degradation is also an important way for m6A modification to exert biological functions [40], indicating that some genes with hypermethylated m6A-modified peaks might participate in the development of ARHL by altering RNA levels.
According to the joint analysis of MeRIP-Seq and RNA-Seq data, we found 262 genes with differential expression and differential m6A methylation simultaneously. Subsequently, the above differentially expressed genes were analysed via GO and KEGG analyses, and the main enriched biological processes were metabolic processes, while the main enriched pathway was the PI3K-AKT signalling pathway. Notably, in our study, most of the genes enriched in the PI3K-AKT signalling pathway are positive regulators of this pathway, including COL1A2, GHR, IL7R, ITGA6, LPAR1, and THBS4 genes. However, the expression of these genes was significantly downregulated in 12 m mice compared to the control group (Table A3). The PI3K-AKT signalling pathway serves an important role in regulating cell proliferation, differentiation, apoptosis, and senescence [41]. Previous studies have found that inhibition of the PI3K-AKT signalling pathway can induce cell senescence, and directly increase the expression of p21 [42,43]. Additionally, activation of this pathway has also been found to reduce hair cell loss and hearing damage in cisplatin-induced and noise-induced hearing loss [44,45]. These results suggested that the alteration of the PI3K-AKT signalling pathway may also play an important role in ARHL. In the following, compared with the control group, we further found a total of 86 genes were hypermethylated, 162 genes were hypomethylated, and 14 genes exhibited both hypermethylation and hypomethylation. The m6A was one of the important RNA modifications that affect RNA stability (half-life) [46]. In 162 hypomethylated genes, changes in the m6A levels of many genes have been found to affect their mRNA stability, such as APC gene (Table A3). Wang et al. found that METTL3 accelerated the degradation of APC mRNA in a YTHDF-dependent manner by up-regulating the m6A modification level of APC mRNA, thereby leading to oesophageal squamous cell carcinoma cell proliferation and tumour formation in mice [47]. Moreover, changes in the m6A levels of many of the above-hypermethylated genes have also been reported to affect their mRNA stability, such as SREBF1 gene (Table A3). Tang et al. demonstrated that upregulation of FTO levels might promote liver steatosis by increasing the stability of SREBF1 mRNA via demethylating its m6A sites [48]. However, whether the m6A methylation changes of these genes are related to the development of ARHL still needs further research in the future.
Among 262 genes with differential expression and differential m6A methylation simultaneously, some functional genes may be potential targets for future therapies, such as EGFR gene (Table A3). EGFR, a member of the epidermal growth factor receptor (HER) family, plays a key role in the regulation of cell proliferation and differentiation [49]. Previous studies have provided evidence that EGFR was required for the regenerative proliferation of auditory supporting cells in mammals and birds [50], and impaired glycolysis promotes alcohol exposure-induced apoptosis in HEI-OC1 cells by inhibiting EGFR signalling [51]. Interestingly, recent research has found that reader proteins IGF2BP3 or YTHDF2 also can directly bind to the m6A modification site of EGFR mRNA to change its stability [52,53]. However, the m6 A modification of EGFR mRNA in ARHL is still poorly understood and its role in ARHL is the next step in our research plan. High-throughput sequencing platforms may produce sequencing errors originating from the sequencing technology itself [54]. We validated differentially expressed genes in the combined MeRIP-Seq and mRNA-Seq analyses by MeRIP-qPCR and qRT-PCR, respectively, the results of which were consistent with those of high-throughput sequencing, further indicating the reliability of high-throughput sequencing in this study.

5. Conclusions

This study confirms that changes in m6A modification are strongly associated with ARHL and is the first to identify an m6A transcriptome pattern in ARHL mice. This study reveals the potential relationship between m6A methylation and mRNA expression in ARHL and provides new insights for further study of the pathogenesis of ARHL and potential therapeutic strategies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biom13101537/s1.

Author Contributions

All of the authors contributed to the conception of the article. Conceived of and designed the study: W.Y.; Performed the experiments and contributed materials: M.F., Y.H. and J.Z.; Analysed the data: W.Y., X.Z., T.Y. and Z.C.; Wrote the manuscript: M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was funded by Major Programs of Chongqing Science and Health Union (No2022DBXM006). National Natural Science Foundation of China (No. 81873702) and Natural Science Foundation of Chongqing (No.2022NSCQ-MSX2839, Nos. CSTB2022NSCQ-MSX0553, 2023jsqd180). Young and Middle-aged Medical Excellence Team of Chongqing. Postdoctoral Science Foundation Project of Chongqing Natural Science Foundation (No. cstc2021jcyj-bshx0026).

Institutional Review Board Statement

The animal study protocol was approved by the Ethics Committee of Chongqing Medical University (IACUC-CQMU-2022-0028; data of approval: 31 December 2022).

Informed Consent Statement

Not applicable.

Data Availability Statement

All of the data generated or analyzed during this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Age-related hearing loss, ARHL; Wilms tumour suppressor-1-associated protein, WTAP; methyltransferase-like 3, METTL3; ALKB homologous protein 5, ALKBH5; Fat mass and obesity associated protein, FTO; methylated RNA immunoprecipitation sequencing, MeRIP-Seq; RNA sequencing, RNA-Seq; quantitative real-time PCR, qRT–PCR; Auditory Brainstem Response, ABR; biological process, BP; cellular component, CC; molecular function, MF; week-old, w; month-old, m; the number of animals, n; the number of experiments, N.

Appendix A

Table A1. Primers used for qRT–PCR.
Table A1. Primers used for qRT–PCR.
GeneForward and Reverse Primer
WTAPF: 5′ GAAGGAGACACGACAGCAGTTGG 3′
R: 5′ GCTTGTGACCTCTGCCTGATCTAC 3′
METTL3F: 5′ CGCTGCCTCCGATGTTGATCTG 3′
R: 5′ CTGACTGACCTTCTTGCTCTGCTG 3′
FTOF: 5′ ATGAAGACGCTGTGCCACTGTG 3′
R: 5′ CACGTTGTAGGCTGCTCTGCTC 3′
ALKBH5F: 5′ GCAAGGTGAAGAGCGGCATCC 3′
R: 5′GTCCACCGTGTGCTCGTTGTAC 3′
RAPGEF6F: 5′ AGCTGTCAGGACCAAAGTCG 3′
R: 5′ CGCTCAGGTGTTGCCTTGAT 3′
BIRC6F: 5′ TGCAGAGGCTCAAAGTAGCC 3′
R: 5′ GGAACAGAATGCCCATTGGA 3′
RPS6KA3F: 5′ GACAGCGCGGAGAATGGACA 3′
R: 5′ CCGAACACGGTCTCGAACTT 3′
SH2D1B1F: 5′ TGTGCCTCTGTGTCTCGTTTA 3′
R: 5′ TTGTCCAATGTCCTTGTCTTCA 3′
β-ACTINF: 5′ AGTGTGACGTTGACATCCGT 3′
R: 5′ TGCTAGGAGCCAGAGCAGTA 3′
Table A2. Primers used for MeRIP-qPCR.
Table A2. Primers used for MeRIP-qPCR.
GeneForward and Reverse Primer
RAPGEF6F: 5′ CAGGGGATATGGAGCAGGC 3′
R: 5′ CAGGCCTCCGTCACTAAATAG 3′
BIRC6F: 5′ GCCTCCTCTCTCTTCCTTGG 3′
R: 5′ AATCAGCACCCCAGTAGTCA 3′
RPS6KA3F: 5′ GGCGTCCTCCTTTATACAATGC 3′
R: 5′ GTCCTTTGCTGTGTCTGAAAC 3′
SH2D1B1F: 5′ TGACCAAGCGAGAGTGTGAA 3′
R: 5′ GAGACACAGAGGCACAGGG 3′
Table A3. The information on some of the 262 genes.
Table A3. The information on some of the 262 genes.
IDGeneM6a Fold Changem6A
p Value		Gene
Log FC		Gene
p Value		Change
ENSMUSG00000029661COL1A2231.21.38075 × 10−10−1.60599020.04229598Hypo-down
ENSMUSG00000055737GHR11.74.29725 × 10−11−3.02140160.03698555Hypo-down
ENSMUSG00000003882IL7R5.18.2386 × 10−8−8.28040950.03681321Hypo-down
ENSMUSG00000027111ITGA62.91.19754 × 10−7−5.60422210.00835887Hyper-down
ENSMUSG00000038668LPAR14.44.55964 × 10−7−2.51374380.0377013Hypo-down
ENSMUSG00000021702THBS41044.72.7696 × 10−12−3.56060820.02895484Hypo-down
ENSMUSG00000005871APC7.71.8278 × 10−8−9.40207050.01765319Hypo-down
ENSMUSG00000020538SREBF12.32.8006 × 10−101.685438980.04796545Hyper-up
ENSMUSG00000020122EGFR6.32.51291 × 10−6−3.40247140.04096643Hyper-down

References

  1. De Iorio, M.L.; Rapport, L.J.; Wong, C.G.; Stach, B.A. Characteristics of Adults with Unrecognized Hearing Loss. Am. J. Audiol. 2019, 28, 384–390. [Google Scholar] [CrossRef]
  2. Makary, C.A.; Shin, J.; Kujawa, S.G.; Liberman, M.C.; Merchant, S.N. Age-Related Primary Cochlear Neuronal Degeneration in Human Temporal Bones. J. Assoc. Res. Otolaryngol. 2011, 12, 711–717. [Google Scholar] [CrossRef]
  3. Chavant, M.; Kapoula, Z. Presbycusis and the Aging of Eye Movement: Common Attention Mechanisms. Brain Sci. 2022, 12, 107. [Google Scholar] [CrossRef]
  4. Parham, K.; McKinnon, B.J.; Eibling, D.; Gates, G.A. Challenges and Opportunities in Presbycusis. Otolaryngol. Head Neck Surg. 2011, 144, 491–495. [Google Scholar] [CrossRef]
  5. Tu, N.C.; Friedman, R.A. Age-Related Hearing Loss: Unraveling the Pieces. Laryngoscope Investig. Otolaryngol. 2018, 3, 68–72. [Google Scholar] [CrossRef]
  6. Locher, H.; Frijns, J.H.; van Iperen, L.; de Groot, J.C.; Huisman, M.A.; de Sousa Lopes, S.M.C. Neurosensory Development and Cell Fate Determination in the Human Cochlea. Neural Dev. 2013, 8, 20. [Google Scholar] [CrossRef]
  7. Choi, J.E.; Ahn, J.; Moon, I.J. Associations Between Age-Related Hearing Loss and Dietaryassessment Using Data From Korean National Health Andnutrition Examination Survey. Nutrients 2021, 13, 1230. [Google Scholar] [CrossRef]
  8. Keithley, E.M. Pathology and Mechanisms of Cochlear Aging. J. Neurosci. Res. 2020, 98, 1674–1684. [Google Scholar] [CrossRef]
  9. Jiang, X.; Liu, B.; Nie, Z.; Duan, L.; Xiong, Q.; Jin, Z.; Yang, C.; Chen, Y. The Role of M6a Modification in the Biological Functions and Diseases. Signal Transduct. Target. Ther. 2021, 6, 74. [Google Scholar] [CrossRef]
  10. Zheng, Y.; Cong, S.; Ding, Y.; Yuan, B.; Wang, H.-Q.; Zhang, J.-B.; Zhang, W.-D.; Guo, H.-X. Comprehensive Analysis of Differences in N6-Methyladenosine RNA Methylomes in the Rat Adenohypophysis after Gnrh Treatment. FASEB J. 2022, 36, e22204. [Google Scholar]
  11. Liang, J.; Sun, J.; Zhang, W.; Wang, X.; Xu, Y.; Peng, Y.; Zhang, L.; Xiong, W.; Liu, Y.; Liu, H. Novel Insights Into the Roles of N(6)-Methyladenosine (M(6)a) Modification and Autophagy in Human Diseases. Int. J. Biol. Sci. 2023, 19, 705–720. [Google Scholar] [CrossRef]
  12. Gao, R.; Ye, M.; Liu, B.; Wei, M.; Ma, D.; Dong, K. M6a Modification: A Double-Edged Sword in Tumor Development. Front. Oncol. 2021, 11, 679367. [Google Scholar] [CrossRef]
  13. Dong, G.; Yu, J.; Shan, G.; Su, L.; Yu, N.; Yang, S. N6-Methyladenosine Methyltransferase METTL3 Promotes Angiogenesis and Atherosclerosis by Upregulating the JAK2/STAT3 Pathway Via M6a Reader IGF2BP1. Front. Cell Dev. Biol. 2021, 9, 731810. [Google Scholar] [CrossRef]
  14. Chen, X.; Yu, C.; Guo, M.; Zheng, X.; Ali, S.; Huang, H.; Zhang, L.; Wang, S.; Huang, Y.; Qie, S.; et al. Down-Regulation of M6a mRNA Methylation is Involved in Dopaminergic Neuronal Death. ACS Chem. Neurosci. 2019, 10, 2355–2363. [Google Scholar] [CrossRef]
  15. Shafik, A.M.; Zhang, F.; Guo, Z.; Dai, Q.; Pajdzik, K.; Li, Y.; Kang, Y.; Yao, B.; Wu, H.; He, C.; et al. N6-Methyladenosine Dynamics in Neurodevelopment and Aging, and its Potential Role in Alzheimer’s Disease. Genome Biol. 2021, 22, 17. [Google Scholar] [CrossRef]
  16. Tang, Z.; Cao, J.; Yao, J.; Fan, X.; Zhao, J.; Zhao, M.; Duan, Q.; Han, B.; Duan, S. KDM1A-Mediated Upregulation of METTL3 Ameliorates Alzheimer’s Disease Via Enhancing Autophagic Clearance of P-Tau through M6a-Dependent Regulation of STUB1. Free. Radic. Biol. Med. 2023, 195, 343–358. [Google Scholar] [CrossRef]
  17. Wang, E.; Li, Y.; Li, H.; Liu, Y.; Ming, R.; Wei, J.; Du, P.; Li, X.; Zong, S.; Xiao, H. METTL3 Reduces Oxidative Stress-Induced Apoptosis in Presbycusis by Regulating the N6-Methyladenosine Level of SIRT1 mRNA. Neuroscience 2023, 521, 110–122. [Google Scholar] [CrossRef]
  18. Teichert, M.; Liebmann, L.; Hübner, C.A.; Bolz, J. Homeostatic Plasticity and Synaptic Scaling in the Adult Mouse Auditory Cortex. Sci. Rep. 2017, 7, 17423. [Google Scholar] [CrossRef]
  19. Kasemsiri, P.; Yimtae, K.; Thanawirattananit, P.; Israsena, P.; Noymai, A.; Laohasiriwong, S.; Vatanasapt, P.; Siripaopradith, P.; Kingkaew, P. Effectiveness of a Programable Body-Worn Digital Hearing Aid for Older Adults in a Developing Country: A Randomized Controlled Trial with a Cross-Over Design. BMC Geriatr. 2021, 21, 437. [Google Scholar] [CrossRef]
  20. Wang, J.; Puel, J.-L. Presbycusis: An Update On Cochlear Mechanisms and Therapies. J. Clin. Med. 2020, 9, 218. [Google Scholar] [CrossRef]
  21. Li, Q.; Zang, Y.; Sun, Z.; Zhang, W.; Liu, H. Long Noncoding RNA GM44593 Attenuates Oxidative Stress From Age-Related Hearing Loss by Regulating Mir-29B/WNK1. Bioengineered 2022, 13, 573–582. [Google Scholar] [CrossRef] [PubMed]
  22. Xu, J.; Zheng, J.; Shen, W.; Ma, L.; Zhao, M.; Wang, X.; Tang, J.; Yan, J.; Wu, Z.; Zou, Z.; et al. Elevated SLC26A4 Gene Promoter Methylation is Associated with the Risk of Presbycusis in Men. Mol. Med. Rep. 2017, 16, 347–352. [Google Scholar] [CrossRef]
  23. Bouzid, A.; Smeti, I.; Dhouib, L.; Roche, M.; Achour, I.; Khalfallah, A.; Gibriel, A.A.; Charfeddine, I.; Ayadi, H.; Lachuer, J.; et al. Masmoudi. Down-Expression of P2RX2, KCNQ5, ERBB3 and SOCS3 through DNA Hypermethylation in Elderly Women with Presbycusis. Biomarkers 2018, 23, 347–356. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, H.H.R.; Wijesinghe, P.; Nunez, D.A. Micrornas in Acquired Sensorineural Hearing Loss. J. Laryngol. Otol. 2019, 133, 650–657. [Google Scholar] [CrossRef] [PubMed]
  25. Xia, H.; Wu, Y.; Zhao, J.; Cheng, C.; Lin, J.; Yang, Y.; Lu, L.; Xiang, Q.; Bian, T.; Liu, Q. N6-Methyladenosine-Modified Circsav1 Triggers Ferroptosis in Copd through Recruiting YTHDF1 to Facilitate the Translation of IREB2. Cell Death Differ. 2023, 30, 1293–1304. [Google Scholar] [CrossRef]
  26. Wang, Y.; Chen, Y.; Liang, J.; Jiang, M.; Zhang, T.; Wan, X.; Wu, J.; Li, X.; Chen, J.; Sun, J.; et al. METTL3-Mediated M6a Modification of HMGA2 mRNA Promotes Subretinal Fibrosis and Epithelial-Mesenchymal Transition. J. Mol. Cell Biol. 2023. [Google Scholar] [CrossRef]
  27. Zhou, L.; Tang, R.; Li, X.; Tian, S.; Li, B.; Qin, G. N(6)-Methyladenosine RNA Modification Regulates Strawberry Fruit Ripening in an ABA-Dependent Manner. Genome Biol. 2021, 22, 168. [Google Scholar] [CrossRef]
  28. Zheng, Y.; Li, Y.; Ran, X.; Wang, D.; Zheng, X.; Zhang, M.; Yu, B.; Sun, Y.; Wu, J. METTL14 Mediates the Inflammatory Response of Macrophages in Atherosclerosis through the NF-KAPPAB/Il-6 Signaling Pathway. Cell. Mol. Life Sci. 2022, 79, 311. [Google Scholar] [CrossRef]
  29. Niu, X.; Xu, J.; Liu, J.; Chen, L.; Qiao, X.; Zhong, M. Landscape of N(6)-Methyladenosine Modification Patterns in Human Ameloblastoma. Front. Oncol. 2020, 10, 556497. [Google Scholar] [CrossRef]
  30. Fan, C.; Ma, Y.; Chen, S.; Zhou, Q.; Jiang, H.; Zhang, J.; Wu, F. Comprehensive Analysis of the Transcriptome-Wide M6a Methylation Modification Difference in Liver Fibrosis Mice by High-Throughput M6a Sequencing. Front. Cell Dev. Biol. 2021, 9, 767051. [Google Scholar] [CrossRef]
  31. Gilbert, W.V.; Bell, T.A.; Schaening, C. Messenger RNA Modifications: Form, Distribution, and Function. Science 2016, 352, 1408–1412. [Google Scholar] [CrossRef] [PubMed]
  32. Coling, D.E.; Yu, K.C.; Somand, D.; Satar, B.; Bai, U.; Huang, T.-T.; Seidman, M.D.; Epstein, C.J.; Mhatre, A.N.; Lalwani, A.K. Effect of SOD1 Overexpression on Age- and Noise-Related Hearing Loss. Free. Radic. Biol. Med. 2003, 34, 873–880. [Google Scholar] [CrossRef] [PubMed]
  33. Zhao, J.; Li, G.; Zhao, X.; Lin, X.; Gao, Y.; Raimundo, N.; Li, G.-L.; Shang, W.; Wu, H.; Song, L. Down-Regulation of AMPK Signaling Pathway Rescues Hearing Loss in TFB1 Transgenic Mice and Delays Age-Related Hearing Loss. Aging 2020, 12, 5590–5611. [Google Scholar] [CrossRef] [PubMed]
  34. Moghbeli, M. Molecular Interactions of Mir-338 During Tumor Progression and Metastasis. Cell. Mol. Biol. Lett. 2021, 26, 13. [Google Scholar] [CrossRef]
  35. Sun, J.; Guo, Y.; Fan, Y.; Wang, Q.; Zhang, Q.; Lai, D. Decreased Expression of IDH1 by Chronic Unpredictable Stress Suppresses Proliferation and Accelerates Senescence of Granulosa Cells through Ros Activated MAPK Signaling Pathways. Free. Radic. Biol. Med. 2021, 169, 122–136. [Google Scholar] [CrossRef]
  36. Yang, H.; Yang, X.; Rong, K.; Liang, J.; Wang, Z.; Zhao, J.; Zhang, P.; Li, Y.; Wang, L.; Ma, H.; et al. Eupatilin Attenuates the Senescence of Nucleus Pulposus Cells and Mitigates Intervertebral Disc Degeneration Via Inhibition of the MAPK/NF-KAPPB Signaling Pathway. Front. Pharmacol. 2022, 13, 940475. [Google Scholar] [CrossRef]
  37. Liu, C.; Zheng, Z.; Li, W.; Tang, D.; Zhao, L.; He, Y.; Li, H. Inhibition of KDM5A Attenuates Cisplatin-Induced Hearing Loss Via Regulation of the MAPK/AKT Pathway. Cell. Mol. Life Sci. 2022, 79, 596. [Google Scholar] [CrossRef]
  38. Zhang, Y.; Huang, S.; Dai, X.; Xia, Z.-F.; Xiao, H.; He, X.-L.; Yang, R.; Li, J. SOD2 Alleviates Hearing Loss Induced by Noise and Kanamycin in Mitochondrial DNA4834-Deficient Rats by Regulating PI3K/MAPK Signaling. Curr. Med Sci. 2021, 41, 587–596. [Google Scholar] [CrossRef]
  39. Deutscher, M.P. Degradation of RNA in Bacteria: Comparison of mRNA and Stable RNA. Nucleic Acids Res. 2006, 34, 659–666. [Google Scholar] [CrossRef]
  40. Imam, H.; Kim, G.-W.; Mir, S.A.; Khan, M.; Siddiqui, A. Interferon-Stimulated Gene 20 (ISG20) Selectively Degrades N6-Methyladenosine Modified Hepatitis B Virus Transcripts. PLoS Pathog. 2020, 16, e1008338. [Google Scholar] [CrossRef]
  41. Weng, J.; Chen, M.; Shi, B.; Liu, D.; Weng, S.; Guo, R. Konjac Glucomannan Defends Against High-Fat Diet-Induced Atherosclerosis in Rabbits by Promoting the PI3K/AKT Pathway. Heliyon 2023, 9, e13682. [Google Scholar] [CrossRef]
  42. Zhang, Y.; Lv, Z.; Liu, Y.; Cao, H.; Yang, J.; Wang, B. PIN1 Protects Hair Cells and Auditory Hei-Oc1 Cells Against Senescence by Inhibiting the PI3K/AKT/MTOR Pathway. Oxidative Med. Cell. Longev. 2021, 2021, 9980444. [Google Scholar] [CrossRef] [PubMed]
  43. Courtois-Cox, S.; Williams, S.M.G.; Reczek, E.E.; Johnson, B.W.; McGillicuddy, L.T.; Johannessen, C.M.; Hollstein, P.E.; MacCollin, M.; Cichowski, K. A Negative Feedback Signaling Network Underlies Oncogene-Induced Senescence. Cancer Cell 2006, 10, 459–472. [Google Scholar] [CrossRef] [PubMed]
  44. Zha, D.-J.; Chen, F.-Q.; Fan, B.; Lu, F.; Du, W.-J.; Chen, J.; An, X.-G.; Wang, R.-F.; Li, W.; Song, Y.-L. PTEN Inhibitor Bisperoxovanadium Protects Against Noise-Induced Hearing Loss. Neural Regen. Res. 2023, 18, 1601–1606. [Google Scholar] [CrossRef] [PubMed]
  45. He, Y.; Zheng, Z.; Liu, C.; Li, W.; Zhao, L.; Nie, G.; Li, H. Inhibiting DNA Methylation Alleviates Cisplatin-Induced Hearing Loss by Decreasing Oxidative Stress-Induced Mitochondria-Dependent Apoptosis Via the LRP1-PI3K/AKT Pathway. Acta Pharm. Sin. B 2022, 12, 1305–1321. [Google Scholar] [CrossRef]
  46. Chen, J.; Fang, X.; Zhong, P.; Song, Z.; Hu, X. N6-Methyladenosine Modifications: Interactions with Novel RNA-Binding Proteins and Roles in Signal Transduction. RNA Biol. 2019, 16, 991–1000. [Google Scholar] [CrossRef]
  47. Wang, W.; Shao, F.; Yang, X.; Wang, J.; Zhu, R.; Yang, Y.; Zhao, G.; Guo, D.; Sun, Y.; Wang, J.; et al. METTL3 Promotes Tumour Development by Decreasing APC Expression Mediated by APC mRNA N(6)-Methyladenosine-Dependent YTHDF Binding. Nat. Commun. 2021, 12, 3803. [Google Scholar] [CrossRef]
  48. Tang, Z.; Sun, C.; Yan, Y.; Niu, Z.; Li, Y.; Xu, X.; Zhang, J.; Wu, Y.; Li, Y.; Wang, L.; et al. Aberrant Elevation of FTO Levels Promotes Liver Steatosis by Decreasing the M6a Methylation and Increasing the Stability of SREBF1 and CHREBP mRNAs. J. Mol. Cell Biol. 2023, 14, mjac061. [Google Scholar] [CrossRef]
  49. Dave, B.; Migliaccio, I.; Gutierrez, M.C.; Wu, M.-F.; Chamness, G.C.; Wong, H.; Narasanna, A.; Chakrabarty, A.; Hilsenbeck, S.G.; Huang, J.; et al. Loss of Phosphatase and Tensin Homolog or Phosphoinositol-3 Kinase Activation and Response to Trastuzumab or Lapatinib in Human Epidermal Growth Factor Receptor 2-Overexpressing Locally Advanced Breast Cancers. J. Clin. Oncol. 2011, 29, 166–173. [Google Scholar] [CrossRef]
  50. White, P.M.; Stone, J.S.; Groves, A.K.; Segil, N. EGFR Signaling is Required for Regenerative Proliferation in the Cochlea: Conservation in Birds and Mammals. Dev. Biol. 2012, 363, 191–200. [Google Scholar] [CrossRef]
  51. Kang, H.; Choi, S.J.; Park, K.H.; Lee, C.K.; Moon, J.S. Impaired Glycolysis Promotes Alcoholexposure-Induced Apoptosis in Hei-Oc1 Cells Via Inhibition of EGFR Signaling. Int. J. Mol. Sci. 2020, 21, 476. [Google Scholar] [CrossRef] [PubMed]
  52. Chen, L.-J.; Liu, H.-Y.; Xiao, Z.-Y.; Qiu, T.; Zhang, D.; Zhang, L.-J.; Han, F.-Y.; Chen, G.-J.; Xu, X.-M.; Zhu, J.-H.; et al. IGF2BP3 Promotes the Progression of Colorectal Cancer and Mediates Cetuximab Resistance by Stabilizing EGFR mRNA in an M(6)a-Dependent Manner. Cell Death Dis. 2023, 14, 581. [Google Scholar] [CrossRef] [PubMed]
  53. Zhong, L.; Liao, D.; Zhang, M.; Zeng, C.; Li, X.; Zhang, R.; Ma, H.; Kang, T. YTHDF2 Suppresses Cell Proliferation and Growth Via Destabilizing the EGFR mRNA in Hepatocellular Carcinoma. Cancer Lett. 2019, 442, 252–261. [Google Scholar] [CrossRef] [PubMed]
  54. Zhong, J.; Xie, Y.; Dang, Y.; Zhang, J.; Song, Y.; Lan, D. Use of RNA - Sequencing to Detect Abnormal Transcription of the Collagen Alpha-2 (Vi) Chain Gene that Can Lead to Bethlem Myopathy. Int. J. Mol. Med. 2021, 47, 28. [Google Scholar] [CrossRef] [PubMed]
Biomolecules 13 01537 g001
Figure 1. The change of ABR hearing thresholds and total m6A levels in C57BL/6J mice. (a–e) ABR thresholds were observed in 6 w (N = 39), 3 m (N = 6), 6 m (N = 6), 9 m (N = 6), and 12 m (N = 39) C57BL/6J mice at click, 4, 8, 16, and 32 kHz. (f) Quantification of m6A in total cochlear RNA was determined by colorimetric method in 6 w, 3 m, 6 m, 9 m, and 12 m C57BL/6J mice (N = 3). Abbreviations: w, week-old; m, month-old; ABR, Auditory brainstem response; N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference; ns represents no statistically significant difference. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 1. The change of ABR hearing thresholds and total m6A levels in C57BL/6J mice. (a–e) ABR thresholds were observed in 6 w (N = 39), 3 m (N = 6), 6 m (N = 6), 9 m (N = 6), and 12 m (N = 39) C57BL/6J mice at click, 4, 8, 16, and 32 kHz. (f) Quantification of m6A in total cochlear RNA was determined by colorimetric method in 6 w, 3 m, 6 m, 9 m, and 12 m C57BL/6J mice (N = 3). Abbreviations: w, week-old; m, month-old; ABR, Auditory brainstem response; N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference; ns represents no statistically significant difference. * p < 0.05; ** p < 0.01; *** p < 0.001.
Biomolecules 13 01537 g001
Biomolecules 13 01537 g002
Figure 2. The RNA expression of m6A enzymes in 6 w and 12 m C57BL/6J mice. (a) WTAP, (b) METTL3, (c) FTO, and (d) ALKBH5 in 6 w C57BL/6J mice compared with 12 m C57BL/6J mice (N = 3). Abbreviations: w, week-old; m, month-old; N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 2. The RNA expression of m6A enzymes in 6 w and 12 m C57BL/6J mice. (a) WTAP, (b) METTL3, (c) FTO, and (d) ALKBH5 in 6 w C57BL/6J mice compared with 12 m C57BL/6J mice (N = 3). Abbreviations: w, week-old; m, month-old; N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference. * p < 0.05; ** p < 0.01; *** p < 0.001.
Biomolecules 13 01537 g002
Biomolecules 13 01537 g003
Figure 3. The protein expression of m6A enzymes in 6 w and 12 m C57BL/6J mice. (a) WTAP, (b) METTL3, (c) FTO, and (d) ALKBH5 in 6 w C57BL/6J mice compared with 12 m C57BL/6J mice (N = 3). Abbreviations: w, week-old; m, month-old; N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference. * p < 0.05; ** p < 0.01. Western Blot images can be found in Supplementary Materials.
Figure 3. The protein expression of m6A enzymes in 6 w and 12 m C57BL/6J mice. (a) WTAP, (b) METTL3, (c) FTO, and (d) ALKBH5 in 6 w C57BL/6J mice compared with 12 m C57BL/6J mice (N = 3). Abbreviations: w, week-old; m, month-old; N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference. * p < 0.05; ** p < 0.01. Western Blot images can be found in Supplementary Materials.
Biomolecules 13 01537 g003
Biomolecules 13 01537 g004
Figure 4. Overview of m6A methylation patterns in 6 w mice and 12 m mice. (a) The average number of m6A peaks in the 6 w mice and the 12 m mice. (b) Venn diagram of m6A-modified genes in the 6 w mice and the 12 m mice. (c) Number of peaks per gene. (d) Density of differential m6A peaks along transcripts. Each transcript was divided into five parts: 5′UTR, CDS, 3′UTR, startC, and stopC. (e) Pie charts showing the region of m6A peaks in each group. (f) The most conserved sequence motif of the differential m6A peak region. (g) The distribution patterns of m6A peaks in different chromosomes. (h) The count of m6A peaks in per chromosome. Abbreviations: w, week-old; m, month-old.
Figure 4. Overview of m6A methylation patterns in 6 w mice and 12 m mice. (a) The average number of m6A peaks in the 6 w mice and the 12 m mice. (b) Venn diagram of m6A-modified genes in the 6 w mice and the 12 m mice. (c) Number of peaks per gene. (d) Density of differential m6A peaks along transcripts. Each transcript was divided into five parts: 5′UTR, CDS, 3′UTR, startC, and stopC. (e) Pie charts showing the region of m6A peaks in each group. (f) The most conserved sequence motif of the differential m6A peak region. (g) The distribution patterns of m6A peaks in different chromosomes. (h) The count of m6A peaks in per chromosome. Abbreviations: w, week-old; m, month-old.
Biomolecules 13 01537 g004
Biomolecules 13 01537 g005
Figure 5. (a) Volcano plot representation of microarray data on the differentially expressed m6A methylation genes. (b) Hierarchical cluster analysis of differentially m6A methylation genes.
Figure 5. (a) Volcano plot representation of microarray data on the differentially expressed m6A methylation genes. (b) Hierarchical cluster analysis of differentially m6A methylation genes.
Biomolecules 13 01537 g005
Biomolecules 13 01537 g006
Figure 6. Gene ontology enrichment and pathway analysis of differentially methylated genes. (a–c) The top 10 GO terms significantly enriched for hypermethylated genes. (e–g) The top 10 GO terms significantly enriched for hypomethylated genes. (d) The top 10 KEGG pathways significantly enriched for hypermethylated genes. (h) The top 10 KEGG pathways significantly enriched for hypomethylated genes.
Figure 6. Gene ontology enrichment and pathway analysis of differentially methylated genes. (a–c) The top 10 GO terms significantly enriched for hypermethylated genes. (e–g) The top 10 GO terms significantly enriched for hypomethylated genes. (d) The top 10 KEGG pathways significantly enriched for hypermethylated genes. (h) The top 10 KEGG pathways significantly enriched for hypomethylated genes.
Biomolecules 13 01537 g006
Biomolecules 13 01537 g007
Figure 7. Conjoint analysis of m6A methylation and mRNA expression. (a) The volcano plot showing the differentially expressed mRNA gene. (b) Hierarchical cluster analysis of differentially expressed mRNA genes. (c) Venn diagram of genes with differentially m6A methylation and differentially expressed mRNA. (d) Four quadrant graphs of peaks with differentially m6A methylation and differentially expressed mRNA. (e–g) The top 10 GO terms significantly enriched for genes with differentially m6A methylation and differentially expressed mRNA. (h) The top 10 KEGG pathways significantly enriched for genes with differentially m6A methylation and differentially expressed mRNA. Notes: the purple color circles represent the upregulated genes with hypermethylation, the red color circles represent the downregulated genes with hypermethylation, the blue color circles represent the downregulated genes with hypomethylation, the green color circles represent the upregulated genes with hypomethylation.
Figure 7. Conjoint analysis of m6A methylation and mRNA expression. (a) The volcano plot showing the differentially expressed mRNA gene. (b) Hierarchical cluster analysis of differentially expressed mRNA genes. (c) Venn diagram of genes with differentially m6A methylation and differentially expressed mRNA. (d) Four quadrant graphs of peaks with differentially m6A methylation and differentially expressed mRNA. (e–g) The top 10 GO terms significantly enriched for genes with differentially m6A methylation and differentially expressed mRNA. (h) The top 10 KEGG pathways significantly enriched for genes with differentially m6A methylation and differentially expressed mRNA. Notes: the purple color circles represent the upregulated genes with hypermethylation, the red color circles represent the downregulated genes with hypermethylation, the blue color circles represent the downregulated genes with hypomethylation, the green color circles represent the upregulated genes with hypomethylation.
Biomolecules 13 01537 g007
Biomolecules 13 01537 g008
Figure 8. Verification of differentially m6A methylation level and differentially mRNA expression. (a–d) MeRIP-qPCR showing m6A enrichment in mRNA of RAPGEF6, BIRC6, RPS6KA3, and SH2D1B1 genes (N = 3). (e–h) qRT-PCR showing mRNA expression of RAPGEF6, BIRC6, RPS6KA3, and SH2D1B1 genes (N = 3). Abbreviations: N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference.** p < 0.01; *** p < 0.001.
Figure 8. Verification of differentially m6A methylation level and differentially mRNA expression. (a–d) MeRIP-qPCR showing m6A enrichment in mRNA of RAPGEF6, BIRC6, RPS6KA3, and SH2D1B1 genes (N = 3). (e–h) qRT-PCR showing mRNA expression of RAPGEF6, BIRC6, RPS6KA3, and SH2D1B1 genes (N = 3). Abbreviations: N, the number of experiments. Notes: p < 0.05 represents a statistically significant difference.** p < 0.01; *** p < 0.001.
Biomolecules 13 01537 g008
Table 1. The top 10 hypermethylated and hypomethylated genes among the 262 genes.
Table 1. The top 10 hypermethylated and hypomethylated genes among the 262 genes.
IDGenem6A Fold Changem6A
p-Value		Gene
Log FC		Gene
p-Value		Change
ENSMUSG00000037533RAPGEF6433.21.2812 × 10−14−10.0253740.00213148Hyper-down
ENSMUSG00000024073BIRC6800.75.0835 × 10−14−2.99438060.02636128Hyper-down
ENSMUSG00000058997VWA8201.17.6235 × 10−14−3.70584470.01945395Hyper-down
ENSMUSG00000068036MLLT4353.61.4383 × 10−13−1.94149220.04990249Hyper-down
ENSMUSG00000022607PTK219.91.6866 × 10−13−3.46520580.01968117Hyper-down
ENSMUSG00000032410XRN11434.82.9063 × 10−13−10.0238740.00555607Hyper-down
ENSMUSG00000003847NFAT58.63.6304 × 10−13−2.35166770.03805174Hyper-down
ENSMUSG00000032555TOPBP111.88.4203 × 10−13−7.29464740.01763448Hyper-down
ENSMUSG00000024672MS4A76.72.8421 × 10−121.833057830.04776993Hyper- up
ENSMUSG00000029186PI4K2B3.52.8892 × 10−12−10.1258560.02227473Hyper-down
ENSMUSG00000067149IGJ312.91.4276 × 10−138.490156020.00112192Hypo-up
ENSMUSG00000031309RPS6KA324.951.8979 × 10−13−5.24279370.00759645Hypo-down
ENSMUSG00000096334SH2D1B1412.23.1109 × 10−136.187696860.02234852Hypo-up
ENSMUSG00000030231PLEKHA5167.23.5733 × 10−13−6.92625120.00946439Hypo-down
ENSMUSG00000024900CPT1A524.43.6441 × 10−13−2.5806430.01844615Hypo-down
ENSMUSG00000025278FLNB93.243.6648 × 10−13−3.70842020.03199615Hypo-down
ENSMUSG00000030811FBXL19470.73.9417 × 10−13−9.49879690.03882934Hypo-down
ENSMUSG00000024749TMC1543.75.0403 × 10−13−8.47348060.03470489Hypo-down
ENSMUSG00000061080LSAMP554.15.7905 × 10−139.72095040.00629816Hypo-up
ENSMUSG00000022297FZD6698.66.9979 × 10−13−9.892310.0045972Hypo-down
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Feng, M.; Zhou, X.; Hu, Y.; Zhang, J.; Yang, T.; Chen, Z.; Yuan, W. Comprehensive Transcriptomic Profiling of m6A Modification in Age-Related Hearing Loss. Biomolecules 2023, 13, 1537. https://doi.org/10.3390/biom13101537

AMA Style

Feng M, Zhou X, Hu Y, Zhang J, Yang T, Chen Z, Yuan W. Comprehensive Transcriptomic Profiling of m6A Modification in Age-Related Hearing Loss. Biomolecules. 2023; 13(10):1537. https://doi.org/10.3390/biom13101537

Chicago/Turabian Style

Feng, Menglong, Xiaoqing Zhou, Yaqin Hu, Juhong Zhang, Ting Yang, Zhiji Chen, and Wei Yuan. 2023. "Comprehensive Transcriptomic Profiling of m6A Modification in Age-Related Hearing Loss" Biomolecules 13, no. 10: 1537. https://doi.org/10.3390/biom13101537

APA Style

Feng, M., Zhou, X., Hu, Y., Zhang, J., Yang, T., Chen, Z., & Yuan, W. (2023). Comprehensive Transcriptomic Profiling of m6A Modification in Age-Related Hearing Loss. Biomolecules, 13(10), 1537. https://doi.org/10.3390/biom13101537

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop