Long Non-Coding RNAs in Multifactorial Diseases: Another Layer of Complexity.

Multifactorial diseases such as cancer, cardiovascular conditions and neurological, immunological and metabolic disorders are a group of diseases caused by the combination of genetic and environmental factors. High-throughput RNA sequencing (RNA-seq) technologies have revealed that less than 2% of the genome corresponds to protein-coding genes, although most of the human genome is transcribed. The other transcripts include a large variety of non-coding RNAs (ncRNAs), and the continuous generation of RNA-seq data shows that ncRNAs are strongly deregulated and may be important players in pathological processes. A specific class of ncRNAs, the long non-coding RNAs (lncRNAs), has been intensively studied in human diseases. For clinical purposes, lncRNAs may have advantages mainly because of their specificity and differential expression patterns, as well as their ideal qualities for diagnosis and therapeutics. Multifactorial diseases are the major cause of death worldwide and many aspects of their development are not fully understood. Recent data about lncRNAs has improved our knowledge and helped risk assessment and prognosis of these pathologies. This review summarizes the involvement of some lncRNAs in the most common multifactorial diseases, with a focus on those with published functional data.

1. Introduction

Multifactorial diseases can be defined as a group of conditions caused by the combination of genetic and environmental factors. Categorized within this group there is a variety of diseases, including cardiovascular conditions, neurological, immunological and metabolic disorders and cancer [1]. From a public health perspective, expenses on multifactorial diseases are at the top of the list of expenditures upon disease treatment in the United States ($2.10 trillion) [2] and are among the leading causes of death worldwide. These data justify the numerous studies focusing on better understanding their molecular basis to improve diagnosis and treatment options.

The Human Genome Project and high-throughput RNA sequencing (RNA-seq) technologies have revealed a surprisingly small number of protein-coding regions in humans. In fact, from our whole genome, less than 2% of the total nucleotide sequence corresponds to protein-coding genes, although most of the human genome is transcribed [3,4,5,6]. The other transcripts include: ribosomal RNA (rRNA); transfer RNA (tRNA); small nuclear RNA (snRNA); small nucleolar RNA (snoRNA); small regulatory RNAs including microRNAs (miRNAs) and piwi-RNAs (piRNA); and long non-coding RNAs (lncRNAs), including here the class of transcribed ultraconserved regions (T-UCRs) [7,8,9,10,11]. The continuous generation of RNA-seq data has made it possible to investigate normal and diseased tissues [12], with the observation that non-coding RNAs (ncRNAs) are strongly deregulated and may be important players in pathological processes [13,14].

From all ncRNA classes, miRNAs and lncRNAs have received major attention [15,16]. miRNAs have been extensively studied in the last decade and a great amount of literature is available, showing their role in many disorders [17,18,19,20].

Unlike miRNAs, lncRNAs comprise a highly diverse group of RNAs [21] with multiple functions and that often exhibit tissue-specific expression patterns [22,23,24]. By definition, lncRNAs are molecules longer than 200 nucleotides that are not translated into functional proteins [25]. Transcriptome studies have revealed that most transcripts belong to the class of lncRNAs [26] and that more molecules are yet to be identified [27]. They are classified into five groups—as sense, antisense, bidirectional, intronic, and intergenic lncRNAs—according to their relative position to a protein-coding gene [28].

In general, RNA molecules can establish specific interactions with DNA, RNA or proteins, by base pairing and/or by their tertiary structure, then forming complexes with distinct biological functions. In the nucleus, lncRNAs act to regulate gene activity through chromatin remodeling and/or by recruiting transcription factors and enhancers. Their association with splicing factors and exportation complexes also interferes with gene regulation, changing the stability and availability of mRNAs in specific circumstances [29]. In the cytoplasm, lncRNAs can act in different ways. They may target mRNA transcripts and modulate their stability, by recruiting proteins that recognize or, conversely, hiding these elements from the degradation machinery. They can also block translation by double stranding with mRNA or even promote cap-independent translation. In this way, long non-coding RNAs can regulate their neighboring genomic environment (in cis) or in distant sites of action (in trans) [30]. Genes of lncRNAs can host miRNA sequences and, in this way, these molecules may be expressed coordinately. Finally, lncRNAs can act as sponges (or competing endogenous RNAs, ceRNAs) of miRNAs and proteins, preventing them from binding to their usual targets [31]. More recently, new studies have shed light on the encoding role of lncRNAs, including the description of functional micropeptides [32,33,34,35].

LncRNAs are essential for both pluripotency in embryonic stem cells [36] and the formation and physiology of specialized cell types [37,38,39]. They are also associated with disease onset and progression, with increasing evidence for their differential expression in a variety of pathological conditions [19,28,40,41,42,43,44,45,46,47,48,49,50]. Among all multifactorial diseases, lncRNAs have already been well described in cancer [9,45,51,52,53,54,55,56,57,58]. Here, we will summarize the involvement of some lncRNAs in cardiovascular, autoimmune and neurodegenerative diseases and in the most common psychiatric disorders. This review focuses on lncRNAs with published functional data.

2. Long Non-Coding RNAs in Complexes Diseases

2.1. Cardiovascular Diseases

Cardiovascular diseases (CVD) are the most common cause of death worldwide and contribute substantially to health and economic burdens globally. Some factors contribute to cardiovascular disease development such as unhealthy behaviors (smoking, physical inactivity, inappropriate nutrition, obesity, abusive use of alcohol), some health conditions (hypertension, high cholesterol, metabolic syndromes and diabetes mellitus) and genetics factors (alleles, mutation, methylation pattern, chromosomal abnormalities) [59,60].

CVD are a group of disorders that involves heart and blood vessel diseases that are related to atherosclerosis in most cases [60,61]. Atherosclerosis is a chronic disease with inflammatory manifestation in which plaques (made of fatty substances, cholesterol, cellular waste products, fibrin and calcium) build up on walls of arteries making it harder for blood to flow [61]. Therefore, this condition is the primary origin of various cardiovascular diseases: angina, myocardial infarction, heart failure, heart attack, ischemic stroke, coronary artery diseases, and other cardiovascular problems [62].

The main events in atherosclerosis involve lipid metabolism, inflammatory response, cell proliferation and apoptosis, adhesion and migration [63]. This multifactorial process is attributed to the interactions of several key components accompanied by fatty plaque formation, a cascade of inflammatory response, proliferation of smooth muscle cells, and migration of monocytes [64]. Consequently, a cascade of immune-inflammatory responses triggers wound healing and repair.

Among the genetic contribution to CVD, lncRNAs may participate in atherosclerosis regulation but several of them remain functionally uncharacterized [65,66,67,68,69,70,71,72]. LncRNAs influence the transcriptional landscape from epigenetic chromosome modification and remodeling to the regulation of transcriptional and post-transcriptional processes, and regulate the functions of endothelial cells (ECs), smooth muscle cells (SMCs) and macrophages, and vascular inflammation and metabolism, suggesting the involvement of lncRNAs in the manifestation and progression of CVD [73].

The lncRNA ANRIL (CDKN2B antisense RNA 1—CDKN2B-AS1) regulates p16 and p15 expression affecting proliferation, development and progression of atherosclerosis [74,75]. ANRIL can interact with polycomb repressive complex 2 (PRC 2) components such as SUZ12 and CBX7, promoting the H3K27me3, resulting in the epigenetic repression of INK4α. These interactions can influence p14, p15 and p16 expression [75,76,77]. There is evidence showing that polymorphisms in ANRIL can interfere in its interaction with p14, 15, and p16, affecting the expression of these proteins and, in consequence, affecting the PRC complex role, therefore influencing progression to some diseases and contributing, for example, to coronary artery disease (CAD) susceptibility [78,79]. This lncRNA is located in 9p21.3 and a risk allele is carried by 75% of the European population that confers risk for coronary atherosclerosis [80], ischemia stroke [81], intracranial aneurysm [82,83], myocardial infarction [84], abdominal aortic aneurysms [85] and peripheral vascular disease [86]. A meta-analysis study demonstrated that this locus is significantly associated with early onset of heart disease [87].

The myosin heavy chain-associated RNA transcripts, also known as MyHeart (MHRT), is a cardioprotective antisense lncRNA encoded by the murine myosin heavy chain 7 (MYH7) and is suggested to have a conserved biological function [65]. BRG1-HDAC-PARP (SMARCA4-HDAC9-PARP1) chromatin complex inhibits MYH6 and MYH7 (MHRT) transcription. In the healthy heart, MHRT binds to BRG1 through its helicase domain, preventing BRG1 from repressing its target gene MYH6 which results in the suppression of cardiac hypertrophy and heart failure [65,88,89]. Increased levels of stress stimuli have been found to downregulate MHRT through the activation of BRG1-HDAC-PARP chromatin repressor complex, which exacerbates cardiac hypertrophy. Conversely, MHRT inhibits BRG1 and increases the levels of genes involved in cardiomyopathy [65,88,90,91]. This lncRNA was also shown to be a novel predictive biomarker of heart failure [92].

Cardiac hypertrophy-associated transcript (CHAST) is a pro-hypertrophic lncRNA that was recently described as a cis-regulator of PLEKHM1 located on the opposite strand [71]. This lncRNA may have an important role in cardiac remodeling and heart failure despite not being a cardiac-specific lncRNA, but its cell type-specific expression is sensitive to pressure overload pro-hypertrophic stimuli in vitro and in vivo. The CHAST gene has a binding site for the pro-hypertrophic NFAT (nuclear factor of activated T cells) transcription factor which negatively regulates PLEKHM1, a protein with autophagy regulatory activity [93]. In this way, CHAST can promote the inhibition of this process and its overexpression can drive cardiomyocytes hypertrophy. This lncRNA was found to be upregulated in a murine model and in patients with aortic stenosis [71,93]. However, the mechanism of CHAST/PLEKHM1 interaction it is still unknown. It seems that CHAST affects pathways of cardiac muscle morphogenesis and Wnt signaling, an additional driver of cardiac hypertrophy, but the mechanism is still unknown [94]. Bioinformatics analysis shows that CHAST can interact with vinculin, laminins and calmodulin, which have already been related to some cardiac diseases [95,96,97]. Inhibition of CHAST shows some therapeutic relevance to prevent or attenuate cardiac remodeling before and after induction of hypertrophy, respectively [71].

LIPCAR (MT-LIPCAR—mitochondrially encoded long non-coding cardiac-associated RNA) is a mitochondrial lncRNA associated with a higher risk of mortality in patients with heart failure [98]. In a study measuring plasma RNAs in patients experiencing heart remodeling following myocardial infarct (MI), and in healthy individuals by a global transcriptomic technology, LIPCAR was detected as downregulated during the initial stages of MI, whereas it was upregulated during the later stages. Patients with chronic heart failure had elevated levels of LIPCAR as compared to patients with continuous heart remodeling at one year after MI. Among a group of 344 individuals with systolic heart failure, LIPCAR levels were indeed associated with future cardiovascular health, suggesting that LIPCAR is a novel biomarker of cardiac remodeling and predicts future death in patients with heart failure [66].

CARMEN (CARMN—Cardiac mesoderm enhancer associated ncRNA) is an enhancer-associated lncRNA that has trans-repressive roles during cardiac differentiation. Ounzain et al. showed, by RNA immunoprecipitation, that this lncRNA interacts with SUZ12 and EZH2, two components of the PRC2, but the molecular mechanisms of this interaction are currently unknown. They found that the inhibition of cardiac specification and differentiation in cardiac precursor cells occurs when CARMEN is knocked down through the ablation of gene expression of cardiac transcription factors and differentiation markers, including GATA4, NKX2.5, TBX5, MYH6, MYH7, and TNNI. During pathological remodeling in the human heart, CARMEN expression is activated and necessary for maintaining cardiac identity in differentiated cardiomyocytes [68].

Cardiac hypertrophy occurs because of a heart damage compensation mechanism, but its molecular basis is still being characterized. Some evidence has demonstrated that myeloid differentiation primary response 88 (MYD88) gene has an impact on cardiac pathology [99,100]. This gene was characterized as a target of miR-489, which prevents hypertrophic responses by regulating the expression of MYD88 [67]. In the same work, Wang and colleagues identified a lncRNA named cardiac hypertrophy related factor (CHRF) that acts as an endogenous sponge for miR-489, regulating cardiac hypertrophy. CHRF induces hypertrophic responses in vitro and increases cardiomyocyte apoptosis in vivo, whereas knockdown attenuates cardiac hypertrophy in the angiotensin II-induced mouse model of heart failure [101].

A combination of increasing expression of genes related to muscle contraction, such as MYH7 and ACTA1, and genes implicated in the control of extracellular fluid volume, such as ANF, contribute to the development of cardiac hypertrophy. The lncRNA named CHAER (cardiac-hypertrophy-associated epigenetic regulator) associates with PRC2, preventing H3 lysine 27 methylation at the promoter regions of those genes. Wang et al. demonstrated that CHAER-PRC2 interaction is induced following stress stimulation, and subsequent methylation at the promoter regions of genes involved in cardiac hypertrophy contributes to the development of cardiac hypertrophy and pathological remodeling. The authors also suggest that this interaction could serve as a potential therapeutic target for the treatment of this kind of cardiovascular disease [72].

It has been verified that polymorphisms in the lncRNA MIAT confer risk of myocardial infarct [78]. In a recent publication, the same group demonstrated that MIAT is upregulated in the myocardium of diabetic rats and may act as a ceRNA to increase Dapk2 expression by sponging miR-22-3p, contributing to cardiomyocyte apoptosis [102]. Studies showed that mice that overexpress miR-150 in the heart are resistant to cardiac hypertrophy through downregulation of serum response factor (SRF). Conversely, the loss of function of miR-150 produced the opposite effects [103]. It has been determined that over-expression of Miat in H9C2 cells reduces mir-150 expression, suggesting that this lncRNA acts as a sponge for miR-150 during hypertrophy development [104]. In an AngII-induced cardiac hypertrophy mouse model, lncRNA Miat was significantly increased and necessary for cardiac hypertrophy. In addition, Yan et al. found that lncRNA Miat regulates microvascular dysfunction by functioning as a ceRNA in rat models [105]. The regulation of gene expression by recruiting miRNAs seems to be well characterized for this lncRNA, at least in animal models. In humans, it was demonstrated that MIAT is overexpressed in Chagas disease patients with inflammatory dilated chronic cardiomyopathy, compared to non-inflammatory dilated cardiomyopathy and controls, suggesting that this lncRNA may be associated with the increased severity of this pathology, operating by yet unknown mechanisms [106].

The lncRNA H19 may participate in the progression of CAD since it is expressed in human atherosclerotic plaques and in a rat model of carotid artery injury [107,108], and that hyperhomocysteinemia, an independent risk factor for CAD, could increase the expression of H19 in aorta and vascular smooth muscle cells [109]. H19 is upregulated in pathological cardiac hypertrophy and heart failure [110], and may act through different mechanisms. First, it may protect cardiomyocytes from phenylephrine-induced hypertrophy through miR-675, mapped in the H19 locus. The authors demonstrated that its overexpression upregulates the expression level of miR-675 in cardiomyocytes, which in turn inhibits hypertrophic growth of these cells. Furthermore, a pro-hypertrophic factor, Ca/calmodulin-dependent protein kinase IIδ (CaMKIId), a direct target of miR-675, partially mediated the effect of H19 on cardiomyocyte hypertrophy [69]. Second, the downregulation of H19 inhibits miR-19b, a microRNA that was previously demonstrated to promote cell proliferation and suppress apoptosis of mouse embryonic carcinoma cells, controlling Sox6 expression which, in turn, promotes cell proliferation and inhibits cell apoptosis during late-stage cardiac differentiation [111,112]. Third, the contribution of lncRNA H19 to cardiac fibroblast proliferation and fibrosis can be in part through repression of DUSP5/ERK1/2 [113]. Cardiac fibrosis depends on cardiac fibroblast activation. The main pathway involved in this activation and cell proliferation is MAPK/ERK1/2 and its regulation and repression by DUSP5 [114]. Tao et al. observed that H19 is upregulated in activated cardiac fibroblast and fibrosis tissue as well as p-ERK1/2, although the upregulation of ERK1/2 is often preceded by the dysregulation of DUSP5. So, they investigated the possible relationship between these two phenomena and found that overexpression of H19 is correlated with DUSP5 downregulation and p-ERK1/2 upregulation in cardiac fibroblast [113]. Finally, Pan suggested that H19 regulates atherosclerosis via MAPK and NF-kB signaling pathway inducing cell proliferation and reducing apoptosis of vascular smooth muscle cells [115].

The lncPPARδ has been indicated to be a good biomarker for CAD [116]. LncPPARδ regulates the expression of neighboring protein-coding gene PPARδ which in turn acts on its targets ADRP and ANGPTL4. Activation of PPARδ increases cholesterol efflux from the macrophages in the lesions and decreases trans-endothelial migration of leukocyte/monocytes into the arterial wall, thus reducing the atherosclerotic lesion size [117], indicating its important role in the process of CAD development.

Besides its roles in cancer [118,119,120], UCA1 inhibited the expression of miR-1 and protected the cardiomyocyte against H2O2-induced apoptosis mediated by caspase3 through downregulation of the tumor suppressor gene p27, providing new insights into the roles of lncRNAs in the development of ischemia/reperfusion-induced injury [121]. UCA1 was proposed as a novel biomarker for acute myocardial infarction since its plasma levels are significantly decreased 12 h after myocardial infarction but elevated after 72 h in comparison with healthy controls. This lncRNA may become a classical marker for cardiac injury such as cardiac troponin I or creatine kinase, improving the diagnose [122].

2.2. Autoimmune Diseases

Autoimmune diseases are a class of various multifactorial disorders, characterized by the exacerbated response of the immunological system against self-antigens. Several lncRNAs from loci associated with autoimmune diseases were found enriched in immune cell types. Moreover, a co-expression pattern of these lncRNAs with protein-coding genes is suggestive of a common pathway in autoimmunity [123].

The Foxp3 long intergenic non-coding RNA (FLICR) is a lncRNA neighboring the transcription factor FOXP3. FLICR negatively regulates FOXP3, decreasing its levels in regulatory T (Treg) cells, mainly in conditions of IL-2 deficiency. By promoting the decreased activity of Treg cells, FLICR is associated with the induction of autoimmunity [124]. Furthermore, the lncRNA GAS5 has been found in microarray analysis to be differentially expressed in the autoimmune diseases rheumatoid arthritis, systemic lupus erythematosus, multiple sclerosis and sarcoidosis, as well as in some infectious diseases, indicating the GAS5 role in immune functions and regulation [125]. LncRNAs involved in some autoimmune diseases are discussed in the next sections.

