Skip to main content
  • Research article
  • Open access
  • Published:

Constructing an integrated genetic and epigenetic cellular network for whole cellular mechanism using high-throughput next-generation sequencing data

Abstract

Background

Epigenetics has been investigated in cancer initiation, and development, especially, since the appearance of epigenomics. Epigenetics may be defined as the mechanisms that lead to heritable changes in gene function and without affecting the sequence of genome. These mechanisms explain how individuals with the same genotype produce phenotypic differences in response to environmental stimuli. Recently, with the accumulation of high-throughput next-generation sequencing (NGS) data, a key goal of systems biology is to construct networks for different cellular levels to explore whole cellular mechanisms. At present, there is no satisfactory method to construct an integrated genetic and epigenetic cellular network (IGECN), which combines NGS omics data with gene regulatory networks (GRNs), microRNAs (miRNAs) regulatory networks, protein-protein interaction networks (PPINs), and epigenetic regulatory networks of methylation using high-throughput NGS data.

Results

We investigated different kinds of NGS omics data to develop a systems biology method to construct an integrated cellular network based on three coupling models that describe genetic regulatory networks, protein–protein interaction networks, microRNA (miRNA) regulatory networks, and methylation regulation. The proposed method was applied to construct IGECNs of gastric cancer and the human immune response to human immunodeficiency virus (HIV) infection, to elucidate human defense response mechanisms. We successfully constructed an IGECN and validated it by using evidence from literature search. The integration of NGS omics data related to transcription regulation, protein-protein interactions, and miRNA and methylation regulation has more predictive power than independent datasets. We found that dysregulation of MIR7 contributes to the initiation and progression of inflammation-induced gastric cancer; dysregulation of MIR9 contributes to HIV-1 infection to hijack CD4+ T cells through dysfunction of the immune and hormone pathways; dysregulation of MIR139-5p, MIRLET7i, and MIR10a contributes to the HIV-1 integration/replication stage; dysregulation of MIR101, MIR141, and MIR152 contributes to the HIV-1 virus assembly and budding mechanisms; dysregulation of MIR302a contributes to not only microvesicle-mediated transfer of miRNAs but also dysfunction of NF-κB signaling pathway in hepatocarcinogenesis.

Conclusion

The coupling dynamic systems of the whole IGECN can allow us to investigate genetic and epigenetic cellular mechanisms via omics data and big database mining, and are useful for further experiments in the field of systems and synthetic biology.

Background

With advances in molecular biology technologies for whole genome sequencing, expression profiling, and high-throughput experiments, large amounts of biological data covering various biological levels have emerged [1–3]. These kinds of ‘omics’ data include genetic sequences (genomics), microarray-based genome-wide expression profiles (transcriptomics), protein abundance data (proteomics), and microRNA (miRNA) and methylation data, and they provide an unprecedented view of cellular components and their cellular mechanisms in biological systems [4]. With the large amount of genomic/transcriptomic and epigenetic data that has accumulated, a more complex cellular function understanding of living organisms is possible. Computational systems biology techniques able to combine these large, heterogeneous genomic/transcriptomic and epigenetic data sets will provide useful tools to gain additional systems biology insights into cellular mechanism under a specific biological condition, for example, cancer cells or infected cells [5].

Many studies have focused on gene regulatory networks (GRNs) and protein–protein interaction networks (PPINs). Classical graph algorithms have been used to integrate cellular networks of protein–protein and protein–DNA interactions. Well-known simple and complex regulatory circuits have been identified and many putative regulatory circuits have been discussed. In [6], this study includes searches for composite network motifs, which consist of both transcription regulation and protein–protein interactions (PPIs) that recur significantly more often than expected in random networks. In [7], integrated networks comprise transcriptional and PPI data. Interspecific analysis has shown that several types of network motifs are not subject to any particular evolutionary pressure to preserve the corresponding interaction patterns of important biological functions. Studies have also integrated protein–protein and protein–DNA interactions to infer signaling regulatory pathways, but they focus only on pathways that explain gene expression changes in response to gene knockouts [8]. Recently, an integrated cellular network of transcription regulation and PPIs was introduced by Wang and Chen [9, 10]. They use gene-expression data at multiple time points to prune and combine candidate gene regulatory and signaling networks obtained from genome-scale data. An integrated and focused network for a specific condition of interest is then obtained. The transcriptional network is characterized as a dynamical system in which the expression of a target gene is modeled as a function of the regulatory effect of its corresponding transcription factors (TFs) and mRNA degradation. The modeling of a signaling/protein interaction network accounts for the activity of neighboring loci in the network. Genomic/transcriptomic and high-throughput methods have successfully identified many GRNs and PPINs. However, to explore cellular mechanism, we need more genomic data, such as epigenetic regulation data.

In real cellular systems, the expression of protein-coding genes is controlled by a complex network of genetic and epigenetic regulatory interactions. In addition to the abovementioned genetic regulation mechanisms, epigenetic regulation via DNA methylation is also an important regulatory mechanism in several cellular processes. Recently, high-throughput next-generation sequencing (NGS) data containing information on mRNAs, microRNAs (miRNAs), and methylation has been generated. Therefore, many compelling biological questions center on how regulation and interactions among genes, proteins, and epigenetic regulators give rise to specific cellular mechanisms. To address this problem, we proposed a method to construct an integrated cellular network that can explain specific cellular mechanisms under genetic and epigenetic regulation in response to specific biological conditions, based on the coupling of stochastic dynamic models. Recently, systems biology and computational biology methods have been widely employed to develop stochastic dynamic models that describe biological functions from a dynamic systems perspective [11–23]. Dynamic models to construct an integrated genetic and epigenetic cellular network (IGECN) not only provide a quantitative description of the integrated cellular network, but also predict the cellular mechanism of the network in response to various conditions, gene knockouts, treatments with external agents, etc [24].

In this study, we integrated omics data, including NGS [25], mRNA and miRNA expression [26], RNA sequencing (RNA-seq) [27], PPIs [28], transcription regulation interaction [29–32], miRNA-target gene association [33–37], and gene ontology (GO) data (http://geneontology.org/) to construct a candidate IGECN. A schematic diagram of the candidate IGECN is shown in Fig. 1. The candidate IGECN mainly consisted of three sub-networks. The first was the candidate PPIN, which included candidate PPIs in signal transduction pathways and metabolic pathways; the second was the GRN, which described transcription regulation; the third was the candidate miRNA regulatory network. Epigenetic DNA methylation was also considered, which involves the modification of DNA to influence mRNA transcription. The TFs in Fig. 1 are at the interface between the PPIN and GRN, and genes are at the interface between the GRN, PPIN, and miRNA and methylation regulatory networks.

Fig. 1
figure 1

Schematic diagram of candidate integrated genetic and epigenetic cellular networks (IGECN). The candidate IGECN was constructed using omics data and database mining. The candidate IGECN mainly consisted of three sub-networks: the candidate protein-protein interaction network (PPIN) including signal transduction pathways and metabolic pathways, the candidate gene regulatory network (GRN) in which transcriptional regulation occurs, and the candidate miRNA regulatory network. The epigenetic regulation of DNA methylation was considered to influence gene transcription. Transcription factors (TFs) represent the interface between PPINs and GRNs, and genes represent the interface between GRNs, miRNA regulatory networks, PPINs, and methylation regulatory networks

Since the candidate IGECN was constructed using omics data and big database mining, for data collected under various experimental conditions (including some that were gathered in vitro), it must contain many false positive interactions and regulations in the candidate PPIN and GRN, respectively. Therefore, it is necessary to prune these false positives, which may not explain real-world cellular mechanisms of biological organisms. In this situation, it is more appealing to evaluate these interactions and regulations in the candidate IGECN for investigating whole cellular mechanisms of a real organism under specific conditions (for example, HIV infections or cancer) using gene expression profile data, protein expression profile data, and miRNA expression profile data via dynamic interaction and regulation models. The false positive interactions and regulations can be pruned based on these real-time data to construct an accurate IGECN for some specific cellular mechanisms.

In this study, three coupling stochastic dynamic models were proposed to develop a candidate IGECN that describes the interplay between the GRN, PPI network, and miRNA network and accounts for methylation regulation in the GRN. After system parameter estimation for coupling stochastic dynamic models of candidate IGECN by these time-profile data, the Akaike information criterion (AIC) and the p-value statistical method [10] were employed to detect the system order (the number of interactions and regulations) of coupling stochastic dynamic models, and then prune the false positive regulations and interactions, which are not significant and out of system order, to obtain the real IGECN for some specific cellular mechanisms in response to some specific biological condition. By the interplay of genetic transcription regulation, PPIs, and miRNA regulation and DNA methylation of the identified IGECN, the cellular machanism can be elucidated comprehensively.

To investigate host-pathogen interspecies cellular interactions and regulations in the host cells is a very important research topic in attempts to understand pathogenic and defensive mechanisms in the HIV infection process and help control the clinic pathogenic infections. In this study, the IGECN is constructed for investigating host-pathogen crosstalks in HIV-infected cells by gene expression profile and miRNA profile of NGS data. Some significant cross-talk hubs in the interspecies IGECN can be extracted as significant crosstalk network marker, from which we could get more insight into the pathogenic mechanism of HIV virus and the defense mechanism of human. Further, the IGECN is also constructed for cancer cells and compared with normal cells to find the significant network marker to investigate some significant cues from the genetic and epigenetic regulatory mechanisms using NGS data.

Materials and methods for constructing IGECN

Data selection and processing in gastric cancer and liver cancer

Kim and colleagues [26] proposed microarray data, including gene expression data and miRNA expression data, 1 pair of data for normal and tumor tissues of stage I gastric cancer in humans, 10 pairs of data for normal and tumor tissues of stage II gastric cancer in humans, and 15 pairs of data for normal and tumor tissues of stage III gastric cancer in humans. According to The Cancer Genome Atlas (TCGA; https://tcga-data.nci.nih.gov/tcga/), the mRNA and miRNA expression data including 17 dataset for normal tissues and 163 dataset for tumor tissues of stage I liver cancer, 7 dataset for normal tissues and 82 dataset for tumor tissues of stage II liver cancer, 8 dataset for normal tissues and 80 dataset for tumor tissues of stage III liver cancer, and 1 dataset for normal tissues and 6 dataset for tumor tissues of stage IV liver cancer were available. The candidate associations for miRNA-gene regulatory associations are available at the TargetScan [38–41], TF-gene regulatory associations are available at the Human Transcriptional Regulation Interactions database (HTRIdb) [29] and Integrated Transcription Factor Platform (ITFP) [32], and PPIs can be found at The Biological General Repository for Interaction Datasets (BioGRID) [28] and available in the manually curated human signaling network in Edwin Wang Lab website (http://www.cancer-systemsbiology.org/) [42]. In order to integrate the above databases, we use gene symbols in NCBI human gene database to denote the standard names of genes in this study. We used Matlab’s text-file and string manipulation tools in text mining. The interaction candidates were integrated and unified from different databases.

Data selection and processing in CD4+ T cells of human with HIV-1 infection

Mohammadi and colleagues [25] reported gene expression data and miRNA expression data in HIV-1-infected cluster of differentiation 4+ (CD4+) T cells and non-infected (normal) cells to analyze host-pathogen interactions at early (or reverse transcription), middle (integration/replication), and late (virus assembly/budding) infection stages.

Construction of the dynamic models in IGECN

Our goal was to construct an IGECN in which the transcription regulation mechanisms, miRNA regulations, PPIs, and epigenetic regulation via DNA3 methylation are integrated for investigating whole cellular mechanisms under specific biological conditions. A flow chart describing the proposed method to construct the IGECN is shown in Fig. 2. Several kinds of NSG omics data and information derived from databases are integrated, including microarray gene expression data, PPI data, miRNA expression data, and miRNA-target gene association data, as the input for the proposed IGECN construction method. From these data, the candidate GRN, PPIN, and miRNA regulatory network using the broad TF-gene regulation pool, PPI pool, and miRNA regulatory pool, respectively, under various biological conditions were retrieved to construct a candidate IGECN. In this study, NGS data including gene expression and miRNA profiles were used to validate these integrated TF-gene regulation data, PPIs, and miRNA regulation data in the candidate IGECN. Therefore, stochastic dynamic models were used to identify interactions and regulatory mechanisms in the candidate IGCEN by gene expression and miRNA expression profiles using NGS data under a specific biological condition of interest, for example, cancer or infections.

Fig. 2
figure 2

Flowchart (a) and block diagram (b) of the proposed method to integrate genetic and epigenetic cellular networks

In the transcriptional regulation sub-network of the candidate IGECN, the candidate GRN was depicted as a stochastic dynamic system in which the regulations of TFs and miRNAs are the inputs and gene expression levels of target genes are outputs (Fig. 1). For target gene i in the candidate GRN, the stochastic dynamic model is described by the following stochastic discrete dynamic equation:

$$ {x}_i\left(t+1\right)={x}_i(t)+{\displaystyle \sum_{j=1}^{J_i}{a}_{ij}{y}_j(t)}-{\displaystyle \sum_{l=1}^{L_i}{c}_{il}{x}_i(t){s}_l(t)}-{\lambda}_i{x}_i(t)+{k}_i+{w}_i(t), $$
(1)

where x i(t + 1) denotes the mRNA expression level of the ith target gene at time t + 1, a ij represents the regulatory ability of the jth TF to the ith target gene (with a positive or negative value indicating activation or repression), the regulatory ability a ij is also influenced by the DNA methylations of the corresponding gene [43], y j (t) represents the regulation function of the jth TF binding to the target gene i, λ i indicates the degradation effect from the present time t to time t + 1, k i represents the basal level (which is also dependent on the methylation of gene i), s l (t) denotes the expression level of the lth miRNA, the parameter c il represents the rate of miRNA-mRNA coupling (and is dependent on the energy of RNA-RNA binding [44]), and w i (t) denotes the stochastic noise due to modeling residue and fluctuation in the expression profile of the target gene. The biological interpretation of equation (1) is that the mRNA expression x i (t + 1) of the target gene i at time t + 1 is a function of the present mRNA expression x i (t) and mRNA expression due to transcriptional regulation of J i TFs that bind to the target gene at time t, minus the inhibitory regulation of L i miRNAs at time t, minus the mRNA due to the degradation over time, plus the basal level of mRNA expression, and stochastic noise due to measurement and random fluctuations in the expression profile of target gene i.

In the candidate miRNA regulatory sub-network, the stochastic dynamic model of the lth miRNA can be described by the following stochastic discrete time dynamic equation [45]:

$$ {s}_l\left(t+1\right)={s}_l(t)-{\displaystyle \sum_{i=1}^{I_l}{c}_{il}{s}_l(t){x}_i(t)}-{r}_l{s}_l(t)+{k}_l+{\varpi}_l(t), $$
(2)

where the degradation of miRNA, i.e. \( {\displaystyle {\sum}_{i=1}^{I_l}{c}_{il}{s}_l(t){x}_i(t)} \), is accompanied by mRNA degradation events induced by the miRNA in (1). That is, the lth miRNA induces the degradation of I l mRNAs at present. r l denotes the self-degradation rate, k l denotes the production rate of the lth miRNA, and Ï– l (t) denotes the stochastic noise.

In the PPI sub-network, the candidate PPIs were depicted as a system in which expressions of mRNA and TFs are input/output of the system (Fig. 1). For a target protein j in the candidate PPI sub-network, the stochastic dynamic model of protein expression level is given as follows:

$$ {y}_j\left(t+1\right)={y}_j(t)-{\displaystyle \sum_{m=1}^{M_j}{b}_{jm}{y}_j(t){y}_m(t)}+{\alpha}_j{x}_j(t)-{\beta}_j{y}_j(t)+{h}_j+{\zeta}_j(t) $$
(3)

where y j (t) represents the expression level at time t of the target protein j, b jm denotes the interaction ability between the mth interactive protein y m (t) and the target protein y j (t), α j represents the effect of protein synthesis from mRNA, x j (t) denotes the mRNA expression level of the jth protein, h j represents the basal expression level, and ζ j (t) is the stochastic noise. The rate of formation of the protein complex y j (t)y m (t) is proportional to the product of the concentration of each proteins [46, 47], i.e., it is proportional to the probability of molecular collisions between two proteins; thus, the protein interaction complex was modeled as nonlinear multiplication scheme.

The biological interpretation of equation (3) is that the expression level of the target protein j at time t + 1 is a function of the protein expression level at time t, minus protein complex interactions with M j proteins, plus the translation effect, minus the present protein degradation, plus the basal protein level from other sources and stochastic interactions. It is noted that there is no directionality for interacting proteins in the PPI sub-network.

The regulatory effects and interactions between genes, proteins, miRNA, and DNA methylation in the candidate IGECN are described in the following. Some TFs y j (t) at the end of candidate signal transduction pathways (PPINs) regulate target genes according to the regulation function a ij y j (t) shown in equation (1), and mRNAs of target genes are also negatively regulated by some miRNAs s l (t). The regulated genes influence the corresponding protein expression level via translation from mRNA x j (t) to protein α j x j (t), as described in (3). An mRNA degradation event, -c il x i (t)s l (t), described in (1), induced by miRNA s l (t), also leads to the degradation of miRNA, -c il s l (t)x i (t), in (2). Epigenetic DNA methylation also influences the regulation parameter a ij and basal level k i in (1) and the translation rate α j from mRNA to protein in (3). The interplay among genes, proteins, miRNAs, and DNA methylation can be seen in (1)- (3) and Fig. 1, in which TFs are at the interface between PPINs and GRNs, and genes are at the interface between GRNs, miRNA, PPINs, and DNA methylation.

Based on the above stochastic dynamic models (1)- (3), the candidate GRN can be linked through the regulatory parameters a ij in (1) between genes and their possible regulatory TFs (Fig. 1), and the candidate PPIN can be linked through the interaction parameters b jm in (3) between possible interacting proteins. The candidate GRN and miRNA regulatory network are linked through the miRNA-mRNA coupling parameters c il in (1) and (2). Since the candidate IGECN only indicated potential TF-gene regulatory effects, PPIs, and miRNA-mRNA couplings based on data collected by database mining, they should be confirmed by gene expression profile, protein expression profile, and miRNA expression profile in NGS data. The regulatory effects a ij and c il in (1), interaction b jm in (3) and regulation c il in (2) , which represent the edges of the candidate IGECN, should be identified and validated by empirical gene, protein, and miRNA expression data. These regulatory and interactive parameters were evaluated using temporal data by solving the constrained least square parameter estimation problem. Owing to the limited number of temporal data points, to avoid overfitting when identifying the parameters in (1–3), the cubic spline method [10] was also used to interpolate extra time points for gene expression data. The gene regulatory parameters a ij and c il were identified gene by gene (and PPI parameters b jm protein by protein), so that the candidate integrated network identification process was not limited by the number of genes, proteins, and miRNAs in the candidate IGECN. In other words, the regulatory parameters a ij and c il in (1) were first identified for target gene i and then for target gene i + 1, i + 2, etc. The models in (1–3) only represent a small fraction of real biology. For example, many post-transcriptional regulators (RNA-binding proteins) as well as epigenetic regulations beyond methylation were involved. The model uncertainty of real biological system and the fluctuation of expression data were involved in the stochastic noises, w i (t) in (1), ϖ l (t) in (2), and ζ j (t) in (3). A model order detection method such as AIC and statistical assessment like the Student’s t-test [10] were used to prune the false-positive regulatory effects and interactions and then detect the true interaction and regulation number of each in the candidate IGECN. AIC [10, 48, 49] was applied to detect the real regulations and interactions by pruning the false positive regulations and interactions in the candidate GRN (1), candidate miRNA regulatory network (2), and the candidate PPI network (3). AIC, which includes both estimated residual error and model complexity in one statistics, quantifies the relative goodness of fit of a model The minimization achieved in AIC will indicate the real network with the true model order. By applying AIC detection method, the insignificant parameters, a ij , c il and b jm , out of network order are considered as false positives and pruned from candidate IGECNs. We then used Student’s t-test to evaluate the significance of each regulatory parameter in (1–3) [10]. In this way, the candidate GRN, candidate PPIN, and candidate miRNA regulatory network based on large-scale data mining were pruned using the temporal data, leading to the construction of an accurate GRN, PPIN, and miRNA regulatory network. Based on the interactive effects between GRNs and signal transduction pathways (or PPINs) via TFs, and among GRNs, miRNA regulatory networks, and PPINs via coupling genes, these three networks were coupled to consist of the IGECN.

We proposed a general method to construct the IGECN using large-scale database mining and system identification for cellular functions. In the following, two examples are given to illustrate the application of the IGECN system to investigate cellular mechanisms under different biological conditions. In the first example, an IGECN was constructed to investigate cellular mechanisms under cancer conditions and in the second example, an IGECN was constructed to investigate pathogenic and defensive cellular mechanisms under HIV infection conditions.

Identification of the regulatory parameter aij from the ith TF to the jth target gene, interaction bjm between the jth and mth proteins, and coupling rate cil between the lth miRNA and the ith mRNA in IGECN

After constructing the coupling dynamic models (1)-(3) in the candidate IGECN, the regulative and interactive parameters have to be identified using NGS data. The strategy is to identify the IGECN gene by gene, protein by protein, and miRNA by miRNA. Before applying the identification method, we first examine the dynamic models carefully. In (1), the basal level k i should be always non-negative, because the gene expressions are always non-negative. Because the parameters in (1) have certain constraints, the regulatory parameters were identified by solving the constrained least square parameter estimation.

The mRNA model (1) was rewritten as the following linear regression form

$$ \begin{array}{l}{x}_i\left(t+1\right)=\left[\begin{array}{cccccccc}\hfill {y}_1(t)\hfill & \hfill \cdots \hfill & \hfill {y}_{J_i}(t)\hfill & \hfill {x}_i(t)\hfill & \hfill -{s}_1(t){x}_i(t)\hfill & \hfill \cdots \hfill & \hfill -{s}_{L_i}(t){x}_i(t)\hfill & \hfill 1\hfill \end{array}\right]\left[\begin{array}{c}\hfill {a}_{i1}\hfill \\ {}\hfill \vdots \hfill \\ {}\hfill {a}_{i{J}_i}\hfill \\ {}\hfill \left(1-{\lambda}_i\right)\hfill \\ {}\hfill {c}_{i1}\hfill \\ {}\hfill \vdots \hfill \\ {}\hfill {c}_{i{L}_i}\hfill \\ {}\hfill {k}_i\hfill \end{array}\right]+{w}_i(t)\\ {}\kern3em ={\phi}_i(t){\theta}_i^1+{w}_i(t)\end{array} $$
(4)

where ϕ i (t) denotes the regression vector which can be obtained from the processing above and θ 1 i is the parameter vector of target gene i to be estimated. In order to avoid overfitting in the parameter estimation process, when identifying the regulatory parameters, the cubic spline method [50–52] was also used to interpolate extra time points for y i (t), x i (t) and s i (t) from NGS data.

The mRNA model (4) at different time points can be arranged as follows.

$$ \left[\begin{array}{c}\hfill {x}_i\left({t}_2\right)\hfill \\ {}\hfill {x}_i\left({t}_3\right)\hfill \\ {}\hfill \vdots \hfill \\ {}\hfill {x}_i\left({t}_L\right)\hfill \end{array}\right]=\left[\begin{array}{c}\hfill {\phi}_i\left({t}_1\right)\hfill \\ {}\hfill {\phi}_i\left({t}_2\right)\hfill \\ {}\hfill \vdots \hfill \\ {}\hfill {\phi}_i\left({t}_{L-1}\right)\hfill \end{array}\right]{\theta}_i^1+\left[\begin{array}{c}\hfill {w}_i\left({t}_1\right)\hfill \\ {}\hfill {w}_i\left({t}_3\right)\hfill \\ {}\hfill \vdots \hfill \\ {}\hfill {w}_i\left({t}_{L-1}\right)\hfill \end{array}\right] $$
(5)

where L is the number of expression time points of NGS data after cubic spline interpolation.

For simplicity, we define the symbols X i , Φ i, and W i to express (5) as follows.

$$ {X}_i={\Phi}_i{\theta}_i^1+{W}_i $$
(6)

The constrained least square parameter estimation problem of θ 1 i is showed as follows.

$$ \underset{\theta_i^1}{ \min }{\left\Vert {\Phi}_i{\theta}_i^1-{X}_i\right\Vert}_2^2\kern0.75em \mathrm{subject}\ \mathrm{t}\mathrm{o}\kern0.5em A{\theta}_i^1\le b $$
(7)

where and b = [0…0]T that give the constraints to force the miRNA inhibition -c il to be always non-positive, and the basal level k i to be always non-negative in (1), i.e. -c il  ≤ 0, and k i  ≥ 0.The constrained least square problem was solved using the active set method for quadratic programming [53, 54].

Similarly, the miRNA model (2) was rewritten in the following regression form.

$$ \begin{array}{l}{s}_l\left(t+1\right)=\left[\begin{array}{ccccc}\hfill {s}_l(t)\hfill & \hfill -{s}_l(t){x}_1(t)\hfill & \hfill \cdots \hfill & \hfill -{s}_l(t){x}_{I_l}(t)\hfill & \hfill 1\hfill \end{array}\right]\left[\begin{array}{c}\hfill 1-{r}_l\hfill \\ {}\hfill {c}_{1l}\hfill \\ {}\hfill \vdots \hfill \\ {}\hfill {c}_{I_ll}\hfill \\ {}\hfill {k}_l\hfill \end{array}\right]+{\varpi}_l(t)\\ {}\kern2.75em ={\psi}_l(t){\theta}_l^2+{\varpi}_l(t)\end{array} $$
(8)

where ψ l (t) denotes the regression vector and θ 2 l indicates the parameter vector to be estimated. By cubic spline method, at different time points, (8) can be represented as the following equation.

$$ {S}_l={\Psi}_l{\theta}_l^2+{W}_l $$
(9)

The parameter identification problem is then expressed as follows.

$$ \underset{\theta_i^2}{ \min }{\left\Vert {\Psi}_l{\theta}_l^2-{S}_l\right\Vert}_2^2\kern0.5em \mathrm{subject}\ \mathrm{t}\mathrm{o}\kern0.5em B{\theta}_l^1\le b $$
(10)

where \( B=\begin{array}{cc}\hfill \left[\begin{array}{c}\hfill 0\hfill \\ {}\hfill 0\hfill \\ {}\hfill 0\hfill \end{array}\right.\hfill & \hfill \overset{I_l+1}{\overbrace{\left.\begin{array}{ccc}\hfill -1\hfill & \hfill 0\hfill & \hfill 0\hfill \\ {}\hfill 0\hfill & \hfill \ddots \hfill & \hfill 0\hfill \\ {}\hfill 0\hfill & \hfill 0\hfill & \hfill -1\hfill \end{array}\right]}}\hfill \end{array} \), and b = [0…0]T that give the constraint to force the miRNA degradation -c il to be always non-positive, and the basal level k l to be always non-negative in (2), i.e. -c il  ≤ 0, and k l  ≥ 0. Finally, the protein model (3) uses the same way like above to make sure –b jm  ≤ 0, and α j , h j  ≥ 0. According to large-scale measurement of protein activities, 73 % of the variance in protein abundance could be illustrated by mRNA abundance [55]. mRNA expression levels were frequently used to substitute for the protein expression levels. If the simultaneously measured genome-wide protein expression data and mRNA expression data in each cancer stage are available, the general models in (1–3) can also be applied to identify the real network in cancer. After the parameter identification problem was solved, we can identify the IGECN. For example, we identified the regulatory parameter a RBM5,ETS1 =-0.22 from the TF ETS1 to the target gene RBM5, interaction parameter b PGR,NCOA2 =0.06 between the 2 proteins PGR and NCOA2, and coupling rate c USP15,MIR381 =0.001 between the miRNA MIR381 and the mRNA USP15 in gastric cancer cells.

By using Student’s t-test [56], the p-values for the estimated parameters, including the regulatory ability a ij in (1), the miRNA-mRNA coupling rate c il in (1) and (2), the protein interaction ability b jm in (3), were calculated to determine the significance of the parameters. By using one-way ANOVA, we calculated the p-value for each gene to determine the significance of expression change between normal and tumor cells in cancers, and between normal and HIV-infected human cells. The human DNA methylation profiles of liver cancers (including 41 normal and 369 tumor cells), and gastric cancer (including 2 normal and 396 tumor cells) were available in TCGA. A One-Way ANOVA was also used to calculate the p-value of each gene to determine the significance of DNA methylation change between normal and tumor cells in each cancer.

Results

The construction of IGECN is very useful to investigate the whole cellular mechanisms of specific biological conditions by big database mining and NGS data. The following two examples were given to illustrate how to construct IGECN to investigate the roles of genetic and epigenetic methylation and miRNA dysregulations in gastric tumorigenesis and HIV infection. Based on these significant core network biomarkers in IGECN, potential multiple drug targets were also proposed for the development of therapeutic treatment.

  1. (I)

    Construction of IGECN to investigate significant cellular mechanisms in gastric tumorigenesis and hepatocarcinogenesis

Cancer is a complex and heterogeneous disease that is highly robust against therapeutic interventions. It is a constellation of diverse and evolving disorders characterized by the uncontrolled proliferation of cells that may eventually lead to fatal dysfunction of the cellular system. Although some cancer subtypes can be cured by early diagnosis and specific treatment, no effective treatment has been established for many subtypes. Due to the complex, heterogeneous, and evolving nature of cancer, it is appealing to construct an IGECN to gain insight into the carcinogenesis process using large-scale data mining and systematic genetic modeling and various types of omics and temporal data.

Based on the importance of epigenetic modifications in cancer, we identified IGECNs using raw data collected at various stages of cancer as the temporal information for gastric tumorigenesis (Fig. 3) and hepatocarcinogenesis (Fig. 4), demonstrating large differences between miRNA regulations of normal and tumor tissues in gastric tumorigenesis. The results (Fig. 3) show that two cell cycle genes, ubiquitin C (UBC) and ataxia telangiectasia mutated (ATM), and one chromosomal translocation gene, CREB binding protein (CREBBP), play important roles in gastric tumorigenesis via genetic and epigenetic regulation. ATM is an important cell cycle checkpoint kinase [57] , and chromosomal translocation plays a critical role in gastric tumorigenesis [58]. We identified that 7 genes have basal level difference between cancer and normal cells, including RUNX1 (P-value < 1.4☓10-3), ARHGDIB (P-value < 1.6☓10-2), PSME1 (P-value < 4☓10-2), VAMP5 (P-value < 2☓10-2), POLG (P-value < 1.3☓10-2), GZMB (P-value < 1.44☓10-2), and RBM5 (P-value < 8.2☓10-3), which are mainly due to DNA methylation of the corresponding genes. In order to validate our finding, we used the epigenomics data of human gastric cancer in TCGA. According to DNA methylation profiles of gastric cancer cells and normal gastric cells from TCGA, we found that 5 genes have methylation change, including RUNX1 (P-value < 3.56☓10-1), ARHGDIB (P-value < 1.81☓10-1), PSME1 (P-value < 2.59☓10-1), GZMB (P-value < 1.79☓10-1), and RBM5 (P-value < 3.94☓10-1), and 2 genes have small methylation change including VAMP5 (P-value < 8.61☓10-1), POLG (P-value < 8.38☓10-1) between gastric cancer cells and normal gastric cells (Fig. 3).

Fig. 3
figure 3

The integrated genetic and epigenetic cellular network for investigating the cellular mechanisms of normal and tumor tissues in gastric cancer. The associations include the networks in tumor cells (dotted lines), normal cells (dashed lines), and both (solid lines). The bold lines indicate large genetic and epigenetic regulatory effects or genetic expression. Dysregulation of MIR7 contributes to initiation and progression of inflammation-induced gastric cancer, and PML, a potential drug target, overcomes platinum resistance in platinum-based chemotherapy in gastric cancer

Fig. 4
figure 4

The integrated genetic and epigenetic cellular network for investigating the cellular mechanisms of normal and tumor tissues in liver cancer. The associations include the networks in tumor cells (dotted lines), normal cells (dashed lines), and both (solid lines). The bold lines indicate large genetic and epigenetic regulatory effects or genetic expression. Dysregulations of MIR141 and MIR335 contributes to dysfunction of DNA damage response and cell proliferation in human HCC, respectively

Moreover, we also identified that 5 genes have basal level difference between Hepatocellular Carcinomas (HCCs) and normal cells, including ELAVL1 (P-value < 3.69☓10-5), JMJD1C (P-value < 3.59☓10-11), GADD45A (P-value < 7.32☓10-11), EIF2AK3 (P-value < 5☓10-5), and GSN (P-value < 1.19☓10-2), which is mainly due to DNA methylation of the corresponding genes (Fig. 4). We used the epigenomics data of HCCs in TCGA to support the results. According to DNA methylation profiles of HCCs and normal liver cells from TCGA, we found that 5 gene have methylation changes, including ELAVL1 (P-value < 5.7☓10-2), JMJD1C (P-value < 1.4☓10-2), GADD45A (P-value < 9.6☓10-2), EIF2AK3 (P-value < 1.85☓10-6), and GSN (P-value < 8.71☓10-4).

  1. (II)

    Construction of an IGECN for investigating pathogenic and defensive mechanisms of HIV-1-infected cells

Human HIV was first identified in 1983 and rapidly emerged as a virus responsible for a pandemic. At the beginning of this century, the HIV epidemic received global attention, and methods to halt and reverse the AIDS epidemic became an important United National Millennium topic. Currently, approximately 35 million people are living with HIV worldwide. With the advent of highly active antiretrovial therapies, HIV now can be managed as a chronic disease, but developing a cure remains a significant challenge. Therefore, the investigation of interspecific host-pathogen regulatory mechanisms and interactions in cellular genetic and epigenetic systems is a very important research topic to improve understanding of pathogenic and defensive mechanisms in the HIV infection process and to help control clinical pathogenic infections. In this section, using two-sided time-profile HIV-human high-throughput sequence data (NGS), RT-PCR data, miRNA data, and other omic data, the interspecific host-pathogen IGECN was constructed for the HIV infection process. To construct this model, the coupling dynamic equations (1)-(3) should be modified to consider the pathogen in addition to the host species as follows:

$$ \begin{array}{c}\hfill {x}_i\left(t+1\right)={x}_i(t)+{\displaystyle \sum_{j=1}^{J_i}{a}_{ij}{y}_j(t)}-{\displaystyle \sum_{l=1}^{L_i}{c}_{il}{x}_i(t){s}_l(t)}+{\displaystyle \sum_{o=1}^{O_i}{d}_{io}{h}_o(t)}-{\lambda}_i{x}_i(t)+{k}_i+{w}_i(t),\hfill \\ {}\hfill {y}_j\left(t+1\right)={y}_j(t)-{\displaystyle \sum_{m=1}^{M_j}{b}_{jm}{y}_j(t){y}_m(t)}-{\displaystyle \sum_{o=1}^{O_j}{g}_{jo}{y}_j(t){h}_o(t)}+{\alpha}_j{x}_j(t)-{\beta}_j{y}_j(t)+{h}_j+{\zeta}_j(t),\hfill \\ {}\hfill {h}_o\left(t+1\right)=h{}_o(t)-{\displaystyle \sum_{j=1}^{J_o}{g}_{jo}{h}_o(t){y}_j(t)}-{\displaystyle \sum_{l=1}^{L_o}{m}_{ol}{h}_o(t){s}_l(t)}-{\kappa}_o{h}_o(t)+{\delta}_o+{\xi}_o(t),\hfill \\ {}\hfill {s}_l\left(t+1\right)={s}_l(t)-{\displaystyle \sum_{i=1}^{I_l}{c}_{il}{s}_l(t){x}_i(t)}-{\displaystyle \sum_{o=1}^{O_l}{m}_{ol}{s}_l(t){h}_o(t)}-{r}_l{s}_l(t)+{k}_l+{\varpi}_l(t),\hfill \end{array} $$
(11)

where h o (t) denotes expression of the oth pathogen protein, and the crosstalk effect of pathogen proteins h o (t) are considered. After system identification of coupling IGECN between host and pathogen in (11) by NGS data at different infection stages, we found that 3 genes had differences in their basal levels between cancer and normal cells at the early infection stage, TAF5 (P-value < 2.4☓10-1), DDX3X (P-value < 2☓10-1), and CELF1 (P-value < 9.47☓10-2); 2 genes had basal level differences between cancer and normal cells at the middle stage, TAF5 (P-value < 5.8☓10-2), and CELF1 (P-value < 2.07☓10-7); 2 genes had basal level differences between cancer and normal cells at the late stage, ZNF451 (P-value < 9☓10-2), and CELF1 (P-value < 1.6☓10-3). These genes with different basal levels at different infection stages are mainly due to DNA methylation of stage-specific genes. Previous evidence for DNA methylation, including TAF5 [59] and DDX3X [60] methylation in individuals with HIV infection, ZNF451 methylation in individuals with Epstein–Barr virus (EBV) infection [61], and CELF1 methylation in individuals with Moloney leukemia virus infection [62], supports our findings. However, DNA methylation of genes are highly tissue-specific. If the methylation profiles of genes in HIV-infected and mock cells are available in the future, the results can be specifically supported.

By applying the raw expression data and the candidate associations in cancer cells to our model (4), we identified the IGECNs at the three infection stages (Figs. 5, 6 and 7). In the following section, we discuss the contribution of miRNA dysregulation to gastric cancer, HIV-1 infection progress, and liver cancer with respect to biological functions. The results have applications to the development of therapeutic treatments.

Fig. 5
figure 5

The integrated genetic and epigenetic cellular network for investigating the cellular mechanisms of normal and HIV-infected T cells at the first HIV-1 infection stage, i.e., the reverse transcription stage. The associations include those in HIV-infected cells (dotted lines), normal cells (dashed lines), and both cell types (solid lines). The bold lines indicate large genetic and epigenetic regulatory effects or genetic expression. Dysregulation of the signaling cascade from MIR9 to SUMO3 contributes to HIV-1 infection, hijacking CD4+ T cells via dysfunction of the immune and hormone pathways at the early infection stage

Fig. 6
figure 6

The integrated genetic and epigenetic cellular network for investigating the cellular mechanisms of normal and HIV-infected T cells at the intermediate HIV infection stage, i.e., the integration/replication stage. The associations include those in HIV-infected cells (dotted lines), normal cells (dashed lines), and both cell types (solid lines). The bold lines indicate large genetic and epigenetic regulatory effects or genetic expression. Dysregulation of MIR139-5p, and MIRLET7i contributes to dysfunction of the immune response and viral replication via DNA hypermethylation

Fig. 7
figure 7

The integrated genetic and epigenetic cellular network for investigating the cellular mechanisms of normal and HIV-infected T cells at the late HIV infection stage, i.e., the virus assembly/budding stage. The associations include those in HIV-infected cells (dotted lines), normal cells (dashed lines), and both cell types (solid lines). The bold lines indicate large genetic and epigenetic regulatory effects or genetic expression. Dysregulation of MIR101, MIR141, and MIR152 to the HIV-1 Gag protein contributes to HIV-1 budding and release via DNA hypermethylation, ubiquitin transfer, and endoplasmic reticulum-associated degradation at the late infection stage

Briefly, dysregulation of MIR7 contributes to the initiation and progression of inflammation-induced gastric cancer; dysregulation of MIR9 contributes to HIV-1 infection to hijack CD4+ T cells through dysfunction of the immune and hormone pathways; dysregulation of MIR139-5p, MIRLET7i, and MIR10a contributes to the HIV-1 integration/replication stage through DNA hypermethylation and immune system dysfunction; dysregulation of MIR101, MIR141, and MIR152 contributes to the HIV-1 virus assembly/budding stage through DNA hypermethylation, ubiquitin transfer, and endoplasmic reticulum-associated degradation; dysregulation of MIR302a contributes to not only microvesicle-mediated transfer of miRNAs but also dysfunction of NF-κB signaling pathway in hepatocarcinogenesis.

Discussion

  1. (a)

    Gastric Cancer

MIR7 contributes to the initiation and progression of inflammation-induced gastric cancer

Until now, it has been suggested that miRNAs play important roles in regulating proliferation and apoptosis to promote the invasion, metastasis, and angiogenesis of gastric cancer cells [63–65]. Although a few connections between genetic-and-epigenetic regulatory mechanisms in gastric cancer cells can be considered independently using biological experiments [63–66], the precise biological functions of miRNAs in gastric cancer cells and miRNAs affecting gastric cancer progression need to be determined by comparing the IGECN with that of normal cells.

In this study, we identified the IGECN including the largest differences in miRNA regulation between normal and cancer cells (Fig. 3). We identified that MIR7 represses SMG1 (p-value < 10-16) in cancer cells. It has been reported that both MIR7 and SMG1 are involved in the initiation and progression of inflammation-induced gastric cancer [67, 68]. The significant expression changes of MIR7 (p-value < 0.151) and SMG1 (p-value < 0.009) between normal gastric cells and gastric cancer cells also supported the result that MIR7 and SMG1 play important role in Human gastric carcinogenesis.

We identified that the expression change of MIR381 (p-value < 0.83) contributes to the expression changes of RNA-binding protein 5 (RBM5) (p-value < 0.008) and ubiquitin specific peptidase 15 (USP15) (p-value < 6☓10-4) via significant miRNA regulations (p-value < 10-16) in cancer cells. The activities of RBM5 and USP15 were also significantly activated in gastric cancer. It has been reported that the reduced RBM5 and USP15 expressions result in tumor cell proliferation [69, 70] and induce tumor cell apoptosis [71], respectively. It has also been suggested that dysregulation of two proteins, IL4 and CD40, leads to dysfunction of cell proliferation in mantle cell lymphoma [72], and dysregulation of MIR381 leads to dysfunction of proliferation in esophageal cancer [73]. Thus, we proposed that dysregulation of miRNAs in gastric cancer contributes substantially to the dysfunction of proliferation and apoptosis in gastric cancer via the ubiquitin-proteasome system [74]. TCGA offers a multilayered view of genomics and epigenomics data of many human cancer types. According to DNA methylation profiles of gastric cancer cells and normal gastric cells from TCGA, the methylation changes of GZMB (p-value < 0.18), ARHGDIB (p-value < 0.18), PSME1 (p-value < 0.26), RUNX1 (p-value < 0.36), RBM5 (p-value < 0.39), POLG (p-value < 0.84), and VAMP5 (p-value < 0.86) were found by the IGECN. Although the methylation changes of POLG and VAMP5 were insignificant, the gene expression changes of POLG (p-value < 0.01) and VAMP5 (p-value < 0.02) were still significant. Therefore, we suggested that other regulators may exist in regulating the expressions of POLG and VAMP5 in gastric cancer.

PML, a potential drug target, overcomes platinum resistance in platinum-based chemotherapy in gastric cancer

Although platinum-based chemotherapy has been used to treat most malignant diseases, platinum resistance still limits its efficacy. The nuclear receptor coactivator family (NCOAs) and the p160 steroid receptor coactivator family (SRCs) in gastric cancer cells exhibit significant differences in interactions observed in the IGECN than in normal cells. It has been proposed that SRC3 plays an important role in intrinsic resistance to platinum therapy in cancers [75]. The results revealed that dysregulation of NCOA2 is transduced by signaling cascades via promyelocytic leukemia (PML) and RUNX1, which is transcriptionally regulated by v-ets erythroblastosis virus E26 oncogene homolog 1 (ETS1). It has been suggested that ETS1 plays an important role in regulating gastric cancer cell invasion, metastasis, and angiogenesis [66, 76]. Because RUNX1 is essential for the development of human hematopoietic stem cells [77], we proposed that PML is a potential drug target that can overcome platinum resistance in the platinum-based chemotherapy in gastric cancer. It has also been suggested that PML mediates glioblastoma resistance to the mammalian target in rapamycin-targeted therapies [78]. Two pharmaceutical drugs, arsenic trioxide and retinoin, potentially target PML to overcome platinum resistance in gastric cancer.

  1. (b)

    Hepatocellular Carcinoma (HCC)

Dysregulations of DNA methylation and MIR141 contribute to dysfunction of DNA damage response in HCC

HCC frequently occurs in the context of chronic disease that promotes DNA damage and chromosomal aberrations [79]. Recently, it has been found that the overexpressed MIR141 results in extensive DNA damage [80], and JMJD1C and TP73 are the components of DNA-damage response pathway with implications for tumorigenesis [81, 82]. The accumulated alterations in the DNA-damage response pathway could trigger hepatocarcinogenesis [83, 84]. In this study, we identified that the 3 genetic and epigenetic regulations, including the gain of the miRNA regulation of MIR141 to JMJD1C (p-value < 0.789) and the genetic regulation of TP73 to GADD45A (p-value < 0.037) in cancer cells, and DNA methylations of JMJD1C and GADD45A, led to not only changes of gene expression profiles of JMJD1C (p-value < 7.3☓10-11) and GADD45A (p-value < 7.3x☓10-11) but also changes of DNA methylation profiles of JMJD1C (p-value < 0.014) and GADD45A (p-value < 0.096) between human normal liver cells and liver cancer cells (Fig. 4). The changes also resulted in the gain of PPI between JMJD1C and GADD45A (p-value < 0.48) in hepatocarcinogenesis. Therefore, we suggested that the dysregulations of DNA methylation and MIR141 give rise to dysfunction of DNA damage response in HCC.

Dysregulation of MIR335 contributes human HCC cell proliferation

Recently, it has been proposed that down-regulation of PTPMT1 is sufficient to trigger cancer cell death [85], and TRIM27 participates in regulating cell proliferation in the development of cancer [86]. We identified that the change in the miRNA expression profile of MIR335 (p-value < 0.808) results in the significant change in the gene expression profile of PTPMT1 (p-value < 8.4☓10-4) between normal liver cells and liver cancer cells (Fig. 4). The expressions of TRIM27 and PTPMT1 were significantly activated in HCC. The significant genetic alterations of TRIM27 (p-value < 8.4☓10-4) and PTPMT1, and their interaction (p-value < 0.32) in HCC promote aberrant proliferation. Therefore, we suggested that dysregulation of MIR335 contributes human HCC cell proliferation.

Dysregulation of MIR302a contributes to not only microvesicle-mediated transfer of miRNAs but also dysfunction of NF-κB signaling pathway in hepatocarcinogenesis

It has been identified that microvesicles transfer miRNA between cells and alter biological functions by affecting signaling pathways in hepatocarcinogenesis [87]. Microvesicles are enriched with maladjusted miRNAs, which could regulate zinc finger proteins. We identified that the change in the miRNA expression profile of MIR302a (p-value < 0.734) leads to the change in the gene expression profile of teashirt zinc finger homeobox 3 (TSHZ3) (p-value < 0.189), also known as zinc finger protein 537 (ZNF537), between normal liver cells and liver cancer cells (Fig. 4). The expression change of TSHZ3 further contributes to the change in the gene expression profile of TNF receptor-associated factor 2 (TRAF2) (p-value < 1.73☓10-11), which mediates the dysfunction of the nuclear factor-kappa B (NF-κB) signaling pathway in hepatocarcinogenesis. The dysregulated NF-κB signaling pathway results in HCC angiogenesis and metastasis [88]. Suppression of TRAF2 inhibits the proliferation of several cancer cells [89], and the oncogene in epithelial cancers, TRAF2, contributing to the dysregulated NF-κB signaling pathway has been also identified [89]. Therefore, dysregulation of MIR302a contributes to not only microvesicle-mediated transfer of miRNAs but also dysfunction of NF-κB signaling pathway in hepatocarcinogenesis.

Dysregulation of the hormone receptors, ESR2 and AR, contributes to hepatocarcinogenesis

The most recent study proposed that menopause, which results from growing estrogen deficiency and physiological aging, increases risk of immune system disorders, mitochondrial dysfunction, cellular senescence, and imbalance between antioxidant formation and oxidative stress [90]. The physiologic and biochemical changes in those functions have a direct effect on liver function and mediate the development of liver disease. In this study, we identified two hormone receptors, estrogen receptor 2 (ESR2) and androgen receptor (AR), which were involved in multiple changes in intracellular signaling cascades. The gene expression profile of ESR2 showed the change (p-value < 0.091) between young and old normal liver cells, and the change (p-value < 0.063) between young normal liver cells and young liver cancer cells (Fig. 4). The relatively small change in the gene expression profile of ESR2 (p-value < 0.45) between old normal liver cells and old liver cancer cells implicates that the gene expression change of ESR2 could be necessary but not sufficient for hepatocarcinogenesis. The significant change in the gene expression profile of AR not only between young normal liver cells and young liver cancer cells (p-value < 2.46☓10-6) but also between old normal liver cells and old liver cancer cells (p-value < 6.13☓10-7) implicates that the gene expression change of AR may be required for hepatocarcinogenesis. The result can be supported by the most recent study in hepatocarcinogenesis [91].

  1. (c)

    HIV-1 Infection

The signaling cascade from MIR9 to SUMO3 contributes to HIV-1 infection to hijack CD4+ T cells through dysfunction of the immune and hormone pathways at early stage

Owing to the available HIV-1 mRNA gene expression data at middle and late infection stages and the lack of HIV-1 protein expression data, the impact of HIV-1 proteins on the host cells cannot be identified, directly. In this study, we identified the IGECNs at early, middle, and late infection stages (Figs. 5, 6 and 7), including the largest differences in the miRNA regulation and HIV-1 protein interactions at the middle and late stages between normal and infected cells. At early infection stage (Fig. 5), we identified that the expression changes of MIR9 (p-value < 0.448) and general transcription factor IIi (GTF2I) (p-value < 0.96) contribute to the expression change of small ubiquitin-like modifier 3 (SUMO3) (p-value < 0.016) via respectively regulating zinc finger protein 131 (ZNF131) (p-value < 9.3☓10-14) and DEAD (Asp-Glu-Ala-Asp) box helicase 3, X-linked (DDX3X) (p-value < 1☓10-16). By comparing miRNA expression profiles between the wildtype HIV-1 infected human CEMx174 lymphocytes and Tat RNA silencing suppressor in HIV-1 infected human CEMx174 lymphocytes [92], Tat RNA silencing suppressor contributed to the expression change of MIR9 (p-value < 0.035). It has been also proposed that Tat can induce MIR9 to control inflammatory responses [93], and Tat can also control GTF2I expression during HIV infection [94]. Additionally, it has been also reported that SUMO modification of ZNF131 suppressed estrogen signaling [95], and HIV transcription was suppressed by inducing estrogen signaling [96]. Therefore, we suggested that the HIV infection leads to dysregulation of immune and hormone pathways to hijack CD4+ T cells.

Dysregulation of MIR34c-5p and MIR449a contributes to the frequency of latent HIV infection

Moreover, we also identified that the expression changes of MIR34c-5p (p-value < 0.33), and MIR449a (p-value < 0.50) contributed to the expression changes of TATA binding protein-associated factor 5 (TAF5) (p-value < 0.25) via miRNA regulations (p-value < 1☓10-16). It has been suggested that the genetic and epigenetic regulations in TATA protein binding may control the frequency of latent HIV infection [97]. Therefore, we suggested that genetic and epigenetic regulations in TAF5 also contribute to control HIV infection at early stage.

Dysregulation of MIR139-5p contributes to HIV-1 replication via DNA hypermethylation

At the intermediate HIV infection stage, we identified that the expression change of MIR139-5p (p-value < 0.23) contributes to the expression change of DDX3X (p-value < 0.14) via miRNA regulation (p-value < 1☓10-16). The expression of DDX3X was also activated at the intermediate HIV infection stage. Therefore, we suggested that DDX3X could regulate HIV replication through MIR139-5p regulations [98, 99].

Dysregulation of MIRLET7i contributes to dysfunction of the immune response and viral replication

It has been reported that the HIV-1 Vif/Vpr protein is essential for viral replication [100, 101]. We determined that the Vif/Vpr protein is regulated by MIRLET7i (p-value < 1☓10-16) and interacts with the human STAT1 protein (p-value < 1☓10-16) (Fig. 6). Therefore, the expression change of MIRLET7i (p-value < 0.079) leads to the expression change of Vif/Vpr protein (p-value < 0.318) during the intermediate HIV infection stage. It was supported by the previous observations that MIRLET7i has high homology to HIV-1 sequences (90–100 %) [102] and the immune response of T cells induced by the Vpr-STAT1 interaction was involved in regulating HIV-1 replication [100]. Moreover, by comparing miRNA expression profiles between the wildtype HIV-1 infected human CEMx174 lymphocytes and Vif/Vpr RNA silencing suppressor in HIV-1 infected human CEMx174 lymphocytes [92], Vif/Vpr RNA silencing suppressor contributes to the expression change of MIR7i (p-value < 0.107). Therefore, we proposed that the miRNA regulation of MIRLET7i to the Vif/Vpr protein is involved in the induction of the immune response and viral replication (Fig. 6).

Dysregulation of MIR10a contributes to dysfunction of immune response in CD4+ T cells and viral replication

At the second infection stage, we determined that MIR10a regulates HIV-1 proteins including Gag/Pol (p-value < 1☓10-16), which significantly interacts with small subunit (SSU) processome component, homolog (yeast) (KRR1) (p-value < 1☓10-16), also known as HIV-1 Rev binding protein 2 (HRB2) (Fig. 6). Therefore, the expression change of MIR10a (p-value < 0.86) contributes to the expression change of Gag/Pol (p-value < 0.36). It has been reported that KRR1, an evolutionarily conserved gene, is essential for viral replication [103]. Additionally, MIR10a is involved in the regulation of immune response in CD4+ T cells [104] and coxsackievirus group B type 3 (CVB3) replication [105]. Moreover, it has been shown that the copy number of loci belonging to the MIRLET7 family and MIR10a is correlated with HIV-1 viral load during viral infection [106]. Therefore, we proposed that dysregulation of MIR10a is involved in HIV-1 replication via post-transcriptional regulation of HIV-1.

Dysregulation of MIR101, MIR141, and MIR152 to HIV-1 budding and releasing through DNA hypermethylation, ubiquitin transfer, and endoplasmic reticulum-associated degradation at the virus assembly/budding infection stage

The final steps of HIV-1 replication are virus assembly and budding. It has been proposed that the HIV-1 Gag protein is involved in the regulation of these steps [107, 108]. At the late infection stage, we determined that regulation of MIR101, MIR141, and MIR152 to Gag (p-value < 1☓10-16) results in dysfunction of the infected cells through a signaling cascade of 5 proteins, PRMT1, inositol 1,4,5-trisphosphate receptor type 1 (ITPR1), autocrine motility factor receptor (AMFR), tripartite motif containing 25 (TRIM25), and ubiquitin-conjugating enzyme E2N (UBE2N) (Fig. 7). The 4 genes were also highly expressed in HIV-1-infected T cells. Therefore, the expression changes of MIR101, MIR141, and MIR152 (p-value < 0.037, p-value < 0.108, and p-value < 0.126, respectively) contributed to the expression change of Gag (p-value < 0.090). We proposed that DNA methylation-associated proteins including PRMT1 and UBE2N contributes to virus release from HIV-1-infected cells to dysregulated uninfected cells through aberrant hypermethylation of host target genes [109]. It has been also suggested that MIR101 [110], and MIR152 [111] induce aberrant DNA hypermethylation in the hepatitis B virus infection. Moreover, ubiquitin transfer is required for efficient HIV-1 release during the final step of viral replication. Therefore, we suggested that the UBE2N-associated ubiquitin transfer is induced by a signaling cascade leading to dysregulation of MIR101, MIR141, and MIR152. Additionally, we also found that ITPR1 is induced to regulate HIV-1 infectivity through endoplasmic reticulum-associated degradation [112].

Conclusions

In this study, we proposed a new method to construct IGECNs for describing cellular mechanisms by system dynamic modeling and large-scale database mining of NGS data. Furthermore, we applied the method to identify differences in the IGECNs between normal and cancer cells to investigate carcinogenesis, and between HIV-1-infected and non-infected cells to investigate pathogenic and defensive mechanisms at the reverse transcription stage, integration/replication stage, and virus assembly/budding stage based on omics databases and NGS data. Database mining provided all possible candidates for genetic and epigenetic regulatory effects and interactions in the IGECNs related to cellular functions. We used AIC and statistical assessment to prune the false-positive regulatory effects and interactions by applying the dynamic coupling model to NGS data. Finally, we compared the major differences in IGECNs between different biological conditions to identify mechanistic differences to explore the evolution of cellular mechanisms. Based on these network comparisons, we examined how genetic and epigenetic regulation mechanism affects gastric tumorigenesis, hepatocarcinogenesis, and the progression of HIV-1 infections in CD4+ T cells. Therefore, the proposed IGECN construction method allowed us to unravel cellular network mechanisms from genetic and epigenetic mechanisms using omics databases and NGS data. A cancer hallmark network framework has been proposed 9 cancer hallmark networks, including angiogenesis-inducing network, mutating network, dedifferentiation network, EMT (epithelial-mesenchymal transition) network, genome duplication network, immune-escaping network, survival network, metabolic network, and Stroma-network [113]. By comparing the hallmark network difference between normal and cancer cells, cancer clonal evolution and clinical phenotypes, which will have impact on diagnosis and personalized treatment, can be predicted.

Abbreviations

IGECN:

integrated genetic and epigenetic cellular network

NGS:

next-generation sequencing

GRNs:

gene regulatory networks

PPINs:

protein-protein interaction networks

miRNAs:

microRNAs

AIDS:

acquired immune deficiency syndrome

HIV:

human immunodeficiency virus

TFs:

transcription factors

RNA-seq:

RNA sequencing

GO:

gene ontology

AIC:

Akaike information criterion

GSEA:

gene set enrichment analysis

UBC:

ubiquitin C

HCC:

hepatocellular carcinoma

USP15:

ubiquitin specific peptidase 15

SRCs:

p160 steroid receptor coactivator family

NCOAs:

nuclear receptor coactivator family

PML:

promyelocytic leukemia

ETS1:

v-ets erythroblastosis virus E26 oncogene homolog 1

SUMOs:

small ubiquitin-related modifier family

PRMT1:

protein arginine methyltransferase 1

TCGA:

the cancer genome atlas

References

  1. Blair JD, Price EM. Illuminating Potential Technical Artifacts of DNA-Methylation Array Probes. Am J Hum Genet. 2012;91(4):760–2.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  2. Chatr-aryamontri A, Breitkreutz BJ, Heinicke S, Boucher L, Winter A, Stark C, et al. The BioGRID interaction database: 2013 update. Nucleic Acids Res. 2013;41(D1):D816–23.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  3. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  4. Broen JCA, Radstake TRDJ, Rossato M. The role of genetics and epigenetics in the pathogenesis of systemic sclerosis. Nat Rev Rheumatol. 2014;10(11):671–81.

    Article  PubMed  CAS  Google Scholar 

  5. Ke XY, Cortina-Borja M, Silva BC, Lowe R, Rakyan V, Balding D. Integrated analysis of genome-wide genetic and epigenetic association data for identification of disease mechanisms. Epigenetics. 2013;8(11):1236–44.

    Article  PubMed  CAS  Google Scholar 

  6. Yeger-Lotem E, Sattath S, Kashtan N, Itzkovitz S, Milo R, Pinter RY, et al. Network motifs in integrated cellular networks of transcription-regulation and protein-protein interaction. Proc Natl Acad Sci U S A. 2004;101(16):5934–9.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  7. Lin CC, Chen YJ, Chen CY, Oyang YJ, Juan HF, Huang HC. Crosstalk between transcription factors and microRNAs in human protein interaction network. BMC Syst Biol. 2012;6:18.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  8. Herrgard MJ, Lee BS, Portnoy V, Palsson BO. Integrated analysis of regulatory and metabolic networks reveals novel regulatory mechanisms in Saccharomyces cerevisiae. Genome Res. 2006;16(5):627–35.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  9. Przytycka TM, Kim YA. Network integration meets network dynamics. BMC Biol. 2010;8:48.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  10. Wang YC, Chen BS. Integrated cellular network of transcription regulations and protein-protein interactions. BMC Syst Biol. 2010;4:20.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  11. Hsu HJ, Lee CF, Kaunas R. A dynamic stochastic model of frequency-dependent stress fiber alignment induced by cyclic stretch. Plos One. 2009;4(3):e4853.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  12. Chen LZ, Nguang SK, Chen X. Modelling and optimization of biotechnological processes : artificial intelligence approaches. Berlin: Springer; 2006.

    Google Scholar 

  13. Chen LZ, Nguang SK, Chen XD. On-line identification and optimization of feed rate profiles for high productivity fed-batch culture of hybridoma cells using genetic algorithms. ISA Trans. 2002;41(4):409–19.

    Article  PubMed  Google Scholar 

  14. Chuang CH, Lin CL, Chang YC, Jennawasin T, Chen PK. Design of synthetic biological logic circuits based on evolutionary algorithm. IET Syst Biol. 2013;7(4):89–105.

    Article  PubMed  Google Scholar 

  15. Huang CH, Jan HY, Lin CL, Lee CS. System identification: DNA computing approach. ISA Trans. 2009;48(3):254–63.

    Article  PubMed  Google Scholar 

  16. Chen BS, Li CW. Robust Observer-Based Tracking Control of Hodgkin-Huxley Neuron Systems Under Environmental Disturbances. Neural Comput. 2010;22(12):3143–78.

    Article  PubMed  Google Scholar 

  17. Chen BS, Li CW. On the Interplay between Entropy and Robustness of Gene Regulatory Networks. Entropy. 2010;12(5):1071–101.

    Article  CAS  Google Scholar 

  18. Chen BS, Li CW. On the Noise-Enhancing Ability of Stochastic Hodgkin-Huxley Neuron Systems. Neural Comput. 2010;22(7):1737–63.

    Article  PubMed  Google Scholar 

  19. Chen BS, Li CW. Measuring information flow in cellular networks by the systems biology method through microarray data. Frontiers in Plant Science. 2015;6:390.

    PubMed Central  PubMed  Google Scholar 

  20. Chen BS, Tsai KW, Li CW. Using Nonlinear Stochastic Evolutionary Game Strategy to Model an Evolutionary Biological Network of Organ Carcinogenesis Under a Natural Selection Scheme. Evol Bioinforma. 2015;11:155–78.

    Article  Google Scholar 

  21. Chen BS, Wong SW, Li CW. On the Calculation of System Entropy in Nonlinear Stochastic Biological Networks. Entropy. 2015;17(10):6801–33.

    Article  Google Scholar 

  22. Li CW, Chen BS. Identifying Functional Mechanisms of Gene and Protein Regulatory Networks in Response to a Broader Range of Environmental Stresses. Comp Funct Genomics. 2010;2010:408705.

  23. Li CW, Lo CC, Chen BS. Robust sensorimotor control of human arm model under state-dependent noises, control-dependent noises and additive noises. Neurocomputing. 2015;167:61–75.

    Article  Google Scholar 

  24. Guintivano J, Brown T, Newcomer A, Jones M, Cox O, Maher BS, et al. Identification and Replication of a Combined Epigenetic and Genetic Biomarker Predicting Suicide and Suicidal Behaviors. Am J Psychiatr. 2014;171(12):1287–96.

    Article  PubMed  Google Scholar 

  25. Mohammadi P, Desfarges S, Bartha I, Joos B, Zangger N, Munoz M, et al. 24 Hours in the Life of HIV-1 in a T Cell Line. Plos Pathogens. 2013;9(1):e1003161.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  26. Kim SY, Park C, Kim HJ, Park J, Hwang J, Kim JI, et al. Deregulation of Immune Response Genes in Patients With Epstein-Barr Virus-Associated Gastric Cancer and Outcomes. Gastroenterology. 2015;148(1):137–47.

  27. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  28. Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M. BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006;34:D535–9.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  29. Bovolenta LA, Acencio ML, Lemke N. HTRIdb: an open-access database for experimentally verified human transcriptional regulation interactions. BMC Genomics. 2012;13:405.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  30. Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1 alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73.

    Article  PubMed  CAS  Google Scholar 

  31. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  32. Zheng GY, Tu K, Yang Q, Xiong Y, Wei CC, Xie L, et al. ITFP: an integrated platform of mammalian transcription factors. Bioinformatics. 2008;24(20):2416–7.

    Article  PubMed  CAS  Google Scholar 

  33. Griffiths-Jones S. The microRNA Registry. Nucleic Acids Res. 2004;32:D109–11.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  34. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ. miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008;36:D154–8.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  35. Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. The role of site accessibility in microRNA target recognition. Nat Genet. 2007;39(10):1278–84.

    Article  PubMed  CAS  Google Scholar 

  36. Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39:D152–7.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  37. Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(D1):D68–73.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  38. Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.

  39. Friedman RC, Farh KKH, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009;19(1):92–105.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  40. Garcia DM, Baek D, Shin C, Bell GW, Grimson A, Bartel DP. Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs. Nat Struct Mol Biol. 2011;18(10):1139–46.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  41. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007;27(1):91–105.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  42. Cui QH, Ma Y, Jaramillo M, Bari H, Awan A, Yang S, et al. A map of human cancer signaling. Mol Syst Biol. 2007;3:152.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  43. Huang YZ, Sun JJ, Zhang LZ, Li CJ, Womack JE, Li ZJ, et al. Genome-wide DNA Methylation Profiles and Their Relationships with mRNA and the microRNA Transcriptome in Bovine Muscle Tissue (Bos taurine). Sci Rep. 2014;4:1–17.

    Google Scholar 

  44. Muckstein U, Tafer H, Hackermuller J, Bernhart SH, Stadler PF, Hofacker IL. Thermodynamics of RNA-RNA binding. Bioinformatics. 2006;22(10):1177–82.

    Article  PubMed  CAS  Google Scholar 

  45. Nandi A, Vaz C, Bhattacharya A, Ramaswamy R. miRNA-regulated dynamics in circadian oscillator models. BMC Syst Biol. 2009;3:45.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  46. Ferreira AEN, Freire AMJP, Voit EO. A quantitative model of the generation of N-epsilon-(carboxymethyl)lysine in the Maillard reaction between collagen and glucose. Biochem J. 2003;376:109–21.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  47. Voit EO, Ferreira AEN. Computational analysis of biochemical systems : a practical guide for biochemists and molecular biologists. New York: Cambridge University Press; 2000.

    Google Scholar 

  48. Akaike H. New Look at Statistical-Model Identification. Ieee Transactions on Automatic Control. 1974;Ac19(6):716–23.

    Article  Google Scholar 

  49. Johansson R. System modeling and identification. Englewood Cliffs: Prentice Hall; 1993.

    Google Scholar 

  50. Bar-Joseph Z, Gerber G, Gifford DK, Jaakkola TS, Simon I. A new approach to analyzing gene expression time series data. In: Proceedings of the sixth annual international conference on Computational biology. Washington, DC: ACM; 2002. p. 39–48.

    Google Scholar 

  51. De Boor C. A practical guide to splines : with 32 figures. Revth ed. New York: Springer; 2001.

    Google Scholar 

  52. Faires JD, Burden RL. Numerical methods. 2nd ed. Pacific Grove: Brooks/Cole Pub. Co.; 1998.

    Google Scholar 

  53. Coleman TF, Hulbert LA. A Direct Active Set Algorithm for Large Sparse Quadratic Programs with Simple Bounds. Math Program. 1989;45(3):373–406.

    Article  Google Scholar 

  54. Gill PE, Murray W, Wright MH. Practical optimization. London; New York: Academic Press; 1981.

    Google Scholar 

  55. Lu P, Vogel C, Wang R, Yao X, Marcotte EM. Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol. 2007;25(1):117–24.

    Article  PubMed  CAS  Google Scholar 

  56. Seber GAF, Lee AJ. Linear regression analysis. 2nd ed. Hoboken: Wiley-Interscience; 2003.

    Book  Google Scholar 

  57. Abraham RT. Cell cycle checkpoint signaling through the ATM and ATR kinases. Genes Dev. 2001;15(17):2177–96.

    Article  PubMed  CAS  Google Scholar 

  58. McMahon M. RAF translocations expand cancer targets. Nat Med. 2010;16(7):749–50.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  59. Wilhelm E, Doyle MC, Nzaramba I, Magdzinski A, Dumais N, Bell B. CTGC motifs within the HIV core promoter specify Tat-responsive pre-initiation complexes. Retrovirology. 2012;9:62.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  60. Sharma D, Bhattacharya J. Evolutionary Constraints Acting on DDX3X Protein Potentially Interferes with Rev-Mediated Nuclear Export of HIV-1 RNA. Plos One. 2010;5(3):e9613.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  61. Kang M-S, Kieff E. Epstein–Barr virus latent genes. Exp Mol Med. 2015;47(1):e131.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  62. Liu L, Ouyang M, Rao JN, Zou T, Xiao L, Chung HK, et al. Competition between RNA-binding proteins CELF1 and HuR modulates MYC translation and intestinal epithelium renewal. Mol Biol Cell. 2015;26(10):1797–810.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  63. Xie LM, Zhang ZW, Tan ZQ, He RF, Zeng X, Xie YJ, et al. microRNA-124 inhibits proliferation and induces apoptosis by directly repressing EZH2 in gastric cancer. Mol Cell Biochem. 2014;392(1-2):153–9.

    Article  PubMed  CAS  Google Scholar 

  64. Zhao XN, Li XW, Yuan HQ. microRNAs in gastric cancer invasion and metastasis. Frontiers in Bioscience-Landmark. 2013;18:803–10.

    Article  CAS  Google Scholar 

  65. Chang S, Gao L, Yang Y, Tong D, Guo B, Liu L, et al. miR-145 mediates the antiproliferative and gene regulatory effects of vitamin D3 by directly targeting E2F3 in gastric cancer cells. Oncotarget. 2015;6(10):7675–85.

    Article  PubMed Central  PubMed  Google Scholar 

  66. Zheng LD, Pu JR, Qi T, Qi M, Li D, Xiang X, et al. miRNA-145 Targets v-ets Erythroblastosis Virus E26 Oncogene Homolog 1 to Suppress the Invasion, Metastasis, and Angiogenesis of Gastric Cancer Cells. Mol Cancer Res. 2013;11(2):182–93.

    Article  PubMed  CAS  Google Scholar 

  67. Cao DH, Jiang J, Tsukamoto T, Liu RM, Ma L, Jia ZF, et al. Canolol Inhibits Gastric Tumors Initiation and Progression through COX-2/PGE2 Pathway in K19-C2mE Transgenic Mice. Plos One. 2015;10(3):e0120938.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  68. Roberts TL, Ho U, Luff J, Lee CS, Apte SH, MacDonald KPA, et al. Smg1 haploinsufficiency predisposes to tumor formation and inflammation. Proc Natl Acad Sci U S A. 2013;110(4):E285–94.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  69. Ishida J. Clinicopathological study of the amount of LUCA15 expressed in gastric cancer tissue. Hokkaido Igaky Zasshi. 2006;81(5):451–9.

    CAS  Google Scholar 

  70. Peng J, Valeshabad AK, Li QF, Wang Y. Differential expression of RBM5 and KRAS in pancreatic ductal adenocarcinoma and their association with clinicopathological features. Oncology Letters. 2013;5(3):1000–4.

    PubMed Central  PubMed  CAS  Google Scholar 

  71. Zou Q, Jin J, Hu HB, Li HYS, Romano S, Xiao YC, et al. USP15 stabilizes MDM2 to mediate cancer-cell survival and inhibit antitumor T cell responses. Nat Immunol. 2014;15(6):562–70.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  72. Guidoboni M, Zancai P, Cariati R, Rizzo S, Dal Col J, Pavan A, et al. Retinoic acid inhibits the proliferative response induced by CD40 activation and interleukin-4 in mantle cell lymphoma. Cancer Res. 2005;65(2):587–95.

    PubMed  CAS  Google Scholar 

  73. Zhou S, Ye WG, Ren J, Shao QJ, Qi YH, Liang J, et al. MicroRNA-381 increases radiosensitivity in esophageal squamous cell carcinoma. American Journal of Cancer Research. 2015;5(1):267–77.

    PubMed Central  PubMed  Google Scholar 

  74. Mani A, Gelmann EP. The ubiquitin-proteasome pathway and its role in cancer. J Clin Oncol. 2005;23(21):4776–89.

    Article  PubMed  CAS  Google Scholar 

  75. Palmieri C, Gojis O, Rudraraju B, Stamp-Vincent C, Wilson D, Langdon S, et al. Expression of steroid receptor coactivator 3 in ovarian epithelial cancer is a poor prognostic factor and a marker for platinum resistance. Br J Cancer. 2013;108(10):2039–44.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  76. Nishida T, Terashima M, Fukami K, Yamada Y. PIASy controls ubiquitination-dependent proteasomal degradation of Ets-1. Biochem J. 2007;405:481–8.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  77. Lacaud G, Gore L, Kennedy M, Kouskoff V, Kingsley P, Hogan C, et al. Runx1 is essential for hematopoietic commitment at the hemangioblast stage of development in vitro. Blood. 2002;100(2):458–66.

    Article  PubMed  CAS  Google Scholar 

  78. Iwanami A, Gini B, Zanca C, Matsutani T, Assuncao A, Nael A, et al. PML mediates glioblastoma resistance to mammalian target of rapamycin (mTOR)-targeted therapies. Proc Natl Acad Sci U S A. 2013;110(11):4339–44.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  79. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global Cancer Statistics. Ca-a Cancer Journal for Clinicians. 2011;61(2):69–90.

    Article  PubMed  Google Scholar 

  80. Itahana Y, Neo SH, Itahana K. miR-141, a new player, joins the senescence orchestra. Cell Cycle. 2013;12(23):3586–7.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  81. Watanabe S, Watanabe K, Akimov V, Bartkova J, Blagoev B, Lukas J, et al. JMJD1C demethylates MDC1 to regulate the RNF8 and BRCA1-mediated chromatin response to DNA breaks. Nat Struct Mol Biol. 2013;20(12):1425–33.

    Article  PubMed  CAS  Google Scholar 

  82. Zheng TS, Wang JB, Song XA, Meng XZ, Pan SH, Jiang HC, et al. Nutlin-3 cooperates with doxorubicin to induce apoptosis of human hepatocellular carcinoma cells through p53 or p73 signaling pathways. J Cancer Res Clin Oncol. 2010;136(10):1597–604.

    Article  PubMed  CAS  Google Scholar 

  83. Lamy E, Herz C, Lutz-Bonengel S, Hertrampf A, Marton MR, Mersch-Sundermann V. The MAPK Pathway Signals Telomerase Modulation in Response to Isothiocyanate-Induced DNA Damage of Human Liver Cancer Cells. Plos One. 2013;8(1):e53240.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  84. Yang SF, Chang CW, Wei RJ, Shiue YL, Wang SN, Yeh YT. Involvement of DNA Damage Response Pathways in Hepatocellular Carcinoma. Biomed Research International. 2014;2014:153867.

    PubMed Central  PubMed  Google Scholar 

  85. Niemi NM, Lanning NJ, Westrate LM, MacKeigan JP. Downregulation of the Mitochondrial Phosphatase PTPMT1 Is Sufficient to Promote Cancer Cell Death. Plos One. 2013;8(1):e53803.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  86. Hatakeyama S. TRIM proteins and cancer. Nat Rev Cancer. 2011;11(11):792–804.

    Article  PubMed  CAS  Google Scholar 

  87. Li H, Sun L, Chen X, Xiong W, Hu D, Jie S. Microvesicle microRNA profiles and functional roles between chronic hepatitis B and hepatocellular carcinoma. Clin Transl Oncol. 2014;16(3):315–21.

    Article  PubMed  CAS  Google Scholar 

  88. Song RP, Song HW, Liang YJ, Yin DL, Zhang H, Zheng TS, et al. Reciprocal Activation Between ATPase Inhibitory Factor 1 and NF-kappa B Drives Hepatocellular Carcinoma Angiogenesis and Metastasis. Hepatology. 2014;60(5):1659–73.

    Article  PubMed  CAS  Google Scholar 

  89. Shen RR, Zhou AY, Kim E, O’Connell JT, Hagerstrand D, Beroukhim R, et al. TRAF2 is an NF-kappa B-activating oncogene in epithelial cancers. Oncogene. 2015;34(2):209–16.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  90. Brady CW. Liver disease in menopause. World J Gastroenterol. 2015;21(25):7613–20.

    Article  PubMed Central  PubMed  Google Scholar 

  91. Tian Y, Xie X, Lin Y, Tan G, Zhong W. Androgen receptor in hepatocarcinogenesis: Recent developments and perspectives (Review). Oncology Letters. 2015;9(5):1983–8.

    PubMed Central  PubMed  Google Scholar 

  92. Hayes AM, Qian SM, Yu LB, Boris-Lawrie K. Tat RNA silencing suppressor activity contributes to perturbation of lymphocyte miRNA by HIV-1. Retrovirology. 2011; 8:36.

  93. Yao H, Ma R, Yang L, Hu G, Chen XF, Duan M, et al. MiR-9 promotes microglial activation by targeting MCPIP1. Nat Commun. 2014;5:4386.

    PubMed Central  PubMed  CAS  Google Scholar 

  94. Taylor HE, Linde ME, Khatua AK, Popik W, Hildreth JEK. Sterol Regulatory Element-Binding Protein 2 Couples HIV-1 Transcription to Cholesterol Homeostasis and T Cell Activation. J Virol. 2011;85(15):7699–709.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  95. Oh Y, Chung KC. Small Ubiquitin-like Modifier (SUMO) Modification of Zinc Finger Protein 131 Potentiates Its Negative Effect on Estrogen Signaling. J Biol Chem. 2012;287(21):17517–29.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  96. Szotek EL, Narasipura SD, Al-Harthi L. 17 beta-Estradiol inhibits HIV-1 by inducing a complex formation between beta-catenin and estrogen receptor alpha on the HIV promoter to suppress HIV transcription. Virology. 2013;443(2):375–83.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  97. Miller-Jensen K, Skupsky R, Shah PS, Arkin AP, Schaffer DV. Genetic Selection for Context-Dependent Stochastic Phenotypes: Sp1 and TATA Mutations Increase Phenotypic Noise in HIV-1 Gene Expression. Plos Computational Biology. 2013;9(7):e1003135.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  98. Mahboobi SH, Javanpour AA, Mofrad MRK. The Interaction of RNA Helicase DDX3 with HIV-1 Rev-CRM1-RanGTP Complex during the HIV Replication Cycle. Plos One. 2015;10(2):e0112969.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  99. Yedavalli VSRK, Jeang KT. Rev-ing up post-transcriptional HIV-1 RNA expression. RNA Biol. 2011;8(2):195–9.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  100. Zahoor MA, Xue GG, Sato H, Murakami T, Takeshima S, Aida Y. HIV-1 Vpr Induces Interferon-Stimulated Genes in Human Monocyte-Derived Macrophages. Plos One. 2014;9(8):e106418.

    Article  PubMed Central  PubMed  Google Scholar 

  101. Zhang H, Pomerantz RJ, Dornadula G, Sun Y. Human immunodeficiency virus type 1 Vif protein is an integral component of an mRNP complex of viral RNA and could be involved in the viral RNA folding and packaging process. J Virol. 2000;74(18):8252–61.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  102. Kanak M, Alseiari M, Balasubramanian P, Addanki K, Aggarwal M, Noorali S, et al. Triplex-forming MicroRNAs Form Stable Complexes With HIV-1 Provirus and Inhibit its Replication. Appl Immunohistochem Mol Morphol. 2010;18(6):532–45.

    Article  PubMed  Google Scholar 

  103. Panaro MA, Mitolo V, Cianciulli A, Cavallo P, Mitolo CI, Acquafredda A. The HIV-1 rev binding family of proteins: The dog proteins as a study model. Endocrine Metabolic & Immune Disorders-Drug Targets. 2008;8(1):30–46.

    Article  CAS  Google Scholar 

  104. Kelada S, Sethupathy P, Okoye IS, Kistasis E, Czieso S, White SD, et al. miR-182 and miR-10a Are Key Regulators of Treg Specialisation and Stability during Schistosome and Leishmania-associated Inflammation. Plos Pathogens. 2013;9(6):e1003451.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  105. Tong L, Lin LX, Wu S, Guo ZW, Wang TY, Qin Y, et al. MiR-10a(star) up-regulates coxsackievirus B3 biosynthesis by targeting the 3D-coding sequence. Nucleic Acids Res. 2013;41(6):3760–71.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  106. Abdel-Mohsen M, Deng XT, Danesh A, Liegler T, Jacobs ES, Rauch A, et al. Role of MicroRNA Modulation in the Interferon-alpha/Ribavirin Suppression of HIV-1 In Vivo. Plos One. 2014;9(10):e109220.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  107. Grigorov B, Arcanger F, Roingeard P, Darlix JL, Muriaux D. Assembly of infectious HIV-1 in human epithelial and T-lymphoblastic cell lines. J Mol Biol. 2006;359(4):848–62.

    Article  PubMed  CAS  Google Scholar 

  108. De Clercq E. Current lead natural products for the chemotherapy of human immunodefiency virus (HIV) infection. Med Res Rev. 2000;20(5):323–49.

    Article  PubMed  Google Scholar 

  109. Luzzi A, Morettini F, Gazaneo S, Mundo L, Onnis A, Mannucci S, et al. HIV-1 Tat induces DNMT over-expression through microRNA dysregulation in HIV-related non Hodgkin lymphomas. Infectious Agents and Cancer. 2014;9:41.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  110. Wei XF, Xiang TX, Ren GS, Tan C, Liu R, Xu X, et al. miR-101 is down-regulated by the hepatitis B virus x protein and induces aberrant DNA methylation by targeting DNA methyltransferase 3A. Cell Signal. 2013;25(2):439–46.

    Article  PubMed  CAS  Google Scholar 

  111. Huang JF, Wang Y, Guo YJ, Sun SH. Down-Regulated MicroRNA-152 Induces Aberrant DNA Methylation in Hepatitis B Virus-Related Hepatocellular Carcinoma by Targeting DNA Methyltransferase 1. Hepatology. 2010;52(1):60–70.

    Article  PubMed  CAS  Google Scholar 

  112. Casini A, Olivieri M, Vecchi L, Burrone OR, Cereseto A. Reduction of HIV-1 Infectivity through Endoplasmic Reticulum-Associated Degradation-Mediated Env Depletion. J Virol. 2015;89(5):2966–71.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  113. Wang E, Zaman N, Mcgee S, Milanese JS, Masoudi-Nejad A, O’Connor-McCourt M. Predictive genomics: A cancer hallmark network framework for predicting tumor clinical phenotypes using genome sequencing data. Semin Cancer Biol. 2015;30:4–12.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

The work supported by the Ministry of Science and Technology of Taiwan under grant No. MOST 103-2745-E-007-001-ASP.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Bor-Sen Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

BSC & CWL: Methodology development, conception and design, data analysis and interpretation, manuscript writing. Both authors read and approved the final manuscript.

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

Chen, BS., Li, CW. Constructing an integrated genetic and epigenetic cellular network for whole cellular mechanism using high-throughput next-generation sequencing data. BMC Syst Biol 10, 18 (2016). https://doi.org/10.1186/s12918-016-0256-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12918-016-0256-5

Keywords