REVIEW article

Front. Cardiovasc. Med., 19 February 2019

Sec. Cardiovascular Genetics and Systems Medicine

Volume 6 - 2019 | https://doi.org/10.3389/fcvm.2019.00009

Multi-Omics Approaches to Study Long Non-coding RNA Function in Atherosclerosis

  • 1. Center for Public Health Genomics, University of Virginia, Charlottesville, VA, United States

  • 2. Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, VA, United States

  • 3. School of Nursing, University of Virginia, Charlottesville, VA, United States

  • 4. Data Science Institute, University of Virginia, Charlottesville, VA, United States

  • 5. Department of Biomedical Engineering, University of Virginia, Charlottesville, VA, United States

  • 6. Department of Public Health Sciences, University of Virginia, Charlottesville, VA, United States

Abstract

Atherosclerosis is a complex inflammatory disease of the vessel wall involving the interplay of multiple cell types including vascular smooth muscle cells, endothelial cells, and macrophages. Large-scale genome-wide association studies (GWAS) and the advancement of next generation sequencing technologies have rapidly expanded the number of long non-coding RNA (lncRNA) transcripts predicted to play critical roles in the pathogenesis of the disease. In this review, we highlight several lncRNAs whose functional role in atherosclerosis is well-documented through traditional biochemical approaches as well as those identified through RNA-sequencing and other high-throughput assays. We describe novel genomics approaches to study both evolutionarily conserved and divergent lncRNA functions and interactions with DNA, RNA, and proteins. We also highlight assays to resolve the complex spatial and temporal regulation of lncRNAs. Finally, we summarize the latest suite of computational tools designed to improve genomic and functional annotation of these transcripts in the human genome. Deep characterization of lncRNAs is fundamental to unravel coronary atherosclerosis and other cardiovascular diseases, as these regulatory molecules represent a new class of potential therapeutic targets and/or diagnostic markers to mitigate both genetic and environmental risk factors.

Introduction

Despite intensive research into the underlying pathogenesis of atherosclerosis/coronary artery disease (CAD), this disease still remains a significant public health burden. Atherosclerosis is a complex disease involving both environmental and genetic risk factors resulting in plaque formation and inflammation in the vessel wall. A number of cellular responses have been proposed to contribute to disease progression, including endothelial cell dysfunction, vascular smooth muscle cell phenotypic switching, macrophage/foam cell activation/necrosis/defective efferocytosis, and defective lipid/lipoprotein metabolism (1, 2). Much of the research to date has focused on the role of protein coding genes in atherosclerosis leading to the identification of a number of proteins described as key drivers. However, our understanding of the causal mechanisms of this disease remains limited, likely due to our incomplete functional knowledge of the non-coding genome. It is now well-established that more than 90% of disease-associated variants reside in non-coding regions, once considered “junk sequences” (3). An increasing number of studies, especially those utilizing high-throughput sequencing technologies, have shown that a number of non-coding RNAs (ncRNAs) are differentially regulated during disease, supporting potential functional roles for these molecules (4).

The formation of the GENCODE project (part of the ENCODE project) and other large-scale initiatives such as NONCODE have revealed 75% of the genome is transcribed, yet only 2% encodes for protein, suggesting alternative functional roles for ncRNA transcripts (5, 6). Since the development of RNA-seq and other high-throughput sequencing assays, ncRNAs are now appreciated as key regulators of gene expression (7). A number of the essential types of non-coding RNAs are already well-characterized such as transfer, spliceosomal, and ribosomal RNAs. Aside from these housekeeping RNA elements, the remaining types of non-coding RNAs are subdivided into two classes, small and long non-coding RNAs, based on their size. Examples of small non-coding RNAs include microRNAs (miRNA), small nucleolar RNAs (snoRNA), and PIWI-interacting RNAs (piRNA) (8). These elements can act as positive or negative regulators of gene expression and generally exert their influence through complementary base pairing to their target transcript 3′ or 5′ untranslated regions (9).

Long non-coding RNAs (lncRNAs) represent a heterogeneous class of non-coding RNAs that includes transcripts >200 nucleotides, which lack functional protein coding ability (10). Within this lncRNA class, they are also classified based on their genomic location and broadly encompass enhancer-related RNAs (eRNAs), transcribed ultraconserved RNAs, intronic RNAs, long intergenic RNAs (lincRNAs), and natural antisense transcripts (NATs) (10). Contrary to canonical linear lncRNAs, a distinct group of lncRNAs are known as circular RNAs (circRNAs) due to their circular structure, which often results from backsplicing (11). By participating in both transcriptional and post-transcriptional stages, lncRNAs modulate gene expression through multiple distinct mechanisms. Further insight into these regulatory mechanisms will facilitate a better understanding of disease biology and identify additional viable targets for therapeutic intervention or diagnostics. Here, we present an overview of various lncRNAs relevant to atherosclerosis and highlight next-generation sequencing approaches to systematically investigate lncRNA function, as well as the ongoing challenges in this exciting field.

Mechanisms of Long Non-coding RNA Function

LncRNAs are a heterogeneous class of ncRNAs (>200 nucleotides in length) that do not contain a functional open reading frame (12). LncRNAs can be encoded on either the sense or antisense DNA strand and may be located within a protein coding gene or in the intergenic regions (12). Similar to mRNAs, lncRNAs are transcribed by RNA polymerase II. Many transcripts are polyadenylated, multi-exonic, undergo RNA splicing, and contain a 5' cap. Often, their active promoters are marked with H3K4me3 and gene bodies have H3K36me3 histone modifications (13). Unlike protein coding genes, lncRNAs are not translated into protein and thereby lack functional initiation and termination codons (8). They are expressed at a much lower levels relative to their protein coding counterparts and lack robust evolutionary conservation. Despite this low level of conservation, the expression pattern of lncRNAs has been shown to be relatively cell/tissue specific (1416).

The mechanism of action of these regulatory factors were categorized into four broad groups proposed by Wang and Chang (17). Signaling lncRNAs represent a class which exhibit a high degree of spatial and temporal specificity that serve a role in signal transduction. Once transcribed, these signaling lncRNAs have effector functions in activating appropriate downstream pathways in response to a stimulus. Additionally, their presence may indicate a particular developmental cell state, condition, or overall transcriptional activity (17). Another mechanism in which lncRNAs exert their regulatory function is by acting as decoy molecules to limit the availability of RNA binding factors to interact with their partners. By impairing the ability of chromatin remodelers, transcription factors, and miRNA from binding to their target genes, decoy lncRNAs can inhibit downstream effector functions (17). It is noteworthy that miRNAs also have the ability to target lncRNAs directly and thereby influence transcriptional regulation and vascular functions (1820). Aided by their ability to bind protein as well as base pair with target sequences, guide lncRNAs are responsible for localizing transcriptional regulators to specific regions. Similar to the role of guide lncRNAs, scaffold lncRNAs use their protein binding ability to provide a surface to mediate protein-protein interactions (17). In these various ways, lncRNAs represent a distinct class of regulatory elements to modulate transcriptional activities.

Overview of LncRNAs Involved in Atherosclerosis

Functional Studies of LncRNAs in Atherosclerosis

Depending on the cell types involved, lncRNAs play myriad roles in diverse atherosclerotic processes in the vessel wall including cell proliferation, migration, differentiation, apoptosis, and inflammation. They also play important roles in the regulation of cholesterol and lipid metabolism. Pertinent cell types include smooth muscle cells, endothelial cells, macrophages, and hepatocytes (Figure 1). A more comprehensive overview of key lncRNAs in atherogenic processes is given in Table 1. Perhaps the most well-studied in atherosclerosis is CDKN2B-AS1 or ANRIL (Antisense Non-coding RNA in the INK4 Locus), which acts in several cell types relevant to CAD (21, 22, 81). ANRIL acts as a guide lncRNA to localize polycomb repressive complex (PRC) at target promoters through a direct interaction with its subunits, CBX7 or SUZ12. PRC then adds H3K27me3 modifications to this region to repress transcription (22). Loss of function studies have suggested ANRIL acts in cis in order to regulate transcription of the nearby tumor suppressors, CDKN2A and CDKN2B (82, 83). Consistent with this mechanism of action, ANRIL expression correlates with a more proliferative phenotype in endothelial cells and vascular smooth muscle cells (VSMC) (22, 84, 85). In addition to acting in cis, ANRIL also acts in trans (via Alu elements) to regulate other genes that participate in proatherogenic pathways (22). Since ANRIL is not well-conserved in mice, in vivo functional studies have been challenging (86). A more complete overview of the atherogenic roles of ANRIL RNA species has been recently documented in a review by Holdt and Teupser (87).

Figure 1

Table 1

NameFull nameInitial detection method(s)Cell types/tissues relevant to CADFunction relevant to CADReferences
ANRILA new large antisense non-coding RNALong range PCR and nucleotide sequencingSMC, EC, MacCell cycle regulation
Regulation of apoptosis
(2124)
GAS5Growth arrest specific 5Characterization of cDNA librarySMC, EC, MacRegulates apoptosis
Regulates autophagy
Regulates smooth muscle differentiation
(2528)
HIF1a-AS1Hypoxia-inducible factor 1-alpha Antisense RNA 1RIP, qRT-PCRSMCRegulation of VSMC apoptosis(2932)
LincRNA-p21Long intergenic non-coding RNA at p21 locusNimblegen lincRNA tiling microarray platformSMC, MacRegulation of cell proliferation
Regulation of apoptosis
(33, 34)
MALAT1Metastasis associated lung adenocarcinoma transcript 1Identified in numerous physiological processesSMC, EC, MacRegulation of cell proliferation
Regulation of inflammation
Regulation of autophagy
(3540)
MEG3Maternally expressed 3Characterization of cDNA libraryEC, SMCRegulates endothelial cell proliferation
Inhibits angiogenesis
Regulates smooth muscle proliferation/migration
(4143)
NEAT1Nuclear paraspeckle assembly transcript 1Affymetrix expression arraySMC, MacRegulation of VSMC phenotypic switching
Macrophage polarization
(4447)
NEXN-AS1Nexilin antisense transcript 1lncRNA arrayEC, SMCRegulation of adhesion molecules/monocyte recruitment
Regulation of inflammatory cytokines
(48)
CHROMECholesterol homeostasis regulator of miRNA expressionCharacterization of transcript proximal to CAD and plasma HDL-C associated locusMac, LiverRegulation of cholesterol homeostasis (49)
AK098656AK098656lncRNA arraySMCRegulation of VSMC phenotypic switching(50)
CASC15Cancer susceptibility 15Quantification of copy number gains in metastatic melanomaSMCVSMC stiffness(51, 52)
H19H19cDNA characterizationSMCRegulation of VSMC differentiation
Let 7 miRNA sponge
(5355)
MIR122HG/ Lnc-Ang362 MicroRNA 122 Host GeneRNA-seqSMCRegulation of VSMC proliferation(56)
MYOSLIDMYOcardin-induced Smooth muscle LncRNA, Inducer of DifferentiationRNA-seqSMCRegulation of SMC differentiation(57)
PACERP50-associated COX-2 extragenic RNAANalysis of ChIP data, RT-qPCRSMCRegulation of COX-2 expression
VSMC stiffness
(51, 58)
SENCRSmooth muscle and Endothelial cell-enriched migration/differentiation-associated long Non-coding RNARNA-seq SMCRegulation of myocardin, SMC contractile gene program(59)
SMILRSmooth muscle-induced lncRNA enhances replicationRNA-seqSMCRegulation of VSMC proliferation(60)
CLDN10-AS1Claudin 10 antisense transcript 1lncRNA microarrayECRegulation of endothelial signaling(61)
CTC-459I6.1RASGRF2 antisense RNA 1lncRNA microarrayECRegulation of endothelial signaling(61)
GATA6-ASGATA6 antisenseRNA-seqECRegulation of endothelial signaling
Effects on endothelial- mesenchymal transition
(62)
LEENElncRNA that enhances eNOS expressionRNA-seq, chromosome conformation captureECRegulation of eNOS and endothelial function (63)
MIATMyocardial infarction associated transcriptIsolation of specific cDNA, RACEECRegulation of angiogenesis(64, 65)
NRONNon-protein coding RNA, repressor of NFATScreen of cDNA library to identify conserved lncRNAsshRNA screen to identify NFAT regulatorsECRegulation of angiogenesis
Regulation of vascular development
(6668)
sONE/ ATG9BAutophagy related 9BcDNA characterizationECRegulation of eNOS(69)
STEELSpliced transcript endothelial-enriched lncRNACustom lncRNA microarrayECAngiogenesis(70)
TIE1-ASEndothelial-specific receptor tyrosine kinase 1 antisenseDetection in EST libraries, RACEECRegulation of vascular development(71, 72)
Dnm3osDynamin-3 opposite strand/antisense RNAIsolation of clones, RACEMacRegulation of macrophage inflammation(73, 74)
lncRNA-Mirt2lncRNA Myocardial infarction associated transcriptlncRNA expression microarrayMacNegative regulator of inflammation(75)
MeXIS Macrophage-expressed LXR-induced sequenceRNA-seqMacRegulation of cholesterol metabolism(76)
APOA1-ASApolipoprotein A1 antisense transcriptIdentification in EST library, RACE LiverRegulation of APOA1 expression(77)
LeXisLiver-expressed LXR-induced sequenceRNA-seqLiverRegulation of cholesterol metabolism(78)
lncLSTRlncRNA Liver-Specific Triglyceride RegulatorCharacterization of imprinted genes
Reannotation of Affymetrix Mouse Genome 430 2.0 Array, search for liver-enriched transcripts
LiverRegulation of triglyceride metabolism(79)
TRIBALTribbles homolog 1-associated locusIdentification in EST library, RACELiverRegulation triglyceride metabolism(80)

