Modeling SAGE tag formation and its effects on data interpretation within a Bayesian framework

<p>Abstract</p> <p>Background</p> <p>Serial Analysis of Gene Expression (SAGE) is a high-throughput method for inferring mRNA expression levels from the experimentally generated sequence based tags. Standard analyses of SAGE data, however, ignore the fact that the probability of generating an observable tag varies across genes and between experiments. As a consequence, these analyses result in biased estimators and posterior probability intervals for gene expression levels in the transcriptome.</p> <p>Results</p> <p>Using the yeast <it>Saccharomyces cerevisiae </it>as an example, we introduce a new Bayesian method of data analysis which is based on a model of SAGE tag formation. Our approach incorporates the variation in the probability of tag formation into the interpretation of SAGE data and allows us to derive exact joint and approximate marginal posterior distributions for the mRNA frequency of genes detectable using SAGE. Our analysis of these distributions indicates that the frequency of a gene in the tag pool is influenced by its mRNA frequency, the cleavage efficiency of the anchoring enzyme (AE), and the number of informative and uninformative AE cleavage sites within its mRNA.</p> <p>Conclusion</p> <p>With a mechanistic, model based approach for SAGE data analysis, we find that inter-genic variation in SAGE tag formation is large. However, this variation can be estimated and, importantly, accounted for using the methods we develop here. As a result, SAGE based estimates of mRNA frequencies can be adjusted to remove the bias introduced by the SAGE tag formation process.</p

Production of phi mesons at mid-rapidity in sqrt(s_NN) = 200 GeV Au+Au collisions at RHIC

We present the first results of meson production in the K^+K^- decay channel from Au+Au collisions at sqrt(s_NN) = 200 GeV as measured at mid-rapidity by the PHENIX detector at RHIC. Precision resonance centroid and width values are extracted as a function of collision centrality. No significant variation from the PDG accepted values is observed. The transverse mass spectra are fitted with a linear exponential function for which the derived inverse slope parameter is seen to be constant as a function of centrality. These data are also fitted by a hydrodynamic model with the result that the freeze-out temperature and the expansion velocity values are consistent with the values previously derived from fitting single hadron inclusive data. As a function of transverse momentum the collisions scaled peripheral.to.central yield ratio RCP for the is comparable to that of pions rather than that of protons. This result lends support to theoretical models which distinguish between baryons and mesons instead of particle mass for explaining the anomalous proton yield.Comment: 326 authors, 24 pages text, 23 figures, 6 tables, RevTeX 4. To be submitted to Physical Review C as a regular article. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Genetic dissection of an amygdala microcircuit that gates conditioned fear

The role of different amygdala nuclei (neuroanatomical subdivisions) in processing Pavlovian conditioned fear has been studied extensively, but the function of the heterogeneous neuronal subtypes within these nuclei remains poorly understood. Here we use molecular genetic approaches to map the functional connectivity of a subpopulation of GABA-containing neurons, located in the lateral subdivision of the central amygdala (CEl), which express protein kinase C-δ (PKC-δ). Channelrhodopsin-2-assisted circuit mapping in amygdala slices and cell-specific viral tracing indicate that PKC-δ^+ neurons inhibit output neurons in the medial central amygdala (CEm), and also make reciprocal inhibitory synapses with PKC-δ^− neurons in CEl. Electrical silencing of PKC-δ^+ neurons in vivo suggests that they correspond to physiologically identified units that are inhibited by the conditioned stimulus, called Cel_(off) units. This correspondence, together with behavioural data, defines an inhibitory microcircuit in CEl that gates CEm output to control the level of conditioned freezing

Search for pair-produced long-lived neutral particles decaying to jets in the ATLAS hadronic calorimeter in ppcollisions at √s=8TeV

The ATLAS detector at the Large Hadron Collider at CERN is used to search for the decay of a scalar boson to a pair of long-lived particles, neutral under the Standard Model gauge group, in 20.3fb−1of data collected in proton–proton collisions at √s=8TeV. This search is sensitive to long-lived particles that decay to Standard Model particles producing jets at the outer edge of the ATLAS electromagnetic calorimeter or inside the hadronic calorimeter. No significant excess of events is observed. Limits are reported on the product of the scalar boson production cross section times branching ratio into long-lived neutral particles as a function of the proper lifetime of the particles. Limits are reported for boson masses from 100 GeVto 900 GeV, and a long-lived neutral particle mass from 10 GeVto 150 GeV

Measurement of the cross-section of high transverse momentum vector bosons reconstructed as single jets and studies of jet substructure in pp collisions at √s = 7 TeV with the ATLAS detector

This paper presents a measurement of the cross-section for high transverse momentum W and Z bosons produced in pp collisions and decaying to all-hadronic final states. The data used in the analysis were recorded by the ATLAS detector at the CERN Large Hadron Collider at a centre-of-mass energy of √s = 7 TeV;{\rm Te}{\rm V}$and correspond to an integrated luminosity of$4.6\;{\rm f}{{{\rm b}}^{-1}}$. The measurement is performed by reconstructing the boosted W or Z bosons in single jets. The reconstructed jet mass is used to identify the W and Z bosons, and a jet substructure method based on energy cluster information in the jet centre-of-mass frame is used to suppress the large multi-jet background. The cross-section for events with a hadronically decaying W or Z boson, with transverse momentum${{p}_{{\rm T}}}\gt 320\;{\rm Ge}{\rm V}$and pseudorapidity$|\eta |\lt 1.9$, is measured to be${{\sigma }_{W+Z}}=8.5\pm 1.7$ pb and is compared to next-to-leading-order calculations. The selected events are further used to study jet grooming techniques

