Clinical trial NCT04533282: Epitranscriptomic Blood Biomarkers for Coronary Artery Disease - A Prospective Cohort Study (IHD-EPITRAN)
Official title: Epitranscriptomic Biomarkers for Ischemic Heart Disease (IHD-EPITRAN) - A Prospective Cohort Study
Study type: OBSERVATIONAL
Overall status: RECRUITING
Conditions: Ischemic Heart Disease
Coronary Artery Disease
Aortic Valve Stenosis
Acute Myocardial Infarction
Keywords: Ischemic heart disease
Coronary artery disease
Epitranscriptomics
Blood Biomarkers
Biomarkers
Direct long-read nanopore sequencing
RNA editing
Lead sponsor: Hospital District of Helsinki and Uusimaa
Sponsor class: OTHER
Collaborators: name: University of Helsinki
name: Tays Heart Hospital
name: Middle East Technical University
name: Folkhälsan Researech Center
name: Karolinska Institutet
name: University of Tartu
First posted: 2020-08-31
Enrollment: 200
Sex: ALL
Minimum age (years): 18
Maximum age (years): 80
Standard ages: ADULT
OLDER_ADULT
Accepts healthy volunteers: no
Observational model: COHORT
Time perspective: PROSPECTIVE
Sampling method: NON_PROBABILITY_SAMPLE
Study population: Cohort I; Acute IHD (STEMI+PCI):

Patients presenting with ST-elevation myocardial infarction to the CCU and admitted for PCI Intervention.

Cohort II; Chronic IHD (elective CABG):

Stable IHD patients with angina pectoris or dyspnea destined to undergo an elective CABG operation based on preceding angiography.

Cohort III; Aortic valve stenosis (elective AVR):

Stable patients with angina pectoris or dyspnea destined to undergo an elective open surgery AVR due to aortic valve calcification and stenosis without recorded comorbid IHD in preceding angiography.

Cohort IV; IHD-negative controls (coronary CT):

Patients referred to Heart Unit's coronary CT angiogram in order to investigate the possibility of IHD for symptoms such as pressing chest pain or dyspnea provoked by exertion with negative results.

Brief summary: Despite advancements in medical care, ischemic heart disease (IHD) remains the leading global cause of death. IHD develops through lipid accumulation into the coronary arteries with subsequent formation of larger atherogenic plaques. During myocardial infarction (MI), a plaque ruptures and subsequent occlusion leads to a death of the heart muscle. The tissue is rapidly replaced with a scar, which may later lead to heart failure (HF).

Optimally, disease biomarkers are analyzed from blood, provide insight into the disease progression and aid the evaluation of therapy efficacy. Unfortunately, no optimal biomarkers have been identified for IHD. The vast but uncounted number of patients with undiagnosed IHD, benefitting from an early diagnosis, underscore the dire need for an IHD biomarker.

Epitranscriptomics, the study of posttranscriptional modifications on RNA, has recently been properly re-established. This expanding field is uncovering a new layer of regulation, controlling processes ranging from cell division to cell death.

Over 170 modiﬁcations have been identiﬁed as posttranscriptional marks in RNA species. These modiﬁcations influence RNA metabolism, including export, stability, and translation. One the most common and intensively studied RNA modiﬁcation is the N6-methyladenosine (m6A), the abundance and effects of which are determined by the interplay between its writers, readers and erasers.

Recent findings suggest a local dysregulation of the m6A dynamics in the myocardium, coalescing in signalling pathway and contractility related RNA transcripts during hypertrophy, MI and HF. While these early reports have focused on the myocardium, the role of the m6A in the circulation during IHD remains unexplored.

We hypothesize the IHD pathophysiology to be reflected in the epitranscriptome of the circulating RNA.

The objective of the IHD-EPITRAN is to identify new IHD biomarkers via cohort comparison of the blood epitranscriptomes from patients with: (1) MI related with coronary angioplasty, (2) IHD treated with elective coronary artery bypass grafting, (3) aortic valve stenosis treated with valve replacement and (4) IHD-healthy controls verified with computerized tomography imaging. The RNA fractionation is followed by the quantitative modifications analysis with mass spectrometry. Ultimately, nanopore RNA sequencing with simultaneous m6A identification in their native sequences is carried out using recently published artificial intelligence-based algorithm.

