Utility of a PI3K/mTOR Inhibitor (NVP-BEZ235) for Thyroid Cancer Therapy

Background

We assessed the utility of the dual PI3K/mTOR inhibitor NVP-BEZ235 (BEZ235) as single agent therapy and in combination with conventional chemotherapy for thyroid cancer.

Methodology/Principal Findings

Eight cell lines from four types of thyroid cancer (papillary, follicular, anaplastic, medullary) were studied. The cytotoxicity of BEZ235 and five conventional chemotherapeutic agents alone and in combination was measured using LDH assay. Quantitative western blot assessed expression of proteins associated with cell cycle, apoptosis and signaling pathways. Cell cycle distribution and apoptosis were measured by flow cytometry. Murine flank anaplastic thyroid cancers (ATC) were treated with oral BEZ235 daily. We found that BEZ235 effectively inhibited cell proliferation of all cancer lines, with ATC exhibiting the greatest sensitivity. BEZ235 consistently inactivated signaling downstream of mTORC1. BEZ235 generally induced cell cycle arrest at G0/G1 phase, and also caused apoptosis in the most sensitive cell lines. Baseline levels of p-S6 ribosomal protein (Ser235/236) and p27 correlated with BEZ235 sensitivity. Growth of 8505C ATC xenograft tumors was inhibited with BEZ235, without any observed toxicity. Combination therapy of BEZ235 and paclitaxel consistently demonstrated synergistic effects against ATC in vitro.

Conclusions

BEZ235 as a single therapeutic agent inhibits thyroid cancer proliferation and has synergistic effects in combination with paclitaxel in treating ATC. These findings encourage future clinical trials using BEZ235 for patients with this fatal disease.     

Midostaurin preferentially attenuates proliferation of triple-negative breast cancer cell lines through inhibition of Aurora kinase family

Background

Breast cancer is classified into three subtypes by the expression of biomarker receptors such as hormone receptors and human epidermal growth factor receptor 2. Triple-negative breast cancer (TNBC) expresses none of these receptors and has an aggressive phenotype with a poor prognosis, which is insensitive to the drugs that target the hormone receptors and human epidermal growth factor receptor 2. It is, thus, required to develop an effective therapeutic reagent to treat TNBC.

Results

The study using a panel of 19 breast cancer cell lines revealed that midostaurin, a multi-target protein kinase inhibitor, suppresses preferentially the growth of TNBC cells comparing with non-TNBC cells. Clustering analysis of the drug activity data for the panel of cancer cell lines predicted that midostaurin shares the target with Aurora kinase inhibitors. Following studies indicated that midostaurin attenuates the phosphorylation reaction mediated by Aurora kinase in the cells and directly inhibits this protein kinase in vitro, and that this reagent induces apoptosis accompanying accumulation of 4N and 8N DNA cells in TNBC cells.

Conclusion

Midostaurin suppresses the proliferation of TNBC cells among the breast cancer cell lines presumably through the inhibition of the Aurora kinase family. The precise study of midostaurin on cell growth will contribute to the development of the drug for the treatment of TNBC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-015-0150-2) contains supplementary material, which is available to authorized users.     

Suppression of HBV by Tenofovir in HBV/HIV Coinfected Patients: A Systematic Review and Meta-Analysis

Background

Hepatitis B coinfection is common in HIV-positive individuals and as antiretroviral therapy has made death due to AIDS less common, hepatitis has become increasingly important. Several drugs are available to treat hepatitis B. The most potent and the one with the lowest risk of resistance appears to be tenofovir (TDF). However there are several questions that remain unanswered regarding the use of TDF, including the proportion of patients that achieves suppression of HBV viral load and over what time, whether suppression is durable and whether prior treatment with other HBV-active drugs such as lamivudine, compromises the efficacy of TDF due to possible selection of resistant HBV strains.

Methods

A systematic review and meta-analysis following PRISMA guidelines and using multilevel mixed effects logistic regression, stratified by prior and/or concomitant use of lamivudine and/or emtricitabine.

Results

Data was available from 23 studies including 550 HBV/HIV coinfected patients treated with TDF. Follow up was for up to seven years but to ensure sufficient power the data analyses were limited to three years. The overall proportion achieving suppression of HBV replication was 57.4%, 79.0% and 85.6% at one, two and three years, respectively. No effect of prior or concomitant 3TC/FTC was shown. Virological rebound on TDF treatment was rare.

