Research Paper Volume 10, Issue 4 pp 775—788

Identifying circRNA-associated-ceRNA networks in the hippocampus of Aβ1-42-induced Alzheimer's disease-like rats using microarray analysis

Zhe Wang1, , Panpan Xu1, , Biyue Chen1, , Zheyu Zhang2, , Chunhu Zhang2, , Qiong Zhan3, , Siqi Huang2, , Zi-an Xia2, , Weijun Peng1, ,

  • 1 Department of Integrated Traditional Chinese and Western Medicine, The Second Xiangya Hospital, Central South University, Changsha, Hunan 410011, China
  • 2 Department of Integrated Traditional Chinese and Western Medicine, Xiangya Hospital, Central South University, Changsha, Hunan 410008, China
  • 3 Department of Neurology, The Second Xiangya Hospital, Central South University, Changsha, Hunan 410011, China

Received: March 17, 2018       Accepted: April 20, 2018       Published: April 27, 2018      

https://doi.org/10.18632/aging.101427
How to Cite

Copyright: © 2018 Wang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Alzheimer’s disease (AD) is the most common form of dementia worldwide. Accumulating evidence indicates that non-coding RNAs are strongly implicated in AD-associated pathophysiology. However, the role of these ncRNAs remains largely unknown. In the present study, we used microarray analysis technology to characterize the expression patterns of circular RNAs (circRNAs), microRNAs (miRNAs), and mRNAs in hippocampal tissue from Aβ1-42-induced AD model rats, to integrate interaction data and thus provide novel insights into the mechanisms underlying AD. A total of 555 circRNAs, 183 miRNAs and 319 mRNAs were identified to be significantly dysregulated (fold-change ≥ 2.0 and p-value < 0.05) in the hippocampus of AD rats. Quantitative real-time polymerase chain reaction (qRT-PCR) was then used to validate the expression of randomly-selected circRNAs, miRNAs and mRNAs. Next, GO and KEGG pathway analyses were performed to further investigate ncRNAs biological functions and potential mechanisms. In addition, we constructed circRNA-miRNA and competitive endogenous RNA (ceRNA) regulatory networks to determine functional interactions between ncRNAs and mRNAs. Our results suggest the involvement of different ncRNA expression patterns in the pathogenesis of AD. Our findings provide a novel perspective for further research into AD pathogenesis and might facilitate the development of novel therapeutics targeting ncRNAs.

Introduction

Alzheimer’s disease (AD), the most common cause of dementia worldwide, is becoming more prevalent due to the aging population, and represents one of the grand challenges to health care systems [1]. Although substantial progress has been made in the identification of disease-related molecular and cellular processes over the last decade, the molecular mechanisms that underlie the pathogenesis of AD remain largely unknown [2], and none of the pharmacological treatments presently available for AD are able to stop or slow down the progression of this disease [3].Therefore, further investigation of the underlying disease mechanisms are urgently required in order to better understand AD and to facilitate the development of effective therapeutic strategies.

Currently, accumulating evidence indicates that noncoding RNAs (ncRNAs), especially microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), are strongly implicated in AD-associated pathophysiology, including amyloid-β plaques and neurofibrillary tangles, synaptic loss and neuronal death [4,5]. In a previous study, we investigated the hippocampal expression patterns of dysregulated lncRNAs in a rat model of AD using microarray analysis and demonstrated that lncRNAs contributed to the pathogenesis of AD [6]. To further investigate the regulatory role of ncRNAs in AD, we focused upon circRNAs, a class of ncRNAs that are highly expressed in the mammalian brain [7,8], and can regulate transcriptional and post-transcriptional gene expression [9]. Unlike linear RNA, circRNA is formed with covalently closed continuous loops without 5’-3’ polarity and a poly(A) tail, and might function as microRNA sponges to modulate the expression of parental genes through the competing endogenous RNA (ceRNA) network [10]. Recent studies have provided evidence that the circRNA-associated ceRNA network may play a crucial role in many disease processes, including AD [11]. For example, Lukiw et al. [12,13] demonstrated that circRNA-7(ciRS-7) acted as a natural miRNA sponge for miRNA-7 and regulated the expression of ubiquitin-conjugating enzyme E2A (UBE2A) and the epidermal growth factor receptor (EGFR) in sporadic AD hippocampal brain. Zhang et al. [14] characterized circRNA-associated-ceRNA networks in the cerebral cortex of senescence-accelerated mouse prone 8(SAMP8). Additionally, Zhong et al. [15] described the expression of circRNAs in different ages of SAMP8 mice. However, the potential role of circRNAs in the pathogenesis of AD is still in its infancy and has yet to be characterized, particularly the role of circRNA-associated-ceRNA networks in the hippocampus of AD.

In the present study, we used a microarray analysis approach to identify differentially-expressed profiles of circRNAs, miRNAs and mRNAs in the hippocampus of Aβ1-42-induced AD model rats, which is a useful experimental animal model of AD which emphasizes the inflammatory component of the disease pathology, and strongly complements the use of transgenic animal models in advancing our understanding of AD [16]. Next we performed Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. We also constructed a circRNA-associated-ceRNA network. Findings arising from this study will expand our understanding of the potential role of the circRNA-associated-ceRNA network involved in AD pathogenesis and therapeutic strategy.

Results

Altered circRNA, miRNA and mRNA profiles in the AD rat hippocampus

The expression of circRNAs, miRNAs, and mRNAs in the hippocampus of AD rats were profiled using microarrays. Using a threshold of fold change (FC) ≥ 2.0 and a p-value < 0.05, 555 circRNAs, 183 miRNAs, and 319 mRNAs were significantly altered in the AD rat hippocampus compared to sham controls. The results showed that 444 circRNAs, 93 miRNAs, and 173 mRNAs were up-regulated, while 111 circRNAs, 90 miRNAs, and 146 mRNAs were down-regulated. The most up-regulated were circRNA_008964, miR-133a-5p, and BGLAP with FCs of 16.97, 355.67, and 22.96, respectively, whereas circRNA_017759, miR-551b-3p, and IL-1RN were the most down-regulated, with corresponding FCs of 6.32, 257.63 and 14.21. The top 20 up-regulated and 20 down-regulated circRNAs, miRNAs, and mRNAs in the AD group are listed in Tables 1-3. Hierarchical clustering and volcano plot visualization (Fig. 1 and 2) suggested that the expression level of circRNAs, miRNAs and mRNAs were distinguishable and variable.

Table 1. Top 40 differently expressed circRNAs in microarray analysis.

circRNAp-valueFDRFold changeRegulationchromstrandcircRNA_typeGene symbol
rno_circRNA_0089640.0006060.09436516.9714upchr2+exonicPpp3ca
mmu_circRNA_320070.0021110.10806916.57521upchr18-sense overlappingRGD1308601
rno_circRNA_0045600.0186090.18284210.47134upchr14-exonicRGD1305110
rno_circRNA_0172950.0376750.22527110.14399upchr9-exonicClk1
rno_circRNA_0163550.0292810.20713410.00513upchr8-sense overlappingBirc2
rno_circRNA_0012850.011660.165069.170892upchr1+exonicPdcd11
rno_circRNA_0111190.0438190.2364278.989032upchr3-exonicExt2
mmu_circRNA_335920.0232370.1964178.459upchr3-exonicAcvr1
rno_circRNA_0037240.0445370.2375037.782604upchr12-exonicRGD1309762
rno_circRNA_0106840.0164630.178357.75411upchr3-exonicGolga1
rno_circRNA_0110540.0343310.2162667.619158upchr3+intergenic
rno_circRNA_0075990.0238290.1973547.284635upchr18+sense overlappingZfp516
rno_circRNA_0071080.0281130.2051017.274158upchr18+exonicWdr33
mmu_circRNA_248280.0191810.1846867.260869upchr6+sense overlappingTssc1
rno_circRNA_0055470.0154650.1755147.112599upchr15+exonicMipep
rno_circRNA_0009670.0147180.1733477.044118upchr1+exonicCpsf7
rno_circRNA_0150730.0179180.1828427.039097upchr7-exonicPolr3b
rno_circRNA_0034630.0134530.1701966.770216upchr12+exonicSlc46a3
rno_circRNA_0020710.0084680.1490266.626886upchr10+exonicTekt5
hsa_circRNA_1020430.005760.1381316.473515upchr10+exonicAcaca
rno_circRNA_0177590.0028990.1118286.320172downchrX+intergenic
mmu_circRNA_368140.010920.1618143.519062downchr5+exonicRad23b
rno_circRNA_0078790.0170540.1796283.370089downchr19-exonicNr3c2
rno_circRNA_0012350.0426610.2349353.322749downchr1-exonicGot1
rno_circRNA_0055600.009720.1543092.990798downchr15+exonicDpysl2
rno_circRNA_0139400.000040.0665852.94782downchr6+exonicPtprn2
rno_circRNA_0139410.0002490.0776012.893453downchr6+exonicPtprn2
rno_circRNA_0086230.0420570.2337772.8228downchr2+sense overlappingAsh1l
rno_circRNA_0139810.0295540.2073812.812952downchr6-exonicNrxn1
rno_circRNA_0130260.0444910.237472.785261downchr5+exonicRere
rno_circRNA_0139420.0001820.0776012.733282downchr6+exonicPtprn2
rno_circRNA_0130250.0419070.2337772.704128downchr5+sense overlappingRere
mmu_circRNA_317940.0018290.1054662.669903downchr18+exonicCamk2a
rno_circRNA_0156120.0295450.2073812.653884downchr7-exonicTrps1
mmu_circRNA_368130.0269850.2036132.590337downchr5+exonicRad23b
rno_circRNA_0130220.0411580.232682.585987downchr5+exonicRere
rno_circRNA_0139430.0088760.152.525526downchr6+exonicPtprn2
rno_circRNA_0111800.0035190.1196022.509546downchr4-exonicMagi2
rno_circRNA_0117310.0121780.1673162.50498downchr4-exonicGrin2b
rno_circRNA_0141660.0170490.1796282.501832downchr6+exonicDtnb

Table 2. Top 40 differently expressed miRNAs in microarray analysis.

NameFCP-valueFDRRegulationSham -1Sham -2Sham-3AD-1AD-2AD-3
rno-miR-133a-5p355.67781.49E-060.000597up0.0300160.0242630.0259.6768199.5321748.988871
rno-miR-133a-3p180.92550.0009180.023498up0.1390210.296360.27812552.4619340.351336.27822
rno-miR-1b121.90970.0179590.081167up0.355450.1109190.20781337.4737730.8156513.89984
rno-miR-133b-3p112.21120.0003520.016336up0.2559240.5892550.512559.6395945.4417447.2655
rno-miR-152-5p77.690730.0315240.108617up0.0063190.0103990.0251.6548221.0695650.516693
rno-miR-378a-5p57.079560.000380.01635up0.0537120.0849220.0406252.9001693.3843483.947536
rno-miR-208b-3p50.910590.0175940.080291up0.0205370.0259970.0140631.1590521.3965220.529412
rno-miR-144-3p45.249515.72E-086.88E-05up0.1042650.2201040.21258.1641298.1686967.960254
rno-miR-48637.367580.0003670.016336up0.14060.2166380.2531257.6785118.7095656.419714
rno-miR-499-5p33.396830.0038520.038921up0.3491310.2218370.193759.7715749.9965225.771065
rno-miR-10b-5p33.323010.0009770.023498up0.0315960.0259970.0359381.2199661.0747830.82194
rno-miR-155-5p33.288752.36E-070.000189up0.0173780.0155980.0109380.495770.4747830.491256
rno-miR-3561-3p30.310980.0097890.060993up0.0221170.0086660.0140630.5668360.5269570.265501
rno-miR-3573-5p27.408150.0039710.039624up0.0205370.0034660.0109380.3587140.2173910.381558
rno-miR-378a-3p23.93810.0112890.064185up0.932072.3916811.07656327.2267327.8921750.21622
rno-miR-223-3p18.232470.0135430.071138up0.2654030.2530330.3140633.9086297.3286963.941176
rno-miR-3557-3p16.942050.0490430.137085up0.0142180.0034660.006250.2182740.1252170.062003
rno-miR-199a-5p16.613418.89E-060.002256up0.3064770.1767760.2078133.5888323.9286963.963434
rno-miR-451-5p15.920290.000890.023498up1.3917852.1577121.81406333.786828.1860923.41653
rno-miR-380-5p15.600350.0040150.039631up0.0473930.0693240.1015630.8917091.4817391.031797
rno-miR-551b-3p0.0038810.0015510.027567down5.9194318.8232249.1968750.0338410.040.019078
rno-miR-153-3p0.0047664.37E-086.88E-05down7.6635077.9081467.7593750.0253810.0365220.049285
rno-miR-539-5p0.0051030.0006880.020933down1.5071091.6499131.3750.008460.006957N/A
rno-miR-376a-3p0.0054570.0001880.012127down3.928915.072794.806250.0186130.0313040.025437
rno-miR-124-3p0.0074932.61E-050.003452down92.0932180.7625680.406250.5922170.5582610.747218
rno-miR-136-3p0.0079820.0002260.013604down2.0205372.5649912.6031250.0253810.020870.011129
rno-miR-129-5p0.0082430.0419110.124593down2.5197476.675912.7906250.0473770.0434780.007949
rno-miR-9a-3p0.0098720.0003540.016336down41.5007956.2651648.034380.4483930.5426090.448331
rno-miR-410-3p0.0100020.0026310.032583down1.9194311.4436741.7828130.016920.017391N/A
rno-miR-153-5p0.0117050.0032470.035819down0.5071090.5823220.68750.0016920.012174N/A
rno-miR-409a-5p0.0121510.003480.036711down0.360190.2980940.4093750.0016920.006957N/A
rno-miR-138-5p0.0151860.0140030.072249down8.2922593.9358754.7578130.0812180.0956520.081081
rno-miR-369-3p0.0153710.0040610.039771down1.2590841.8561532.2921880.0253810.0243480.033386
rno-miR-4950.0158329.38E-060.002256down1.2654031.3778161.4234380.0236890.0295650.011129
rno-miR-129-2-3p0.0203830.0050480.043957down3.3617691.8596192.3468750.0524530.0747830.027027
rno-miR-376b-3p0.0234080.0011050.025237down1.9620851.9532062.7296880.0795260.0521740.023847
rno-miR-127-5p0.0281050.009530.06067down1.2322272.1074521.131250.0609140.0504350.014308
rno-miR-132-5p0.028830.0355420.115982down1.142184.0883883.131250.0896790.0782610.073132
rno-miR-496-3p0.0293830.0190140.082393down0.502371.3154250.9484380.0253810.040.015898
rno-miR-35420.0297910.0297990.10517down0.0758290.1525130.1171880.0033840.003478N/A

Table 3. Top 40 differently expressed mRNAs in microarray analysis.

Genbank AccessionGene SymbolP-valueFDRFold ChangeRegulationAD-1AD-2AD-3Sham -1Sham -2Sham -3
NM_013414Bglap0.000550.02640622.96703up8.0520339.0514688.5751234.180264.551223.382667
NM_001270665Tnnt30.0281490.12876312.35968up5.8061895.2759918.6429433.2113193.2090042.422092
NM_012587Ibsp0.0010960.0314245.771831up6.2528165.2791575.6475493.3847293.163913.043795
NM_013122Igfbp20.0389930.153085.754793up11.066089.30309712.162068.3008288.4275178.228595
NM_012605Mylpf0.0267810.1252975.641142up7.4663027.1291579.4621925.4646035.4796665.625421
NM_057104Enpp20.0384890.151985.07149up12.4611210.6819713.244729.6466469.64582810.0681
NM_001000438Olr12140.0022170.0393444.895564up5.1856194.6429265.0700822.3608333.2412772.422092
NM_031511Igf20.0490390.1754234.757985up13.6567911.6663614.3154110.7884310.9109811.1881
XR_086340LOC3014440.0146150.0904224.340022up10.7889310.891279.7461958.4489618.9326087.691718
XM_001059752Vom2r453.47E-050.0164.296108up4.3667944.3800624.6724772.3608332.3273172.422092
NM_001108140Cd3e0.0035430.0473364.22607up4.7389034.6570024.9623162.3608332.4021693.357269
NM_001107541Art10.0430710.1625623.875367up4.7572094.7242556.2576833.7994083.6546482.422092
NM_031703Aqp30.0085240.0685043.685033up3.968044.7698044.8469252.3608332.3905173.188386
BC078859Mpz0.0289460.1303313.67334up7.9995987.8975566.2860765.4219195.4266665.703367
NM_019212Acta10.0475680.172133.632821up6.6376626.6670138.6222185.4968245.3446615.502137
NM_030838Slco1a50.0317830.1371863.400608up8.4330277.2967419.1473056.4442436.4373626.69809
NM_001107564Ano10.0012390.032453.397198up8.8169459.0671168.6959156.8000287.0445667.442346
XM_227107RGD15618410.0143160.0893963.391897up10.0780310.320659.0682678.3931948.0271467.760331
NM_001109599Pou2af10.0081570.0671193.171449up3.9347134.7885683.7049842.3608332.6499132.422092
NM_001004129Stfa2l10.0116190.079943.137655up7.9001528.1287247.2930165.5701286.4591836.343522
NM_022194Il1rn0.0367530.14845114.21197down5.0168144.9923115.1593869.62517210.565316.465138
XR_361775LOC1009103670.0045830.05281910.16713down3.4687414.3079652.3218736.8759656.7161666.54397
BC098733LOC3627950.0324480.1388889.708935down6.8443716.6756957.02145610.9133311.354128.112008
NM_017210Dio30.0026690.0427826.519498down6.2318446.6210635.5032289.0682679.0627218.339429
XM_006255043LOC6894530.0498740.1772245.455306down5.0155995.8415143.1285027.6754957.217996.435111
XM_006227100LOC5014670.0302320.1337515.431563down5.40696.0831513.7945818.0356317.57486.998302
NM_199253Pcsk90.0021330.0390565.204765down2.4265242.6658532.9462935.6046215.0127884.56076
NM_130748Slc38a40.0009380.0302924.934402down4.5874824.6379614.4756717.3545936.7838756.471272
NM_031512Il1b0.0343290.1433024.617557down4.0374014.3647453.8747767.5025996.2632315.132481
NM_001014221LOC3633370.0394230.1540474.587823down8.4737129.351027.10960510.8374410.834519.855817
XR_349427LOC1025572060.0431260.1627114.261659down7.1662777.9467255.9045769.5120479.4774238.302353
XM_006224554Tmco5a0.0145830.0903734.222635down5.5903174.6651176.3819157.7476077.6741067.450068
XM_006255374LOC6851830.0114310.0793583.885825down4.7192625.0915993.7787116.8746886.4233156.166232
NM_001014221LOC3633370.0429260.1623223.775706down9.38010610.536038.45947411.8147611.3977610.91333
NM_001127377LOC6806630.0013380.0331243.748348down3.411883.8284323.5327975.7664425.6291985.096233
NM_031561Cd360.0213480.1109653.577859down4.2270854.4626584.1763515.9027047.0716495.40903
XM_006222689LOC6917120.0171040.0983593.482986down4.7236325.2009543.8460466.7337756.4814395.956392
NM_053587S100a90.020890.1098583.409311down5.4651034.9430093.9675656.1991666.6372696.847683
NM_001014091Ccdc330.0009360.0302923.293413down2.4775672.6023662.9261514.6755794.2536474.235609
NM_053822S100a80.0067070.0624173.266306down5.7682225.8369434.8620417.008027.2132297.368936
Heat map of the top 40 differentially-expressed circRNAs (A), miRNAs (B), and mRNAs (C) in AD hippocampal tissue. The data are depicted as matrices in which each row represents one circRNA, miRNA, or mRNA and each column represents one of the hippocampal samples. Relative circRNA, miRNA, or mRNA expression is depicted according to the color scale shown at the top. Red and green represent high and low relative expression, respectively; -2.0, 0, and 2.0 are fold-changes in the corresponding spectrum. The magnitude of deviation from the median is represented by color saturation.

Figure 1. Heat map of the top 40 differentially-expressed circRNAs (A), miRNAs (B), and mRNAs (C) in AD hippocampal tissue. The data are depicted as matrices in which each row represents one circRNA, miRNA, or mRNA and each column represents one of the hippocampal samples. Relative circRNA, miRNA, or mRNA expression is depicted according to the color scale shown at the top. Red and green represent high and low relative expression, respectively; -2.0, 0, and 2.0 are fold-changes in the corresponding spectrum. The magnitude of deviation from the median is represented by color saturation.

Volcano plot of differentially-expressed circRNAs (A), miRNAs (B), and mRNAs (C) between AD and sham hippocampal samples. Volcano plots were constructed using fold-change values and p-values. The vertical lines correspond to 2.0-fold up- and down-regulation between normal and AD samples (N vs.D), and the horizontal lines represent p-values. Red plot points represent differentially-expressed circRNAs with statistical significance.

Figure 2. Volcano plot of differentially-expressed circRNAs (A), miRNAs (B), and mRNAs (C) between AD and sham hippocampal samples. Volcano plots were constructed using fold-change values and p-values. The vertical lines correspond to 2.0-fold up- and down-regulation between normal and AD samples (N vs.D), and the horizontal lines represent p-values. Red plot points represent differentially-expressed circRNAs with statistical significance.

Expression profile validation

To validate the accuracy and reliability of the microarray profiling data, some transcripts, including four circRNAs (rno_circRNA_001555, rno_circRNA_010684, rno_circRNA_01398, and rno_circRNA_017759), four miRNAs (rno-miR-181a-2-3p, rno-miR-124-3p, rno-miR-136-3p, and rno-miR-206-3p), and four mRNAs (IGF2, IGFBP2, S100a8, and IGF1) were randomly selected for quantitative real-time polymerase chain reaction (qRT-PCR) analysis in nine samples including those used for microarray analysis. As shown in Fig. 3, the microarray data were consistent with the qRT-PCR results in terms of the expression levels of the validated ncRNAs and mRNAs.

qRT-PCR-validated ncRNA and mRNA expression changes. Microarray validation by qRT-PCR. Expression levels of four circRNAs, four miRNAs, and four mRNAs were determined by qRT-PCR. Each assay was performed at least three times. *p

Figure 3. qRT-PCR-validated ncRNA and mRNA expression changes. Microarray validation by qRT-PCR. Expression levels of four circRNAs, four miRNAs, and four mRNAs were determined by qRT-PCR. Each assay was performed at least three times. *p < 0.05.

GO and KEGG pathway analyses of differentially-expressed mRNAs

GO and KEGG pathway analyses allow for the functional annotation of differentially-expressed mRNAs. GO analysis indicated that the most enriched mRNAs correlated with the extracellular region part (GO: 0044421) of the cellular component analysis (Figure 4a). Meanwhile, the majority of genes were related to toll-like receptor binding (GO:0035325) in the molecular functions and regulation of biological quality (GO:0065008) aspect of the biological processes analyses (Figure 4b and 4c). KEGG pathway analysis indicated that 10 KEGG pathways (p < 0.05) were associated with dysregulated mRNAs involved in neuroactive ligand-receptor interaction, AMP-activated protein kinase (AMPK) signaling pathway, longevity regulating pathway-multiple species, fatty acid elongation, inflammatory mediator regulation of transient receptor potential (TRP) channels, p53 signaling pathway, hematopoietic cell lineage, adipocytokine signaling pathway, focal adhesion and osteoclast differentiation (Figure 4d).

Gene ontology (GO) enrichment and pathway analysis of differentially- expressed mRNAs showing the most significantly enriched (-log10, p-value) GO terms of differentially-expressed mRNAs according to (A) cellular components, (B) molecular functions, and (C) biological processes. (D) The top ten enrichment scores (-log10, p-value) of significantly enriched KEGG pathways.

Figure 4. Gene ontology (GO) enrichment and pathway analysis of differentially- expressed mRNAs showing the most significantly enriched (-log10, p-value) GO terms of differentially-expressed mRNAs according to (A) cellular components, (B) molecular functions, and (C) biological processes. (D) The top ten enrichment scores (-log10, p-value) of significantly enriched KEGG pathways.

Prediction of circRNA-miRNA interactions

To determine the function of circRNA, interactions between circRNAs and their target miRNAs were theoretically predicted by conserved seed-matching sequences. To find the potential miRNA target, two confirmed circRNAs (rno_circRNA_008964 and rno_circRNA_017759) were selected, and circRNA-miRNA interaction was predicted using Arraystar's miRNA target prediction software based on the TargetScan and miRanda databases. The potential miRNA targets of rno_circRNA_008964 included rno-miR-216b-5p, rno-miR-181d-5p, rno-miR-337-5p, rno-miR-497-3p, and rno-miR-181b-5p (Fig. 5a). For rno_circRNA_007879, the potential miRNA targets included rno-miR-702-5p, rno-miR-3547, rno-miR-329-5p, rno-miR-203b-5p and rno-miR-3576 (Fig. 5).

Example of detailed annotation for circRNA-miRNA interactions. (A) rno

Figure 5. Example of detailed annotation for circRNA-miRNA interactions. (A) rno_circRNA_008964. (B) rno_circRNA_017759. 8mer: bases from number 2 to 8 matched perfectly, and number 1 base is A; 7mer-m8: bases from number 2 to 8 matched perfectly, and number 1 base is not A; 6mer: bases from number 2 to 7 matched perfectly, and number 1 base is not A; offset 6mer: bases from number 3 to 8 matched perfectly; imperfect match: there is imperfect base match from number 2 to 7; M: circRNA-miRNA interaction can be predicted by miRanda; T: circRNA-miRNA interaction can be predicted by TargetScan.

Construction of a circRNA-miRNA regulatory network

A circRNA-miRNA regulatory network, based on the microarray results, was constructed containing 245 circRNAs, 144 miRNAs and 279 relationships (Fig. 6). We observed that one circRNA could regulate multiple miRNAs in different ways, while one miRNA could be regulated by multiple circRNAs. For example, circRNA_34441 was co-related to four dysregulated miRNAs and MiR-153-5p was co-related with 25 dysregulated circRNAs. Thus, there appears to be a complex circRNA-miRNA regulatory network involved in the pathogenesis of AD.

CircRNA-miRNA network analysis. A total of 245 circRNAs and 144 miRNAs containing 279 relationships were selected to generate a network map. The circRNA-miRNA co-expression network was constructed with Cytoscape V3.5.0 software. Within this network, V-shaped nodes represent circRNA and octagon nodes represent miRNA. Red and green represent up- and down-regulation, respectively. Node size represents p-values (larger nodes for more significant p-values).

Figure 6. CircRNA-miRNA network analysis. A total of 245 circRNAs and 144 miRNAs containing 279 relationships were selected to generate a network map. The circRNA-miRNA co-expression network was constructed with Cytoscape V3.5.0 software. Within this network, V-shaped nodes represent circRNA and octagon nodes represent miRNA. Red and green represent up- and down-regulation, respectively. Node size represents p-values (larger nodes for more significant p-values).

Construction of a circRNA-miRNA-mRNA regulatory network

