Next Article in Journal
DNA Hypomethylation May Contribute to Metabolic Recovery of Frozen Wood Frog Brains
Next Article in Special Issue
Epigenomic Approaches for the Diagnosis of Rare Diseases
Previous Article in Journal
The Role of Genetics, Epigenetics, and the Environment in ASD: A Mini Review
Previous Article in Special Issue
Epigenetic Signatures of Centrosomes Are Novel Targets in Cancer Diagnosis: Insights from an Analysis of the Cancer Genome Atlas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Global m6A RNA Methylation in SARS-CoV-2 Positive Nasopharyngeal Samples in a Mexican Population: A First Approximation Study

by
Jorge Luis Batista-Roche
1,
Bruno Gómez-Gil
2,
Gertrud Lund
3,
César Alejandro Berlanga-Robles
4 and
Alejandra García-Gasca
1,*
1
Molecular and Cellular Biology, Centro de Investigación en Alimentación y Desarrollo, Mazatlán 82112, Mexico
2
Microbial Genomics, Centro de Investigación en Alimentación y Desarrollo, Mazatlán 82112, Mexico
3
Department of Genetic Engineering, CINVESTAV Irapuato Unit, Irapuato 36821, Mexico
4
Environmental Management, Centro de Investigación en Alimentación y Desarrollo, Mazatlán 82112, Mexico
*
Author to whom correspondence should be addressed.
Epigenomes 2022, 6(3), 16; https://doi.org/10.3390/epigenomes6030016
Submission received: 28 April 2022 / Revised: 11 June 2022 / Accepted: 20 June 2022 / Published: 29 June 2022

Abstract

:
The Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2) is the causal agent of COVID-19 (Coronavirus Disease-19). Both mutation and/or recombination events in the SARS-CoV-2 genome have resulted in variants that differ in transmissibility and severity. Furthermore, RNA methylation of the N6 position of adenosine (m6A) is known to be altered in cells infected with SARS-CoV-2. However, it is not known whether this epitranscriptomic modification differs across individuals dependent on the presence of infection with distinct SARS-CoV-2 variants, the viral load, or the vaccination status. To address this issue, we selected RNAs (n = 60) from SARS-CoV-2 sequenced nasopharyngeal samples (n = 404) of 30- to 60-year-old outpatients or hospitalized individuals from the city of Mazatlán (Mexico) between February 2021 and March 2022. Control samples were non-infected individuals (n = 10). SARS-CoV-2 was determined with real-time PCR, viral variants were determined with sequencing, and global m6A levels were determined by using a competitive immunoassay method. We identified variants of concern (VOC; alpha, gamma, delta, omicron), the variant of interest (VOI; epsilon), and the lineage B.1.1.519. Global m6A methylation differed significantly across viral variants (p = 3.2 × 10−7). In particular, we found that m6A levels were significantly lower in the VOC delta- and omicron-positive individuals compared to non-infected individuals (p = 2.541236 × 10−2 and 1.134411 × 10−4, respectively). However, we uncovered no significant correlation between global m6A levels and viral nucleocapsid (N) gene expression or age. Furthermore, individuals with complete vaccination schemes showed significantly lower m6A levels than unvaccinated individuals (p = 2.6 × 10−4), and differences in methylation levels across variants in unvaccinated individuals were significant (p = 3.068 × 10−3). These preliminary results suggest that SARS-CoV-2 variants show differences in global m6A levels.

1. Introduction

