- Open Access
Shared genetic variants between serum levels of high-density lipoprotein cholesterol and wheezing in a cohort of children from Cyprus
Italian Journal of Pediatrics volume 42, Article number: 67 (2016)
In a cohort of children in Cyprus, we recently reported low levels of high density lipoprotein cholesterol (HDL-C) to be associated with asthma. We examined whether genetic polymorphisms that were previously linked individually to asthma, obesity, or HDL-C are associated with both asthma and HDL-C levels in the Cyprus cohort.
We assessed genotypes frequencies in current-wheezers (n = 190) and non-asthmatic controls (n = 671) and HDL-C levels across several genotypes. Binary logistic regression models were used to assess the effect of genotypes on wheezing risk and examined whether this effect is carried out through changes of HDL–C.
Of the 16 polymorphisms tested, two polymorphisms TNFa rs3093664 and PRKCA rs9892651 presented significant differences in genotype distribution among current-wheezers and controls. Higher HDL-C levels were noted in carriers of genotype GG of polymorphism TNFa rs3093664 that was protective for wheezing Vs AG and AA genotypes (65.3 Vs 51.8 and 53.3 mg/dl, p-value < 0.001 and p-value for trend = 0.028). In polymorphism PRKCA rs9892651, HDL-C levels were lower in carriers of CC and TC genotypes that were more frequent in current-wheezers Vs TT genotype (52.2 and 52.7 Vs 55.2 mg/dl, p-value = 0.042 and p-value for trend = 0.02). The association of TNFa rs3093664 with wheezing is partly mediated by its effect on HDL-C whereas association of PRKCA rs9892651 with wheezing appeared to be independent of HDL-C.
We found evidence that two SNPs located in different genetic loci, are associated with both wheezing and HDL-C levels, although more studies in other populations are needed to confirm our results.
High density lipoprotein cholesterol (HDL-C) is widely known to have systemic and vascular anti-inflammatory properties  and, as shown in vitro, also anti-inflammatory effects consistent with a protective role in asthma . In a cohort of 3982 children in Cyprus, we recently reported low levels of HDL-C at age 11–12 years to be associated with active asthma five years later in adolescence, suggesting a potential role of this lipoprotein in the pathogenesis of pediatric asthma . Obesity is also widely accepted to be related to asthma  as well as with low serum HDL-C levels  and consequently it may play a confounding or mediating role in the asthma – HDL-C association.
Numerous individual studies confirmed the implication of genetic factors in the development of asthma , obesity  or low HDL-C . In addition, epidemiological studies in twins have suggested that the association of asthma and allergy with obesity or low HDL-C levels could, at least in part, be explained by the presence of shared genetic factors . In the case of asthma and obesity, overlapping genomic regions for common susceptibility genes on chromosomes 5, 6, 11 and 12 have been identified by linkage analysis studies . Specifically, some of the most known candidate genes and genetic loci for asthma and allergy [6, 11] have already been linked  to the pathogenesis of obesity such as the tumor necrosis factor alpha (TNF-a) gene , the protein kinase C alpha (PRKCA) gene  and β2‐adrenergic receptor (ADRB2) . To date however, no common genetic loci were reported to be associated with asthma and HDL-C levels.
In this study, we selected 16 single nucleotide polymorphisms (SNPs) located in genes that were linked to asthma and/or obesity or HDL-C and examined in the Cyprus cohort of children a) whether they were associated with both asthma and HDL-C levels and b) explored, in the SNPs that were associated with both conditions, whether the association with asthma is mediated through changes in the levels of HDL-C.
The participants of this study were selected from a cohort of 3982 children who participated in two large school-based health surveys in Cyprus (Fig. 1). The first survey spanned years 2001–2003 and involved all children (n = 19,849) attending the 6th grade from all primary schools across Cyprus and focused on nutrition and physical fitness. The second survey was conducted in year 2007, followed-up 3982 of the subjects who took part in the first survey (20.1 %) and focused on respiratory health (International Study of Asthma and Allergies in Childhood – ISAAC) and risk factors for asthma. In a case-control design, we invited all those from the 3982 children who on the second survey (ISAAC questionnaire) reported wheezing in the past 12 months (n = 297) to participate in this study along with a triplicate number of healthy controls (n = 932).
Current wheezers (CUWH) were named those who reported wheezing in the past 12 months. For the purpose of performing a sensitivity analysis, the case definition was further refined to Active Asthma (ACAS), if also there was report of diagnosis of asthma ever. Controls were selected amongst the 3982 adolescents who did not report any wheezing or asthma ever (Never Wheezers Never Asthmatics - NWNA). NWNA were selected using a stratified random sampling approach in order to increase the probability of selection of children at the extremes of Body Mass Index (BMI) change between childhood and adolescence, in line with the scope of another study on the relation of adiposity and serum lipids with asthma . The full scale of BMI change (and BMI) was represented in the final sample, just in slightly different proportions compared to the general population i.e. there was slight over-representation of participants whose BMI increased or decreased between childhood and adolescence and, equivalently slight under-representation of children whose BMI remained stable over this period. In all analyses normalized sampling weights based on inverse probability of selection from the original population were used.
In the asthmatic groups, information was also obtained regarding the presence of exercise induced wheezing and admissions ever in hospital for an asthma attack. IgE sensitization was assessed with skin prick tests (SPT) to 8 common aeroallergens (Greer, USA), performed and interpreted according to the GA2LEN recommendations . Children with a positive SPT to any of the 8 allergens tested were considered sensitized to IgE. Asthmatic children had measurements of forced expiratory volume in one second (FEV1) (Vitallograph, UK) and values were expressed as percent of predicted .
The adolescents had anthropometric measurements and assessments of serum total cholesterol (TC), low density lipoprotein cholesterol (LDL-C), HDL-C and triglycerides (TG) after a 12-h fast. We calculated BMI (kg/m2) and expressed it as age- and gender-specific z-scores based on the United States’ Centers for Disease Control 2000 growth charts . Lipids were determined with the ILab 600 analyzer using ILab Chemistry systems (Instrumentation Laboratory SpA, Monza, Milano, Italy). The assays’ coefficients of variation were 2.08 % for TC, 1.53 % for HDL-C, 1.15 % for LDL-C and 2.45 % for TG.
SNPs selection, DNA extraction and genotyping
Using the NCBI PubMed Entrez database, a total of 16 SNPs in ten genes that were previously reported [12–14, 19–30] to be associated individually with asthma, obesity or HDL-C were selected for genotyping (Table 1).
Venous blood specimens were collected from participants in EDTA treated tubes that were transported on ice to the laboratory. DNA was extracted from peripheral blood lymphocytes using the standard phenol extraction method, and stored at −70 °C until genotyping. Genotyping was carried out using TaqMan Real-time PCR assay. Thirty ng of genomic DNA was used for each reaction and the end-point PCR was carried out on the ABI primer 7900 HT (Applied Biosystems, Foster City, CA).
Differences in participant characteristics between study-groups were investigated in a pairwise manner (NWNA versus CUWH and NWNA versus ACAS) using chi square test in the case of categorical variables and t-test in the case of continuous variables. Furthermore, chi-squared testing was used for Hardy-Weinberg (HW) equilibrium determination and for the investigation of genotypes and allele frequencies associations with study groups. Mean HDL-C levels were compared across genotypes using the one-way ANOVA test and p values and p values for trend were reported. For SNPs that presented significant effects, two binary logistic regression models were used a) to assess the effect of genotypes on asthma risk and b) examine whether this effect is carried out through the change of HDL–C, after adjusting for age and sex (Model 1), age, sex and BMI (Model 2) and age, sex, BMI and HDL-C levels (Model 3). The genotype with the most protective effect for asthma was chosen as the reference in all cases. Binary logistic regression analysis was performed only among NWNA and CUWH as the ACAS group was very small and for most SNPs, at least one cell had no cases. Analysis was performed with SPSS statistical package, version 20.0 (IBM, SPSS Inc., Chicago, IL).
A total of 190 subjects with CUWH (of whom 69 ACAS) and 671 NWNA participated in the study. This corresponds to 64 and 72 % response rates respectively of the targeted population for recruitment (Fig. 1). Participants were aged 17 years (SD 0.6) and the three groups did not differ significantly in terms of age, gender and BMI (Table 2). Mean HDL-C levels in the ACAS and CUWH were significantly lower than in the controls (47.9 and 49.1 Vs 53.4 mg/dl; p < 0.001 for both) (Table 2). Mean serum levels of TC, LDL-C and TG were similar in the three groups (results not shown in detail). The frequency of atopic sensitization was 52.0 % in ACAS and 41.4 % in CUWH (p-value = 0.073), whereas FEV1 was normal and not different between them (97.6 % Vs 99.1 % of predicted; p-value = 0.332) indicating that they were clinically stable. ACAS were experiencing more frequently exercise induced wheezing (83.6 % vs 67.6 %, p-value < 0.001) and lifetime hospitalizations for asthma attack (18.8 % vs 9.0 %, p-value = 0.003) than the CUWH (Table 2). Results of minor allele frequencies and HW equilibrium for genotype distribution are presented in Additional file 1: Table S1.
SNPs genotype distribution in CUWH Vs NWNA and HDL-C levels
Distribution of the genotypes of the 16 polymorphisms in CUWH Vs NWNA and the corresponding HDL-C levels can be found in Additional file 2: Table S2. Two polymorphisms, TNFa rs3093664 and PRKCA rs9892651 presented statistically significant and ADRB2 rs1800888 nearly significant differences in genotype distribution among NWNA and CUWH and are presented separately at Table 3.
The larger differences were exhibited in the genotype distribution of TNFa rs3093664 polymorphism, where the rare GG genotype appeared to offer a protective effect with a frequency of 4.1 % Vs 1.1 % in NWNA and CUWH respectively (p-value = 0.03, p-value for trend = 0.912). Significantly higher HDL-C levels were noted in subjects with GG genotype as opposed to those in AG and AA genotype carriers (65.3 mg/dl Vs 51.8 mg/dl and 53.3 mg/dl, p-value < 0.001 and p-value for trend = 0.028) (Table 3). The genotypes of polymorphism PRKCA rs9892651 were also found to be differentially distributed with the frequency of genotype CC to be higher in CUWH (19.2 %) compared to NWNA (14.3 %) (p-value = 0.141, p-value for trend = 0.048). HDL-C levels were also significantly lower in subjects with the CC and TC genotype compared to TT genotype (52.2 mg/dl and 52.7 mg/dl Vs 55.2 mg/dl, p-value = 0.042 and p-value for trend = 0.02).
The only other SNP that displayed nearly significant differences in genotype frequencies was ADRB2 rs1800888, where the heterozygous genotype CT had a higher frequency in NWNA compared to CUWH (5 % Vs 2.2 %, p-value = 0.076) thus exhibiting a protective effect compared with the homozygote CC genotype. HDL-C levels in the protective CT genotype carriers were not significantly higher at the 5 % level to those in CC, probably due to the small size of the CT genotype group (57.24 mg/dl Vs 53.28 mg/dl, p-value = 0.092).
SNPs genotype distribution in ACAS Vs NWNA and HDL-C levels
As a sensitivity analysis, the same investigations were performed between the more strictly defined asthmatic phenotype ACAS and NWNA. In general, the magnitude and direction of the relationships of the three SNPs genotypes in ACAS and NWNA remained consistent with the findings of the NWNA Vs CUWH analysis although these were not always significant probably due to the smaller size of the ACAS study-group. A detailed description of these findings is presented with Additional file 3: Table S3. Allelic associations of the 16 tested SNPs with study groups were not statistically different and are presented in Additional file 4: Table S4.
Associations of SNPs genotypes with CUWH and the mediating role of HDL-C levels
Associations of rs3093664, rs9892651 and rs1800888 genotypes with CUWH were explored in binary logistic regression models adjusted for age and sex (Model 1) age, sex and BMI (Model 2) and with the inclusion of HDL-C levels (Model 3) (Table 4). The estimates in all comparisons between Models 1 and 2 were very similar, thus excluding any confounding or mediating effect of BMI. In particular, for the TNF-a rs3093664 polymorphism, the Odds Ratios (ORs) of having CUWH for genotypes AA and AG against GG were 3.35 (95 % CI 0.78, 14.35) and 4.92 (95 % CI 1.11, 21.8) respectively (Model 2). After further adjusting for HDL-C (Model 3), the calculated ORs of CUWH were attenuated by approximately 40 % (OR = 2.44, 95 % CI 0.56, 10.58 for genotype AA and OR = 3.43 95 % CI 0.76, 15.53 for genotype AG), thus indicating that the protective effect of GG for CUWH is mediated to an important extent through HDL-C levels.
In the case of PRKCA rs9892651, the ORs for genotype CC carriers to have CUWH compared to TT (OR = 1.63, 95 % CI 0.99, 2.69) and CT (OR = 1.26, 95 % CI 0.85, 1.85) were slightly short of significance in Model 2. In Model 3, these estimates were only moderately attenuated for TT (OR = 1.52, 95 % CI 0.91, 2.53) and CT (OR = 1.17, 95 % CI 0.79, 1.73). Thus PRKCA rs9892651 association with wheezing appeared to be independent of HDL-C levels.
In this case-control study, we examined the effect of 16 genetic variants on asthma phenotypes and identified two SNPs that were simultaneously associated with the wheezing and HDL-C levels. Although many studies have examined separately the relationship of genetic loci with respiratory outcomes, only few have searched for common genetic variants between asthma and metabolic markers [12, 13, 21, 31]. Skaaby et al., studied loss-of-function mutations in the filaggrin gene, which are known to be associated with asthma and allergies and higher vitamin D, and found higher vitamin D status to be associated with higher HDL-C levels , suggesting an indirect potential genetic link between asthma and allergies with high HDL-C levels. In our cohort, we have not examined polymorphisms in the filaggrin gene, but to the best of our knowledge, this is the only study so far which, in the same population, assessed simultaneously the direct association of selected SNPs with wheezing and HDL-C levels. Interestingly, in our study the direction of the association of wheezing and HDL-C levels is in the opposite direction from the one implied in the Skaaby et al. study. More importantly, our results provide preliminary evidence that the association of TNFa rs3093664 with wheezing is mediated in part by its effect on HDL-C levels.
Evidence for underlying mechanisms of action
TNF-a gene codes for a highly pro-inflammatory cytokine  involved in many pathophysiological mechanisms of asthma [34, 35]. TNF-a is also known to affect lipid metabolism and lipid serum levels [36, 37]. Specific polymorphisms of the TNFa gene have been reported to affect HDL-C levels  and have been proposed as risk factors towards obesity, metabolic syndrome and type 2 diabetes mellitus . Our findings point towards a novel protective effect of genotype GG of TNFa rs3093664 against wheezing but also at the same time suggest a strong association of this genotype with higher levels of HDL-C. To date, only few studies have examined the association of this particular polymorphism with asthma  or any metabolic condition . Although TNF-a serum levels were not assessed in this study, previous reports have documented that TNF-a serum levels are inversely associated with serum HDL-C in the general population , whereas administration of anti TNF-a agents in patients with rheumatic diseases has been shown to result in elevated HDL-C levels [42, 43]. Coupled with these observations, our results suggest that in this small group of carriers of the GG genotype we may have a compromised TNF-a pro-inflammatory activity, manifesting as significantly lower risk for wheezing. At the same time, the direct protective effect of genotype GG for wheezing is further enhanced by the higher HDL-C levels that result from the reduced activity of TNF-a and its inflammatory properties .
Enzyme PKC-a has been proposed to be involved in asthma pathogenesis , in adipogenesis  and lipid mediated insulin resistance . The potential pleiotropic effect of PRKCA gene on asthma and BMI was demonstrated by different SNPs influence on asthma or obesity [13, 21], including rs9892651. Our findings also support a pleiotropic effect of PRKCA rs9892651 as the minor C allele appeared to be more prevalent among current wheezers and also associated with lower HDL-C levels, although the effect on the risk of wheezing appeared to be very modest.
Polymorphisms of ADRB2 gene have been studied extensively in asthma but also in lipolysis of the adipose tissue . We examined three of the most well studied polymorphisms of ADRB2 gene, the common rs1042714 and rs1042713 and the rarer rs1800888 [48, 49]. Our study provides evidence in support of a modest protective role for wheezing of T allele of the rare polymorphism rs1800888 but also for a modest association with higher HDL-C levels. This is the first study which evaluated the effects of rs1800888 on HDL-C levels since previous studies mainly examined cardiovascular outcomes with conflicting results [50–52].
The use of self-reported questionnaire data obtained from a young population to define wheezing and select cases and controls is not as accurate as the use of clinical parameters. We do have clinical phenotypic features of asthma severity in the selected cases populations, as depicted in previous studies , but restriction of the analyses to those with features of asthma severity would have reduced the numbers substantially and weaken any possible statistical associations. We think however that any misclassification or recall bias in defining the disease status is unlikely to have been influenced by the genotypes of the examined polymorphisms and thus confound the recovered associations. Furthermore, as five SNPs deviated from HW equilibrium, a random subset (10 %) of samples was genotyped in duplicate for quality control and reproducibility of the genotype calls was confirmed. Populations may deviate from Hardy-Weinberg equilibrium due to selection pressure, structured populations, in-migration, fluctuations in small populations, and if mutations are recent. It might be expected that selection pressure might account for at least some disequilibrium in the disease associated SNPs identified in this study, but other SNPs that did not associate with this particular disease also showed disequilibrium. Even if other disease processes might account for those, this is also a small island population that has experienced multiple historic and recent in-migrations, and which is also a structured population. The number of SNPs studied was limited and although we cannot exclude that other genetic loci that were not tested might also be implicated in the association between wheezing and HDL, this does not affect the validity of the recovered associations. On the other hand, we did investigate several genes/genotypes and although we had some a priori evidence to support our hypotheses, the associations we found do not hold after adjusting for multiple testing by the Bonferroni method, which underlines the need of confirmation of our results in other populations. Finally, this study remains a genetically-focused report without functional phenotypic data on the analyzed SNP variants and HDL-related metabolites. Nevertheless, our novel findings call for further studies to explore the physiological mechanisms behind these associations.
In conclusion, this study provides a first indication that two SNPs located in different genetic loci, TNFa rs3093664 and PRKCA rs9892651, are associated with both the wheezing and HDL-C levels. The association of TNFa rs3093664 with wheezing was mediated in part by its effect on HDL-C levels whereas the associations of PRKCA rs9892651 with wheezing appeared to be independent of HDL-C levels. Further studies are needed to confirm these findings and reveal the functional mechanisms implicated in the interactions of HDL with wheezing and genetic loci expression.
ACAS, active asthmatics; ADRB2, β2‐adrenergic receptor; CUWH, current wheezers; HDL-C, high density lipoprotein cholesterol; NWNA, never wheezers never asthmatics; OR, odds ratio; PRKCA, protein kinase C alpha; SNPs, single nucleotide polymorphisms; TNF-a, tumor necrosis factor alpha
Barter PJ, Nicholls S, Rye KA, Anantharamaiah GM, Navab M, Fogelman AM. Antiinflammatory properties of HDL. Circ Res. 2004;95(8):764–72.
Otera H, Ishida T, Nishiuma T, Kobayashi K, Kotani Y, Yasuda T, Kundu RK, Quertermous T, Hirata K, Nishimura Y. Targeted inactivation of endothelial lipase attenuates lung allergic inflammation through raising plasma HDL level and inhibiting eosinophil infiltration. Am J Physiol Lung Cell Mol Physiol. 2009;296(4):L594–602.
Yiallouros P, Savva S, Kolokotroni O, Behbod B, Zeniou M, Economou M, Chadjigeorgiou C, Kourides Y, Tornaritis M, Lamnisos D. Low serum high‐density lipoprotein cholesterol in childhood is associated with adolescent asthma. Clin Exp Allergy. 2012;42(3):423–32.
Salome C, Marks G. Sex, asthma and obesity: an intimate relationship? Clin Exp Allergy. 2011;41(1):6–8.
Eisenmann JC, Welk GJ, Ihmels M, Dollman J. Fatness, fitness, and cardiovascular disease risk factors in children and adolescents. Med Sci Sports Exerc. 2007;39(8):1251–6.
Vercelli D. Discovering susceptibility genes for asthma and allergy. Nat Rev Immunol. 2008;8(3):169–82.
Day FR, Loos RJ. Developments in obesity genetics in the era of genome-wide association studies. J Nutrigenet Nutrigenomics. 2011;4(4):222–38.
Weissglas-Volkov D, Pajukanta P. Genetic causes of high and low serum HDL-cholesterol. J Lipid Res. 2010;51(8):2032–57.
Becker KG. The common variants/multiple disease hypothesis of common complex genetic disorders. Med Hypotheses. 2004;62(2):309–17.
Litonjua AA, Gold DR. Asthma and obesity: common early-life influences in the inception of disease. J Allergy Clin Immunol. 2008;121(5):1075–84.
Ober C, Hoffjan S. Asthma genetics 2006: the long and winding road to gene discovery. Genes Immun. 2006;7(2):95–100.
Castro-Giner F, Kogevinas M, Imboden M, de Cid R, Jarvis D, Machler M, Berger W, Burney P, Franklin KA, Gonzalez JR, Heinrich J, Janson C, Omenaas E, Pin I, Rochat T, Sunyer J, Wjst M, Anto JM, Estivill X, Probst-Hensch NM. Joint effect of obesity and TNFA variability on asthma: two international cohort studies. Eur Respir J. 2009;33(5):1003–9.
Murphy A, Tantisira KG, Soto-Quirós ME, Avila L, Klanderman BJ, Lake S, Weiss ST, Celedón JC. < i > PRKCA</i>: A Positional Candidate Gene for Body Mass Index and Asthma. Am J Hum Genet. 2009;85(1):87–96.
Rankinen T, Zuberi A, Chagnon YC, Weisnagel SJ, Argyropoulos G, Walts B, Pérusse L, Bouchard C. The human obesity gene map: the 2005 update. Obesity. 2006;14(4):529–644.
Yiallouros PK, Savva SC, Kolokotroni O, Dima K, Zerva A, Kouis P, Bousquet J, Middleton N. Asthma: The Role of Low High-Density-Lipoprotein Cholesterol in Childhood and Adolescence. Int Arch Allergy Immunol. 2014;165(2):91–9.
Bousquet J, Heinzerling L, Bachert C, Papadopoulos NG, Bousquet PJ, Burney PG, Canonica GW, Carlsen KH, Cox L, Haahtela T, Lodrup Carlsen KC, Price D, Samolinski B, Simons FER, Wickman M, Annesi-Maesano I, Baena-Cagnani CE, Bergmann KC, Bindslev-Jensen C, Casale TB, Chiriac A, Cruz AA, Dubakiene R, Durham SR, Fokkens WJ, Gerth-van-Wijk R, Kalayci O, Kowalski ML, Mari A, Mullol J, Nazamova-Baranova L, O’Hehir RE, Ohta K, Panzner P, Passalacqua G, Ring J, Rogala B, Romano A, Ryan D, Schmid-Grendelmeier P, Todo-Bom A, Valenta R, Woehrl S, Yusuf OM, Zuberbier T, Demoly P. Practical guide to skin prick tests in allergy to aeroallergens. Allergy. 2012;67(1):18–24.
Knudson RJ, Lebowitz MD, Holberg CJ, Burrows B. Changes inthe normal maximal expiratory flow-volume curve with growth and aging. Am Rev Respir Dis. 1983;127:725–34.
Kuczmarski RJ, Ogden CL, Guo SS, Grummer-Strawn LM, Flegal KM, Mei Z, Wei R, Curtin LR, Roche AF, Johnson CL: 2000 CDC Growth Charts for the United States: methods and development. Vital Health Stat 11 2002, (246):1–190.
Gudbjartsson DF, Bjornsdottir US, Halapi E, Helgadottir A, Sulem P, Jonsdottir GM, Thorleifsson G, Helgadottir H, Steinthorsdottir V, Stefansson H. Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction. Nat Genet. 2009;41(3):342–7.
Bottini E, Bergamaschi A, Magrini A, et al. Allergy and ACP1 genetic polymorphism. Allergy Asthma Proc. 2007;28:87–92. OceanSide Publications, Inc.
Melén E, Himes BE, Brehm JM, Boutaoui N, Klanderman BJ, Sylvia JS, Lasky-Su J. Analyses of shared genetic factors between asthma and obesity in children. J Allergy Clin Immunol. 2010;126(3):631–7. e8.
Cui L, Jia J, Ma C, Li S, Wang Y, Guo X, Li Q, Yu H, Liu W, Gao L. IL-13 polymorphisms contribute to the risk of asthma: A meta-analysis. Clin Biochem. 2012;45(4–5):285–8.
Hall IP, Blakey JD, Al Balushi KA, Wheatley A, Sayers I, Pembrey ME, Ring SM, McArdle WL, Strachan DP. β2-adrenoceptor polymorphisms and asthma from childhood to middle age in the British 1958 birth cohort: a genetic association study. Lancet. 2006;368(9537):771–9.
Park HS, Shin ES, Lee JE. Genotypes and haplotypes of < i > β</i > < sub > 2</sub > −adrenergic receptor and parameters of the metabolic syndrome in Korean adolescents. Metab Clin Exp. 2008;57(8):1064–70.
Garenc C, Perusse L, Chagnon YC, Rankinen T, Gagnon J, Borecki IB, Leon AS, Skinner JS, Wilmore JH, Rao DC, Bouchard C. The alpha 2-adrenergic receptor gene and body fat content and distribution: the HERITAGE Family Study. Mol Med. 2002;8(2):88–94.
Joubert BR, Reif DM, Edwards SW, Leiner KA, Hudgens EE, Egeghy P, Gallagher JE, Hubal EC. Evaluation of genetic susceptibility to childhood allergy and asthma in an African American urban population. BMC Med Genet 2011, 12:25-2350-12-25.
Thomas PS. Tumour necrosis factor-&agr;: The role of this multifunctional cytokine in asthma. Immunol Cell Biol. 2001;79(2):132–40.
Joffe YT, Van der Merwe L, Evans J, Collins M, Lambert E, September A, Goedecke J. The tumor necrosis factor-α gene −238 G > A polymorphism, dietary fat intake, obesity risk and serum lipid concentrations in black and white South African women. Eur J Clin Nutr. 2012;66:1295–302.
Ichihara S, Yamada Y. Genetic factors for human obesity. Cell Mol Life Sci. 2008;65(7–8):1086–98.
Halapi E, Gudbjartsson DF, Jonsdottir GM, Bjornsdottir US, Thorleifsson G, Helgadottir H, Williams C, Koppelman GH, Heinzmann A, Boezen HM. A sequence variant on 17q21 is associated with age at onset and severity of asthma. Eur J Hum Genet. 2010;18(8):902–8.
Szczepankiewicz A, Breborowicz A, Sobkowiak P, Popiel A. Are genes associated with energy metabolism important in asthma and BMI? J Asthma. 2009;46(1):53–8.
Skaaby T, Husemoen LL, Martinussen T, Thyssen JP, Melgaard M, Thuesen BH, Pisinger C, Jørgensen T, Johansen JD, Menné T, Carlsen B, Szecsi PB, Stender S, Fenger RV, Fenger M, Linneberg A. Vitamin D status, filaggrin genotype, and cardiovascular risk factors: a Mendelian randomization approach. PLoS One. 2013;8(2):e57647. doi:10.1371/journal.pone.0057647.
Desai D, Brightling C. Cytokines and cytokine-specific therapy in asthma. Adv Clin Chem. 2012;57:59.
Thomas PS, Yates DH, Barnes PJ. Tumor necrosis factor-alpha increases airway responsiveness and sputum neutrophilia in normal human subjects. Am J Respir Crit Care Med. 1995;152(1):76–80.
Brightling C, Berry M, Amrani Y. Targeting TNF-α: a novel therapeutic approach for asthma. J Allergy Clin Immunol. 2008;121(1):5–10.
van Sijl AM, Peters MJL, Knol DL, de Vet RHC, Sattar N, Ben Dijkmans BAC, Smulders YM, Nurmohamed MT. The effect of TNF-alpha blocking therapy on lipid levels in rheumatoid arthritis: a meta-analysis. Semin Arthritis Rheum. 2011;41(3):393–400.
Glass CK, Olefsky JM. Inflammation and lipid signaling in the etiology of insulin resistance. Cell Metab. 2012;15(5):635–45.
Fontaine-Bisson B, El-Sohemy A. Genetic polymorphisms of tumor necrosis factor-alpha modify the association between dietary polyunsaturated fatty acids and plasma high-density lipoprotein-cholesterol concentrations in a population of young adults. J Nutrigenet Nutrigenomics. 2008;1(5):215–23.
Sobti RC, Kler R, Sharma YP, Talwar KK, Singh N. Risk of obesity and type 2 diabetes with tumor necrosis factor-α 308G/A gene polymorphism in metabolic syndrome and coronary artery disease subjects. Mol Cell Biochem. 2012;360(1–2):1–7.
Boraska V, Rayner NW, Groves CJ, Frayling TM, Diakite M, Rockett KA, Kwiatkowski DP, Day-Williams AG, McCarthy MI, Zeggini E. Large-scale association analysis of TNF/LTA gene region polymorphisms in type 2 diabetes. BMC Med Genet. 2010;11(1):69.
Yamagishi S, Adachi H, Matsui T, Nakamura K, Enomoto M, Fukami A, Otsuka M, Kumagae S, Nanjo Y, Imaizumi T. Decreased High-Density Lipoprotein Cholesterol Level is an Independent Correlate of Circulating Tumor Necrosis Factor-? in a General Population. Clin Cardiol. 2009;32(9):E29–32.
Yeh K, Lee C, Chang C, Lin Y, Huang J. Lipid Profiles Alter from Pro-Atherogenic into Less Atherogenic and Proinflammatory in Juvenile Idiopathic Arthritis Patients Responding to Anti TNF-α Treatment. PLoS One. 2014;9(3):e90757.
Di Minno MN, Ambrosino P, Peluso R, Di Minno A, Lupoli R, Dentali F. Lipid profile changes in patients with rheumatic diseases receiving a treatment with TNF-α blockers: A meta-analysis of prospective studies. Ann Med. 2014;46:73–83.
McGraw DW, Elwing JM, Fogel KM, Wang WC, Glinka CB, Mihlbachler KA, Rothenberg ME, Liggett SB. Crosstalk between Gi and Gq/Gs pathways in airway smooth muscle regulates bronchial contractility and relaxation. J Clin Invest. 2007;117(5):1391–8.
Artemenko Y, Gagnon A, Aubin D, Sorisky A. Anti‐adipogenic effect of PDGF is reversed by PKC inhibition. J Cell Physiol. 2005;204(2):646–53.
Schmitz‐Peiffer C. Protein kinase C and lipid‐induced insulin resistance in skeletal muscle. Ann N Y Acad Sci. 2002;967(1):146–57.
Arner P. Genetic variance and lipolysis regulation: implications for obesity. Ann Med. 2001;33(8):542–6.
Jalba MS, Rhoads GG, Demissie K. Association of Codon 16 and Codon 27 β2‐Adrenergic Receptor Gene Polymorphisms with Obesity: A Meta‐analysis. Obesity. 2008;16(9):2096–106.
Contopoulos-Ioannidis DG, Manoli EN, Ioannidis J. Meta-analysis of the association of β2-adrenergic receptor polymorphisms with asthma phenotypes. J Allergy Clin Immunol. 2005;115(5):963–72.
Thomsen M, Dahl M, Tybjærg‐Hansen A, Nordestgaard B. β2‐adrenergic receptor Thr164IIe polymorphism, blood pressure and ischaemic heart disease in 66 750 individuals. J Intern Med. 2012;271(3):305–14.
Thomsen M, Dahl M, Tybjærg-Hansen A, Nordestgaard BG. β2-adrenergic receptor Thr164Ile polymorphism, obesity, and diabetes: comparison with FTO, MC4R, and TMEM18 polymorphisms in more than 64,000 individuals. J Clin Endocrinol Metab. 2012;97(6):E1074–9.
Schurks M, Kurth T, Ridker PM, Buring JE, Zee RY. Association between polymorphisms in the beta2-adrenergic receptor gene with myocardial infarction and ischaemic stroke in women. Thromb Haemost. 2009;101(2):351–8.
Montella S, Baraldi E, Cazzato S, Aralla R, Berardi M, Brunetti LM, Cardinale F, Cutrera R, de Benedictis FM, di Palmo E, Di Pillo S, Fenu G, La Grutta S, Lombardi E, Piacentini G, Santamaria F, Ullmann N, Rusconi F, Italian Pediatric Severe Asthma Network (IPSAN) on behalf of the Italian Society of Pediatric Respiratory Diseases (SIMRI). Severe asthma features in children: a case–control online survey. Ital J Pediatr. 2016;42:9. doi:10.1186/s13052-016-0217-z.
The baseline health survey was under the auspices of the Ministries of Health, and Education and Culture of Cyprus.
The baseline health survey was funded by the Ministry of Health and Marfin Popular Bank (local bank). The follow-up evaluation was partly funded by (i) Merck Sharpe and Dohme Medical School Grant and (ii) Electricity Authority of Cyprus Grant. The assessments and work for this study were supported by the European Regional Development Fund and the Republic of Cyprus through the Cyprus Research Promotion Foundation (Project YΓEIA/ΔYΓEIA/0308(BIE)/22). The supporting sources had no involvement in study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the report for publication.
Availability of data and material
The dataset supporting the conclusions of this article is available in the Figshare repository. https://dx.doi.org/10.6084/m9.figshare.3462782.v1.
PY conceived, designed the study, secured the funding and prepared the first draft of the manuscript. PK performed the statistical analysis under NM supervision and contributed towards the preparation of the first draft of the manuscript. OK coordinated the field study. SY performed the genetic analysis under the guidance of PZ. KD and AZ performed the serum lipids assessment. SS, DP, NM and PZ contributed towards the interpretation of the findings. All authors read and approved the final manuscript.
This work was carried out at the Cyprus International Institute for Environmental & Public Health in Association with Harvard School of Public Health, Cyprus University of Technology, Limassol, Cyprus.
The author(s) declare that they have no competing interests.
Ethics approval and consent to participate
The study has been performed in accordance with the declaration of Helsinki and was approved by the Cyprus National Bioethics Committee (ref: EEBK/EP/2009/18). Written informed consent was obtained from each participant and their guardian.
Hardy – Weinberg equilibrium and MAF results. (DOCX 15 kb)
Genotype distribution of all 16 SNPs in Current Wheezers (CUWH) Vs Controls (NWNA) and HDL-C levels. (DOCX 28 kb)
SNPs genotype distribution in Active Asthmatics (ACAS) Vs Controls (NWNA) and HDL-C levels. (DOCX 16 kb)
Allelic comparisons for 16 SNPs in Never Wheezers Never Asthmatics (NWNA) Vs Current Wheezers (CUWH) and in Never Wheezers Never Asthmatics (NWNA) Vs Active Asthmatics (ACAS). (DOCX 21 kb)