Skip to main content

scdNet: a computational tool for single-cell differential network analysis

Abstract

Background

Single-cell RNA sequencing (scRNA-Seq) is an emerging technology that has revolutionized the research of the tumor heterogeneity. However, the highly sparse data matrices generated by the technology have posed an obstacle to the analysis of differential gene regulatory networks.

Results

Addressing the challenges, this study presents, as far as we know, the first bioinformatics tool for scRNA-Seq-based differential network analysis (scdNet). The tool features a sample size adjustment of gene-gene correlation, comparison of inter-state correlations, and construction of differential networks. A simulation analysis demonstrated the power of scdNet in the analyses of sparse scRNA-Seq data matrices, with low requirement on the sample size, high computation efficiency, and tolerance of sequencing noises. Applying the tool to analyze two datasets of single circulating tumor cells (CTCs) of prostate cancer and early mouse embryos, our data demonstrated that differential gene regulation plays crucial roles in anti-androgen resistance and early embryonic development.

Conclusions

Overall, the tool is widely applicable to datasets generated by the emerging technology to bring biological insights into tumor heterogeneity and other studies. MATLAB implementation of scdNet is available at https://github.com/ChenLabGCCRI/scdNet.

Background

Single-cell sequencing is a developing technology that enables a close look into the heterogeneity and clonal evolution of cancer cells. While many methods have been designed to analyze single-cell DNA-Seq data [1], the analysis of scRNA-Seq data remains challenging due to high sparsity that prevents direct applications of methods originally developed for microarray and bulk RNA sequencing. On the other hand, the analysis of gene regulatory networks is a widely used approach to realize the signaling and interactions among genes. Recent studies have successfully applied correlation onto the inference of gene regulatory networks by using scRNA-Seq data [2, 3]. However, realizing that tumor cells are highly heterogeneous, network topologies may be massively changed between cells of different cellular states [4]. Recently, a computational method was proposed to investigate the change in mean absolute distances of a gene to others [5]. However, the computational method for studying individual gene pairs of which regulatory strengths alter between conditions was only carried out in the bulk RNA sequencing data [6,7,8].

Addressing this research need, here we developed a comprehensive bioinformatics tool for single cell-based differential network analysis, namely scdNet. It features two main functions: i) gene correlation analysis out of highly sparse data matrices and ii) differential network analysis between cellular states. Performance of scdNet was tested by simulated datasets. We further applied the tool to scRNA-Seq datasets of CTCs and early-stage mouse embryos for differential networks associated with anti-androgen resistance of prostate cancer and early embryonic development.

Methods

Transformation of intra-state gene-gene correlation

Suppose a scRNA-Seq dataset EG × K = {e(g, k)} contains read counts normalized by DESeq2 [9] of the g-th genes in the k-th single cells (g [1, G] and k [1, K]). We z-transformed E with respect to genes to eliminate biases: \( {\boldsymbol{E}}_{\boldsymbol{G}\times \boldsymbol{K}}^{\boldsymbol{z}}=\left\{\frac{e\left(g,k\right)-{\mu}_g}{\sigma_g}\right\} \), where μg and σg are mean and standard deviation values of gene g across K cells. For the sparsity of scRNA-Seq data, we adopted a sample size adjustment and comparison of Pearson correlation coefficients among states. Within a cellular state n (n {0, 1}), gene-gene correlation coefficients were computed into the correlation matrix \( {\boldsymbol{C}}_{\boldsymbol{G}\times \boldsymbol{G}}^{\boldsymbol{n}} \):

$$ {\boldsymbol{C}}^{\boldsymbol{n}}\left(i,j\right)=\rho \left(\boldsymbol{E}\left(i,{\boldsymbol{x}}_{\boldsymbol{i},\boldsymbol{j}}\right),\boldsymbol{E}\left(j,{\boldsymbol{x}}_{\boldsymbol{i},\boldsymbol{j}}\right)\right), $$
(1)

