1. Introduction

In the central dogma of molecular biology, genetic information flows from DNA to RNA and then to proteins. In addition to the canonical nucleobases, plenty of chemical modifications have been identified presence in genomic DNA and RNA [1, 2]. These modifications do not alter the sequence context of nucleic acids, but change their biochemical and physical properties, and eventually control or regulate the spatiotemporal expression of genes [3]. Up to date, over 150 different types of chemical modifications have been found existence in nucleic acids [4, 5]. The modifications in nucleic acids play sophisticated roles in annotating genetic information [6, 7].

So far, seven modifications have been discovered in genomic DNA ofeukaryotes[7-9].5-Methylcytosine(m⁵C)is the most abundant and important epigenetic mark in DNA [10].m⁵C modification in DNA is dynamic and reversible and involves in numerous physiological processes[10].In2009, 5-hydroxymethylcytosine(hm⁵C), anoxidized derivative of m⁵C, was discovered in mammalian genome and TET (ten-eleven translocation) proteins were identified to be able to convert m⁵C to hm⁵C [11, 12]. Subsequent studies revealed that TET proteins can further oxidize hm⁵C to 5-formylcytosine (f⁵C) and 5-carboxycytosine (ca⁵C) [13, 14]. With these discovered modifications, these studies uncovered the mechanism of active DNA demethylation in mammals. Moreover, it was recently reported that DNA adenine methylation (N⁶-methyladenine, m⁶A) and the hydroxylation derivative of m⁶A were also present in genomic DNA of eukaryotes [9, 15-18]. These adenine modifications in DNA show crucial roles in regulating expression of genes [8].

In addition to DNA modifications, RNA molecules also carry various levels of modifications [19]. More than 150 different kinds of modifications exist in cellular RNA [2]. These modifications exert essential and critical functions in a great number of cellular processes in eukaryotic organisms, including affecting the stability, maturation, and secondary structures of RNA [19, 20]. Besides maintaining the higher structure and catalytic activity, RNA modifications also function as new information carrier in regulating cell physiology [19, 20].

Deciphering the biological roles of modifications in nucleic acids relies on the identification and location analysis of these modifications[21]. Methodsfor analyzing modifications in nucleic acids have tremendously advanced in the last decade [21-29]. Detection of nucleic acid modifications that relies on enzymatic digestion followed by liquid chromatography [30], thin layerchromatography (TLC) [31], and capillary electrophoresis [32], generally only provide the overall contents of modifications. Location information of these modifications in nucleic acids is essential and critical to reveal the functions of modifications. To map modifications in nucleic acids, some masss pectrometry(MS)-based analytical strategies have been developed [33]. Recent progress in next-generation sequencing (NGS) in conjugationwith immunoprecipitation, chemical reaction, or enzyme-mediated mutation allow genome-wide or transcriptome-wide mapping of modifications, which greatly revolutionize the field of nucleic acid modifications [34].

Herein, we review the established methods and the breakthrough of the techniques for location analysis of nucleic acid modifications. Fig. 1 is the summary of the analytical methods for location analysis of modifications in nucleic acids. We discuss the principles, applications, advantages and drawbacks of these methods. We believe that the advancement of techniques and methods will promote uncovering the functional roles of nucleic acid modifications.

2. Location analysis of modifications in nucleic acids 2.1. Mass spectrometry

Mass spectrometry has superior capability in identification of compounds [35-39]. The rapid advances in MS-based analytical methods greatly promote the study of nucleic acid modifications and it becomes increasing important platform in the research field of nucleic acid modifications [40-46]. Tandem MS (MS/MS) is commonly performed to analyze the sequence of nucleic acids and modifications in nucleic acids [47]. Particularly, liquid chromatography-electrospray ionization/mass spectrometry (LC-ESI/MS), matrix-assisted laser desorption/ionization mass spectrometry (MALDI/MS), and top-down mass spectrometry (top-down MS) are the most widely used platforms for MS-base location analysis of modifications in nucleic acids (Fig. 2).

