Module networks revisited: computational assessment and prioritization of model predictions

The solution of high-dimensional inference and prediction problems in computational biology is almost always a compromise between mathematical theory and practical constraints such as limited computational resources. As time progresses, computational power increases but well-established inference methods often remain locked in their initial suboptimal solution. We revisit the approach of Segal et al. (2003) to infer regulatory modules and their condition-specific regulators from gene expression data. In contrast to their direct optimization-based solution we use a more representative centroid-like solution extracted from an ensemble of possible statistical models to explain the data. The ensemble method automatically selects a subset of most informative genes and builds a quantitatively better model for them. Genes which cluster together in the majority of models produce functionally more coherent modules. Regulators which are consistently assigned to a module are more often supported by literature, but a single model always contains many regulator assignments not supported by the ensemble. Reliably detecting condition-specific or combinatorial regulation is particularly hard in a single optimum but can be achieved using ensemble averaging.Comment: 8 pages REVTeX, 6 figure

Changes in undergraduate student alcohol consumption as they progress through university

BACKGROUND: Unhealthy alcohol use amongst university students is a major public health concern. Although previous studies suggest a raised level of consumption amongst the UK student population there is little consistent information available about the pattern of alcohol consumption as they progress through university. The aim of the current research was to describe drinking patterns of UK full-time undergraduate students as they progress through their degree course. METHOD: Data were collected over three years from 5895 undergraduate students who began their studies in either 2000 or 2001. Longitudinal data (i.e. Years 1–3) were available from 225 students. The remaining 5670 students all responded to at least one of the three surveys (Year 1 n = 2843; Year 2 n = 2219; Year 3 n = 1805). Results: Students reported consuming significantly more units of alcohol per week at Year 1 than at Years 2 or 3 of their degree. Male students reported a higher consumption of units of alcohol than their female peers. When alcohol intake was classified using the Royal College of Physicians guidelines [1] there was no difference between male and females students in terms of the percentage exceeding recommended limits. Compared to those who were low level consumers students who reported drinking above low levels at Year 1 had at least 10 times the odds of continuing to consume above low levels at year 3. Students who reported higher levels of drinking were more likely to report that alcohol had a negative impact on their studies, finances and physical health. Consistent with the reduction in units over time students reported lower levels of negative impact during Year 3 when compared to Year 1. CONCLUSION: The current findings suggest that student alcohol consumption declines over their undergraduate studies; however weekly levels of consumption at Year 3 remain high for a substantial number of students. The persistence of high levels of consumption in a large population of students suggests the need for effective preventative and treatment interventions for all year groups

Simple correction improving long-term reproducibility of HPLC-MS

Chromatographic peak areas in long series of HPLC-MS experiments often vary, which decreases reproducibility and may cause bias in the results. It was found that the sensitivitiy of various components change differently; in our case variability is in the order of 20-40%; and it is most likely due to changing conditions in ESI ionization. The most often used peak area correction methods do not take this effect into account. The change in peak areas can be well described by a polynomial function; we found that a 4th order polynomial is most often suitable. We suggest a simple correction algorithm based on polynomial fitting. When the experiments were inherently well reproducible, this correction improved reproducibility from 12% to 3% (on average for various components). When random errors were larger, this improvement was less significant (15% to 12% in nano-ESI), but nevertheless essential in order to avoid possible bias in the results

Perinatal risk factors for neonatal encephalopathy: an unmatched case-control study

OBJECTIVE: Neonatal encephalopathy (NE) is the third leading cause of child mortality. Preclinical studies suggest infection and inflammation can sensitise or precondition the newborn brain to injury. This study examined perinatal risks factor for NE in Uganda. DESIGN: Unmatched case-control study. SETTING: Mulago National Referral Hospital, Kampala, Uganda. METHODS: 210 term infants with NE and 409 unaffected term infants as controls were recruited over 13 months. Data were collected on preconception, antepartum and intrapartum exposures. Blood culture, species-specific bacterial real-time PCR, C reactive protein and placental histology for chorioamnionitis and funisitis identified maternal and early newborn infection and inflammation. Multivariable logistic regression examined associations with NE. RESULTS: Neonatal bacteraemia (adjusted OR (aOR) 8.67 (95% CI 1.51 to 49.74), n=315) and histological funisitis (aOR 11.80 (95% CI 2.19 to 63.45), n=162) but not chorioamnionitis (aOR 3.20 (95% CI 0.66 to 15.52), n=162) were independent risk factors for NE. Among encephalopathic infants, neonatal case fatality was not significantly higher when exposed to early neonatal bacteraemia (OR 1.65 (95% CI 0.62 to 4.39), n=208). Intrapartum antibiotic use did not improve neonatal survival (p=0.826). After regression analysis, other identified perinatal risk factors (n=619) included hypertension in pregnancy (aOR 3.77), male infant (aOR 2.51), non-cephalic presentation (aOR 5.74), lack of fetal monitoring (aOR 2.75), augmentation (aOR 2.23), obstructed labour (aOR 3.8) and an acute intrapartum event (aOR 8.74). CONCLUSIONS: Perinatal infection and inflammation are independent risk factors for NE in this low-resource setting, supporting a role in the aetiological pathway of term brain injury. Intrapartum antibiotic administration did not mitigate against adverse outcomes. The importance of intrapartum risk factors in this sub-Saharan African setting is highlighted

RNA polymerase II stalling promotes nucleosome occlusion and pTEFb recruitment to drive immortalization by Epstein-Barr virus

Epstein-Barr virus (EBV) immortalizes resting B-cells and is a key etiologic agent in the development of numerous cancers. The essential EBV-encoded protein EBNA 2 activates the viral C promoter (Cp) producing a message of ~120 kb that is differentially spliced to encode all EBNAs required for immortalization. We have previously shown that EBNA 2-activated transcription is dependent on the activity of the RNA polymerase II (pol II) C-terminal domain (CTD) kinase pTEFb (CDK9/cyclin T1). We now demonstrate that Cp, in contrast to two shorter EBNA 2-activated viral genes (LMP 1 and 2A), displays high levels of promoter-proximally stalled pol II despite being constitutively active. Consistent with pol II stalling, we detect considerable pausing complex (NELF/DSIF) association with Cp. Significantly, we observe substantial Cp-specific pTEFb recruitment that stimulates high-level pol II CTD serine 2 phosphorylation at distal regions (up to +75 kb), promoting elongation. We reveal that Cp-specific pol II accumulation is directed by DNA sequences unfavourable for nucleosome assembly that increase TBP access and pol II recruitment. Stalled pol II then maintains Cp nucleosome depletion. Our data indicate that pTEFb is recruited to Cp by the bromodomain protein Brd4, with polymerase stalling facilitating stable association of pTEFb. The Brd4 inhibitor JQ1 and the pTEFb inhibitors DRB and Flavopiridol significantly reduce Cp, but not LMP1 transcript production indicating that Brd4 and pTEFb are required for Cp transcription. Taken together our data indicate that pol II stalling at Cp promotes transcription of essential immortalizing genes during EBV infection by (i) preventing promoter-proximal nucleosome assembly and ii) necessitating the recruitment of pTEFb thereby maintaining serine 2 CTD phosphorylation at distal regions

Barriers to Scientific Contributions: The Author’s Formula

Recently I completed a review of the empirical research on scientific journals (Armstrong 1982). This review provided evidence for an “author’s formula,” a set of rules that authors can use to increase the likelihood and speed of acceptance of their manuscripts. Authors should: (1) not pick an important problem, (2) not challenge existing beliefs, (3) not obtain surprising results, (4) not use simple methods, (5) not provide full disclosure, and (6) not write clearly. Peters & Ceci (P&C) are obviously ignorant of the author’s formula. In their extension of the Kosinski study (Ross 1979; 1980), they broke most of the rules

The WEBT Campaign on the Blazar 3C279 in 2006

The quasar 3C279 was the target of an extensive multiwavelength monitoring campaign from January through April 2006, including an optical-IR-radio monitoring campaign by the Whole Earth Blazar Telescope (WEBT) collaboration. In this paper we focus on the results of the WEBT campaign. The source exhibited substantial variability of optical flux and spectral shape, with a characteristic time scale of a few days. The variability patterns throughout the optical BVRI bands were very closely correlated with each other. In intriguing contrast to other (in particular, BL Lac type) blazars, we find a lag of shorter- behind longer-wavelength variability throughout the RVB ranges, with a time delay increasing with increasing frequency. Spectral hardening during flares appears delayed with respect to a rising optical flux. This, in combination with the very steep IR-optical continuum spectral index of ~ 1.5 - 2.0, may indicate a highly oblique magnetic field configuration near the base of the jet. An alternative explanation through a slow (time scale of several days) acceleration mechanism would require an unusually low magnetic field of < 0.2 G, about an order of magnitude lower than inferred from previous analyses of simultaneous SEDs of 3C279 and other FSRQs with similar properties.Comment: Accepted for publication in Ap

Bayesian Centroid Estimation for Motif Discovery

Biological sequences may contain patterns that are signal important biomolecular functions; a classical example is regulation of gene expression by transcription factors that bind to specific patterns in genomic promoter regions. In motif discovery we are given a set of sequences that share a common motif and aim to identify not only the motif composition, but also the binding sites in each sequence of the set. We present a Bayesian model that is an extended version of the model adopted by the Gibbs motif sampler, and propose a new centroid estimator that arises from a refined and meaningful loss function for binding site inference. We discuss the main advantages of centroid estimation for motif discovery, including computational convenience, and how its principled derivation offers further insights about the posterior distribution of binding site configurations. We also illustrate, using simulated and real datasets, that the centroid estimator can differ from the maximum a posteriori estimator.Comment: 24 pages, 9 figure

Bayesian integration of isotope ratio for geographic sourcing of castor beans

pre-printRecent years have seen an increase in the forensic interest associated with the poison ricin, which is extracted from the seeds of the Ricinus communis plant. Both light element (C, N, O, and H) and strontium (Sr) isotope ratios have previously been used to associate organic material with geographic regions of origin. We present a Bayesian integration methodology that can more accurately predict the region of origin for a castor bean than individual models developed independently for light element stable isotopes or Sr isotope ratios. Our results demonstrate a clear improvement in the ability to correctly classify regions based on the integrated model with a class accuracy of 60.9 ± 2.1% versus 55.9 ± 2.1% and 40.2 ± 1.8% for the light element and strontium (Sr) isotope ratios, respectively. In addition, we show graphically the strengths and weaknesses of each dataset in respect to class prediction and how the integration of these datasets strengthens the overall model

Early cranial ultrasound findings among infants with neonatal encephalopathy in Uganda: an observational study.

BACKGROUND: In sub-Saharan Africa, the timing and nature of brain injury and their relation to mortality in neonatal encephalopathy (NE) are unknown. We evaluated cranial ultrasound (cUS) scans from term Ugandan infants with and without NE for evidence of brain injury. METHODS: Infants were recruited from a national referral hospital in Kampala. Cases (184) had NE and controls (100) were systematically selected unaffected term infants. All had cUS scans <36 h reported blind to NE status. RESULTS: Scans were performed at median age 11.5 (interquartile range (IQR): 5.2-20.2) and 8.4 (IQR: 3.6-13.5) hours, in cases and controls respectively. None had established antepartum injury. Major evolving injury was reported in 21.2% of the cases vs. 1.0% controls (P < 0.001). White matter injury was not significantly associated with bacteremia in encephalopathic infants (odds ratios (OR): 3.06 (95% confidence interval (CI): 0.98-9.60). Major cUS abnormality significantly increased the risk of neonatal death (case fatality 53.9% with brain injury vs. 25.9% without; OR: 3.34 (95% CI: 1.61-6.95)). CONCLUSION: In this low-resource setting, there was no evidence of established antepartum insult, but a high proportion of encephalopathic infants had evidence of major recent and evolving brain injury on early cUS imaging, suggesting prolonged or severe acute exposure to hypoxia-ischemia (HI). Early abnormalities were a significant predictor of death