Search

Abstract Background Although smoking behavior is known to affect body mass index (BMI), the potential for smoking to influence genetic associations with BMI is largely unexplored. Methods As part of the ‘Population Architecture using Genomics and Epidemiology (PAGE)’ Consortium, we investigated interaction between genetic risk factors associated with BMI and smoking for 10 single nucleotide polymorphisms (SNPs) previously identified in genome-wide association studies. We included 6 studies with a total of 56,466 subjects (16,750 African Americans (AA) and 39,716 European Americans (EA)). We assessed effect modification by testing an interaction term for each SNP and smoking (current vs. former/never) in the linear regression and by stratified analyses. Results We did not observe strong evidence for interactions and only observed two interactions with p-values TMEM18, the risk allele (C) was associated with BMI only among AA females who were former/never smokers (β = 0.018, p = 0.002), vs. current smokers (β = 0.001, p = 0.95, pinteraction = 0.10). For rs9939609/FTO, the A allele was more strongly associated with BMI among current smoker EA females (β = 0.017, p = 3.5x10-5), vs. former/never smokers (β = 0.006, p = 0.05, pinteraction = 0.08). Conclusions These analyses provide limited evidence that smoking status may modify genetic effects of previously identified genetic risk factors for BMI. Larger studies are needed to follow up our results. Clinical Trial Registration NCT00000611

Abstract Background We sought to examine whether ε4 carrier status modifies the relation between body mass index (BMI) and HDL. The National Heart, Lung, and Blood Institute Family Heart Study included 657 families with high family risk scores for coronary heart disease and 588 randomly selected families of probands in the Framingham, Atherosclerosis Risk in Communities, and Utah Family Health Tree studies. We selected 1402 subjects who had ε4 carrier status available. We used generalized estimating equations to examine the interaction between BMI and ε4 allele carrier status on HDL after adjusting for age, gender, smoking, alcohol intake, mono- and poly-unsaturated fat intake, exercise, comorbidities, LDL, and family cluster. Results The mean (standard deviation) age of included subjects was 56.4(11.0) years and 47% were male. Adjusted means of HDL for normal, overweight, and obese BMI categories were 51.2(± 0.97), 45.0(± 0.75), and 41.6(± 0.93), respectively, among 397 ε4 carriers (p for trend < 0.0001) and 53.6(± 0.62), 51.3(± 0.49), and 45.0(± 0.62), respectively, among 1005 non-carriers of the ε4 allele (p-value for trend < 0.0001). There was no evidence for an interaction between BMI and ε4 status on HDL(p-value 0.39). Conclusion Our findings do not support an interaction between ε4 allele status and BMI on HDL.

Abstract Background The ferlin gene family possesses a rare and identifying feature consisting of multiple tandem C2 domains and a C-terminal transmembrane domain. Much currently remains unknown about the fundamental function of this gene family, however, mutations in its two most well-characterised members, dysferlin and otoferlin, have been implicated in human disease. The availability of genome sequences from a wide range of species makes it possible to explore the evolution of the ferlin family, providing contextual insight into characteristic features that define the ferlin gene family in its present form in humans. Results Ferlin genes were detected from all species of representative phyla, with two ferlin subgroups partitioned within the ferlin phylogenetic tree based on the presence or absence of a DysF domain. Invertebrates generally possessed two ferlin genes (one with DysF and one without), with six ferlin genes in most vertebrates (three DysF, three non-DysF). Expansion of the ferlin gene family is evident between the divergence of lamprey (jawless vertebrates) and shark (cartilaginous fish). Common to almost all ferlins is an N-terminal C2-FerI-C2 sandwich, a FerB motif, and two C-terminal C2 domains (C2E and C2F) adjacent to the transmembrane domain. Preservation of these structural elements throughout eukaryotic evolution suggests a fundamental role of these motifs for ferlin function. In contrast, DysF, C2DE, and FerA are optional, giving rise to subtle differences in domain topologies of ferlin genes. Despite conservation of multiple C2 domains in all ferlins, the C-terminal C2 domains (C2E and C2F) displayed higher sequence conservation and greater conservation of putative calcium binding residues across paralogs and orthologs. Interestingly, the two most studied non-mammalian ferlins (Fer-1 and Misfire) in model organisms C. elegans and D. melanogaster, present as outgroups in the phylogenetic analysis, with results suggesting reproduction-related divergence and specialization of species-specific functions within their genus. Conclusions Our phylogenetic studies provide evolutionary insight into the ferlin gene family. We highlight the existence of ferlin-like proteins throughout eukaryotic evolution, from unicellular phytoplankton and apicomplexan parasites, through to humans. We characterise the preservation of ferlin structural motifs, not only of C2 domains, but also the more poorly characterised ferlin-specific motifs representing the DysF, FerA and FerB domains. Our data suggest an ancient role of ferlin proteins, with lessons from vertebrate biology and human disease suggesting a role relating to vesicle fusion and plasma membrane specialization.

Abstract Background Transcription factor 7-like 2 (TCF7L2) has emerged as a consistently replicated susceptibility gene for type 2 diabetes, however, whether the TCF7L2 gene also has similar effects on the retinal microvasculature is less clear. We therefore aimed to investigate the association between the transcription factor 7-like 2 (TCF7L2) rs7903146 polymorphism and retinal microvascular phenotypes in the Atherosclerosis Risk in Communities (ARIC) Study (1993-1995). Methods This was a population-based, cross-sectional study of 10,320 middle-aged African American (n = 2,199) and Caucasian (n = 8,121) men and women selected from four United States communities to examine the association between TCF7L2 rs7903146 polymorphism and retinal microvascular signs (retinopathy, focal arteriolar narrowing, arteriovenous nicking, arteriolar and venular calibers). Photographs on one randomly selected eye were graded for presence of retinal microvascular signs and used to measure retinal vessel calibres. Results After adjusting for age, sex, study center, mean arterial blood pressure, total serum cholesterol, triglycerides, and other covariates, few associations of TCF7L2 rs7903146 and retinal microvascular signs were noted. TCF7L2 rs7903146 T risk allele was significantly associated with focal arteriolar narrowing in Caucasians with hypertension [odds ratio (OR)CT vs. CC (95% CI) = 1.25 (1.09-1.44); ORTT vs. CC = 1.56 (1.18-2.06); P = 0.002] and in Caucasians without diabetes [OR CT vs. CC = 1.18 (1.06-1.32); OR TT vs. CC = 1.40 (1.12, 1.75); P = 0.003]. No significant association of the TCF7L2 rs7903146 polymorphism and retinal vascular signs was noted among African American individuals. Conclusions TCF7L2 rs7903146 is not consistently associated with retinal microvascular signs. However, we report an association between the TCF7L2 rs7903146 polymorphism and focal arteriolar narrowing in Caucasians with hypertension or without diabetes. Further research in other large, population-based studies is needed to replicate these findings.

Abstract Background Osteoblasts are considered to primarily arise from osseous progenitors within the periosteum or bone marrow. We have speculated that cells from local soft tissues may also take on an osteogenic phenotype. Myoblasts are known to adopt a bone gene program upon treatment with the osteogenic bone morphogenetic proteins (BMP-2,-4,-6,-7,-9), but their osteogenic capacity relative to other progenitor types is unclear. We further hypothesized that the sensitivity of cells to BMP-2 would correlate with BMP receptor expression. Methods We directly compared the BMP-2 sensitivity of myoblastic murine cell lines and primary cells with osteoprogenitors from osseous tissues and fibroblasts. Fibroblasts forced to undergo myogenic conversion by transduction with a MyoD-expressing lentiviral vector (LV-MyoD) were also examined. Outcome measures included alkaline phosphatase expression, matrix mineralization, and expression of osteogenic genes (alkaline phosphatase, osteocalcin and bone morphogenetic protein receptor-1A) as measured by quantitative PCR. Results BMP-2 induced a rapid and robust osteogenic response in myoblasts and osteoprogenitors, but not in fibroblasts. Myoblasts and osteoprogenitors grown in osteogenic media rapidly upregulated Bmpr-1a expression. Chronic BMP-2 treatment resulted in peak Bmpr-1a expression at day 6 before declining, suggestive of a negative feedback mechanism. In contrast, fibroblasts expressed low levels of Bmpr-1a that was only weakly up-regulated by BMP-2 treatment. Bioinformatics analysis confirmed the presence of myogenic responsive elements in the proximal promoter region of human and murine BMPR-1A/Bmpr-1a. Forced myogenic gene expression in fibroblasts was associated with a significant increase in Bmpr-1a expression and a synergistic increase in the osteogenic response to BMP-2. Conclusion These data demonstrate the osteogenic sensitivity of muscle progenitors and provide a mechanistic insight into the variable response of different cell lineages to BMP-2.