2.1.1. LC-ESI/MS

As nucleic acids typically are long polymers and have extremely large molecular weights, they generally are first enzymatically digested to small fragments followed by LC-ESI/MS analysis [48-54]. RNase T1 and RNase A are widely used in the digestion of RNA to form small oligonucleotide fragments, which were then subjected to LC-ESI/MS [55, 56]. Tandem MS (MS/MS) analysis can generate sequence information, including the sites of modifications on oligonucleotides [57, 58]. With this analytical strategy, the sites of various modifications in RNA were identified, including 3-methylpseudouridine (m³U) in 23S rRNA of E. coli [59], 5-methyluridine (m⁵U) and pseudouridine (ψ) in tmRNA of E. coli [60], 2-thiouridine (s²U) in mitochondria tRNA of L. tarentolae [61], 2-methylthio cyclic N⁶-threonylcarbamoyladenosine (ms²ct⁶A) at position 37 of tRNA from Bacillus subtilis, plants and Trypanosoma brucei [62]. LC-ESI/MS is only able to analyze short fragments of DNA and RNA, which limits its application on the genome-wide or transcriptome-wide mapping of modifications.

2.1.2. MALDI/MS

In addition to LC-ESI/MS, MALDI/MS is also frequently used for the analysis of nucleic acid fragments due to its capability on detection of relatively large molecular weight compounds. Prior to MALDI/MS analysis, RNA was also need to be fragmented by specific RNase, which are then analyzed by MALDI-MS/MS to locate the modifications [63]. The MALDI/MS was successfully employed in screening multiple modifications in rRNA [64].

ψ is a mass-silent modification and the normal MALDI/MS analysis cannot locate ψ in RNA. To address this issue, chemical labeling in combination with MALDI/MS analysis was developed for locating ψ modification in RNA [65]. In this method, 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide (CMC) was used to selectively label ψ and the formed derivative of ψ carrying a 252 Da mass tag that can be easily identified by subsequent MALDI-MS/MS analysis [65]. In addition, the specific cyanoethylation of N1 in ψ by acrylonitrile was also utilized to locate ψ in RNA [66]. Cyanoethylated fragments can be readily identified by comparing the MS spectra of untreated and acrylonitrile-treated samples since the addition of one acrylonitrile leads to a mass increase of 53 Da. With this analytical strategy, ψ in E. coli tRNA^TyrII was identified. Cyanoethylation is a useful complement to the CMC labeling on ψ location analysis.

2.1.3. Top-down MS

An alternative strategy for the location analysis of modifications in nucleic acids is through the top-down MS analysis of the intact RNA without enzymatic digestion. For example, top-down MS was successfully used for the identification of the presence of 2'-O-methylation on the terminal of microRNA in plant [67]. With the top-down MS analysis, 89% sequence coverage in a single collisionally activated dissociation spectrum of E. coli tRNA^Val was obtained and many modifications in tRNA also can be located [68]. Very recently, top-down MS was used to characterize RNA methylations (m⁶A, m⁵C, m³U, and m⁵U) in the mixture of either isomers of RNA or non-isomeric RNA forms [69]. The top-down MS-based analysis relies on advanced mass spectrometers. Right now only pure RNA samples can be analyzed by top-down MS. Thus, this method still needs great improvement for locating modifications in nucleic acids from biological samples.

2.2. Sequencing

Sequencing can provide comprehensive information of nucleotides in DNA and RNA. In the last several years, many powerful strategies based on sequencing in conjugation with a variety of approaches have been developed and established for genomewide and transcriptome-wide mapping of a variety of modifications in nucleic acids.

2.2.1. Immunoprecipitation-mediated sequencing