Interpretation

TDF suppresses HBV to undetectable levels in the majority of HBV/HIV coinfected patients with the proportion fully suppressed continuing to increase during continuous treatment. Prior treatment with 3TC/FTC does not compromise efficacy of TDF treatment. The use of combination treatment with 3TC/FTC offers no significant benefit over TDF alone.     

miR-26a and its Target CKS2 Modulate Cell Growth and Tumorigenesis of Papillary Thyroid Carcinoma

Background

While many studies have shown that levels of miR-26a are lower in papillary thyroid carcinoma (PTC), the role and mechanism of miR-26a in PTC are unclear.

Method

We used database searches to select potential mRNA targets of miR-26a. Anti-miR-26a, miR-26a mimic, siRNA for CKS2 and their effects on cell growth, cell-cycle distribution and colony formation were evaluated. We also evaluate the over-expressed miR-26a in TPC-1 cells in severe combined immune-deficient mice. We used luciferase reporter assays, real-time PCR and western blot analysis to measure the expression and activity of miR-26a, CKS2, and related factors such as cyclin B1, cyclin A, cdk1, bcl-xl and Akt. Finally, we measured the relationship between the levels of miR-26a and CKS2 in PTC and normal thyroid tissues.

Results

Relative to normal thyroid tissues, miR-26a is consistently down-regulated in TPC specimens, and CKS2 was identified as a potential target. Up-regulated miR-26a expression or down-regulated CKS2 expression in TPC-1 and CGTH W3 cells lines caused G2 phase-arrest. Decreased miR-26a expression or increased CKS2 expression could have inverse function on PTC cell lines. CyclinB1, cyclinA, bcl-xl and AKt are indirectly regulated by miR-26a in a CKS2-dependent manner. Finally, CKS2 is overexpressed in PTC specimens relative to normal thyroid tissue, and a significant inverse correlation exists between miR-26a and CKS2 expression in clinical PTC specimens.

Conclusion

Our data indicate that miR-26a functions as a growth-suppressive miRNA in PTC, and that its suppressive effects are mediated mainly by repressing CKS2 expression.     

Tumorigenesis of Papillary Thyroid Cancer Is Not BRAF-Dependent in Patients with Acromegaly

Introduction

Several studies have reported a high frequency of papillary thyroid cancer (PTC) in patients with acromegaly. The aim of this study was to determine the prevalence and predictors of thyroid cancer in patients with acromegaly and to investigate the frequency of the BRAF ^V600E mutation in PTC patients with and without acromegaly.

Materials and Methods

We conducted a retrospective study of 60 patients with acromegaly. Thyroid ultrasonography (US) and US-guided fine needle aspiration were performed on nodules with sonographic features of malignancy. We selected 16 patients with non-acromegalic PTC as a control group. The BRAF ^V600E mutation was analyzed in paraffin-embedded surgical specimens of PTC by real-time polymerase chain reaction, and tumor specimens from patients with PTC were stained immunohistochemically with an antibody against insulin-like growth factor-1 receptor β (IGF-1Rβ).

Results

Thyroid cancer was found in 15 (25.0%) patients. No differences in age, sex, initial growth hormone (GH) and IGF-1 percentage of the upper limit of normal values or treatment modalities were observed between patients with and without PTC. Acromegaly was active in 12 of 15 patients at the time of PTC diagnosis; uncontrolled acromegaly had a significantly higher frequency in the PTC group (60%) than in the non-PTC group (28.9%) (p = 0.030). The BRAF ^V600E mutation was present in only 9.1% (1/11) of PTC patients with acromegaly, although 62.5% (10/16) of control patients with PTC had the mutation (p = 0.007). IGF-1Rβ immunostaining showed moderate-to-strong staining in all malignant PTC cells in patients with and without acromegaly. Significantly less staining for IGF-1Rβ was observed in normal adjacent thyroid tissues of PTC patients with acromegaly compared with those without (p = 0.014).

Conclusion

The prevalence of PTC in acromegalic patients was high (25%). An uncontrolled hyperactive GH-IGF-1 axis may play a dominant role in the development of PTC rather than the BRAF ^V600E mutation in patients with acromegaly.     

A survey of tandem repeat instabilities and associated gene expression changes in 35 colorectal cancers

Background

Colorectal cancer is a major contributor to cancer morbidity and mortality. Tandem repeat instability and its effect on cancer phenotypes remain so far poorly studied on a genome-wide scale.

Results

Here we analyze the genomes of 35 colorectal tumors and their matched normal (healthy) tissues for two types of tandem repeat instability, de-novo repeat gain or loss and repeat copy number variation. Specifically, we study for the first time genome-wide repeat instability in the promoters and exons of 18,439 genes, and examine the association of repeat instability with genome-scale gene expression levels. We find that tumors with a microsatellite instable (MSI) phenotype are enriched in genes with repeat instability, and that tumor genomes have significantly more genes with repeat instability compared to healthy tissues. Genes in tumor genomes with repeat instability in their promoters are significantly less expressed and show slightly higher levels of methylation. Genes in well-studied cancer-associated signaling pathways also contain significantly more unstable repeats in tumor genomes. Genes with such unstable repeats in the tumor-suppressor p53 pathway have lower expression levels, whereas genes with repeat instability in the MAPK and Wnt signaling pathways are expressed at higher levels, consistent with the oncogenic role they play in cancer.

Conclusions

Our results suggest that repeat instability in gene promoters and associated differential gene expression may play an important role in colorectal tumors, which is a first step towards the development of more effective molecular diagnostic approaches centered on repeat instability.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1902-9) contains supplementary material, which is available to authorized users.     

An Ep-ICD Based Index Is a Marker of Aggressiveness and Poor Prognosis in Thyroid Carcinoma

Background

Nuclear accumulation of the intracellular domain of epithelial cell adhesion molecule (Ep-ICD) in tumor cells was demonstrated to predict poor prognosis in thyroid carcinoma patients in our earlier study. Here, we investigated the clinical significance of Ep-ICD subcellular localization index (ESLI) in distinguishing aggressive papillary thyroid carcinoma (PTC) from non-aggressive cases.

Methods

Using domain specific antibodies against the intracellular (Ep-ICD) and extracellular (EpEx) domains of epithelial cell adhesion molecule, 200 archived tissues from a new cohort of patients with benign thyroid disease as well as malignant aggressive and non aggressive PTC were analyzed by immunohistochemistry (IHC). ESLI was defined as sum of the IHC scores for accumulation of nuclear and cytoplasmic Ep-ICD and loss of membranous EpEx; ESLI = [Ep-ICD_nuc + Ep-ICD_cyt + loss of membranous EpEx].

Results

For the benign thyroid tissues, non-aggressive PTC and aggressive PTC, the mean ESLI scores were 4.5, 6.7 and 11 respectively. Immunofluorescence double staining confirmed increased nuclear Ep-ICD accumulation and decreased membrane EpEx expression in aggressive PTC. Receiver-operating characteristic (ROC) curve analysis showed an area under the curve (AUC) of 0.841, 70.2% sensitivity and 83.9% specificity for nuclear Ep-ICD for differentiating aggressive PTC from non-aggressive PTC. ESLI distinguished aggressive PTC from non-aggressive cases with improved AUC of 0.924, 88.4% sensitivity and 85.5% specificity. Our study confirms nuclear accumulation of Ep-ICD and loss of membranous EpEx occurs in aggressive PTC underscoring the potential of Ep-ICD and ESLI to serve as diagnostic markers for aggressive PTC. Kaplan Meier survival analysis revealed significantly reduced disease free survival (DFS) for ESLI positive (cutoff >10) PTC (p<0.05), mean DFS = 133 months as compared to 210 months for patients who did not show positive ESLI.

Conclusion

ESLI scoring improves the identification of aggressive PTC and thereby may serve as a useful index for defining aggressiveness and poor prognosis among PTC patients.     

Comparative DNA methylome analysis of endometrial carcinoma reveals complex and distinct deregulation of cancer promoters and enhancers

Background

Aberrant DNA methylation is a hallmark of many cancers. Classically there are two types of endometrial cancer, endometrioid adenocarcinoma (EAC), or Type I, and uterine papillary serous carcinoma (UPSC), or Type II. However, the whole genome DNA methylation changes in these two classical types of endometrial cancer is still unknown.

Results