Detailed description: Background - Ischemic Heart Disease

With a global prevalence of 126.5 million people and a yearly mortality of 8.9 million, ischemic heart disease (IHD) is the leading cause of death (GBD Collaborators 2017). IHD develops as a result of an ongoing atherogenesis, a process of lipid buildup into the walls of the coronary arteries. Eventually, the accumulation of cholesterol and calcium form deposits, i.e. plaques, that narrow down the vessel lumina. At the onset of myocardial infarction (MI), such a plaque ruptures, which leads to anoxic death of the myocardium relying on the blood supply of the vessel. In the short term, the heart's function is rescued by the formation of a stiff scar. Over time, the pumping function typically deteriorates with ensuing heart failure (HF).

Optimally, disease biomarkers are analyzed from blood and provide both insight into disease progression, severity and aid the evaluation of therapy efficacy. Unfortunately, optimal biomarkers for IHD have not been identified to date. The vast but uncounted number of those living with undiagnosed IHD, benefitting from an early diagnosis, underscore the need for an IHD biomarker.

Background - Epitranscriptomics

Epitranscriptomics, the field of the posttranscriptional modifications of RNA, analogous to DNA epigenetics, has recently gained more scientific remark (Saletore 2012). This expanding field is uncovering a new layer of biological regulation that controls processes ranging from cell proliferation to death.

Although more than 170 RNA modifications have been identified in RNA species (Yang 2018), adenosine methylation at the nitrogen-6 position (m6A) stands out as one the most common and intensively studied (Liu 2020). Following the initial findings of m6A methylated mRNAs in 1970s (Desrosiers 1974), epitranscriptomics has only reignited during the last decade with the development of reliable in vivo methods to assess m6A (Dominissini 2012). These early immunoprecipitation-based methodologies (i.e. meRIP-seq) facilitated the discoveries of the specific m6A punctuation patterns and the specifically denoted writers (methyltransferases; METTL3/14/WTAP1 complex and METTL16), readers (binding proteins; e.g. YTH-protein family), and erasers (demethylases; ALKBH5; Zaccara 2019).

Multiple readers mediating the m6A effects from the YTH domain family have been identified. YTHDF1-3 promote mRNA translation, degradation or both, respectively. YTHDC1-2, eIF2, METTL3 and ribosomes recognize m6A as well and promote translation and cytoplasmic compartmentalization (Zaccara 2019).

Changes in m6A levels have been associated with numerous pathologies, including cancers, cardiovascular and neurological disorders (Frye 2016). Studies using both knockout and overexpression of the m6A writers and erasers have revealed their role in driving immune reactivity, cell proliferation, migration and apoptosis (Delaunay 2019).

Background - Epitranscriptomics \& heart

A recently published set of studies suggest a role for both the m6A and A-to-I editing (adenosine-to-inosine, another common modification diversifying the transcriptome by altering base pairing) during hypertension (Jain 2018), neovascularization (Kwast 2019), myocardial hypertrophy (Dorn 2019, Kmietczyk 2019), ischemia (Mathiyalagan 2019) and HF (Berulava 2020). While a body of evidence indicate epitranscriptomics as a contributor in cardiac health and disease, all of the initial studies have focused on the myocardium.

Rationale, objectives and significance

We hypothesize that the pathophysiology of IHD is reflected in the epitranscriptome of the circulating RNA. For example, ribosomal RNA has a long half-life and its modifications in the coronary circulation can respond to the state of the heart.

The objectives of the IHD-ERPITRAN project are: (1) to establish screening protocols for epitranscriptomic biomarkers, (2) provide epitranscriptomic insight into the IHD pathophysiology, (3) identify a set of IHD biomarker candidates and (4) open avenues for therapeutic development that were previously disregarded due to lack of research and methodological limitations. For example, our research group has reported a discovery of small molecule ligands for RNA methyltransferases (Selberg 2019).

Since epitranscriptomics has been suggested to play a regulatory role in heart and the IHD-EPITRAN is among the first to evaluate epitranscriptomic blood biomarkers in IHD, the project holds keys for a scientific breakthrough with expectable global clinical benefits.

Methods