where (i, j)  [1, G] × [1, G] and xi, j = {s| E(i, s) ≠ 0} ∩ {t| E(j, t) ≠ 0} for samples of the state. Number of samples used to calculate the correlation was stored in the sample size matrix \( {\boldsymbol{S}}_{\boldsymbol{G}\times \boldsymbol{G}}^{\boldsymbol{n}} \):

$$ {\boldsymbol{S}}^{\boldsymbol{n}}\left(i,j\right)=\parallel {\boldsymbol{x}}_{\boldsymbol{i},\boldsymbol{j}}\parallel . $$
(2)

We applied the Fisher transformation \( \mathcal{F} \), as described in our previous studies [6,7,8], to convert the correlation coefficients into a sample size-free domain and termed as the interaction matrix \( {\boldsymbol{I}}_{\boldsymbol{G}\times \boldsymbol{G}}^{\boldsymbol{n}} \):

$$ {\boldsymbol{I}}^{\boldsymbol{n}}\left(i,j\right)=\mathcal{F}\left({\boldsymbol{C}}^{\boldsymbol{n}}\left(i,j\right),{\boldsymbol{S}}^{\boldsymbol{n}}\left(i,j\right)\right)=\frac{\sqrt{{\boldsymbol{S}}^{\boldsymbol{n}}\left(i,j\right)-3}}{2}\ln \frac{1+{\boldsymbol{C}}^{\boldsymbol{n}}\left(i,j\right)}{1-{\boldsymbol{C}}^{\boldsymbol{n}}\left(i,j\right)}. $$
(3)

Elements of the interaction matrix followed the standard normal distribution.

Inter-state differential network analysis

We analyzed the changes of interaction strengths of each gene-gene pair, say i and j, between two cellular states:

$$ \mathbf{\Delta }\boldsymbol{I}\left(i,j\right)=\left|{\boldsymbol{I}}^{\mathbf{1}}\left(i,j\right)\right|-\left|{\boldsymbol{I}}^{\mathbf{0}}\left(i,j\right)\right| $$
(4)

and tested for statistical significance by the cumulative distribution function (CDF) [8]:

$$ F\left(\mathbf{\Delta }\boldsymbol{I}\left(i,j\right)\right)=\frac{1}{2}+\mathit{\operatorname{erf}}\left(\frac{\mathbf{\Delta }\boldsymbol{I}\left(i,j\right)}{2}\right)-\frac{1}{2}\mathit{\operatorname{sgn}}\left(\mathbf{\Delta }\boldsymbol{I}\left(i,j\right)\right)\bullet {\left[\mathit{\operatorname{erf}}\left(\frac{\mathbf{\Delta }\boldsymbol{I}\left(i,j\right)}{2}\right)\right]}^2 $$
(5)

, where erf() and sgn() are the Gauss error function and sign function, respectively. Given the exact CDF of the changes in the interaction values, P-values are directly assessed and time-consuming permutation or simulation tests can be avoided. Gene pairs with significant changes were defined as differential pairs and merged into a differential network.

Visualization and functional annotation analysis of gene regulatory networks

Networks of identified genes and their dynamic interactions across cellular states were visualized by the Cytoscape software (version 3.5.1) [10], with nodes and edges denoting genes and differential gene regulations, respectively. To investigate the functional relevance of cellular state-modulated differential networks, top hub genes of the network were analyzed for the associations with Gene Ontology (GO) terms of molecular functions, cellular components, and biological processes by the Database for Annotation, Visualization and Integrated Discovery (DAVID) [11].

Results

Model overview

scdNet is devised to analyze differential gene regulatory networks associated with cellular states at the single cell level. Fig. 1 shows the flowchart of the tool. Briefly, scRNA-Seq data were preprocessed and normalized to eliminate inter-cell biases, and non-informative genes with a coefficient of variation < 0.25 in either state were eleminated. Within each group of cells, gene-gene correlation coefficients were calculated with an exclusion of zeros, and transformed to a sample size independent domain by the Fisher transformation to eliminate sample size related biases. Normalized correlation coefficients were compared between groups of cells and the changes in correlation were statistically tested in the Fisher domain. Significantly changed gene-gene pairs were merged into a differential network. Visualization and functional annotation analyses were performed to realize the biological relevance of such dynamic network. MATLAB implementation of scdNet is available at https://github.com/ChenLabGCCRI/scdNet.