2.2.1. Rheumatoid Arthritis

Rheumatoid arthritis (RA) is the most common inflammatory arthritis and is a major cause of disability. In RA, CD4+ T cells, B cells and macrophages infiltrate the synovium and sometimes organize into discrete lymphoid aggregates with germinal centers [126]. The understanding of RA pathogenesis is still a challenge but may involve complex mechanisms of synovial inflammation and lncRNA dysregulation.

H19 overexpression has been reported in synovial tissue from patients with RA [127]. HOTAIR is another lncRNA that is overexpressed in RA, specifically in peripheral blood mononuclear cells (PBMCs) and exosomes purified from patient’s serum, suggesting that HOTAIR may be used as a potential biomarker for RA diagnosis [128]. The C5T1lncRNA (C5-OT1—C5 3’UTR overlapping transcript 1) is transcribed from the region between TRAF1 and C5 genes previously associated with RA pathogenesis. It is localized predominantly in the nucleus and its expression correlates with C5 mRNA levels in many tissues [129].

Analysis of T cells from RA patients and controls demonstrated some deregulated lncRNAs, such as increased expression of LOC100652951 and LOC100506036 [130]. Furthermore, an investigation of lncRNA expression in rheumatoid arthritis fibroblast-like synoviocytes (RAFLS) following the administration of quercetin found, among other things, differential expression of MALAT1, HOTAIR, MEG3, CBR3-AS1, GAS5 and YIYA. In this study, MALAT1 was found to be upregulated approximately 4-fold and knockdown of MALAT1 inhibited apoptosis in RAFLS and led to the activation of the PI3K/AKT pathway [131].

2.2.2. Type 1 Diabetes

Type 1 diabetes is a chronic disease in which genetic predisposition, coupled with environmental influence predominantly in early life, induces pancreatic beta cell autoimmunity, eventually resulting in progressive apoptotic destruction of these cells. Some research has been done on the involvement of lncRNAs in diabetes pathogenesis. A microarray analysis of a murine insulin-secreting cell line (MIN6) before and after exposition to proinflammatory cytokines (IL-1β, TNF-α, and IFN-γ) reported many lncRNAs differentially expressed in MIN6, many of which were modified by cytokine treatment: 467 transcripts were upregulated and 219 were downregulated. Four upregulated lncRNAs that displayed the most significant expression changes were evaluated in prediabetic non-obese diabetic (NOD) mice islet cells, where they were also overexpressed. These results allowed the conclusion that lncRNA overexpression sensitizes beta cells to apoptosis, probably contributing to their failure during the initial phases of the disease [132].

In addition, a similar study was performed by Sun et al., who conducted a microarray analysis in MIN6 cells exposed to the same proinflammatory cytokines from the previous study, resulting in 444 upregulated and 279 downregulated lncRNAs. A lncRNA-mRNA co-expression network identified interactions between the differentially expressed lncRNAs and mRNAs. The most enriched pathways among upregulated lncRNAs were: cytokine–cytokine receptor interaction, the NFkB signaling pathway, and toll-like receptor signaling. Regarding the downregulated lncRNAs, the most enriched were: metabolic pathways, viral carcinogenesis and protein processing in the endoplasmic reticulum [133]. Together, these results are indicative of an important role of lncRNAs in cytokine-mediated beta cell dysfunction.

2.2.3. Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is characterized by the production of pathogenic autoantibodies against nuclear autoantigens, pointed as the primary cause of tissue damage. There are diverse and heterogeneous clinical presentations in SLE patients, in which different tissues can be affected. The SLE autoimmunity arises by means of complex mechanisms involving every key facet of the immune system [134]. Purified monocytes from SLE patients and controls evidenced 29 and 60 lncRNAs upregulated and downregulated in SLE, respectively [135]. The expression of four lincRNAs was measured in PBMCs from SLE patients, patients with rheumatoid arthritis and subjects without autoimmune diseases, and it was found that linc0949 and linc0597 expression was significantly decreased in patients with SLE compared to patients with RA and control subjects. Besides, the level of linc0949 was reduced in patients with SLE who had the presence of cumulative organ damage and lupus nephritis. For linc1992 and linc3995, no significant differential expression was found [136]. More recently, the contribution of lncRNA NEAT1 for SLE pathogenesis was investigated. The NEAT1 expression was abnormally increased in monocytes from SLE patients. Additionally, it was demonstrated that NEAT1 expression was induced by lipopolysaccharide via p38 activation. Moreover, NEAT1 silencing experiments with target-specific antisense oligonucleotides resulted in significantly lower expression of a group of chemokines and cytokines, including IL-6 and CXCL10, highlighting the NEAT1 role in the SLE autoimmune response [137].

2.2.4. Psoriasis

Psoriasis is a common inflammatory disease in which skin cells multiply faster than normal, forming scales and red plaques patching the skin. Among the lncRNAs with a role in psoriasis pathogenesis, PRINS (psoriasis-associated non-protein coding RNA induced by stress) is the most well-stablished. PRINS regulates the G1P3 gene, which has an anti-apoptotic effect and is overexpressed in hyper proliferative lesional and non-lesional psoriatic epidermis. Thus, PRINS deregulation leads to a higher expression of G1P3, which may contribute to psoriasis as a result of keratinocyte proliferation in psoriatic lesion development [138]. In psoriatic patients, the uninvolved epidermis has a higher expression of PRINS, compared both healthy and psoriatic lesioned epidermis. As it is known that PRINS plays a protective role in cellular stress response, this differential expression suggests that PRINS might be involved in the susceptibility to psoriasis, and not in the psoriatic lesion maintenance [139]. In the functional studies of Szegedi et al., the phosphoprotein nucleophosmin (NPM) was found in HaCaT and NHEK cell lysates as a potential PRINS-direct interacting molecule. The authors revealed by immunohistochemistry that in psoriatic-lesioned samples, NPM is overexpressed in dividing cells of the basal layers compared with healthy and non-lesioned samples. It was also found that PRINS prevents the translocation of NPM from the nucleolus to the nucleoplasm of ultraviolet-B (UV-B) irradiated cultured keratinocytes, an incident that occurs in certain cells as a stress response after skin UV-B irradiation [140].

