Long non-coding RNAs in immune regulation and their potential as therapeutic targets

Expression profile of long-noncoding RNAs MIR31HG, NKILA, and PACER in systemic lupus erythematosus patients

OBJECTIVES

Growing evidence suggests that systemic lupus erythematosus (SLE), an organ-damaging systemic autoimmune illness, may be influenced by long-noncoding RNAs (lncRNAs). This study aimed to assess the relative expression of lncRNAs (MIR31HG, NKILA, and PACER) in patients with SLE to evaluate their role in the disease.

DESIGN AND METHODS

This study involved 70 patients with SLE and 70 apparently healthy control subjects. The expression levels of lnc-MIR31HG, NKILA, and PACER were quantified using real-time PCR.

RESULTS

Lnc-MIR31HG, NKILA, and PACER were significantly upregulated in SLE cases compared to controls (P < 0.001). ROC curve analysis revealed a 91.43 % sensitivity of PACER for the diagnosis of SLE at a cutoff point of > 1.46, followed by NKILA with 90 % sensitivity at a cutoff point of > 1.16, and MIR31HG with 85.71 % sensitivity at a cutoff point of > 1.43. MIR31HG had the highest sensitivity for the diagnosis of lupus nephritis (86.67 %) at a cutoff point of > 7.19, then NKILA with 80 % sensitivity at a cutoff point of > 8.12, and finally PACER expression with 73.33 % sensitivity at a cutoff point of > 18.19. Moreover, MIR31HG and NKILA revealed a significant correlation with albumin/creatinine ratio, estimated glomerular filtration rate, and the SLEDAI score. Regression analysis revealed the potential roles of MIR31HG, NKILA, and PACER expression as predictors for SLE.

CONCLUSION

An upregulated lncRNA panel (MIR31HG, NKILA, and PACER) could play a role in the pathogenesis and, hence, the predispositiontoSLE. MIR31HG and NKILA can serve as prognostic markers significantly linked with disease activity.

Dysregulation of LncRNAs ANRIL, MALAT1, and LINC00305 in Coronary Slow Flow Patients: Implications for Inflammation and Endothelial Dysfunction

Coronary Slow Flow (CSF) is observed in individuals who experience delayed blood supply in the coronary arteries. Inflammation and endothelial dysfunction may play a role in the etiology and development of CSF. The current investigation aimed to compare the expression of specific long noncoding RNAs (lncRNAs) associated with endothelial dysfunction and inflammation in CSF patients. This case‒control study enrolled 72 CSF patients and 71 healthy individuals. Blood samples were collected, and serum marker levels were measured. The expression levels of lncRNAs ANRIL, MALAT1, and LINC00305 in peripheral blood mononuclear cells (PBMCs) were assessed using real-time Polymerase Chain Reaction (PCR). All statistical analyses were performed using SPSS 22, with the significance level set at P < 0.05. The study revealed that the relative expression of MALAT1 and LINC00305 was significantly lower in the CSF group (p < 0.01), whereas ANRIL was expressed at higher levels (p < 0.0001). The areas under the ROC curves (AUCs) for MALAT1, LINC00305, and ANRIL were 0.64, 0.66, and 0.75, respectively. Notably, the expression level of LINC00305 exhibited an inverse correlation with CSF incidence (OR: 0.83, p: 0.008) in contrast to that of ANRIL (OR: 1.43, p < 0.0001). Additionally, compared to those in the control group, the average BMI, WBC, RBC, Hb, LDH, LDL, FBS, and percentage of neutrophils in the CSF group were significantly greater (p< 0.05). lncRNA ANRIL is upregulated in CSF patients, whereas MALAT1 and LINC00305 are downregulated. Dysregulation of ANRIL, MALAT1, and LINC00305 may serve as diagnostic and predictive factors for CSF leakage.

Long Non-Coding RNA TUG1 Gene Polymorphism and TUG1 Expression Level as Molecular Biomarkers of Systemic Lupus Erythematosus and Lupus Nephritis