Nucleic acid modifications generally are present in lowabundance. Therefore, the enrichment of the fragments that carry modifications will benefit the subsequent sequencing for mapping modifications. Immunoprecipitation-mediated enrichment of nucleic acid fragments using antibodies or affinity binding proteins has been demonstrated to be very useful for mapping modifications in nucleic acids [70]. The drawback of immunoprecipitation-mediated sequencing is that this strategy generally cannot offer location information of modifications at single-base resolution.

Many immunoprecipitation-mediated sequencing, such as methylated DNA binding domain sequencing (MBD-seq) [71], methylation DNA immunoprecipitation sequencing (MeDIP-seq) [72], MIRA-seq [73], MethylCap-seq [74], and m⁶A DIP-seq [16, 75, 76] have been established for the comprehensive mapping of m⁵C and m⁶A in genomic DNA. In addition, hm⁵C DNA immunoprecipitation-sequencing using antibodies against hm⁵C [77, 78], or antibody against cytosine 5-methylenesulfonate (CMS) [79] were also developed to map hm⁵C in genome. These mapping studies demonstrated that hm⁵C was mainly enriched in enhancers, promoters, and gene bodies, revealing a critical role of hm⁵C in regulation of gene expressions. Similarly, immunoprecipitation-mediated enrichment of f⁵C followed by NGS realized the genome-wide mapping of f⁵C [80].

In addition to DNA modifications, immunoprecipitation of RNA fragments carrying modifications followed by NGS (RIP-seq) offers a promising platform for transcriptome-wide mapping of RNA modifications. In recent years, m⁶A-specific RNA immunoprecipitation coupled with NGS (MeRIP-seq) was developed for transcriptome-wide mapping of m⁶A, which provided an effective approach for locating m⁶A in mRNA and revealed the regulatory roles of m⁶A [81-83]. However, MeRIP-seq only can map m⁶A in RNA at a resolution of approximately 100-200 nt. Recently, crosslinking and immunoprecipitation (miCLIP) in conjugation NGS analysis was developed for mapping m⁶A at single-base-resolution [84]. Moreover, the m¹A-ID-seq [85, 86] and ac⁴C-RIP-seq [87] for transcriptome-wide mapping of m¹A and N⁴-acetylcytidine (ac⁴C), respectively, based on immunoprecipitation-mediated sequencing were also established. Collectively, these approaches allow comprehensive location analysis of modifications and provide valuable platforms for functional studies of modifications in nucleic acids.

2.2.2. Chemical reaction-mediated sequencing

The bisulphite-converted strategy is the gold standard for mapping m⁵C in DNA. Bisulphite sequencing can provide single-base resolution of m⁵C in genomic DNA [88]. After the bisulphite treatments, the genome is largely composed of only three nucleotides (A, G, and T). Thus, methods relying on bisulphite sequencing need extensive bioinformatics analysis.

Inadditionto m⁵C inDNA, oxidative bisulphite sequencing(oxBS-seq)strategy was developed for mappinghm⁵C inDNA atsingle-base resolution [89]. The oxBS-seq approach is based on the specific oxidation of hm⁵C to form f⁵C by potassium perruthenate (KRuO₄), which enables conversion of the formed f⁵C to uracil by bisulphite treatment. The hm⁵C would therefore read as thymine, whereas m⁵C would still read as cytosine in DNA treated with KRuO₄ by subsequent bisulphite sequencing. By comparing the sequencing results obtained by oxBS-seq and traditional bisulphite sequencing, the hm⁵C sites can be obtained at single-baseresolution.Morerecently, a method named hmC-CATCH that was based on selective hm⁵C oxidation by potassium ruthenate (K₂RuO₄), chemical labeling and subsequentC-to-Ttransition during PCR, was developed for mapping hm⁵C in genomes (Fig. 3) [90]. This method required only nanoscale genomic DNA, which could greatly benefit the hm⁵C mapping analysis with limited biological and clinical samples. In addition, peroxotungstate-mediated oxidation and sequencing analysis was developed for mapping hm⁵C at single-base resolution [91]. In this strategy, the oxidized productof hm⁵C, trihydroxylated thymine, can pair with dATP. Bycomparing the sequencing results before and after the peroxotungstate oxidation, hm⁵C sites in DNA could be specifically located (Fig. 4).