Here we described complete genome-wide DNA methylome maps of EAC, UPSC, and normal endometrium by applying a combined strategy of methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylation-sensitive restriction enzyme digestion sequencing (MRE-seq). We discovered distinct genome-wide DNA methylation patterns in EAC and UPSC: 27,009 and 15,676 recurrent differentially methylated regions (DMRs) were identified respectively, compared with normal endometrium. Over 80% of DMRs were in intergenic and intronic regions. The majority of these DMRs were not interrogated on the commonly used Infinium 450K array platform. Large-scale demethylation of chromosome X was detected in UPSC, accompanied by decreased XIST expression. Importantly, we discovered that the majority of the DMRs harbored promoter or enhancer functions and are specifically associated with genes related to uterine development and disease. Among these, abnormal methylation of transposable elements (TEs) may provide a novel mechanism to deregulate normal endometrium-specific enhancers derived from specific TEs.

Conclusions

DNA methylation changes are an important signature of endometrial cancer and regulate gene expression by affecting not only proximal promoters but also distal enhancers.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-868) contains supplementary material, which is available to authorized users.     

Simplification to Abacavir/Lamivudine + Atazanavir Maintains Viral Suppression and Improves Bone and Renal Biomarkers in ASSURE,a Randomized,Open Label,Non-Inferiority Trial

Objective

Simplification of antiretroviral therapy in patients with suppressed viremia may minimize long-term adverse effects. The study’s primary objective was to determine whether abacavir/lamivudine + atazanavir (ABC/3TC+ATV) was virologically non-inferior to tenofovir/emtricitabine + atazanavir/ritonavir (TDF/FTC+ATV/r) over 24 weeks in a population of virologically suppressed, HIV-1 infected patients.

Design

This open-label, multicenter, non-inferiority study enrolled antiretroviral experienced, HIV-infected adults currently receiving a regimen of TDF/FTC+ATV/r for ≥6 months with no history of virologic failure and whose HIV-1 RNA had been ≤75 copies/mL on 2 consecutive measurements including screening. Patients were randomized 1∶2 to continue current treatment or simplify to ABC/3TC+ATV.

Methods

The primary endpoint was the proportion of patients with HIV-RNA<50 copies/mL at Week 24 by the Time to Loss of Virologic Response (TLOVR) algorithm. Secondary endpoints included alternative measures of efficacy, adverse events (AEs), and fasting lipids. Exploratory endpoints included inflammatory, coagulation, bone, and renal biomarkers.

Results

After 24 weeks, ABC/3TC+ATV (n = 199) was non-inferior to TDF/FTC+ATV/r (n = 97) by both the primary analysis (87% in both groups) and all secondary efficacy analyses. Rates of grade 2–4 AEs were similar between the two groups (40% vs 37%, respectively), but an excess of hyperbilirubinemia made the rate of grade 3–4 laboratory abnormalities higher in the TDF/FTC+ATV/r group (30%) compared with the ABC/3TC+ATV group (13%). Lipid levels were stable except for HDL cholesterol, which increased significantly in the ABC/3TC+ATV group. Bone and renal biomarkers improved significantly between baseline and Week 24 in patients taking ABC/3TC+ATV, and the difference between groups was significant at Week 24. No significant changes occurred in any inflammatory or coagulation biomarker within or between treatment groups.

Conclusions

After 24 weeks, simplification to ABC/3TC+ATV from TDF/FTC+ATV/r maintained viral suppression was well-tolerated, and led to improvements in bone and renal biomarkers and HDL cholesterol.

Trial Registration

ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT01102972","term_id":"NCT01102972"}}NCT01102972 GlaxoSmithKline Clinical Study Register #113734     

Hashimoto's Thyroiditis: Similar and Dissimilar Characteristics in Neighboring Areas. Possible Implications for the Epidemiology of Thyroid Cancer

Context

Medical centers worldwide report an increased frequency of Hashimoto''s thyroiditis (HT) and thyroid cancer (TC), two environmentally influenced diseases. In Sicily, data on HT are available for the province of Messina (1975–2005); data on TC are available for the whole island (2002–2004), with the volcanic province of Catania having the highest incidence.

Objective

To replicate in Catania, on comparable years, the HT data of Messina.

Design, Methods, Setting

Review of the clinical records of patients in years 1995–2005 to compare presentation and yearly changes of HT. During 1995–2005, records were computer stored in the Endocrine Divisions of the University Hospitals of Catania and Messina, two tertiary referral centers.

Results