List of long non-coding RNAs with functional relevance in coronary artery disease cell types/tissues.

Color scheme: Gray, lncRNAs associated with multiple cell types/tissues. Purple, lncRNAs associated with smooth muscle cells (SMC). Green, lncRNAs associated with endothelial cells (EC), Yellow, lncRNAs associated with macrophages (Mac), Blue, lncRNAs associated with cholesterol metabolism in liver. RACE: Rapid amplification of cDNA ends, EST: Expressed sequence tag.

The ubiquitously expressed and evolutionarily conserved lncRNA, MALAT1, is decreased in atherosclerotic plaques (88, 89), and reduced MALAT1 levels in hematopoietic cells promotes atherosclerosis and inflammation in mice in vivo (89). In contrast, MALAT1 knockdown in VSMC and EC results in cell cycle arrest and reduced proliferation (35, 36). Other lncRNAs (e.g., lincRNA-p21, HIF1α-AS, circular ANRIL, and GAS5) have been implicated in cell death/apoptotic pathways through various mechanisms (23, 29, 33, 90), as described elsewhere (91, 92). In particular, lincRNA-p21 was shown to be reduced in CAD patients and mouse models of atherosclerosis, and regulates p53-dependent smooth muscle cell proliferation and apoptosis (33). Many lncRNAs have established immune and inflammatory roles. For example, heterozygous MALAT1-deficient ApoE-/- mice have increased inflammation and atherosclerosis (93). LncRNA Dnm3os (dynamin 3 opposite strand) is upregulated in diabetic induced macrophages and regulates nucleolin and epigenetic mediated inflammatory responses (73). Finally, in human macrophages treated with oxidized LDL, HOTAIR regulates oxidative stress and inflammation (94).

There are several lncRNAs with regulatory roles in lipid and cholesterol metabolism. CHROME (cholesterol homeostasis regulator of miRNA expression) is a lncRNA upregulated in carotid plaques, which regulates cholesterol homeostasis in primates in liver and macrophages by inhibiting miRNAs, such as miR-33 (49). NEAT1 promotes pro-atherogenic functions in THP-1 human macrophage cells such as increased ox-LDL lipid accumulation and inflammation by serving as a sponge of miR-342-3p target (44). Finally, differential expression of TRIBAL, APOA1-AS, and lncLSTR is linked to defects in lipid metabolic pathways, mainly in the liver. TRIBAL (TRIB1 associated locus) regulates Trib1 mRNA stability through mitogen activated kinase, consistent with Trib1 regulation (80, 95, 96). Increased TRIBAL expression stabilizes Trib1 expression and upregulates fatty acid oxidative pathways (80). Likewise, lncLSTR (liver-specific triglyceride regulator) regulates plasma triglyceride clearance by modulating apolipoprotein C2 (APOC2) levels and lipoprotein lipase activity (79). APOA1-AS regulates cholesterol levels through epigenetic modulation of APOA1, a protein involved in the cholesterol efflux pathway (77).

LncRNAs With Genetic Associations in Atherosclerosis

Genome wide association studies (GWAS) have linked genetic variation at the ANRIL locus (9p21.3) to many complex phenotypes including CAD, stroke, type 2 diabetes and multiple cancers (97). In addition to ANRIL, the 9p21.3 locus encodes three tumor suppressor proteins: CDKN2A, CDKN2B, and MTAP. Despite each being attractive candidates underlying the locus association with various diseases, several studies report CAD risk polymorphisms associated with ANRIL expression (21, 98). However, association studies of 9p21.3 genotype with ANRIL expression remain complex due to the numerous linear and circular ANRIL forms (23).

Other less studied lncRNAs have been identified from genetic studies of CAD or related traits and may play critical roles in atherosclerosis. For instance, genetic variation in the imprinted lncRNA H19, involved in embryonic development (99) and oncogenesis (100), was associated with CAD and ischemic stroke in Chinese populations (101, 102). H19 was initially shown to be re-expressed in smooth muscle cells in human and rodent atherosclerotic plaques (103), and promotes VSMC proliferation by acting as a let-7a miRNA sponge to upregulate cyclin D (104). However, a recent study revealed endothelial cell restricted expression in human atherosclerotic plaques and a role in endothelial cell aging by suppressing STAT3 signaling (105), similar to lncRNA MEG3 (106). Another endothelial cell lncRNA, MIAT (Myocardial Infarction Associated Transcript), was previously associated with myocardial infarction in a large genetic study of a Japanese population (64). MIAT is upregulated in atherosclerosis plaques (88), and regulates microvascular dysfunction by acting as a competing endogenous RNA (65). Another lncRNA associated with CAD through large-scale GWAS is known as TARID (TCF21 antisense RNA inducing promoter demethylation) (107). TARID was identified as an eQTL target gene in human coronary artery smooth muscle cells (108), and molecular studies suggest this lncRNA guides GADD45A mediated DNA demethylation and inactivation of TCF21 (109), a known tumor suppressor and vascular wall transcription factor associated with CAD (110113). Yet, functional studies of TARID both in VSMCs and in vivo are needed to elucidate its potential role in atherosclerosis. With larger GWAS sample sizes, and complementary eQTL colocalization (114), and transcriptome-wide association studies (TWAS) (115), it is anticipated that even more lncRNAs will be identified with genetic association evidence.

Application of Transcriptomics to Identify and Study LncRNAs Relevant to CAD

General Considerations for Transcriptomics Studies of LncRNAs

While traditional methods to profile lncRNA transcriptomes have relied on microarrays or serial analysis of gene expression (SAGE), these approaches have largely been replaced with the decreasing costs and greater output achieved by RNA-seq (116). In general, RNA-seq provides greater sensitivity and specificity to detect a broad range of ncRNA transcripts, novel isoforms, and interactions between ncRNAs (117). Nonetheless there are some important considerations when designing and conducting RNA-seq based lncRNA screening experiments. For instance, given that lncRNAs are approximately 10X less abundant than mRNAs on average, the basal expression of a typical lncRNA is < 5 fragments per kilobase of transcript per million mapped reads (FPKM) (118). Thus, it is highly recommended to obtain deeper sequencing per sample (~100X read depth) than a typical RNA-seq experiment. Also, while up to 50% of lncRNAs appear to be poly-adenylated (119) and would be detected with mRNA library preparation kits, a more comprehensive landscape of lncRNAs, other ncRNAs, including eRNAs, would require total RNA [poly(A) and non-poly(A)], ribosomal RNA depletion methods of purification. Distinguishing lncRNA transcripts from mRNA transcripts from short-read sequencing data remains a challenge, however deeper, paired-end and stranded sequencing should improve identification of lncRNAs (120). Also, careful study design is needed to ensure sufficient power to detect differentially expressed and transcript-specific lncRNAs, when using standard count-based tools (120). Since many lncRNAs are tissue and cell-specific (121123), it is also worth considering the effects of diluting weak signals from bulk populations of cells, as well as specific environmental contexts that may regulate lncRNA transcript levels. Below, we summarize recent findings of RNA-seq based lncRNA discoveries in specific cell types relevant to CAD/atherosclerosis.

Transcriptomics of Vascular Smooth Muscle Cell Function

The first VSMC lncRNA discovered via RNA-seq was Lnc-Ang362 (HG-MIR222), which is upregulated in rat aortic smooth muscle cells upon stimulation with angiotensin II (56). Lnc-Ang362 promotes VSMC proliferation and is the host-transcript for both miR-221 and miR-222. Bell et al. conducted RNA-seq in human coronary artery smooth muscle cells and identified 31 previously unidentified lncRNAs (59). Notably, one of these was Smooth muscle and Endothelial cell-enriched migration/differentiation-associated long Non-coding RNA (SENCR), which is located antisense to the FLI1 gene. SENCR functionally promotes a contractile smooth muscle phenotype and inhibits migration (59). In a follow-up study, RNA-seq was performed in human coronary artery smooth muscle cells to examine the effect of myocardin (MYOCD) overexpression (57). MYOCD is a potent co-factor that binds with serum response factor (SRF) to activate an array of smooth muscle-specific genes that maintain smooth muscle cell differentiation (124128). Over 100 lncRNAs were differentially expressed, one of which was identified as MYOcardin-induced Smooth muscle LncRNA, Inducer of Differentiation (MYOSLID). Functional studies demonstrated that MYOSLID, a direct transcriptional target of MYOCD/SRF, promotes smooth muscle differentiation and inhibits proliferation (57).

Yu et al. used RNA-seq to compare transcriptomes of coronary and aortic smooth muscle cells subjected to both normal and pathological aortic stiffness, a subclinical risk factor for CAD and various aortic diseases (51). Only two of the top 20 ranked differentially expressed lncRNAs have been studied to date: CASC15 and PACER (RP5-973M2.2). These lncRNAs regulate expression of protein-coding genes in cis and PACER activates COX2 expression (52, 58, 129). Analysis of RNA-seq data highlighted the lncRNA MALAT1 as a key regulator of VSMC stiffness-induced proliferation and migration. Although MALAT1 was originally described as an endothelial lncRNA, MALAT1 regulates the phenotyping switching of VSMCs via activation of the autophagy pathway (36). Using RNA-seq in human smooth muscle cells Ballantyne et al. identified over 300 differentially expressed lncRNAs upon platelet-derived growth factor and interleukin-1 alpha stimulation. The novel lncRNA, Smooth Muscle-Induced LncRNA enhances Replication (SMILR) identified from this study enhances smooth muscle cell proliferation and has increased expression in unstable atherosclerotic plaques (60). The lncRNA NEAT1 (nuclear paraspeckle assembly transcript 1) has recently been implicated in promoting the phenotypic switching of VSMCs (45). RNA-seq demonstrated NEAT1 silencing increases the mRNA levels of numerous critical smooth muscle cell marker genes. Finally, to identify lncRNAs key in smooth muscle cell differentiation, Lim et al. combined and queried diverse RNA-seq datasets from Gene Expression Omnibus (GEO). Dozens of lncRNAs with no previous evidence for roles in VSMC differentiation were identified in this analysis that warrant further investigation, either as cis transcriptional regulators or suppressing miRNA function (130).

The development of custom lncRNA arrays has been applied to identify lncRNAs involved in various processes critical in atherosclerosis. One example is a microarray analysis which identified 580 lncRNAs differentially expressed upon exposure of human aortic smooth muscle cells to cyclic mechanical stretch (131). Another example is identification of AK098656, predominantly expressed in VSMCs, also upregulated in hypertensive patients and involved in promoting a synthetic smooth muscle cell phenotype (50).

Transcriptomics of Endothelial Cell Function

Although not all lncRNAs have a poly(A) tail, Michalik et al. performed deep sequencing of poly(A)-selected RNA in human umbilical vein endothelial cells (HUVECs) and found over half of total RNA composed of non-coding RNA, many of which are lncRNAs (35). This study focused on five lncRNAs with high endothelial expression and strong conservation between mice and humans: MALAT1, linc00493, maternally expressed 3 (MEG3), taurine upregulated gene 1 (TUG1), and linc00657. MALAT1 and MEG3 are strongly upregulated in response to hypoxia while linc006757 are TUG1 are moderately upregulated. In regards to angiogenesis, MALAT1 promotes angiogenesis and induces a switch of endothelial cells from a migratory cell phenotype to a proliferative cell phenotype (132). Huang et al. postulated exosomal MALAT1 from oxidized LDL (oxLDL) treated endothelial cells (HUVECs) promotes macrophage polarization toward the M2 phenotype (133). MEG3 was shown to interact with epigenetic modifiers, to inhibit angiogenesis and contribute to age-related endothelial dysfunction (106, 134, 135).

In another study Miao et al. conducted RNA-seq profiling of endothelial cells subjected to both physiological and pathological flow for various time points (63). They identified and characterized LEENE (lncRNA that enhances eNOS expression) as a lncRNA highly correlated with endothelial nitric oxide synthase (eNOS) expression levels, which is downregulated upon pathological flow (63). Several lncRNAs characterized in smooth muscle cells also have functional significance in endothelial cells. For instance, the SMC lncRNA SENCR regulates the differentiation of pluripotent cells into endothelial cells and promotes angiogenesis in HUVECs (136).

Custom lncRNA microarrays have identified endothelial cell enriched lncRNAs (70) and endothelial lncRNAs differentially expressed in response to specific treatments (e.g., oxidized LDL) (61). These arrays have revealed new candidate lncRNAs in atherosclerosis including spliced-transcript endothelial-enriched lncRNA (STEEL) (70), NEXN-AS1 (48), CLDN10-AS1, and CTC-459I6.1 (61). STEEL facilitates the transcriptional stimulation of both eNOS and Kruppel-like factor (KLF2) and in vivo promotes angiogenesis (70). STEEL expression is decreased in HUVECs exposed to “atheroprotective” flow and expression is increased in HUVECs exposed to disturbed “atheroprone” flow.

Transcriptomics of Macrophage Function and Inflammation

In macrophages, LXR activation promotes cholesterol efflux through activation of target genes such as Abca1 during the formation of HDL. To investigate the regulation of LXR-dependent transcription in macrophages, a recent study conducted large-scale transcriptional profiling of mouse peritoneal macrophages in response to the LXR agonist GW3965. LXR activation stimulated transcription of an array of lncRNAs, of which MeXis was among the strongest induced (76). MeXis is well-conserved in mice and was shown to amplify the LXR-dependent expression of Abca1 in vivo and promote cholesterol efflux in macrophages (76). Loss of MeXis in Ldlr−/− mice was shown to accelerate atherosclerosis through impaired Abca1 expression in macrophages and resulted in decreased cholesterol efflux (76). ATAC-seq in peritoneal macrophages demonstrated decreased chromatin accessibility across the Abca1 locus in response to loss of MeXis. Querying the MeXis interactome through mass spectrometry revealed protein interactions with the nuclear receptor coactivator DDX17. Either directly or indirectly through one of its interacting targets, MeXis represents a potential therapeutic target to regulate macrophage cholesterol efflux.