Moreover, a strategy of TET-assisted bisulphite sequencing (TAB-seq) was also developed for single-base resolution mapping of hm⁵C [92]. In TAB-seq appraoch, hm⁵C is glycosylated by β-glucosyltransferase (β-GT), which prevented from TET-mediated oxidation. Then TET enzyme treatment causes the oxidation of m⁵C and f⁵C to ca⁵C, which were read as thymine in subsequent bisulphite sequencing. The remaining cytosine signals are then from the protected hm⁵C. The advantage of TAB-seq is that there is no need for the comparison with traditional bisulphite sequencing and therefore hm⁵C can be identified directly.

Recently, a reduced bisulphite sequencing (redBS-seq) method was established to map f⁵C in DNA at single-base resolution [93].This redBS-seq is based on a selective reduction of f⁵C to hm⁵C followed by bisulphite sequencing. By combining redBS-seq and oxBS-seq, the genome-wide mapping of m⁵C, hm⁵C and f⁵C in mouse embryonic stem cells was achieved. In addition, Yi's group developed a bisulfite-free method for genome-wide analysis of f⁵C based on 1, 3-indandione labeling of f⁵C by and subsequent C-to-T transition in sequencing [94]. This strategy shows no noticeable DNA degradation and provides single-base resolution mapping of 5f C. As for ca⁵C, a chemical modification-assisted bisulphite sequencing (CAB-seq) was developed for locating ca⁵C in DNA at single-base resolution [95]. The carboxyl group of ca⁵C was labeled by the compound containing primary amine group, which endows the resistance of labeled ca⁵C to deamination by bisulphite treatment. Therefore, ca⁵C can be distinguished from other nucleobases during bisulphite sequencing.

Similarly, bisulfite-mediated sequencing was also reported for the location analysis of m⁵C in both rRNA and tRNA of Drosophila [96], m⁵C in mRNA of human cells [97] and in mRNA of S. solfataricus [98]. Since bisulphite treatment may result in substantial degradation of nucleic acids and RNA sample is more prone to degradation than DNA, the bisulphite treatment reaction need to be carefully optimized while analysis of RNA samples.

Alkaline hydrolysis followed by NGS analysis (RiboMeth-seq) was developed to map 2′-O-methylation sites in yeast rRNA [99]. Under alkaline condition, the ribose sugar with 2'-OH can be cleaved; however, the 2′-O-methylated ribose is resistant to alkaline cleavage. The differential properties of 2′-O-methylated and 2'-OH ribose were then utilized to map 2′-O-methylation in rRNA [100]. More recently, periodate oxidation reaction was also utilized to establish a method of Nm-seq for transcriptome-wide mapping of 2′-O-methylation at single-base-resolution by virtue of the differential reactivity of 2′-O-methylated and 2'0-OH nucleotides to periodate oxidation [101]. In addition, a method of thiouridine-to-cytidine-sequencing (TUC-seq) was developed for mapping 4-thiouridine (s⁴U) in RNA [102]. The reaction of s⁴U with osmium tetroxide (OsO₄) and ammonia can result in conversion of s⁴U into cytosine, which was then utilized for a comparative sequencing analysis. For untreated RNA, s⁴U is read as T while for OsO₄-treated RNA, s⁴U is read as C. This TUC-seq enabled the location analysis of s⁴U in RNA at single-base resolution.

