Epidemiological study provide evidence that some diseases tend to develop together more often that it would be expected by chance, suggesting shared pathogenetic mechanisms . Various terms are used to describe this phenomenon, such as co-morbidity, syntropy, multi-morbidity. Examples of common comorbid diseases include allergic disorders (bronchial asthma, allergic rhinitis, atopic dermatitis) , autoimmune phenotypes of so-called “the kaleidoscope of autoimmunity” , cardiovascular diseases continuum .
Importantly, such combinations are common not only in an individual, but in their relatives, thus implying hereditary factors presented by pleiotropic, co-regulated syntropic genes involved in shared pathways .
Bronchial asthma, a common disease of airways, is often co-morbid with cardiovascular disorders , , , especially in women , and late-onset patients ; however, the increased rates of hypertension in asthma patients are observed in the vast majority of the studied populations .
Mechanisms underlying this common co-morbidity remain unknown. Among the discussed reasons of the development of hypertension in asthma patients are the decrease of lung function, systemic inflammation, impairment of smooth muscle tonus regulation, vascular remodelling, adverse effects of therapy, as well as genetic factors .
Systems biology provides ground for understanding molecular basis of pathogenesis of co-morbid diseases and offers novel view on diseases classification. Rzhetsky et al. built a phenotypic network of diseases including 657 using 1.5 mln clinical records based on the co-morbidity between the diseases . Both positive and negative correlations between the diseases was identified, e.g. schizophrenia and bipolar disorder (positive correlation) or aorta aneurism and schizophrenia, breast cancer and bipolar disorder (negative correlations). Another approach that used both data on disorders and genes from OMIM catalogue allowed building a network of 867 diseases with genes associated with at least two disorders . Protein interaction networks can also be used to determine molecular causes of co-morbidity .
We hypothesise that co-morbidity between asthma and hypertension is largely based on underlying genetic variation. Recently, we used a bioinformatics approach to prioritise genes that are likely important for the co-morbidity based on the analysis of associative networks . To provide further support for the hypothesis, in the current study we set out to carry out the analysis of association between asthma, hypertension and their comorbidity and functional SNPs linked to variation in gene expression (eQTL) for the genes established in the previous study.
2 Materials and Methods
2.1 Study Groups
A total of 587 individuals from West Siberia (Tomsk Region, Russia) were studied (Table 1). Participants were included in the study in case they were not related to each other and were of self-reported European ancestry. Recruitment was carried out in clinical facilities of the city of Tomsk and the city of Seversk. Diagnoses of bronchial asthma and essential hypertension were defined based on clinical and laboratory investigations of the participants in accordance with WHO accepted criteria. All patients with arterial hypertension exhibited not clinical or laboratory symptoms and signs of secondary hypertension. Individuals with the presence of systemic complex diseases such as type 1 diabetes; obesity III-IV degree and cancer were excluded from analysis. Control group comprised healthy individuals with normal arterial blood pressure and without bronchial asthma.
The study was approved by the Bioethics Committee of the Research Institute of Medical Genetics; a full verbal explanation of the investigation was given and written informed consent was obtained from all participants.
2.2 Genetic Analyses
Total genomic DNA was isolated from peripheral blood by phenol/chloroform extraction method. A total of 96 eQTL SNP were selected and genotyped by MALDI-TOF mass-spectrometry using Sequenom MassARRAY iPLEX (USA). Genotypes were established by automatic analysis using MassARRAY Typer 4 software.
2.3 Statistical Analysis
Statistical analysis was carried out in R programming environment. The analysis of associations between the SNPs and diseases was carried out using logistic regression adjusting for age and sex. Additive genetic model was tested. The Hardy-Weinberg equilibrium (HWE) was tested using Fisher’s exact test. The significance threshold was set at p < 0.05.
3 Results and Discussion
We studied association of eQTL SNPs of genes (TLR4, CAT, IL10, CST3, ICAM1, IRF6, AKT1, NFKB1, PNP, SELL, CCL5, IL2RB, IDS, FOS, NT5C2, BHLHE40 and etc.) selected by bioinformatics analysis  with bronchial asthma, hypertension and their combination. A total of 96 SNPs were studied. All of the SNPs met the expectations under HWE. For six SNPs (rs17350806, rs4987367, rs4743691, rs2284368, rs7548533, and rs4743692) genotyping rate was <80 %; therefore, they were excluded from the analysis.
3.1 Polymorphisms associated with bronchial asthma
Three polymorphisms were found to be associated with isolated bronchial asthma (Table 2). Two of them (rs1927914 and rs1928298) are located in a proximity to TLR4 and are in strong LD (r2 = 0.97) . They influence TLR4 expression in various cell types including blood cells and oesophagus smooth muscles , , . The gene encodes a pattern-recognition receptor essential in innate immunity by transducing microbial signals, such as LPS of gram-negative bacteria. The role of innate immunity in bronchial asthma pathogenesis is not fully understood; however, in patients with neutrophilic asthma, TLR4 along with related protein TLR2, are highly activated . These receptors are also thought to be linked to fatal asthma exacerbations . In the literature, there are data about the impact of TLR family genes polymorphisms on susceptibility to various phenotypes of bronchial asthma , , .
Another gene found to be associated with bronchial asthma in the current study, SERPINA1, encodes a serine protease inhibitor. The chromosomal region where the gene resides (14q32.13) is associated with the levels of IgE, bronchial hyperreactivity and other signs of allergic inflammation in different populations , , . Gene expression profiling of epithelial cells stimulated by key allergy cytokines IL-4 and IL-13 showed most remarkable increase of expression of serpin family genes SERPINB3 и SERPINB4 . Their role in the pathogenesis of bronchial asthma is not known, but they can exhibit protective effects via inhibiting endogenous proteases associated with inflammatory response.
3.2 Polymorphisms associated with hypertension
A polymorphism rs11065987 associated with hypertension in the current study have been found to be associated in genome-wide association studies (GWAS) with a number of phenotypes related to the risks of cardiovascular diseases, including lipid levels ,  and BMI . This variant is considered as the marker of the locus of susceptibility to ischemic stroke and ischemic heart disease .
Another variant associated with hypertension in the current study resides in IL2-receptor gene IL2RB, which takes part in T-cell mediated immune response. Previously, IL2RB gene expression was found to be correlated with diastolic blood pressure levels in a large sample of 7017 individuals .
The rs11191582 SNP associated with hypertension in the current study is located in the gene-reach region near CYP17A1-CNNM2-NT5C2. The region contains a remarkable number of regulatory SNPs associated with the development of cardiovascular disorders in GWAS studies . The SNP rs11191582 is a cis-eQTL for NT5C2, USMG5, and hsa-mir-1307 genes that are abundantly expressed in the whole blood .
Another polymorphism associated with hypertension in the current study, rs11669386, is located in an intron of the caspase recruitment domain family member 8 gene (CARD8). The gene is involved in apoptosis, inflammation and NF-kappaB signaling. CARD8 along with NLRP3 controls activity of caspase-1 and is involved in caspase-1-mediated activation of IL-1β . In turn, both caspase-1 and IL-1β are crucial players in the development of coronary pathology in mouse model of Kawasaki disease .
3.3 Polymorphisms associated with the comorbidity between bronchial asthma and hypertension
The polymorphism rs1010461 associated with the comorbidity, is located in intron 1 of angiogenin gene (ANG) and is an eQTL for ANG and RNASE4 genes with shared promoter. The proteins encoded by ANG and RNASE4 genes have 38.7 % similarity and belong to ribonuclease family (RNase). Their major function is induction of neovascularization . Angiogenic activity of these genes may be important for interpretation of their associations with the comorbidity. Angiogenin is induced in inflammation and exhibit wound-healing effects and microbicide activity . In patients with arterial hypertension, the decrease of angiogenin and fibroblasts growth factor levels is observed along with the increase of the levels of endostatin, vascular endothelial grown factor, and IL-8, suggesting a disbalance of pro-/anti-angiogenic factors in this condition . Increased vascularisation promotes bronchial asthma pathogenesis and correlates with several signs of clinical course of the disease including bronchial hyperreactivity and airflow limitation . The levels of angiogenin and vascular endothelial grown factor are increased during asthma exacerbation and decreased followed by corticosteroid treatment . This raises a possibility that this is corticosteroid treatment that promotes the development of hypertension in genetically predisposed asthma patients.
The comorbidity between bronchial asthma and hypertension was associated with three SNPs known to be eQTLs for TLR4 gene according to GTEx database: two LD-bound polymorphisms rs7038716 and rs7026297 as well as intergenic rs7025144 SNP.
Finally, rs2022318 SNP, an eQTL for CAT gene was found to be associated with the comorbidity between asthma and hypertension. Catalase encoded by the gene is an important antioxidant enzyme preventing harmful effects of active oxygen species. As a result, polymorphisms of this gene are widely studied in diseases though to be linked with oxidative stress. Promoter SNPs for this gene were found to be associated with cardiometabolic diseases , , bronchial asthma , type 1 diabetes and other diseases .
The important observation from the current study is the difference in SNPs associated with isolated bronchial asthma and hypertension or their comorbidity, thus suggesting different genetic background of the comorbidity. On the other hand, considering associations of TLR4 SNPs, this gene variation may be important both for the development of isolated bronchial asthma and comorbidity of bronchial asthma and hypertension. The most pronounced impact on the development of comorbidity between bronchial asthma and hypertension can be attributed to the TLR4, CAT and ANG/RNASE4 genes, suggesting the importance of inflammation, neovascularisation, and oxidative stress in co-pathogenesis of both diseases. Importantly, according to our bioinformatics analysis, TLR4 and CAT genes have the highest priority in the development of bronchial asthma and hypertension comorbidity . Consistently, we observed associations of these genes with the comorbidity, thus providing evidence for productivity in the use of bioinformatics approaches for the planning of human association studies.
Currently not enough of GWAS are published aimed at the study of comorbidity. These works aimed to study psychiatric diseases , . Our study, unlike GWAS, focuses on specific SNPs of genes found using bioinformatics approaches that can cause comorbidity of asthma and hypertension.
This study has limitations, including no correction for multiple testing, small sample sizes, wide age range of the patients, and no correction for the use of therapy. Thus, the results of the current study must be treated with caution and considered as a starting point for further development of studies in this direction.
Future large-scale studies, such as GWAS, are needed to elucidate the molecular mechanisms of asthma and hypertension comorbidity; however, this will require large samples to evaluate the effects of genetic variants with moderate and small effects.
This study was supported by 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” (# 90335). Project of the Integrated Basic Research Program of the SB RAS II.1 “Reconstruction, computer analysis and modelling of the structural and functional organization of biomedically significant gene networks” (# 0324-2018-0021; # 0550-2018-0003).
Anaya JM, Corena R, Castiblanco J, Rojas-Villarraga A, Shoenfeld Y. The kaleidoscope of autoimmunity: multiple autoimmune syndromes and familial autoimmunity. Expert Rev Clin Immunol 2007;3:623–35. PubMedCrossrefGoogle Scholar
Dzau V, Braunwald E. Resolved and unresolved issues in the prevention and treatment of coronary artery disease: a workshop consensus statement. Am Heart J 1991;121:1244–63. PubMedCrossrefGoogle Scholar
Johnson M, Nriagu J, Hammad A, Savoie K, Jamil H. Asthma, environmental risk factors, and hypertension among Arab Americans in metro Detroit. J Immigr Minor Health 2010;12:640–51. PubMedCrossrefWeb of ScienceGoogle Scholar
Lee HS, Park YM, Han K, Pekler G, Lee SS, Yoo S, et al. Sex-specific association between asthma and hypertension in nationally representative young Korean adults. Sci Rep 2017;7:15667. Web of SciencePubMedCrossrefGoogle Scholar
Onufrak S, Abramson J, Vaccarino V. Adult-onset asthma is associated with increased carotid atherosclerosis among women in the Atherosclerosis Risk in Communities (ARIC) study. Atherosclerosis 2007;195:129–37. CrossrefPubMedWeb of ScienceGoogle Scholar
Ferguson S, Teodorescu MC, Gangnon RE, Peterson AG, Consens FB, Chervin RD, et al. Factors associated with systemic hypertension in asthma. Lung 2014;192:675–83. Web of SciencePubMedCrossrefGoogle 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(Suppl 1):15. Web of SciencePubMedCrossrefGoogle Scholar
Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res 2012;40(Database issue):D930–4. Web of SciencePubMedCrossrefGoogle Scholar
Westra HJ, Peters MJ, Esko T, Yaghootkar H, Schurmann C, Kettunen J, et al. Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat Genet 2013;45:1238–43. CrossrefPubMedWeb of ScienceGoogle Scholar
Lappalainen T, Sammeth M, Friedländer MR, ’t Hoen PA, Monlong J, Rivas MA, et al. Transcriptome and genome sequencing uncovers functional variation in humans. Nature 2013;501:506–11. PubMedWeb of ScienceCrossrefGoogle Scholar
Ferreira DS, Annoni R, Silva LF, Buttignol M, Santos AB, Medeiros MC, et al. Toll-like receptors 2, 3 and 4 and thymic stromal lymphopoietin expression in fatal asthma. Clin Exp Allergy 2012;42:1459–71. PubMedWeb of ScienceCrossrefGoogle Scholar
Zhang Q, Qian FH, Zhou LF, Wei GZ, Jin GF, Bai JL, et al. Polymorphisms in toll-like receptor 4 gene are associated with asthma severity but not susceptibility in a Chinese Han population. J Investig Allergol Clin Immunol 2011;21:370–7. Google Scholar
Tizaoui K, Kaabachi W, Hamzaoui K, Hamzaoui A. Association of single nucleotide polymorphisms in toll-like receptor genes with asthma risk: a systematic review and meta-analysis. Allergy Asthma Immunol Res 2015;7:130–40. Web of ScienceCrossrefGoogle Scholar
Lau MY, Dharmage SC, Burgess JA, Win AK, Lowe AJ, Lodge C, et al. The interaction between farming/rural environment and TLR2, TLR4, TLR6 and CD14 genetic polymorphisms in relation to early- and late-onset asthma. Sci Rep 2017;7:43681. CrossrefPubMedWeb of ScienceGoogle Scholar
Mansur AH, Bishop DT, Markham AF, Morton NE, Holgate ST, Morrison JF. Suggestive evidence for genetic linkage between IgE phenotypes and chromosome 14q markers. Am J Respir Crit Care Med 1999;159:1796–802. PubMedCrossrefGoogle Scholar
Ober C, Tsalenko A, Parry R, Cox NJ. A second-generation genomewide screen for asthma-susceptibility alleles in a founder population. Am J Hum Genet 2000;67:1154–62. Google Scholar
Malerba G, Patuzzo C, Trabetti E, Lauciello MC, Galavotti R, Pescollderungg L, et al. Chromosome 14 linkage analysis and mutation study of 2 serpin genes in allergic asthmatic families. J Allergy Clin Immunol 2001;107:654–8. PubMedCrossrefGoogle Scholar
Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet 2013;45:1274–83. Web of ScienceCrossrefPubMedGoogle Scholar
Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature 2010;466:707–13. Web of SciencePubMedCrossrefGoogle Scholar
Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature 2015;518:197–206. PubMedCrossrefWeb of ScienceGoogle Scholar
Dichgans M, Malik R, König IR, Rosand J, Clarke R, Gretarsdottir S, et al. Shared genetic susceptibility to ischemic stroke and coronary artery disease: a genome-wide analysis of common variants. Stroke 2014;45:24–36. PubMedCrossrefWeb of ScienceGoogle Scholar
Huan T, Esko T, Peters MJ, Pilling LC, Schramm K, Schurmann C, et al. A meta-analysis of gene expression signatures of blood pressure and hypertension. PLoS Genet 2015;11:e1005035. CrossrefWeb of SciencePubMedGoogle Scholar
Cheema AN, Rosenthal SL, Ilyas Kamboh M. Proficiency of data interpretation: identification of signaling SNPs/specific loci for coronary artery disease. Database (Oxford) 2017;2017. doi: 10.1093/database/bax078. PubMedWeb of ScienceGoogle Scholar
Lee Y, Schulte DJ, Shimada K, Chen S, Crother TR, Chiba N, et al. Interleukin-1β is crucial for the induction of coronary artery inflammation in a mouse model of Kawasaki disease. Circulation 2012;125:1542–50. Web of ScienceCrossrefGoogle Scholar
Li S, Sheng J, Hu JK, Yu W, Kishikawa H, Hu MG, et al. Ribonuclease 4 protects neuron degeneration by promoting angiogenesis, neurogenesis, and neuronal survival under stress. Angiogenesis 2013;16:387–404. Web of ScienceCrossrefPubMedGoogle Scholar
Marek-Trzonkowska N, Kwieczyńska A, Reiwer-Gostomska M, Koliński T, Molisz A, Siebert J. Arterial hypertension is characterized by imbalance of pro-angiogenic versus anti-angiogenic factors. PLoS One 2015;10:e0126190. CrossrefPubMedWeb of ScienceGoogle Scholar
Chetta A, Zanini A, Torre O, Olivieri D. Vascular remodelling and angiogenesis in asthma: morphological aspects and pharmacological modulation. Inflamm Allergy Drug Targets 2007;6:41–5. PubMedCrossrefGoogle Scholar
Abdel-Rahman AM, el-Sahrigy SA, Bakr SI. A comparative study of two angiogenic factors: vascular endothelial growth factor and angiogenin in induced sputum from asthmatic children in acute attack. Chest 2006;129:266–71. PubMedCrossrefGoogle Scholar
Zhou XF, Cui J, DeStefano AL, Chazaro I, Farrer LA, Manolis AJ, et al. Polymorphisms in the promoter region of catalase gene and essential hypertension. Dis Markers 2005;21:3–7. PubMedCrossrefGoogle Scholar
Polonikov AV, Ivanov VP, Solodilova MA, Kozhuhov MA, Panfilov VI. Tobacco smoking, fruit and vegetable intake modify association between -21A>T polymorphism of catalase gene and risk of bronchial asthma. J Asthma 2009;46:217–24. CrossrefPubMedWeb of ScienceGoogle Scholar
Zhou H, Polimanti R, Yang BZ, Wang Q, Han S, Sherva R, et al. Genetic Risk Variants Associated With Comorbid Alcohol Dependence and Major Depression. JAMA Psychiatry 2017;74:1234–41. Web of ScienceCrossrefPubMedGoogle Scholar
Edwards AC, Aliev F, Bierut LJ, Bucholz KK, Edenberg H, Hesselbrock V, et al. Genome-wide association study of comorbid depressive syndrome and alcohol dependence. Psychiatr Genet 2012;22:31–41. CrossrefPubMedWeb of ScienceGoogle Scholar
About the article
Published Online: 2018-12-10
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, 20180052, ISSN (Online) 1613-4516, DOI: https://doi.org/10.1515/jib-2018-0052.
©2018, Elena Yu. Bragina 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