Differential expression analysis with global network adjustment

<p>Background: Large-scale chromosomal deletions or other non-specific perturbations of the transcriptome can alter the expression of hundreds or thousands of genes, and it is of biological interest to understand which genes are most profoundly affected. We present a method for predicting a gene’s expression as a function of other genes thereby accounting for the effect of transcriptional regulation that confounds the identification of genes differentially expressed relative to a regulatory network. The challenge in constructing such models is that the number of possible regulator transcripts within a global network is on the order of thousands, and the number of biological samples is typically on the order of 10. Nevertheless, there are large gene expression databases that can be used to construct networks that could be helpful in modeling transcriptional regulation in smaller experiments.</p> <p>Results: We demonstrate a type of penalized regression model that can be estimated from large gene expression databases, and then applied to smaller experiments. The ridge parameter is selected by minimizing the cross-validation error of the predictions in the independent out-sample. This tends to increase the model stability and leads to a much greater degree of parameter shrinkage, but the resulting biased estimation is mitigated by a second round of regression. Nevertheless, the proposed computationally efficient “over-shrinkage” method outperforms previously used LASSO-based techniques. In two independent datasets, we find that the median proportion of explained variability in expression is approximately 25%, and this results in a substantial increase in the signal-to-noise ratio allowing more powerful inferences on differential gene expression leading to biologically intuitive findings. We also show that a large proportion of gene dependencies are conditional on the biological state, which would be impossible with standard differential expression methods.</p> <p>Conclusions: By adjusting for the effects of the global network on individual genes, both the sensitivity and reliability of differential expression measures are greatly improved.</p&gt

A semi-parametric approach to estimate risk functions associated with multi-dimensional exposure profiles: application to smoking and lung cancer

A common characteristic of environmental epidemiology is the multi-dimensional aspect of exposure patterns, frequently reduced to a cumulative exposure for simplicity of analysis. By adopting a flexible Bayesian clustering approach, we explore the risk function linking exposure history to disease. This approach is applied here to study the relationship between different smoking characteristics and lung cancer in the framework of a population based case control study

Invasive meningococcal disease in three siblings with hereditary deficiency of the 8th component of complement: Evidence for the importance of an early diagnosis

Deficiency of the eighth component of complement (C8) is a very rare primary immunodeficiency, associated with invasive, recurrent infections mainly caused by Neisseria species. We report functional and immunochemical C8 deficiency diagnosed in three Albanian siblings who presented with severe meningococcal infections at the age of 15 years, 4 years and 17 months, respectively. The youngest suffered serious complications (necrosis of fingers and toes requiring amputation). METHODS: Functional activity of the classical, alternative and mannose-binding lectin complement pathways was measured in serum from the 3 siblings and their parents (37-year-old woman and 42-year-old man). Forty healthy subjects (20 males and 20 females aged 4-38 years) served as normal controls. Serum complement factors were measured by haemolytic assays and immunoblotting. Sequence DNA analysis of the C8B gene was performed. RESULTS: Analyses of the three complement pathways revealed no haemolytic activity and also absence of C8beta in serum samples from all three siblings. The genetic analysis showed that the three siblings were homozygous for the p.Arg428* mutation in the C8B gene on chromosome 1p32 (MIM 120960). The parents were heterozygous for the mutation and presented normal complement activities. A 2-year follow-up revealed no further infective episodes in the siblings after antibiotic prophylaxis and meningococcal vaccination. CONCLUSIONS: Complement deficiencies are rare and their occurrence is often underestimated. In presence of invasive meningococcal infection, we highlight the importance of complement screening in patients and their relatives in order to discover any genetic defects which would render necessary prophylaxis to prevent recurrent infections and severe complications

Search For Heavy Pointlike Dirac Monopoles

We have searched for central production of a pair of photons with high transverse energies in $p\bar p$ collisions at $\sqrt{s} = 1.8$ TeV using $70 pb^{-1}$ of data collected with the D\O detector at the Fermilab Tevatron in 1994--1996. If they exist, virtual heavy pointlike Dirac monopoles could rescatter pairs of nearly real photons into this final state via a box diagram. We observe no excess of events above background, and set lower 95% C.L. limits of $610, 870, or 1580 GeV/c^2$ on the mass of a spin 0, 1/2, or 1 Dirac monopole.Comment: 12 pages, 4 figure

Limits on WWZ and WW\gamma couplings from p\bar{p}\to e\nu jj X events at \sqrt{s} = 1.8 TeV

We present limits on anomalous WWZ and WW-gamma couplings from a search for WW and WZ production in p-bar p collisions at sqrt(s)=1.8 TeV. We use p-bar p -> e-nu jjX events recorded with the D0 detector at the Fermilab Tevatron Collider during the 1992-1995 run. The data sample corresponds to an integrated luminosity of 96.0+-5.1 pb^(-1). Assuming identical WWZ and WW-gamma coupling parameters, the 95% CL limits on the CP-conserving couplings are -0.33<lambda<0.36 (Delta-kappa=0) and -0.43<Delta-kappa<0.59 (lambda=0), for a form factor scale Lambda = 2.0 TeV. Limits based on other assumptions are also presented.Comment: 11 pages, 2 figures, 2 table