Adenosine-to-inosine (A-to-I) editing is prevalent in RNA [103]. However, inosine cannot be directly detected by reverse transcription due to its base-pair with cytidine. Specifically chemical labeling with modifications may change their reverse transcription (RT) signature by RT arrest or misincorporation. In this respect, inosine chemical erasing (ICE) approach was established to map inosines in mRNA based on cyanoethylation reaction [104]. In ICE-seq, acrylonitrile was used to cyanoethylate inosine to form N1-cyanoethylinosine (ce¹I), which can compromise the base-pairing with cytidine and induce the RT stop. This unique signature can be obtained by subsequent NGS analysis. With this ICE-seq, 5072 inosine sites in mRNA of human brain tissue were identified [104], which provides a powerful approach to decipher the function of A-to-I editing. Similarly, the borohydride-based reduction of ac⁴C in RNA can induce the RT misincorporation, which provided a unique chemical signature and was successfully applied for the location analysis of ac⁴C in RNA (Fig. 5) [105]. In addition, an approach of N₃-CMC-enriched ψ sequencing (CeU-seq) was also developed to map ψ in mRNA of human cells [106]. In CeU-seq, N₃-CMC was used to specifically react with ψ on RNA followed by subsequent enrichment and NGS analysis. The CMC-ψ complex can cause the RT stop one base downstream of ψ, which can provide transcriptomewide mapping of ψ at single-base-resolution. Similarly, iodoacetamide (IAA) was used to specifically react with s⁴U to form alkylated s⁴U that preferably pairs with G than A and therefore causes the T-to-C mutation [107]. This signature of IAA labeled s⁴U can be used for map s⁴U in RNA at single-base resolution.

2.2.3. Enzyme-mediated sequencing