Long non-coding RNA (lncRNA) TUG1 acts as a proto-oncogene, allowing the proliferation of tumor cells, and it has been related to inflammation. Therefore, we aimed in this study to investigate for the first time the role of TUG1 gene polymorphism and the TUG1 level as biomarkers in systemic lupus erythematosus (SLE) and their link to lupus nephritis 145 SLE. A total of 145 healthy controls were subjected to clinical and laboratory evaluation. The disease activity was assessed by the SLE disease activity index (SLEDAI) score. SLE patients were divided into two subgroups according to the presence of lupus nephritis. The TUG1 gene polymorphisms rs5749201 and rs886471 were determined by Sanger sequencing, and TUG1 expression was assessed by qRT-PCR. There was a significant increase in the risk of SLE AA, TA, dominant genotypes, and the A allele of rs5749201 (p < 0.001) by 4.9-, 10.1-, 6.5-, and 2.5-fold in comparison to the relative control. GG and TG, dominant genotypes and the G allele of rs886471 (p < 0.01) increased the risk by 5.09-, 11.9-, 6.5-, and 2.6-fold. AA, A allele, dominant and recessive rs5749201genotypes increased the risk of lupus nephritis by 16.6-, 7.4-, 7.1-, and 12.2-fold, respectively (p < 0.05). GG, dominant and recessive genotypes, and the G allele of rs886471 increased the risk of lupus nephritis by 17.04-, 7.8-, 9.4-, and 6.08-fold, respectively (p < 0.05). Additionally, the AG haplotype increased the risk of SLE and lupus nephritis by 2.7- and 7.8-fold, respectively. The AA rs5749201 and GG rs886471 variants are significantly associated with more severe disease (p < 0.001). TUG1 expression was significantly higher in SLE than in the control and in the lupus nephritis than in non-lupus nephritis cases (p < 0.05). Interestingly, AA rs5749201 and GG rs886471 were significantly associated with higher TUG1 levels (p < 0.001). It was also found that AA rs5749201 and high SLEDAI were predictors of lupus nephritis. Overall, our findings illustrated for the first time that TUG1 gene rs5749201 and rs886471 variants were associated with increased risk of SLE, more severe disease, and lupus nephritis, and the TUG1 level could be used as a diagnostic biomarker of SLE and lupus nephritis.

A genome-wide CRISPR activation screen identifies SCREEM a novel SNAI1 super-enhancer demarcated by eRNAs

The genome is pervasively transcribed to produce a vast array of non-coding RNAs (ncRNAs). Long non-coding RNAs (lncRNAs) are transcripts of >200 nucleotides and are best known for their ability to regulate gene expression. Enhancer RNAs (eRNAs) are subclass of lncRNAs that are synthesized from enhancer regions and have also been shown to coordinate gene expression. The biological function and significance of most lncRNAs and eRNAs remain to be determined. Epithelial to mesenchymal transition (EMT) is a ubiquitous cellular process that occurs during cellular migration, homeostasis, fibrosis, and cancer-cell metastasis. EMT-transcription factors, such as SNAI1 induce a complex transcriptional program that coordinates the morphological and molecular changes associated with EMT. Such complex transcriptional programs are often subject to coordination by networks of ncRNAs and thus can be leveraged to identify novel functional ncRNA loci. Here, using a genome-wide CRISPR activation (CRISPRa) screen targeting ∼10,000 lncRNA loci we identified ncRNA loci that could either promote or attenuate EMT. We discovered a novel locus that we named SCREEM (SNAI1 cis-regulatory eRNAs expressed in monocytes). The SCREEM locus contained a cluster of eRNAs that when activated using CRISPRa induced expression of the neighboring gene SNAI1, driving concomitant EMT. However, the SCREEM eRNA transcripts themselves appeared dispensable for the induction of SNAI1 expression. Interestingly, the SCREEM eRNAs and SNAI1 were co-expressed in activated monocytes, where the SCREEM locus demarcated a monocyte-specific super-enhancer. These findings suggest a potential role for SNAI1 in monocytes. Exploration of the SCREEM-SNAI axis could reveal novel aspects of monocyte biology.

Predicting prognosis and immunotherapy response among colorectal cancer patients based on a tumor immune microenvironment-related lncRNA signature

Long non-coding RNAs (lncRNAs) remodel the tumor immune microenvironment (TIME) by regulating the functions of tumor-infiltrating immune cells. It remains uncertain the way that TIME-related lncRNAs (TRLs) influence the prognosis and immunotherapy response of colorectal cancer (CRC). Aiming at providing survival and immunotherapy response predictions, a CRC TIME-related lncRNA signature (TRLs signature) was developed and the related potential regulatory mechanisms were explored with a comprehensive analysis on gene expression profiles from 97 immune cell lines, 61 CRC cell lines and 1807 CRC patients. Stratifying CRC patients with the TRLs signature, prolonged survival was observed in the low-risk group, while the patients in the high-risk group had significantly higher pro-tumor immune cells infiltration and higher immunotherapy response rate. Through the complex TRLs-mRNA regulation network, immunoregulation pathways and immunotherapy response pathways were found to be differently activated between the groups. In conclusion, the CRC TRLs signature is capable of making prognosis and immunotherapy response predictions, which may find application in stratifying patients for immunotherapy in the bedside.