Study cohorts - 200 patients are recruited into four cohorts, counting 50 patients each. First cohort consists of patients presenting to Heart Unit with an acute MI revascularized with percutaneous coronary intervention (PCI), a cohort that gives valuable insight to the way the acute IHD phenotype is mirrored to the blood epitranscriptomes. Second, the main study cohort of the project will be formed from patients with stable IHD phenotype treated with coronary artery bypass grafting (CABG) and offers critical information in respect to chronic ischemia and coronary atherogenesis, both ongoing processes in the silent, i.e. undiagnosed, IHD. Third, patients destined for an aortic valve replacement (AVR) therapy due to stenosis with no IHD provide the positive control group with a differing cardiac pathology, and, fourth, together with the healthy non-IHD patients verified with the coronary computed tomography (CT), form the control cohorts of the project.

Study samples - The IHD-EPITRAN uncovers blood cellular and cell-free RNA epitranscriptomes from the twice collected TEMPUS(TM) and EDTA blood samples, respectively. In addition, for MI and CABG cohorts, the epitranscriptomes from the right atrial appendage tissue samples are deciphered for a reference point. Finally, a set of supplementary HF and IHD biomarkers: NT-proBNP, hsCRP, sST2 and TMAO (Aimo 2019, Ahmadmehrabi 2017 and Tibaut 2019) are measured.

RNA methods - The RNA isolation and fractionation with in-lab-validated methods is followed by the quantitative analysis of seven base modifications utilizing an UHPLC-triple quadrupole liquid chromatography-tandem mass spectrometry (LC-MS/MS) system (Selberg 2019). Following the quantitative analysis, anti-m6A antibody-based RNA pre-selection rich in m6A marks in tandem with sequencing, i.e. meRIP-seq, is performed. Ultimately, after the UHPLC-LC-MS/MS and meRIP-seq, direct nanopore long-read sequencing enabling identification of the m6A marks in their native sequences is carried out using the novel and recently developed algorithm by Liu et al. (2019). The algorithm bases on the recognition of the characteristic disruptions of electrical current by the m6A flowing through specific protein nanopores.

DNA sequencing - In order to identify the A-to-I RNA editing events, a whole genome next generation sequencing is performed for all study participants to enable a matched DNA-to-RNA base comparison (Park et al., 2012).

Bioinformatics - In order to interpret the massive datasets yielded from the sequencing, bioinformatics protocols are developed in collaboration with the bioinformatics teams from Folkhälsan (Karolinska Institutet, Dr. Shintaro Katayma) and Middle East Technical University (professor Nurcan Tuncbag). The sequences are evaluated against known transcript libraries.

Clinical methods - The echocardiography gives insight into the functional state of the heart and is performed with pre-specified parameters both during the initial hospitalization and 3 months after either PCI, CABG or AVR. These cardiologists' appointments include the severity assessment of angina pectoris and exertional dyspnea with CCS and NYHA gradings, respectively. In addition, morbidity self-assessment with the standardized Short Form 36 Health Survey is included. Moreover, a patient record system search including hospital admissions, deaths and medication changes associable to cardiovascular health are collected 6 months prior and after operation.

Power analysis - The cohort sizes were determined with RNAseqPS web tool designed for evaluating statistical power for RNAseq experiments (Guo et al., 2014). Parameter values used (false discovery rate \[FDR\] \<0.05, total number of genes for testing 20 000, predicted prognostic genes 1500, minimum fold change threshold 2 for differential expression and average read counts of 10 for prognostic genes) were derived from an applicable report concerning epitranscriptomics during HF (Berulava et al., 2020). Prognostic gene dispersion value of 0,215 was applied based on the article focusing on the issue by Yoon \& Nam 2017, which analysed dispersions from ten publicly available RNA-seq datasets, four of which reflected unrelated replicates applicable for the IHD-EPITRAN with the dispersions ranging between 0.15 and 0.28.

Based on the power analysis described above, n=25 per cohort provides sufficient power (P≥0,95). To fluently tackle possible preanalytical errors, perform subsequent validation and follow-up analyses, the cohorts were doubled.

Patient registries