Our group and others recently reported an enzyme-mediated sequencing strategy for mapping hm⁵C in DNA [108, 109]. Previous report demonstrated that cytosine deaminases of APOBEC3A (apolipoproteinBmRNA-editing catalytic polypeptide-like 3 A)can deaminate cytosine to form uracil [110]. We further identified that APOBEC3A can effectively deaminate C, m⁵C, and hm⁵C, but show no observable deamination activity toward glycosylated hm⁵C (β-glucosyl-5-hydroxymethyl-2'-deoxycytidine, ghm⁵C). With this unique property of APOBEC3A, we established the APOBEC3 A-mediated deamination sequencing (AMD-seq) for location analysis of hm⁵C in DNA at single-base resolution (Fig. 6). In AMD-seq, APOBEC3A deaminate C and m⁵C in DNA to form U and T; while ghm⁵C is resistant to APOBEC3 Adeamination and will still read as C during sequencing (Fig. 7). Thus, the detected remaining C in DNA by sequencing could only come from original hm⁵C, which provides the single-base resolution analysis of hm⁵C. Since the AMD-seq analytical strategy does not need bisulphite treatment, the significant degradation of DNA can be avoided, which benefit the highlysensitive mappingof hm⁵C in DNAeven at single cell level.

2.2.4. Reverse transcription-mediated sequencing

The RT-mediated arrest or misincorporation in cDNA sequences can be analyzed by subsequent NGS. For example, AlkBfacilitated RNA methylation sequencing (ARM-seq) utilized the property that m¹A, m¹G, and m³C can cause RT-arrest to map their sites in tRNA (Fig. 8) [111]. In ARM-seq, the full-length cDNA can be obtained using AlkB-treated RNA since AlkB removed the methyl group of m¹A, m¹G, and m³C. However, the native RNA without AlkB treatment only produced truncated cDNA due to RTarrest caused by these modifications. Thus, the comparison of NGS results in treated and untreated samples can provide the modification sites in RNA. Recently, an antibody-independent method to map m⁶A at single-nucleotide resolution by 4SedTTP incorporation was established [112]. 4SedTTP can form stable base pair with adenine but not with m⁶A. Therefore, the use of the 4SedTTP instead of TTP in reverse transcription led to the RT stop at m⁶A sites, which was then utilized to locate m⁶A in RNA. In addition, an engineered RT-active KlenTaq DNA polymerase variant was employed for mapping m⁶A in RNA [113]. The KlenTaq DNA polymerase variant exhibited increased misincorporation opposite m⁶A compared to adenine, which allowed the direct identification of m⁶A sites from the sequencing data. However, the RT-mediated sequencing for mapping modifications in RNA can be largely affected by the RNA higher structures. Thus, other methods are essential to confirm the results obtained by RT-mediated sequencing analysis.

2.3. Endonuclease digestion-mediated location analysis

Some endonucleases can recognize and cleave specific modifications in RNA. The produced patterns can offer the location information of modifications by subsequent PCR, TLC, or gel electrophoresis analysis. Glyoxal can form stable complex with G, but not inosine, which therefore was used to map inosine sites in RNA [114]. Upon glyoxal treatment, the formed glyoxalG-borate complex is resistant to RNase T1 cleavage and only inosine sites were cleaved by RNase T1. The resulting fragments were then analyzed by gel electrophoresis, which provides a way to locate inosines in RNA. In addition, the RNase H cleavage in combination with gel electrophoresis analysis has been successfully employed in the location analysis of ψ in RNA [115], 2′-O-methylguanosine in 18S rRNA and 2′-O-methyladenosine in 28S rRNA of Xenopus [116].

Endonuclease digestion-mediated location analysis generally requires relatively large amount of purified RNA, which causes the analysis of modifications in individual mRNA or lncRNA challenging. Recently, a method of site-specific cleavage and radioactivelabeling followed by ligation-assisted extraction and thin-layer chromatography (SCARLET) was developed to map the sites of m⁶A in mRNA and lncRNA [117]. In the SCARLET method, RNase H digestion, ligation, and TLC analysis were sequentially carried out. With the SCARLET method, the m⁶A sites in two human lncRNA and three mRNA were successfully identified. SCARLET method requires no expensive instruments and can be easily applicable for mapping various modifications in RNA.

2.4. Reverse transcription-mediated location analysis

A number of modifications in RNA can lead to RT arrest or misincorporation of dNTP, which can be used as characterized signature for locating modifications in RNA by subsequent PCR or gel electrophoresis analysis [118]. For example, some modifications in RNA were specifically mapped by RT-mediated arrest or misincorporation, such as m¹G in E. coli 23S rRNA [119], N⁶, N⁶- dimethyladenosine in yeast 18S rRNA [120], 2′-O-methylation in mouse piwi-interacting RNA (piRNA) and in yeast and human rRNA [121].

Chemical reagents that can selectively react with certain modification can be used to label the modification and then enhance the signature of RT-mediated arrest or misincorporation. In this respect, CMC (N-cyclohexyl-N′-β-(4-methylmorpholinium) ethylcarbodiimide p-tosylate) was utilized to modify ψ, U and G [122]. After labeling, only CMC-ψ complex is stable, which therefore induced the RT arrest at ψ. Application of this method enabled the identification of four ψ sites in E. coli 23S RNA.

2.5. Single molecule real-time detection

The innovative third-generation sequencing technologies of single molecule real-time detection (SMRT) can realize the direct location analysis of m⁵C and m⁶A in DNA [123-125]. In SMRT sequencing, DNA polymerases catalyze the incorporation of nucleotides into DNA. The discrimination of different nucleobase in DNA can be obtained by the unique kinetic signatures of each nucleotide. In addition to m⁵C and m⁶A, chemical labeling-mediated SMRT sequencing also achieved the mapping of hm⁵C in genomic DNA of mouse embryonic stem cells at single-base resolution [126].

In addition to DNA, modifications in RNA also can be mapped by SMRT sequencing [127]. The kinetics of reverse transcriptase that synthesizes cDNA can be monitored to obtain the sequence content and sites of modifications in RNA. The kinetic signature of reverse transcriptase on the m⁶A is obviously different with other nucleobases, which was then successfully used to locate m⁶A [127]. However, the location analysis of modifications in RNA by SMRT currently has only been achieved in synthetic RNA, and great effort is still needed to realize the location analysis of modifications in cellular RNA by SMRT.

2.6. Nanopore

Nanopore sequencing has emerged as a promising approach in location analysis of modifications in nucleic acids. The characterized ionic current while molecules pass through the nanoscale pore can be measured and utilized to differentiate modified nucleosides from normal nucleosides [128]. Nanopore sequencing is able to readily distinguish m⁵C from cytosine in DNA [129, 130]. m⁵C and hm⁵C also can be easily discriminated in DNA molecules by nanopore sequencing [131]. Moreover, chemical labeling of hm⁵C can improve the performance of nanopore sequencing for locating hm⁵C [132, 133]. Modifications in RNA, such as m⁶A, m⁵C, and inosine also can be distinguished by nanopore sequencing (Fig. 9) [134]. The nanopore sequencing shows to be promising and is expected to have a transformative impact on nucleic acid modifications study.

2.7. Cryo-electron microscopy

The fast advancement of cryo-electron microscopy technology makes it become an important tool for deciphering the structures of biological specimens [135]. Recently, over 130 chemical modifications were visualized in rRNA of HeLa cells by cryoelectron microscopy analysis (Fig. 10) [136]. The high resolution of approximately 2.5Å allowed the distinct determination of the chemical structures of the rRNA. In addition, cryo-electron microscopy can provide detailed structures of the modifications in rRNA, which led to the discovery of many new types and sites of modifications in rRNA. It is believed that the powerful cryoelectron microscopy may play increasing important role and promote the in-depth investigation of nucleic acid modifications.

2.8. Nanomaterials

Ligase chain reaction in combination with gold nanoparticle (AuNP)-based colorimetric assay was developed for the location analysis of m⁵C in DNA at one-base resolution [137]. In this strategy, DNA was treated with bisulfite, where the m⁵C in DNA remained unchanged and cytosine in DNA was converted to uracil. Then the bisulfite-treated DNA was used as the template for the ligation of the two pieces of complementary strands. The ligation position was designed at the site of m⁵C. Thus, the ligated products can be detected by subsequent ligase chain reaction and AuNPs aggregation induced red-to-blue color change, which allowed the visual detection of m⁵C in DNA with the naked eye. In addition to AuNPs, the quantum dot nanomaterials also was employed in conjugation with ligase chain reaction for the location analysis of m⁵C in DNA [138]. This kind of nanomaterial-based detection generally can determine m⁵C in DNA at single site with high detection limit.

3. Conclusions and perspectives

Innovations in technologies and analytical methods have greatly promoted the location study of nucleic acid modifications. In this review, we discussed the principles, advantages and drawbacks of the analytical strategies for mapping nucleic acid modifications and we summarized the important features of these methods in Table 1. Although MS exhibits particular capability on the discovery of new modifications, MS-based analysis may continue to play an important role in the location analysis of nucleic acid modifications. The fast advances of NGS-based methods are becoming the routine tools in study of nucleic acids. Further investigation of a variety of chemical reactions to be coupled with NGS is a promising direction for mapping modifications in nucleic acids. Many chemical reactions can be explored and new reagents can be synthesized to selectively label specific modifications in DNA or RNA.

Innovative third-generation sequencing technologies show extraordinary promising on location analysis of modifications. SMRT and nanopore sequencing provide the useful platforms for direct mapping modifications occurring in nucleic acids without additional chemical or enzymatic treatment. Nanopore sequencing has been well applied in location study of DNA modifications; while the application of third-generation sequencing technologies for locating modifications in RNA from real biological samples may still take some time.

Identification of modifications and mapping their sites in nucleic acids remain challenging. The advancement of technologies and new methods may lead to the discovery of novel modifications. By understanding the property of each analytical method, researchers can select the most suitable strategy for mapping these modifications in nucleic acids. And we expect that new advancement of techniques will promote the mapping study of modifications in nucleic acids and eventually improve our understanding of these modifications.

Acknowledgments

The authors thank the National Key R&D Program of China (No.2017YFC0906800) and the National Natural Science Foundation of China (Nos. 21672166, 21635006, 21721005, 21728802) for the financial support.