Search for direct pair production of the top squark in all-hadronic final states in proton-proton collisions at s√=8 TeV with the ATLAS detector

The results of a search for direct pair production of the scalar partner to the top quark using an integrated luminosity of 20.1fb−1 of proton–proton collision data at √s = 8 TeV recorded with the ATLAS detector at the LHC are reported. The top squark is assumed to decay via t˜→tχ˜01 or t˜→ bχ˜±1 →bW(∗)χ˜01 , where χ˜01 (χ˜±1 ) denotes the lightest neutralino (chargino) in supersymmetric models. The search targets a fully-hadronic final state in events with four or more jets and large missing transverse momentum. No significant excess over the Standard Model background prediction is observed, and exclusion limits are reported in terms of the top squark and neutralino masses and as a function of the branching fraction of t˜ → tχ˜01 . For a branching fraction of 100%, top squark masses in the range 270–645 GeV are excluded for χ˜01 masses below 30 GeV. For a branching fraction of 50% to either t˜ → tχ˜01 or t˜ → bχ˜±1 , and assuming the χ˜±1 mass to be twice the χ˜01 mass, top squark masses in the range 250–550 GeV are excluded for χ˜01 masses below 60 GeV

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

SUMOylation of the Forkhead Transcription Factor FOXL2 Promotes Its Stabilization/Activation through Transient Recruitment to PML Bodies

International audienceBACKGROUND: FOXL2 is a transcription factor essential for ovarian development and maintenance. It is mutated in the genetic condition called Blepharophimosis Ptosis Epicantus inversus Syndrome (BPES) and in cases of isolated premature ovarian failure. We and others have previously shown that FOXL2 undergoes several post-translational modifications. METHODS AND PRINCIPAL FINDINGS: Here, using cells in culture, we show that interference with FOXL2 SUMOylation leads to a robust inhibition of its transactivation ability, which correlates with a decreased stability. Interestingly, FOXL2 SUMOylation promotes its transient recruitment to subnuclear structures that we demonstrate to be PML (Promyelocytic Leukemia) Nuclear Bodies. Since PML bodies are known to be sites where post-translational modifications of nuclear factors take place, we used tandem mass spectrometry to identify new post-translational modifications of FOXL2. Specifically, we detected four phosphorylated, one sulfated and three acetylated sites. CONCLUSIONS: By analogy with other transcription factors, we propose that PML Nuclear Bodies might transiently recruit FOXL2 to the vicinity of locally concentrated enzymes that could be involved in the post-translational maturation of FOXL2. FOXL2 acetylation, sulfation, phosphorylation as well as other modifications yet to be discovered might alter the transactivation capacity of FOXL2 and/or its stability, thus modulating its global intracellular activity

U-Shaped Relation between Plasma Oxytocin Levels and Behavior in the Trust Game

10.1371/journal.pone.0051095PLoS ONE712

Common genetic variation in IGF1, IGFBP-1, and IGFBP-3 in relation to mammographic density: a cross-sectional study

INTRODUCTION: Mammographic density is one of the strongest risk factors for breast cancer and is believed to represent epithelial and stromal proliferation. Because of the high heritability of breast density, and the role of the insulin-like growth factor (IGF) pathway in cellular proliferation and breast development, we examined the association between common genetic variation in this pathway and mammographic density. METHODS: We conducted a cross-sectional analysis among controls (n = 1,121) who were between the ages of 42 and 78 years at mammography, from a breast cancer case-control study nested within the Nurses' Health Study cohort. At the time of mammography, 204 women were premenopausal and 917 were postmenopausal. We genotyped 29 haplotype-tagging SNPs demonstrated to capture common genetic variation in IGF1, IGF binding protein (IGFBP)-1, and IGFBP-3. RESULTS: Common haplotype patterns in three of the four haplotype blocks spanning the gene encoding IGF1 were associated with mammographic density. Haplotype patterns in block 1 (p = 0.03), block 3 (p = 0.009), and block 4 (p = 0.007) were associated with mammographic density, whereas those in block 2 were not. None of the common haplotypes in the three haplotype blocks spanning the genes encoding IGFBP-1/IGFBP-3 were significantly associated with mammographic density. Two haplotype-tagging SNPs in IGF1, rs1520220 and rs2946834, showed a strong association with mammographic density. Those with the homozygous variant genotype for rs1520220 had a mean percentage mammographic density of 19.6% compared with those with the homozygous wild-type genotype, who had a mean percentage mammographic density of 27.9% (p for trend < 0.0001). Those that were homozygous variant for rs2946834 had a mean percentage mammographic density of 23.2% compared with those who were homozygous wild-type with a mean percentage mammographic density of 28.2% (p for trend = 0.0004). Permutation testing demonstrated that results as strong as these are unlikely to occur by chance (p = 0.0005). CONCLUSION: Common genetic variation in IGF1 is strongly associated with percentage mammographic density