To explore the molecular mechanism of ncRNAs, a circRNA-miRNA-mRNA regulatory network was constructed with 140 circRNAs as decoys, 140 miRNAs as centres, and 20 mRNAs as targets based on the microarray data (Fig. 7). Because there are binding sites between circRNAs and miRNAs, circRNAs could indirectly regulate miRNA target genes by competitively binding to miRNA as a miRNA sponge. For instance, we focused on miR-7a-5p, circRNA_101834 and circRNA_004690 that could regulate AQP3 expression by competing miRNA response elements (MREs) of miR-7a-5p. These data suggested that circRNAs harbor MREs and play pivotal regulatory roles in AD.

CircRNA-miRNA-mRNA regulatory network. The regulatory network consists of 140 circRNAs, 140 miRNAs, and 20 genes with 503 relationships. V-shaped, octagon, and diamond nodes represent circRNA, miRNA, and mRNA, respectively. Red and green represent up- and down-regulation, respectively. Node size represents p-values (larger nodes for more significant p-values).

Figure 7. CircRNA-miRNA-mRNA regulatory network. The regulatory network consists of 140 circRNAs, 140 miRNAs, and 20 genes with 503 relationships. V-shaped, octagon, and diamond nodes represent circRNA, miRNA, and mRNA, respectively. Red and green represent up- and down-regulation, respectively. Node size represents p-values (larger nodes for more significant p-values).

Discussion

To the best of our knowledge, this is the first integrated microarray analysis of circRNA, miRNA and mRNA expression profiles in the hippocampus of Aβ1-42-induced AD model rats. With FC ≥ 2.0 and p-value < 0.05 thresholds, 444 up-regulated and 111 down-regulated circRNAs, 93 up-regulated and 90 down-regulated miRNAs, and 173 up-regulated and 146 down-regulated mRNAs showed significant differential expression between the AD and control groups. These transcripts are thought to be associated with the pathogenesis of AD. For instance, the S100A8 precedes Aβ plaque formation [17], IGFBP-2 drives AD neurodegeneration [18], miR-146a-5p facilitates neuroinflammation in AD pathogenesis [19], and miR-132-3p contributes to tau hyper-phosphorylation [20]. Our qRT-PCR validation showd that the qRT-PCR results and microarray data were consistent, indicating that the latter were reliable.

To better understand the biological functions and potential mechanisms of mRNAs in the pathogenesis of AD, we performed GO and KEGG pathway analysis. Among the GO terms found in this study, the extracellular region part (GO:0044421), toll-like receptor binding (GO:0035325), and the regulation of biological quality (GO:0065008) have been reported to play important roles in AD development. Remarkable among the KEGG pathways were the AMPK and p53 signaling pathways, both known to potentially mediate AD pathogenesis [2123].

Increasing lines of evidence indicate that circRNAs can function as endogenous sponges to influence miRNA activity, thereby regulating other transcripts containing miRNA-binding sites [24,25]. Therefore, we examined circRNA-miRNA interactions and found that the majority of circRNAs contained one or more miRNA binding sites based on sequence analyses. The association of miRNAs with AD indicated that circRNAs might have a regulatory role in AD. For example, rno_circRNA_003295 is potentially able to interact with rno-miR-206-3p, rno_circRNA_002441 with miR-144-3p, and rno_circRNA_012846 with rno-miR-10a-5p. Moreover, we constructed a circRNA-miRNA regulatory network to investigate potential relationships among the circRNAs and miRNAs. Our results indicated that circRNA-miRNA regulatory networks might play important roles in the pathogenesis of AD. However, beyond acting as a miRNA sponge, circRNAs might also influence mRNA transcriptional levels by interacting with the Pol II complex in the nucleus [26], or by directly binding to RNA-binding proteins and RNA-associated proteins to form RNA-protein complexes [27]. Furthermore, circRNAs can act on gene expression trans-functionally by competing with pre-mRNA splicing machinery [28], and endogenous circRNAs can be used as templates to encode proteins [29]. Therefore, the biological function of circRNA in the pathogenesis of AD requires further investigation.

Perturbations in competing endogenous RNA (ceRNA) regulatory networks (ceRNETs) including mRNAs, miRNAs, and circRNAs have been proposed to play critical roles in the pathogenesis of human disease [3032]. For example, the circRNA-7/miRNA-7/UBE2A signaling pathway is significantly dysregulated in AD [33,34]. The application of ceRNET analysis may provide a novel view of interplay between ncRNAs and mRNAs, thereby offering some insight into molecular pathways. Our results provide the first assessment of ceRNETs in AD and suggest that specific ceRNETs are involved in disease etiology and pathogenesis.

There are certain limitations to our study which should be considered when interpreting our findings. Firstly, the sample size was limited, which could have resulted in under- or over-estimation of the numbers of altered circRNAs, miRNAs, and mRNAs. Larger sample sizes are needed to confirm our findings. Secondly, our results refer only to hippocampal samples from Aβ1-42-induced AD rats. Further research with other AD models, using additional sample sources such as blood and cerebrospinal fluid, is now needed to more accurately capture the pathophysiology of AD. Thirdly, our analyses were performed using gene expression microarrays with limited dynamic range that lack the ability to identify novel features. RNA-sequence technology should be utilized to unravel previously inaccessible transcriptome complexities. Finally, because the functions of circRNAs and miRNAs remain largely unknown, the interpretation of our data was not straightforward. Thus, we only predicted the functions of differentially-expressed mRNAs. Future studies that overcome such limitations are now highly warranted.

In conclusion, we identified dysregulated expression profiles of circRNAs, miRNA, and mRNAs in the hippocampus of AD rats, and present an innovative data integration analysis of circRNAs, miRNAs and mRNAs. Our results indicate that ncRNAs may interact to regulate the expression of their target protein-genes involved in the pathogenesis of AD. The next step is to validate and expand these findings in future studies, which might ultimately enable us to fully elucidate the mechanisms underlying AD.

Materials and Methods

Ethics statement

All animal protocols were approved by the Central South University (Changsha, China) Institutional Review Board and were performed in compliance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals. This investigation was conducted in accordance with appropriate ethical standards and the Declaration of Helsinki, as well as national and international guidelines.

Brain samples

All hippocampal tissue samples were from brain tissues obtained in our previous study [6]. Briefly, 20 adult male Sprague Dawley rats (250 ± 30 g) were randomly divided into AD (n = 10) and sham (n = 10) groups. We performed intracerebroventricular injections of Aβ1-42 oligomers into the cerebral ventricles to induce a validated AD model, as previously described [35,36]. Briefly, the animals were anaesthetized and placed in a stereotactic frame, then the Aβ1–42 oligomers were injected bilaterally into the lateral ventricles through a stainless-steel cannula.

RNA extraction