SARS-CoV-2 is an enveloped, spheroidal-shaped virus that causes the disease known as COVID-19 in humans. It has a capsid that contains the structural proteins spike (S), membrane (M), envelope (E), and the nucleocapsid protein (N). Its genome is a single, positive stranded ribonucleic acid (RNA) that behaves as messenger RNA (mRNA) [1]. The N protein is a crucial structural component of SARS-CoV-2 that participates in the virion assembly and is often used as a diagnostic marker of viral infection [2]. In a vaccine setting, the N protein induces SARS-CoV-2-specific T cell proliferation and cytotoxic activity and promotes long-lasting T cell immunity [3].
The main receptor and route of entry of SARS-CoV-2 into the host cells is the angiotensin-converting enzyme 2 (ACE2), which binds with the S protein of SARS-CoV-2 and mediates the fusion of viral and cell membranes. ACE2 is found in the membranes of different human cells such as pneumocytes, enterocytes, proximal tubular cells of the kidney, and endothelial cells of veins and arteries [4]. Various mutations in the viral genome have given rise to SARS-CoV-2 variants, which differ not only in nucleotide sequences but also in transmissibility and severity [5,6].
Following the discovery of RNA methylation at N6- position of adenosine (m6A), the most prevalent modification of RNA in mammalian cells [7,8], this modification has also been detected in several viral genomes [9]. Most m6A is found in the consensus sequence motif DRACH (where A* denotes the methylated adenosine, D denotes A, G or U, R denotes A and G, and H denotes A, C or U) and mRNAs may contain from one to up to 20 m6A sites or more [10,11]. Genomic and viral levels of m6A are dynamically regulated by the activity of methyltransferase enzymes (i.e., “writers” such as METTL3, METTL14 that transfer methyl groups to nitrogen 6 of the adenosine) and demethylases (“erasers” such as fat mass and obesity-associated protein (FTO) and a ketoglutarate dependent dioxygenase homolog 5 (ALKBH5)). Once methylated, m6A is recognized by “reader” proteins (e.g., YTH N6-Methyladenosine RNA Binding Protein domain-containing proteins YTHDC1 and C2, YTHDF1, F2 and F3) that influence the fate of m6A RNA, either by guiding mRNAs to ribosomes for translation or to other complexes for their degradation [11,12].
It was recently reported that SARS-CoV-2 RNA, as well as negative-stranded viral RNAs, are dynamically methylated in vitro, as is the host cell RNA [13]. In human hepatocarcinoma cells (Huh7), SARS-CoV-2 genomic RNA is gradually methylated during infection and m6A occurs more frequently towards the 3′ region of the genome [14]. In addition, SARS-CoV-2 proteins interact with the host m6A methylation machinery to modulate viral replication. For example, an increase in METTL3 expression 48 h post-infection in Vero E6 cells was positively associated with SARS-CoV-2 replication [15]. In accordance, depletion of METTL3 leads to a reduction in SARS-CoV-2 replication [13,15,16], although the opposite has also been reported [14]. Furthermore, knockdown of specific erasers (ALKBH5) and readers (YTHDF1, YTHDF2, YTHDF3) have also been shown to impact on viral replication [14,16].
Recent studies show that SARS-CoV-2 infection changes the host cell m6A methylome in vitro, promoting differential expression of host genes [13], and in vivo, altering m6A modification levels in lymphocytes from peripheral blood samples by increased expression of the m6A methyltransferase RNA-binding motif protein 15 (RBM15) [17] or decreased expression of METTL3 in epithelial cells of bronchoalveolar lavage fluid of COVID-19 patients [13].
To date, only a few studies have probed for alterations in global m6A levels in a limited number of COVID-19 positive individuals (n = 2–20) [13,17,18]. These studies found increased levels of m6A in peripheral blood samples [17] or in epithelial cells of bronchoalveolar lavage fluid [13] in infected individuals relative to non-infected individuals that are associated with increased expression of RBM15, which encodes for a protein that facilitates m6A by guiding METTL3 to the target RNA sequence [17]. To our knowledge, no studies have compared m6A levels across individuals infected with distinct COVID-19 variants. Therefore, we measured global levels of this modification in nasopharyngeal samples from adult (30- to 60-year-old) male and female patients infected with different SARS-CoV-2 variants between February 2021 and March 2022.

2. Results

2.1. Viral Variants Detected in Mazatlán, Mexico during the COVID-19 Pandemic (February 2021–March 2022)

Analysis of nasopharyngeal samples from a total of 404 individuals (219 males and 185 females) revealed the presence of distinct SARS-CoV-2 variants during the second, third, and fourth waves of the COVID-19 pandemic in Mazatlán, Mexico from February 2021 to March 2022 (Figure 1). During the second wave, epsilon and B.1.1.519 lineage variants predominated in the population, peaking between February and March 2021. During the third wave, peaking between June and August, the delta variant dominated, although the alpha and gamma variants were also detected. Finally, the fourth wave, peaking between January and February, was dominated by the omicron variant.

2.2. Global m6A Levels and Nucleocapsid (N) Gene Expression in Nasopharyngeal Samples

Global m6A levels were determined with 60 randomly selected sequenced samples of each of the following variants: omicron (BA.1), delta (B.1.617.2), gamma (P.1), alpha (B.1.1.7), and epsilon (B.1.429), as well as the B.1.1.519 lineage (Figure 2A). Table 1 shows the average age, sex distribution, and vaccination status of individuals in each variant group. Vaccination data were available for 64 of the 70 samples analyzed; 8 and 4 samples were from patients with a complete or partial vaccination schemes, respectively, while 52 samples came from unvaccinated patients.
A Kruskal–Wallis test revealed significant differences in global m6A according to viral variant (p = 3.2 × 10−7). In particular, delta and omicron-positive patients showed significantly lower m6A levels than the control group (SARS-CoV-2 negative) (p = 2.541236 × 10−2 and 1.134411 × 10−4 for delta and omicron respectively).
To understand whether the differences in variant m6A levels were related to viral load, we compared expression levels of the viral nucleocapsid (N) gene across variants. Only delta and B.1.1.519 variants showed significant differences in N expression levels (Kruskal-Wallis, p = 4.1 × 10−2; Dunn, p = 2.4 × 10−2) (Figure 2B). However, we found no significant correlation, neither between m6A levels and relative expression of the N gene of all, or only unvaccinated individuals (r = −0.06, p = 0.63 and r = 0.004, p = 0.9803, respectively) (Figure 3), nor between m6A and age (r = 0.15, p = 0.21).
Finally, we evaluated whether patient vaccination status affected global m6A RNA methylation levels. Samples from patients with complete vaccination schemes showed significantly lower levels of RNA methylation compared to unvaccinated patients (Kruskal-Wallis, p = 2.6 × 10−4; Dunn test, p = 2.5 × 10−4) (Figure 4A). Importantly, samples from vaccinated patients were obtained during the third and fourth waves, dominated by delta and omicron variants, respectively. To exclude that vaccination status affected the observed differences between variants, we performed a Kruskal-Wallis analysis restricted to unvaccinated individuals. As previously observed, m6A levels differed significantly across variants (p = 3.068 × 10−3), with omicron and some gamma-infected patients showing the lowest methylation levels (Figure 4B). The delta variant was not represented in unvaccinated individuals.