Abstract Background Mannose-binding lectin (MBL) forms an integral part of the innate immune system. Persistent, subclinical infections and chronic inflammatory states are hypothesized to contribute to the pathogenesis of atherosclerosis. MBL gene (MBL2) variants with between 12 to 25% allele frequency in Caucasian and other populations, result in markedly reduced expression of functional protein. Prospective epidemiologic studies, including a nested, case-control study from the present population, have demonstrated the ability of MBL2 genotypes to predict complications of atherosclerosis,. The genetic control of MBL2 expression is complex and genetic background effects in specific populations are largely unknown. Methods The Strong Heart Study is a longitudinal, cohort study of cardiovascular disease among American Indians. A subset of individuals genotyped for the above mentioned case-control study were selected for analysis of circulating MBL levels by double sandwich ELISA method. Mean MBL levels were compared between genotypic groups and multivariate regression was used to determine other independent factors influencing MBL2 expression. Results Our results confirm the effects of variant structural (B, C, and D) and promoter (H and Y) alleles that have been seen in other populations. In addition, MBL levels were found to be positively associated with male gender and hemoglobin A1c levels, but inversely related to triglyceride levels. Correlation was not found between MBL and other markers of inflammation. Conclusion New data is presented concerning the effects of known genetic variants on MBL levels in an American Indian population, as well as the relationship of MBL2 expression to clinical and environmental factors, including inflammatory markers.

Abstract Background In a previous study of the Hypertension Genetic Epidemiology Network (HyperGEN) we have shown that metabolic syndrome (MetS) risk factors were moderately and significantly associated with echocardiographic (ECHO) left ventricular (LV) phenotypes. Methods The study included 1,393 African Americans and 1,133 whites, stratified by type 2 diabetes mellitus (DM) status. Heritabilities of seven factor scores based on the analysis of 15 traits were sufficiently high to pursue QTL discovery in this follow-up study. Results Three of the QTLs discovered relate to combined MetS-ECHO factors of "blood pressure (BP)-LV wall thickness" on chromosome 3 at 225 cM with a 2.8 LOD score, on chromosome 20 at 2.1 cM with a 2.6 LOD score; and for "LV wall thickness" factor on chromosome 16 at 113.5 with a 2.6 LOD score in whites. The remaining QTLs include one for a "body mass index-insulin (BMI-INS)" factor with a LOD score of 3.9 on chromosome 2 located at 64.8 cM; one for the same factor on chromosome 12 at 91.4 cM with a 3.3 LOD score; one for a "BP" factor on chromosome 19 located at 67.8 cM with a 3.0 LOD score. A suggestive linkage was also found for "Lipids-INS" with a 2.7 LOD score located on chromosome 11 at 113.1 cM in African Americans. Of the above QTLs, the one on chromosome 12 for "BMI-INS" is replicated in both ethnicities, (with highest LOD scores in African Americans). In addition, the QTL for "LV wall thickness" on chromosome 16q24.2-q24.3 reached its local maximum LOD score at marker D16S402, which is positioned within the 5th intron of the cadherin 13 gene, implicated in heart and vascular remodeling. Conclusion Our previous study and this follow-up suggest gene loci for some crucial MetS and cardiac geometry risk factors that contribute to the risk of developing heart disease.

Abstract Background Body fat mass distribution and deposition are determined by multiple environmental and genetic factors. Obesity is associated with insulin resistance, hyperinsulinemia, and type 2 diabetes. We previously identified evidence for genotype-by-diabetes interaction on obesity traits in Strong Heart Family Study (SHFS) participants. To localize these genetic effects, we conducted genome-wide linkage scans of obesity traits in individuals with and without type 2 diabetes, and in the combined sample while modeling interaction with diabetes using maximum likelihood methods (SOLAR 2.1.4). Methods SHFS recruited American Indians from Arizona, North and South Dakota, and Oklahoma. Anthropometric measures and diabetes status were obtained during a clinic visit. Marker allele frequencies were derived using maximum likelihood methods estimated from all individuals and multipoint identity by descent sharing was estimated using Loki. We used variance component linkage analysis to localize quantitative trait loci (QTLs) influencing obesity traits. We tested for evidence of additive and QTL-specific genotype-by-diabetes interactions using the regions identified in the diabetes-stratified analyses. Results Among 245 diabetic and 704 non-diabetic American Indian individuals, we detected significant additive gene-by-diabetes interaction for weight and BMI (P < 0.02). In analysis accounting for QTL-specific interaction (P < 0.001), we detected a QTL for weight on chromosome 1 at 242 cM (LOD = 3.7). This chromosome region harbors the adiponectin receptor 1 gene, which has been previously associated with obesity. Conclusion These results suggest distinct genetic effects on body mass in individuals with diabetes compared to those without diabetes, and a possible role for one or more genes on chromosome 1 in the pathogenesis of obesity.

Abstract Background The pathogenesis of atherosclerosis involves both hemostatic and inflammatory mechanisms. Fibrinogen is associated with both risk of thrombosis and inflammation. A recent meta-analysis showed that risk of coronary heart disease may increase 1.8 fold for 1 g/L of increased fibrinogen, independent of traditional risk factors. It is known that fibrinogen levels may be influenced by demographic, environmental and genetic factors. Epidemiologic and candidate gene studies are available; but few genome-wide linkage studies have been conducted, particularly in minority populations. The Strong Heart Study has demonstrated an increased incidence of cardiovascular disease in the American Indian population, and therefore represents an important source for genetic-epidemiological investigations. Methods The Strong Heart Family Study enrolled over 3,600 American Indian participants in large, multi-generational families, ascertained from an ongoing population-based study in the same communities. Fibrinogen was determined using standard technique in a central laboratory and extensive additional phenotypic measures were obtained. Participants were genotyped for 382 short tandem repeat markers distributed throughout the genome; and results were analyzed using a variance decomposition method, as implemented in the SOLAR 2.0 program. Results Data from 3535 participants were included and after step-wise, linear regression analysis, two models were selected for investigation. Basic demographic adjustments constituted model 1, while model 2 considered waist circumference, diabetes mellitus and postmenopausal status as additional covariates. Five LOD scores between 1.82 and 3.02 were identified, with the maximally adjusted model showing the highest score on chromosome 7 at 28 cM. Genes for two key components of the inflammatory response, i.e. interleukin-6 and "signal transducer and activator of transcription 3" (STAT3), were identified within 2 and 8 Mb of this 1 LOD drop interval respectively. A LOD score of 1.82 on chromosome 17 between 68 and 93 cM is supported by reports from two other populations with LOD scores of 1.4 and 1.95. Conclusion In a minority population with a high prevalence of cardiovascular disease, strong evidence for a novel genetic determinant of fibrinogen levels is found on chromosome 7 at 28 cM. Four other loci, some of which have been suggested by previous studies, were also identified.

Abstract Background Data accumulated from mouse studies and in vitro studies of human arteries support the notion that soluble intercellular adhesion molecule-1 (sICAM-1) and monocyte chemoattractant protein-1 (MCP-1) play important roles in the inflammation process involved in atherosclerosis. However, at the population level, the utility of sICAM-1 and MCP-1 as biomarkers for subclinical atherosclerosis is less clear. In the follow-up exam of the NHLBI Family Heart Study, we evaluated whether plasma levels of sICAM-1 and MCP-1 were associated with coronary artery calcification (CAC), a measure of the burden of coronary atherosclerosis. Methods CAC was measured using the Agatston score with multidetector computed tomography. Information on CAC and MCP-1 was obtained in 2246 whites and 470 African Americans (mean age 55 years) without a history of coronary heart disease (CHD). Information on sICAM-1 was obtained for white participants only. Results In whites, after adjustment for age and gender, the odds ratios (ORs) of CAC (CAC > 0) associated with the second, third, fourth, and fifth quintiles of sICAM-1 compared to the first quintile were 1.22 (95% confidence interval [CI]: 0.91–1.63), 1.15 (0.84–1.58), 1.49 (1.09–2.05), and 1.72 (1.26–2.36) (p = 0.0005 for trend test), respectively. The corresponding ORs for the second to fifth quintiles of MCP-1 were 1.26 (0.92–1.73), 0.99 (0.73–1.34), 1.42 (1.03–1.96), and 2.00 (1.43–2.79) (p < 0.0001 for trend test), respectively. In multivariable analysis that additionally adjusted for other CHD risk factors, the association of CAC with sICAM-1 and MCP-1 was attenuated and no longer statistically significant. In African Americans, the age and gender-adjusted ORs of CAC associated with the second and third tertiles of MCP-1 compared to the first tertile were 1.16 (0.64–2.08) and 1.25 (0.70–2.23) (p = 0.44 for trend test), respectively. This result did not change materially after additional adjustment for other CHD risk factors. Test of race interaction showed that the magnitude of association between MCP-1 and CAC did not differ significantly between African Americans and whites. Similar results were obtained when CAC ≥ 10 was analyzed as an outcome for both MCP-1 and sICAM-1. Conclusion This study suggests that sICAM-1 and MCP-1 are biomarkers of coronary atherosclerotic burden and their association with CAC was mainly driven by established CHD risk factors.

Abstract Background We report longitudinal changes in the metabolic syndrome (MetS) in 2,458 participants from 480 families in the Family Heart Study. Participants were examined between 1994–96 (FHS-T1) and 2002–03 (FHS-T2), about 7.4 years apart. Additionally, the impact of medication on estimates of MetS prevalence, and associations of MetS with prevalent coronary heart disease (CHD) and type 2 diabetes (T2D) were studied. Methods Three definitions for MetS prevalence were considered. One represented the original (o) National Cholesterol Education Program (NCEP) MetS criteria. Two others considered the confounding of medications effects, respectively (m) lipid medications constituted a categorical diagnostic criterion for lipids variables, and (c) lipids and blood pressure variables were corrected with average clinical trials medications effects. Logistic regression of MetS on CHD and T2D, as well as the trend analysis of MetS by age, were performed. Results MetS increased from 17.1% in FHS-T1(o) to 28.8% in FHS-T2(o); from 19.7% in FHS-T1(m) to 42.5% in FHS-T2(m); and from 18.4% in FHS-T1(c) to 33.6% in FHS-T2(c). While we observed adverse changes in all risk factors, the greatest increase was for waist circumference (25%). The percentages of MetS were about 2 to almost 3 times higher in ages 50 years and older than in younger ages. The odds of having prevalent CHD were about 2.5 times higher in the subjects classified with MetS than without. Conclusion MetS percentages increased noticeably longitudinally and cross-sectionally with older age. These conclusions were reached with and without considering medication use, but correcting risk factors for medications use affects the MetS prevalence estimates. As found in other studies, MetS was associated with increased odds for prevalent CHD.

Changes in gene expression form a key component of the molecular mechanisms by which plants adapt and respond to environmental stresses. There is compelling evidence for the role of stimulus-specific Ca(2+) signatures in plant stress responses. However, our understanding of how they orchestrate the differential expression of stress-induced genes remains fragmentary. We have undertaken a global study of changes in the Arabidopsis transcriptome induced by the pollutant ozone in order to establish a robust transcriptional response against which to test the ability of Ca(2+) signatures to encode stimulus-specific transcriptional information. We show that the expression of a set of co-regulated ozone-induced genes is Ca(2+)-dependent and that abolition of the ozone-induced Ca(2+) signature inhibits the induction of these genes by ozone. No induction of this set of ozone-regulated genes was observed in response to H(2)O(2), one of the reactive oxygen species (ROS) generated by ozone, or cold stress, which also generates ROS, both of which stimulate changes in [Ca(2+)](cyt). These data establish unequivocally that the Ca(2+)-dependent changes in gene expression observed in response to ozone are not simply a consequence of an ROS-induced increase in [Ca(2+) ](cyt) per se. The magnitude and temporal dynamics of the ozone, H(2)O(2) , and cold Ca(2+) signatures all differ markedly. This finding is consistent with the hypothesis that stimulus-specific transcriptional information can be encoded in the spatiotemporal dynamics of complex Ca(2+) signals in plants., (© 2012 The Authors. The Plant Journal © 2012 Blackwell Publishing Ltd.)