Catania is outnumbered by Messina (742 vs. 3,409 HT patients). Similar were the linear increase in the yearly number of HT patients, rates of thyroid dysfunctions though with different proportions of subclinical and overt hypothyroidism, and rates of positiveness for TgAb or TPOAb. Different were age and its yearly trend; gender distribution and rates of the sonography variants, though yearly trends were similar.

Conclusion

The HT epidemics is smaller in Catania, with changes in presentation overlapping partially those in Messina. Whatever environmental factors might be involved, they (and/or their intensity) were not necessarily the same in these provinces. Intriguingly, the expected number of TC in HT patients with thyroid nodules in Catania is congruent with that of the general population of this province, but it is far less than in the Messina province. Thus, TC and HT incidences could be influenced by distinct environmental factors.     

An Integrated Approach to Uncover Driver Genes in Breast Cancer Methylation Genomes

Background

Cancer cells typically exhibit large-scale aberrant methylation of gene promoters. Some of the genes with promoter methylation alterations play “driver” roles in tumorigenesis, whereas others are only “passengers”.

Results

Based on the assumption that promoter methylation alteration of a driver gene may lead to expression alternation of a set of genes associated with cancer pathways, we developed a computational framework for integrating promoter methylation and gene expression data to identify driver methylation aberrations of cancer. Applying this approach to breast cancer data, we identified many novel cancer driver genes and found that some of the identified driver genes were subtype-specific for basal-like, luminal-A and HER2+ subtypes of breast cancer.

Conclusion

The proposed framework proved effective in identifying cancer driver genes from genome-wide gene methylation and expression data of cancer. These results may provide new molecular targets for potential targeted and selective epigenetic therapy.     

MiR-886-3p regulates cell proliferation and migration, and is dysregulated in familial non-medullary thyroid cancer

Background

The molecular basis and characteristics of familial non-medullary thyroid cancer are poorly understood. In this study, we performed microRNA (miRNA) profiling of familial and sporadic papillary thyroid cancer tumor samples.

Methodology/Principal Findings

Genome wide miRNA profiling of sporadic and familial papillary thyroid cancer was performed. Differentially expressed miRNAs were validated by quantitative RT-PCR. Ectopic expression of miR-886-3p in thyroid cancer lines was performed to identify pathways targeted by the miRNA, as well as, to determine its effect on tumor cell biology. We found four differentially expressed miRNAs between familial and sporadic papillary thyroid cancer tumor samples. MiR-886-3p and miR-20a were validated to be differentially expressed by 3- and 4-fold, respectively. Pathway analysis of genome-wide expression data on cells overexpressing miR-886-3p and target prediction analysis showed genes involved in DNA replication and focal adhesion pathways were regulated by miR-886-3p. Overexpression of miR-886-3p in thyroid cancer cell lines significantly inhibited cellular proliferation, the number and size of spheroids and cellular migration. Additionally, overexpression of miR-886-3p increased the number of cells in S phase.

Conclusions/Significance

Our findings for the first time suggest that miR-886-3p plays an important role in thyroid cancer tumor cell biology and regulates genes involved in DNA replication and focal adhesion. Thus, miR-886-3p may play a role in the initiation and or progression of papillary thyroid cancer.     

An Appraisal of Proliferation and Apoptotic Markers in Papillary Thyroid Carcinoma: An Automated Analysis

Introduction

Proliferation and apoptosis are opposing processes by which the cell numbers are kept in a delicate balance, essential for tissue homeostasis, whereas uncontrolled growth of cells is a hallmark of cancer. Papillary thyroid cancer (PTC) is the commonest type of thyroid cancer, with some PTC following an indolent course, whereas the other ones are more aggressive.

Aim

To evaluate respective contribution of proliferation and apoptosis in the tumorigenesis of PTC by automated analysis.

Materials and Methods

We investigated the immunolabeling of phosphorylated histone H3 (pHH3), cyclin D1, active caspase-3, and bcl-2 in thirteen cases each of metastatic PTC, follicular variant of PTC (FVPTC), papillary microcarcinoma (PMC) and well differentiated tumor of uncertain malignant potential (WDT-UMP). FVPTC cases comprised seven encapsulated and six unencapsulated cases.

Results