3. Discussion

RNA modifications such as m6A play an important role in different physiological processes and human diseases. In the innate immune response, m6A-modfied mRNAs are essential for the translation of co-stimulatory molecules that allow dendritic cell maturation and subsequent T cell activation [19]. Also, the immunosuppressive functions of regulatory T cells are maintained by m6A-driven degradation of the suppressor of cytokine signaling (SOCS) transcripts [8]. In addition, RNAs from several viral families acquire m6A modifications as a common strategy to evade innate host immunity [20].
It has been reported that m6A deficiency in the vesicular stomatitis virus (VSV) genome triggers significantly higher levels of type I interferon [20]. However, some m6A sites of viral RNAs favor the evasion of the pattern recognition receptors of the cellular innate immune response and promote viral replication [21]. Conversely, certain transcripts of some RNA viruses and adenoviruses contain m6A-methylated sites that negatively influence infection [14]. Currently, there is a controversy in the literature regarding the effects of m6A methyltransferases on SARS-CoV-2 replication. Apparently, depending on the cell type, depletion of m6A methyltransferases can decrease viral load [13] or increase viral replication and the percentage of SARS-CoV-2 infected cells [14].
To our knowledge, this is the first report showing that global m6A levels of nasopharyngeal RNA samples of patients infected with SARS-CoV-2 differ among viral variants. Interestingly, the two most contagious variants (delta and omicron)—with omicron being the most contagious globally—[22] showed the lowest methylation levels. These differences are unlikely to reflect variation in disease severity across variants, as COVID-19 patients with beta and delta variants generally are at higher risk of developing severe disease compared to patients with alpha, gamma [23], and omicron variant(s), whose mutations have been suggested to contribute to the host immune escape [24]. Interestingly, DRACH motifs are highly conserved among SARS-CoV-2 variants [25], thus the contribution of DRACH motifs to differential methylation levels among variants remains unclear and require further research.
From a public health perspective, identifying a reliable marker of severe COVID-19 risk is of utmost importance. Unfortunately, we have no information regarding disease severity of the individuals analyzed in this study since most samples were from outpatients. In addition, we found no correlation between viral load and m6A levels across variants. However, previous studies have reported a higher m6A level in lung or peripheral blood from patients with moderate or severe COVID-19 disease compared to healthy individuals [13,17,18,26]. Interestingly, in peripheral blood these changes were associated with increased expression of RBM15 [17]. Moreover, Qiu et al. [26] proposed a predictive “m6A score” to quantify and model the m6A pattern in blood leukocytes for each COVID-19 patient based on m6A levels and nine selected differentially expressed genes (m6A-DEGs) mostly related to the immune response; in this model, patients displaying higher (protective) scores showed a better prognosis related to T-cell activation compared to patients with lower scores; in addition to clinical prognosis, the model may predict the possibility of contracting COVID-19 in patients infected with SARS-CoV-2, as well as the detection of SARS-CoV-2 carriers [26].
Our study also revealed significantly lower levels of m6A in samples from fully vaccinated compared to unvaccinated individuals (n = 8 and 52, respectively). Most unvaccinated individuals were positive for alpha, gamma, epsilon, and B.1.1.519 variants, whereas delta and omicron were the dominant variants in vaccinated individuals. Importantly, we show that the variation in m6A levels across variants, in particular the low levels in delta and omicron, could not be explained by vaccination status alone. While changes in DNA methylation levels have been reported after influenza vaccination [27,28], there is no information regarding RNA methylation after vaccination. Potential vaccination-mediated m6A is intriguing and deserves further investigation.
The delta variant showed significant differences in the mean levels of N expression when compared with the B.1.1.519 variant, indicating that the delta variant had the highest viral loads. Similarly, additional studies show that patients infected with the delta variant presented increased viral loads compared to other SARS-CoV-2 variants [29,30]. Conversely, a recent study that compared a larger number of individuals positive for alpha or delta to other variants (n = 36 and 41, respectively), neither found significant differences between mean viral load across variants, nor between vaccinated and unvaccinated participants (37 vaccinated and 136 unvaccinated) [31]. In our study, no significant differences in N gene expression were observed among VOC- and VOI-positive samples, consistent with the previous report.
The results presented in this study are preliminary and should be interpreted with caution due to the following limitations. First, although all patients enrolled in this study were Mexican and living in the city of Mazatlán, genetic background, health and immune conditions, as well as lifestyle, greatly differ across individuals. Furthermore, we received samples from different hospitals which may differ in the sampling process and/or the social/economic status of the patients. Second, global m6A data represented the sum of methylation of both human and viral RNAs; therefore, additional studies are needed to determine their relative contributions. Third, due to the small amounts of RNA obtained from each nasopharyngeal sample, we were unable to perform technical duplicates of m6A content. Fourth, since we only had access to amplification data of the viral N and human Rp genes, calculations of viral N expression were performed relative to Rp. Fifth, increasing the sample size of analyzed individuals is important in order to improve the reliability of the conclusions presented.
The port of Mazatlán is a very touristic location, receiving thousands of visitors from Mexico and other countries around the world every year. Rapid changes of the dominant variants over the sampling period could be due to the crowded environment allowing the virus to spread. In addition, high mutation and recombination rates of RNA viruses, as well as asymptomatic infections, could also contribute to the rapid diversification of SARS-CoV-2 [32] even among fully vaccinated individuals. Despite the limitations of this study (which can be addressed in future studies), our results suggest that viral variants modulate genome and/or viral m6A levels differentially (Figure 5). Further research is needed to fully unravel and understand the molecular basis of these differences and its relevance in viral infection and transmission, as well as the implications in human immunity and health.

4. Materials and Methods

4.1. Samples

This study used 404 SARS-CoV-2 PCR-positive sequenced samples from the National Epidemiological Surveillance Consortium (see Figure 1), from which 60 samples (plus 10 samples from SARS-CoV-2 negative patients) were randomly selected for further analysis. Patients <30 or >60 years old, as well as SARS-CoV-2 PCR-positive samples with Cts <30 were excluded.

4.2. Identification and Sequencing of SARS-CoV-2 Positive Samples

Nasopharyngeal samples were taken by health professionals from different hospitals in Mazatlán (Mexico) and transported to the Molecular Diagnostic Laboratory at the Research Centre for Nutrition and Development (LDM-CIAD, official authorization number: DGE-DDYR-DSAT-04471-2020) for SARS-CoV-2 detection by real-time PCR (2019-nCov CDC® Integrated DNA Technologies, Coralville, IA, USA). Relative expression levels of the viral N gene were determined for each sample as follows [33]:
Δ C T = C T N 1 + C T N 2 2 C T R p
where N1 and N2 are regions of the viral N gene and Rp is the human RNAse P gene (Supplementary File S1).
RNA samples (n = 404) with higher N levels (CT < 30) were selected for sequencing through the National Epidemiological Surveillance Consortium (COVIGEN). Complementary DNAs (cDNAs) were synthesized from total RNAs using the GoScript™ Reverse Transcription System (Promega, Madison, WI, USA). Libraries were prepared for sequencing using the COVID-Seq test (Illumina, San Diego, CA, USA), and the MiniSeqTM sequencing system.

4.3. Global RNA Methylation Assay

Sequenced SARS-CoV-2-positive RNA samples (n = 60; 10 per variant) were randomly selected from adult (30–60 years old) patients (men and women). In addition, 10 negative samples were used as controls. RNA concentration and purity were measured with a spectrophotometer (Nanodrop, Thermo Fisher Scientific, Waltham, MA, USA Thermo). Because of the small amounts of RNA in nasopharyngeal samples, all samples were vacuum-lyophilized (LABCONCO, Kansas City, MO, USA) at −56 °C and 0.123 mBar, for 48 h. The samples were then resuspended in nuclease-free water for a final concentration of 60 ng/µL. The methylation assay was performed with the EpiQuikTM m6A RNA Methylation Quantification Kit (Epigentek, Farmingdale, NY, USA), according to the manufacturer’s instructions. Absorbance at 450 nm was measured in a microplate reader (EPOCH2, BioTek, Thermo Fisher Scientific, Waltham, MA, USA). The standard curve was obtained and the relative methylation level m6A (%) was determined according to the following formula:
m 6 A   ( % ) = ( O D s a m p l e O D N C ) ÷ S ( O D P C O D N C ) ÷ P × 100
where ODsample is the optical density of the sample, ODNC is the optical density of the negative control, ODPC is the optical density of the positive control, S is the amount of RNA (ng) applied to the well, and P is the amount of RNA (ng) applied to the positive control well (Supplementary File S2).

4.4. Statistical Analysis

Normality and homoscedasticity were evaluated with Shapiro-Wilk and Leneve tests, respectively. Kruskal–Wallis and post-hoc Dunn tests were performed with relative m6A RNA methylation and N gene expression data to detect differences among variants. Spearman correlations between relative m6A methylation and relative expression levels of the N gene or age were also done. In addition, differences in m6A methylation levels among vaccination schemes were analyzed by Kruskal–Wallis (unbalanced) and Dunn tests. The significance level was 0.05. All analyses were performed in RStudio 4.0.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/epigenomes6030016/s1, Supplementary File S1: Real-time PCR data; Supplementary File S2: m6A RNA methylation data.

Author Contributions

Conceptualization, A.G.-G. and J.L.B.-R.; methodology, A.G.-G., J.L.B.-R., B.G.-G. and G.L.; formal analysis, J.L.B.-R., C.A.B.-R. and B.G.-G.; investigation, A.G.-G. and J.L.B.-R.; resources, A.G.-G. and B.G.-G.; data curation, J.L.B.-R., B.G.-G. and A.G.-G.; writing—original draft preparation, A.G.-G. and J.L.B.-R.; writing—review and editing, A.G.-G., J.L.B.-R., B.G.-G., G.L. and C.A.B.-R.; visualization, A.G.-G., J.L.B.-R. and B.G.-G.; supervision, A.G.-G.; project administration, B.G.-G. and A.G.-G.; funding acquisition, B.G.-G. and A.G.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the LDM-COVID project 2017710752 awarded to B.G.G. and A.G.-G., and the project “Vigilancia Genómica del Virus SARS-CoV-2 (CoViGen) en México-2022” (PP-F003) from the National Council for Science and Technology (CONACyT) awarded to B.G.G. The APC was funded by project 2017710752.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Ethics Committee for Research at CIAD (CONBIOETICA-26-CEI-001-20200122CEI/028-1/2021).