In association studies, the lncRNA PSORS1C3, located in PSORS1 (in 6p21.3) and near to HLA-C gene region, was also related to psoriasis susceptibility. The PSORS1C3 alleles −79C, −26C and +246 were strongly associated with psoriasis in a Caucasian population [141], and the allele 582A was associated with psoriasis vulgaris in Chinese patients [142].

Furthermore, several other lncRNAs were found as differentially expressed in psoriatic subjects in studies using large-scale analysis [143,144,145,146], and many of them are co-expressed with immune-related genes [144]. These data demonstrate the involvement of lncRNAs in this disease and raise the question of their role in other autoimmune diseases of the skin.

2.2.5. Other Autoimmune Diseases

Studies on lncRNA roles in other autoimmune diseases have also been done. In Sjogren’s Syndrome, increased expression of TMEVPG1 lncRNA (IFNG-AS1—IFNG antisense RNA 1) was detected in CD4+ T lymphocytes and positively correlated with Th1 cell proportion among CD4+ T cells [147]. In Kawasaki disease, the lncRNA PINC (pregnancy induced noncoding RNA) was found overexpressed in human umbilical vascular endothelial cells (HUVECs). It was also reported that knockdown of PINC enhanced the viability of HUVECs treated with TNF-α, increased the expression of anti-apoptotic genes and, on the other hand, reduced the expression of apoptotic genes [148]. In addition, THRIL lncRNA expression was significantly lower in Kawasaki disease patients during the acute phase than in the post-treatment convalescent phase [149]. In autoimmune thyroid disease, a single nucleotide variant Ex9b-SNP10 located in the promoter region of B cell-specific lncRNA SAS-ZFAT was associated with increased risk for disease [150]. For celiac disease (CeD), it has been demonstrated that Lnc13 (lncRNA13) regulates a subset of CeD-related inflammatory genes, that may account for inflammation observed, contributing to its pathogenesis [151].

2.3. Neurodegenerative Disorders

Neurodegenerative disorders are characterized by neuronal damage and reduction of brain capacity, ability to think, memory and motor coordination [152]. Many of these diseases have an accumulation of misfolded and different toxic proteins, which contribute to neuronal death [153]. In addition to the role of many genes involved in regulation, production, and degradation of these malfunctioning proteins inside the neurons, part of the neurodegeneration is caused by the involvement of other cells and inflammation components [154]. LncRNA deregulation has also been included as part of these processes.

2.3.1. Alzheimer’s Disease

In Alzheimer’s disease (AD), mutations in one of three genes—the amyloid precursor protein (APP) gene and the presenilin protein 1 and 2 genes (PS1 and PS2), which encode subunits of gamma-secretase—contribute to the development of an early onset disease [155], and the allele E4 of the APOE gene increases considerably the risk of development of AD [156]. Toxic proteins in AD are the result of APP processing in amyloid beta peptides (Aβ) and hyperphosphorylation of microtubule-associated TAU protein. Aβ peptides are exported and accumulated outside the neuron, interfering with synapses. The accumulation of abnormal TAU protein causes toxic tangles inside the neuron and facilitates cell death [153]. The process of Aβ deposition is connected to the development of AD symptoms and it has been reported that lncRNA regulates this processing pathway.

BACE1 antisense (BACE1-AS) is an lncRNA which regulates its host gene-encoding β-site APP-cleaving enzyme 1 (BACE1) [157]. In AD, its upregulation drives an increase of Aβ production, by increasing the stability of BACE1 mRNA through the formation of RNA duplexes. Through the silencing of BACE1-AS in SH-SY5Y human neuroblastoma cells treated with exogenous Aβ, a significant reduction of cytotoxic effect was observed and these cells maintained their normal state, suggesting that BACE1-AS is critical in AD development [158].

Another lncRNA that regulates its host gene is 17A, located in the intron 3 of the GPR51 gene coding for the GABA B2 receptor, whose mRNA was found upregulated in the brain of AD. 17A expression is induced in response to inflammatory stimuli and is reduced by the treatment with non-steroidal anti-inflammatory drugs, such as diclofenac. They also demonstrated that 17A controls the alternative splicing of GPR51 mRNA; decreases transcription of the canonical form of GABA-B R2, impairing GABA-B signaling; and enhances the secretion of Aβ peptide in SH-SY5Y cells treated with IL-1α [159].

Deregulated cholesterol metabolism is a risk factor for the development of AD, interfering with Aβ production. The lncRNA n336934 is upregulated in AD brains and associated with the cholesterol homeostasis pathway [160]. In the same study, the authors also identified n341006 downregulated in AD patients brain and associated with protein ubiquitination and proteasome system pathway, affecting the toxic Aβ degradation.

Magistri et al. investigated the expression of RNAs in the brain from AD patients and age-matched controls. Four lncRNAs (HAO2-AS, EBF3-AS, AD-linc1 and AD-linc2) were reported as having a differential expression. AD-linc2 was not detected in the normal brain, though this lncRNA was found upregulated in specific regions of the brain (superior frontal gyrus, entorhinal cortex, hippocampus, and cerebellum) of AD patients. The authors investigated if the differential expression was caused by Aβ accumulation. Indeed, neuronal damage and the level of AD-linc1 were increased in cells treated with exogenous Aβ [161].

2.3.2. Parkinson’s Disease

Parkinson’s disease (PD) is characterized by the accumulation of misfolded α-synuclein protein, which forms aggregates called Lewy bodies affecting the neuron viability [162]. Clinically PD is characterized by stilly shacking, bradykinesia, muscular rigidity and abnormal posture [163]. In a few percentage of cases, PD behaves as a monogenic condition, with mutations in one of the genes involved in the accumulation of α-synuclein and coding for components of the ubiquitin-proteasome system: α-synuclein (SNCA/PARK1), leucine-rich repeat kinase 2 (LRRK2/PARK8), ubiquitin carboxyl-terminal esterase L1 (UCHL1), Parkinson protein 2 (Parkin/PARK2), PTEN-induced kinase 1 (PINK1/PARK6), PD 7 (DJ-1/PARK7) and PD 9 (ATP13A2/PARK9) [152,164].

The α-synuclein pathway is regulated by the lncRNA MALAT1, which is highly expressed in neurons and regulates a subset of genes involved in the nervous system development and in α-synuclein production. The level of MALAT1 was increased in SH-SY5Y cells after exposure to MPP (1-methyl-4-phenylpyridinium—that induces neurodegeneration to dopamine neurons) and decreased when treated with the immunosuppressant β-asarone (cis-2,4,5-trimethoxy-1-allyl phenyl—a neuroprotective drug). Additionally, MALAT1 overexpression upregulated α-synuclein expression, while inhibition of MALAT1 downregulated its expression, by a suggested mechanism in which MALAT1 bound α-synuclein protein enhancing its stability [165].

Kraus et al. observed lncRNAs (MALAT1, HUC1 and 2, lincRNA-p21, SNHG1 and NEAT1) deregulated in different parts of the brain. HUC 1 and 2 were found downregulated, while the others were upregulated in PD compared with controls. MALAT1, lincRNA-p21 and small nucleolar RNA host gene 1 (SNHG1) may be involved in apoptosis pathways and the molecular pathogenesis of PD may be associated with activation of apoptotic events by the p53 pathway [166].

Another lncRNA that was found to be overexpressed in PD patients is AL049437. The reduction in AL049437 level increased cell viability, mitochondrial transmembrane potential, mitochondrial mass, and tyrosine hydroxylase secretion. The lncRNA AK021630, downregulated in the PD brain, is involved in the same pathways, while its overexpression increases cell survival [167].

Mutations in the multifunctional protein LRRK are the major cause of familial and sporadic PD. LRRK2 may interact with the α-synuclein and tau pathway, contributing to PD development, although the mechanism has not yet been defined [168]. Two lncRNAs located in the LRRK2 locus, UC001lva.4 and AC079630, were found to be deregulated in cerebrospinal fluid (CSF) from patients with PD, however LRRK2 expression was unaltered. The roles of these two lncRNA are still unclear [169]. The lncRNA HOTAIR is also important to the development of PD. In an animal model, both HOTAIR and LRRK2 were upregulated and an overexpression of HOTAIR increased the stability and level of LRRK2 mRNA, while the knockout of HOTAIR enhanced the cell viability in SH-SY5Y cells [170].

2.3.3. Huntington’s Disease

LncRNAs are also studied in Huntington’s disease (HD) models. HD is a progressive brain disorder that causes uncontrolled movements, emotional problems, and loss of thinking ability (cognition) caused by trinucleotide (CAG) repeat expansion in the huntingtin (HTT) gene, resulting in a toxic huntingtin protein [153]. The first study with lncRNA in HD was performed by Johnson and colleagues. They investigated the differential expression of REST (RE1-silencing transcription factor) targets, affected by mutant huntingtin and found two lncRNAs (HAR1F and HAR1R), located in the HAR1 locus, downregulated in the striatum of HD patients [171].

In an animal model for HD, the LINC00035, also called ABHD11-AS1 (abhydrolase domain containing 11 antisense) in humans and Abhd11os in mice, was investigated by Francelle and colleagues. They showed that the level of Abhd11os was decreased in HD mouse models, that Abhd11os overexpression protects against neuronal toxicity, and that its knockdown exacerbates the toxicity caused by the mutant huntingtin [172].

2.3.4. Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease of motor neurons. About 10% of ALS cases are inherited, out of which 20% are caused by mutations in the gene encoding superoxide dismutase 1 (SOD1), which aggregates and deposits inside the neuron, causing death [153].

In ALS, NEAT has been found to be predominantly expressed in spinal motor neurons in an early phase of the ALS pathological process. NEAT1 is modestly expressed in the human brain and the transcript NEAT1_2 is barely detectable in the central nervous system, while it was found in the brain of ALS patients. NEAT1_2 acts as a scaffold and changes the functions of TDP-43 and FUS/TLS proteins already associated with ALS and involved in RNA processing and stress granule formation [173].

2.4. Psychiatric Disorders

Psychiatric disorders include several mental illnesses which involve psychological and/or neurological disabilities. To date, it is known that certain lncRNAs are involved in the pathophysiology of common psychiatric diseases.

Copy number variations (CNVs), especially deletions upstream the NRXN1 gene, were identified to affect the lncRNA AK127244, segregating with psychiatric conditions in affected individuals of a family [174]. Furthermore, translocations in the disrupted in schizophrenia (DISC) genome locus, including DISC1 (which encodes the scaffold protein DISC1) and DISC2, a lncRNA antisense to DISC1, have been linked with various psychiatric disorders, such as schizophrenia (SZ), autistic spectrum disorders (ASD), major depression disorder (MDD), among others [174,175]. DISC2 is involved in SZ pathogenesis possibly by regulating DISC1 expression [175,176].

2.4.1. Schizophrenia

Schizophrenia (SZ) is a mental disorder characterized by hallucinations, delusions and loss of reality. Several lncRNAs were found to be deregulated in SZ. In a study with murine cortical neurons, neuronal depolarization was induced, and then transcriptional changes were evaluated at different time points in response to neuronal activity, evidencing that MIAT was the most downregulated lncRNA [177]. MIAT binds to splicing factors QKI and serine/arginine-rich splicing factor 1 (SFSF1), and its inhibition leads to increased levels of splice variants of DISC1 and ERBB4, both upregulated genes in post-mortem schizophrenic brains [178]. Furthermore, in an association study in a Chinese population, the MIAT tag SNPs rs1894720 and rs4274 were associated with paranoid SZ [179].