Individual participant identification information data collected during the IHD-EPITRAN will be mainly stored in electrical form on HUS's and Tays's network hard drives that are protected by role-based access rights. Two separate registries are created for the IHD-EPITRAN: (1) A Key Registry, which contains all participant identification information and links between pseudonymization codes that are used in the (2) Research Registry (fully pseudonymized) that contains all other collected study information from the participants. Antti Vento, the director of the IHD-EPITRAN, has the accession rights to the Key Registry. The registries are stored for 10 years to enable possible later follow-up projects.

Processing description documents for both registries, alongside self assessment document for the risks of the study, have been created for the IHD-EPITRAN in Finnish. These documents has been reviewed and approved by the Ethics Board of the Hospital District of Helsinki and Uusimaa (Dnr. HUS/1211/2020).

Research group

Docent Antti Vento is the director of both the IHD-EPITRAN project and the HUS Heart and Lung Center. The interdisciplinary team composed of clinicians, basic researchers and analytical experts enables the IHD-EPITRAN to synergistically perform recruitment, sample collection and preparation followed by the bioinformatic analyses of the epitranscriptomic landscape of circulating RNA in IHD. Docents Mika Laine, Pasi Karjalainen, Helena Rajala and Satu Suihko, all cardiologists from HUS heart station, will carry out the CT patients' recruitment, imaging analyses, treatment of the MI patients and, thirdly, MI, CABG and AVR patients' control appointments. Heart and vascular surgeon Kari Teittinen from the HUS Heart and Lung Center, together with the cardiothoracic surgeon Jahangir Khan PhD and professor Jari Laurikka, both from Tays Heart Hospital, carry out the recruitment of the CABG and AVR patients and their operations. Research nurses Kati Oksaharju and Kati Helleharju, from HUS Heart and Lung Center and Tays Heart Hospital, respectively, coordinate patient recruitment and, in close collaboration with laborant Lahja Eurajoki, the sample collection, transport and prompt storage after freezing. The project will be carried out in the laboratory of docent Esko Kankuri and professor Eero Mervaala (University of Helsinki, CardioReg Group, Department of Pharmacology). B.M. Vilbert Sikorski and MSc Daria Blokhina, will collaboratively perform RNA isolation, fractionation, UHPLC-LC-MS/MS analysis and protocol article formulation. Additionally, Sikorski has collaboratively planned study protocol with the other IHD-EPITRAN co-authors, applied study permissions and continues to coordinate information exchange between collaborators and carries out results article writing with the co-authors.

Ethics

This study will be conducted in compliance with the protocol approved by the HUS Ethics Board (Dnr. HUS/1211/2020). Respective regional and national boards will evaluate the protocol subsequently for the multicenter sample collection. Care will be taken so that each study participant has enough time to familiarize himself/herself with the study protocol. Thereafter, an informed consent discussion will be held and written informed consent obtained from the participant.

Eligibility criteria: Inclusion Criteria:

1. Cohort I, STEMI + PCI:

   1. Earlier PCIs and silent infarctions eligible.
   2. ECG confirmed STEMI with Troponin I elevation and pressing chest pain.
   3. ECG-indicated local damage correlates with recorded dyskinesia in TTE.
   4. During acute PCI and angiography, only one clear occlusion.
   5. Successful initial coronary artery reperfusion during PCI.
2. Cohort II, Chronic IHD + elective CABG:

   1. Chronic and either CCS or NYHA II-IV symptoms for at least one month.
   2. First and elective operation. Only heart operation to be performed.
   3. In transthoracic echocardiogram (TTE):

      * No indication of cardiomyopathy other than ischemic.
      * No pathological remodelling (valves, ventricles and atrias).
      * No clear indication of significant heart failure (i.e. LVEF \> 25%)
3. Cohort III, elective aortic replacement therapy (AVR) for stenosis:

   1. Chronic and either CCS or NYHA II-IV symptoms for at least one month.
   2. Operated as an open heart surgery (either prosthetic or biovalves)
   3. No signs of IHD in coronary angiography.
   4. Both bicuspid and tricuspid valves eligible.
4. Cohort IV, IHD-negative healthy controls defined by coronary CT:

   1. Computerized tomography angiogram results are categorized as negative for coronary artery disease.
   2. No known heart disease.

Exclusion Criteria:

* Condition that limits life expectancy.
* Combination procedures (i.e. CABG+valve).
* Chronic renal insufficiency (KDIGO scale Pt-GFR \< 45/min).
* Active inflammatory/infectious process.
* Known disease affecting either blood or bone marrow.
* Structural or functional congenital heart disease.
* Recorded atrial fibrillation.
* Other comorbidities in poor clinical control (i.e. uncontrolled severe hypertension \>170-180/100 and for diabetes HbA1c \> 60 mmol/l).
* Insulin treated diabetes.

Primary outcomes:
measure: Blood leukocyte RNA's epitranscriptomic changes specifically attributable for IHD
timeFrame: 2020-2023
Secondary outcomes:
measure: Blood cell-free RNA's epitranscriptomic alterations specifically attributable for IHD.
timeFrame: 2020-2023
Other outcomes:
Locations:
city: Helsinki
country: Finland
facility: Hospital District of Helsinki and Uusimaa, Helsinki University Hospital, Heart and Lung Center & Cardiac Unit
state: Uusimaa
Study officials:
affiliation: Helsinki University Central Hospital, Heart and Lung Center
name: Antti E Vento, Docent
role: STUDY_DIRECTOR
affiliation: University of Helsinki, Faculty of Medicine, Department of Pharmacology
name: Esko Kankuri, Docent
role: PRINCIPAL_INVESTIGATOR
References:
citation: Aimo A, Januzzi JL Jr, Bayes-Genis A, Vergaro G, Sciarrone P, Passino C, Emdin M. Clinical and Prognostic Significance of sST2 in Heart Failure: JACC Review Topic of the Week. J Am Coll Cardiol. 2019 Oct 29;74(17):2193-2203. doi: 10.1016/j.jacc.2019.08.1039.
citation: Ahmadmehrabi S, Tang WHW. Gut microbiome and its role in cardiovascular diseases. Curr Opin Cardiol. 2017 Nov;32(6):761-766. doi: 10.1097/HCO.0000000000000445.
citation: Berulava T, Buchholz E, Elerdashvili V, Pena T, Islam MR, Lbik D, Mohamed BA, Renner A, von Lewinski D, Sacherer M, Bohnsack KE, Bohnsack MT, Jain G, Capece V, Cleve N, Burkhardt S, Hasenfuss G, Fischer A, Toischer K. Changes in m6A RNA methylation contribute to heart failure progression by modulating translation. Eur J Heart Fail. 2020 Jan;22(1):54-66. doi: 10.1002/ejhf.1672. Epub 2019 Dec 17.
citation: Delaunay S, Frye M. RNA modifications regulating cell fate in cancer. Nat Cell Biol. 2019 May;21(5):552-559. doi: 10.1038/s41556-019-0319-0. Epub 2019 May 2.
citation: Desrosiers R, Friderici K, Rottman F. Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells. Proc Natl Acad Sci U S A. 1974 Oct;71(10):3971-5. doi: 10.1073/pnas.71.10.3971.
citation: Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, Cesarkas K, Jacob-Hirsch J, Amariglio N, Kupiec M, Sorek R, Rechavi G. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012 Apr 29;485(7397):201-6. doi: 10.1038/nature11112.
citation: Dorn LE, Lasman L, Chen J, Xu X, Hund TJ, Medvedovic M, Hanna JH, van Berlo JH, Accornero F. The N6-Methyladenosine mRNA Methylase METTL3 Controls Cardiac Homeostasis and Hypertrophy. Circulation. 2019 Jan 22;139(4):533-545. doi: 10.1161/CIRCULATIONAHA.118.036146.
citation: Frye M, Jaffrey SR, Pan T, Rechavi G, Suzuki T. RNA modifications: what have we learned and where are we headed? Nat Rev Genet. 2016 Jun;17(6):365-72. doi: 10.1038/nrg.2016.47. Epub 2016 May 3.
citation: GBD 2017 Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018 Nov 10;392(10159):1789-1858. doi: 10.1016/S0140-6736(18)32279-7. Epub 2018 Nov 8. Erratum In: Lancet. 2019 Jun 22;393(10190):e44. doi: 10.1016/S0140-6736(19)31047-5.
citation: Jain M, Mann TD, Stulic M, Rao SP, Kirsch A, Pullirsch D, Strobl X, Rath C, Reissig L, Moreth K, Klein-Rodewald T, Bekeredjian R, Gailus-Durner V, Fuchs H, Hrabe de Angelis M, Pablik E, Cimatti L, Martin D, Zinnanti J, Graier WF, Sibilia M, Frank S, Levanon EY, Jantsch MF. RNA editing of Filamin A pre-mRNA regulates vascular contraction and diastolic blood pressure. EMBO J. 2018 Oct 1;37(19):e94813. doi: 10.15252/embj.201694813. Epub 2018 Aug 7.
citation: Kmietczyk V, Riechert E, Kalinski L, Boileau E, Malovrh E, Malone B, Gorska A, Hofmann C, Varma E, Jurgensen L, Kamuf-Schenk V, Altmuller J, Tappu R, Busch M, Most P, Katus HA, Dieterich C, Volkers M. m6A-mRNA methylation regulates cardiac gene expression and cellular growth. Life Sci Alliance. 2019 Apr 9;2(2):e201800233. doi: 10.26508/lsa.201800233. Print 2019 Apr.
citation: van der Kwast RVCT, Quax PHA, Nossent AY. An Emerging Role for isomiRs and the microRNA Epitranscriptome in Neovascularization. Cells. 2019 Dec 25;9(1):61. doi: 10.3390/cells9010061.
citation: Liu J, Li K, Cai J, Zhang M, Zhang X, Xiong X, Meng H, Xu X, Huang Z, Peng J, Fan J, Yi C. Landscape and Regulation of m6A and m6Am Methylome across Human and Mouse Tissues. Mol Cell. 2020 Jan 16;77(2):426-440.e6. doi: 10.1016/j.molcel.2019.09.032. Epub 2019 Oct 29.
citation: Liu H, Begik O, Lucas MC, Ramirez JM, Mason CE, Wiener D, Schwartz S, Mattick JS, Smith MA, Novoa EM. Accurate detection of m6A RNA modifications in native RNA sequences. Nat Commun. 2019 Sep 9;10(1):4079. doi: 10.1038/s41467-019-11713-9.
citation: Mathiyalagan P, Adamiak M, Mayourian J, Sassi Y, Liang Y, Agarwal N, Jha D, Zhang S, Kohlbrenner E, Chepurko E, Chen J, Trivieri MG, Singh R, Bouchareb R, Fish K, Ishikawa K, Lebeche D, Hajjar RJ, Sahoo S. FTO-Dependent N6-Methyladenosine Regulates Cardiac Function During Remodeling and Repair. Circulation. 2019 Jan 22;139(4):518-532. doi: 10.1161/CIRCULATIONAHA.118.033794.
citation: Saletore Y, Meyer K, Korlach J, Vilfan ID, Jaffrey S, Mason CE. The birth of the Epitranscriptome: deciphering the function of RNA modifications. Genome Biol. 2012 Oct 31;13(10):175. doi: 10.1186/gb-2012-13-10-175.
citation: Selberg S, Blokhina D, Aatonen M, Koivisto P, Siltanen A, Mervaala E, Kankuri E, Karelson M. Discovery of Small Molecules that Activate RNA Methylation through Cooperative Binding to the METTL3-14-WTAP Complex Active Site. Cell Rep. 2019 Mar 26;26(13):3762-3771.e5. doi: 10.1016/j.celrep.2019.02.100.
citation: Tibaut M, Caprnda M, Kubatka P, Sinkovic A, Valentova V, Filipova S, Gazdikova K, Gaspar L, Mozos I, Egom EE, Rodrigo L, Kruzliak P, Petrovic D. Markers of Atherosclerosis: Part 1 - Serological Markers. Heart Lung Circ. 2019 May;28(5):667-677. doi: 10.1016/j.hlc.2018.06.1057. Epub 2018 Oct 4.
citation: Zaccara S, Ries RJ, Jaffrey SR. Reading, writing and erasing mRNA methylation. Nat Rev Mol Cell Biol. 2019 Oct;20(10):608-624. doi: 10.1038/s41580-019-0168-5. Epub 2019 Sep 13.
citation: Park E, Williams B, Wold BJ, Mortazavi A. RNA editing in the human ENCODE RNA-seq data. Genome Res. 2012 Sep;22(9):1626-33. doi: 10.1101/gr.134957.111.