Informed Consent Statement

Patient consent was waived because of the health emergency and retrospective nature of the study. All health providers were required to inform the patients that their samples could have one out of three different destinations: (1) destruction; (2) storage (in case the Mexican authority needs the sample), or (3) research and sequencing. All patients used in this study were identified with a code provided by the health provider, without sharing any personal information.

Data Availability Statement

Data supporting reported results are provided as supplementary files.

Acknowledgments

The authors would like to thank Julissa Enciso-Ibarra, Alejandra Díaz-Sánchez, Karla Aguilar-Rendón, Alexia Audelo-Fajardo, Eunice Barraza-Ortega, and Víctor Álvarez-López for technical assistance in processing (RNA extraction, real-time PCR, and sequencing) COVID samples. We thank Crisantema Hernández for lending the microplate reader, and Nayely Leyva-López, and Marian Mirabent-Casals for technical assistance with the methylation assay. We are also grateful to Daniel Fregoso-Rueda for logistic and administrative support, and Carlos F. Arias and Blanca Taboada from the CoViGen project. We thank all anonymous reviewers for their constructive comments to this document.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mariano, G.; Farthing, R.J.; Lale-Farjat, S.L.M.; Bergeron, J.R.C. Structural characterization of SARS-CoV-2: Where we are, and where we need to be. Front. Mol. Biosci. 2020, 7, 605236. [Google Scholar] [CrossRef] [PubMed]
  2. Bai, Z.; Cao, Y.; Liu, W.; Li, J. The SARS-CoV-2 nucleocapsid protein and its role in viral structure, biological functions, and a potential target for drug or vaccine mitigation. Viruses 2021, 13, 1115. [Google Scholar] [CrossRef] [PubMed]
  3. Le Bert, N.; Tan, A.T.; Kunasegaran, K.; Tham, C.Y.L.; Hafezi, M.; Chia, A.; Chng, M.H.Y.; Lin, M.; Tan, N.; Linster, M.; et al. SARS-CoV-2-specific T cell immunity in cases of COVID-19 and SARS, and uninfected controls. Nature 2020, 584, 457–462. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, Y.; Liu, Q.; Guo, D. Emerging coronaviruses: Genome structure, replication, and pathogenesis. J. Med. Virol. 2020, 92, 418–423. [Google Scholar] [CrossRef]
  5. Davies, N.G.; Abbott, S.; Barnard, R.C.; Jarvis, C.I.; Kucharski, A.J.; Munday, J.D.; Pearson, C.A.B.; Russell, T.W.; Tully, D.C.; Washburne, A.D.; et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science 2021, 372, eabg3055. [Google Scholar] [CrossRef]
  6. Bager, P.; Wohlfahrt, J.; Fonager, J.; Rasmussen, M.; Albertsen, M.; Michaelsen, T.Y.; Møller, C.H.; Ethelberg, S.; Legarth, R.; Button, M.; et al. Risk of hospitalisation associated with infection with SARS-CoV-2 lineage B.1.1.7 in Denmark: An observational cohort study. Lancet Infect. Dis. 2021, 21, 1507–1517. [Google Scholar] [CrossRef]
  7. Desrosiers, R.; Friderici, K.; Rottman, F. Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells. Proc. Natl. Acad. Sci. USA 1974, 71, 3971–3975. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, M.; Wong, C.M. The emerging roles of N6-methyladenosine (m6A) deregulation in liver carcinogenesis. Mol. Cancer 2020, 19, 44. [Google Scholar] [CrossRef] [Green Version]
  9. Manners, O.; Baquero-Perez, B.; Whitehouse, A. m6A: Widespread regulatory control in virus replication. Biochim. Biophys. Acta Gene Regul. Mech. 2019, 1862, 370–381. [Google Scholar] [CrossRef]
  10. Fu, Y.; Dominissini, D.; Rechavi, G.; He, C. Gene expression regulation mediated through reversible m6A RNA methylation. Nat. Rev. Genet. 2014, 15, 293–306. [Google Scholar] [CrossRef]
  11. Ke, S.; Pandya-Jones, A.; Saito, Y.; Fak, J.J.; Vågbø, C.B.; Geula, S.; Hanna, J.H.; Black, D.L.; Darnell, J.E., Jr.; Darnell, R.B. m6A mRNA modifications are deposited in nascent pre-mRNA and are not required for splicing but do specify cytoplasmic turnover. Genes Dev. 2017, 31, 990–1006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Zaccara, S.; Ries, J.R.; Jaffrey, R.S. Reading, writing and erasing mRNA methylation. Nat. Rev. Mol. Cell Biol. 2019, 20, 608–624. [Google Scholar] [CrossRef] [PubMed]
  13. Li, N.; Hui, H.; Bray, B.; Gonzalez, G.M.; Zeller, M.; Anderson, K.G.; Knight, R.; Smith, D.; Wang, Y.; Carlin, A.F.; et al. METTL3 regulates viral m6A RNA modification and host cell innate immune responses during SARS-CoV-2 infection. Cell Rep. 2021, 35, 109091. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, J.; Xu, Y.P.; Li, K.; Ye, Q.; Zhou, H.Y.; Sun, H.; Li, X.; Yu, L.; Deng, Y.Q.; Li, R.T.; et al. The m6A methylome of SARS-CoV-2 in host cells. Cell Res. 2021, 31, 404–414. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, X.; Hao, H.; Ma, L.; Zhang, Y.; Hu, X.; Chen, Z.; Liu, D.; Yuan, J.; Hu, Z.; Guan, W. Methyltransferase-like 3 modulates severe acute respiratory syndrome coronavirus-2 RNA N6-methyladenosine modification and replication. mBio 2021, 12, e01067-21. [Google Scholar] [CrossRef]
  16. Burgess, H.M.; Depledge, D.P.; Thompson, L.; Srinivas, K.P.; Grande, R.C.; Vink, E.I.; Abebe, J.S.; Blackaby, W.P.; Hendrick, A.; Albertella, M.R.; et al. Targeting the m6A RNA modification pathway blocks SARS-CoV-2 and HCoV-OC43 replication. Genes Dev. 2021, 35, 1005–1019. [Google Scholar] [CrossRef]
  17. Meng, Y.; Zhang, Q.; Wang, K.; Zhang, X.; Yang, R.; Bi, K.; Chen, W.; Diao, H. RBM15-mediated N6-methyladenosine modification affects COVID-19 severity by regulating the expression of multitarget genes. Cell Death Dis. 2021, 12, 732. [Google Scholar] [CrossRef]
  18. An, S.; Xie, Z.; Liao, Y.; Jiang, J.; Dong, W.; Yin, F.; Li, W.X.; Ye, L.; Lin, J.; Liang, H. Systematic analysis of clinical relevance and molecular characterization of m6A in COVID-19 patients. Genes Dis. 2022, in press. [Google Scholar] [CrossRef]
  19. Liu, C.; Yang, Z.; Li, R.; Wu, Y.; Chi, M.; Gao, S.; Sun, X.; Meng, X.; Wang, B. Potential roles of N6-methyladenosine (m6A) in immune cells. J. Transl. Med. 2021, 19, 251. [Google Scholar] [CrossRef]
  20. Lu, M.; Xue, M.; Wang, H.T.; Kairis, E.L.; Ahmad, S.; Wei, J.; Zhang, Z.; Liu, Q.; Zhang, Y.; Gao, Y.; et al. Nonsegmented negative-sense RNA viruses utilize N6-methyladenosine (m6A) as a common strategy to evade host innate immunity. J. Virol. 2021, 95, e01939-20. [Google Scholar] [CrossRef]
  21. Brocard, M.; Ruggieri, A.; Locker, N. m6A RNA methylation, a new hallmark in virus-host interactions. J. Gen. Virol. 2017, 98, 2207–2214. [Google Scholar] [CrossRef] [PubMed]
  22. Statement on Omicron Sublineage BA.2. Available online: https://www.who.int/news/item/22-02-2022-statement-on-omicron-sublineage-ba.2 (accessed on 29 March 2022).
  23. Da Silva, S.; de Lima, S.C.; da Silva, R.C.; Kohl, A.; Pena, L. Viral load in COVID-19 patients: Implications for prognosis and vaccine efficacy in the context of emerging SARS-CoV-2 variants. Front. Med. 2022, 8, 836826. [Google Scholar] [CrossRef] [PubMed]
  24. McCallum, M.; Czudnochowski, N.; Rosen, L.E.; Zepeda, S.K.; Bowen, J.E.; Walls, A.C.; Hauser, K.; Joshi, A.; Stewart, C.; Dillen, J.R.; et al. Structural basis of SARS-CoV-2 Omicron immune evasion and receptor engagement. Science 2022, 375, 864–868. [Google Scholar] [CrossRef] [PubMed]
  25. Campos, J.H.C.; Maricato, J.T.; Braconi, C.T.; Antoneli, F.; Janini, L.M.R.; Briones, M.R.S. Direct RNA sequencing reveals SARS-CoV-2 m6A sites and possible differential drach motif methylation among variants. Viruses 2021, 13, 2108. [Google Scholar] [CrossRef]
  26. Qiu, X.; Hua, X.; Li, Q.; Zhou, Q.; Chen, J. m6A regulator-mediated methylation modification patterns and characteristics of immunity in blood leukocytes of COVID-19 patients. Front. Immunol. 2021, 12, 774776. [Google Scholar] [CrossRef]
  27. Gensous, N.; Franceschi, C.; Blomberg, B.B.; Pirazzini, C.; Ravaioli, F.; Gentilini, D.; Di Blasio, A.M.; Garagnani, P.; Frasca, D.; Bacalini, M.G. Responders and non-responders to influenza vaccination: A DNA methylation approach on blood cells. Exp. Gerontol. 2018, 105, 94–100. [Google Scholar] [CrossRef]
  28. Zimmermann, M.T.; Oberg, A.L.; Grill, D.E.; Ovsyannikova, I.G.; Haralambieva, I.H.; Kennedy, R.B.; Poland, G.A. System-wide associations between DNA methylation, gene expression, and humoral immune response to influenza vaccination. PLoS ONE 2016, 11, e0152034. [Google Scholar] [CrossRef]
  29. Luo, C.H.; Morris, C.P.; Sachithanandham, J.; Amadi, A.; Gaston, D.C.; Li, M.; Swanson, N.J.; Schwartz, M.; Klein, E.Y.; Pekosz, A.; et al. Infection with the SARS-CoV-2 Delta variant is associated with higher recovery of infectious virus compared to the Alpha variant in both unvaccinated and vaccinated individuals. Clin. Infect. Dis. 2021, 18, ciab986. [Google Scholar] [CrossRef]
  30. Ong, S.; Chiew, C.J.; Ang, L.W.; Mak, T.M.; Cui, L.; Toh, M.; Lim, Y.D.; Lee, P.H.; Lee, T.H.; Chia, P.Y.; et al. Clinical and virological features of SARS-CoV-2 variants of concern: A retrospective cohort study comparing B.1.1.7 (Alpha), B.1.315 (Beta), and B.1.617.2 (Delta). Clin. Infect. Dis. 2021, ciab721. [Google Scholar] [CrossRef]
  31. Kissler, S.M.; Fauver, J.R.; Mack, C.; Tai, C.G.; Breban, M.I.; Watkins, A.E.; Samant, R.M.; Anderson, D.J.; Metti, J.; Khullar, G.; et al. Viral dynamics of SARS-CoV-2 variants in vaccinated and unvaccinated persons. N. Engl. J. Med. 2021, 385, 2489–2491. [Google Scholar] [CrossRef]
  32. Duarte, C.M.; Ketcheson, D.I.; Eguíluz, V.M.; Agustí, S.; Fernández-Gracia, J.; Jamil, T.; Laiolo, E.; Gojobori, T.; Alam, I. Rapid evolution of SARS-CoV-2 challenges human defenses. Sci. Rep. 2022, 12, 6457. [Google Scholar] [CrossRef] [PubMed]
  33. Rao, X.; Huang, X.; Zhou, Z.; Lin, X. An improvement of the 2ˆ(-delta delta CT) method for quantitative real-time polymerase chain reaction data analysis. Biostat. Bioinform. Biomath. 2013, 3, 71–85. [Google Scholar]