RNA-seq and lncRNA arrays have identified a number of other macrophage lncRNAs that could represent novel CAD targets. Zhang et al. performed deep RNA sequencing of human monocyte-derived macrophages as well as M1 activated (via interferon gamma and lipopolysaccharide stimulation) and M2 activated (via interleukin 4 stimulation) macrophages (137). This study identified 861 previously unannotated lincRNAs, most of which are not syntenic in mouse. Furthermore, the lncRNA expression profile is dramatically shifted upon M1 activation, supporting the inflammatory nature of atherosclerosis. Similarly, 109 unannotated CD14+ monocyte lincRNAs were highlighted upon exposure to inflammatory stress in vivo (138). Other recent array studies highlighted the macrophage lncRNAs Dnm3os amd Mirt2. Dnm3os is upregulated in bone marrow derived macrophages in diabetic mice compared to controls and is higher in monocytes in human type 2 diabetic patients compared to controls (73). Dnm3os alters global histone modifications in macrophages and upregulates various immune-response and inflammatory genes. LncRNA-Mirt2 is strongly induced by LPS, a toll-like receptor 4 (TLR4) ligand where it acts as a negative feedback inflammatory regulator (75).

Transcriptomics of Cholesterol Metabolism and Hepatocyte Function

Liver X receptors (LXRs) are nuclear factor transcription factors that are important mediators of lipid and cholesterol metabolism. LXR targets include the ABC family of transporters, ApoE, LPL, and SREBP (139, 140). Liver-specific LXR alpha knockout mice develop increased cholesterol levels and atherosclerosis (141). Sallam et al. performed genome-wide transcriptional profiling of primary mouse hepatocytes upon stimulation with an LXR agonist (78). The strongest induced gene was a non-coding RNA termed LeXis (liver-expressed LXR-induced sequence) that lies adjacent to the Abca1 gene. LeXis regulates several genes with roles in cholesterol biosynthesis, subsequently altering both liver and plasma cholesterol levels. Mass spectrometry was used to characterize the LeXis interactome and revealed binding to RALY, a ribonucleoprotein that acts a transcriptional cofactor in regulation of cholesterol biosynthetic genes (78). In the context of atherosclerosis, adenoviral overexpression of LeXis in the liver reduces atherosclerosis in a familial hypercholesterolemia mouse model (142). As discussed above, CHROME is another LXR-regulated lncRNA involved in cholesterol homeostasis (49). CHROME was first identified through a combination of genetic association studies for premature CAD and HDL-C and microarray based expression profiling in human atherosclerotic plaques (49). RNA-seq of control and CHROME shRNA treated HepG2 hepatocytes revealed downstream pathways affected, including the LXR pathway, bile acid metabolism, cholesterol excretion and fatty-acid β-oxidation pathways (49).

Application of Novel Genomic Technologies to Detect and Study LncRNA Functions

Novel Sequencing Technologies to Discover and Annotate Long Non-coding RNAs

Although next-generation sequencing has resulted in the identification of thousands of lncRNAs in the genome, many of these lncRNAs remain poorly characterized and annotated. It is often unclear where transcription begins and which exons are present in a particular isoform. Since lncRNAs are often expressed at lower levels compared to protein-coding genes, current transcriptomic data is unable to provide comprehensive mapping/characterization of isoforms. However, new sequencing technologies allow for better characterization due to longer read lengths, higher sensitivity, and higher accuracy. Techniques such as Iso-Seq (Pacific Biosystems) offer long-read sequencing using single-molecule, real-time (SMRT) sequencing, in which the sequence of a full-length transcript is captured in a single read (143). Despite these benefits, these single-molecule sequencers yield higher error rates compared with short read sequencing technologies (e.g., Illumina). Nanopore technologies such as the MinION instrument (Oxford Nanopore Technologies) also allow single cDNA molecules to be sequenced without the need for amplification, providing sufficient read lengths to cover the full-length non-coding RNA, and results in less bias than other long-read approaches (144). This technique passes nucleic acids through an orifice 10−9 m in diameter, where instrumental electric current changes are utilized to decipher the identity of each nucleotide (145).

Since lncRNAs are typically less abundant than protein coding genes (usually one order of magnitude less), they remain a challenge to study in bulk transcriptomic datasets. To improve the detection and annotation of lncRNAs, a method known as RACE (Rapid Amplification of cDNA Ends)-Seq was developed (146), however this approach was limited by its low-throughput. Later a technique called RNA CaptureSeq was developed to enrich for long non-coding RNAs (147). RNA CaptureSeq employs an array of oligonucleotide probes to capture select genes of interest, which can be applied to pull-down lncRNAs of interest (148, 149). More recently the GENCODE consortium improved upon RNA CaptureSeq by developing RNA Capture Long Seq (RNA CLS) with the goal of annotating lncRNAs with much higher confidence (150). RNA CLS overcomes the short-read length hurdle of RNA CaptureSeq by first capturing lncRNAs and then integrating with long-read sequencing.

DNA-Based LncRNA Interactions

Despite their low abundance, lncRNAs are known to function through specific molecular interactions with other RNA species and RNA binding proteins. Several high-throughput methods are now available to uncover the genomic DNA sequences that lncRNAs interact with and likely regulate (Figure 2A). Chromatin Isolation by RNA Purification (ChIRP-Seq) is a well-established technique to study lncRNA-chromatin interactions through RNA/chromatin crosslinking, purification using biotinylated antisense oligonucleotides, followed by high-throughput sequencing (151, 152). Domain-specific ChIRP (dChIRP) is a variation of ChIRP that can characterize lncRNA function and architecture at the RNA domain level (153). dChIRP can not only investigate lncRNA-chromatin interactions but also pairwise lncRNA-RNA and lncRNA-protein interactions.

Figure 2

Capture hybridization analysis of RNA targets (CHART) (154, 155) is a similar method to experimentally determine where lncRNAs target and localize in the genome (Figure 2A). In the CHART protocol chromatin is crosslinked and lncRNAs subsequently hybridized to biotinylated C-oligos. After bead immobilization of lncRNA/DNA complexes, sequencing is conducted to identify lncRNA binding DNA regions.

GRID-seq (global RNA interactions with DNA by deep sequencing) is a new unbiased method to capture global RNA-interactions (Figure 2A) that can be applied to investigate lncRNA-DNA interactions in cell lines relevant to atherosclerosis (156). This GRID-seq technique uses a bivalent linker consisting of double-stranded DNA and single-stranded RNA to link RNAs with DNA in nuclei that have been fixed. Finally, MARGI (mapping RNA-genome interactions) is a high-throughput method that can be performed in vivo or on cells and reveal the genomic target sites of lncRNAs (157).

Protein-Based LncRNA Interactions

ChIRP-MS is an adaptation of the ChIRP protocol and used to characterize the interacting proteome for a lncRNA (158). ChIRP-MS has identified protein interactors for lncRNAs such as LeXis, MeXis, and AK098656 (50, 76, 78). lncRNA pull-down followed by mass spectrometry has been conducted for several lncRNAs with potential roles in CAD such as circANRIL (23), STEEL (70), MALAT1 (159, 160), Dnm3os (73), lncLSTR (79), and GATA6-AS (62). Numerous additional methods exist to decipher the proteins binding lncRNAs. RAP-MS uses ultraviolet light to crosslink direct RNA-protein interactions (161). UV-C crosslinking immunoprecipitation (CLIP) is another powerful technique to interrogate direct protein-RNA interactions and many variations have been adapted based on the implementation of high-throughput sequencing (Figure 2B). These include iCLIP, PAR-CLIP, HITS-CLIP, irCLIP, and eCLIP (162166). High-throughput sequencing of RNA isolation by crosslinking immunoprecipitation (HITS-CLIP) was developed as genome-wide means to interrogate RNA-protein interactions in vivo (164).

TRIBE (targets of RNA-binding proteins identified by editing) is designed for identifying RNA molecules that bind to RNA binding proteins (RBP) (Figure 2B) (167). Advantages of TRIBE include application to in vivo samples, ability to performed on a small number of cells, and no need for antibodies in the procedure. The TRIBE protocol couples an RNA editing enzyme to the RBP and RNA targets that have been edited are identified via next-generation sequencing (TRIBE-seq). HyperTRIBE extends upon the TRIBE procedure by introducing a hyperactive mutation into the RNA editing enzyme, which improves the RNA editing efficiency and reduces the sequence bias of editing (168).

RNA-Based LncRNA Interactions

RNA-centric RNA antisense purification (RAP) is a general approach to identify and study lncRNA functions (Figure 2C). This method uses long capture probes (120 nucleotides) tiled across an entire RNA sequence to pull down lncRNAs, followed by stringent wash conditions to reduce non-specific binding (169). There are now next-generation sequencing derived methodologies that have been established to better define RNA-RNA interactions. LIGR-seq (LIGation of interacting RNA followed by high-throughput sequencing) can capture base-paired RNA-RNA interactions (Figure 2C) (170). In LIGR-seq, RNA duplexes are cross linked with the psoralen derivative 4'-aminomethyltrioxalen (AMT) along with UV irradiation at 365 nm, and RNase R is added to digest linear and structural RNAs. This step enriches for AMT-crosslinked RNA-RNA duplexes that are subsequently subjected to next-generation sequencing. Though LIGR-seq does not work well for small RNAs such as microRNA (miRNA), it should be able to uncover novel dynamic and long range interactions between lncRNAs and other RNA molecules. Various other methods have been developed to study the RNA interactome for lncRNAs with functional relevance in atherosclerosis including PARIS (Psoralen Analysis of RNA Interactions and Structures) (171), SPLASH (172), and MARIO (173). These techniques all can provide valuable information because many lncRNA sequence and structural motifs act as functional scaffolds in the assembly of RNA-protein complexes (17). However, it should be noted that many of these assays could be biased toward capturing stable interactions, while more transient and stimulation specific interactions may require some enrichment steps.

In situ Hybridization-Based Methods

A critical consideration when interrogating a given lncRNA function, is identifying its endogenous tissue and cellular localization. While many lncRNAs are expected to be cytosolic and contribute to post-transcriptional, translational or post-translational gene regulation, nuclear lncRNAs could participate in transcriptional regulation, chromatin structure or mRNA export mechanisms (17). RNA Fluorescence in situ Hybridization (FISH) has been a traditional method to identify the subcellular localization of RNA within cells, however it lacks sensitivity for lowly expressed lncRNAs. Single-molecule RNA FISH (smFISH) is a quantitative technique that provides the sensitivity to detect these lncRNAs and measures absolute transcript levels by using multiple short probes per target RNA (174). However, given that smFISH relies heavily on the optical detection of a limited number of fluorophores, it is restricted in its multiplexing capacity. Attempts to overcome this issue include implementation of combinatorial labeling by spectral barcodes and the incorporation of sequential hybridizations (seqFISH) using different colored probes in each hybridization round (175, 176). In seqFISH individual transcripts are imaged as different colored dots and quantified by counting the number of dots. Multiplex error-robust combinatorial labeling (merFISH) is an in situ targeted approach that utilizes two-step labeling and the detection of binary barcodes assigned to specific targets. This is accomplished by several rounds of hybridization, imaging, and cleavage of fluorophores from probes conjugated to readout sequences that interchange each cycle. Hybridization to readout sequences by the merFISH technique is much less time consuming than methods that utilize hybridization directly to target RNAs (177, 178). RNA SPOTS (sequential probing of targets) follows the same rationale as merFISH, except that it is used in vitro instead of in situ (179).

While still in the nascent stage, the emergence of spatial transcriptomics facilitates integration of RNA-seq expression data with spatial locations of RNA molecules in individual tissue sections (180). In this procedure, fixed tissue samples are annealed to regionally barcoded reverse transcription primers. Following reverse transcription, RNA-seq followed by computational reconstruction allows the two-dimensional localization and quantification of RNA molecules (180). This barcoded method has already been applied to spatially resolve gene expression in the human adult heart (181). While this procedure was originally developed to study mRNAs, it shows promise for the spatial resolution of lncRNAs, given the increased sensitivity and ability to identify context-specific expression profiles. One consideration for atherosclerosis FISH experiments, is that heterogeneous cell types in lesions may be impacted differently by various fixation and hybridization conditions, so careful titration of reagents is recommended.

Other LncRNA Functions