Improving nutrient use efficiency of crop plants, especially at low input, is essential to ensure sustainable food production in the future. In order to address the genetic basis of nutrient use efficiency in a model system, growth of Arabidopsis ecotypes at normal and low nitrogen (N) supply was compared. The ecotypes differed significantly in the extent of growth reduction in limiting conditions. The fresh weight of Shahdara and Ws grown at 1mM nitrate was reduced by 30% compared to control, whereas Col-0 and Ga-0 were almost unaffected. Total N content was reduced in all ecotypes by 10-30%. The capacity to store nitrate correlated with the tolerance to low N; in Shahdara and Ws, but not in Col-0 and Ga-0, nitrate content on low N was significantly reduced compared to control nutrition. The mRNA levels for genes of nitrate uptake and assimilation were only moderately affected by the treatment. The transcript levels of nitrate reductase NIA1 and nitrite reductase were higher in the ecotypes tolerant to low N (Col-0 and Ga-0) with normal N nutrition but on low N they were reduced to a much higher extent than the sensitive ecotypes (Shahdara and Ws). It seems that a higher capacity to keep nitrate reserves at low N, perhaps due to the ability to turn down nitrate reduction rate, is responsible for a better tolerance of Col-0 and Ga-0 to low N supply.

Sulfur-containing compounds play an important role in plant stress defense; however, only a little is known about the molecular mechanisms of regulation of sulfate assimilation by stress. Using known Arabidopsis (Arabidopsis thaliana) mutants in signaling pathways, we analyzed regulation of the key enzyme of sulfate assimilation, adenosine 5'-phosphosulfate reductase (APR), by salt stress. APR activity and mRNA levels of all three APR isoforms increased 3-fold in roots after 5 h of treatment with 150 mm NaCl. The regulation of APR was not affected in mutants deficient in abscisic acid (ABA) synthesis and treatment of the plants with ABA did not affect the mRNA levels of APR isoforms, showing that APR is regulated by salt stress in an ABA-independent manner. In mutants deficient in jasmonate, salicylate, or ethylene signaling, APR mRNA levels were increased upon salt exposure similar to wild-type plants. Surprisingly, however, APR enzyme activity was not affected by salt in these plants. The same result was obtained in mutants affected in cytokinin and auxin signaling. Signaling via gibberellic acid, on the other hand, turned out to be essential for the increase in APR mRNA by salt treatment. These results demonstrate an extensive posttranscriptional regulation of plant APR and reveal that the sulfate assimilation pathway is controlled by a complex network of multiple signals on different regulatory levels.

Most agronomic traits of importance, whether physiological (such as nutrient use efficiency) or developmental (such as flowering time), are controlled simultaneously by multiple genes and their interactions with the environment. Here, we show that variation in sulfate content between wild Arabidopsis thaliana accessions Bay-0 and Shahdara is controlled by a major quantitative trait locus that results in a strong interaction with nitrogen availability in the soil. Combining genetic and biochemical results and using a candidate gene approach, we have cloned the underlying gene, showing how a single-amino acid substitution in a key enzyme of the assimilatory sulfate reduction pathway, adenosine 5'-phosphosulfate reductase, is responsible for a decrease in enzyme activity, leading to sulfate accumulation in the plant. This work illustrates the potential of natural variation as a source of new alleles of known genes, which can aid in the study of gene function and metabolic pathway regulation. Our new insights on sulfate assimilation may have an impact on sulfur fertilizer use and stress defense improvement.

Stomatal responses to air pollutants are complex, varying among species and with concentration, environmental conditions and age. In general, short-term exposure to sulphur dioxide (SO 2 ) promotes stomatal opening, whereas longer-term exposure can cause partial stomatal closure. By contrast, the effects of oxides of nitrogen (NO x ) are often small or insignificant. The effects of ozone, and oxidative stress, are equally complex. Short-term exposure to ozone stimulates a rapid reduction in stomatal aperture, whilst longer-term exposure causes stomatal responses to become sluggish. The response of stomata to abscisic acid (ABA) has been shown to be slower in plants exposed to a combination of SO 2 and NO 2 suggesting an adverse effect on guard cell ABA signal transduction. In addition, ozone causes a reduction in stomatal closure under drought conditions. There is an increasing body of evidence to suggest that air pollutants and oxidative stresses can have a marked effect on the Ca 2+ homeostasis of guard cells and the intracellular machinery responsible for stomatal movements. Here we discuss the effects of air pollutants on stomatal responses and their possible effects on Ca 2+ based signalling in guard cells focusing on the effects of ozone and oxidative stress.