Research article Special Issues

Caputo-Hadamard fractional differential equations with nonlocal fractional integro-differential boundary conditions via topological degree theory

  • This article aims to prove the existence and uniqueness of solutions to a nonlinear boundary value problem of fractional differential equations involving the Caputo-Hadamard fractional derivative with nonlocal fractional integro-differential boundary conditions. The concerned results are obtained employing topological degree for condensing maps via a priori estimate method and the Banach contraction principle fixed point theorem. Besides, two illustrative examples are presented.

    Citation: Choukri Derbazi, Hadda Hammouche. Caputo-Hadamard fractional differential equations with nonlocal fractional integro-differential boundary conditions via topological degree theory[J]. AIMS Mathematics, 2020, 5(3): 2694-2709. doi: 10.3934/math.2020174

    Related Papers:

    [1] Ping Yan, Gerardo Chowell . Modeling sub-exponential epidemic growth dynamics through unobserved individual heterogeneity: a frailty model approach. Mathematical Biosciences and Engineering, 2024, 21(10): 7278-7296. doi: 10.3934/mbe.2024321
    [2] Jun Zhai, Bilin Xu . Research on meme transmission based on individual heterogeneity. Mathematical Biosciences and Engineering, 2021, 18(5): 5176-5193. doi: 10.3934/mbe.2021263
    [3] Sarafa A. Iyaniwura, Musa Rabiu, Jummy F. David, Jude D. Kong . Assessing the impact of adherence to Non-pharmaceutical interventions and indirect transmission on the dynamics of COVID-19: a mathematical modelling study. Mathematical Biosciences and Engineering, 2021, 18(6): 8905-8932. doi: 10.3934/mbe.2021439
    [4] Huiping Zang, Shengqiang Liu, Yi Lin . Evaluations of heterogeneous epidemic models with exponential and non-exponential distributions for latent period: the Case of COVID-19. Mathematical Biosciences and Engineering, 2023, 20(7): 12579-12598. doi: 10.3934/mbe.2023560
    [5] Sara Y. Del Valle, J. M. Hyman, Nakul Chitnis . Mathematical models of contact patterns between age groups for predicting the spread of infectious diseases. Mathematical Biosciences and Engineering, 2013, 10(5&6): 1475-1497. doi: 10.3934/mbe.2013.10.1475
    [6] Liqiong Pu, Zhigui Lin . A diffusive SIS epidemic model in a heterogeneous and periodically evolvingenvironment. Mathematical Biosciences and Engineering, 2019, 16(4): 3094-3110. doi: 10.3934/mbe.2019153
    [7] Pengfei Liu, Yantao Luo, Zhidong Teng . Role of media coverage in a SVEIR-I epidemic model with nonlinear incidence and spatial heterogeneous environment. Mathematical Biosciences and Engineering, 2023, 20(9): 15641-15671. doi: 10.3934/mbe.2023698
    [8] Yongli Cai, Yun Kang, Weiming Wang . Global stability of the steady states of an epidemic model incorporating intervention strategies. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1071-1089. doi: 10.3934/mbe.2017056
    [9] Ahmed Alshehri, Saif Ullah . A numerical study of COVID-19 epidemic model with vaccination and diffusion. Mathematical Biosciences and Engineering, 2023, 20(3): 4643-4672. doi: 10.3934/mbe.2023215
    [10] Shubham Chaudhry, Gauri Agrawal, Maia Martcheva, A. K. Misra . Modeling the impact of temperature on the dynamics of carrier-dependent infectious diseases with control strategies. Mathematical Biosciences and Engineering, 2025, 22(7): 1722-1750. doi: 10.3934/mbe.2025063
  • This article aims to prove the existence and uniqueness of solutions to a nonlinear boundary value problem of fractional differential equations involving the Caputo-Hadamard fractional derivative with nonlocal fractional integro-differential boundary conditions. The concerned results are obtained employing topological degree for condensing maps via a priori estimate method and the Banach contraction principle fixed point theorem. Besides, two illustrative examples are presented.


    With the advent of technologies allowing for large-scale, high throughput data, a much clearer understanding of the genomic mechanisms behind gene regulation have been gained. The scientists found that there are unexpected far more noncoding RNAs comparing with protein-coding genes and, and these noncoding regions play important roles in determining the complexity observed in the human genome [1,2]. Within these noncoding regions, long noncoding RNAs (lncRNAs), which are functionally defined as noncoding regions of RNA that are at least 200 base-pairs in length, have attracted lots of attention. Certain lncRNAs appear to act locally, while others have more distal regulatory effects, even acting across multiple chromosomes [3]. Many studies have identified specific functions of particular lncRNAs, including embryonic mechanisms, cell cycle functions, innate immunity, and disease processes. However, there are still thousands of lncRNAs have no identified functions [1,3,4,5,6]. Some studies have been performed that produce relatively few numbers of lncRNA functions [7], and have shown that the function of lncRNAs is highly cell-type-specific: one lncRNA may inhibit particular genes in one type of cell while promoting the same gene in another. This phenomenon makes it even more difficult to identify lncRNA functions on a large scale. Due to this specificity, researchers propose that future lncRNA studies should be performed on specific cell types to identify particular regulatory mechanisms.

    One of the most prominent and intriguing applications of lncRNA regulatory investigation comes from cancer studies [8,9]. It has been shown that lncRNAs appear to have high connectivity with numerous diseases, especially cancer. Because of the highly cell type-specific nature of lncRNA regulatory functions and the irregularity of cancer cell genetic information, studying lncRNA regulation in specific cancer types may provide promising insight into specific genomic regulations of common cancer cells. In a few documented cases, specific lncRNAs have been shown to be significantly differentially expressed in specific cancer types, such as prostate cancer and breast cancer [1]. For these reasons, it seems appropriate to further investigate lncRNA-gene interactions in particular cancer cells.

    The wealth of gene expression datasets available provides an opportunity to computationally identify co-expressed gene modules(CEMs), each of which is defined as a highly structured expression pattern on a specific gene set [10,11]. These CEMs tend to be functionally related or co-regulated by the same transcriptional regulatory signals (e.g., transcription factors, lncRNA and so on) under a specific condition or in a particular disease cell type. Overall, successful derivation of the CEMs may grant a higher-level interpretation of large-scale gene expression data, improve functional annotation of condition-specific gene activities, facilitate inference of gene regulatory relationships, hence, provide a better mechanism level understanding of complex diseases.

    The computational identification of CEMs can be solved by a biclustering approach [12], which is a two-dimensional data mining technique that simultaneously identifies co-expressed genes under a subset of conditions. a high proportion of enriched biclusters on real datasets. Within this study, we try to identify new lncRNA-gene interactions and transcription factor-lncRNA partnerships from cancer RNA-seq data using a biclustering approach. The biclustering method will allow for the identification of particular expression patterns across multiple datasets, indicating networks of lncRNA and gene interactions. This developed method will also provide a framework for future lncRNA interaction studies. We applied this method on two sets of TCGA breast cancer RNA-seq data to generated CEMs based on known lncRNA-gene interactions. Then, the predicted CEMs are linked to lncRNA by a statistic p-value and the new lncRNA-gene relationship are generated. The evaluation on the predicted results showed that the pipeline can find some target genes for given lncRNA, and meanwhile the performance still has some space to be improved. We further conducted a TF motif analysis on the predicted CEMs and provide potential regulation cooperation between TFs and lncRNAs. The related original data with codes, results and supplementary data can be downloaded on https://github.com/IvesG/sGavin.git.

    Two sets of TCGA (The Cancer Genome Atlas) breast cancer RNA-seq data, one from the normal cell (referred as normal data) and the other from tumor cell (referred as tumor data) were downloaded from https://portal.gdc.cancer.gov/. The normal and tumor data consist of 113 and 1091 samples, respectively. And of the 113 normal samples, 112 of them are from the same patient among the tumors. Both datasets contain 60,483 genes, among which there are 19,824 protein-coding genes and 7,399 long intergenic noncoding RNAs (lincRNAs) genes. The RNA-seq data are all Upper Quartile normalized FPKM (UQ-FPKM) values.

    A total of 1,081 experimentally validated lncRNA-associated regulatory entries were downloaded from LncReg [13], describing the comprehensive regulatory relationships among 258 lncRNAs and 571 genes. All these relationships were manually collected from PubMed with focus on the data generated by laboratory methods, and can be categorized into up/down/active/inactive based on regulatory relationships or transcription/post-transcription/translation/post-translation based on regulatory mechanisms.

    As we focus on lncRNA-gene interactions, the relationships downloaded from LncReg were filtered to retain only relationships describing genes regulated by lncRNAs with specified species information (constrained to Homo sapiens and Mus musculus), resulting 925 relationships in total for the downstream analysis, covering interactions between 309 unique human genes and 103 human lncRNAs, as well as between 199 mouse genes and 100 mouse lncRNAs. It is noteworthy that these 925 relationships include 28 post-transcriptional regulations, 41 post-translational regulations, 714 transcriptional regulations, 23 translational regulations, 1 transcriptional & translational regulation and 118 unspecified relationships.

    As the table from LncReg [13] only provides gene symbols, while the RNA-seq dataset uses Ensembl ID as gene's identifiers, we use Ensembl BioMart [14] to match gene symbols with Ensembl IDs for all the genes and lncRNAs. Then we got orthologous genes between mouse and human also using BioMart; we found orthologous human genes for all 199 mouse genes, and 38 overlapped with original human genes. For convenience, we recorded human genes, mouse genes that don't overlap with human genes, human lncRNAs and mouse lncRNAs that don't overlap with human lncRNAs as HG, MG, HL, and ML, respectively.

    We combined the normal and tumor RNA-seq dataset together, then extracted expression values for all the HG, MG, HL, ML, protein-coding genes (PC, the remaining protein-coding genes except HG and MG) and lincRNAs (linc, the remaining lincRNA except HL and ML). Taking the genes as rows and the conditions as columns, we obtained the RNA-seq expression matrix on which biclustering will be performed to detect CEMs.

    QUBIC is a biclustering analysis tool designed for co-expression analyses of genes based on their gene-expression patterns under multiple conditions. The software can generally identify all statistically significant groups, or biclusters, of genes with similar expression patterns under at least a specific number of experimental conditions, which tend to be more sensitive and more specific than other biclustering tools [15]. We use a quantile-based discretization method of QUBIC to generate a qualitative representing matrix for the RNA-seq expression matrix. Then we extracted the rows of known lncRNA regulated HG and MG from this representing matrix as seed 1 and HG, MG, HL, and ML rows as seed 2. Next bi-clustering analysis was performed on these two seeds to predict co-expressed gene modules (CEMs) in the qualitative representing matrix, respectively.

    For an identified CEMs, we calculated the P-value of a bicluster enriched with genes regulated by a lncRNA using the hypergeometric function [16],

    Pr(r|N,K,n)=(Kr)(NKnr)(Nn)

    where r is the number of genes in a CEMs (with size n) that regulated by certain lncRNA, N is the total number of known lncRNA regulated genes in the whole genome, K is the number of genes regulated by that lncRNA in the whole genome.

    We assumed that, if the known target genes of a given lncRNA are highly covered by a CEM with a significant p-value, the other genes in this CEM have high possibilities regulated by the given lncRNA. Thus, we used the smallest P-value for all possible lncRNAs as the p-value of the current bicluster and the relationships between lncRNA and genes in the bicluster are predicted.

    To evaluate the performance of the new methods on the prediction of new relationships between lncRNA and genes, we randomly separate seed2 into two parts with equal size named seedpart1 and seedpart2, for multiple times. Then bi-clustering analysis will be performed on seedpart2 to predict co-expressed gene modules (CEMs). For seedpart1 we find its part which is covered by co-expressed gene modules (CEMs) from seedpart2. We calculate the cover ratios by the size of seedpart1 to be divided by the size of the covered part by CEMs generated from seedpart2. Also, we calculate the p-values for the coverage rates to present the statistical significance of them.

    We choose several significant CEMs with sizes or conditions below 100, to conducted TF motif analysis. The promoter regions of the corresponding genes are inputted into the sub-routine findMotifs.pl of Homer [17], respectively. The script findMotifs.pl can firstly search for the upstream promoter sequences of a certain length automatically, and then perform motif finding on the promoters. For each run of findMotifs.pl on the datasets, we let the program output at most 5 top-ranking motifs, i.e. there will be up to 5 motifs discovered by findMotifs.pl for each CEMs. To evaluate the validity of the discovered motifs, findMotifs.pl automatically compares the similarity between the discovered motif profiles and the motif profiles archived in JASPAR [18] v2018 (http://jaspar.genereg.net/) under its default parameter setting. For each discovered motif having similarity with at least one motif archived in JASPAR, we present its motif logo as well as the information of its most similar motif in JASPAR.

    All the known interactions between lncRNAs and genes are showcased in Figure 1A. The related data can be download from https://github.com/IvesG/sGavin.git data/LncReg0419 and more details are written in data/readme.txt. In figure 1A, dark-blue nodes represent LncRNAs, light-blue nodes represent proteins, pink edges represent interactions documented in Homo, green edges represent interactions documented in Mus, orange edges represent interactions documented in both Homo and Mus. Meanwhile, there are some labels on the edges, categorized based on regulatory mechanisms including PTL (post-translational regulation), TC (transcriptional regulation), PTC (post-transcriptional regulation), TL (translational regulation), and NS (not sure). The distribution above is displayed in Figure 1B and nearly three fourth (714/925) of them are identified at the transcriptional level. Other labels on the edges are categorized based on regulatory relationships including down, up, active and inactive. The distribution above is displayed in Figure 1C. The down relationships (575) are more than up relationships (308), and the proportion of active/inactive is scarce (4.5%).

    Figure 1.  Correlation Analysis of the known interactions of genes and lncRNA. (A) The network of known interactions between genes and lncRNA. The distribution of the regulatory mechanisms of the network. (C) The distribution of the regulatory relationships of the network, (D & E) The number of lncRNA that regulate genes in the scatter diagram and bar chart respectively.

    Figure 1E showed the distribution of a number of genes regulated by each lncRNA. It can be found that most lncRNA (~78%) regulate less than 5 genes. To show the specific details of the number of genes regulated by each lncRNA, Figure 1D is made, each point in the Figure 1D reflect the number of lncRNA (horizontal coordinate) that regulate certain number of genes (longitudinal coordinates) e.g. the point with coordinate (4, 14) in Figure 1D indicate that there are 14 lncRNA and each of them regulate 4 genes. The lncRNA that regulate more genes in Figure 1D belongs to the more concentrated parts in Figure 1(A).

    With the quantile-based discretization method and biclustering analysis, there are some co-expressed gene modules (CEMs) are found. The details of the way we identify CEMs are showcased in Figure 2C. Figure 2A shows the number of co-expressed gene modules(CEMs) we have got from seed1 and seed2 processed by max(min [19])-based (QUBIC1.0, [15,20]) and KL-based bi-clustering analysis (QUBIC2.0 [21]) respectively. And the distributions of numbers of genes and conditions for each CEM can be found in Figure 2D and Figure 2E. For instance, label seed1_1genes represent that Qubic1.0 is performed on seed1. To better illustrate the distributions, we have constrained the number of each gene and the number of each condition below 100 (in Figure 2B). It is found that the KL-based biclustering method tends to generate CEMs that contain fewer genes while more conditions than max(min [22])-based biclustering method from Figure 2D and Figure 2E. The difference between the results (whether the distribution of genes sizes and condition sizes or the number of CEMs that we have predicted) we get from seed1 and seed2 is subtle.

    Figure 2.  Correlation Analysis of the predicted co-expressed gene modules. (A) The number of CEMs that we have obtained. (B) The number of CEMs with size (both numbers of conditions and genes) constrained below 100. (C) The flow-process diagram describing the way we find CEMs. (D) The distribution of numbers of conditions of CEMs. (E) The distribution of the number of genes of CEMs.

    The proportion of CEMs that have significant P-values (below a pre-selected P-value cutoff) as well as proportions of the number of unique enriched lncRNA in each bicluster that belong to certain categories (i.e., number of lncRNA = 2, 3 or > 3) are calculated and shown in Figure 3A and Figure 3B. In the figures, Seed1_qubic1 represent the proportions from the results obtained using quantile discretization and using max(min [23])-based biclustering on seed 1, seed1_qubic2 represent using quantile discretization and using KL-based biclustering on seed1. In Figure 3A, it can be found that most CEMs have P-value more than 0.001 and seed1_qubic1 seems to have more significant P-value. Constrain P-value below 0.00001 and there is barely CEMs remained (less than 7%). In Figure 3B the majority of (more than70%) CEMs are with enriched lncRNA more than 3 and especially most (around 85%) of CEMs from seed1_qubic1.

    Figure 3.  Correlation Analysis of the enriched lncRNA and P-value in CEMs. (A)Proportions of CEMs that significantly enriched with lncRNAs and proportions of the number of enriched lncRNA for seed1. (B)Proportions of CEMs that significantly enriched with lncRNAs and proportions of the number of enriched lncRNA for seed2.

    For validation, we separated the HG + MG genes into two parts randomly and equally for 10 times and obtained 10 cover ratios correspondingly to check the accuracy of the previously predicted genes. The results of our validation are calculated and shown in Table 1. From Table 1 it can be found that all of the cover ratios are under 25% and the average ratio is 16.25%. We further calculated the p-value of the coverage rates. The results indicated that even the coverage date has a lot of space to be improved, the statistical significance of them are acceptable.

    Table 1.  Groups refer to genes that we extracted.
    Group 1 2 3 4 5 6 7 8 9 10
    Ratio 14.00% 19.60% 14.00% 15.70% 14.50% 21.30% 9.80% 24.30% 14.00% 15.30%
    P-value 7.58E-07 5.14E-14 7.58E-07 8.11E-09 2.56E-07 1.24E-16 5.37E-03 1.27E-21 7.58E-07 2.65E-08

     | Show Table
    DownLoad: CSV

    Since lncRNA plays an important role in regulation, they should have cooperation with transcription factor [23,24]. Thus we conduct the analysis about the DNA binding sites of related to CEMs [19,22,25]. As described in the Method section, we choose five CEMs to conducted TF motif analysis. The corresponding gene list files each containing 3, 6, 14, 20 and 21 genes. The predicted motifs and the comparison between them and JASPAR motifs are listed in Table 2, along with the function of the target TFs. In Table 2, the second column has the name of lncRNA related with this CEMs and the p-value of their correlations; The third column contains the motif consensus by Homer; The fourth column provides TF names of the most similar motifs in JASPAR, along with the similarity scores in the fifth column. These TFs may have cooperation with corresponding lncRNAs. In the first column, all the P-value of the LncRNA from the CEM is below 0.01 and the least P-value is from LncRNA HOTAIR. The supplementary table S1 with more details, including and the logo of discovered motifs and the functions of corresponding TFs, can be downloaded by visiting the GitHub link.

    Table 2.  Comparison between discovered motifs and JASPAR motifs.
    lncRNA (p-value) Homer motifs JASPAR TFs Scores
    1 FOXCUT
    3.6e-3
    AACCAVTTHDCG TFCP2 0.64
    TCCTATCACACR MEIS2 0.62
    TTTTHAAAGGGG CHR 0.67
    ARTGGTTGTWGA FOXJ2 0.58
    GCAATCTCGC IRF4 0.66
    2 ANCR
    1.1e-3
    AGGGTGACAG SPZ1 0.80
    GGTATCTTAC GATA5 0.64
    CTCATAGGAG GCM1 0.65
    TAAGTGAAAG PRDM1 0.86
    CTTTTGGAAC CHR 0.65
    3 250-280
    2.2e-4
    WYTRTCTTTGCG RXR 0.61
    TCTTACGG ELK1 0.71
    GGCAAGGA SD 0.76
    GAGGTATGTT TEAD1 0.70
    TGCCGGGAGCGT POL 0.64
    4 HOXD-AS1
    6.1e-3
    CTCGAGTAGG PB0114 0.63
    GCCCCCTGCA PB0076 0.74
    ACGYMYATKYCC GFY 0.59
    AGCGGGTT PH 0.68
    AGGCGCCGCGCC SP1 0.69
    5 HOTAIR
    5e-6
    TGGCGCAGCGCG PB 0.67
    GTACAACTTT PB 0.66
    CMTSTGTCWCYK NeuroG2 0.66
    GTGATCCATT RHOXF1 0.68
    GGTMGRRGTGMW TBX20 0.58

     | Show Table
    DownLoad: CSV

    In order to further evaluate the biological significance of the identified CEMs, we tested the enrichment of the genes in each CEM in Gene ontology terms and KEGG pathways using clusterProfiler package of R project BioConductor under q-value cutoff 0.05, of which the description of the GO terms and KEGG pathways that the CEMs are enriched in are presented in Table 3 and Table 4 respectively. And the supplementary table S2 with more details, including original and adjusted P-value, proportion of the matched genes, gene's ID, etc., can be downloaded on GitHub link.

    Table 3.  Gene ontology information of selected CEMs.
    LncRNA ID Description q-value
    FOXCUT GO:0033613 transmembrane receptor protein tyrosine kinase activity 1.2937E-02
    GO:0033613 activating transcription factor binding 1.2937E-02
    GO:0019199 transmembrane receptor protein kinase activity 1.2937E-02
    GO:0001085 RNA polymerase Ⅱ transcription factor binding 2.0767E-02
    250-280 GO:0003735 structural constituent of ribosome 3.1600E-06
    GO:0003729 mRNA binding 1.1140E-03
    GO:0008483 transaminase activity 1.1140E-03
    GO:0048027 mRNA 5'-UTR binding 1.1140E-03
    GO:0016769 transferase activity, transferring nitrogenous groups 1.1140E-03
    GO:0045182 translation regulator activity 1.5826E-03
    GO:0030170 pyridoxal phosphate binding 1.5826E-03
    GO:0070279 vitamin B6 binding 1.5826E-03
    GO:0019843 rRNA binding 1.5826E-03
    GO:0019842 vitamin binding 3.1903E-03
    HOXD-AS1 GO:0004714 transmembrane receptor protein tyrosine kinase activity 1.2937E-02
    GO:0033613 activating transcription factor binding 1.2937E-02
    GO:0019199 transmembrane receptor protein kinase activity 1.2937E-02
    GO:0001085 RNA polymerase Ⅱ transcription factor binding 2.0767E-02
    HOTAIR GO:0005109 frizzled binding 1.9097E-03
    GO:0001227 transcriptional repressor activity, RNA polymerase Ⅱ transcription regulatory region sequence-specific binding 1.9097E-03
    GO:0001664 G-protein coupled receptor binding 1.9686E-03
    GO:0001078 transcriptional repressor activity, RNA polymerase Ⅱ core promoter proximal region sequence-specific binding 1.1196E-02
    GO:0008201 heparin binding 1.2885E-02
    GO:0005539 glycosaminoglycan binding 1.8790E-02
    GO:1901681 sulfur compound binding 1.9016E-02
    GO:0045236 CXCR chemokine receptor binding 2.1785E-02
    GO:0008301 DNA binding, bending 2.1785E-02
    GO:0001223 transcription coactivator binding 2.1785E-02
    GO:0042813 Wnt-activated receptor activity 2.1785E-02
    GO:0035198 miRNA binding 2.2807E-02
    GO:0017147 Wnt-protein binding 2.5258E-02
    GO:1990841 promoter-specific chromatin binding 2.5258E-02
    GO:0000982 transcription factor activity, RNA polymerase Ⅱ core promoter proximal region sequence-specific binding 2.5258E-02
    GO:0001221 transcription cofactor binding 2.5587E-02

     | Show Table
    DownLoad: CSV
    Table 4.  KEGG pathway information of selected CEMs.
    LncRNA ID Description q-value
    FOXCUT hsa05216 Thyroid cancer 8.8587E-03
    hsa04510 Focal adhesion 8.8587E-03
    hsa05205 Proteoglycans in cancer 8.8587E-03
    hsa05218 Melanoma 1.7683E-02
    hsa05214 Glioma 1.7683E-02
    hsa04151 PI3K-Akt signaling pathway 1.8745E-02
    hsa05215 Prostate cancer 1.8745E-02
    hsa01522 Endocrine resistance 1.8745E-02
    hsa04919 Thyroid hormone signaling pathway 2.2317E-02
    hsa04152 AMPK signaling pathway 2.2317E-02
    hsa04068 FoxO signaling pathway 2.4442E-02
    hsa04550 Signaling pathways regulating pluripotency of stem cells 2.5129E-02
    hsa05224 Breast cancer 2.5129E-02
    hsa04218 Cellular senescence 2.7471E-02
    hsa05225 Hepatocellular carcinoma 2.7471E-02
    hsa04530 Tight junction 2.7471E-02
    250-280 hsa03010 Ribosome 2.9900E-05
    hsa01210 2-Oxocarboxylic acid metabolism 3.7322E-03
    hsa00220 Arginine biosynthesis 3.7322E-03
    hsa00250 Alanine, aspartate and glutamate metabolism 4.7849E-03
    hsa01230 Biosynthesis of amino acids 7.9156E-03
    HOXD-AS1 hsa05216 Thyroid cancer 8.8587E-03
    hsa04510 Focal adhesion 8.8587E-03
    hsa05205 Proteoglycans in cancer 8.8587E-03
    hsa05218 Melanoma 1.7683E-02
    hsa05214 Glioma 1.7683E-02
    hsa04151 PI3K-Akt signaling pathway 1.8745E-02
    hsa05215 Prostate cancer 1.8745E-02
    hsa01522 Endocrine resistance 1.8745E-02
    hsa04919 Thyroid hormone signaling pathway 2.2317E-02
    hsa04152 AMPK signaling pathway 2.2317E-02
    hsa04068 FoxO signaling pathway 2.4442E-02
    hsa04550 Signaling pathways regulating pluripotency of stem cells 2.5129E-02
    hsa05224 Breast cancer 2.5506E-02
    hsa04218 Cellular senescence 2.7471E-02
    hsa05225 Hepatocellular carcinoma 2.7471E-02
    hsa04530 Tight junction 2.7471E-02
    HOTAIR hsa04310 Wnt signaling pathway 5.7295E-03

     | Show Table
    DownLoad: CSV

    Within this study, we have developed a method for elucidating lncRNA-gene and transcription factor-lncRNA interactions using a biclustering approach. The method was performed on 2 breast cancer RNA-seq datasets from TCGA. The bicluster method allows for the identification of particular expression patterns across multiple datasets, indicating networks of lncRNA and gene interactions. The developed method will also provide a way for future lncRNA interaction studies. Certainly, the predict performance still far from satisfactory, which is not unexpected since we only used RNA-Seq data. Actually, the interaction mechanism between lncRNA and genes are far more complex, and more data should be involved if we want to capture the whole picture of them. We are planning to include some other data, like proteomics and chromatin accessibility information, to improve the prediction. Besides, the evaluation on the relationship between lncRNA and predicted CEMs also has the potential to be improved, e.g. calculating the adjusted P-value or overall P-value in place of the original P-values used in this study. In view of the application, we will work on more specific examples of the regulatory functions of some particular lncRNAs and identify some hypothesized mechanisms of these regulatory functions. Also, the further analysis of the difference of lncRNA related genes between tumor and normal samples could provide more information for studying the process and mechanism of cancer occurrence and development, e.g. determination of the stage of developed tumors, which will be our concern in the future research.

    This work was supported by the National Nature Science Foundation of China (NSFC) [61772313 and 61432010], Young Scholars Program of Shandong University [YSPSDU, 2015WLJH19], the Innovation Method Fund of China (2018IM020200), and Shanghai Municipal Science and Technology Major Project (2018SHZDZX01) and ZHANGJIANGLAB. Qin Ma's work was supported by an R01 Award from the National Institute of General Medical Sciences of the National Institutes of Health [GM131399-01]. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by the National Science Foundation [ACI-1548562]. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health and the National Science Foundation.

    All authors declare no conflicts of interest in this paper.



    [1] S. Abbas, M. Benchohra, N. Hamidi, et al. Caputo-Hadamard fractional differential equations in Banach spaces, Fract. Calc. Appl. Anal., 21 (2018), 1027-1045. doi: 10.1515/fca-2018-0056
    [2] G. Adomian, G. E. Adomian, Cellular systems and aging models, Comput. Math. Appl., 11 (1985), 283-291. doi: 10.1016/0898-1221(85)90153-1
    [3] R. P. Agarwal, D. O'Regan, Toplogical degree theory and its applications, Tylor and Francis, 2006.
    [4] B. Ahmad and A. Alsaedi, Nonlinear fractional differential equations with nonlocal fractional integro-differential boundary conditions, Bound. Value Probl., 2012 (2012), 124.
    [5] B. Ahmad, A. Alsaedi and S. K. Ntouyas, Nonlinear Langevin equations and inclusions involving mixed fractional order derivatives and variable coefficient with fractional nonlocal-terminal conditions, AIMS Mathematics, 4 (2019), 626-647. doi: 10.3934/math.2019.3.626
    [6] B. Ahmad, M. Alghanmi, S. K. Ntouyas, et al. A study of fractional differential equations and inclusions involving generalized Caputo-type derivative equipped with generalized fractional integral boundary conditions, AIMS Mathematics, 4 (2019), 26-42.
    [7] A. Ali, B. Samet, K. Shah, et al. Existence and stability of solution to a toppled systems of differential equations of non-integer order, Bound. Value Probl., 2017 (2017), 16.
    [8] N. Ali, K. Shah, D. Baleanu, et al. Study of a class of arbitrary order differential equations by a coincidence degree method, Bound. Value Probl., 2017 (2017), 111.
    [9] A. Ardjouni, A. Djoudi, Positive solutions for nonlinear Caputo-Hadamard fractional differential equations with integral boundary conditions, Open J. Math. Anal., 3 (2019), 62-69.
    [10] A. Ardjouni, Positive solutions for nonlinear Hadamard fractional differential equations with integral boundary conditions, AIMS Mathematics, 4 (2019) 1101-1113.
    [11] Y. Arioua, N. Benhamidouche, Boundary value problem for Caputo-Hadamard fractional differential equations, Surveys in Mathematics and its Applications, 12 (2017), 103-115.
    [12] M. B. Zada, K. Shah and R. A. Khan, Existence theory to a coupled system of higher order fractional hybrid differential equations by topological degree theory, Int. J. Appl. Comput. Math., 4 (2018), 102.
    [13] P. W. Bates, On some nonlocal evolution equations arising in materials science. Nonlinear dynamics and evolution equations, Fields Inst. Commun., 48 (2006), 13-52.
    [14] M. Benchohra, S. Hamani, S. K. Ntouyas, Boundary value problems for differential equations with fractional order and nonlocal conditions, Nonlinear Anal., 71 (2009), 2391-2396. doi: 10.1016/j.na.2009.01.073
    [15] M. Benchohra, S. Bouriah, J. J. Nieto, Existence of periodic solutions for nonlinear implicit Hadamard's fractional differential equations, Rev. R. Acad. Cienc. Exactas Fís. Nat., Ser. A Mat., 112 (2018), 25-35.
    [16] M. Benchohra, Z. Bouteffal, J. Henderson, et al. Measure of noncompactness and fractional integro-differential equations with state-dependent nonlocal conditions in Fréchet spaces, AIMS Mathematics, 5 (2020), 15-25. doi: 10.3934/math.2020002
    [17] W. Benhamida, J. R. Graef, S. Hamani, Boundary value problems for Hadamard fractional differential equations with nonlocal multi-point boundary conditions, Fract. Differ. Calc., 8 (2018), 165-176.
    [18] W. Benhamida, S. Hamani, J. Henderson, Boundary Value Problems For Caputo-Hadamard Fractional Differential Equations, Adv. Theory Nonlinear Anal. Appl., 2 (2018), 138-145.
    [19] P. L. Butzer, A. A. Kilbas, J. J. Trujillo, Composition of Hadamard-type fractional integration operators and the semigroup property, J. Math. Anal. Appl., 269 (2002), 387-400. doi: 10.1016/S0022-247X(02)00049-5
    [20] P. L. Butzer, A. A. Kilbas, J. J. Trujillo, Fractional calculus in the Mellin setting and Hadamardtype fractional integrals, J. Math. Anal. Appl., 269 (2002), 1-27. doi: 10.1016/S0022-247X(02)00001-X
    [21] P. L. Butzer, A. A. Kilbas, J. J. Trujillo, Mellin transform analysis and integration by parts for Hadamard-type fractional integrals, J. Math. Anal. Appl., 270 (2002), 1-15. doi: 10.1016/S0022-247X(02)00066-5
    [22] A. Bednarz, L. Byszewski, On abstract nonlocal Cauchy problem, Czasopismo Techniczne, (2015), 11-17.
    [23] L. Byszewski, Theorem about existence and uniqueness of continuous solutions of nonlocal problem for nonlinear hyperbolic equation, Appl. Anal., 12 (1991), 173-180.
    [24] L. Byszewski, Existence and uniqueness of a classical solution to a functional-differential abstract nonlocal Cauchy problem, J. Appl. Math. Stochastic Anal., 12 (1999), 91-97. doi: 10.1155/S1048953399000088
    [25] D. Chergui, T. E. Oussaeif, M. Ahcene, Existence and uniqueness of solutions for nonlinear fractional differential equations depending on lower-order derivative with non-separated type integral boundary conditions, AIMS Mathematics, 4 (2019), 112-133. doi: 10.3934/Math.2019.1.112
    [26] K. Deimling, Nonlinear Functional Analysis, Springer, Berlin, Germany, 1985.
    [27] Y. Y. Gambo, F. Jarad, D. Baleanu, et al. On Caputo modification of the Hadamard fractional derivatives, Adv. Differ. Equ., 2014 (2014), 10.
    [28] J. R. Graef, N. Guerraiche and S. Hamani, Boundary value problems for fractional differential inclusions with Hadamard type derivatives in Banach spaces, Stud. Univ. Babes-Bolyai Math., 62 (2017), 427-438. doi: 10.24193/subbmath.2017.4.02
    [29] A. Guezane-Lakoud, R. Khaldi, Solvability of a fractional boundary value problem with integral condition, Nonlinear Analysis, 75 (2012), 2692-2700. doi: 10.1016/j.na.2011.11.014
    [30] J. Hadamard, Essai sur l'etude des fonctions donnees par leur developpment de Taylor, J. Mat. Pure Appl. Ser., 8 (1892), 101-186.
    [31] R. Hilfer, Application of fractional calculus in physics, New Jersey: World Scientific, 2001.
    [32] F. Isaia, On a nonlinear integral equation without compactness, Acta Math. Univ. Comenian. (N.S.), 75 (2006), 233-240.
    [33] F. Jarad, D. Baleanu, A. Abdeljawad, Caputo-type modification of the Hadamard fractional derivatives, Adv. Differ. Equ., 2012 (2012), 142.
    [34] R. A. Khan and K. Shah, Existence and uniqueness of solutions to fractional order multi-point boundary value problems, Communications in Applied Analysis, 19 (2015), 515-526.
    [35] A. A. Kilbas, Hadamard-type fractional calculus, J. Korean Math. Soc., 38 (2001), 1191-1204.
    [36] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and Applications of Fractional Differential Equations, North-Holland Mathematics Sudies, Elsevier Science, 2006.
    [37] Z. Laadjal, Q.-H. Ma, Existence and uniqueness of solutions for nonlinear Volterra-Fredholm integro-differential equation of fractional order with boundary conditions, Math. Meth. Appl. Sci., 2019.
    [38] Q. Ma, R. Wang, J. Wang, et al. Qualitative analysis for solutions of a certain more generalized two-dimensional fractional differential system with Hadamard derivative, Appl. Math. Comput., 257 (2015), 436-445.
    [39] K. S. Miller, B. Ross, An Introduction to Fractional Calculus and Fractional Differential Equations, wiley, New YorK, 1993.
    [40] K. B. Oldham, Fractional differential equations in electrochemistry, Adv. Eng. Softw., 41 (2010), 9-12. doi: 10.1016/j.advengsoft.2008.12.012
    [41] I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, 1993.
    [42] J. Sabatier, O. P. Agrawal, J. A. T. Machado, Advances in Fractional Calculus-Theoretical Developments and Applications in Physics and Engineering, Dordrecht: Springer, 2007.
    [43] K. Shah, A. Ali and R. A. Khan, Degree theory and existence of positive solutions to coupled systems of multi-point boundary value problems, Bound. Value Probl., 2016 (2016), 43.
    [44] K. Shah and R. A. Khan, Existence and uniqueness results to a coupled system of fractional order boundary value problems by topological degree theory, Numer. Funct. Anal. Optim., 37 (2016), 887-899. doi: 10.1080/01630563.2016.1177547
    [45] K. Shah, W. Hussain, P. Thounthong, et al. On nonlinear implicit fractional differential equations with integral boundary condition involving p-Laplacian operator without compactness, Thai J. Math., 16 (2018), 301-321.
    [46] K. Shah and W. Hussain, Investigating a class of nonlinear fractional differential equations and its Hyers-Ulam stability by means of topological degree theory, Numer. Funct. Anal. Optim., 40 (2019), 1355-1372. doi: 10.1080/01630563.2019.1604545
    [47] M. Shoaib, K. Shah, R. Ali Khan, Existence and uniqueness of solutions for coupled system of fractional differential equation by means of topological degree method, Journal Nonlinear Analysis and Application, 2018 (2018), 124-135.
    [48] V. E. Tarasov, Fractional Dynamics: Application of Fractional Calculus to Dynamics of Particles, Fields and Media, Springer, Heidelberg & Higher Education Press, Beijing, 2010.
    [49] J. Tariboon, A. Cuntavepanit, S. K. Ntouyas, et al. Separated boundary value problems of sequential Caputo and Hadamard fractional differential equations, J. Funct. Space. Appl., 2018 (2018), 1-8.
    [50] J. Wang, Y. Zhou and W. Wei, Study in fractional differential equations by means of topological degree methods, Numer. Funct. Anal. Optim., 33 (2012), 216-238. doi: 10.1080/01630563.2011.631069
  • This article has been cited by:

    1. Gerardo Chowell, Amna Tariq, James M. Hyman, A novel sub-epidemic modeling framework for short-term forecasting epidemic waves, 2019, 17, 1741-7015, 10.1186/s12916-019-1406-6
    2. Ping Yan, Gerardo Chowell, 2019, Chapter 6, 978-3-030-21922-2, 183, 10.1007/978-3-030-21923-9_6
    3. Ping Yan, Gerardo Chowell, 2019, Chapter 9, 978-3-030-21922-2, 317, 10.1007/978-3-030-21923-9_9
    4. Christiana Tsiligianni, Aristeides Tsiligiannis, Christos Tsiliyannis, A stochastic inventory model of COVID-19 and robust, real-time identification of carriers at large and infection rate via asymptotic laws, 2023, 304, 03772217, 42, 10.1016/j.ejor.2021.12.037
    5. Ping Yan, Gerardo Chowell, Modeling sub-exponential epidemic growth dynamics through unobserved individual heterogeneity: a frailty model approach, 2024, 21, 1551-0018, 7278, 10.3934/mbe.2024321
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4851) PDF downloads(582) Cited by(18)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog