Asthma is a chronic inflammatory disease of the respiratory tract, the main characteristics of which are the hypersensitivity of the airways to various stimuli and reversible airflow obstruction . More than 300 million patients worldwide are suffering from this disease . Another widespread disease in the world is hypertension, which is characterized by a persistent increase in blood pressure from 140/90 mmHg. Art. and higher . Annually about nine million people die due to complications associated with hypertension . A number of studies have shown that asthma is related to an increase in the incidence of cardiovascular diseases , . In this regard, the combination of asthma and arterial hypertension (comorbidity) in a large percentage of clinical cases is not surprising. Prevalence of hypertension in asthma patients is high with odds ratio = 1.66 , , p-value < 0.00001 . This may not be a simple coincidence, because standard asthma mechanisms have been implicated in key processes of initiation of hypertension. Good examples of this are the establishment of the role of arachidonic acid-leukotriene B4 production in spontaneously hypertensive rat  or the value of Th17 and IL17 in arterial hypertension . Therapeutic control of allergic inflammation leads to better control of blood pressure . The important role of the transcription factor STAT3, which is involved in signal transduction for a variety of cytokines and takes an active part in allergic inflammation ,  and vascular remodeling , is well-known. On the basis on this information, it can be assumed that the imbalance between the pro-inflammatory and anti-inflammatory factors in the body creates the optimal situation for the realization of a hereditary predisposition to both asthma and arterial hypertension. In addition to the important role played by immune responses and control of inflammation in the pathogenesis of asthma and hypertension, it is also possible to assume the involvement of other mechanisms underlying the comorbidity of these diseases. For example, β-adrenoblockers and ACE inhibitors, which for a long time were contraindicated for patients with bronchial asthma due to possible bronchospasm, are widely used in the treatment of hypertension. Polymorphism of β-adrenoreceptor genes is associated with the risk of hypertension and bronchial asthma , . It is assumed that mutations in the SLC26A4 gene can play an important role in the pathogenesis of bronchial asthma and hypertension, and therefore in the comorbidity of these diseases , , . However, despite the accumulated data on the participation of specific genes simultaneously in asthma and hypertension, the mechanisms of interaction of these two diseases and the formation of their comorbid state are still poorly understood. Previously, we performed the prioritization of genes that are associated with asthma and hypertension simultaneously, to find the candidates which can potentially be involved in the mechanisms of comorbidity of these two diseases .
In the present work, a prioritization of an extended set of genes, which are associated either with asthma or with hypertension, was carried out. In addition, a new criterion for prioritization, based on the assessment of cross-talk centrality, previously used by us for the prioritization of neuronal apoptosis genes in the associative gene network of Parkinson’s disease , was introduced. When assessing the quality of prioritization, this criterion was of the best significance. Among the most significant genes were cytokines (CXCL10, CD40LG, IL1A, IL13 and IL15), toll-like receptor TLR2, death domain receptor FAS, transcription factor STAT3, complement factor C3 and amyloid beta precursor protein APP. An analysis of the overrepresentation of Gene Ontology biological processes for the top 100 highest-priority genes has shown the potential importance of the response to hypoxia, the inflammatory and immune responses and regulation of the apoptotic process in the mechanism of interaction of asthma and hypertension. In addition, we analyzed the dynamics of interest in genes by the frequency of their mention in scientific publications. It appeared that genes simultaneously associated with asthma and hypertension are enriched with genes that have a high positive dynamic of interest in them from researchers. This fact may indicate the involvement of these genes in many mechanisms of manifestation of various phenotypic traits that can determine the importance of their role in the development of comorbidity of diseases. An analysis of the most priority genes also showed significant enrichment with genes that have a high dynamic of interest in them. This may indicate their involvement in multiple phenotypic signs and, accordingly, in the comorbidity of asthma and hypertension. The genes identified as a result of prioritization can be used in the planning of experiments on genotyping and search for associations with the comorbid state of asthma and hypertension.
The sets of genes associated with asthma and hypertension were extracted from the ANDSystem knowledge base , . The set of genes associated with asthma included 749 genes. The set of genes associated with hypertension included 706 genes. The number of genes simultaneously presented in the first and second sets was 202. A test set of genes, used for prioritization, was obtained by combining sets of asthma and hypertension genes, minus the 202 genes gained by the intersection of these two sets. The total number of genes in the test set was 1051.
The over-represented Gene Ontology biological processes were identified using the DAVID 6.8 service  with the following parameters: the organism – "Homo sapiens", Gene_Ontology – "GOTERM_BP_DIRECT".
To prioritize the genes we used four criteria, which we previously discussed in the analysis of asthma, hypertension and Parkinson’s disease , . Criteria 1 and 2 were based on the well-known methods of prioritization: ToppGene  (https://toppgene.cchmc.org/prioritization.jsp) and Endeavor  version 3.71 (https://endeavour.esat.kuleuven.be/), respectively. The ToppGene and Endeavor programs allow one to perform a ranking of a test set of genes by a training set of genes, according to specific criteria characterizing the proximity of genes from the test set to the genes from the training sample. The genetic information (co-localization in the genome), the functional properties of genes (participation in the same GO categories), etc., can be used by these programs as properties of the vertices of the graph of protein-protein networks. As a training set for each of these programs, a list of genes associated with asthma and hypertension, described above, was provided, while a list of 1051 genes, which is also described above, was used as a test set.
For the ToppGene in the Training parameters section, all features were selected. The "Rank" indicator of the ToppGene program was used to rank the genes. In the case of the Endeavor system, all settings were set by default. The ranking of genes was carried out according to the "p-value" indicator of the Endeavor program. The lowest rank received genes with the lowest values of "p-value", and the largest rank – genes with the highest values of "p-value".
Criterion 3 reflects the specificity of association of biological processes in which the analyzed gene participates with comorbidity mechanisms. Criterion 3 was calculated by formula N1/N2, where N1 is a number of the Gene Ontology biological processes, overrepresented for a set of genes associated with asthma and hypertension; N2 is a number of all Gene Ontology biological processes in which the analyzed gene was involved. Information indicating that the gene is involved in the Gene Ontology biological process was obtained from the AmiGO system . Ranks for genes were determined by sorting the list of genes by descending the number of over-represented Gene Ontology biological processes. Thus, the lowest rank was attributed to the gene for which the number of the over-represented Gene Ontology biological processes was the highest.
To calculate criterion 4, an associative gene network which describes the interactions between the analyzed gene/protein and the genes/proteins from the training set was used. The network was automatically reconstructed by ANDSystem. Criterion 4 was based on the use of the cross-talk centrality (CTC), calculated using the “Intelligent Filtration” function of the “ANDVisio” program. Within this criterion, the CTC value was calculated separately for vertices in the gene network corresponding to the gene and its product (protein). Thus, for the gene, two indices (CTCgene and CTCprotein) were determined.
The CTCgene was calculated by the formula:
where Ni – is a number of links of the i-th gene with members of the associative gene network of asthma/hypertension; M – is a number of vertices of the associative gene network of asthma/hypertension. CTCprotein for the i-th protein was determined the same way. The CTCgene and CTCprotein grades were determined by sorting the values in descending order, like in the case of criterion 3.
In the case when several genes had the same value for the criterion under which the ranking was conducted, the rank of these genes was calculated the similar way as the Spearman rank correlation coefficient . Such genes were assigned the same rank, equal to the average arithmetic ranks of these genes, according to their position, in the sorted list of genes.
As the fifth criterion, the average value of the ranks obtained according to the criteria 1–4 was taken like it was done in .
For the reconstruction of the dynamics, the information on investigated genes from PubMed publications for the period from 1990 to 2017 were considered. For each year from this period, the number of publications in which the analyzed gene was mentioned was calculated. The marking of gene names in the texts of publications was carried out by using the ANDSystem. The trend reflecting the dynamics of the number of gene references in the articles was calculated using the Pearson correlation coefficient between the year and the total number of articles published this year in which the gene was mentioned. The enrichment of the gene set with genes with increased interest dynamics was calculated using the hypergeometric distribution realized in the hypergeom.sf function of the scipy.stats package of the Python programming language . Genes with high dynamics of interest were genes for which the Pearson correlation coefficient exceeded 0.9.
3 Results and Discussion
To prioritize the genes for their potential simultaneous involvement in asthma and hypertension 749 genes of asthma and 706 genes of hypertension were considered. The training set of genes consisted of 202 genes, which were presented both in asthma and hypertension lists. Prioritization was carried out by using 4 criteria (see Online Appendix 1). The ranking of genes was carried out according to the value calculated as the mean for the ranks of the genes obtained using each of the 4 criteria. The top 10 of the most significant genes is listed in Table 1. Among them were cytokines (CXCL10, CD40LG, IL1A, IL13 and IL15), toll-like receptor TLR2, death domain receptor FAS, transcription factor STAT3, complement component C3 and amyloid beta precursor protein APP. It is interesting to note that all these genes are genes from the group of asthma. At the same time, the ratio of asthma genes to the genes of hypertension in the top 100 of the most significant genes appeared to be 2:1.
Figure 1 shows that the top 10 genes are relatively well connected to each other inside the associative gene network, reconstructed with ANDSystem (see Online Appendix 2). In particular, the transcription factor STAT3 regulates CXCL10  and Fas . In turn, TLR2, IL-13 and IL-15 regulate the activity of STAT3 , , .
For better understanding which biological processes can underlie the molecular mechanisms of interaction between asthma and hypertension, an analysis of the overrepresentation of Gene Ontology biological processes on two sets of data was conducted: the training set of genes and the top 100 of the highest priority genes from the table in Online Appendix 1 were analyzed. For the set of genes associated with asthma and hypertension simultaneously, there were found 69 significantly over-represented biological processes (p-value < 0.05 with FDR correction). For the top 100 of the most priority genes, there were 42 such processes. Twenty-two biological processes were over-represented simultaneously for these two sets of genes (see Online Appendix 3). For example, the response to hypoxia had the best relevance for the training set of genes (p-value < 10−12). In the second set of genes, the p-value for this process was <10−8, and the most significant process appeared to be the inflammatory response (p-value < 10−19). The intersection of over-represented processes was more than 50% of the total number of all overrepresented processes for a set of top 100 priority genes. Such a high value was partly related to the 3rd criterion of prioritization, which was based on considering the overrepresented biological processes. But since the ranking was performed by using the rank averaged over all 4 criteria, such intersection indirectly confirms the correctness of the prioritization.
It should be noted that the manual analysis of the full-text articles of PubMed Central showed that for genes from the top 10 there was additional information that was not represented in ANDSystem that speaks in favor of their possible connection with hypertension. For example, for the IL1A gene, it is known that in Chinese Han population genetic variations in IL1A are associated with the risk of preeclampsia, which is characterized by hypertension and proteinuria . In addition in Chinese children with type 1 diabetes the prevalence of hypertension in carriers of T allele in the IL1A gene was 4.2-fold higher than for C allele carriers with p-value < 0.001 . The article  stated that the down-regulation of CD40LG might be related to kidney dysfunction in an early phase of hypertension, and  discussed that CD40LG was included in pulmonary hypertension gene subnetwork. Gene STAT3 regulate angiogenesis and vasoconstriction and is involved in the development of portal hypertension . In patients with essential hypertension, IL-15 might be involved in the development of cardiovascular diseases . For gene FAS it was shown that its expression was significantly increased in spontaneously hypertensive rats compared to Wistar-Kyoto rats . In  APP was identified to undergo differential ubiquitination during pulmonary hypertension pathogenesis. Hypomethylation of the TLR2 gene increased the risk of essential hypertension in the Chinese population . It was shown that C3 is highly expressed in mesenchymal tissues in spontaneously hypertensive rats (SHR)  and increased C3 induces salt-sensitive hypertension in SHR . In  was discussed that IL-13 was involved in the induction of pulmonary hypertension. CXCL10 appeared to be involved in the pathogenesis of hypertension  and its level in serum was elevated in untreated essential hypertension .
One of the interesting characteristics of a gene may be the dynamics of increasing researcher’s interest in this gene, which may reflect the prospect of its study. It turned out that the genes from the training set, simultaneously associated with asthma and hypertension, are statistically significantly enriched with genes that have a high positive dynamic of their mention in the scientific publications of PubMed (p-value < 10−15). The total number of genes with a positive growing trend (r > 0.9) was 872 of the 17,868 human genes, presented in the ANDSystem. Among the 202 genes of the training set, it appeared to be 20% of such genes (42 genes). The HMOX1 gene, which encodes heme oxygenase 1, possessed the greatest positive dynamics of interest (r = 0.99). It is interesting to note that the top 100 of most priority genes contained 25% of genes possessing increased positive dynamics of interest. This enrichment was also statistically significant (p-value < 10−11). Analysis of the top 5 genes showed that among them STAT3 belongs to the group of 25% of genes possessing increased positive dynamics of interest (r = 0.96). The differences in the dynamics of interest for genes from the top 5 group are shown in Figure 2.
Gene STAT3 encodes a protein that in response to cytokines and growth factors is phosphorylated and translocated to the nucleus where it starts to act as a transcription activator . STAT3 is involved in many cellular processes including cell growth and apoptosis .
4 Concluding Remarks
It is well known that many diseases are comorbid, which significantly complicates their treatment, and also prevents their proper diagnosis. The study of the molecular mechanisms of comorbid states is an urgent task of the modern biomedicine. Asthma and hypertension are examples of widespread comorbid diseases. In this study, we used prioritization methods, including existing well-known methods (ToppGene and Endeavor), and an approach (cross-talk centrality), based on the analysis of associative gene networks from ANDSystem, to find genes that could potentially be involved in asthma and hypertension simultaneously. Cytokines (CXCL10, CD40LG, IL1A, IL13 and IL15), toll-like receptor TLR2, death domain receptor FAS, transcription factor STAT3, complement factor C3 and amyloid beta precursor protein (APP), appeared to be the most significant. In addition, an analysis of the dynamics of the frequency of mentioning genes in scientific publications was carried out. It turned out that the top 100 of the priority genes are significantly enriched by genes with increased positive dynamics of researcher’s interest in them. This may indicate a promise for further researches on these genes. A further study of the predicted genes, simultaneously involved in asthma and hypertension, can help to illuminate the mechanisms of the comorbid state of these two diseases.
Gene prioritization was carried out with the support of the Volkswagen Stiftung Trilateral Partnerships – Cooperation Projects between Scholars and Scientists from Ukraine, Russia, and Germany “In silico screening and experimental validation of new drug targets for the treatment of co-morbid multifactorial diseases” (90.335).
Development of methods and analysis of the dynamics of interest in the area of genomics and medical genetics on the basis of frequency of mentioning of genes in scientific publications was supported by Federal target program “Research and development in priority areas of development of Russia’s scientific and technological complex for 2014–2020”. Agreement on granting a subsidy as of October 23, 2017 No. 14.601.21.0015 between the Ministry of Education and Science of the Russian Federation and Tomsk NRMC on conducting research project titled: Development of a forecast for the implementation of the priority of scientific and technological development defined in paragraph 20c “personalized medicine, high-tech health care and health saving technologies, including the rational use of medicines (primarily antibacterial)” Scientific and technological development strategy of the Russian Federation. The unique identifier of the project is RFMEFI60117X0015. The state agreement identifier is 0000000007417PE10002.
Camm AJ, Lüscher TF, Serruys PW, eds. The ESC textbook of cardiovascular medicine. Oxford, UK: OXFORD University Press, 2009:1136. Google Scholar
Prosser R, Carleton B, Smith A. The comorbidity burden of the treated asthma patient population in British Columbia. Chronic Dis Inj Can 2010;30(2):46–55. Google Scholar
Tattersall MC, Guo M, Korcarz CE, Gepner AD, Kaufman JD, Liu KJ, et al. Asthma predicts cardiovascular disease events: the multi-ethnic study of atherosclerosis. Arterioscler Thromb Vasc Bio 2015;35(6):1520–5. CrossrefGoogle Scholar
Waki H, Hendy EB, Hindmarch CC, Gouraud S, Toward M, Kasparov S, et al. Excessive leukotriene B4 in nucleus tractus solitarii is prohypertensive in spontaneously hypertensive rats. Hypertension 2013;61(1):194–201. PubMedWeb of ScienceCrossrefGoogle Scholar
Litonjua AA, Tantisira KG, Lake S, Lazarus R, Richter BG, Gabriel S, et al. Polymorphisms in signal transducer and activator of transcription 3 and lung function in asthma. Respir Res 2005;6(1):52. CrossrefPubMedGoogle Scholar
Simeone-Penney MC, Severgnini M, Tu P, Homer RJ, Mariani TJ, Cohn L, et al. Airway epithelial STAT3 is required for allergic inflammation in a murine model of asthma. J Immunol 2007;178(10):6191–9. Web of ScienceCrossrefGoogle Scholar
Liu Z, Qi H, Liu B, Liu K, Wu J, Cao H, et al. Genetic susceptibility to salt-sensitive hypertension in a Han Chinese population: a validation study of candidate genes. Hypertens Res 2017;40(10):876–84. CrossrefWeb of ScienceGoogle Scholar
Dossena S, Bizhanova A, Nofziger C, Bernardinelli E, Ramsauer J, Kopp P, et al. Identification of allelic variants of pendrin (SLC26A4) with loss and gain of function. Cell Physiol Biochem 2011;28(3):467–76. CrossrefPubMedWeb of ScienceGoogle Scholar
Nofziger C, Vezzoli V, Dossena S, Schönherr T, Studnicka J, Nofziger J, et al. STAT6 links IL‐4/IL‐13 stimulation with pendrin expression in asthma and chronic obstructive pulmonary disease. Clin Pharmacol Ther 2011;90(3):399–405. CrossrefPubMedWeb of ScienceGoogle Scholar
Saik OV, Demenkov PS, Ivanisenko TV, Bragina EY, Freidin MB, Goncharova IA, et al. Novel candidate genes important for asthma and hypertension comorbidity revealed from associative gene networks. BMC Med Genomics 2018;11(1):15. Web of SciencePubMedCrossrefGoogle Scholar
Yankina MA, Saik O, Demenkov PS, Khusnutdinova EK, Rogaev E, Lavrik IN, et al. Analysis of the interactions of neuronal apoptosis genes in the associative gene network of Parkinson’s disease. Vavilovskii Zhurnal Genet Selektsii 2018;22(1):153–60. Web of ScienceGoogle Scholar
Demenkov PS, Ivanisenko TV, Kolchanov NA, Ivanisenko VA. ANDVisio: a new tool for graphic visualization and analysis of literature mined associative gene networks in the ANDSystem. In Silico Biol 2012;11(3,4):149–61. Google Scholar
Ivanisenko VA, Saik OV, Ivanisenko NV, Tiys ES, Ivanisenko TV, Demenkov PS, et al. ANDSystem: an Associative Network Discovery System for automated literature mining in the field of biology. BMC Syst Biol 2015;9(2):S2. PubMedCrossrefWeb of ScienceGoogle Scholar
Tranchevent LC, Ardeshirdavani A, ElShal S, Alcaide D, Aerts J, Auboeuf D, et al. Candidate gene prioritization with Endeavour. Nucleic Acids Res 2016;44(W1):W117–21. CrossrefWeb of SciencePubMedGoogle Scholar
Kwok SK, Cho ML, Her YM, Oh HJ, Park MK, Lee SY, et al. TLR2 ligation induces the production of IL-23/IL-17 via IL-6, STAT3 and NF-kB pathway in patients with primary Sjogren’s syndrome. Arthritis Res Ther 2012;14(2):R64. Web of SciencePubMedCrossrefGoogle Scholar
Umeshita-Suyama R, Sugimoto R, Akaiwa M, Arima K, Yu B, Wada M, et al. Characterization of IL-4 and IL-13 signals dependent on the human IL-13 receptor α chain 1: redundancy of requirement of tyrosine residue for STAT3 activation. Int Immunol 2000;12(11):1499–509. PubMedCrossrefGoogle Scholar
Willerslev-Olsen A, Litvinov IV, Fredholm SM, Petersen DL, Sibbesen NA, Gniadecki R, et al. IL-15 and IL-17F are differentially regulated and expressed in mycosis fungoides (MF). Cell Cycle 2014;13(8):1306–12. PubMedCrossrefWeb of ScienceGoogle Scholar
Li J, Liu M, Zong J, Tan P, Wang J, Wang X, et al. Genetic variations in IL1A and IL1RN are associated with the risk of preeclampsia in Chinese Han population. Sci Rep 2014;4:5250. PubMedWeb of ScienceGoogle Scholar
Zhou X, Ca JG, Peng H, Wang JL, Li GM. Association of IL-1alpha gene polymorphism with susceptibility to type 1 diabetes in Chinese children. Genet Mol Res. 2016;15(3). DOI: . CrossrefPubMedGoogle Scholar
Batchu SN, Hughson A, Gerloff J, Fowell DJ, Korshunov VA. Role of Axl in early kidney inflammation and progression of salt-dependent hypertension. Hypertension 2013;62(2):302–9. PubMedCrossrefWeb of ScienceGoogle Scholar
Chun HJ, Bonnet S, Chan SY. Translational advances in the field of pulmonary hypertension. Translating MicroRNA biology in pulmonary hypertension. It will take more than “miR” words. Am J Respir Crit Care Med 2017;195(2):167–78. PubMedCrossrefGoogle Scholar
Wade BE, Zhao J, Ma J, Hart CM, Sutliff RL. Hypoxia-induced alterations in the lung ubiquitin proteasome system during pulmonary hypertension pathogenesis. Pulm Circ 2018:2045894018788267. Web of SciencePubMedGoogle Scholar
Mao S, Gu T, Zhong F, Fan R, Zhu F, Ren P, et al. Hypomethylation of the Toll-like receptor-2 gene increases the risk of essential hypertension. Mol Med Rep 2017;16(1):964–70. CrossrefPubMedWeb of ScienceGoogle Scholar
Lin ZH, Fukuda N, Jin XQ, Yao EH, Ueno T, Endo M, et al. Complement 3 is involved in the synthetic phenotype and exaggerated growth of vascular smooth muscle cells from spontaneously hypertensive rats. Hypertension 2004;44(1):42–7. CrossrefPubMedGoogle Scholar
Negishi E, Fukuda N, Otsuki T, Katakawa M, Komatsu K, Chen L, et al. Involvement of complement 3 in the salt-sensitive hypertension by activation of renal renin-angiotensin system in spontaneously hypertensive rats. Am J Physiol Renal Physiol. 2018;315(6). DOI: 10.1152/ajprenal.00370.2018. Web of ScienceGoogle Scholar
Park SH, Chen WC, Esmaeil N, Lucas B, Marsh LM, Reibman J, et al. Interleukin 13–and interleukin 17A–induced pulmonary hypertension phenotype due to inhalation of antigen and fine particles from air pollution. Pulm Circ 2014;4(4):654–68. Web of SciencePubMedCrossrefGoogle Scholar
Antonelli A, Fallahi P, Ferrari SM, Ghiadoni L, Virdis A, Mancusi C, et al. High serum levels of CXC (CXCL10) and CC (CCL2) chemokines in untreated essential hypertension. Int J Immunopathol Pharmacol 2012;25(2):387–95. Web of SciencePubMedCrossrefGoogle Scholar
The online version of this article offers supplementary material (DOI: https://doi.org/10.1515/jib-2018-0054).
About the article
Published Online: 2018-12-25
Conflict of interest statement: Authors state no conflict of interest. All authors have read the journal’s Publication ethics and publication malpractice statement available at the journal’s website and hereby confirm that they comply with all its parts applicable to the present scientific work.
Citation Information: Journal of Integrative Bioinformatics, Volume 15, Issue 4, 20180054, ISSN (Online) 1613-4516, DOI: https://doi.org/10.1515/jib-2018-0054.
©2018, Olga V. Saik et al., published by Walter de Gruyter GmbH, Berlin/Boston. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0