Genome-wide mapping of 5-hydroxymethylcytosine in three rice cultivars reveals its preferential localization in transcriptionally silent transposable element genes

5-Hydroxymethylcytosine (5hmC), a modified form of cytosine that is considered the sixth nucleobase in DNA, has been detected in mammals and is believed to play an important role in gene regulation. In this study, 5hmC modification was detected in rice by employing a dot-blot assay, and its levels was further quantified in DNA from different rice tissues using liquid chromatography-multistage mass spectrometry (LC-MS/MS/MS). The results showed large intertissue variation in 5hmC levels. The genome-wide profiles of 5hmC modification in three different rice cultivars were also obtained using a sensitive chemical labelling followed by a next-generation sequencing method. Thousands of 5hmC peaks were identified, and a comparison of the distributions of 5hmC among different rice cultivars revealed the specificity and conservation of 5hmC modification. The identified 5hmC peaks were significantly enriched in heterochromatin regions,and mainly located in transposable element (TE) genes, especially around retrotransposons. The correlation analysis of 5hmC and gene expression data revealed a close association between 5hmC and silent TEs. These findings provide a resource for plant DNA 5hmC epigenetic studies and expand our knowledge of 5hmC modification.

Female-specific intergenerational transmission patterns of the human corticolimbic circuitry

Parents have large genetic and environmental influences on offspring’s cognition, behavior, and brain. These intergenerational effects are observed in mood disorders, with particularly robust association in depression between mothers and daughters. No studies have thus far examined the neural bases of these intergenerational effects in humans. Corticolimbic circuitry is known to be highly relevant in a wide range of processes including mood regulation and depression. These findings suggest that corticolimbic circuitry may also show matrilineal transmission patterns. We therefore examined human parent-offspring association in this neurocircuitry, and investigated the degree of association in gray matter volume between parent and offspring. We used voxel-wise correlation analysis in a total of 35 healthy families, consisting of parents and their biological offspring. We found positive associations of regional grey matter volume in the corticolimbic circuit including the amygdala, hippocampus, anterior cingulate cortex, and ventromedial prefrontal cortex between biological mothers and daughters. This association was significantly greater than mother-son, father-daughter, and father-son associations. The current study suggests that the corticolimbic circuitry, which has been implicated in mood regulation, shows a matrilineal specific transmission patterns. Our preliminary findings are consistent with what has been found behaviorally in depression, and may have clinical implications for disorders known to have dysfunction in mood regulation such as depression. Studies such as ours will likely bridge animal work examining gene expression in the brains and clinical symptom-based observations, and provide promising ways to investigate intergenerational transmission patterns in the human brain.

Annular modes and apparent eddy feedbacks in the Southern Hemisphere

Lagged correlation analysis is often used to infer intraseasonal dynamical effects but is known to be affected by non-stationarity. We highlight a pronounced quasi-two-year peak in the anomalous zonal wind and eddy momentum flux convergence power spectra in the Southern Hemisphere, which is prima facie evidence for non-stationarity. We then investigate the consequences of this non-stationarity for the Southern Annular Mode and for eddy momentum flux convergence. We argue that positive lagged correlations previously attributed to the existence of an eddy feedback are more plausibly attributed to non-stationary interannual variability external to any potential feedback process in the mid-latitude troposphere. The findings have implications for the diagnosis of feedbacks in both models and re-analysis data as well as for understanding the mechanisms underlying variations in the zonal wind.

Impact of increasing fruit and vegetables and flavonoid intake on the human gut microbiota

Epidemiological studies have shown protective effects of fruits and vegetables (F&V) in lowering the risk of developing cardiovascular diseases (CVD) and cancers. Plant-derived dietary fibre (non-digestible polysaccharides) and/or flavonoids may mediate the observed protective effects particularly through their interaction with the gut microbiota. The aim of this study was to assess the impact of fruit and vegetable (F&V) intake on gut microbiota, with an emphasis on the role of flavonoids, and further to explore relationships between microbiota and factors associated with CVD risk. In the study, a parallel design with 3 study groups, participants in the two intervention groups representing high-flavonoid (HF) and low flavonoid (LF) intakes were asked to increase their daily F&V intake by 2, 4 and 6 portions for a duration of 6 weeks each, while a third (control) group continued with their habitual diet. Faecal samples were collected at baseline and after each dose from 122 subjects. Faecal bacteria enumeration was performed by fluorescence in situ hybridisation (FISH). Correlations of dietary components, flavonoid intake and markers of CVD with bacterial numbers were also performed. A significant dose X treatment interaction was only found for Clostidium leptum-Ruminococcus bromii/flavefaciens with a significant increase after intake of 6 additional portions in the LF group. Correlation analysis of the data from all 122 subjects independent from dietary intervention indicated an inhibitory role of F&V intake, flavonoid content and sugars against the growth of potentially pathogenic clostridia. Additionally, we observed associations between certain bacterial populations and CVD risk factors including plasma TNF-α, plasma lipids and BMI/waist circumference.

Correlation between process cheese meltability determined by sensory analysis, computer vision method and Olson and Price test

The meltabilities of 14 process cheese samples were determined at 2 and 4 weeks after manufacture using sensory analysis, a computer vision method, and the Olson and Price test. Sensory analysis meltability correlated with both computer vision meltability (R-2 = 0.71, P < 0.001) and Olson and Price meltability (R-2 = 0.69, P < 0.001). There was a marked lack of correlation between the computer vision method and the Olson and Price test. This study showed that the Olson and Price test gave greater repeatability than the computer vision method. Results showed process cheese meltability decreased with increasing inorganic salt content and with lower moisture/fat ratios. There was very little evidence in this study to show that process cheese meltability changed between 2 and 4 weeks after manufacture..

Identifying the critical parameters of a cyanobacterial growth and movement model by using generalised sensitivity analysis

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling of cyanobacteria in freshwaters is an important tool for understanding their population dynamics and predicting the location and timing of the bloom events in lakes and rivers. A new deterministic-mathematical model was developed, which simulates the growth and movement of cyanobacterial blooms in river systems. The model focuses on the mathematical description of the bloom formation, vertical migration and lateral transport of colonies within river environments by taking into account the major factors that affect the cyanobacterial bloom formation in rivers including, light, nutrients and temperature. A technique called generalised sensitivity analysis was applied to the model to identify the critical parameter uncertainties in the model and investigates the interaction between the chosen parameters of the model. The result of the analysis suggested that 8 out of 12 parameters were significant in obtaining the observed cyanobacterial behaviour in a simulation. It was found that there was a high degree of correlation between the half-saturation rate constants used in the model.

Spatial analysis of the error in a model of soil nitrogen

Models of the dynamics of nitrogen in soil (soil-N) can be used to aid the fertilizer management of a crop. The predictions of soil-N models can be validated by comparison with observed data. Validation generally involves calculating non-spatial statistics of the observations and predictions, such as their means, their mean squared-difference, and their correlation. However, when the model predictions are spatially distributed across a landscape the model requires validation with spatial statistics. There are three reasons for this: (i) the model may be more or less successful at reproducing the variance of the observations at different spatial scales; (ii) the correlation of the predictions with the observations may be different at different spatial scales; (iii) the spatial pattern of model error may be informative. In this study we used a model, parameterized with spatially variable input information about the soil, to predict the mineral-N content of soil in an arable field, and compared the results with observed data. We validated the performance of the N model spatially with a linear mixed model of the observations and model predictions, estimated by residual maximum likelihood. This novel approach allowed us to describe the joint variation of the observations and predictions as: (i) independent random variation that occurred at a fine spatial scale; (ii) correlated random variation that occurred at a coarse spatial scale; (iii) systematic variation associated with a spatial trend. The linear mixed model revealed that, in general, the performance of the N model changed depending on the spatial scale of interest. At the scales associated with random variation, the N model underestimated the variance of the observations, and the predictions were correlated poorly with the observations. At the scale of the trend, the predictions and observations shared a common surface. The spatial pattern of the error of the N model suggested that the observations were affected by the local soil condition, but this was not accounted for by the N model. In summary, the N model would be well-suited to field-scale management of soil nitrogen, but suited poorly to management at finer spatial scales. This information was not apparent with a non-spatial validation. (c),2007 Elsevier B.V. All rights reserved.

Predictions of the arrival time of coronal mass ejections at 1AU: an analysis of the causes of errors

Three existing models of Interplanetary Coronal Mass Ejection (ICME) transit between the Sun and the Earth are compared to coronagraph and in situ observations: all three models are found to perform with a similar level of accuracy (i.e. an average error between observed and predicted 1AU transit times of approximately 11 h). To improve long-term space weather prediction, factors influencing CME transit are investigated. Both the removal of the plane of sky projection (as suffered by coronagraph derived speeds of Earth directed CMEs) and the use of observed values of solar wind speed, fail to significantly improve transit time prediction. However, a correlation is found to exist between the late/early arrival of an ICME and the width of the preceding sheath region, suggesting that the error is a geometrical effect that can only be removed by a more accurate determination of a CME trajectory and expansion. The correlation between magnetic field intensity and speed of ejecta at 1AU is also investigated. It is found to be weak in the body of the ICME, but strong in the sheath, if the upstream solar wind conditions are taken into account.

Genetic parameters and evaluations from single- and multiple-trait analysis of dairy cow fertility and milk production

Genetic parameters and breeding values for dairy cow fertility were estimated from 62 443 lactation records. Two-trait analysis of fertility and milk yield was investigated as a method to estimate fertility breeding values when culling or selection based on milk yield in early lactation determines presence or absence of fertility observations in later lactations. Fertility traits were calving interval, intervals from calving to first service, calving to conception and first to last service, conception success to first service and number of services per conception. Milk production traits were 305-day milk, fat and protein yield. For fertility traits, range of estimates of heritability (h(2)) was 0.012 to 0.028 and of permanent environmental variance (c(2)) was 0.016 to 0.032. Genetic correlations (r(g)) among fertility traits were generally high ( > 0.70). Genetic correlations of fertility with milk production traits were unfavourable (range -0.11 to 0.46). Single and two-trait analyses of fertility were compared using the same data set. The estimates of h(2) and c(2) were similar for two types of analyses. However, there were differences between estimated breeding values and rankings for the same trait from single versus multi-trait analyses. The range for rank correlation was 0.69-0.83 for all animals in the pedigree and 0.89-0.96 for sires with more than 25 daughters. As single-trait method is biased due to selection on milk yield, a multi-trait evaluation of fertility with milk yield is recommended. (C) 2002 Elsevier Science B.V. All rights reserved.

Juvenile-mature wood transition in pine: correlation between wood properties and candidate gene expression profiles

There is a strong desire to exploit transcriptomics data from model species for the genetic improvement of non-model crops. Here, we use gene expression profiles from the commercial model Pinus taeda to identify candidate genes implicated in juvenile-mature wood transition in the non-model relative, P. sylvestris. Re-analysis of 'public domain' SAGE data from xylem tissues of P. taeda revealed 283 mature-abundant and 396 juvenile-abundant tags (P < 0.01), of which 70 and 137, respectively matched to genes with known function. Based on sequence similarity, we then isolated 16 putative homologues of genes that in P. taeda exhibited widest divergence in expression between juvenile and mature samples. Candidate expression levels in P. sylvestris were almost invariably differential between juvenile and mature woody tissue samples among two cohorts of five trees collected from the same seed source and selected for genetic uniformity by genetic distance analysis. However, the direction of differential expression was not always consistent with that described in the original P. taeda SAGE data. Correlation was observed between gene expression and juvenile-mature wood anatomical characteristics by OPLS analysis. Four candidates (alpha-tubulin, porin MIP1, lipid transfer protein and aquaporin like protein) apparently had greatest influence on the wood traits measured. Speculative function of these genes in relation to juvenile-mature wood transition is briefly explored. Thus, we demonstrate the feasibility of exploiting SAGE data from a model species to identify consistently differentially expressed candidates in a related non-model species.

Bayesian analysis of an inverse Gaussian correlated frailty model

In survival analysis frailty is often used to model heterogeneity between individuals or correlation within clusters. Typically frailty is taken to be a continuous random effect, yielding a continuous mixture distribution for survival times. A Bayesian analysis of a correlated frailty model is discussed in the context of inverse Gaussian frailty. An MCMC approach is adopted and the deviance information criterion is used to compare models. As an illustration of the approach a bivariate data set of corneal graft survival times is analysed. (C) 2006 Elsevier B.V. All rights reserved.

The use of image analysis to estimate population growth rate in Daphnia magna

1. Population growth rate (PGR) is central to the theory of population ecology and is crucial for projecting population trends in conservation biology, pest management and wildlife harvesting. Furthermore, PGR is increasingly used to assess the effects of stressors. Image analysis that can automatically count and measure photographed individuals offers a potential methodology for estimating PGR. 2. This study evaluated two ways in which the PGR of Daphnia magna, exposed to different stressors, can be estimated using an image analysis system. The first method estimated PGR as the ratio of counts of individuals obtained at two different times, while the second method estimated PGR as the ratio of population sizes at two different times, where size is measured by the sum of the individuals' surface areas, i.e. total population surface area. This method is attractive if surface area is correlated with reproductive value (RV), as it is for D. magna, because of the theoretical result that PGR is the rate at which the population RV increases. 3. The image analysis system proved reliable and reproducible in counting populations of up to 440 individuals in 5 L of water. Image counts correlated well with manual counts but with a systematic underestimate of about 30%. This does not affect accuracy when estimating PGR as the ratio of two counts. Area estimates of PGR correlated well with count estimates, but were systematically higher, possibly reflecting their greater accuracy in the study situation. 4. Analysis of relevant scenarios suggested the correlation between RV and body size will generally be good for organisms in which fecundity correlates with body size. In these circumstances, area estimation of PGR is theoretically better than count estimation. 5. Synthesis and applications. There are both theoretical and practical advantages to area estimation of population growth rate when individuals' reproductive values are consistently well correlated with their surface areas. Because stressors may affect both the number and quality of individuals, area estimation of population growth rate should improve the accuracy of predicting stress impacts at the population level.

GeneMCL in microarray analysis

Accurately and reliably identifying the actual number of clusters present with a dataset of gene expression profiles, when no additional information on cluster structure is available, is a problem addressed by few algorithms. GeneMCL transforms microarray analysis data into a graph consisting of nodes connected by edges, where the nodes represent genes, and the edges represent the similarity in expression of those genes, as given by a proximity measurement. This measurement is taken to be the Pearson correlation coefficient combined with a local non-linear rescaling step. The resulting graph is input to the Markov Cluster (MCL) algorithm, which is an elegant, deterministic, non-specific and scalable method, which models stochastic flow through the graph. The algorithm is inherently affected by any cluster structure present, and rapidly decomposes a graph into cohesive clusters. The potential of the GeneMCL algorithm is demonstrated with a 5730 gene subset (IGS) of the Van't Veer breast cancer database, for which the clusterings are shown to reflect underlying biological mechanisms. (c) 2005 Elsevier Ltd. All rights reserved.

Analysis of molecular determinants of affinity and relative efficacy of a series of R- and S-2-(dipropylamino)tetralins at the 5-HT1A serotonin receptor

1 Factors influencing agonist affinity and relative efficacy have been studied for the 5-HT1A serotonin receptor using membranes of CHO cells expressing the human form of the receptor and a series of R-and S-2-(dipropylamino)tetralins (nonhydroxylated and monohydroxylated (5-OH, 6-OH, 7-OH, 8-OH) species). 2 Ligand binding studies were used to determine dissociation constants for agonist binding to the 5HT(1A) receptor: (a) K-i values for agonists were determined in competition versus the binding of the agonist [H-3]-8-OH DPAT. Competition data were all fitted best by a one-binding site model. (b) K-i values for agonists were also determined in competition versus the binding of the antagonist [H-3]-NAD-199. Competition data were all fitted best by a two-binding site model, and agonist affinities for the higher (K-h) and lower affinity (K-1) sites were determined. 3 The ability of the agonists to activate the 5-HT1A receptor was determined using stimulation of [S-35]-GTPgammaS binding. Maximal effects of agonists (E-max) and their potencies (EC50) were determined from concentration/response curves for stimulation of [S-35]-GTPgammaS binding. 4 K-1/K-h determined from ligand binding assays correlated with the relative efficacy (relative Em) of agonists determined in [S-35]-GTPgammaS binding assays. There was also a correlation between K-1/K-h and K-1/EC50 for agonists determined from ligand binding and [S-35]-GTPgammaS binding assays. 5 Simulations of agonist binding and effect data were performed using the Ternary Complex Model in order to assess the use of K-1/K-h for predicting the relative efficacy of agonists. British Journal of Pharmacology (2003) 138, 1129-1139. doi: 10. 1038/sj.bjp.705085.