Yi Xu1,2, Xingming Jiang1, Yunfu Cui1. 1. Department of Hepatopancreatobiliary Surgery, Second Affiliated Hospital of Harbin Medical University. 2. The Key Laboratory of Myocardial Ischemia, Harbin Medical University, Ministry of Education, Harbin, People's Republic of China.
Abstract
Cholangiocarcinoma (CCA) is one of the most malignant human cancers with increasing incidence worldwide. LncRNAs have emerged as gene regulators and prognostic biomarkers in a variety of neoplasms. PANDAR, a novel cancer-related lncRNA, has been reported to be upregulated in diverse human carcinomas. In this study, we aimed to investigate the clinical significance of lncRNA PANDAR in CCA and explore its functional roles in CCA cells including cell proliferation, apoptosis, migration, invasion and epithelial-to-mesenchymal transition (EMT). The results showed that PANDAR was significantly upregulated in CCA tissue specimens and cell lines, and its high expression was closely associated with lymph node invasion (P=0.004), TNM stage (P=0.034) and postoperative relapse (P=0.006) in patients with CCA. Thus, overexpression of PANDAR could serve as an independent prognostic biomarker of CCA. Furthermore, silencing of PANDAR followed by siRNA significantly inhibited cell proliferation and increased apoptosis in CCA cells. In addition, suppression of PANDAR impaired migration and invasion capacity in vitro partly by affecting EMT. Overall, our findings showed that lncRNA PANDAR serves as a novel prognostic biomarker and therapeutic target for CCA.
Cholangiocarcinoma (CCA) is one of the most malignant humancancers with increasing incidence worldwide. LncRNAs have emerged as gene regulators and prognostic biomarkers in a variety of neoplasms. PANDAR, a novel cancer-related lncRNA, has been reported to be upregulated in diverse humancarcinomas. In this study, we aimed to investigate the clinical significance of lncRNA PANDAR in CCA and explore its functional roles in CCA cells including cell proliferation, apoptosis, migration, invasion and epithelial-to-mesenchymal transition (EMT). The results showed that PANDAR was significantly upregulated in CCA tissue specimens and cell lines, and its high expression was closely associated with lymph node invasion (P=0.004), TNM stage (P=0.034) and postoperative relapse (P=0.006) in patients with CCA. Thus, overexpression of PANDAR could serve as an independent prognostic biomarker of CCA. Furthermore, silencing of PANDAR followed by siRNA significantly inhibited cell proliferation and increased apoptosis in CCA cells. In addition, suppression of PANDAR impaired migration and invasion capacity in vitro partly by affecting EMT. Overall, our findings showed that lncRNA PANDAR serves as a novel prognostic biomarker and therapeutic target for CCA.
Cholangiocarcinoma (CCA) originates from the epithelium of hepatic biliary trees with limited therapeutic options and poor prognosis.1 According to epidemiological data, there is an increasing incidence of CCA worldwide over the past decades, and currently, radical resection is the only treatment option available.2,3 However, due to the early invasive and metastatic characteristics of CCA, most patients are diagnosed at the advanced stages and are ineligible for surgical intervention. In addition, surgical resection is associated with a high rate of recurrence, and the 5-year survival rate is less than 10%. Moreover, CCA is also insensitive to palliative treatments including chemotherapy and radiotherapy.4,5 Therefore, exploring the underlying mechanisms of tumorigenesis and cancer progression, in addition to identifying new therapeutic targets and novel biomarkers associated with diagnosis/prognosis, is very important.With the progression of whole-genome sequencing technology, it has been found that less than 2% of the mammalian genome is in protein-encoded regions and the rest is in noncoding RNAs.6 LncRNAs, imperative members of noncoding RNA family, are transcripts of at least 200 nucleotides in length and have no or limited protein-coding potential.7 Recently, many lncRNAs have been reported to play imperative roles in tumorigenesis, cancer progression and drug resistance.8–11 LncRNAs exert a large number of molecular functions, such as modulation of alternative splicing, chromatin remodeling and RNA metabolism.12–14 They function mostly by binding with miRNAs as sponges or interacting with proteins. However, studies exploring the functional role of lncRNAs in CCA are very scarce.15PANDAR, a novel lncRNA consisting of 1506 nucleotides in length and mapped to 6p21.2, was originally discovered in a p53-dependent manner and interacts with the transcription factor NF-YA. NF-YA is linked to the expression of proapoptotic genes, which has been proved by chromatin immunoprecipitation method.16 Recently, emerging evidence indicates that PANDAR is upregulated and involved in carcinogenesis and progression of various neoplasms including gastric cancer, hepatocellular carcinoma, bladder cancer and colorectal cancer.17–20 However, the biological role and clinical relevance of PANDAR in CCA remains to be elucidated.In this study, we aimed to characterize the expression of PANDAR in CCA and their corresponding adjacent noncancerous tissues and evaluate the association of PANDAR with clinicopathological features of CCA. Furthermore, we explored the functional roles of PANDAR in CCA through examining proliferation, apoptosis, migration, invasion and epithelial-to-mesenchymal transition (EMT) after silencing PANDAR.
Materials and methods
Ethics statement
This study was performed with the approval of the Ethics Review Committee of Harbin Medical University (KY2016-015, approval date, 2016.3.1), and all the recruited patients signed informed consent before participating in this study.
Patients and tissue samples
In this study, we collected 67 pairs of CCA tumor and matched nontumor tissue samples at the Second Affiliated Hospital of Harbin Medical University from 2011 to 2013. Tissue specimens after surgery were collected and immediately stored in liquid nitrogen until required. All recruited patients did not receive chemotherapy or radiotherapy preoperatively, and the diagnoses were confirmed by histological examination.
RNA extraction and quantitative real-time PCR (qRT-PCR)
Total RNA of the tissues and cells was extracted by Trizol reagent (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s recommendation and then converted to cDNA using a reverse transcription kit (Roche, Germany). We performed qRT-PCR by using FastStart Universal SYBR Green Master (Roche Applied Science, Mannheim, Germany) and a BIO-RAD C1000 Thermal Cycler. The following primer sets for PANDAR and GAPDH were used: PANDAR primers forward: 5′-CTGTTAAGGTGGTGGCATTG-3′, reverse: 5′-GGAGGCTCATACTGGCTGAT-3′; GAPDH primers forward: 5′-GGGAGCCAAAAGGGTCAT-3′, reverse: 5′-GAGTCCTTCCACGATACCAA-3′. The relative expression was calculated by delta-delta CT method.
Cell culture and siRNA transfection
HCCC-9810 and RBE cells used in the study were purchased from the Cell Bank of Type Culture of Chinese Academy of Sciences (Shanghai, People’s Republic of China). The other CCA cell lines including QBC939, Huh-28, HuCCT1, KMBC and CCLP-1 and humanintrahepatic biliary epithelial cell (HIBEC) were preserved in our laboratory. The cells were cultivated in DMEM or RPMI-1640 supplemented with 10% FBS (Thermo Fisher Scientific) in a 37°C humidified incubator with 5% CO2. All of the cells were passaged for no more than 6 months. The siRNA specifically targeting PANDAR (si-PANDAR) and a scrambled nucleotide (si-NC) were commercially constructed by GenePharma Corporation (Shanghai, People’s Republic of China). The cells cultured in plates were allowed to grow to half confluence prior to transfection with corresponding siRNA using Lipofectamine 3000 transfection reagent (Thermo Fisher Scientific) in serum-free medium. Flow cytometry was used to detect the fluorescently labeled siRNAs and examine the transfection efficiency. After 48 h, cell lysates were harvested for qRT-PCR and western blot analysis. The target sequences of si-PANDAR are as follows: si-PANDAR-1: 5′-AAUGUGUGCACGUAACAGAUU-3′, si-PANDAR-2: 5′-GGGCAUGUUUUCACAGAGGUU-3′, si-PANDAR-3: 5′-GCAAUCUACAACCUGUCUU-3′.
Cell proliferation assays
Cell proliferation was determined using Cell Counting Kit-8 (CCK-8) (Dojindo, Japan) according to the manufacturer’s instruction. Briefly, QBC939 and RBE cells transfected with siRNA targeting PANDAR or negative control were seeded in 96-well flat-bottomed plates at a density of 1,500 cells/well. At each monitored time (0, 24, 48, 72 and 96 h), 10 μL CCK-8 reagent was added into the corresponding wells and incubated at 37°C for 2 h prior to reading absorbance at 450 nm by using a microplate reader (Tecan, Männedorf, Switzerland). For the clonogenic assay, a certain number of cells transfected with si-PANDAR or si-NC were trypsinized into a single-cell suspension and planted in each well of a 3.5 cm dish and maintained in proper media supplemented with 10% FBS for about 14 days until there were visible colonies formed. After fixation with paraformaldehyde and staining with 0.1% crystal violet (Beyotime, Haimen, People’s Republic of China), number of colonies were counted from three randomly chosen fields under a light microscope.
Cell apoptosis assays
QBC939 and RBE cells transfected with corresponding siRNA in the exponential growth phase were trypsinized and washed twice with cold PBS, followed by resuspension in the provided binding buffer. After being labeled with 5 μL Annexin V and 5 μL propidium iodide (BD Biosciences, San Jose, CA, USA) and incubated at room temperature for 15 min according to the manufacturer’s instructions, the treated cells were then detected by flow cytometry (FACScan, BD Biosciences). Cell apoptosis was also determined by acridine orange/ethidium bromide (AO/EB) double fluorescence staining analysis. AO instantly permeates intact cell membrane and stains the nucleus with green fluorescence. EB stains DNA in orange or red fluorescence when the cell membrane is broken. Thus, AO/EB staining helps to differentiate live cells from apoptotic cells. In brief, after 48 h transfection, cells were stained with prepared AO/EB mixing solution for 5 min (Solarbio, Beijing, People’s Republic of China) before observation under a fluorescence microscope.
Relative caspase activity determination
Caspase-3 and caspase-9 were introduced to evaluate si-PANDAR induced apoptosis in this study. The relative expression of caspase-3 and caspase-9 was measured using the Caspase-3 Activity Kit and Caspase-9 Activity Kit (Solarbio). In brief, proteins were lysed from siRNA-treated cells after 48 h and then were added into 96-well plates with 90 μL of the reaction buffer and caspase substrate (2 mM). After incubation at 37°C for 4 h, caspase activity was measured by a microplate reader at 405 nm.
Wound healing assay
QBC939 and RBE cells transfected with si-PANDAR and their corresponding control cells were plated in 6-well plates and were grown up to a confluency of 70%–80%. Straight wounds were created by scraping the cells with a 200 μL pipette tip and the cell debris was cleaned up. Photographs of wound closure area were captured 0 h and 36 h after wound formation.
Migration and invasion assays
The cell motility assay was performed by using 24-well transwell chamber (Costar, Washington, DC, USA) with or without 40 μL Matrigel (BD Biosciences) coating. Transfected cells (5×104 cells/well) were suspended in 0.2 mL serum-free medium prior to the addition of upper chamber of inserts in the 24-well plate. The lower compartment was filled with 0.6 mL of complete medium. After 24 h incubation, the cells migrated or invaded through membrane were fixed with paraformaldehyde and were stained with 0.1% crystal violet. Three randomly high-power fields were counted. Experiments were performed in triplicate.
Western blot assay
After 48 h transfection, QBC939 and RBE cells were lysed in RIPA buffer with protease inhibitors (Beyotime). Protein concentration was measured by the bicinchoninic acid method (Beyotime). Equal amounts of proteins were electrophoresed on sodium dodecyl sulfatepolyacrylamide gel before being transferred to polyvinylidene fluoride membranes (GE Healthcare, Milwaukee, WI, USA). Afterwards, 5% nonfat milk diluted in Tris-buffered saline containing 0.05% Tween-20 was used to block the membrane. Then, the membranes were incubated with primary antibodies to E-cadherin, N-cadherin, Vimentin, Bax, Bcl-2, MMP-9 or GAPDH (Abcam, Cambridge, MA, USA) overnight, followed by probing with an HRP-conjugated secondary antibody (Cell Signalling Technologies, Danvers, MA, USA) for 2 h. The bound antibodies were detected using an enhanced chemiluminescence kit (Beyotime).
Statistical analysis
All statistical tests were conducted by using GraphPad Prism 5.01 software (GraphPad Software, Inc., La Jolla, CA, USA) and SPSS 19.0 statistical software package (IBM, Armonk, NY, USA). All experimental results from three independent experiments were expressed as mean ± standard deviation (SD) and analyzed using Student’s t-test. Survival analysis was calculated by Kaplan–Meier method with the log-rank test. Fisher’s exact test was applied to evaluate the links between PANDAR and clinicopathological parameters. Survival data were analyzed using univariate and multivariate Cox regression methods. A P-value <0.05 was considered statistically significant.
Results
PANDAR is upregulated in CCA clinical samples and cell lines
First, we explored the relative expression of PANDAR in tumor and adjacent noncancerous tissues of patients with CCA. As shown in Figure 1A, the relative expression of PANDAR in tumor tissues from patients with CCA was remarkably increased compared with the corresponding nontumor tissue by qRT-PCR. The expression of PANDAR in CCA cell lines including KMBC, HCCC-9810, HuCCT1, CCLP-1, Huh-28, QBC939 and RBE was also higher than that of HIBEC (Figure 1B).
Figure 1
Relative expression of PANDAR in CCA and its clinical significance.
Notes: (A) Relative expression of PANDAR in CCA tissues and adjacent nontumor tissues were determined by qRT-PCR; (B) Relative expression of PANDAR in HIBEC and CCA cell lines were determined by qRT-PCR; (C) Kaplan–Meier survival curves showed that overexpression of PANDAR decreased overall survival in patients with CCA. *P<0.05, **P<0.01, ***P<0.001.
PANDAR overexpression is correlated with poor prognosis of patients with CCA
To verify the clinical significance of PANDAR in CCA, we evaluated the correlation of PANDAR expression with the clinicopathological characteristics of patients with CCA. As shown in Table 1, the high PANDAR group patients showed a more frequent lymph node invasion (P=0.004), advanced TNM stage (P=0.034) and more postoperative recurrence (P=0.006). However, the other parameters, such as patients’ gender, age and tumor site, were not associated with the expression of PANDAR. Moreover, Kaplan–Meier overall survival analysis and the log-rank test showed that patients with CCA with higher PANDAR expression had the worse overall survival after surgical procedure (Figure 1C). The univariate Cox regression analyses of overall survival showed that the lymph node invasion (P=0.023), TNM stage (P=0.020), postoperative recurrence (P=0.004) and PANDAR expression (P=0.005) were good prognostic indicators. Furthermore, multivariate Cox regression analyses confirmed the expression of PANDAR as an independent predictor for overall survival in patients with CCA (P=0.041, Table 2).
Table 1
Association between PANDAR expression and clinicopathological characteristics of CCA
To shed light on the roles of PANDAR in CCA, first, three siRNA targeting PANDAR and a scrambled nucleotide were transfected into QBC939 and RBE cells separately because the two cell lines expressed highest levels of PANDAR. Then, transfection analysis by flow cytometry demonstrated that the transfection efficiency of the fluorescence-labeled siRNA and Lipofectamine 3000 transfection reagent mixture reached 70%–90% in the two selected cell lines (Figure 2A). At 48 h after incubation, the relative expression of PANDAR was detected by qRT-PCR, and the results revealed a dramatic decline of PANDAR in all the three siRNAs targeting PANDAR compared to si-NC groups (Figure 2B). Among these, si-PANDAR-2 got the most powerful knockdown efficiency, and thus was chosen for subsequent studies. Compared to the negative control, transfection with PANDAR siRNA resulted in a remarkably restrained QBC939 and RBE cell viability as monitored by CCK-8 proliferation assay (Figure 2C). Moreover, clonogenic assay was conducted to probe the effect of PANDAR on anchor-independent growth of QBC939 and RBE cells. As shown in Figure 2D, CCA cells transfected with si-PANDAR-2 significantly abrogated the formative colony numbers.
Figure 2
Transfection and knockdown efficiency of PANDAR-specific siRNAs and PANDAR depletion retarded cell proliferation and colony formation in CCA cells.
Notes: (A) Transfection efficiency detected by flow cytometry in QBC939 and RBE cells; (B) Relative expression of PANDAR in QBC939 and RBE cells was significantly decreased by three siRNAs specifically targeting PANDAR compared with si-NC; (C) CCK-8 assay showed that silenced PANDAR inhibited cell proliferation of QBC939 and RBE cells; (D) Clonogenic assay showed that silenced PANDAR dramatically impaired the colony formation ability of QBC939 and RBE cells. *P<0.05, **P<0.01.
Abbreviation: CCA, cholangiocarcinoma.
Knockdown of PANDAR increases apoptosis of CCA cells
Flow cytometric analysis clarified that silencing of PANDAR induced apoptosis of CCA cells (Figure 3A). AO/EB double fluorescence staining assay showed similar results subsequently (Figure 3B). Next, we examined the activity of caspase-3 and caspase-9 in response to PANDAR. According to the results, PANDAR knockdown increased the expression of caspase-3 and caspase-9 in CCA cells (Figure 3C). Western blot analysis showed that silencing of PANDAR restrained the expression of Bcl-2 and increased the expression of Bax (Figure 3D).
Figure 3
PANDAR depletion increased apoptosis in CCA cells.
Notes: (A) Flow cytometry for apoptosis analysis showed that silenced PANDAR significantly promoted apoptosis of QBC939 and RBE cells compared to negative control siRNA; (B) AO/EB double fluorescence staining assay showed that silenced PANDAR significantly promoted apoptosis of QBC939 and RBE cells compared to negative control siRNA (×200); (C) Relative expression of caspase-3 and caspase-9 in QBC939 and RBE cells were activated by si-PANDAR-2; (D) Western blot analysis showed that knockdown of PANDAR restrained the expression of Bcl-2 and increased Bax expression in QBC939 and RBE cells. *P<0.05, **P<0.01.
Silenced PANDAR inhibits metastasis in vitro and affects EMT in CCA cells
To estimate the extent of PANDAR function in migration and invasion of CCA, we introduced wound healing assay and Transwell assay for detection. Interestingly, cells transfected with selected siRNA specifically targeting PANDAR led to less wound closure area than si-NC group (Figure 4A). Moreover, Transwell assay demonstrated that the migration capacity as well as invasive potential of the two selected cells was strikingly impaired due to PANDAR depletion (Figure 4B and C). Furthermore, MMP-9 was downregulated after silencing of PANDAR as seen by western blot analysis (Figure 4D). EMT is a vital mechanism associated with cell migration and invasion. Hence, we detected the expression of EMT-related proteins in CCA cells after transfection with corresponding siRNA. According to our results of western blot assay, the expression levels of N-cadherin and Vimentin examined were strikingly decreased, whereas the expression of E-cadherin was boosted followed by si-PANDAR (Figure 4E).
Figure 4
Knockdown of PANDAR impaired migration and invasion capacity in vitro and reversed EMT of CCA cells.
Notes: (A) Wound healing assay showed that silenced PANDAR inhibited migration ability of QBC939 and RBE cells (×50); (B) Transwell migration assay showed that silenced PANDAR inhibited migration capacity of QBC939 and RBE cells (×100); (C) Transwell invasion assay showed that silenced PANDAR inhibited invasion potential of QBC939 and RBE cells (×100); (D) Western blot analysis showed that knockdown of PANDAR restrained the expression of MMP-9 in QBC939 and RBE cells; (E) Knockdown of PANDAR reversed EMT in QBC939 and RBE cells. *P<0.05, **P<0.01.
During the past decades, the worldwide incidence of CCA has increased, especially in Asian countries.2,21,22 In addition, CCA is widely known to be an aggressive malignancy of the biliary duct that shows high mortality rate typically due to its late presentation and limited treatment measures.23,24 Therefore, the detection of specific therapeutic targets as well as identification of new diagnostic and prognostic biomarkers is the need of the hour for patients with CCA. lncRNAs were mistaken for nonfunctional secondary product of transcription for ages. During the past few years, lncRNAs were shown to be vital elements in cell biology and human diseases, especially in various cancers. Numerous evidences indicate that lncRNAs have imperative roles in various cancer-related processes including tumorigenesis, metastasis and angiogenesis.25–27 However, studies regarding the association between lncRNAs and CCA are scarce. Therefore, in this study, we focused on the novel cancer-related lncRNA PANDAR, which has got much attention in recent years and may provide new insights into CCA.In this study, we found that PANDAR expression was upregulated not only in CCA tissue samples but also in CCA cell lines, which is in line with all of the digestive system neoplasms studied to date.17,18,20 We further explored the correlation of expression of PANDAR with clinicopathological features of patients with CCA and the data showed that higher PANDAR expression was closely associated with lymph node invasion, advanced TNM stage and postoperative recurrence. The patients with higher PANDAR expression had statistically shorter survival time than those with lower PANDAR expression (by Kaplan–Meier analysis). Further multivariate analysis confirmed that PANDAR serves as an independent indicator for patients’ overall survival. Taken together, the results of this study suggests that PANDAR is an oncogenic player or a predictor of overall survival for CCA.Recent studies indicated that PANDAR was upregulated and led to cell growth in multiple cancers, such as gastric cancer, hepatocellular carcinoma, bladder carcinoma, colon cancer and breast cancer.17–20,28 A recent study indicated that PANDAR acts as a tumor-promoter in breast cancer by regulating G1/S transition. PANDAR-mediated promotion of cell proliferation occurs via regulation of p16INK4A.28 In order to understand the biological functions of PANDAR in CCA, we examined the differences in cell proliferation after suppressing PANDAR expression by using CCK-8 and colony formation assays. The results showed that cell viability and colony formation capacity were both abrogated in the PANDAR knockdown groups. Apoptosis signaling can be activated by receptor-mediated pathway and mitochondrial pathway and both of the pathways are correlated with the regulation of caspase-3, caspase-9, Bcl-2 and Bax. In our study, flow cytometric analysis as well as AO/EB staining assay showed that PANDAR acts as a defender against cancer apoptosis, and apoptosis-related factors (caspase-3 and caspase-9) were both activated analyzed by caspases detection. Furthermore, Bcl-2 was downregulated after PANDAR knockdown, whereas Bax was activated followed by silencing of PANDAR. These results suggest that Bcl-2/caspase-3 pathway act by mediating apoptosis induced by the downregulation of PANDAR.Next, we performed Transwell and western blot assays to explore whether expression of PANDAR affects migration capacity and invasive potential of CCA cells in vitro and investigate its associations with MMP-9. Transwell assay showed that siRNA targeting PANDAR significantly reduced the migration and invasion of CCA cells. MMP-9 has been proved to facilitate the invasion capacity of cancer cells by damaging the histological barrier around the extracellular matrix.29 In this study, we found that MMP-9 was decreased followed by the knockdown of PANDAR. These results demonstrate that PANDAR facilitates invasion potential partly by upregulating MMP-9 expression. EMT, a major mechanism involved in cell invasiveness, was proven to be involved in tumor metastasis, and the characteristics of EMT included the lost E-cadherin and enhanced N-cadherin and Vimentin expression.30–32 It has been reported that LINC01133 inhibits EMT and metastasis in colorectal carcinoma via binding to SRSF6 as a target mimic.33 In addition, overexpressed MALAT1 in tongue cancer activated EMT and restrained cell apoptosis by regulating Wnt/beta-catenin signaling pathway.34 For CCA, E-cadherin and N-cadherin switch promotes cell progression through TGF-beta-induced EMT, which suggested a target for illuminating the mechanisms of invasion and metastasis in CCA.35 In this study, western blot results showed reversed EMT markers after knockdown of PANDAR, which indicates that PANDAR is involved in mechanisms related to the facilitation of CCA cell migration and invasion, and impaired the potential of metastasis in vitro partly by regulating the EMT pathway.In our study, there are still some limitations to come to an incontrovertible conclusion. CCA is a multifactorial, multistep, aggressive disease with various genes activated or silenced. Moreover, the molecular mechanisms of lncRNAs seem to be highly diverse and act at miscellaneous levels to manipulate cancer development and progression. Hence, further studies are needed to unravel the precise mechanism behind the altered expression of PANDAR in CCA.
Conclusion
We first confirmed that PANDAR was statistically overexpressed in CCA and therefore might be a potential indicator of a negative prognosis for patients with CCA. Moreover, downregulation of PANDAR inhibited the proliferation, migration, invasion and EMT of CCA cells and facilitate its apoptosis. In summary, these results indicate that PANDAR plays oncogenic roles in CCA and serves as a potential prognostic biomarker and therapeutic target for treatment of CCA.
Authors: Michelle L DeOliveira; Steven C Cunningham; John L Cameron; Farin Kamangar; Jordan M Winter; Keith D Lillemoe; Michael A Choti; Charles J Yeo; Richard D Schulick Journal: Ann Surg Date: 2007-05 Impact factor: 12.969
Authors: P Carninci; T Kasukawa; S Katayama; J Gough; M C Frith; N Maeda; R Oyama; T Ravasi; B Lenhard; C Wells; R Kodzius; K Shimokawa; V B Bajic; S E Brenner; S Batalov; A R R Forrest; M Zavolan; M J Davis; L G Wilming; V Aidinis; J E Allen; A Ambesi-Impiombato; R Apweiler; R N Aturaliya; T L Bailey; M Bansal; L Baxter; K W Beisel; T Bersano; H Bono; A M Chalk; K P Chiu; V Choudhary; A Christoffels; D R Clutterbuck; M L Crowe; E Dalla; B P Dalrymple; B de Bono; G Della Gatta; D di Bernardo; T Down; P Engstrom; M Fagiolini; G Faulkner; C F Fletcher; T Fukushima; M Furuno; S Futaki; M Gariboldi; P Georgii-Hemming; T R Gingeras; T Gojobori; R E Green; S Gustincich; M Harbers; Y Hayashi; T K Hensch; N Hirokawa; D Hill; L Huminiecki; M Iacono; K Ikeo; A Iwama; T Ishikawa; M Jakt; A Kanapin; M Katoh; Y Kawasawa; J Kelso; H Kitamura; H Kitano; G Kollias; S P T Krishnan; A Kruger; S K Kummerfeld; I V Kurochkin; L F Lareau; D Lazarevic; L Lipovich; J Liu; S Liuni; S McWilliam; M Madan Babu; M Madera; L Marchionni; H Matsuda; S Matsuzawa; H Miki; F Mignone; S Miyake; K Morris; S Mottagui-Tabar; N Mulder; N Nakano; H Nakauchi; P Ng; R Nilsson; S Nishiguchi; S Nishikawa; F Nori; O Ohara; Y Okazaki; V Orlando; K C Pang; W J Pavan; G Pavesi; G Pesole; N Petrovsky; S Piazza; J Reed; J F Reid; B Z Ring; M Ringwald; B Rost; Y Ruan; S L Salzberg; A Sandelin; C Schneider; C Schönbach; K Sekiguchi; C A M Semple; S Seno; L Sessa; Y Sheng; Y Shibata; H Shimada; K Shimada; D Silva; B Sinclair; S Sperling; E Stupka; K Sugiura; R Sultana; Y Takenaka; K Taki; K Tammoja; S L Tan; S Tang; M S Taylor; J Tegner; S A Teichmann; H R Ueda; E van Nimwegen; R Verardo; C L Wei; K Yagi; H Yamanishi; E Zabarovsky; S Zhu; A Zimmer; W Hide; C Bult; S M Grimmond; R D Teasdale; E T Liu; V Brusic; J Quackenbush; C Wahlestedt; J S Mattick; D A Hume; C Kai; D Sasaki; Y Tomaru; S Fukuda; M Kanamori-Katayama; M Suzuki; J Aoki; T Arakawa; J Iida; K Imamura; M Itoh; T Kato; H Kawaji; N Kawagashira; T Kawashima; M Kojima; S Kondo; H Konno; K Nakano; N Ninomiya; T Nishio; M Okada; C Plessy; K Shibata; T Shiraki; S Suzuki; M Tagami; K Waki; A Watahiki; Y Okamura-Oho; H Suzuki; J Kawai; Y Hayashizaki Journal: Science Date: 2005-09-02 Impact factor: 47.728
Authors: K Araki; T Shimura; H Suzuki; S Tsutsumi; W Wada; T Yajima; T Kobayahi; N Kubo; H Kuwano Journal: Br J Cancer Date: 2011-11-08 Impact factor: 7.640
Authors: B Luzón-Toro; R M Fernández; J M Martos-Martínez; M Rubio-Manzanares-Dorado; G Antiñolo; S Borrego Journal: Sci Rep Date: 2019-10-07 Impact factor: 4.379
Authors: Dino Bekric; Daniel Neureiter; Markus Ritter; Martin Jakab; Martin Gaisberger; Martin Pichler; Tobias Kiesslich; Christian Mayr Journal: J Clin Med Date: 2020-04-22 Impact factor: 4.241