Fig. 1
figure 1

Flowchart of the proposed method. The method is designed to analyze differential gene regulatory networks from scRNA-Seq data. It features two functions: i) measurement of sample size corrected gene-gene correlation for each state to handle the sparse data matrices and ii) statistical inference of the changes in correlation across cellular states. The identified differential gene-gene pairs were subject to network and functional annotation analyses

Assessment of model performace – simulation design

We simulated scRNA-Seq datasets to test the performance of the proposed method in analyzing gene-gene correlations out of highly sparse scRNA-Seq data. Since we z-transformed the sequencing data, the synthetic datasets were generated by randomly sampling the standard normal distribution. We note that the performance of the comparison of interaction matrices (Eqs. 45) has been described in our previous papers [7, 8]. Thus, the simulation was simply focused on the comparison of intra-group correlations against zero (uncorrelated). In each simulation scenario, 10,000 gene pairs were generated, of which 20% were defined as correlated (covariance=θ), and the remaining 80% as independent (covariance = 0). Four parameters were tuned in the simulation: i) covariance of correlated gene pairs (θ, 0, 0.3, 0.7, and 1.0), ii) relative power of Gaussian noises added to the original signals to mimic sequencing errors (ε, 0, 0.05, 0.1, 0.25, 0.50, 0.75, 0.9, and 1.0), iii) number of single cells (K, 10, 20, 50,100, and 200), and iv) proportions of low-signal elements to be eliminated (representing zeroes of scRNA-Seq) (τ, 0, 0.25, 0.50, 0.75. 0.9, and 0.95). Yielded Fisher-transformed scores were compared to zero; gene pairs with Bonferroni-adjusted P-value< 0.05 were called as significant. Performance was evaluated by accuracy, sensitivity, specificity, and time consumption. The simulation processes were performed on MATLAB.

Assessment of model performace by four parameters

Gene-gene Covariance (θ)

With ε, K, and τ set at 0.05, 50, 0.25, respectively, we tested the performance of scdNet across different gene-gene covariance θ. As shown in Fig. 2a, the increase in θ greatly rose sensitivity from to 0.14 to 1.00 (θ=0 and 1.0) and slightly increased accuracy from 0.70 to 0.88. On the other hand, specificity seemed to be independent of θ (range, 0.84–0.85).

Fig. 2
figure 2

Simulation analysis for performance assessment with respect to gene-gene covariance and relative power of noise. a Performance of the tool on a dataset composed of 20% of correlated samples with a covariance ranging from 0 to 1 and 80% of uncorrelated samples. b Performance of the tool when Gaussian noises of different power relative to the original signals were added

Relative Power of Noises (ε)

We then tested the effect of sequencing noises. Simulation settings were identical as described above, while θ was fixed at 0.7. At a general level of sequencing errors (ε=0.1), scdNet achieved high performance (accuracy, sensitivity, and specificity = 0.86, 0.90, and 0.85). While accuracy and specificity were not sensitive to ε, the sensitivity dropped to 0.40 when the noise was extremely high (Fig. 2b).

Number of Single Cells (K) and Proportion of Low-signal Elements (τ)

Due to limitations in budget and specimens, scRNA-Seq data are typically of limited sample sizes. To realize the effects of sample size and sparsity of scRNA-Seq data, we jointly analyzed the two factors. As shown in Table 1, the proposed tool achieved generally favorable performance in regardless of settings of K and τ when either group of cells had 10 or more non-zero values.

Table 1 Simulation analysis on the number of single cells and proportions of low-signal elements

Time Consumption

Facilitated by the exact probability function (Eq. 5), scdNet is of remarkably high computation efficiency. An inference of 10,000 gene pairs described above took an average of 1.3 s (std., 0.5) on a personal computer.

Application to a prostate cancer dataset of CTCs

We then applied scdNet to a real dataset derived from prostate CTCs. CTCs is an emerging technology of liquid biopsies that allows minimally invasive assessment and prediction of metastasis and treatment outcomes of advanced prostate cancer [12, 13]. Here we studied on the resistance of enzalutamide, a second-generation anti-androgen medication which is used in the treatment of prostate cancer [14, 15] while its resistance mechanisms remain vague. Thus, we utilized a public dataset of 169 scRNA-Seq of prostate CTCs (Gene Expression Omnibus accession number: GSE67980) [16]. We normalized raw read counts by DESeq2. The dataset was found very sparse. Among 21,696 unique genes, each cell carried an average of ~ 76.4% genes (16,573, std., 2293) with no sequencing reads. We analyzed 77 samples isolated from 13 patients, of which 41 progressed on enzalutamide (hereafter referred to as the enzalutamide-resistant group) and 36 were enzalutamide-naïve. Out of ~ 4.7 million transcriptome-wide gene pairs, we set a stringent criterion to identify the most significant subset of pairs that exhibited significant changes between the two groups (Bonferroni adjusted P < 1 × 10− 5). In total, 6023 and 10,670 pairs of genes were correlated with each other specifically in enzalutamide-resistant and -naïve groups, respectively, involving 2735 genes. These gene pairs formed a highly intertwined network (Fig. 3a); on average, each gene was connected to 12.2 partner genes. We note that only 2.8% of these genes were differentially expressed (with t-test P < 0.05) between the two groups of cells, confirming that our analysis was not biased by differential expression. (Fig. 3b). Table 2 tabulates the top 10 hub genes of the network. ENOSF1, a mitochondria enzyme that have been shown as a serum biomarker for gastric cancer [17], was the top hub with 150 differential pairs (Fig. 3c). Consistently, functional annotation analysis of the top 100 hub genes in the network also highlighted the association between mitochondria-related biological functions and anti-androgen resistance (Table 3). Thus, our data indicate differential gene regulatory networks, at least partially, govern mitochondria functions and play an essential role in anti-androgen resistance of prostate cancer.

Fig. 3
figure 3

Differential gene networks associated with anit-androgen resistance of prostate cancer. a Differential gene regulatory network. The network was constructed by merging differential gene pairs between enzalutamide-resistant or -naïve CTCs, with nodes and edges representing genes and differential correlations, respectively. Top hub genes are labeled with gene symbols. b Venn diagram of genes involved in the differential network and those differentially expressed between two groups of cells. c Subnetwork of the top hub gene ENOSF1

Table 2 Top hub genes in the differential gene network between enzalutamide-resistant and -naïve CTCs
Table 3 Gene Ontology terms associated with top 100 hub genes of the enzalutamide resistance-modulated differential network

Application to a dataset of early mouse embryos

We also applied scdNet to study the differential network associated with early development of mouse embryos. We preprocessed scRNA-Seq data of 32 single cells at the 8-cell stage and 16 cells at the 2-cell stage (ArrayExpress accession number E-MTAB-3321 [18]) as described in the CTC study. The data were also highly sparse (68.6% sparsity). Using identical criteria, we identified 11,245 gene regulatory pairs specifically shown in 2-cell (2915 pairs) and 8-cell (8330) embryos composed of 3998 genes (Fig. 4a). Again, the proportion of these genes to be differentially expressed was lower than expectation (39.1% compared to an expectation of 49.0%; Fig. 4b). The top 100 hubs genes were significantly associated with cell division and differentiation-related functions, such as cytoskeleton and chromatin assembly, and the ribosome and mitochondria, which are known to regulate early development of mouse embryos [19, 20] (Table 4).

Fig. 4
figure 4

Differential gene networks associated with early development of mouse embryos. a Differential gene regulatory network of differential gene pairs identified by comparing the 2-cell and 8-cell stages of mouse embryonic development. Top hub genes are labeled with gene symbols. b Venn diagram of genes involved in the differential network and those differentially expressed between two groups of cells

Table 4 Gene Ontology terms associated with top 100 hub genes of the early embryonic development-modulated differential network