Figure 1. Percentage of genomes of SARS-CoV-2 variants sequenced in Mazatlán, Mexico, from February 2021 to March 2022 (n = 404). Delta (AY. other) = AY.100, AY.103, AY.113, AY.119, AY.122.4, AY.25, AY.3, AY.39, AY.43, AY.44. AY.53. Omicron (other) = BA.1.1.16, BA.1.1.2, BA.1.1.8, BA.1.13, BA.1.14, BA.1.17.
Figure 1. Percentage of genomes of SARS-CoV-2 variants sequenced in Mazatlán, Mexico, from February 2021 to March 2022 (n = 404). Delta (AY. other) = AY.100, AY.103, AY.113, AY.119, AY.122.4, AY.25, AY.3, AY.39, AY.43, AY.44. AY.53. Omicron (other) = BA.1.1.16, BA.1.1.2, BA.1.1.8, BA.1.13, BA.1.14, BA.1.17.
Epigenomes 06 00016 g001
Figure 2. (A) Global m6A levels of RNAs extracted from human nasopharyngeal samples, tested for SARS-CoV-2 by real-time PCR based on viral variant (n = 10 each). Control: SARS-CoV-2 negative samples (Kruskal-Wallis p = 3.2 × 10−7); different letters indicate significant differences in methylation levels. (B) Relative expression of the N gene in RNAs extracted from human nasopharyngeal samples, tested for SARS-CoV-2 by real-time PCR based on viral variant (n = 10 each) (Kruskal-Wallis, p = 0.041; Dunn, p = 0.024); different letters indicate significant differences in gene expression levels. VOC: omicron, delta, gamma, and alpha. VOI: epsilon; B.1.1.519: lineage found in Mexico during the second wave.
Figure 2. (A) Global m6A levels of RNAs extracted from human nasopharyngeal samples, tested for SARS-CoV-2 by real-time PCR based on viral variant (n = 10 each). Control: SARS-CoV-2 negative samples (Kruskal-Wallis p = 3.2 × 10−7); different letters indicate significant differences in methylation levels. (B) Relative expression of the N gene in RNAs extracted from human nasopharyngeal samples, tested for SARS-CoV-2 by real-time PCR based on viral variant (n = 10 each) (Kruskal-Wallis, p = 0.041; Dunn, p = 0.024); different letters indicate significant differences in gene expression levels. VOC: omicron, delta, gamma, and alpha. VOI: epsilon; B.1.1.519: lineage found in Mexico during the second wave.
Epigenomes 06 00016 g002
Figure 3. Scatter plot between relative expression of the N gene and m6A levels (Spearman’s rank correlation, rho = −0.063, p = 0.63).
Figure 3. Scatter plot between relative expression of the N gene and m6A levels (Spearman’s rank correlation, rho = −0.063, p = 0.63).
Epigenomes 06 00016 g003
Figure 4. (A) Global m6A levels in of RNAs extracted from human nasopharyngeal samples, tested for SARS-CoV-2 by real-time PCR with different vaccination schemes (n = 64). Partial and complete vaccination refers to the application of one or two doses, respectively, of authorized vaccines in Mexico (AZD1222, CoronaVac, BNT162b2, Ad5-nCOV, and mRNA-1273 from the companies AstraZeneca, Sinovac, Pfizer-BioNTech, CanSinoBio, and Moderna, respectively) (Kruskal-Wallis, p = 2.6 × 10−4; Dunn test, p = 2.5 × 10−4); n = 8, 4, and 52 for complete, partial, and unvaccinated individuals, respectively. (B) Global m6A levels of RNAs extracted from human nasopharyngeal samples from unvaccinated patients, tested for SARS-CoV-2 by real-time PCR based on viral variant (n = 52). Control: SARS-CoV-2 negative samples (Kruskal–Wallis p = 3.068 × 10−3); different letters indicate significant differences in methylation levels.
Figure 4. (A) Global m6A levels in of RNAs extracted from human nasopharyngeal samples, tested for SARS-CoV-2 by real-time PCR with different vaccination schemes (n = 64). Partial and complete vaccination refers to the application of one or two doses, respectively, of authorized vaccines in Mexico (AZD1222, CoronaVac, BNT162b2, Ad5-nCOV, and mRNA-1273 from the companies AstraZeneca, Sinovac, Pfizer-BioNTech, CanSinoBio, and Moderna, respectively) (Kruskal-Wallis, p = 2.6 × 10−4; Dunn test, p = 2.5 × 10−4); n = 8, 4, and 52 for complete, partial, and unvaccinated individuals, respectively. (B) Global m6A levels of RNAs extracted from human nasopharyngeal samples from unvaccinated patients, tested for SARS-CoV-2 by real-time PCR based on viral variant (n = 52). Control: SARS-CoV-2 negative samples (Kruskal–Wallis p = 3.068 × 10−3); different letters indicate significant differences in methylation levels.
Epigenomes 06 00016 g004
Figure 5. Schematic representation of m6A levels (grey triangle) among different SARS-CoV-2 variants from February 2021 to March 2022. Circle size represents number of cases per variant (smaller circles = less cases; larger circles = more cases) in the studied population.
Figure 5. Schematic representation of m6A levels (grey triangle) among different SARS-CoV-2 variants from February 2021 to March 2022. Circle size represents number of cases per variant (smaller circles = less cases; larger circles = more cases) in the studied population.
Epigenomes 06 00016 g005
Table 1. Characteristics of patients (average age, sex distribution, and vaccination status) included in each variant group. The Control group refers to SARS-CoV-2 negative patients.
Table 1. Characteristics of patients (average age, sex distribution, and vaccination status) included in each variant group. The Control group refers to SARS-CoV-2 negative patients.
VariantAverage Age ± SDSex Distribution (n = 70)Vaccination Status (n = 64)
MalesFemales
Control46 ± 755Not vaccinated (n = 10)
Omicron42 ± 855Complete (n = 8)
Not vaccinated (n = 2)
Delta39 ± 1055Partial (n = 4)
Gamma44 ± 755Not vaccinated (n = 10)
Alpha49 ± 855Not vaccinated (n = 10)
Epsilon45 ± 1055Not vaccinated (n = 10)
B.1.1.51947 ± 1155Not vaccinated (n = 10)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Batista-Roche, J.L.; Gómez-Gil, B.; Lund, G.; Berlanga-Robles, C.A.; García-Gasca, A. Global m6A RNA Methylation in SARS-CoV-2 Positive Nasopharyngeal Samples in a Mexican Population: A First Approximation Study. Epigenomes 2022, 6, 16. https://doi.org/10.3390/epigenomes6030016

AMA Style

Batista-Roche JL, Gómez-Gil B, Lund G, Berlanga-Robles CA, García-Gasca A. Global m6A RNA Methylation in SARS-CoV-2 Positive Nasopharyngeal Samples in a Mexican Population: A First Approximation Study. Epigenomes. 2022; 6(3):16. https://doi.org/10.3390/epigenomes6030016

Chicago/Turabian Style

Batista-Roche, Jorge Luis, Bruno Gómez-Gil, Gertrud Lund, César Alejandro Berlanga-Robles, and Alejandra García-Gasca. 2022. "Global m6A RNA Methylation in SARS-CoV-2 Positive Nasopharyngeal Samples in a Mexican Population: A First Approximation Study" Epigenomes 6, no. 3: 16. https://doi.org/10.3390/epigenomes6030016

APA Style

Batista-Roche, J. L., Gómez-Gil, B., Lund, G., Berlanga-Robles, C. A., & García-Gasca, A. (2022). Global m6A RNA Methylation in SARS-CoV-2 Positive Nasopharyngeal Samples in a Mexican Population: A First Approximation Study. Epigenomes, 6(3), 16. https://doi.org/10.3390/epigenomes6030016

Article Metrics

Back to TopTop