Proliferation, as assessed by pHH3 and cyclin D1 immunolabeling, was increased in all PTC variants, including the putative precursor lesion WDT-UMP, compared to normal thyroid tissue. pHH3 was immunolabeled in more cells of metastatic PTC than of PMC and of encapsulated FVPTC. Surprisingly, metastatic PTC and unencapsulated FVPTC also demonstrated more cleaved caspase-3 immunolabeled cells than the other types. In contrast, increased expression of bcl-2 protein was seen in normal thyroid areas, encapsulated FVPTC and PMC as compared to metastatic PTC. Metastatic PTC shows higher proliferation than other types of PTC but unexpectedly also higher apoptotic levels. Similar results were also seen with unencapsulated FVPTC, thus suggesting that unencapsulated FVPTC has a potential for adverse outcome. Bcl-2 was immunolabeled in a low percentage of cells in WDT-UMP.

Conclusions

The expression of the proliferative protein pHH3 together with the apoptotic marker cleaved caspase-3 may indicate an aggressive behaviour of PTC and loss of apoptosis inhibition by bcl-2 protein can further amplify the role of these proteins in tumor progression. Both cyclin D1 and bcl-2 could prove to be interesting markers of PTC precursor lesions. Automated/digital image quantification approach helps in refining the diagnostic accuracy.     

MiRNA Expression May Account for Chronic but Not for Acute Regulation of mRNA Expression in Human Thyroid Tumor Models

Background

For thyroid tumorigenesis, two main human in vitro models are available: primary cultures of human thyrocytes treated with TSH or EGF/serum as models for autonomous adenomas (AA) or papillary thyroid carcinomas (PTC) respectively, and human thyroid tumor derived cell lines. Previous works of our group have assessed properties of those models, with a special emphasis on mRNA regulations. It is often assumed that miRNA may be one of the primary events inducing these mRNA regulations.

Methods

The purpose of this study was to investigate the representativity of those models to study microRNA regulations and their relation with mRNA expression. To achieve this aim, the miRNA expressions profiles of primary cultures treated with TSH or EGF/serum and of 6 thyroid cancer cell lines were compared to the expression profiles of 35 tumor tissues obtained by microarrays.

Results

Our data on primary cultures have shown that the TSH or EGF/serum treatment did not greatly modify the microRNA expression profiles, which is contrary to what is observed for mRNA expression profiles, although they still evolved differently according to the treatment. The analysis of miRNA and mRNA expressions profiles in the cell lines has shown that they have evolved into a common, dedifferentiated phenotype, closer to ATC than to the tumors they are derived from.

Conclusions

Long-terms TSH or EGF/serum treatments do not mimic AA or PTC respectively in terms of miRNA expression as they do for mRNA, suggesting that the regulations of mRNA expression induced by these physiological agents occur independently of miRNA. The general patterns of miRNA expression in the cell lines suggest that they represent a useful model for undifferentiated thyroid cancer. Mirna probably do not mediate the rapid changes in gene expression in rapid cell biology regulation.     

Inferring the perturbed microRNA regulatory networks from gene expression data using a network propagation based method

Background

MicroRNAs (miRNAs) are a class of endogenous small regulatory RNAs. Identifications of the dys-regulated or perturbed miRNAs and their key target genes are important for understanding the regulatory networks associated with the studied cellular processes. Several computational methods have been developed to infer the perturbed miRNA regulatory networks by integrating genome-wide gene expression data and sequence-based miRNA-target predictions. However, most of them only use the expression information of the miRNA direct targets, rarely considering the secondary effects of miRNA perturbation on the global gene regulatory networks.

Results

We proposed a network propagation based method to infer the perturbed miRNAs and their key target genes by integrating gene expressions and global gene regulatory network information. The method used random walk with restart in gene regulatory networks to model the network effects of the miRNA perturbation. Then, it evaluated the significance of the correlation between the network effects of the miRNA perturbation and the gene differential expression levels with a forward searching strategy. Results show that our method outperformed several compared methods in rediscovering the experimentally perturbed miRNAs in cancer cell lines. Then, we applied it on a gene expression dataset of colorectal cancer clinical patient samples and inferred the perturbed miRNA regulatory networks of colorectal cancer, including several known oncogenic or tumor-suppressive miRNAs, such as miR-17, miR-26 and miR-145.

Conclusions

Our network propagation based method takes advantage of the network effect of the miRNA perturbation on its target genes. It is a useful approach to infer the perturbed miRNAs and their key target genes associated with the studied biological processes using gene expression data.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-255) contains supplementary material, which is available to authorized users.