Discussion

Recent advances in scRNA-Seq technology have revolutionized the investigation of tumor heterogeneity [21,22,23] and construction of the cell atlas in human [24] and mouse [25]. This study addresses the unmet demand for studying condition-specific gene regulatory network using scRNA-Seq data. We proposed and implemented a novel bioinformatics algorithm, scdNet, for systematic identification of gene pairs with regulatory strengths significantly changed between two groups of single cells. We adopted a sample-size correction transformation on correlation coefficients to cope with the sparsity of scRNA-Seq data. Using simulated datasets, we demonstrated the tolerance of scdNet to gene-gene covariances, relative power of noises, and number of single cells, as well as its great computational efficiency. We also applied the method to study two real-world datasets. CTC is a minimally invasive liquid biopsy strongly indicated for investigation of metastasis [26], prediction of treatment response [12, 13], and risk assessment of cancers [27, 28]. Our work is a unique extension into the underlying mechanisms of CTCs in treatment response of metastatic prostate cancer. By comparing CTCs obtained from patients responsive and naïve to a second-generation anti-androgen therapy, we constructed an intertwined gene regulatory network. Our data are in line with previous in vitro studies in the critical role of mitochondria and oxidative stress in the development of hormone-refractory prostate cancer [29,30,31]. In a non-cancer setting, scdNet also identified biologically meaningful results related to early development of mouse embryos. Overall, we demonstrated the feasibility of our method using simulated and real datasets. We expect the method to be widely applicable to different studies of biomedicine with the emerging applications of scRNA-Seq.

Conclusions

Here we present a novel bioinformatics tool, namely scdNet, for a fast and comprehensive inference of differential gene regulatory networks out of scRNA-Seq data. Performance and computation efficiency of scdNet were demonstrated by simulation analysis. Applying the tool to a dataset of prostate cancer, we showed the involvement of mitochondria-related biological functions in anti-androgen resistance. We also illuminated crucial biological functions regulating early development of mouse embryos. Taken together, our data suggest wide applications of scdNet in exploring differential networks out of the rapidly increasing scRNA-Seq studies.

Abbreviations

CTC:

Circulating tumor cell

GO:

Gene Ontology

scdNet:

scRNA-Seq-based differential network analysis

scRNA-Seq:

Single-cell RNA sequencing

References

  1. Gawad C, Koh W, Quake SR. Single-cell genome sequencing: current state of the science. Nat Rev Genet. 2016;17(3):175–88.

    Article  CAS  Google Scholar 

  2. Crow M, Paul A, Ballouz S, Huang ZJ, Gillis J. Exploiting single-cell expression to characterize co-expression replicability. Genome Biol. 2016;17:101.

    Article  Google Scholar 

  3. Wang J, Xia S, Arand B, Zhu H, Machiraju R, Huang K, Ji H, Qian J. Single-cell co-expression analysis reveals distinct functional modules, co-regulation mechanisms and clinical outcomes. PLoS Comput Biol. 2016;12(4):e1004892.

    Article  Google Scholar 

  4. Ideker T, Krogan NJ. Differential network biology. Mol Syst Biol. 2012;8:565.

    Article  Google Scholar 

  5. Wang Y, Wu H, Yu T. Differential gene network analysis from single cell RNA-seq. J Genet Genomics. 2017;44(6):331–4.

    Article  Google Scholar 

  6. Chiu YC, Hsiao TH, Wang LJ, Chen Y, Chuang EY. Analyzing differential regulatory networks modulated by continuous-state genomic features in Glioblastoma Multiforme. IEEE/ACM Trans Comput Biol Bioinform. 2016.

  7. Chiu YC, Wang LJ, Lu TP, Hsiao TH, Chuang EY, Chen Y. Differential correlation analysis of glioblastoma reveals immune ceRNA interactions predictive of patient survival. BMC Bioinformatics. 2017;18(1):132.

    Article  Google Scholar 

  8. Hsiao TH, Chiu YC, Hsu PY, Lu TP, Lai LC, Tsai MH, Huang TH, Chuang EY, Chen Y. Differential network analysis reveals the genome-wide landscape of estrogen receptor modulation in hormonal cancers. Sci Rep. 2016;6:23035.

    Article  CAS  Google Scholar 

  9. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  Google Scholar 

  10. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  Google Scholar 

  11. Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  Google Scholar 

  12. Antonarakis ES, Lu C, Wang H, Luber B, Nakazawa M, Roeser JC, Chen Y, Mohammad TA, Chen Y, Fedor HL, et al. AR-V7 and resistance to enzalutamide and abiraterone in prostate cancer. N Engl J Med. 2014;371(11):1028–38.

    Article  Google Scholar 

  13. de Bono JS, Scher HI, Montgomery RB, Parker C, Miller MC, Tissing H, Doyle GV, Terstappen LW, Pienta KJ, Raghavan D. Circulating tumor cells predict survival benefit from treatment in metastatic castration-resistant prostate cancer. Clin Cancer Res. 2008;14(19):6302–9.

    Article  CAS  Google Scholar 

  14. Scher HI, Fizazi K, Saad F, Taplin ME, Sternberg CN, Miller K, de Wit R, Mulders P, Chi KN, Shore ND, et al. Increased survival with enzalutamide in prostate cancer after chemotherapy. N Engl J Med. 2012;367(13):1187–97.

    Article  CAS  Google Scholar 

  15. Taneja SS. Re: increased survival with enzalutamide in prostate cancer after chemotherapy. J Urol. 2013;189(1):123–4.

    Article  Google Scholar 

  16. Miyamoto DT, Zheng Y, Wittner BS, Lee RJ, Zhu H, Broderick KT, Desai R, Fox DB, Brannigan BW, Trautwein J, et al. RNA-Seq of single prostate CTCs implicates noncanonical Wnt signaling in antiandrogen resistance. Science. 2015;349(6254):1351–6.

    Article  CAS  Google Scholar 

  17. Yang J, Xiong X, Wang X, Guo B, He K, Huang C. Identification of peptide regions of SERPINA1 and ENOSF1 and their protein expression as potential serum biomarkers for gastric cancer. Tumour Biol. 2015;36(7):5109–18.

    Article  CAS  Google Scholar 

  18. Goolam M, Scialdone A, Graham SJL, Macaulay IC, Jedrusik A, Hupalowska A, Voet T, Marioni JC, Zernicka-Goetz M. Heterogeneity in Oct4 and Sox2 targets biases cell fate in 4-cell mouse embryos. Cell. 2016;165(1):61–74.

    Article  CAS  Google Scholar 

  19. Engel W, Zenzes MT, Schmid M. Activation of mouse ribosomal RNA genes at the 2-cell stage. Hum Genet. 1977;38(1):57–63.

    Article  CAS  Google Scholar 

  20. Mitchell M, Schulz SL, Armstrong DT, Lane M. Metabolic and mitochondrial dysfunction in early mouse embryos following maternal dietary protein intervention. Biol Reprod. 2009;80(4):622–30.

    Article  CAS  Google Scholar 

  21. Chung W, Eum HH, Lee HO, Lee KM, Lee HB, Kim KT, Ryu HS, Kim S, Lee JE, Park YH, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun. 2017;8:15081.

    Article  CAS  Google Scholar 

  22. Karaayvaz M, Cristea S, Gillespie SM, Patel AP, Mylvaganam R, Luo CC, Specht MC, Bernstein BE, Michor F, Ellisen LW. Unravelling subclonal heterogeneity and aggressive disease states in TNBC through single-cell RNA-seq. Nat Commun. 2018;9(1):3588.

    Article  Google Scholar 

  23. Fan J, Lee HO, Lee S, Ryu DE, Lee S, Xue C, Kim SJ, Kim K, Barkas N, Park PJ, et al. Linking transcriptional and genetic tumor heterogeneity through allele analysis of single-cell RNA-seq data. Genome Res. 2018;28(8):1217–27.

    Article  CAS  Google Scholar 

  24. Rozenblatt-Rosen O, Stubbington MJT, Regev A, Teichmann SA. The human cell atlas: from vision to reality. Nature. 2017;550(7677):451–3.

    Article  CAS  Google Scholar 

  25. Han X, Wang R, Zhou Y, Fei L, Sun H, Lai S, Saadatpour A, Zhou Z, Chen H, Ye F, et al. Mapping the mouse cell atlas by microwell-Seq. Cell. 2018;172(5):1091–107 e1017.

    Article  CAS  Google Scholar 

  26. Aceto N, Bardia A, Miyamoto DT, Donaldson MC, Wittner BS, Spencer JA, Yu M, Pely A, Engstrom A, Zhu H, et al. Circulating tumor cell clusters are oligoclonal precursors of breast cancer metastasis. Cell. 2014;158(5):1110–22.

    Article  CAS  Google Scholar 

  27. Moreno JG, Miller MC, Gross S, Allard WJ, Gomella LG, Terstappen LW. Circulating tumor cells predict survival in patients with metastatic prostate cancer. Urology. 2005;65(4):713–8.

    Article  Google Scholar 

  28. Rack B, Schindlbeck C, Juckstock J, Andergassen U, Hepp P, Zwingers T, Friedl TW, Lorenz R, Tesch H, Fasching PA, et al. Circulating tumor cells predict survival in early average-to-high risk breast cancer patients. J Natl Cancer Inst. 2014;106(5).

  29. Lee EC, Tenniswood MP. Emergence of metastatic hormone-refractory disease in prostate cancer after anti-androgen therapy. J Cell Biochem. 2004;91(4):662–70.

    Article  CAS  Google Scholar 

  30. Shiota M, Yokomizo A, Naito S. Oxidative stress and androgen receptor signaling in the development and progression of castration-resistant prostate cancer. Free Radic Biol Med. 2011;51(7):1320–8.

    Article  CAS  Google Scholar 

  31. Shiota M, Yokomizo A, Naito S. Pro-survival and anti-apoptotic properties of androgen receptor signaling by oxidative stress promote treatment resistance in prostate cancer. Endocr Relat Cancer. 2012;19(6):R243–53.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors greatly appreciate the brilliant and constructive inputs from reviewers and participants of the International Conference on Intelligent Biology and Medicine (ICIBM 2018).

Funding

Publication charges for this article have been funded by Taichung Veterans General Hospital. This research was also supported by the Ministry of Science and Technology, Taiwan (MOST104-2314-B-038-044) to YJS; the National Health Research Institutes, Taiwan (NHRI-EX107-10710BC) to THH; the NCI Cancer Center Shared Resources (NIH-NCI P30CA54174), NIH (CTSA 1UL1RR025767–01), and CPRIT (RP160732) to YC; and San Antonio Life Science Institute (SALSI Postdoctoral Research Fellowship 2018) to YCC. The funding sources had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

MATLAB implementation of scdNet is available at https://github.com/ChenLabGCCRI/scdNet.

About this supplement

This article has been published as part of BMC Systems Biology Volume 12 Supplement 8, 2018: Selected articles from the International Conference on Intelligent Biology and Medicine (ICIBM) 2018: systems biology. The full contents of the supplement are available online at https://bmcsystbiol.biomedcentral.com/articles/supplements/volume-12-supplement-8.

Author information

Authors and Affiliations

Authors

Contributions

YCC, THH, LJW, YC, and YJS conceived the study. YCC, THH, and YC designed the model. YCC and LJW performed data analysis. YCC, THH, YC, and YJS interpreted the data. YCC, THH, LJW, YC, and YJS wrote and approved the final version of paper.

Corresponding authors

Correspondence to Yidong Chen or Yu-Hsuan Joni Shao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chiu, YC., Hsiao, TH., Wang, LJ. et al. scdNet: a computational tool for single-cell differential network analysis. BMC Syst Biol 12 (Suppl 8), 124 (2018). https://doi.org/10.1186/s12918-018-0652-0

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/s12918-018-0652-0

Keywords