Towards the clinical implementation of pharmacogenetics in bipolar disorder.

BackgroundBipolar disorder (BD) is a psychiatric illness defined by pathological alterations between the mood states of mania and depression, causing disability, imposing healthcare costs and elevating the risk of suicide. Although effective treatments for BD exist, variability in outcomes leads to a large number of treatment failures, typically followed by a trial and error process of medication switches that can take years. Pharmacogenetic testing (PGT), by tailoring drug choice to an individual, may personalize and expedite treatment so as to identify more rapidly medications well suited to individual BD patients.DiscussionA number of associations have been made in BD between medication response phenotypes and specific genetic markers. However, to date clinical adoption of PGT has been limited, often citing questions that must be answered before it can be widely utilized. These include: What are the requirements of supporting evidence? How large is a clinically relevant effect? What degree of specificity and sensitivity are required? Does a given marker influence decision making and have clinical utility? In many cases, the answers to these questions remain unknown, and ultimately, the question of whether PGT is valid and useful must be determined empirically. Towards this aim, we have reviewed the literature and selected drug-genotype associations with the strongest evidence for utility in BD.SummaryBased upon these findings, we propose a preliminary panel for use in PGT, and a method by which the results of a PGT panel can be integrated for clinical interpretation. Finally, we argue that based on the sufficiency of accumulated evidence, PGT implementation studies are now warranted. We propose and discuss the design for a randomized clinical trial to test the use of PGT in the treatment of BD

The Dijet Mass Spectrum and a Search for Quark Compositeness in bar{p}p Collisions at sqrt{s} = 1.8 TeV

Using the DZero detector at the 1.8 TeV pbarp Fermilab Tevatron collider, we have measured the inclusive dijet mass spectrum in the central pseudorapidity region |eta_jet| < 1.0 for dijet masses greater than 200 Gev/c^2. We have also measured the ratio of spectra sigma(|eta_jet| < 0.5)/sigma(0.5 < |eta_jet| < 1.0). The order alpha_s^3 QCD predictions are in good agreement with the data and we rule out models of quark compositeness with a contact interaction scale < 2.4 TeV at the 95% confidence level.Comment: 11 pages, 4 figures, 2 tables, submitted to Phys. Rev. Let

Zgamma Production in pbarp Collisions at sqrt(s)=1.8 TeV and Limits on Anomalous ZZgamma and Zgammagamma Couplings

We present a study of Z +gamma + X production in p-bar p collisions at sqrt{S}=1.8 TeV from 97 (87) pb^{-1} of data collected in the eegamma (mumugamma) decay channel with the D0 detector at Fermilab. The event yield and kinematic characteristics are consistent with the Standard Model predictions. We obtain limits on anomalous ZZgamma and Zgammagamma couplings for form factor scales Lambda = 500 GeV and Lambda = 750 GeV. Combining this analysis with our previous results yields 95% CL limits |h{Z}_{30}| < 0.36, |h{Z}_{40}| < 0.05, |h{gamma}_{30}| < 0.37, and |h{gamma}_{40}| < 0.05 for a form factor scale Lambda=750 GeV.Comment: 17 Pages including 2 Figures. Submitted to PR

Identification of candidate tumour suppressor genes frequently methylated in renal cell carcinoma

Promoter region hyermethylation and transcriptional silencing is a frequent cause of tumour suppressor gene (TSG) inactivation in many types of human cancers. Functional epigenetic studies, in which gene expression is induced by treatment with demethylating agents, may identify novel genes with tumour-specific methylation. We used high-density gene expression microarrays in a functional epigenetic study of 11 renal cell carcinoma (RCC) cell lines. Twenty-eight genes were then selected for analysis of promoter methylation status in cell lines and primary RCC. Eight genes (BNC1, PDLIM4, RPRM, CST6, SFRP1, GREM1, COL14A1 and COL15A1) showed frequent (30% of RCC tested) tumour-specific promoter region methylation. Hypermethylation was associated with transcriptional silencing. Re-expression of BNC1, CST6, RPRM and SFRP1 suppressed the growth of RCC cell lines and RNA interference knock-down of BNC1, SFRP1 and COL14A1 increased the growth of RCC cell lines. Methylation of BNC1 or COL14A1 was associated with a poorer prognosis independent of tumour size, stage or grade. The identification of these epigenetically inactivated candidate RCC TSGs can provide insights into renal tumourigenesis and a basis for developing novel therapies and biomarkers for prognosis and detection. © 2010 Macmillan Publishers Limited.Published versio

A Measurement of the W Boson Mass

We report a measurement of the W boson mass based on an integrated luminosity of 82 pb$^{-1}$ from \ppbar collisions at $\sqrt{s}=1.8$ TeV recorded in 1994--1995 by the \Dzero detector at the Fermilab Tevatron. We identify W bosons by their decays to $e\nu$ and extract the mass by fitting the transverse mass spectrum from 28,323 W boson candidates. A sample of 3,563 dielectron events, mostly due to Z to ee decays, constrains models of W boson production and the detector. We measure \mw=80.44\pm0.10(stat)\pm0.07(syst)~GeV. By combining this measurement with our result from the 1992--1993 data set, we obtain \mw=80.43\pm0.11 GeV.Comment: 11 pages, 5 figure