Finally, there are various omics methods that can define the dynamics of lncRNA transcription, stability and RNA modifications. Nascent RNA sequencing analysis, including global nuclear run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq) assays, could enable comprehensive detection of transient RNA transcriptional events for multiple RNA species, including mRNA, lncRNA, and eRNA (182). While most transcriptomic datasets capture steady-state levels of lncRNA transcripts, they do not provide direct insights into the stability of lncRNAs. BRIC-seq (5'-bromo-uridine immunoprecipitation chase-deep sequencing analysis) is a method that pulse-labels endogenous RNAs and employs next-generation sequencing to measure RNA decay over time (183). Total RNAs (including lncRNAs) can be isolated from cells at desired time points under various cell-specific perturbations to facilitate functional analysis of lncRNA stability (e.g., lncRNA related to CAD). For example, direct measurements of lncRNA stability in response to CRISPR based loss/gain of gene function or drug treatments could be examined.

Another technique, ICE-seq (inosine chemical erasing coupled with sequencing) (184) represents a promising approach to globally identify lncRNA adenosine to inosine modifications (e.g., in the context of atherosclerosis). Adenosine to inosine (A-to-I) RNA editing is the most abundant form of RNA editing in humans and results from adenosine deaminase acting on RNA (ADAR). A-to-I editing is common to all lncRNAs and affects lncRNA function through altered stability and target recognition (185, 186). A-to-I RNA editing of mRNA has already been demonstrated to have important functional consequences in atherosclerosis. For example, A-to-I editing of cathepsin S mRNA (CTSS) is associated with cathepsin S levels in patients with atherosclerosis. Treatment of endothelial cells with inflammatory cytokines or exposure to hypoxia was shown to induce cathepsin S RNA editing and gene/protein expression (187).

Novel Computational Tools for LncRNA Annotation and Functional Prediction

Genomic annotation of lncRNA sequences requires defining the precise genomic coordinates of lncRNA exons and their respective transcription start sites. LncRNA annotation also involves functional annotation with respect to predicted biological mechanisms, subcellular localization, and affected cell types/tissues. While lncRNAs share some similarities with mRNAs such as transcript length and splicing structure (188), proper identification and characterization of specific long non-coding transcripts still remains a challenge. Unlike mRNAs, lncRNAs often exhibit lower stability, lower abundance, less splicing and greater nuclear localization (189).

With the widespread application of high-throughput sequencing technologies, both automated and manual methods have been adopted to properly define lncRNA sequences from RNA-seq data. Automated annotation generates a larger catalog of lncRNAs and harnesses a transcriptome assembly consisting of two distinct strategies. In one automated approach, reads are first aligned to the reference genome to reveal all the possible splicing events which are subsequently assembled into transcripts (190, 191). In another automated approach, transcripts are built de novo from experimental reads and later aligned to a particular reference genome. Fu et al. (192) used both short and long sequencing reads to demonstrate superior sensitivity of transcript assembly and isoform annotation accuracy with the de novo approach. Automated assembly is fast as it does not require wet-lab based characterization, and it is considerable cheaper than the manual approach (144).

Although it produces a smaller catalog of lncRNAs compared to the automated method, manual annotation produces higher quality lncRNA transcript sequences and thus improves functional characterization. The widely adopted GENCODE project annotation of lncRNAs utilizes a manual curation approach (193), and integrates different sources of data together with computational analyses to generate a transcript model. cDNA and expressed sequence tag (EST) sequences deposited in publicly available databases are typically the starting point for manually annotating lncRNA transcripts. These are integrated with Cap Analysis of Gene Expression (5′-CAGE) and poly(A) position profiling by sequencing (3P-seq) to characterize 5′ and 3′ ends, respectively. These manually annotated transcripts are then mapped to reference genomes and assigned exon and splice site locations (119). The RefSeq (Reference Sequence) project also implements manual annotation of long non-coding RNAs that are integrated with automated methods (194). Manually annotated lncRNAs can be further divided into subclasses such as intergenic lncRNAs (lincRNAs), antisense lncRNAs, and intronic lncRNAs. As cDNA annotation depends on the availability of full length transcripts, manual annotation focuses primarily on genomic annotation. As a result the manual approach produces a more comprehensive set of pseudogenes and alternatively spliced transcripts (193).

Another comprehensive database established in 2016 is NONCODE that dedicates itself to collecting lncRNAs through integration with other databases (e.g., RefSeq and Ensembl) and exhaustive annotation. Compared to these other databases, NONCODE has collected more lncRNA transcripts (excluding tRNAs and rRNAs) and provides unique annotations of lncRNAs (e.g., RNA secondary structure, expression in exosomes, associations between lncRNA and disease) (6, 195). NONCODE also provides lncRNAs for over 15 species including mouse, zebrafish, and C-elegans. The latest version of NONCODE (v5), which also captures lncRNAs from the literature, consists of nearly 550,000 annotated lncRNAs (195).

There are now an array of computational tools to annotate the sequences and functions of the expanding catalog of lncRNAs, as described in Table 2. Existing computational methods for lncRNA identification include those that require a reference genome and those that are reference-free. Examples of methods requiring a reference genome include UClncR, lncScore (205), COME (206), and lncRScan-SVM (207). Reference-free methods to identify lncRNAs from RNA-seq data include LncADeep (198), lncRNAnet (197, 208), FEElnc (197), longdist (204), lncRNA-MFDL (209), and CPC2 (210). Many of these tools employ artificial intelligence algorithms (e.g., machine learning, deep learning) in order to distinguish lncRNAs from their protein-coding transcript counterparts.

Table 2

Tool name and reference*Model#Code availableApplicationPerformance results
AnnoLnc (196)Statistical approachhttp://annolnc.cbi.pku.edu.cnAnnotation of human lncRNAs.Not reported.
FEELnc (197)Random foresthttps://github.com/tderrien/FEELncAnnotation of lncRNAs.High classification power (AUC = 0.97).
LncADeep (198)Deep belief network built as a stack of restricted Boltzmann machineshttps://github.com/cyang235/LncADeepIdentification and functional annotation for lncRNAsWith 10-fold cross validation, average sensitivity of 98.1% and specificity of 97.2% and an average harmonic mean of 97.7%
LncFunTK (199)Statistical approachhttps://github.com/zhoujj2013/lncfuntkTo integrate ChIP-seq, CLIP-seq and RNA-seq data to predict, prioritize and annotate lncRNA functions.Calculates a Functional Information Score (FIS) to quantitatively predict functional importance.
lncLocator (200)Ensemble of support vector machine and random forest classifiers.http://www.csbio.sjtu.edu.cn/bioinf/lncLocator/To predict lncRNA subcellular localizations.Accuracy of 59% for prediction.
PennDiff (201)Regression-based statistical approachhttps://github.com/tigerhu15/PennDiffTo detect differential transcript isoforms from RNA-seq dataBased on both annotations (RefSeq and Ensembl), estimates from PennDiff have Spearman correlation coefficients of 0.87 and 0.76, respectively.
SEEKR (202)Statistical approachhttps://github.com/CalabreseLabPrediction of lncRNA subcellular localization, protein interactorsLncRNAs of related function have similar k-mer profiles, despite linear sequence similarity
UClncR (203)Statistical approachhttp://bioinformaticstools.mayo.edu/research/UClncRPerforms transcript assembly, prediction of lncRNA candidates in bulk RNA-seq data, quantification and annotation both known and novel lncRNA candidates.For lincRNA prediction, UClncR reported 66 “novel” lincRNA transcripts and 12 lncRNAs overlapping with nearby genes (the recall rate of 90.7%).
A support vector machine based method to distinguish long non-coding RNAs from protein transcripts (204)Support vector machinehttps://github.com/hugowschneider/longdist.pyTo distinguish lncRNAs from protein coding transcripts.98.21% accuracy in classifying long non-coding RNAs from protein coding transcripts.

Comparison of selected computational tools.

*

Three of the publications have not been constructed into available tools but rather represent a framework for analysis.

#

Model type does not include preprocessing which may or may not including alignment of protein-coding regions. The link is provided if the code is available otherwise the column is marked with an “X”. AUC, area under the curve.

Unlike protein functions that can be inferred from protein-coding sequences, it is more difficult to infer lncRNA function from RNA sequences. Zhou et al. developed a tool, lncFunTK that calculates a Functional Information Score (FIS) to quantitatively measure the functional importance of a lncRNA (199), based on the top Gene Ontology and inferred regulatory networks for lncRNAs and their neighboring genes. Another tool, FEELnc, annotates lncRNA function by evaluating neighboring genes to predict both lncRNA function and mRNA partners (197). Given that lncRNA function often depends on subcellular localization, the lncLocator tool predicts five lncRNA categories: nucleus, cytoplasm, cytosol, exosome, and ribosome (200). LncADeep provides enriched pathways and functional modules for lncRNA functional annotation by integrating KEGG and Reactome Pathway databases in a deep learning framework (211). A novel method for lncRNA classification is SEEKR, which counts lncRNA k-mer frequencies from nucleotide sequences, which may be correlated with lncRNA localization or protein binding (202).

Like other ncRNAs, lncRNA functional annotation also depends on accurate secondary structure prediction. There are several computational tools to predict lncRNA structures such as pysster (212), RNAfold (213), RNAstructure (214), and UNAfold (215). Other tools are available to study the evolutionary conservation of lncRNAs, including EvoFold (216), Evolinc (217), RNAz (218), and SISSIz (219). Since lncRNAs are generally poorly conserved, there are strategies to examine remnants of protein coding sequences in these lncRNAs (220). Identifying evolutionarily conserved lncRNA structures and binding domains may provide clues to predict lncRNA function for follow-up experimentation.

Conclusion

The emergence of RNA-seq and other omics technologies in the past decade have catalyzed the identification of a plethora of novel lncRNAs. To date, more than 30 lncRNAs with functional relevance to CAD have been characterized (Table 1), yet numerous lncRNAs remain to be studied in greater detail that are linked to endothelial, smooth muscle, macrophage, and lipid traits. With the growing number of CAD GWAS candidate loci harboring lncRNAs, and improved fine-mapping and annotation approaches, there is an opportunity to functionally dissect these regions to develop novel strategies to target non-coding genomic risk factors. As outlined in this review, a multi-faceted approach is likely required to successfully prioritize and study these lncRNAs, which may include implementation of long-read and high-depth sequencing, improved computational tools, coupled with orthogonal high-throughput experimental validation assays. Careful consideration of the lower abundance, context-specific expression of lncRNAs, and thoughtful study designs may improve chances of success in these multi-omics assays. However, it should also be noted that in many cases, more traditional and lower throughput approaches would be equally appropriate to characterize a given lncRNA, thus reducing the overall costs and required expertise.

For conserved lncRNAs with predicted roles in altering CAD pathogenesis, loss of function studies can be performed in animal models, such as the mouse (ApoE−/− or LDLR−/− backgrounds) or zebrafish. However, with the majority of human lncRNAs being poorly conserved across species, they may be better suited to studies in primary human cells or induced pluripotent stem cell (iPSC) derived vascular cells. In the context of CAD and other cardiometabolic disorders, genetic manipulation of lncRNAs via antisense oligonucleotides (221) or CRISPR/Cas9 to either delete (23, 222, 223) or activate/repress lncRNA expression, may lead to the identification of specific lncRNA binding partners, subcellular localization and functional insights relevant to CAD. lncRNA discovery/annotation can be further improved by integrating these genetic perturbations with high-dimensional transcriptomic and epigenomic assays (e.g., RNA-seq, ATAC-seq and ChIP-seq) to mark lncRNA promoters, decipher RNA polymerase and transcription factor binding, and reveal the dynamics of lncRNA regulatory activities. Single-cell based assays may also shed light on cell-specific markers and dynamics of lncRNAs across lineages (224, 225). Unraveling the complexity of lncRNA function in the setting of atherosclerosis may hold the key to delineate causal disease-associated pathways. In this regard it will also be important to determine whether lncRNAs operate synergistically, serve redundant and/or compensatory roles with other dysregulated lncRNAs and/or mRNAs associated with CAD.

Statements

Author contributions

AT and CM conceived of the manuscript. AT, DW, MK, CD, MP, and CM wrote the manuscript.

Funding

This work was supported by National Institutes of Health (NIH) grants R00 HL125912 (CM) and F31 NR017821 (CD) and a Leducq Foundation Transatlantic Network of Excellence award (CM). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    WeberCNoelsH. Atherosclerosis: current pathogenesis and therapeutic options. Nat Med. (2011) 17:141022. 10.1038/nm.2538

  • 2.

    TabasIGarcía-CardeñaGOwensGK. Recent insights into the cellular biology of atherosclerosis. J Cell Biol. (2015) 209:1322. 10.1083/jcb.201412052

  • 3.

    St. LaurentGVyatkinYKapranovP. Dark matter RNA illuminates the puzzle of genome-wide association studies. BMC Med. (2014) 12:97. 10.1186/1741-7015-12-97

  • 4.

    LiuYZhengLWangQHuY-W. Emerging roles and mechanisms of long noncoding RNAs in atherosclerosis. Int J Cardiol. (2017) 228:57082. 10.1016/j.ijcard.2016.11.182

  • 5.

    DjebaliSDavisCAMerkelADobinALassmannTMortazaviAet al. Landscape of transcription in human cells. Nature (2012) 489:1018. 10.1038/nature11233

  • 6.

    ZhaoYLiHFangSKangYWuWHaoYet al. NONCODE 2016: an informative and valuable data source of long non-coding RNAs. Nucleic Acids Res. (2016) 44:D2038. 10.1093/nar/gkv1252

  • 7.

    PalazzoAFLeeES. Non-coding RNA: what is functional and what is junk?Front Genet. (2015) 6:2. 10.3389/fgene.2015.00002

  • 8.

    BhartiyaDScariaV. Genomic variations in non-coding RNAs: Structure, function and regulation. Genomics (2016) 107:5968. 10.1016/j.ygeno.2016.01.005

  • 9.

    CostaMCLeitãoALEnguitaFJ. Biogenesis and mechanism of action of small non-coding RNAs: insights from the point of view of structural biology. Int J Mol Sci. (2012) 13:1026895. 10.3390/ijms130810268

  • 10.

    BhatSAAhmadSMMumtazPTMalikAADarMAUrwatUet al. Long non-coding RNAs: mechanism of action and functional utility. Noncoding RNA Res. (2016) 1:4350. 10.1016/j.ncrna.2016.11.002

  • 11.

    BartschDZirkelAKurianL. Characterization of Circular RNAs (circRNA) Associated with the Translation Machinery. Methods Mol Biol. (2018). 1724:15966. 10.1007/978-1-4939-7562-4_13

  • 12.

    ShiXSunMLiuHYaoYSongY. Long non-coding RNAs: a new frontier in the study of human diseases. Cancer Lett. (2013) 339:15966. 10.1016/j.canlet.2013.06.013

  • 13.

    GuttmanMAmitIGarberMFrenchCLinMFFeldserDet al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature (2009) 458:2237. 10.1038/nature07672

  • 14.

    PangKCFrithMCMattickJS. Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function. Trends Genet. (2006) 22:15. 10.1016/j.tig.2005.10.003

  • 15.

    QuZAdelsonDL. Evolutionary conservation and functional roles of ncRNA. Front Genet. (2012) 3:205. 10.3389/fgene.2012.00205

  • 16.

    DiederichsS. The four dimensions of noncoding RNA conservation. Trends Genet. (2014) 30:1213. 10.1016/j.tig.2014.01.004

  • 17.

    WangKCChangHY. Molecular mechanisms of long noncoding RNAs. Mol Cell. (2011) 43:90414. 10.1016/j.molcel.2011.08.018

  • 18.

    NatarelliLGeißlerCCsabaGWeiYZhuMdi FrancescoAet al. miR-103 promotes endothelial maladaptation by targeting lncWDR59. Nat Commun. (2018) 9:2645. 10.1038/s41467-018-05065-z

  • 19.

    SantovitoDEgeaVWeberC. Small but smart: MicroRNAs orchestrate atherosclerosis development and progression. Biochim Biophys Acta (2016) 1861(12 Pt B):207586. 10.1016/j.bbalip.2015.12.013

  • 20.

    FeinbergMWMooreKJ. Microrna regulation of atherosclerosis. Circ Res. (2016) 118:70320. 10.1161/CIRCRESAHA.115.306300

  • 21.

    JarinovaOStewartAFRRobertsRWellsGLauPNaingTet al. Functional analysis of the chromosome 9p21.3 coronary artery disease risk locus. Arterioscler Thromb Vasc Biol. (2009) 29:16717. 10.1161/ATVBAHA.109.189522

  • 22.

    HoldtLMHoffmannSSassKLangenbergerDScholzMKrohnKet al. Alu elements in ANRIL non-coding RNA at chromosome 9p21 modulate atherogenic cell functions through trans-regulation of gene networks. PLoS Genet. (2013) 9:e1003588. 10.1371/journal.pgen.1003588

  • 23.

    HoldtLMStahringerASassKPichlerGKulakNAWilfertWet al. Circular non-coding RNA ANRIL modulates ribosomal RNA maturation and atherosclerosis in humans. Nat Commun. (2016) 7:12429. 10.1038/ncomms12429

  • 24.

    PasmantELaurendeauIHéronDVidaudMVidaudDBiècheI. Characterization of a germ-line deletion, including the entire INK4/ARF locus, in a melanoma-neural system tumor family: identification of ANRIL, an antisense noncoding RNA whose expression coclusters with ARF. Cancer Res. (2007) 67:39639. 10.1158/0008-5472.CAN-06-2004

  • 25.

    SchneiderCKingRMPhilipsonL. Genes specifically expressed at growth arrest of mammalian cells. Cell (1988) 54:78793. 10.1016/S0092-8674(88)91065-3

  • 26.

    CocciaEMCicalaCCharlesworthACiccarelliCRossiGBPhilipsonLet al. Regulation and expression of a growth arrest-specific gene (gas5) during growth, differentiation, and development. Mol Cell Biol. (1992) 12:351421. 10.1128/MCB.12.8.3514

  • 27.

    ChenLYangWGuoYChenWZhengPZengJet al. Exosomal lncRNA GAS5 regulates the apoptosis of macrophages and vascular endothelial cells in atherosclerosis. PLoS ONE (2017) 12:e0185406. 10.1371/journal.pone.0185406

  • 28.

    TangRZhangGWangY-CMeiXChenS-Y. The long non-coding RNA GAS5 regulates transforming growth factor β (TGF-β)-induced smooth muscle cell differentiation via RNA Smad-binding elements. J Biol Chem. (2017) 292:142708. 10.1074/jbc.M117.790030

  • 29.

    WangSZhangXYuanYTanMZhangLXueXet al. BRG1 expression is increased in thoracic aortic aneurysms and regulates proliferation and apoptosis of vascular smooth muscle cells through the long non-coding RNA HIF1A-AS1 in vitro. Eur J Cardiothorac Surg. (2015) 47:43946. 10.1093/ejcts/ezu215

  • 30.

    BaranelloLBertozziDFogliMVPommierYCapranicoG. DNA topoisomerase I inhibition by camptothecin induces escape of RNA polymerase II from promoter-proximal pause site, antisense transcription and histone acetylation at the human HIF-1alpha gene locus. Nucleic Acids Res. (2010) 38:15971. 10.1093/nar/gkp817

  • 31.

    BertozziDIurlaroRSordetOMarinelloJZaffaroniNCapranicoG. Characterization of novel antisense HIF-1α transcripts in human cancers. Cell Cycle (2011) 10:318997. 10.4161/cc.10.18.17183

  • 32.

    ZhaoYFengGWangYYueYZhaoW. Regulation of apoptosis by long non-coding RNA HIF1A-AS1 in VSMCs: implications for TAA pathogenesis. Int J Clin Exp Pathol. (2014) 7:764352.

  • 33.

    WuGCaiJHanYChenJHuangZ-PChenCet al. LincRNA-p21 regulates neointima formation, vascular smooth muscle cell proliferation, apoptosis, and atherosclerosis by enhancing p53 activity. Circulation (2014) 130:145265. 10.1161/CIRCULATIONAHA.114.011675

  • 34.

    HuarteMGuttmanMFeldserDGarberMKoziolMJKenzelmann-BrozDet al. A large intergenic noncoding RNA induced by p53 mediates global gene repression in the p53 response. Cell (2010) 142:40919. 10.1016/j.cell.2010.06.040

  • 35.

    MichalikKMYouXManavskiYDoddaballapurAZörnigMBraunTet al. Long noncoding RNA MALAT1 regulates endothelial cell function and vessel growth. Circ Res. (2014) 114:138997. 10.1161/CIRCRESAHA.114.303265

  • 36.

    SongT-FHuangL-WYuanYWangH-QHeH-PMaW-Jet al. LncRNA MALAT1 regulates smooth muscle cell phenotype switch via activation of autophagy. Oncotarget. (2018) 9:441126. 10.18632/oncotarget.23230

  • 37.

    LuoJ-HRenBKeryanovSTsengGCRaoUNMMongaSPet al. Transcriptomic and genomic analysis of human hepatocellular carcinomas and hepatoblastomas. Hepatology (2006) 44:101224. 10.1002/hep.21328

  • 38.

    GuoFLiYLiuYWangJLiYLiG. Inhibition of metastasis-associated lung adenocarcinoma transcript 1 in CaSki human cervical cancer cells suppresses cell proliferation and invasion. Acta Biochim Biophys Sin (Shanghai) (2010) 42:2249. 10.1093/abbs/gmq008

  • 39.

    KoshimizuTFujiwaraYSakaiNShibataKTsuchiyaH. Oxytocin stimulates expression of a noncoding RNA tumor marker in a human neuroblastoma cell line. Life Sci. (2010) 86(11–12):45560. 10.1016/j.lfs.2010.02.001

  • 40.

    GuffantiAIaconoMPelucchiPKimNSoldàGCroftLJet al. A transcriptional sketch of a primary human breast cancer by 454 deep sequencing. BMC Genomics (2009) 10:163. 10.1186/1471-2164-10-163

  • 41.

    MiyoshiNWagatsumaHWakanaSShiroishiTNomuraMAisakaKet al. Identification of an imprinted gene, Meg3/Gtl2 and its human homologue MEG3, first mapped on mouse distal chromosome 12 and human chromosome 14q. Genes Cells (2000) 5:21120. 10.1046/j.1365-2443.2000.00320.x

  • 42.

    WuZHeYLiDFangXShangTZhangHet al. Long noncoding RNA MEG3 suppressed endothelial cell proliferation and migration through regulating miR-21. Am J Transl Res. (2017) 9:332635.

  • 43.

    ZhuBGongYYanGWangDQiaoYWangQet al. Down-regulation of lncRNA MEG3 promotes hypoxia-induced human pulmonary artery smooth muscle cell proliferation and migration via repressing PTEN by sponging miR-21. Biochem Biophys Res Commun. (2018) 495:212532. 10.1016/j.bbrc.2017.11.185

  • 44.

    WangLXiaJ-WKeZ-PZhangB-H. Blockade of NEAT1 represses inflammation response and lipid uptake via modulating miR-342–3 p in human macrophages THP-1 cells. J Cell Physiol. (2019) 234:531926. 10.1002/jcp.27340

  • 45.

    AhmedASIDongKLiuJWenTYuLXuFet al. Long noncoding RNA NEAT1 (nuclear paraspeckle assembly transcript 1) is critical for phenotypic switching of vascular smooth muscle cells. Proc Natl Acad Sci USA. (2018) 115:E86607. 10.1073/pnas.1803725115

  • 46.

    HutchinsonJNEnsmingerAWClemsonCMLynchCRLawrenceJBChessA. A screen for nuclear transcripts identifies two linked noncoding RNAs associated with SC35 splicing domains. BMC Genomics (2007) 8:39. 10.1186/1471-2164-8-39

  • 47.

    ImamuraKImamachiNAkizukiGKumakuraMKawaguchiANagataKet al. Long noncoding RNA NEAT1-dependent SFPQ relocation from promoter region to paraspeckle mediates IL8 expression upon immune stimuli. Mol Cell (2014) 53:393406. 10.1016/j.molcel.2014.06.013

  • 48.

    HuY-WGuoF-XXuY-JLiPLuZ-FMcVeyDGet al. Long non-coding RNA NEXN-AS1 mitigates atherosclerosis by regulating the actin-binding protein NEXN. J Clin Invest. (2018). 10.1172/JCI98230 [Epub ahead of print].

  • 49.

    HennessyEJvan SolingenCScacalossiKROuimetMAfonsoMSPrinsJet al. The long noncoding RNA CHROME regulates cholesterol homeostasis in primates. Nat Metab. (2018). 1:98110. 10.1038/s42255-018-0004-9

  • 50.

    JinLLinXYangLFanXWangWLiSet al. AK098656, a novel vascular smooth muscle cell-dominant long noncoding RNA, promotes hypertension. Hypertension (2018) 71:26272. 10.1161/HYPERTENSIONAHA.117.09651

  • 51.

    YuCKXuTAssoianRKRaderDJ. Mining the stiffness-sensitive transcriptome in human vascular smooth muscle cells identifies long noncoding RNA stiffness regulators. Arterioscler Thromb Vasc Biol. (2018) 38:16473. 10.1161/ATVBAHA.117.310237

  • 52.

    LessardLLiuMMarzeseDMWangHChongKKawasNet al. The CASC15 long intergenic noncoding RNA locus is involved in melanoma progression and phenotype switching. J Invest Dermatol. (2015) 135:246474. 10.1038/jid.2015.200

  • 53.

    FeilRWalterJAllenNDReikW. Developmental control of allelic methylation in the imprinted mouse Igf2 and H19 genes. Development (1994) 120:293343.

  • 54.

    FeinbergAPKalikinLMJohnsonLAThompsonJS. Loss of imprinting in human cancer. Cold Spring Harb Symp Quant Biol. (1994) 59:35764. 10.1101/SQB.1994.059.01.040

  • 55.

    KallenANZhouX-BXuJQiaoCMaJYanLet al. The imprinted H19 lncRNA antagonizes let-7 microRNAs. Mol Cell (2013) 52:10112. 10.1016/j.molcel.2013.08.027

  • 56.

    LeungATracCJinWLantingLAkbanyASætromPet al. Novel long noncoding RNAs are regulated by angiotensin II in vascular smooth muscle cells. Circ Res. (2013) 113:26678. 10.1161/CIRCRESAHA.112.300849

  • 57.

    ZhaoJZhangWLinMWuWJiangPTouEet al. MYOSLID is a novel serum response factor-dependent long noncoding rna that amplifies the vascular smooth muscle differentiation program. Arterioscler Thromb Vasc Biol. (2016) 36:208899. 10.1161/ATVBAHA.116.307879

  • 58.

    KrawczykMEmersonBM. p50-associated COX-2 extragenic RNA (PACER) activates COX-2 gene expression by occluding repressive NF-κB complexes. Elife (2014) 3:e01776. 10.7554/eLife.01776

  • 59.

    BellRDLongXLinMBergmannJHNandaVCowanSLet al. Identification and initial functional characterization of a human vascular cell-enriched long noncoding RNA. Arterioscler Thromb Vasc Biol. (2014) 34:124959. 10.1161/ATVBAHA.114.303240

  • 60.

    BallantyneMDPinelKDakinRVeseyATDiverLMackenzieRet al. Smooth muscle enriched long noncoding RNA (SMILR) regulates cell proliferation. Circulation (2016) 133:205065. 10.1161/CIRCULATIONAHA.115.021019

  • 61.

    SinghKKMatkarPNPanYQuanAGuptaVTeohHet al. Endothelial long non-coding RNAs regulated by oxidized LDL. Mol Cell Biochem. (2017) 431:13949. 10.1007/s11010-017-2984-2

  • 62.

    NeumannPJaéNKnauAGlaserSFFouaniYRossbachOet al. The lncRNA GATA6-AS epigenetically regulates endothelial gene expression via interaction with LOXL2. Nat Commun. (2018) 9:237. 10.1038/s41467-017-02431-1

  • 63.

    MiaoYAjamiNEHuangT-SLinF-MLouC-HWangY-Tet al. Enhancer-associated long non-coding RNA LEENE regulates endothelial nitric oxide synthase and endothelial function. Nat Commun. (2018) 9:292. 10.1038/s41467-017-02113-y

  • 64.

    IshiiNOzakiKSatoHMizunoHSaitoSTakahashiAet al. Identification of a novel non-coding RNA, MIAT, that confers risk of myocardial infarction. J Hum Genet. (2006) 51:108799. 10.1007/s10038-006-0070-9

  • 65.

    YanBYaoJLiuJ-YLiX-MWangX-QLiY-Jet al. lncRNA-MIAT regulates microvascular dysfunction by functioning as a competing endogenous RNA. Circ Res. (2015) 116:114356. 10.1161/CIRCRESAHA.116.305510

  • 66.

    WillinghamATOrthAPBatalovSPetersECWenBGAza-BlancPet al. A strategy for probing the function of noncoding RNAs finds a repressor of NFAT. Science (2005) 309:15703. 10.1126/science.1115901

  • 67.

    TangYWoLChaiH. [Effects of noncoding RNA NRON gene regulation on human umbilical vein endothelial cells functions]. Zhonghua Xin Xue Guan Bing Za Zhi. (2013) 41:24550. 10.3760/cma.j.issn.0253-3758.2013.03.015

  • 68.

    TsaoH-WTaiT-STsengWChangH-HGrenninglohRMiawS-Cet al. Ets-1 facilitates nuclear entry of NFAT proteins and their recruitment to the IL-2 promoter. Proc Natl Acad Sci USA. (2013) 110:1577681. 10.1073/pnas.1304343110

  • 69.

    RobbGBCarsonARTaiSCFishJESinghSYamadaTet al. Post-transcriptional regulation of endothelial nitric-oxide synthase by an overlapping antisense mRNA transcript. J Biol Chem. (2004) 279:3798296. 10.1074/jbc.M400271200

  • 70.

    ManHSJSukumarANLamGCTurgeonPJYanMSKuKHet al. Angiogenic patterning by STEEL, an endothelial-enriched long noncoding RNA. Proc Natl Acad Sci USA. (2018) 115:24016. 10.1073/pnas.1715182115

  • 71.

    LiKBlumYVermaALiuZPramanikKLeighNRet al. A noncoding antisense RNA in tie-1 locus regulates tie-1 function in vivo. Blood (2010) 115:1339. 10.1182/blood-2009-09-242180

  • 72.

    ChowdhuryTAKocejaCEisa-BeygiSKleinstiverBPKumarSNLinC-Wet al. Temporal and Spatial Post-Transcriptional Regulation of Zebrafish tie1 mRNA by long noncoding RNA during brain vascular assembly. Arterioscler Thromb Vasc Biol. (2018) 38:156275. 10.1161/ATVBAHA.118.310848

  • 73.

    DasSReddyMASenapatiPStapletonKLantingLWangMet al. Diabetes mellitus-induced long noncoding RNA Dnm3os regulates macrophage functions and inflammation via nuclear mechanisms. Arterioscler Thromb Vasc Biol. (2018) 38:180620. 10.1161/ATVBAHA.117.310663

  • 74.

    LoebelDAFTsoiBWongNTamPPL. A conserved noncoding intronic transcript at the mouse Dnm3 locus. Genomics (2005) 85:7829. 10.1016/j.ygeno.2005.02.001

  • 75.

    DuMYuanLTanXHuangDWangXZhengZet al. The LPS-inducible lncRNA Mirt2 is a negative regulator of inflammation. Nat Commun. (2017) 8:2049. 10.1038/s41467-017-02229-1

  • 76.

    SallamTJonesMThomasBJWuXGillilandTQianKet al. Transcriptional regulation of macrophage cholesterol efflux and atherogenesis by a long noncoding RNA. Nat Med. (2018) 24:30412. 10.1038/nm.4479

  • 77.

    HalleyPKadakkuzhaBMFaghihiMAMagistriMZeierZKhorkovaOet al. Regulation of the apolipoprotein gene cluster by a long noncoding RNA. Cell Rep. (2014) 6:22230. 10.1016/j.celrep.2013.12.015

  • 78.

    SallamTJonesMCGillilandTZhangLWuXEskinAet al. Feedback modulation of cholesterol metabolism by the lipid-responsive non-coding RNA LeXis. Nature (2016) 534:1248. 10.1038/nature17674

  • 79.

    LiPRuanXYangLKiesewetterKZhaoYLuoHet al. A liver-enriched long non-coding RNA, lncLSTR, regulates systemic lipid metabolism in mice. Cell Metab. (2015) 21:45567. 10.1016/j.cmet.2015.02.004

  • 80.

    DouvrisASoubeyrandSNaingTMartinukANikpayMWilliamsAet al. Functional analysis of the TRIB1 associated locus linked to plasma triglycerides and coronary artery disease. J Am Heart Assoc. (2014) 3:e000884. 10.1161/JAHA.114.000884

  • 81.

    CongrainsAKamideKOhishiMRakugiH. ANRIL: molecular mechanisms and implications in human health. Int J Mol Sci. (2013) 14:127892. 10.3390/ijms14011278

  • 82.

    KotakeYNakagawaTKitagawaKSuzukiSLiuNKitagawaMet al. Long non-coding RNA ANRIL is required for the PRC2 recruitment to and silencing of p15(INK4B) tumor suppressor gene. Oncogene (2011) 30:195662. 10.1038/onc.2010.568

  • 83.

    YapKLLiSMuñoz-CabelloAMRaguzSZengLMujtabaSet al. Molecular interplay of the noncoding RNA ANRIL and methylated histone H3 lysine 27 by polycomb CBX7 in transcriptional silencing of INK4a. Mol Cell. (2010) 38:66274. 10.1016/j.molcel.2010.03.021

  • 84.

    MotterleAPuXWoodHXiaoQGorSNgFLet al. Functional analyses of coronary artery disease associated variation on chromosome 9p21 in vascular smooth muscle cells. Hum Mol Genet. (2012) 21:40219. 10.1093/hmg/dds224

  • 85.

    ZhouXHanXWittfeldtASunJLiuCWangXet al. Long non-coding RNA ANRIL regulates inflammatory responses as a novel component of NF-κB pathway. RNA Biol. (2016) 13:98108. 10.1080/15476286.2015.1122164

  • 86.

    ViselAZhuYMayDAfzalVGongEAttanasioCet al. Targeted deletion of the 9p21 non-coding coronary artery disease risk interval in mice. Nature (2010) 464:40912. 10.1038/nature08801

  • 87.

    HoldtLMTeupserD. Long Noncoding RNA ANRIL: Lnc-ing genetic variation at the chromosome 9p21 locus to molecular mechanisms of atherosclerosis. Front Cardiovasc Med. (2018) 5:145. 10.3389/fcvm.2018.00145

  • 88.

    ArslanSBerkanÖLalemTÖzbilümNGökselSKorkmazÖet al. Long non-coding RNAs in the atherosclerotic plaque. Atherosclerosis (2017) 266:17681. 10.1016/j.atherosclerosis.2017.10.012

  • 89.

    CremerSMichalikKMFischerAPfistererLJaéNWinterCet al. Hematopoietic deficiency of the long non-coding RNA MALAT1 promotes atherosclerosis and plaque inflammation. Circulation (2018). [Epub ahead of print]. 10.1161/CIRCULATIONAHA.117.029015

  • 90.

    KoldemirOÖzgürEGezerU. Accumulation of GAS5 in exosomes is a marker of apoptosis induction. Biomed Rep. (2017) 6:35862. 10.3892/br.2017.848

  • 91.

    AryalBRotllanNFernández-HernandoC. Noncoding RNAs and atherosclerosis. Curr Atheroscler Rep. (2014) 16:407. 10.1007/s11883-014-0407-3

  • 92.

    ZhouTDingJWangXZhengX. Long noncoding RNAs and atherosclerosis. Atherosclerosis (2016) 248:5161. 10.1016/j.atherosclerosis.2016.02.025

  • 93.

    GastMRauchBHNakagawaSHaghikiaAJasinaAHaasJet al. Immune system-mediated atherosclerosis caused by deficiency of long noncoding RNA MALAT1 in ApoE-/- mice. Cardiovasc Res. (2018) 115:30214. 10.1093/cvr/cvy202

  • 94.

    LiuJHuangG-QKeZ-P. Silence of long intergenic noncoding RNA HOTAIR ameliorates oxidative stress and inflammation response in ox-LDL-treated human macrophages by upregulating miR-330–5 p. J Cell Physiol. (2019) 234:513442. 10.1002/jcp.27317

  • 95.

    SoubeyrandSNaingTMartinukAMcPhersonR. ERK1/2 regulates hepatocyte Trib1 in response to mitochondrial dysfunction. Biochim Biophys Acta. (2013) 1833:340514. 10.1016/j.bbamcr.2013.10.001

  • 96.

    Kiss-TothEBagstaffSMSungHYJozsaVDempseyCCauntJCet al. Human tribbles, a protein family controlling mitogen-activated protein kinase cascades. J Biol Chem. (2004) 279:427038. 10.1074/jbc.M407732200

  • 97.

    HannouSAWoutersKPaumelleRStaelsB. Functional genomics of the CDKN2A/B locus in cardiovascular and metabolic disease: what have we learned from GWASs?Trends Endocrinol Metab. (2015) 26:17684. 10.1016/j.tem.2015.01.008

  • 98.

    HoldtLMBeutnerFScholzMGielenSGäbelGBergertHet al. ANRIL expression is associated with atherosclerosis risk at chromosome 9p21. Arterioscler Thromb Vasc Biol. (2010) 30:6207. 10.1161/ATVBAHA.109.196832

  • 99.

    ZhouJXuJZhangLLiuSMaYWenXet al. Combined single-cell profiling of lncRNAs and functional screening reveals that H19 is pivotal for embryonic hematopoietic stem cell development. Cell Stem Cell. (2019) 24:28598.e5. 10.1016/j.stem.2018.11.023

  • 100.

    YoshimuraHMatsudaYYamamotoMKamiyaSIshiwataT. Expression and role of long non-coding RNA H19 in carcinogenesis. Front Biosci. (2018) 23:61425. 10.2741/4608

  • 101.

    GaoWZhuMWangHZhaoSZhaoDYangYet al. Association of polymorphisms in long non-coding RNA H19 with coronary artery disease risk in a Chinese population. Mutat Res. (2015) 772:1522. 10.1016/j.mrfmmm.2014.12.009

  • 102.

    ZhuRLiuXHeZ. Long non-coding RNA H19 and MALAT1 gene variants in patients with ischemic stroke in a northern Chinese Han population. Mol Brain (2018) 11:58. 10.1186/s13041-018-0402-7

  • 103.

    KimDKZhangLDzauVJPrattRE. H19, a developmentally regulated gene, is reexpressed in rat vascular smooth muscle cells after injury. J Clin Invest. (1994) 93:35560. 10.1172/JCI116967

  • 104.

    SunWLvJDuanLLinRLiYLiSet al. Long noncoding RNA H19 promotes vascular remodeling by sponging let-7a to upregulate the expression of cyclin D1. Biochem Biophys Res Commun. (2019) 508:103842. 10.1016/j.bbrc.2018.11.185

  • 105.

    HofmannPSommerJTheodorouKKirchhofLFischerALiYet al. Long non-coding RNA H19 regulates endothelial cell aging via inhibition of STAT3 signalling. Cardiovasc Res. (2019) 115:23042. 10.1093/cvr/cvy206

  • 106.

    BoonRAHofmannPMichalikKMLozano-VidalNBerghäuserDFischerAet al. Long noncoding RNA meg3 controls endothelial cell aging and function: implications for regenerative angiogenesis. J Am Coll Cardiol. (2016) 68:258991. 10.1016/j.jacc.2016.09.949

  • 107.

    van der HarstPVerweijN. Identification of 64 novel genetic loci provides an expanded view on the genetic architecture of coronary artery disease. Circ Res. (2018) 122:43343. 10.1161/CIRCRESAHA.117.312086

  • 108.

    LiuBPjanicMWangTNguyenTGloudemansMRaoAet al. Genetic regulatory mechanisms of smooth muscle cells map to coronary artery disease risk loci. Am J Hum Genet. (2018) 103:37788. 10.1016/j.ajhg.2018.08.001

  • 109.

    ArabKParkYJLindrothAMSchäferAOakesCWeichenhanDet al. Long noncoding RNA TARID directs demethylation and activation of the tumor suppressor TCF21 via GADD45A. Mol Cell. (2014) 55:60414. 10.1016/j.molcel.2014.06.031

  • 110.

    AcharyaABaekSTHuangGEskiocakBGoetschSSungCYet al. The bHLH transcription factor Tcf21 is required for lineage-specific EMT of cardiac fibroblast progenitors. Development (2012) 139:213949. 10.1242/dev.079970

  • 111.

    MillerCLAndersonDRKunduRKRaiesdanaANürnbergSTDiazRet al. Disease-related growth factor and embryonic signaling pathways modulate an enhancer of TCF21 expression at the 6q23.2 coronary heart disease locus. PLoS Genet. (2013) 9:e1003652. 10.1371/journal.pgen.1003652

  • 112.

    MillerCLHaasUDiazRLeeperNJKunduRKPatlollaBet al. Coronary heart disease-associated variation in TCF21 disrupts a miR-224 binding site and miRNA-mediated regulation. PLoS Genet. (2014) 10:e1004263. 10.1371/journal.pgen.1004263

  • 113.

    NurnbergSTChengKRaiesdanaAKunduRMillerCLKimJBet al. Coronary artery disease associated transcription factor TCF21 regulates smooth muscle precursor cells that contribute to the fibrous cap. PLoS Genet. (2015) 11:e1005155. 10.1371/journal.pgen.1005155

  • 114.

    HormozdiariFvan de BuntMSegrèAVLiXJooJWJBilowMet al. Colocalization of GWAS and eQTL signals detects target genes. Am J Hum Genet. (2016) 99:124560. 10.1016/j.ajhg.2016.10.003

  • 115.

    GusevAKoAShiHBhatiaGChungWPenninxBWJHet al. Integrative approaches for large-scale transcriptome-wide association studies. Nat Genet. (2016) 48:24552. 10.1038/ng.3506

  • 116.

    KashiKHendersonLBonettiACarninciP. Discovery and functional analysis of lncRNAs: Methodologies to investigate an uncharacterized transcriptome. Biochim Biophys Acta. (2016) 1859:315. 10.1016/j.bbagrm.2015.10.010

  • 117.

    KukurbaKRMontgomerySB. RNA sequencing and analysis. Cold Spring Harb Protoc. (2015) 2015:95169. 10.1101/pdb.top084970

  • 118.

    FreedmanJEMianoJM, National Heart, Lung, and Blood Institute Workshop Participants*. Challenges and opportunities in linking long noncoding rnas to cardiovascular, lung, and blood diseases. Arterioscler Thromb Vasc Biol. (2017) 37:215. 10.1161/ATVBAHA.116.308513

  • 119.

    DerrienTJohnsonRBussottiGTanzerADjebaliSTilgnerHet al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. (2012) 22:177589. 10.1101/gr.132159.111

  • 120.

    ConesaAMadrigalPTarazonaSGomez-CabreroDCerveraAMcPhersonAet al. A survey of best practices for RNA-seq data analysis. Genome Biol. (2016) 17:13. 10.1186/s13059-016-0881-8

  • 121.

    CabiliMNTrapnellCGoffLKoziolMTazon-VegaBRegevAet al. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. (2011) 25:191527. 10.1101/gad.17446611

  • 122.

    GuilSSolerMPortelaACarrèreJFonallerasEGómezAet al. Intronic RNAs mediate EZH2 regulation of epigenetic targets. Nat Struct Mol Biol. (2012) 19:66470. 10.1038/nsmb.2315

  • 123.

    TsoiLCIyerMKStuartPESwindellWRGudjonssonJETejasviTet al. Analysis of long non-coding RNAs highlights tissue-specific expression patterns and epigenetic profiles in normal and psoriatic skin. Genome Biol. (2015) 16:24. 10.1186/s13059-014-0570-4

  • 124.

    Ackers-JohnsonMTalasilaASageAPLongXBotIMorrellNWet al. Myocardin regulates vascular smooth muscle cell inflammatory activation and disease. Arterioscler Thromb Vasc Biol. (2015) 35:81728. 10.1161/ATVBAHA.114.305218

  • 125.

    ChenJKitchenCMStrebJWMianoJM. Myocardin: a component of a molecular switch for smooth muscle differentiation. J Mol Cell Cardiol. (2002) 34:134556. 10.1006/jmcc.2002.2086

  • 126.

    WangD-ZLiSHockemeyerDSutherlandLWangZSchrattGet al. Potentiation of serum response factor activity by a family of myocardin-related transcription factors. Proc Natl Acad Sci USA. (2002) 99:1485560. 10.1073/pnas.222561499

  • 127.

    DuKLIpHSLiJChenMDandreFYuWet al. Myocardin is a critical serum response factor cofactor in the transcriptional program regulating smooth muscle cell differentiation. Mol Cell Biol. (2003) 23:242537. 10.1128/MCB.23.7.2425-2437.2003

  • 128.

    YoshidaTSinhaSDandréFWamhoffBRHoofnagleMHKremerBEet al. Myocardin is a key regulator of CArG-dependent transcription of multiple smooth muscle marker genes. Circ Res. (2003) 92:85664. 10.1161/01.RES.0000068405.49081.09

  • 129.

    RussellMRPenikisAOldridgeDAAlvarez-DominguezJRMcDanielLDiamondMet al. CASC15-S is a tumor suppressor lncRNA at the 6p22 neuroblastoma susceptibility locus. Cancer Res. (2015) 75:315566. 10.1158/0008-5472.CAN-14-3613

  • 130.

    LimY-HKwonD-HKimJParkWJKookHKimY-K. Identification of long noncoding RNAs involved in muscle differentiation. PLoS ONE (2018) 13:e0193898. 10.1371/journal.pone.0193898

  • 131.

    MantellaL-ESinghKKSandhuPKantoresCRamadanAKhyzhaNet al. Fingerprint of long non-coding RNA regulated by cyclic mechanical stretch in human aortic smooth muscle cells: implications for hypertension. Mol Cell Biochem. (2017) 435(1–2):16373. 10.1007/s11010-017-3065-2

  • 132.

    ZhangXTangXHamblinMHYinK-J. Long Non-Coding RNA Malat1 regulates angiogenesis in Hindlimb Ischemia. Int J Mol Sci. (2018) 19:E1723. 10.3390/ijms19061723

  • 133.

    HuangCHanJWuYLiSWangQLinWet al. Exosomal MALAT1 derived from oxidized low-density lipoprotein-treated endothelial cells promotes M2 macrophage polarization. Mol Med Report. (2018) 18:50915. 10.3892/mmr.2018.8982

  • 134.

    HeCYangWYangJDingJLiSWuHet al. Long noncoding RNA MEG3 negatively regulates proliferation and angiogenesis in vascular endothelial cells. DNA Cell Biol. (2017) 36:47581. 10.1089/dna.2017.3682

  • 135.

    LiuJLiQZhangK-SHuBNiuXZhouS-Met al. Downregulation of the Long Non-Coding RNA Meg3 Promotes Angiogenesis After Ischemic Brain Injury by Activating Notch Signaling. Mol Neurobiol. (2017) 54:817990. 10.1007/s12035-016-0270-z

  • 136.

    BoulberdaaMScottEBallantyneMGarciaRDescampsBAngeliniGDet al. A role for the long noncoding RNA SENCR in commitment and function of endothelial cells. Mol Ther. (2016) 24:97890. 10.1038/mt.2016.41

  • 137.

    ZhangHXueCWangYShiJZhangXLiWet al. Deep RNA sequencing uncovers a repertoire of human macrophage long intergenic noncoding rnas modulated by macrophage activation and associated with cardiometabolic diseases. J Am Heart Assoc. (2017) 6:e007431. 10.1161/JAHA.117.007431

  • 138.

    XueCZhangXZhangHFergusonJFWangYHinkleCCet al. De novo RNA sequence assembly during in vivo inflammatory stress reveals hundreds of unannotated lincRNAs in human blood CD14+ monocytes and in adipose tissue. Physiol Genomics (2017) 49:287305. 10.1152/physiolgenomics.00001.2017

  • 139.

    HongCTontonozP. Liver X receptors in lipid metabolism: opportunities for drug discovery. Nat Rev Drug Discov. (2014) 13:43344. 10.1038/nrd4280

  • 140.

    ParikhMPatelKSoniSGandhiT. Liver X receptor: a cardinal target for atherosclerosis and beyond. J Atheroscler Thromb. (2014) 21:51931. 10.5551/jat.19778

  • 141.

    ZhangYBreevoortSRAngdisenJFuMSchmidtDRHolmstromSRet al. Liver LXRα expression is crucial for whole body cholesterol homeostasis and reverse cholesterol transport in mice. J Clin Invest. (2012) 122:168899. 10.1172/JCI59817

  • 142.

    TontonozPWuXJonesMZhangZSalisburyDSallamT. Long noncoding RNA facilitated gene therapy reduces atherosclerosis in a murine model of familial hypercholesterolemia. Circulation (2017) 136:7768. 10.1161/CIRCULATIONAHA.117.029002

  • 143.

    TilgnerHJahanbaniFBlauwkampTMoshrefiAJaegerEChenFet al. Comprehensive transcriptome analysis using synthetic long-read sequencing reveals molecular co-association of distant splicing events. Nat Biotechnol. (2015) 33:73642. 10.1038/nbt.3242

  • 144.

    Uszczynska-RatajczakBLagardeJFrankishAGuigóRJohnsonR. Towards a complete map of the human long non-coding RNA transcriptome. Nat Rev Genet. (2018) 19:53548. 10.1038/s41576-018-0017-y

  • 145.

    JainMOlsenHEPatenBAkesonM. The Oxford Nanopore MinION: delivery of nanopore sequencing to the genomics community. Genome Biol. (2016) 17:239. 10.1186/s13059-016-1103-0

  • 146.

    LagardeJUszczynska-RatajczakBSantoyo-LopezJGonzalezJMTapanariEMudgeJMet al. Extension of human lncRNA transcripts by RACE coupled with long-read high-throughput sequencing (RACE-Seq). Nat Commun. (2016) 7:12339. 10.1038/ncomms12339

  • 147.

    MercerTRClarkMBCrawfordJBrunckMEGerhardtDJTaftRJet al. Targeted sequencing for gene discovery and quantification using RNA CaptureSeq. Nat Protoc. (2014) 9:9891009. 10.1038/nprot.2014.058

  • 148.

    ClarkMBMercerTRBussottiGLeonardiTHaynesKRCrawfordJet al. Quantitative gene profiling of long noncoding RNAs with targeted RNA sequencing. Nat Methods (2015) 12:33942. 10.1038/nmeth.3321

  • 149.

    BussottiGLeonardiTClarkMBMercerTRCrawfordJMalquoriLet al. Improved definition of the mouse transcriptome via targeted RNA sequencing. Genome Res. (2016) 26:70516. 10.1101/gr.199760.115

  • 150.

    LagardeJUszczynska-RatajczakBCarbonellSPérez-LluchSAbadADavisCet al. High-throughput annotation of full-length long noncoding RNAs with capture long-read sequencing. Nat Genet. (2017) 49:173140. 10.1038/ng.3988

  • 151.

    ChuCQuinnJChangHY. Chromatin isolation by RNA purification (ChIRP). J Vis Exp. (2012) e3912. 10.3791/3912

  • 152.

    ChuCQuKZhongFLArtandiSEChangHY. Genomic maps of long noncoding RNA occupancy reveal principles of RNA-chromatin interactions. Mol Cell (2011) 44:66778. 10.1016/j.molcel.2011.08.027

  • 153.

    QuinnJJIlikIAQuKGeorgievPChuCAkhtarAet al. Revealing long noncoding RNA architecture and functions using domain-specific chromatin isolation by RNA purification. Nat Biotechnol. (2014) 32:93340. 10.1038/nbt.2943

  • 154.

    SimonMDWangCIKharchenkoPVWestJAChapmanBAAlekseyenkoAAet al. The genomic binding sites of a noncoding RNA. Proc Natl Acad Sci USA. (2011) 108:20497502. 10.1073/pnas.1113536108

  • 155.

    SimonMD. Capture hybridization analysis of RNA targets (CHART). Curr Protoc Mol Biol. (2013) Chapter 21:Unit 21.25. 10.1002/0471142727.mb2125s101

  • 156.

    LiXZhouBChenLGouL-TLiHFuX-D. GRID-seq reveals the global RNA-chromatin interactome. Nat Biotechnol. (2017) 35:94050. 10.1038/nbt.3968

  • 157.

    SridharBRivas-AstrozaMNguyenTCChenWYanZCaoXet al. Systematic Mapping of RNA-chromatin interactions in vivo. Curr Biol. (2017) 27:6029. 10.1016/j.cub.2017.01.011

  • 158.

    ChuCZhangQCda RochaSTFlynnRABharadwajMCalabreseJMet al. Systematic discovery of Xist RNA binding proteins. Cell (2015) 161:40416. 10.1016/j.cell.2015.03.025

  • 159.

    YangLLinCLiuWZhangJOhgiKAGrinsteinJDet al. ncRNA- and Pc2 methylation-dependent gene relocation between nuclear structures mediates gene activation programs. Cell (2011) 147:77388. 10.1016/j.cell.2011.08.054

  • 160.

    ChenRLiuYZhuangHYangBHeiKXiaoMet al. Quantitative proteomics reveals that long non-coding RNA MALAT1 interacts with DBC1 to regulate p53 acetylation. Nucleic Acids Res. (2017) 45:994759. 10.1093/nar/gkx600

  • 161.

    McHughCAGuttmanM. RAP-MS: a method to identify proteins that interact directly with a specific RNA molecule in cells. Methods Mol Biol. (2018) 1649:47388. 10.1007/978-1-4939-7213-5_31

  • 162.

    HuppertzIAttigJD'AmbrogioAEastonLESibleyCRSugimotoYet al. iCLIP: protein-RNA interactions at nucleotide resolution. Methods (2014) 65:27487. 10.1016/j.ymeth.2013.10.011

  • 163.

    HafnerMLandthalerMBurgerLKhorshidMHausserJBerningerPet al. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell (2010) 141:12941. 10.1016/j.cell.2010.03.009

  • 164.

    LicatalosiDDMeleAFakJJUleJKayikciMChiSWet al. HITS-CLIP yields genome-wide insights into brain alternative RNA processing. Nature (2008) 456:4649. 10.1038/nature07488

  • 165.

    ZarnegarBJFlynnRAShenYDoBTChangHYKhavariPA. irCLIP platform for efficient characterization of protein-RNA interactions. Nat Methods (2016) 13:48992. 10.1038/nmeth.3840

  • 166.

    Van NostrandELPrattGAShishkinAAGelboin-BurkhartCFangMYSundararamanBet al. Robust transcriptome-wide discovery of RNA-binding protein binding sites with enhanced CLIP (eCLIP). Nat Methods (2016) 13:50814. 10.1038/nmeth.3810

  • 167.

    McMahonACRahmanRJinHShenJLFieldsendALuoWet al. TRIBE: Hijacking an RNA-editing enzyme to identify cell-specific targets of RNA-binding proteins. Cell (2016) 165:74253. 10.1016/j.cell.2016.03.007

  • 168.

    RahmanRXuWJinHRosbashM. Identification of RNA-binding protein targets with HyperTRIBE. Nat Protoc. (2018) 13:182949. 10.1038/s41596-018-0020-y

  • 169.

    EngreitzJLanderESGuttmanM. RNA antisense purification (RAP) for mapping RNA interactions with chromatin. Methods Mol Biol. (2015) 1262:18397. 10.1007/978-1-4939-2253-6_11

  • 170.

    SharmaESterne-WeilerTO'HanlonDBlencoweBJ. Global mapping of human RNA-RNA Interactions. Mol Cell. (2016) 62:61826. 10.1016/j.molcel.2016.04.030

  • 171.

    LuZZhangQCLeeBFlynnRASmithMARobinsonJTet al. RNA duplex map in living cells reveals higher-order transcriptome structure. Cell (2016) 165:126779. 10.1016/j.cell.2016.04.028

  • 172.

    AwJGAShenYWilmASunMLimXNBoonK-Let al. In vivo mapping of Eukaryotic RNA interactomes reveals principles of higher-order organization and regulation. Mol Cell. (2016) 62:60317. 10.1016/j.molcel.2016.04.028

  • 173.

    NguyenTCCaoXYuPXiaoSLuJBiaseFHet al. Mapping RNA-RNA interactome and RNA structure in vivo by MARIO. Nat Commun. (2016) 7:12023. 10.1038/ncomms12023

  • 174.

    FeminoAMFayFSFogartyKSingerRH. Visualization of single RNA transcripts in situ. Science (1998) 280:58590. 10.1126/science.280.5363.585

  • 175.

    ShahSLubeckEZhouWCaiL. seqFISH accurately detects transcripts in single cells and reveals robust spatial organization in the hippocampus. Neuron (2017) 94:752758.e1. 10.1016/j.neuron.2017.05.008

  • 176.

    LubeckECoskunAFZhiyentayevTAhmadMCaiL. Single-cell in situ RNA profiling by sequential hybridization. Nat Methods (2014) 11:3601. 10.1038/nmeth.2892

  • 177.

    MoffittJRZhuangX. RNA Imaging with Multiplexed Error-Robust Fluorescence in situ Hybridization (MERFISH). Meth Enzymol. (2016) 572:149. 10.1016/bs.mie.2016.03.020

  • 178.

    MoffittJRHaoJBambah-MukkuDLuTDulacCZhuangX. High-performance multiplexed fluorescence in situ hybridization in culture and tissue with matrix imprinting and clearing. Proc Natl Acad Sci USA. (2016) 113:1445661. 10.1073/pnas.1617699113

  • 179.

    EngC-HLShahSThomassieJCaiL. Profiling the transcriptome with RNA SPOTs. Nat Methods (2017) 14:11535. 10.1038/nmeth.4500

  • 180.

    StåhlPLSalménFVickovicSLundmarkANavarroJFMagnussonJet al. Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science (2016) 353:7882. 10.1126/science.aaf2403

  • 181.

    AspMSalménFStåhlPLVickovicSFelldinULöflingMet al. Spatial detection of fetal marker genes expressed at low level in adult human heart tissue. Sci Rep. (2017) 7:12941. 10.1038/s41598-017-13462-5

  • 182.

    WangJZhaoYZhouXHiebertSWLiuQShyrY. Nascent RNA sequencing analysis provides insights into enhancer-mediated gene regulation. BMC Genomics (2018) 19:633. 10.1186/s12864-018-5016-z

  • 183.

    TaniHMizutaniRSalamKATanoKIjiriKWakamatsuAet al. Genome-wide determination of RNA stability reveals hundreds of short-lived noncoding transcripts in mammals. Genome Res. (2012) 22:94756. 10.1101/gr.130559.111

  • 184.

    SakuraiMUedaHYanoTOkadaSTerajimaHMitsuyamaTet al. A biochemical landscape of A-to-I RNA editing in the human brain transcriptome. Genome Res. (2014) 24:52234. 10.1101/gr.162537.113

  • 185.

    NishikuraK. A-to-I editing of coding and non-coding RNAs by ADARs. Nat Rev Mol Cell Biol. (2016) 17:8396. 10.1038/nrm.2015.4

  • 186.

    NigitaGVenezianoDFerroA. A-to-I RNA editing: current knowledge sources and computational approaches with special emphasis on non-coding RNA molecules. Front Bioeng Biotechnol. (2015) 3:37. 10.3389/fbioe.2015.00037

  • 187.

    StellosKGatsiouAStamatelopoulosKPerisic MaticLJohnDLunellaFFet al. Adenosine-to-inosine RNA editing controls cathepsin S expression in atherosclerosis by enabling HuR-mediated post-transcriptional regulation. Nat Med. (2016) 22:114050. 10.1038/nm.4172

  • 188.

    GuttmanMRinnJL. Modular regulatory principles of large non-coding RNAs. Nature (2012) 482:33946. 10.1038/nature10887

  • 189.

    SchlackowMNojimaTGomesTDhirACarmo-FonsecaMProudfootNJ. Distinctive patterns of transcription and RNA processing for human lincRNAs. Mol Cell (2017) 65:2538. 10.1016/j.molcel.2016.11.029

  • 190.

    SteijgerTAbrilJFEngströmPGKokocinskiFRGASPConsortiumHubbardTJet al. Assessment of transcript reconstruction methods for RNA-seq. Nat Methods (2013) 10:117784. 10.1038/nmeth.2714

  • 191.

    YouB-HYoonS-HNamJ-W. High-confidence coding and noncoding transcriptome maps. Genome Res. (2017) 27:105062. 10.1101/gr.214288.116

  • 192.

    FuSMaYYaoHXuZChenSSongJet al. IDP-denovo: de novo transcriptome assembly and isoform annotation by hybrid sequencing. Bioinformatics (2018) 34:216876. 10.1093/bioinformatics/bty098

  • 193.

    HarrowJFrankishAGonzalezJMTapanariEDiekhansMKokocinskiFet al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. (2012) 22:176074. 10.1101/gr.135350.111

  • 194.

    O'LearyNAWrightMWBristerJRCiufoSHaddadDMcVeighRet al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. (2016) 44(D1):D73345. 10.1093/nar/gkv1189

  • 195.

    FangSZhangLGuoJNiuYWuYLiHet al. NONCODEV5: a comprehensive annotation database for long non-coding RNAs. Nucleic Acids Res. (2018) 46(D1):D30814. 10.1093/nar/gkx1107

  • 196.

    HouMTangXTianFShiFLiuFGaoG. AnnoLnc: a web server for systematically annotating novel human lncRNAs. BMC Genomics (2016) 17:931. 10.1186/s12864-016-3287-9

  • 197.

    WucherVLegeaiFHédanBRizkGLagoutteLLeebTet al. FEELnc: a tool for long non-coding RNA annotation and its application to the dog transcriptome. Nucleic Acids Res. (2017) 45:e57. 10.1093/nar/gkw1306

  • 198.

    YangCYangLZhouMXieHZhangCWangMDet al. LncADeep: an ab initio lncRNA identification and functional annotation tool based on deep learning. Bioinformatics (2018) 34:382534. 10.1093/bioinformatics/bty428

  • 199.

    ZhouJHuangYDingYYuanJWangHSunH. lncFunTK: a toolkit for functional annotation of long noncoding RNAs. Bioinformatics (2018) 34:34156. 10.1093/bioinformatics/bty339

  • 200.

    CaoZPanXYangYHuangYShenH-B. The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier. Bioinformatics (2018) 34:218594. 10.1093/bioinformatics/bty085

  • 201.

    HuYLinJHuJHuGWangKZhangHet al. PennDiff: detecting differential alternative splicing and transcription by RNA sequencing. Bioinformatics (2018) 34:238491. 10.1093/bioinformatics/bty097

  • 202.

    KirkJMKimSOInoueKSmolaMJLeeDMSchertzerMDet al. Functional classification of long non-coding RNAs by k-mer content. Nat Genet. (2018) 50:147482. 10.1038/s41588-018-0207-8

  • 203.

    SunZNairAChenXProdduturiNWangJKocherJ-P. UClncR: Ultrafast and comprehensive long non-coding RNA detection from RNA-seq. Sci Rep. (2017) 7:14196. 10.1038/s41598-017-14595-3

  • 204.

    SchneiderHWRaiolTBrigidoMMWalterMEMTStadlerPF. A support vector machine based method to distinguish long non-coding RNAs from protein coding transcripts. BMC Genomics (2017) 18:804. 10.1186/s12864-017-4178-4

  • 205.

    ZhaoJSongXWangK. lncScore: alignment-free identification of long noncoding RNA from assembled novel transcripts. Sci Rep. (2016) 6:34838. 10.1038/srep34838

  • 206.

    HuLXuZHuBLuZJ. COME: a robust coding potential calculation tool for lncRNA identification and characterization based on multiple features. Nucleic Acids Res. (2017) 45:e2. 10.1093/nar/gkw798

  • 207.

    SunLLiuHZhangLMengJ. lncRScan-SVM: a tool for predicting long non-coding RNAs using support vector machine. PLoS ONE (2015) 10:e0139654. 10.1371/journal.pone.0139654

  • 208.

    BaekJLeeBKwonSYoonS. lncRNAnet: long Non-coding RNA identification using deep learning. Bioinformatics (2018) 34:388997. 10.1093/bioinformatics/bty418

  • 209.

    FanX-NZhangS-W. lncRNA-MFDL: identification of human long non-coding RNAs by fusing multiple features and using deep learning. Mol Biosyst. (2015) 11:8927. 10.1039/C4MB00650J

  • 210.

    KangY-JYangD-CKongLHouMMengY-QWeiLet al. CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Res. (2017) 45(W1):W126. 10.1093/nar/gkx428

  • 211.

    KanehisaMGotoSFurumichiMTanabeMHirakawaM. KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. (2010) 38(Database issue):D35560. 10.1093/nar/gkp896

  • 212.

    BudachSMarsicoA. pysster: classification of biological sequences by learning sequence and structure motifs with convolutional neural networks. Bioinformatics (2018) 34:30357. 10.1093/bioinformatics/bty222

  • 213.

    LorenzRBernhartSHHöner Zu SiederdissenCTaferHFlammCStadlerPFet al. ViennaRNA Package 2.0. Algorithms Mol Biol. (2011) 6:26. 10.1186/1748-7188-6-26

  • 214.

    ReuterJSMathewsDH. RNAstructure: software for RNA secondary structure prediction and analysis. BMC Bioinformatics (2010) 11:129. 10.1186/1471-2105-11-129

  • 215.

    MarkhamNRZukerM. UNAFold: software for nucleic acid folding and hybridization. Methods Mol Biol. (2008) 453:331. 10.1007/978-1-60327-429-6_1

  • 216.

    PedersenJSBejeranoGSiepelARosenbloomKLindblad-TohKLanderESet al. Identification and classification of conserved RNA secondary structures in the human genome. PLoS Comput Biol. (2006) 2:e33. 10.1371/journal.pcbi.0020033

  • 217.

    NelsonADLDevisettyUKPalosKHaug-BaltzellAKLyonsEBeilsteinMA. Evolinc: a tool for the identification and evolutionary comparison of long intergenic non-coding RNAs. Front Genet. (2017) 8:52. 10.3389/fgene.2017.00052

  • 218.

    GruberARNeuböckRHofackerILWashietlS. The RNAz web server: prediction of thermodynamically stable and evolutionarily conserved RNA structures. Nucleic Acids Res. (2007) 35(Web Server issue):W3358. 10.1093/nar/gkm222

  • 219.

    GesellTWashietlS. Dinucleotide controlled null models for comparative RNA gene prediction. BMC Bioinformatics (2008) 9:248. 10.1186/1471-2105-9-248

  • 220.

    TalyanSAndrade-NavarroMAMuroEM. Identification of transcribed protein coding sequence remnants within lincRNAs. Nucleic Acids Res. (2018) 46:87209. 10.1093/nar/gky608

  • 221.

    ZhangXXueCLinJFergusonJFWeinerALiuWet al. Interrogation of nonconserved human adipose lincRNAs identifies a regulatory role of linc-ADAL in adipocyte metabolism. Sci Transl Med. (2018) 10:aar5987. 10.1126/scitranslmed.aar5987

  • 222.

    Aparicio-PratEArnanCSalaIBoschNGuigoRJohnsonR. DECKO: Single-oligo, dual-CRISPR deletion of genomic elements including long non-coding RNAs. BMC Genomics (2015) 16:846. 10.1186/s12864-015-2086-z

  • 223.

    CovarrubiasSRobinsonEKShapleighBVollmersAKatzmanSHanleyNet al. CRISPR/Cas-based screening of long non-coding RNAs (lncRNAs) in macrophages with an NF-KB reporter. J Biol Chem. (2017) 292:2091120. 10.1074/jbc.M117.799155

  • 224.

    KimDHMarinovGKPepkeSSingerZSHePWilliamsBet al. Single-cell transcriptome analysis reveals dynamic changes in lncRNA expression during reprogramming. Cell Stem Cell. (2015) 16:88101. 10.1016/j.stem.2014.11.005

  • 225.

    LiuSJNowakowskiTJPollenAALuiJHHorlbeckMAAttenelloFJet al. Single-cell analysis of long non-coding RNAs in the developing human neocortex. Genome Biol. (2016) 17:67. 10.1186/s13059-016-0932-1

Summary

Keywords

long noncoding (lnc) RNAs, genomics, cardiovascular disease, atherosclerosis, gene regulation

Citation

Turner AW, Wong D, Khan MD, Dreisbach CN, Palmore M and Miller CL (2019) Multi-Omics Approaches to Study Long Non-coding RNA Function in Atherosclerosis. Front. Cardiovasc. Med. 6:9. doi: 10.3389/fcvm.2019.00009

Received

03 November 2018

Accepted

30 January 2019

Published

19 February 2019

Volume

6 - 2019

Edited by

Daniel Duerschmied, University of Freiburg, Germany

Reviewed by

Donato Santovito, Klinikum der Universität München, Germany; Ville-Petteri Makinen, South Australian Health and Medical Research Institute (SAHMRI), Australia

Updates

Copyright

*Correspondence: Clint L. Miller

This article was submitted to Cardiovascular Genetics and Systems Medicine, a section of the journal Frontiers in Cardiovascular Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics