Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.

Developing a reliable strategy to infer the effective soil hydraulic properties from field evaporation experiments for agro-hydrological models

The Richards equation has been widely used for simulating soil water movement. However, the take-up of agro-hydrological models using the basic theory of soil water flow for optimizing irrigation, fertilizer and pesticide practices is still low. This is partly due to the difficulties in obtaining accurate values for soil hydraulic properties at a field scale. Here, we use an inverse technique to deduce the effective soil hydraulic properties, based on measuring the changes in the distribution of soil water with depth in a fallow field over a long period, subject to natural rainfall and evaporation using a robust micro Genetic Algorithm. A new optimized function was constructed from the soil water contents at different depths, and the soil water at field capacity. The deduced soil water retention curve was approximately parallel but higher than that derived from published pedo-tranfer functions for a given soil pressure head. The water contents calculated from the deduced soil hydraulic properties were in good agreement with the measured values. The reliability of the deduced soil hydraulic properties was tested in reproducing data measured from an independent experiment on the same soil cropped with leek. The calculation of root water uptake took account for both soil water potential and root density distribution. Results show that the predictions of soil water contents at various depths agree fairly well with the measurements, indicating that the inverse analysis is an effective and reliable approach to estimate soil hydraulic properties, and thus permits the simulation of soil water dynamics in both cropped and fallow soils in the field accurately. (C) 2009 Elsevier B.V. All rights reserved.

Genetic polymorphisms in the hypothalamic pathway in relation to subsequent weight change--the DiOGenes study

BACKGROUND: Single nucleotide polymorphisms (SNPs) in genes encoding the components involved in the hypothalamic pathway may influence weight gain and dietary factors may modify their effects. AIM: We conducted a case-cohort study to investigate the associations of SNPs in candidate genes with weight change during an average of 6.8 years of follow-up and to examine the potential effect modification by glycemic index (GI) and protein intake. METHODS AND FINDINGS: Participants, aged 20-60 years at baseline, came from five European countries. Cases ('weight gainers') were selected from the total eligible cohort (n = 50,293) as those with the greatest unexplained annual weight gain (n = 5,584). A random subcohort (n = 6,566) was drawn with the intention to obtain an equal number of cases and noncases (n = 5,507). We genotyped 134 SNPs that captured all common genetic variation across the 15 candidate genes; 123 met the quality control criteria. Each SNP was tested for association with the risk of being a 'weight gainer' (logistic regression models) in the case-noncase data and with weight gain (linear regression models) in the random subcohort data. After accounting for multiple testing, none of the SNPs was significantly associated with weight change. Furthermore, we observed no significant effect modification by dietary factors, except for SNP rs7180849 in the neuromedin β gene (NMB). Carriers of the minor allele had a more pronounced weight gain at a higher GI (P = 2 x 10⁻⁷). CONCLUSIONS: We found no evidence of association between SNPs in the studied hypothalamic genes with weight change. The interaction between GI and NMB SNP rs7180849 needs further confirmation.

Physical activity attenuates the body mass index-increasing influence of genetic variation in the FTO gene 1-3.

BACKGROUND: Intronic variation in the FTO (fat mass and obesity-associated) gene has been unequivocally associated with increased body mass index (BMI; in kg/m(2)) and the risk of obesity in populations of different ethnicity. OBJECTIVE: We examined whether this robust genetic predisposition to obesity can be attenuated by being more physically active. DESIGN: The FTO variant rs1121980 was genotyped in 20,374 participants (39-79 y of age) from the European Prospective Investigation into Cancer and Nutrition-Norfolk Study, an ethnically homogeneous population-based cohort. Physical activity (PA) was assessed with a validated self-reported questionnaire. The interaction between rs1121980 and PA on BMI and waist circumference (WC) was examined by including the interaction term in mixed-effect models. RESULTS: We confirmed that the risk (T) allele of rs1121980 was significantly associated with BMI (0.31-unit increase per allele; P < 0.001) and WC (0.77-cm increase per allele; P < 0.001). The PA level attenuated the effect of rs1121980 on BMI and WC; ie, whereas in active individuals the risk allele increased BMI by 0.25 per allele, the increase in BMI was significantly (P for interaction = 0.004) more pronounced (76%) in inactive individuals (0.44 per risk allele). We observed similar effects for WC (P for interaction = 0.02): the risk allele increased WC by 1.04 cm per allele in inactive individuals but by only 0.64 cm in active individuals. CONCLUSIONS: Our results showed that PA attenuates the effect of the FTO rs1121980 genotype on BMI and WC. This observation has important public health implications because we showed that a genetic susceptibility to obesity induced by FTO variation can be overcome, at least in part, by adopting a physically active lifestyle.

PPARGC1A sequence variation and cardiovascular risk-factor levels: a study of the main genetic effects and gene x environment interactions in children from the European Youth Heart Study.

AIMS/HYPOTHESIS: The PPARGC1A gene coactivates multiple nuclear transcription factors involved in cellular energy metabolism and vascular stasis. In the present study, we genotyped 35 tagging polymorphisms to capture all common PPARGC1A nucleotide sequence variations and tested for association with metabolic and cardiovascular traits in 2,101 Danish and Estonian boys and girls from the European Youth Heart Study, a multicentre school-based cross-sectional cohort study. METHODS: Fasting plasma glucose concentrations, anthropometric variables and blood pressure were measured. Habitual physical activity and aerobic fitness were objectively assessed using uniaxial accelerometry and a maximal aerobic exercise stress test on a bicycle ergometer, respectively. RESULTS: In adjusted models, nominally significant associations were observed for BMI (rs10018239, p = 0.039), waist circumference (rs7656250, p = 0.012; rs8192678 [Gly482Ser], p = 0.015; rs3755863, p = 0.02; rs10018239, beta = -0.01 cm per minor allele copy, p = 0.043), systolic blood pressure (rs2970869, p = 0.018) and fasting glucose concentrations (rs11724368, p = 0.045). Stronger associations were observed for aerobic fitness (rs7656250, p = 0.005; rs13117172, p = 0.008) and fasting glucose concentrations (rs7657071, p = 0.002). None remained significant after correcting for the number of statistical comparisons. We proceeded by testing for gene x physical activity interactions for the polymorphisms that showed nominal evidence of association in the main effect models. None of these tests was statistically significant. CONCLUSIONS/INTERPRETATION: Variants at PPARGC1A may influence several metabolic traits in this European paediatric cohort. However, variation at PPARGC1A is unlikely to have a major impact on cardiovascular or metabolic health in these children.

Genetic variation in the oxytocin receptor (OXTR) gene is associated with Asperger Syndrome

Background Autism Spectrum Conditions (ASC) are a group of neurodevelopmental conditions characterized by impairments in communication and social interaction, alongside unusually repetitive behaviors and narrow interests. ASC are highly heritable and have complex patterns of inheritance where multiple genes are involved, alongside environmental and epigenetic factors. Asperger Syndrome (AS) is a subgroup of these conditions, where there is no history of language or cognitive delay. Animal models suggest a role for oxytocin (OXT) and oxytocin receptor (OXTR) genes in social-emotional behaviors, and several studies indicate that the oxytocin/oxytocin receptor system is altered in individuals with ASC. Previous studies have reported associations between genetic variations in the OXTR gene and ASC. Methods The present study tested for an association between nine single nucleotide polymorphisms (SNPs) in the OXTR gene and AS in 530 individuals of Caucasian origin, using SNP association test and haplotype analysis. Results There was a significant association between rs2268493 in OXTR and AS. Multiple haplotypes that include this SNP (rs2268493-rs2254298, rs2268490-rs2268493-rs2254298, rs2268493-rs2254298-rs53576, rs237885-rs2268490-rs2268493-rs2254298, rs2268490-rs2268493-rs2254298-rs53576) were also associated with AS. rs2268493 has been previously associated with ASC and putatively alters several transcription factor-binding sites and regulates chromatin states, either directly or through other variants in linkage disequilibrium (LD). Conclusions This study reports a significant association of the sequence variant rs2268493 in the OXTR gene and associated haplotypes with AS.

Age and microenvironment outweigh genetic influence on the Zucker rat microbiome

Animal models are invaluable tools which allow us to investigate the microbiome-host dialogue. However, experimental design introduces biases in the data that we collect, also potentially leading to biased conclusions. With obesity at pandemic levels animal models of this disease have been developed; we investigated the role of experimental design on one such rodent model. We used 454 pyrosequencing to profile the faecal bacteria of obese (n = 6) and lean (homozygous n = 6; heterozygous n = 6) Zucker rats over a 10 week period, maintained in mixed-genotype cages, to further understand the relationships between the composition of the intestinal bacteria and age, obesity progression, genetic background and cage environment. Phylogenetic and taxon-based univariate and multivariate analyses (non-metric multidimensional scaling, principal component analysis) showed that age was the most significant source of variation in the composition of the faecal microbiota. Second to this, cage environment was found to clearly impact the composition of the faecal microbiota, with samples from animals from within the same cage showing high community structure concordance, but large differences seen between cages. Importantly, the genetically induced obese phenotype was not found to impact the faecal bacterial profiles. These findings demonstrate that the age and local environmental cage variables were driving the composition of the faecal bacteria and were more deterministically important than the host genotype. These findings have major implications for understanding the significance of functional metagenomic data in experimental studies and beg the question; what is being measured in animal experiments in which different strains are housed separately, nature or nurture?

A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.)

Faba bean (Vicia faba L.) is a globally important nitrogen-fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence-based consensus map was used to explore synteny with the most closely-related crop species, lentil, and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2 cM interval of Vf chromosome 2 which was highly collinear with Mt3. The obvious candidate gene from 77 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene (Mt3g092830, Mt3g092840) controlling anthocyanin production in Medicago and re-sequencing of the Vf orthologue showed a putative causative deletion of the entire 5’ end of the gene.

Genetic dissection of heterosis using epistatic association mapping in a partial NCII mating design

Heterosis refers to the phenomenon in which an F1 hybrid exhibits enhanced growth or agronomic performance. However, previous theoretical studies on heterosis have been based on bi-parental segregating populations instead of F1 hybrids. To understand the genetic basis of heterosis, here we used a subset of F1 hybrids, named a partial North Carolina II design, to perform association mapping for dependent variables: original trait value, general combining ability (GCA), specific combining ability (SCA) and mid-parental heterosis (MPH). Our models jointly fitted all the additive, dominance and epistatic effects. The analyses resulted in several important findings: 1) Main components are additive and additive-by-additive effects for GCA and dominance-related effects for SCA and MPH, and additive-by-dominant effect for MPH was partly identified as additive effect; 2) the ranking of factors affecting heterosis was dominance > dominance-by-dominance > over-dominance > complete dominance; and 3) increasing the proportion of F1 hybrids in the population could significantly increase the power to detect dominance-related effects, and slightly reduce the power to detect additive and additive-by-additive effects. Analyses of cotton and rapeseed datasets showed that more additive-by-additive QTL were detected from GCA than from trait phenotype, and fewer QTL were from MPH than from other dependent variables.