The 1p21.3 loci have been associated with SZ in genome wide association studies (GWAS). Located in this genomic region, miR-137 and (dihydropyrimidine dehydrogenase (DPYD) have both been previously associated with SZ and neurological and psychiatric disorders [180,181]. The lncRNA linc01930 (EU358092) was also characterized in this SZ-associated region, being co-expressed with miR-137 and acting as a regulatory transcript, with alteration of expression in response to psychoactive drug administration [182].

The pas-lncRNA is a lncRNA antisense to the promoter of CSMD1, a well-known gene associated with SZ. Pas-lncRNA expression is co-regulated with CSMD1, specifically in the central nervous system (CNS) tissue, and presents a high expression in peripheral tissues where CSMD1 is downregulated, suggesting a brain-specific promoter activity regulation in CNS and a relevance of these genes in brain function [183].

Some of the symptoms present in the early and acute phases of SZ are depression and anxiety. The presence of these common psychiatric states in SZ is a challenge in disease diagnosis and choice of appropriate treatment, indicating the importance of distinguishing SZ from depression and anxiety. Cui and colleagues (2017) found, in a microarray analysis, lncRNAs NONHSAT089447, NONHSAT021545 and NONHSAT041499 to be aberrantly expressed in SZ patients, besides another nine lncRNAs with also different expression patterns between SZ, depression and anxiety (see topic below), suggesting these lncRNAs as potential biomarkers for the differential diagnosis of these diseases [184]. Accordingly, NONHSAT089447 and NONHSAT041499 were also found to be upregulated in SZ patients in a previous study. Moreover, these lncRNA levels were significantly reduced after psychotropic treatment, and NONHSAT041499 was negatively correlated with better treatment outcomes [185].

2.4.2. Depression and Anxiety Disorders

Major depression disorder (MDD), which refers to symptoms of depression, is a psychiatric condition mainly characterized by sadness, tiredness, guilt and feelings of low self-esteem, among other mental symptoms which vary among MDD patients, such as suicidal thoughts and anxiety conditions. The anxiety disorders include anxiety, phobias, post-traumatic stress and obsessive-compulsive conditions.

As cited in the last topic, some lncRNAs were found as potential biomarkers to distinguish SZ from MDD and anxiety. Among them, the lncRNAs TCONS_00019174, ENST00000566208, NONHSAG045500, ENST00000517573, NONHSAT034045 and NONHSAT142707 were significantly downregulated in MDD patients in comparison with controls, while ENST00000505825, NONHSAG017299 and NONHSAT078768 were upregulated in anxiety patients [184].

Among several genes differentially expressed in the medial prefrontal cortex (mPFC) of a group of fear-conditioned mice, the lncRNA Miat was significantly downregulated. The fear-related anxiety was further enhanced after Miat knockdown. The authors propose that Miat plays a role in anxiety-like behavior modulation through interaction with Bmi1, a key member of the polycomb repressive complex 1, which epigenetically regulates the expression of the schizophrenia-related gene crystallin beta 1 (Crybb1) [186].

In a genome-wide study of DNA methylation patterns in the cortical brain tissues of MDD suicidal individuals and non-psychiatric controls, in MDD suicidal samples was identified a significant hypomethylation upstream of the lncRNA gene PSORS1C3 [187]. Some co-methylated loci modules associated with propensity to suicide, which were enriched in important gene networks related to suicidality and depression, were also found [187]. This might explain the increased risk for depression and suicide in psoriatic patients [188].

In another recent study, the expression of the Linc01108 was found to be significantly higher in MDD patients in comparison with controls, and the rs12526133 SNP present in this gene also associated with the susceptibility to the disease. Moreover, the rs2272260 SNP, present in the lncRNA LINC00998, which is in turn under expressed in patients compared to controls, is also associated to MDD [189].

The role of the brain-derived neurotrophic factor (BDNF) in some psychiatric disorders, such as MDD and SZ, is well stablished [174]. The antisense to the BDNF gene is BDNF-AS. It was shown that BDNF-AS silencing, highly increases BDNF expression, upregulating glial-derived neurotrophic factor (GDNF) and ephrin receptor B2 (EPHB2) and, then, neuronal outgrowth and differentiation [190]. Moreover, it is suggested that BDNF-AS may inhibit, by complementarity, the binding between the 3’UTR of BDNF and miR-124a, a miRNA associated with MDD and neurobehavioral deficits [191,192].

Among the 2007 MDD-differentially expressed lncRNAs found by Liu et al. in a microarray analysis, the lncRNAs located at chr10:874695-874794, chr10:75873456-75873642, and chr3:47048304-47048512 were found to be able to regulate four differentially expressed mRNAs [193]. Moreover, the lncRNAs BACE-AS [157] and BC200 [194] also seem to be involved in depression [174].

2.4.3. Autistic Spectrum Disorders

The autistic spectrum disorders (ASD), which include autism, Asperger’s syndrome and other mental-related diseases, are common mental illnesses diagnosed in childhood, which involve impairments in social interaction and in communication, intellectual disabilities and repetitive behaviors.

In the cluster 5p14.1, to where important common variants associated with ASD map, the moesin pseudogene 1 (MSNP1) and its antisense lncRNA MSNP1-AS may also be found. MSNP1-AS overexpression was found in post-mortem ASD cerebral cortex and in individuals with the rs4307059 ASD risk allele [195,196]. The overexpression of this lncRNA decreased neurite number and length in human neuronal progenitor cells and dysregulated the expression of genes involved in protein synthesis and chromatin remodeling. MSNP1AS knockdown altered the expression of several genes, mainly those involved in chromatin remodeling and immune response [196].

In a study with a 3-year-old male with ASD, a 2Mb deletion in 1q42 was reported. This genomic region contains the coding genes translin-associated factor X (TSNAX) and DISC1, and the lncRNA DISC2, all of them previously associated with psychiatric disorders. Authors suggest that DISC2 haploinsufficiency may lead to a further dysregulation of the remaining DISC1 in ASD patients [197]. Moreover, in another family study, some deletions in PTCHD1AS1 and PTCHD1AS2, antisense lncRNAs which seem to regulate its overlapping coding gene PTCHD1, were associated with ASD in males [198]. Furthermore, analysis of post-mortem prefrontal cortex and cerebellum tissue of ASD patients indicated some lncRNAs as differentially expressed and located around the human imprinted loci C9orf85, SLC4A2, and UBE3A. The UBE3A gene is related with Angelman syndrome, a disorder with some ASD symptoms [199]. LncRNAs antisense to BDNF (BDNF-AS) and SHANK2 (SHANK2-AS), as some lncRNAs located within the HOX loci, were also associated with ASD [200].

3. Conclusions

The discovery of lncRNAs as another intricate layer in multifactorial disease development has generated more questions than answers since they are required for the proper function of higher eukaryote metabolism. Here we summarize the most relevant lncRNAs described until now, focused on those with published functional data (Table 1). The databases of lncRNA sequences in the human genome are still under construction and complete annotation in different cells and tissues is necessary for subsequent work. For clinical purposes, lncRNAs may have advantages mainly because of their specificity and differential expression patterns, as well as their ideal qualities for diagnosis and therapeutics. Besides, the potential to detect distinct circulating lncRNAs between patients and healthy individuals may allow the development of new diagnostic inputs, increasing the accuracy of the tests already available and helping risk assessment and prognosis.

Table 1

List of the most relevant lncRNAs and the associated mechanisms.

LncRNA	Expression	Disease	Mechanism	Reference
ANRIL	Down	Coronary Artery Disease (CAD), Myocardial Infarction	Chromatin remodelling (interacts with PRC1 and PRC2); regulation of CDKN2A and CDKN2B expression.	[74,75,76,77,79]
MHRT	Up	Heart failure, Acute Myocardial Infarction, Hypertrophy, Cardiomyopathy	Chromatin remodelling (cardiac stress); interacts with BRG1 preventing its ligation to targets.	[65,88,89,90,91,92]
CHAST	Up	Heart failure, Aortic Stenosis, Cardiac Hypertrophy.	Cardiac remodelling; negatively regulates Plekhm1.	[71,93,94,95,96,97]
LIPCAR	Down	Initial Myocardial Infarction, Chronic Cardiovascular Disease	Cardiac remodeling; unknown mechanism.	[66,98]
CARMEN	Up	Cardiovascular Disease (Cardiac Hypertrophy)	Chromatin remodelling; interacts with SUZ12 and EZH2 (PRC2).	[68]
CHRF	Up	Cardiac Hypertrophy	Sponge for mir-489.	[67,101]
CHAER (mouse)	Up	Cardiac Hypertrophy	Epigenetic Regulator; inhibits H3L27 methylation.	[72]
MIAT	Up	Cardiac Hypertrophy, Myocardial Infarction	Sponge for mir-150.	[78,102,103,104,105,106]
H19	Up	Atherosclerosis, Heart failure, Cardiac Hypertrophy	Epigenetic Regulator; suppresses miR-675 and miR-19b; repression of DUSP5/ERK1/2; regulating MAPK and NF-kB signalling pathway.	[69,107,108,109,110,111,112,113,114,115]
LncPPARδ	Up	Atherosclerosis, Cardiovascular Disease	Decreases cholesterol efflux and increases migration of leukocyte/monocytes into the arterial wall.	[116]
UCA1	Down/Up (Biphasic)	Myocardial Infarction (after 12 h/72 h)	Unknown.	[122]
H19	Up	Rheumatoid Arthritis	Unknown.	[127]
HOTAIR	Up	Rheumatoid Arthritis	Unknown.	[128]
C5T1lncRNA	Up	Rheumatoid Arthritis	Regulates C5 mRNA expression.	[129]
LOC10065295/LOC100506036	Up	Rheumatoid Arthritis	Unknown.	[130]
MALAT 1	Up	Rheumatoid Arthritis	Promotes the apoptosis of rheumatoid arthritis fibroblast-like synoviocytes and leads to the activation of the PI3K/AKT pathway.	[131]
NEAT1	Up	Systemic Lupus Erythematosus (SLE)	Contributes to expression of a group of chemokines and cytokines, including IL-6 and CXCL10.	[137]
TMEVPG1	Up	Sjogren’s Syndrome	Positively correlates expression with Th1 cell proportion among CD4+ T cells.	[147]
PINC	Up	Kawasaki disease	Promotes expression of apoptosis genes in human umbilical vascular endothelial cells (HUVEC).	[148]
PRINS	Up	Psoriasis	Interacts with NPM protein, which is associated with keratinocytes proliferation.	[139]
BACE1-AS	Up	Alzheimer’s Disease	Increases the stability of BACE1 mRNA. Aβ production.	[157,158]
17A	Up	Alzheimer’s Disease	Alternative splicing of GPR51 mRNA, reducing canonical form of GABAB R2 and impairing GABAB signaling. Aβ production.	[159]
HAO2-AS, EBF3-AS, AD-linc1 and AD-linc2	Up	Alzheimer’s Disease	Neurotoxicity.	[161]
n341006	Down	Alzheimer’s Disease	Acting in the ubiquitin-proteasome system (UPS), affecting the turnover and degradation, is impaired in AD.	[160]
n336934	Up	Alzheimer’s Disease	Cholesterol pathway.	[160]
MALAT1	Up	Parkinson’s Disease	α-synuclein production.	[165]
HOTAIR	Up	Parkinson’s Disease	Increasing the stability and level of LRRK2 mRNA.	[170]
AL049437	Up	Parkinson’s Disease	Apoptosis pathway.	[167]
AK021630	Down	Parkinson’s Disease	Apoptosis pathway.	[167]
p21	Up	Parkinson’s Disease	Targets p53 and H1F1 (apoptosis pathway).	[166]
SNHG1	Up	Parkinson’s Disease	Affects the p53 stability (cellular proliferation).	[166]
HAR1F and HAR1R	Down	Huntington’s Disease	Unknown.	[171]
ABHD11-AS1	Down	Huntington’s Disease	Affects the neuronal toxicity.	[172]
NEAT1_2	Up	Amyotrophic Lateral Sclerosis (ALS)	Changes the TDP-43 and FUS/TLS functions, stress granules formation).	[173]
DISC2	Down	Schizophrenia, Autistic Spectrum Disorders (ASD)	DISC2 regulates DISC1 expression. Translocation in both genes are involved with psychiatric disorders.	[174,175,176]
EU358092	?	Schizophrenia	Located in SZ- associated region 1p21.3 and co-expressed with mir-137; alters response to psychoactive drugs.	[182]
MIAT	Down	Schizophrenia, Fear-related anxiety	Increases pathogenic splice variants of DISC1 and ERBB4. Associated with risk SNPs rs1894720 and rs4274/Possibly interacts with BMI1, regulating Crybb1 expression.	[177,178,179]
BDNF-AS	?	Depression	Repress BDNF expression, and neuronal outgrowth and differentiation.	[191,192]
MSNP1AS	Up	ASD	Decreases neurite number and length in neuronal progenitor cells; dysregulates the expression of genes involved in protein synthesis and chromatin remodelling.	[196]

The group of diseases are classified by colors: red for cardiovascular; yellow for autoimmune and green for neurological diseases; blue for psychiatric disorders. Up and Down are related to expression levels compared with healthy patients. A question mark means “information not known”.