Statistical learning theory for geospatial data. Case study: Aral sea

In recent years there has been an explosive growth in the development of adaptive and data driven methods. One of the efficient and data-driven approaches is based on statistical learning theory (Vapnik 1998). The theory is based on Structural Risk Minimisation (SRM) principle and has a solid statistical background. When applying SRM we are trying not only to reduce training error ? to fit the available data with a model, but also to reduce the complexity of the model and to reduce generalisation error. Many nonlinear learning procedures recently developed in neural networks and statistics can be understood and interpreted in terms of the structural risk minimisation inductive principle. A recent methodology based on SRM is called Support Vector Machines (SVM). At present SLT is still under intensive development and SVM find new areas of application (www.kernel-machines.org). SVM develop robust and non linear data models with excellent generalisation abilities that is very important both for monitoring and forecasting. SVM are extremely good when input space is high dimensional and training data set i not big enough to develop corresponding nonlinear model. Moreover, SVM use only support vectors to derive decision boundaries. It opens a way to sampling optimization, estimation of noise in data, quantification of data redundancy etc. Presentation of SVM for spatially distributed data is given in (Kanevski and Maignan 2004).

Prevalence of hyperuricemia and relation of serum uric acid with cardiovascular risk factors in a developing country.

BACKGROUND: The prevalence of hyperuricemia has rarely been investigated in developing countries. The purpose of the present study was to investigate the prevalence of hyperuricemia and the association between uric acid levels and the various cardiovascular risk factors in a developing country with high average blood pressures (the Seychelles, Indian Ocean, population mainly of African origin). METHODS: This cross-sectional health examination survey was based on a population random sample from the Seychelles. It included 1011 subjects aged 25 to 64 years. Blood pressure (BP), body mass index (BMI), waist circumference, waist-to-hip ratio, total and HDL cholesterol, serum triglycerides and serum uric acid were measured. Data were analyzed using scatterplot smoothing techniques and gender-specific linear regression models. RESULTS: The prevalence of a serum uric acid level >420 micromol/L in men was 35.2% and the prevalence of a serum uric acid level >360 micromol/L was 8.7% in women. Serum uric acid was strongly related to serum triglycerides in men as well as in women (r = 0.73 in men and r = 0.59 in women, p < 0.001). Uric acid levels were also significantly associated but to a lesser degree with age, BMI, blood pressure, alcohol and the use of antihypertensive therapy. In a regression model, triglycerides, age, BMI, antihypertensive therapy and alcohol consumption accounted for about 50% (R2) of the serum uric acid variations in men as well as in women. CONCLUSIONS: This study shows that the prevalence of hyperuricemia can be high in a developing country such as the Seychelles. Besides alcohol consumption and the use of antihypertensive therapy, mainly diuretics, serum uric acid is markedly associated with parameters of the metabolic syndrome, in particular serum triglycerides. Considering the growing incidence of obesity and metabolic syndrome worldwide and the potential link between hyperuricemia and cardiovascular complications, more emphasis should be put on the evolving prevalence of hyperuricemia in developing countries.

A Novel Acute Retroviral Syndrome Severity Score Predicts the Key Surrogate Markers for HIV-1 Disease Progression.

OBJECTIVE: Best long-term practice in primary HIV-1 infection (PHI) remains unknown for the individual. A risk-based scoring system associated with surrogate markers of HIV-1 disease progression could be helpful to stratify patients with PHI at highest risk for HIV-1 disease progression. METHODS: We prospectively enrolled 290 individuals with well-documented PHI in the Zurich Primary HIV-1 Infection Study, an open-label, non-randomized, observational, single-center study. Patients could choose to undergo early antiretroviral treatment (eART) and stop it after one year of undetectable viremia, to go on with treatment indefinitely, or to defer treatment. For each patient we calculated an a priori defined "Acute Retroviral Syndrome Severity Score" (ARSSS), consisting of clinical and basic laboratory variables, ranging from zero to ten points. We used linear regression models to assess the association between ARSSS and log baseline viral load (VL), baseline CD4+ cell count, and log viral setpoint (sVL) (i.e. VL measured ≥90 days after infection or treatment interruption). RESULTS: Mean ARSSS was 2.89. CD4+ cell count at baseline was negatively correlated with ARSSS (p = 0.03, n = 289), whereas HIV-RNA levels at baseline showed a strong positive correlation with ARSSS (p<0.001, n = 290). In the regression models, a 1-point increase in the score corresponded to a 0.10 log increase in baseline VL and a CD4+cell count decline of 12/µl, respectively. In patients with PHI and not undergoing eART, higher ARSSS were significantly associated with higher sVL (p = 0.029, n = 64). In contrast, in patients undergoing eART with subsequent structured treatment interruption, no correlation was found between sVL and ARSSS (p = 0.28, n = 40). CONCLUSION: The ARSSS is a simple clinical score that correlates with the best-validated surrogate markers of HIV-1 disease progression. In regions where ART is not universally available and eART is not standard this score may help identifying patients who will profit the most from early antiretroviral therapy.

Utility scores associated with the world health organisation drinking risk-level classification in alcohol dependence

BACKGROUND: Extensive research exists estimating the effect hazardous alcohol¦use on morbidity and mortality, but little research quantifies the association between¦alcohol consumption and utility scores in patients with alcohol dependence.¦In the context of comparative research, the World Health Organisation (WHO)¦proposed to categorise the risk for alcohol-related acute and chronic harm according¦to patients' average daily alcohol consumption. OBJECTIVES: To estimate utility¦scores associated with each category of the WHO drinking risk-level classification¦in patients with alcohol dependence (AD). METHODS: We used data from¦CONTROL, an observational cohort study including 143 AD patients from the Alcohol¦Treatment Center at Lausanne University Hospital, followed for 12 months.¦Average daily alcohol consumption was assessed monthly using the Timeline Follow-¦back method and patients were categorised according to the WHO drinking¦risk-level classification: abstinent, low, medium, high and very high. Other measures¦as sociodemographic characteristics and utility scores derived from the EuroQoL¦5-Dimensions questionnaire (EQ-5D) were collected every three months.¦Mixed models for repeated measures were used to estimate mean utility scores¦associated with WHO drinking risk-level categories. RESULTS: A total of 143 patients¦were included and the 12-month follow-up permitting the assessment of¦1318 person-months. At baseline the mean age of the patients was 44.6 (SD 11.8)¦and the majority of patients was male (63.6%). Using repeated measures analysis,¦utility scores decreased with increasing drinking levels, ranging from 0.80 in abstinent¦patients to 0.62 in patients with very high risk drinking level (p_0.0001).¦CONCLUSIONS: In this sample of patients with alcohol dependence undergoing¦specialized care, utility scores estimated from the EQ-5D appeared to substantially¦and consistently vary according to patients' WHO drinking level.

Molecular Proxies for Climate Maladaptation in a Long-Lived Tree (Pinus pinaster Aiton, Pinaceae).

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.

Evolution of spirituality and religiousness in chronic schizophrenia or schizo-affective disorders: A 3 years follow-up study

Purpose : Spirituality and religiousness have been shown to be highly prevalent in patients with schizophrenia. Religion can help instil a positive sense of self, decrease the impact of symptoms and provide social contacts. Religion may also be a source of suffering. In this context, this research explores whether religion remains stable over time. Methods : From an initial cohort of 115 out-patients, 80% completed the 3-years follow-up assessment. In order to study the evolution over time, a hierarchical cluster analysis using average linkage was performed on factorial scores at baseline and follow-up and their differences. A sensitivity analysis was secondarily performed to check if the outcome was influenced by other factors such as changes in mental states using mixed models. Results : Religion was stable over time for 63% patients; positive changes occurred for 20% (i.e., significant increase of religion as a resource or a transformation of negative religion to a positive one) and negative changes for 17% (i.e., decrease of religion as a resource or a transformation of positive religion to a negative one). Change in spirituality and/or religiousness was not associated with social or clinical status, but with reduced subjective quality of life and self-esteem; even after controlling for the influence of age, gender, quality of life and clinical factors at baseline. Conclusions : In this context of patients with chronic schizophrenia, religion appeared to be labile. Qualitative analyses showed that those changes expressed the struggles of patients and suggest that religious issues need to be discussed in clinical settings.

Plasma lipoprotein-associated phospholipase A2 activity in Alzheimer's disease, amnestic mild cognitive impairment, and cognitively healthy elderly subjects: a cross-sectional study.

ABSTRACT: INTRODUCTION: Lipoprotein-associated phospholipase A2 (Lp-PLA2) is a circulating enzyme with pro-inflammatory and oxidative activities associated with cardiovascular disease and ischemic stroke. While high plasma Lp-PLA2 activity was reported as a risk factor for dementia in the Rotterdam study, no association between Lp-PLA2 mass and dementia or Alzheimer's disease (AD) was detected in the Framingham study. The objectives of the current study were to explore the relationship of plasma Lp-PLA2 activity with cognitive diagnoses (AD, amnestic mild cognitive impairment (aMCI), and cognitively healthy subjects), cardiovascular markers, cerebrospinal fluid (CSF) markers of AD, and apolipoprotein E (APOE) genotype. METHODS: Subjects with mild AD (n = 78) and aMCI (n = 59) were recruited from the Memory Clinic, University Hospital, Basel, Switzerland; cognitively healthy subjects (n = 66) were recruited from the community. Subjects underwent standardised medical, neurological, neuropsychological, imaging, genetic, blood and CSF evaluation. Differences in Lp-PLA2 activity between the cognitive diagnosis groups were tested with ANOVA and in multiple linear regression models with adjustment for covariates. Associations between Lp-PLA2 and markers of cardiovascular disease and AD were explored with Spearman's correlation coefficients. RESULTS: There was no significant difference in plasma Lp-PLA2 activity between AD (197.1 (standard deviation, SD 38.4) nmol/min/ml) and controls (195.4 (SD 41.9)). Gender, statin use and low-density lipoprotein cholesterol (LDL) were independently associated with Lp-PLA2 activity in multiple regression models. Lp-PLA2 activity was correlated with LDL and inversely correlated with high-density lipoprotein (HDL). AD subjects with APOE-ε4 had higher Lp-PLA2 activity (207.9 (SD 41.2)) than AD subjects lacking APOE-ε4 (181.6 (SD 26.0), P = 0.003) although this was attenuated by adjustment for LDL (P = 0.09). No strong correlations were detected for Lp-PLA2 activity and CSF markers of AD. CONCLUSION: Plasma Lp-PLA2 was not associated with a diagnosis of AD or aMCI in this cross-sectional study. The main clinical correlates of Lp-PLA2 activity in AD, aMCI and cognitively healthy subjects were variables associated with lipid metabolism.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Predicting future distributions of mountain plants under climate change: does dispersal capacity matter?

Many studies have investigated the impacts that climate change could potentially have on the distribution of plant species, but few have attempted to constrain projections through plant dispersal limitations. Instead, most studies published so far have been using the simplification of considering dispersal as either unlimited or null. However, depending on a species' dispersal capacity, landscape fragmentation, and the rate of climatic change, these assumptions can lead to serious over- or underestimation of a species' future distribution. To quantify the discrepancies between unlimited, realistic, and no dispersal scenarios, we carried out projections of future distribution over the 21st century for 287 mountain plant species in a study area of the Western Swiss Alps. For each species, simulations were run for four dispersal scenarios (unlimited dispersal, no dispersal, realistic dispersal and realistic dispersal with long-distance dispersal events) and under four climate change scenarios. Although simulations accounting for realistic dispersal limitations did significantly differ from those considering dispersal as unlimited or null in terms of projected future distribution, using the unlimited dispersal simplification nevertheless provided good approximations for species extinctions under more moderate climate change scenarios. Overall, simulations accounting for dispersal limitations produced, for our mountainous study area, results that were significantly closer to unlimited dispersal than to no dispersal. Finally, analyzing the temporal pattern of species extinctions over the entire 21st century showed that, due to the possibility of a large number of species shifting their distribution to higher elevation, important species extinctions for our study area might not occur before the 2080-2100 time periods.

Influence of weather, rank, and home advantage on football outcomes in the gulf region.

PURPOSE: The objective of this study was to investigate the effects of weather, rank, and home advantage on international football match results and scores in the Gulf Cooperation Council (GCC) region. METHODS: Football matches (n = 2008) in six GCC countries were analyzed. To determine the weather influence on the likelihood of favorable outcome and goal difference, generalized linear model with a logit link function and multiple regression analysis were performed. RESULTS: In the GCC region, home teams tend to have greater likelihood of a favorable outcome (P < 0.001) and higher goal difference (P < 0.001). Temperature difference was identified as a significant explanatory variable when used independently (P < 0.001) or after adjustment for home advantage and team ranking (P < 0.001). The likelihood of favorable outcome for GCC teams increases by 3% for every 1-unit increase in temperature difference. After inclusion of interaction with opposition, this advantage remains significant only when playing against non-GCC opponents. While home advantage increased the odds of favorable outcome (P < 0.001) and goal difference (P < 0.001) after inclusion of interaction term, the likelihood of favorable outcome for a GCC team decreased (P < 0.001) when playing against a stronger opponent. Finally, the temperature and wet bulb globe temperature approximation were found as better indicators of the effect of environmental conditions than absolute and relative humidity or heat index on match outcomes. CONCLUSIONS: In GCC region, higher temperature increased the likelihood of a favorable outcome when playing against non-GCC teams. However, international ranking should be considered because an opponent with a higher rank reduced, but did not eliminate, the likelihood of a favorable outcome.

Direct Effect of Birth Weight on Childhood Blood Pressure: a Causal Mediation Analysis

Background: Numerous studies have shown a negative association between birth weight (BW) and blood pressure (BP) later in life. To estimate the direct effect of BW on BP, it is conventional to condition on current weight (CW). However, such conditioning can induce collider stratification bias in the estimate of the direct effect. Objective: To bound the potential bias due to U, an unmeasured common cause of CW and BP, on the estimate of the (controlled) direct effect of BW on BP. Methods: Data from a school based study in Switzerland were used (N = 4,005; 2,010 B/1,995 G; mean age: 12.3 yr [range: 10.1-14.9]). Measured common causes of BW-BP (SES, smoking, body weight, and hypertension status of the mother) and CW-BP (breastfeeding and child's physical activity and diet) were identified with DAGs. Linear regression models were fitted to estimate the association between BW and BP. Sensitivity analyses were conducted to assess the potential effect of U on the association between BW and BP. U was assumed 1) to be a binary variable that affected BP by the same magnitude in low BWand in normal BW children and 2) to have a different prevalence in low BW children and in normal BW children for a given CW. Results: A small negative association was observed between BW and BP [beta: -0.3 mmHg/kg (95% CI: -0.9 to 0.3)]. The association was strengthened upon conditioning for CW [beta: -1.5 mmHg/kg (95% CI: -2.1 to -0.9)]. Upon further conditioning on common causes of BW-BP and CW-BP, the association did not change substantially [beta: -1.4 mmHg/kg (95% CI: -2.0 to -0.8)]. The negative association could be explained by U only if U was strongly associated with BP and if there was a large difference in the prevalence of U between low BWand normal BW children. Conclusion: The observed negative association between BW and BP upon adjustment for CW was not easily explained by an unmeasured common cause of CWand BP.

Biomarkers of oxidative stress and its association with the urinary reducing capacity in bus maintenance workers.

BACKGROUND: Exposure to particles (PM) induces adverse health effects (cancer, cardiovascular and pulmonary diseases). A key-role in these adverse effects seems to be played by oxidative stress, which is an excess of reactive oxygen species relative to the amount of reducing species (including antioxidants), the first line of defense against reactive oxygen species. The aim of this study was to document the oxidative stress caused by exposure to respirable particles in vivo, and to test whether exposed workers presented changes in their urinary levels for reducing species.METHODS: Bus depot workers (n = 32) exposed to particles and pollutants (respirable PM4, organic and elemental carbon, particulate metal content, polycyclic aromatic hydrocarbons, NOx, O3) were surveyed over two consecutive days. We collected urine samples before and after each shift, and quantified an oxidative stress biomarker (8-hydroxy-2'-deoxyguanosine), the reducing capacity and a biomarker of PAH exposure (1-hydroxypyrene). We used a linear mixed model to test for associations between the oxidative stress status of the workers and their particle exposure as well as with their urinary level of reducing species.RESULTS: Workers were exposed to low levels of respirable PM4 (range 25-71 μg/m3). However, urinary levels of 8-hydroxy-2'-deoxyguanosine increased significantly within each shift and between both days for non-smokers. The between-day increase was significantly correlated (p < 0.001) with the concentrations of organic carbon, NOx, and the particulate copper content. The within-shift increase in 8OHdG was highly correlated to an increase of the urinary reducing capacity (Spearman ρ = 0.59, p < 0.0001).CONCLUSIONS: These findings confirm that exposure to components associated to respirable particulate matter causes a systemic oxidative stress, as measured with the urinary 8OHdG. The strong association observed between urinary 8OHdG with the reducing capacity is suggestive of protective or other mechanisms, including circadian effects. Additional investigations should be performed to understand these observations.

Correlation between microtensile bond strength data and clinical outcome of Class V restorations.

OBJECTIVE: To determine if the results of resin-dentin microtensile bond strength (µTBS) is correlated with the outcome parameters of clinical studies on non-retentive Class V restorations. METHODS: Resin-dentin µTBS data were obtained from one test center; the in vitro tests were all performed by the same operator. The µTBS testing was performed 8h after bonding and after 6 months of storing the specimens in water. Pre-test failures (PTFs) of specimens were included in the analysis, attributing them a value of 1MPa. Prospective clinical studies on cervical restorations (Class V) with an observation period of at least 18 months were searched in the literature. The clinical outcome variables were retention loss, marginal discoloration and marginal integrity. Furthermore, an index was formulated to be better able to compare the laboratory and clinical results. Estimates of adhesive effects in a linear mixed model were used to summarize the clinical performance of each adhesive between 12 and 36 months. Spearman correlations between these clinical performances and the µTBS values were calculated subsequently. RESULTS: Thirty-six clinical studies with 15 adhesive/restorative systems for which µTBS data were also available were included in the statistical analysis. In general 3-step and 2-step etch-and-rinse systems showed higher bond strength values than the 2-step/3-step self-etching systems, which, however, produced higher values than the 1-step self-etching and the resin modified glass ionomer systems. Prolonged water storage of specimens resulted in a significant decrease of the mean bond strength values in 5 adhesive systems (Wilcoxon, p<0.05). There was a significant correlation between µTBS values both after 8h and 6 months of storage and marginal discoloration (r=0.54 and r=0.67, respectively). However, the same correlation was not found between µTBS values and the retention rate, clinical index or marginal integrity. SIGNIFICANCE: As µTBS data of adhesive systems, especially after water storage for 6 months, showed a good correlation with marginal discoloration in short-term clinical Class V restorations, longitudinal clinical trials should explore whether early marginal staining is predictive for future retention loss in non-carious cervical restorations.

Associations of ambulatory blood pressure with urinary caffeine and caffeine metabolite excretions.

Intake of caffeinated beverages might be associated with reduced cardiovascular mortality possibly via the lowering of blood pressure. We estimated the association of ambulatory blood pressure with urinary caffeine and caffeine metabolites in a population-based sample. Families were randomly selected from the general population of Swiss cities. Ambulatory blood pressure monitoring was conducted using validated devices. Urinary caffeine, paraxanthine, theophylline, and theobromine excretions were measured in 24 hours urine using ultrahigh performance liquid chromatography tandem mass spectrometry. We used mixed models to explore the associations of urinary excretions with blood pressure although adjusting for major confounders. The 836 participants (48.9% men) included in this analysis had mean age of 47.8 and mean 24-hour systolic and diastolic blood pressure of 120.1 and 78.0 mm Hg. For each doubling of caffeine excretion, 24-hour and night-time systolic blood pressure decreased by 0.642 and 1.107 mm Hg (both P values <0.040). Similar inverse associations were observed for paraxanthine and theophylline. Adjusted night-time systolic blood pressure in the first (lowest), second, third, and fourth (highest) quartile of paraxanthine urinary excretions were 110.3, 107.3, 107.3, and 105.1 mm Hg, respectively (P trend <0.05). No associations of urinary excretions with diastolic blood pressure were generally found, and theobromine excretion was not associated with blood pressure. Anti-hypertensive therapy, diabetes mellitus, and alcohol consumption modify the association of caffeine urinary excretion with systolic blood pressure. Ambulatory systolic blood pressure was inversely associated with urinary excretions of caffeine and other caffeine metabolites. Our results are compatible with a potential protective effect of caffeine on blood pressure.

Genome-wide profiling of blood pressure in adults and children.

Hypertension is an important determinant of cardiovascular morbidity and mortality and has a substantial heritability, which is likely of polygenic origin. The aim of this study was to assess to what extent multiple common genetic variants contribute to blood pressure regulation in both adults and children and to assess overlap in variants between different age groups, using genome-wide profiling. Single nucleotide polymorphism sets were defined based on a meta-analysis of genome-wide association studies on systolic blood pressure and diastolic blood pressure performed by the Cohort for Heart and Aging Research in Genome Epidemiology (n=29 136), using different P value thresholds for selecting single nucleotide polymorphisms. Subsequently, genetic risk scores for systolic blood pressure and diastolic blood pressure were calculated in an independent adult population (n=2072) and a child population (n=1034). The explained variance of the genetic risk scores was evaluated using linear regression models, including sex, age, and body mass index. Genetic risk scores, including also many nongenome-wide significant single nucleotide polymorphisms, explained more of the variance than scores based only on very significant single nucleotide polymorphisms in adults and children. Genetic risk scores significantly explained ≤1.2% (P=9.6*10(-8)) of the variance in adult systolic blood pressure and 0.8% (P=0.004) in children. For diastolic blood pressure, the variance explained was similar in adults and children (1.7% [P=8.9*10(-10)] and 1.4% [P=3.3*10(-5)], respectively). These findings suggest the presence of many genetic loci with small effects on blood pressure regulation both in adults and children, indicating also a (partly) common polygenic regulation of blood pressure throughout different periods of life.