Total RNA was extracted from each hippocampal tissue sample by soaking in TRIzol Reagent (Invitrogen, Carlsbad, CA, USA) in accordance with the manufacturer’s instructions, which mainly includes homogenization, phase separation, RNA precipitation, RNA wash and RNA solubilization. RNA quantity and quality were measured using a NanoDrop ND-1000 (Thermo Scientific, Waltham, MA, USA), and RNA integrity was assessed by standard denaturing agarose gel electrophoresis.

Microarray analysis

Sample labelling and microarray hybridization for gene expression were performed according to the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology, Santa Clara, CA, USA), which included RNA purification, transcription into fluorescent-labeled cRNA, and hybridization onto the Rat 4x44K Gene Expression Array (Agilent). Finally, the hybridized arrays were washed, fixed, and scanned using the Agilent DNA Microarray Scanner G2505C. RNA labeling and array hybridization for miRNA analysis were conducted following the Exiqon manual (Vedbaek, Denmark). After quality control, the miRCURY™ Hy3™/Hy5™ Power labeling kit (Exiqon) was used for miRNA labeling according to the manufacturer guidelines. After stopping the labeling procedure, Hy3™-labeled samples were hybridized on the Rat miRCURY LNA™ microRNA Array 7th Gen (Exiqon), according to the array manual. The slides were scanned using the Axon GenePix 4000B microarray scanner (Axon Instruments, Foster City, CA, USA). circRNA sample preparation and microarray hybridization were performed based on Arraystar standard protocols (Super RNA Labeling Kit; Arraystar, Rockville, MD, USA). Briefly, total RNAs were digested with Rnase R (Epicentre, Inc., Madison, WI, USA) to remove linear RNAs and enrich for circular RNAs. Then, the enriched circular RNAs were amplified and transcribed into fluorescent-labeled cRNA using a random priming method (Super RNA Labeling Kit). The labeled cRNAs were hybridized onto the Arraystar Rat circRNA Array (8x15K, Arraystar). After washing the slides, the arrays were scanned using the Agilent Scanner G2505C. All data collection was performed using Agilent Feature Extraction software (version 11.0.1.1). KangChen Bio-tech (Shanghai, China) performed all microarray analyses.

Quantitative real-time polymerase chain reaction (qRT-PCR) validation

As previously described [31], total RNA was isolated using TRIzol Reagent, and then reverse-transcribed into cDNA using SuperScript III Reverse Transcriptase (Invitrogen) according to the manufacturer’s instruction. An Applied Biosystems ViiA™7 Real-Time PCR System and 2× PCR Master Mix were used to perform qRT-PCR (Arraystar) in accordance with the manufacturer's instructions. The relative circRNA and mRNA expression levels were calculated using the 2-ΔΔCt method and were normalized to GAPDH as an endogenous reference transcript [37]. miRNA expression levels were normalized to that of U6. The specific primers for each gene are listed in Table 4. Data shown represent the means of three experiments.

Table 4. Primers designed for qRT-PCR validation of candidate circRNAs, miRNAs, and mRNAs.

Forward primerReverse PrimerPrProduct length
(bp)
Tm(°C)
rno_circRNA_0015555’- ATGAGCAATGACTCCCCAGAA-3’5’- GAGAGTATGGTCTGTTGCGTTG-3’6060
rno_circRNA_0106845’- TGGATCTAAAGCAGCTACAGAA-3’5’- CTTTGGTTCCATTCATCCTTAT-3’8260
rno_circRNA_0139815’- CTACCTTGAGCTGCACATACTG-3’5’- TTTGTCCACCACCTTTGCT-3’5860
rno_circRNA_0177595’- GAGTATCCACTGGTGACGACTG-3’5’- AATATGCTGATCTTGTTTTCACC-3’6960
rno-miR-181a-2-3p5’-GGACCACTGACCGTTGAC-3’5’-CAGTGCGTGTCGTGGAG-3’6460
rno-miR-206-3p5’- TGGGGTGGAATGTAAGGAAGT-3’5’-CAGTGCGTGTCGTGGAGT-3’6560
rno-miR-124-3p5’-GGGTAAGGCACGCGGT-3’5’-GTGCGTGTCGTGGAGTCG-3’6160
rno-miR-136-3p5’- GGGGACATCATCGTCTCAAAT -3’5’-CAGTGCGTGTCGTGGAGT-3’6560
IGFBP25’- TCTACTCCCTGCATATCCCCA-3’5’- GGTTCACACACCAGCACTCC-3’10560
IGF25’- GCTTGTTGACACGCTTCAGTT-3’5’- TAGACACGTCCCTCTCGGA-3’17960
IGF15’- GGGCATTGTGGATGAGTGTTG -3’5’- GCTGGGACTTCTGAGTCTTGG -3’14860
S100A85’-GGGAATCACCATGCCCTCTAC-3’5’-GCCCACCCTTATCACCAACAC-3’16860

GO annotations and KEGG pathway analyses

GO annotations and KEGG pathway analyses were performed to investigate the roles of all differentially-expressed mRNAs, as previously described [38,39]. Briefly, GO analysis was applied to elucidate genetic regulatory networks of interest by forming hierarchical categories according to the molecular functions, biological processes, and cellular component aspects of the differentially expressed genes (http://www.geneontology.org). The -log10 (p-value) denotes enrichment scores that represent the significance of GO term enrichment among differentially-expressed genes. KEGG pathway analyses were performed to explore significant pathways associated with the differentially-expressed genes (http://www.genome.jp/kegg/). The -log10 (p-value) denotes an enrichment score for the significance of pathway correlations.

Annotation for circRNA-miRNA interaction

As described previously [40,41], circRNA-miRNA interactions were predicted with Arraystar's home-made miRNA target prediction software based on TargetScan (http://www.targetscan.org/) and miRanda (http://microrna.org/). The top five putative target miRNAs were identified. Then, we constructed a circRNA-miRNA regulatory network using the Cytoscape software V3.5.0 (San Diego, CA, USA).

CircRNA-associated ceRNA network construction

The circRNA-associated ceRNA network was constructed and visually displayed using the Cytoscape software V3.5.0 (San Diego, CA, USA) based on microarray data analysis results, as previous described [42]. Different shapes and colors represent different RNA types and regulated relationships, respectively. Node size was inversely proportion to the p-value.

Statistical analysis

All data were analysed using SPSS version 22.0 software (IBM Corp. Armonk, NY, USA) and presented as mean ± standard error of the mean (SEM). Student’s t-tests were used for comparisons between two groups, whereas one-way analysis of variance was performed for repeated measures. False discovery rates were calculated to correct p-values. Differences with p < 0.05 were considered to be statistically significant. Fold changes (FCs) and Student’s t-tests were used to determine the statistical significance of the microarray results. FC ≥ 2 and p < 0.05 were used as thresholds for designating differentially-expressed ncRNAs and mRNAs.

Acknowledgements

We thank Mr. Xing Xiang from KangChen Bio-tech Inc. for technical advice relating to our bioinformatics analysis.

Conflicts of Interest

The authors declare that they have no competing interests.

Funding

This work was financially supported by the National Natural Science Foundation of China (Nos. 81373705, 81603670, 81601140, 81673951) and the Hunan Provincial Natural Science Foundation of China (No. 2017JJ3459).

References

  • 1. Jia J, Wei C, Chen S, Li F, Tang Y, Qin W, Zhao L, Jin H, Xu H, Wang F, Zhou A, Zuo X, Wu L, et al. The cost of Alzheimer’s disease in China and re-estimation of costs worldwide. Alzheimers Dement. 2018; 14:483–91. https://doi.org/10.1016/j.jalz.2017.12.006 [PubMed]
  • 2. Zádori D, Veres G, Szalárdy L, Klivényi P, Vécsei L. Alzheimer’s Disease: Recent Concepts on the Relation of Mitochondrial Disturbances, Excitotoxicity, Neuroinflammation, and Kynurenines. J Alzheimers Dis. 2018; 62:523–47. https://doi.org/10.3233/JAD-170929 [PubMed]
  • 3. Association As. 2017 Alzheimer’s disease facts and figures. Alzheimers Dement. 2017; 13:325–73. https://doi.org/10.1016/j.jalz.2017.02.001
  • 4. Idda ML, Munk R, Abdelmohsen K, Gorospe M. Noncoding RNAs in Alzheimer’s disease. Wiley Interdiscip Rev RNA. 2018; 9:1–13. https://doi.org/10.1002/wrna.1463 [PubMed]
  • 5. Millan MJ. Linking deregulation of non-coding RNA to the core pathophysiology of Alzheimer’s disease: an integrative review. Prog Neurobiol. 2017; 156:1–68. https://doi.org/10.1016/j.pneurobio.2017.03.004 [PubMed]
  • 6. Yang B, Xia ZA, Zhong B, Xiong X, Sheng C, Wang Y, Gong W, Cao Y, Wang Z, Peng W. Distinct Hippocampal Expression Profiles of Long Non-coding RNAs in an Alzheimer’s Disease Model. Mol Neurobiol. 2017; 54:4833–46. https://doi.org/10.1007/s12035-016-0038-5 [PubMed]
  • 7. You X, Vlatkovic I, Babic A, Will T, Epstein I, Tushev G, Akbalik G, Wang M, Glock C, Quedenau C, Wang X, Hou J, Liu H, et al. Neural circular RNAs are derived from synaptic genes and regulated by development and plasticity. Nat Neurosci. 2015; 18:603–10. https://doi.org/10.1038/nn.3975 [PubMed]
  • 8. Rybak-Wolf A, Stottmeister C, Glažar P, Jens M, Pino N, Giusti S, Hanan M, Behm M, Bartok O, Ashwal-Fluss R, Herzog M, Schreyer L, Papavasileiou P, et al. Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol Cell. 2015; 58:870–85. https://doi.org/10.1016/j.molcel.2015.03.027 [PubMed]
  • 9. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013; 495:384–88. https://doi.org/10.1038/nature11993 [PubMed]
  • 10. Du WW, Zhang C, Yang W, Yong T, Awan FM, Yang BB. Identifying and Characterizing circRNA-Protein Interaction. Theranostics. 2017; 7:4183–91. https://doi.org/10.7150/thno.21299 [PubMed]
  • 11. Rong D, Sun H, Li Z, Liu S, Dong C, Fu K, Tang W, Cao H. An emerging function of circRNA-miRNAs-mRNA axis in human diseases. Oncotarget. 2017; 8:73271–81. https://doi.org/10.18632/oncotarget.19154 [PubMed]
  • 12. Lukiw W, Zhao Y, Rogaev E, Bhattacharjee S. A Circular RNA (circRNA) ciRS-7 in Alzheimer’s disease (AD) targets miRNA-7 trafficking and promotes deficits in the expression of the ubiquitin conjugase (UBE2A) and the epidermal growth factor receptor (EGFR). The FASEB journal. 2016; 30:587.1. https://www.fasebj.org/doi/abs/10.1096/fasebj.30.1_supplement.587.1.
  • 13. Lukiw WJ. Circular RNA (circRNA) in Alzheimer’s disease (AD). Front Genet. 2013; 4:307. https://doi.org/10.3389/fgene.2013.00307 [PubMed]
  • 14. Zhang S, Zhu D, Li H, Li H, Feng C, Zhang W. Characterization of circRNA-Associated-ceRNA Networks in a Senescence-Accelerated Mouse Prone 8 Brain. Mol Ther. 2017; 25:2053–61. https://doi.org/10.1016/j.ymthe.2017.06.009 [PubMed]
  • 15. Huang JL, Qin MC, Zhou Y, Xu ZH, Yang SM, Zhang F, Zhong J, Liang MK, Chen B, Zhang WY, Wu DP, Zhong ZG. Comprehensive analysis of differentially expressed profiles of Alzheimer’s disease associated circular RNAs in an Alzheimer’s disease mouse model. Aging (Albany NY). 2018; 10:253–65. https://doi.org/10.18632/aging.101387 [PubMed]
  • 16. McLarnon JG, Ryu JK. Relevance of abeta1-42 intrahippocampal injection as an animal model of inflamed Alzheimer’s disease brain. Curr Alzheimer Res. 2008; 5:475–80. https://doi.org/10.2174/156720508785908874 [PubMed]
  • 17. Lodeiro M, Puerta E, Ismail MA, Rodriguez-Rodriguez P, Rönnbäck A, Codita A, Parrado-Fernandez C, Maioli S, Gil-Bea F, Merino-Serrais P, Cedazo-Minguez A. Aggregation of the Inflammatory S100A8 Precedes Aβ Plaque Formation in Transgenic APP Mice: Positive Feedback for S100A8 and Aβ Productions. J Gerontol A Biol Sci Med Sci. 2017; 72:319–28. https://doi.org/10.1093/gerona/glw073 [PubMed]
  • 18. Lane EM, Hohman TJ, Jefferson AL, and Alzheimer’s Disease Neuroimaging Initiative. Insulin-like growth factor binding protein-2 interactions with Alzheimer’s disease biomarkers. Brain Imaging Behav. 2017; 11:1779–86. https://doi.org/10.1007/s11682-016-9636-0 [PubMed]
  • 19. Zhang B, Wang A, Xia C, Lin Q, Chen C. A single nucleotide polymorphism in primary-microRNA-146a reduces the expression of mature microRNA-146a in patients with Alzheimer’s disease and is associated with the pathogenesis of Alzheimer’s disease. Mol Med Rep. 2015; 12:4037–42. https://doi.org/10.3892/mmr.2015.3968 [PubMed]
  • 20. Lau P, Bossers K, Janky R, Salta E, Frigerio CS, Barbash S, Rothman R, Sierksma AS, Thathiah A, Greenberg D, Papadopoulou AS, Achsel T, Ayoubi T, et al. Alteration of the microRNA network during the progression of Alzheimer’s disease. EMBO Mol Med. 2013; 5:1613–34. https://doi.org/10.1002/emmm.201201974 [PubMed]
  • 21. Cai Z, Yan LJ, Li K, Quazi SH, Zhao B. Roles of AMP-activated protein kinase in Alzheimer’s disease. Neuromolecular Med. 2012; 14:1–14. https://doi.org/10.1007/s12017-012-8173-2 [PubMed]
  • 22. Checler F, Alves da Costa C. p53 in neurodegenerative diseases and brain cancers. Pharmacol Ther. 2014; 142:99–113. https://doi.org/10.1016/j.pharmthera.2013.11.009 [PubMed]
  • 23. Merlo P, Frost B, Peng S, Yang YJ, Park PJ, Feany M. p53 prevents neurodegeneration by regulating synaptic genes. Proc Natl Acad Sci USA. 2014; 111:18055–60. https://doi.org/10.1073/pnas.1419083111 [PubMed]
  • 24. Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014; 505:344–52. https://doi.org/10.1038/nature12986 [PubMed]
  • 25. Xie L, Mao M, Xiong K, Jiang B. Circular RNAs: A Novel Player in Development and Disease of the Central Nervous System. Front Cell Neurosci. 2017; 11:354. https://doi.org/10.3389/fncel.2017.00354 [PubMed]
  • 26. Zhang Y, Zhang XO, Chen T, Xiang JF, Yin QF, Xing YH, Zhu S, Yang L, Chen LL. Circular intronic long noncoding RNAs. Mol Cell. 2013; 51:792–806. https://doi.org/10.1016/j.molcel.2013.08.017 [PubMed]
  • 27. Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016; 44:2846–58. https://doi.org/10.1093/nar/gkw027 [PubMed]
  • 28. Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, Evantal N, Memczak S, Rajewsky N, Kadener S. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 2014; 56:55–66. https://doi.org/10.1016/j.molcel.2014.08.019 [PubMed]
  • 29. Legnini I, Di Timoteo G, Rossi F, Morlando M, Briganti F, Sthandier O, Fatica A, Santini T, Andronache A, Wade M, Laneve P, Rajewsky N, Bozzoni I. Circ-ZNF609 Is a Circular RNA that Can Be Translated and Functions in Myogenesis. Mol Cell. 2017; 66:22–37.e9. https://doi.org/10.1016/j.molcel.2017.02.017 [PubMed]
  • 30. Guo LL, Song CH, Wang P, Dai LP, Zhang JY, Wang KJ. Competing endogenous RNA networks and gastric cancer. World J Gastroenterol. 2015; 21:11680–87. https://doi.org/10.3748/wjg.v21.i41.11680 [PubMed]
  • 31. Jin X, Feng CY, Xiang Z, Chen YP, Li YM. CircRNA expression pattern and circRNA-miRNA-mRNA network in the pathogenesis of nonalcoholic steatohepatitis. Oncotarget. 2016; 7:66455–67. https://doi.org/10.18632/oncotarget.12186 [PubMed]
  • 32. Chen L, Zhang S, Wu J, Cui J, Zhong L, Zeng L, Ge S. circRNA_100290 plays a role in oral cancer by functioning as a sponge of the miR-29 family. Oncogene. 2017; 36:4551–61. https://doi.org/10.1038/onc.2017.89 [PubMed]
  • 33. Zhao Y, Alexandrov PN, Jaber V, Lukiw WJ. Deficiency in the Ubiquitin Conjugating Enzyme UBE2A in Alzheimer’s Disease (AD) is Linked to Deficits in a Natural Circular miRNA-7 Sponge (circRNA; ciRS-7). Genes (Basel). 2016; 7:1–9. https://doi.org/10.3390/genes7120116 [PubMed]
  • 34. Lukiw WJ, Zhao Y, Pogue AI, Bhattacharjee S, Dua P. Circular RNA (circRNA-7; ciRS-7) impacts microrna-7 trafficking and downregulates the ubiquitin-conjugating enzyme E2A (UBE2A) in sporadic Alzheimer’s disease (AD) brain. Alzheimers Dement. 2015; 11:624. https://doi.org/10.1016/j.jalz.2015.06.878
  • 35. Xia Z, Peng W, Cheng S, Zhong B, Sheng C, Zhang C, Gong W, Cheng S, Li J, Wang Z. Naoling decoction restores cognitive function by inhibiting the neuroinflammatory network in a rat model of Alzheimer’s disease. Oncotarget. 2017; 8:42648–63. https://doi.org/10.18632/oncotarget.17337 [PubMed]
  • 36. Sheng C, Xu P, Zhou K, Deng D, Zhang C, Wang Z. Icariin Attenuates Synaptic and Cognitive Deficits in an Aβ1-42-Induced Rat Model of Alzheimer’s Disease. BioMed Res Int. 2017; 2017:7464872. https://doi.org/10.1155/2017/7464872 [PubMed]
  • 37. Gebhardt FM, Scott HA, Dodd PR. Housekeepers for accurate transcript expression analysis in Alzheimer’s disease autopsy brain tissue. Alzheimers Dement. 2010; 6:465–74. https://doi.org/10.1016/j.jalz.2009.11.002 [PubMed]
  • 38. Chen R, Liu L, Xiao M, Wang F, Lin X. Microarray expression profile analysis of long noncoding RNAs in premature brain injury: A novel point of view. Neuroscience. 2016; 319:123–33. https://doi.org/10.1016/j.neuroscience.2016.01.033 [PubMed]
  • 39. Zhao Z, Bai J, Wu A, Wang Y, Zhang J, Wang Z, Li Y, Xu J, Li X. Co-LncRNA: investigating the lncRNA combinatorial effects in GO annotations and KEGG pathways based on human RNA-Seq data. Database (Oxford). 2015; 2015:2015. https://doi.org/10.1093/database/bav082 [PubMed]
  • 40. Lin SP, Ye S, Long Y, Fan Y, Mao HF, Chen MT, Ma QJ. Circular RNA expression alterations are involved in OGD/R-induced neuron injury. Biochem Biophys Res Commun. 2016; 471:52–56. https://doi.org/10.1016/j.bbrc.2016.01.183 [PubMed]
  • 41. Zhou J, Xiong Q, Chen H, Yang C, Fan Y. Identification of the Spinal Expression Profile of Non-coding RNAs Involved in Neuropathic Pain Following Spared Nerve Injury by Sequence Analysis. Front Mol Neurosci. 2017; 10:91. https://doi.org/10.3389/fnmol.2017.00091 [PubMed]
  • 42. Liu W, Zhang J, Zou C, Xie X, Wang Y, Wang B, Zhao Z, Tu J, Wang X, Li H, Shen J, Yin J. Microarray Expression Profile and Functional Analysis of Circular RNAs in Osteosarcoma. Cell Physiol Biochem. 2017; 43:969–85. https://doi.org/10.1159/000481650 [PubMed]