Bronchial epithelium epithelial-mesenchymal plasticity forms aberrant basaloid-like cells in vitro

Although epithelial-mesenchymal transition (EMT) is a common feature of fibrotic lung disease, its role in fibrogenesis is controversial. Recently, aberrant basaloid cells were identified in fibrotic lung tissue as a novel epithelial cell type displaying a partial EMT phenotype. The developmental origin of these cells remains unknown. To elucidate the role of EMT in the development of aberrant basaloid cells from the bronchial epithelium, we mapped EMT-induced transcriptional changes at the population and single-cell levels. Human bronchial epithelial cells grown as submerged or air-liquid interface (ALI) cultures with or without EMT induction were analyzed by bulk and single-cell RNA-Sequencing. Comparison of submerged and ALI cultures revealed differential expression of 8,247 protein coding (PC) and 1,621 long noncoding RNA (lncRNA) genes and revealed epithelial cell-type-specific lncRNAs. Similarly, EMT induction in ALI cultures resulted in robust transcriptional reprogramming of 6,020 PC and 907 lncRNA genes. Although there was no evidence for fibroblast/myofibroblast conversion following EMT induction, cells displayed a partial EMT gene signature and an aberrant basaloid-like cell phenotype. The substantial transcriptional differences between submerged and ALI cultures highlight that care must be taken when interpreting data from submerged cultures. This work supports that lung epithelial EMT does not generate fibroblasts/myofibroblasts and confirms ALI cultures provide a physiologically relevant system to study aberrant basaloid-like cells and mechanisms of EMT. We provide a catalog of PC and lncRNA genes and an interactive browser (https://bronc-epi-in-vitro.cells.ucsc.edu/) of single-cell RNA-Seq data for further exploration of potential roles in the lung epithelium in health and lung disease.

LINC01123 promotes immune escape by sponging miR-214-3p to regulate B7-H3 in head and neck squamous-cell carcinoma

Numerous studies have shown that long noncoding RNAs (LncRNAs) are involved in the development and immune escape of head and neck squamous-cell carcinoma (HNSCC). However, the specific regulatory mechanisms by which LINC01123 regulates HNSCC and its correlation with immunity remain unclear. Therefore, this study’s primary purpose was to explore the mechanisms by which LINC01123 regulates the immune escape and progression of HNSCC. This study confirmed that LINC01123 is competitively bound to miR-214-3p, and miR-214-3p specifically targets B7–H3. The effects of LINC01123, B7–H3, and miR-214-3p on tumor progression, CD8⁺T-cell-mediated immune response, and the tumorigenicity of HNSCC in vitro and in vivo were examined through the downregulation or upregulation of LINC01123, B7–H3, and miR-214-3p. Our results indicated that LINC01123 and B7–H3 were highly expressed in HNSCC and are associated with poor prognosis in patients. Notably, overexpression of LINC01123 or B7–H3 or downregulation of miR-214-3p inhibited the function of CD8⁺T cells and promoted the progression of HNSCC. Therefore, LINC01123 acts as a miR-214-3p sponge to inhibit the activation of CD8⁺T cells and promote the progression of HNSCC by upregulating B7–H3.

Noncoding RNAs link metabolic reprogramming to immune microenvironment in cancers

Altered metabolic patterns in tumor cells not only meet their own growth requirements but also shape an immunosuppressive microenvironment through multiple mechanisms. Noncoding RNAs constitute approximately 60% of the transcriptional output of human cells and have been shown to regulate numerous cellular processes under developmental and pathological conditions. Given their extensive action mechanisms based on motif recognition patterns, noncoding RNAs may serve as hinges bridging metabolic activity and immune responses. Indeed, recent studies have shown that microRNAs, long noncoding RNAs and circRNAs are widely involved in tumor metabolic rewiring, immune cell infiltration and function. Hence, we summarized existing knowledge of the role of noncoding RNAs in the remodeling of tumor metabolism and the immune microenvironment, and notably, we established the TIMELnc manual, which is a free and public manual for researchers to identify pivotal lncRNAs that are simultaneously correlated with tumor metabolism and immune cell infiltration based on a bioinformatic approach.

Expression profile of long non-coding RNAs in porcine lymphnode response to porcine circovirus type 2 infection

Porcine circovirus type 2 (PCV2) can cause various clinical diseases in pigs, resulting in huge losses for the pig farms all over the world. In order to develop a new strategy to control PCV2, it is essential to understand its mechanisms firstly, especially PCV2 interferes with the host's innate immunity. In the present study, lncRNA and mRNA expression profiles in porcine lymphnode response to PCV2 infection were deeply sequenced and analyzed. 3271 novel lncRNAs were identified in all. 1898 mRNAs and 282 lncRNAs showed differential expression between control and PCV2-infected groups. The bioinformatics analysis including lncRNA-mRNA co-expression network construction, as well as GO and KEGG pathway analysis focused on the DEGs was carried out. The results indicated that lncRNAs might participate in PCV2 infection-induced the pathogenesis of immunosuppression through regulating the host's immune responses, biological regulation, response to stimulus, cellular component organization or biogenesis and metabolism. And these differentially expressed lncRNAs might play important roles in response to PCV2 infection in the host's innate immune system. These findings provided a large-scale survey of dysregulated lncRNAs after PCV2 infection, especially the lncRNAs responded to host's innate immune within the lymphnode. This study will provide a novel insight into the lncRNAs' functions and the possible immunosuppressive mechanism induced by PCV2 infection. However, further research will be required to verify the characteristic function of the dysregulated lncRNAs.

LINC00942 Promotes Tumor Proliferation and Metastasis in Lung Adenocarcinoma via FZD1 Upregulation

BACKGROUND

Long non-coding RNAs (lncRNAs) have been reported to play important roles in the progression of human cancers. Herein, bioinformatic analysis identified that LINC00942 was a highly overexpressed lncRNA in lung adenocarcinoma (LUAD). The present study aimed to explore the roles and possible molecular mechanisms of LINC00942 in LUAD.

METHODS

First, on the basis of TCGA database, the expression and prognosis of LINC00942 were analyzed in LUAD tissues. Then, si-LINC00942 was transfected into A549 and H1299 cells to knockdown the expression of LINC00942. Cell viability was detected by MTT assay. Flow cytometry was used to analyze cell apoptosis. The expressions of PCNA, Bax, Bcl-2, and wnt/β-catenin pathway proteins were detected by western blotting. Dual-luciferase reporter assay was used to evaluate the regulatory relationship between LINC00942 and miR-5006-5p, or miR-5006-5p and FZD1.

RESULTS

We discovered that LINC00942 was up-regulated in LUAD tissues compared with adjacent tissues. Besides, we found the increased LINC00942 expression was associated with poor survival. In addition, silencing of LINC00942 suppressed the proliferation, migration, invasion and facilitated the apoptosis of A549 and H1299 cells. Moreover, silencing of LINC00942 repressed the expression of PCNA, Bcl-2, and enhanced Bax expression in A549 and H1299 cells. Mechanically, LINC00942 exerted its effects via enhancing Wnt signaling. LINC00942 functioned as competing endogenous RNA (ceRNA) by binding to miR-5006-5p, upregulating the expression of FZD1, which was a direct target of miR-5006-5p.

CONCLUSION

Our findings indicated that LINC00942/miR-5006-5p/FZD1 axis played important roles in LUAD growth through enhancing Wnt signaling. LINC00942/miR-5006-5p/FZD1 axis might serve as a potential biomarker and therapeutic target for LUAD treatment.

Long non-coding RNAs mediate cerebral vascular pathologies after CNS injuries

Central nervous system (CNS) injuries are one of the leading causes of morbidity and mortality worldwide, accompanied with high medical costs and a decreased quality of life. Brain vascular disorders are involved in the pathological processes of CNS injuries and might play key roles for their recovery and prognosis. Recently, increasing evidence has shown that long non-coding RNAs (lncRNAs), which comprise a very heterogeneous group of non-protein-coding RNAs greater than 200 nucleotides, have emerged as functional mediators in the regulation of vascular homeostasis under pathophysiological conditions. Remarkably, lncRNAs can regulate gene transcription and translation, thus interfering with gene expression and signaling pathways by different mechanisms. Hence, a deeper insight into the function and regulatory mechanisms of lncRNAs following CNS injury, especially cerebrovascular-related lncRNAs, could help in establishing potential therapeutic strategies to improve or inhibit neurological disorders. In this review, we highlight recent advancements in understanding of the role of lncRNAs and their application in mediating cerebrovascular pathologies after CNS injury.

Silenced lncRNA DDX11-AS1 or up-regulated microRNA-34a-3p inhibits malignant phenotypes of hepatocellular carcinoma cells via suppression of TRAF5

Background

Studies have discussed long noncoding RNA DDX11-AS1 (DDX11-AS1)-mediated downstream mechanism in hepatocellular carcinoma (HCC). The goal of this study was to investigate the regulatory mechanism of DDX11-AS1-mediated microRNA-34a-3p (miR-34a-3p)/tumor necrosis factor receptor-associated factor 5 (TRAF5) axis on HCC cells.

Methods

DDX11-AS1, miR-34a-3p and TRAF5 expression levels in HCC were detected. The correlation of DDX11-AS1, miR-34a-3p and TRAF5 in HCC patients was analyzed by Pearson test. HCC cells were transfected with corresponding plasmid/oligonucleotide, and cell proliferation, migration, invasion, apoptosis and tumor formation ability were detected. Bioinformatics software, dual luciferase report experiment and RNA-pull down experiment analysis were applied to verify the targeting relationship between DDX11-AS1, miR-34a-3p and TRAF5.

Results

Elevated DDX11-AS1 and TRAF5 and reduced miR-34a-3p exhibited in HCC. Silenced DDX11-AS1 or up-regulated miR-34a-3p inhibited the proliferation, migration, invasion, promoted apoptosis of HCC cells and repressed the tumor growth in nude mice. In addition, DDX11-AS1 bound to miR-34a-3p to target TRAF5. Silencing TRAF5 or elevating miR-34a-3p expression mitigated up-regulated DDX11-AS1-mediated promotion of tumor growth.

Conclusion

Silenced DDX11-AS1 or up-regulated miR-34a-3p inhibits HCC cell growth via elevation of TRAF5, which could be of great benefit to find early diagnostic markers for HCC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01847-6.

Identification of an immune-related long noncoding RNA signature that predicts prognosis in breast cancer patients

Aim: The purpose of this study was to identify an immune-related long noncoding RNA (lncRNA) signature that predicts the prognosis of breast cancer. Materials & methods: The expression profiles of breast cancer were downloaded from The Cancer Genome Atlas. Cox regression analysis was used to identify an immune-related lncRNA signature. Results: The five immune-related lncRNAs could be used to construct a breast cancer survival prognosis model. The receiver operating characteristic curve evaluation found that the accuracy of the model for predicting the 1-, 3- and 5-year prognosis of breast cancer was 0.688, 0.708 and 0.686. Conclusion: This signature may have an important clinical significance for improving predictive results and guiding the treatment of breast cancer patients.

Role of Long Non-Coding RNAs in the Chemoresistance of Gastric Cancer: A Systematic Review

PURPOSE

Accumulating evidence demonstrates that long non-coding RNAs (lncRNAs) play a vital role in the chemoresistance of gastric cancer (GC). The present systematic review summarises the emerging role, potential targets or pathways and regulatory mechanisms of lncRNAs involved in chemoresistance and proposes a number of clinical implications of lncRNAs as novel therapeutic targets for GC.

METHODS

Studies on lncRNAs involved in the chemoresistance of GC published until July 2020 in the PubMed and Web of Science databases were systematically reviewed and the expression form, role in chemoresistance, targets or pathways, corresponding drugs and potential mechanisms of relevant lncRNAs were summarised in detail.

RESULTS

A total of 48 studies were included in this systematic review. Amongst these studies, 32 involved single drug resistance and 16 involved in multidrug resistance (MDR). The 48 studies collected described 38 lncRNAs in the drug-resistant cells of GC, including 33 upregulated and 5 downregulated lncRNAs. Cisplatin (DDP) was the most studied drug and lncRNA MALAT1 was the most studied lncRNA related to the chemoresistance of GC. The potential mechanisms of chemoresistance for lncRNAs in GC mainly included, amongst others, reduction of apoptosis, induction of autophagy, repair of DNA damage, promotion of epithelial-mesenchymal transition (EMT) and regulation of the related signalling pathways.

CONCLUSION

LncRNAs play a vital role in the chemoresistance of GC and are novel therapeutic targets for the disease. Detailed chemoresistance mechanisms, translational studies and clinical trials on lncRNAs in GC are urgently needed.

The Role of Long Non-Coding RNAs in Intracranial Aneurysms and Subarachnoid Hemorrhage

Intracranial aneurysms (IAs) represent the most complex and relevant problem of modern neurology and neurosurgery. They serve as one of the main causes of non-traumatic subarachnoid hemorrhage (SAH), causing up to 85% of all cases of intracranial hemorrhage, which is associated with frequent disability and high mortality among patients. Unfortunately, the molecular mechanisms of the development and rupture of IAs are still under study. Long non-coding RNAs (lncRNAs) are non-coding RNAs that typically have a length of more than 200 nucleotides. It is known that lncRNAs regulate many processes, such as transcription, translation, cell differentiation, regulation of gene expression, and regulation of the cell cycle. In recent years, a lot of evidence has established their role in human diseases from oncology to cardiovascular disease. Recent studies have shown that lncRNAs may be involved in the pathogenesis of IAs. The study of lncRNAs and its targets in various pathological conditions of a person is a rapidly developing field, and it is likely that the knowledge obtained from these studies regarding the pathogenesis of intracranial aneurysms will have the potential to use lncRNAs in therapy, as well as in the diagnosis and prediction of high aneurysms risk of rupture.

lncRNA expression profiles related to apoptosis and autophagy in peripheral blood mononuclear cells of patients with rheumatoid arthritis

Long noncoding RNAs (lncRNAs) are >200-bp molecules that do not generally code for proteins. Human lncRNAs have well-characterized roles in gene expression regulation, particularly with regard to protein-coding genes, and their dysregulation has been linked to disease. Here, we set out to investigate changes in the expression of lncRNAs related to apoptosis and autophagy in the peripheral blood mononuclear cells (PBMCs) of rheumatoid arthritis (RA). In addition, we aimed to correlate lncRNA expression profiles with clinical indexes and self-perception of patients (SPP). To this end, we employed RNA sequencing of lncRNAs in PBMCs from three patients with RA and three healthy controls. We used bioinformatics to screen several dysregulated lncRNAs related to apoptosis and autophagy. To validate key lncRNA candidates, we performed quantitative reverse transcriptase-PCR on 20 patients with RA and 20 healthy controls. We found the expression of seven lncRNAs (MAPKAPK5-AS1, ENST00000619282, C5orf17, LINC01189, LINC01006, DSCR9 and MIR22HG) was significantly altered in PBMCs of patients with RA. Receiver operating characteristic curve analysis suggested that MIR22HG [area under the curve (AUC) = 0.846, P = 0.000], DSCR9 (AUC = 0.783, P = 0.005), LINC01189 (AUC = 0.677, P = 0.034), MAPKAPK5-AS1 (AUC = 0.644, P = 0.025) and ENST00000619282 (AUC = 0.636, P = 0.043) are potential biomarkers of RA. Spearman's correlation analysis revealed selected lncRNAs correlated with clinical indexes and SPP. Therefore, we highlight that some lncRNAs related to apoptosis and autophagy may serve as potential biomarkers for diagnosis and monitoring of RA progression, which also correlate with several clinical indexes and SPP.

Identification of Differentially Expressed lncRNAs and mRNAs in Children with Acquired Aplastic Anemia by RNA Sequencing

Background

The effects of long noncoding RNAs (lncRNAs) and their related messenger RNAs (mRNAs) remain unknown in children with acquired aplastic anemia (AA). The aim of this study is to screen key lncRNAs and mRNAs and investigate their potential roles in the pathology of acquired AA in children.

Methods

RNA sequencing was performed to identify differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) between blood samples of acquired AA children and healthy controls. cis-regulation, trans-regulation, competing endogenous (Ce) regulation networks of DElncRNAs and DEmRNAs were constructed. A literature search was performed to identify immune- or hematopoietic-related DElncRNA-DEmRNA pairs, and qPCR was conducted to validate the expression of the immune- or hematopoietic-related DElncRNA and DEmRNA.

Results

60 DElncRNAs and 364 DEmRNAs were identified. 13 DElncRNAs were predicted to have 15 cis-regulated target DEmRNAs, 16 DElncRNAs might have 28 trans-regulated DEmRNAs, and 2 DElncRNAs might have 9 Ce-regulated DEmRNAs. After literature screen and qPCR validation, 6 immune- or hematopoietic-related DElncRNA-DEmRNA pairs in the networks above were identified as key RNAs in the pathology of acquired AA.

Conclusion

This study revealed key lncRNAs in children with acquired AA and proposed their potential functions by predicting their target mRNAs, which lay the foundation for future study of potential effects of lncRNAs in children with acquired AA.