RNA-seq: technical variability and sampling

<p>Abstract</p> <p>Background</p> <p>RNA-seq is revolutionizing the way we study transcriptomes. mRNA can be surveyed without prior knowledge of gene transcripts. Alternative splicing of transcript isoforms and the identification of previously unknown exons are being reported. Initial reports of differences in exon usage, and splicing between samples as well as quantitative differences among samples are beginning to surface. Biological variation has been reported to be larger than technical variation. In addition, technical variation has been reported to be in line with expectations due to random sampling. However, strategies for dealing with technical variation will differ depending on the magnitude. The size of technical variance, and the role of sampling are examined in this manuscript.</p> <p>Results</p> <p>In this study three independent Solexa/Illumina experiments containing technical replicates are analyzed. When coverage is low, large disagreements between technical replicates are apparent. Exon detection between technical replicates is highly variable when the coverage is less than 5 reads per nucleotide and estimates of gene expression are more likely to disagree when coverage is low. Although large disagreements in the estimates of expression are observed at all levels of coverage.</p> <p>Conclusions</p> <p>Technical variability is too high to ignore. Technical variability results in inconsistent detection of exons at low levels of coverage. Further, the estimate of the relative abundance of a transcript can substantially disagree, even when coverage levels are high. This may be due to the low sampling fraction and if so, it will persist as an issue needing to be addressed in experimental design even as the next wave of technology produces larger numbers of reads. We provide practical recommendations for dealing with the technical variability, without dramatic cost increases.</p

Uncovering hidden genetic variation in photosynthesis of field‐grown maize under ozone pollution

Ozone is the most damaging air pollutant to crops, currently reducing Midwest US maize production by up to 10%, yet there has been very little effort to adapt germ‐ plasm for ozone tolerance. Ozone enters plants through stomata, reacts to form reactive oxygen species in the apoplast and ultimately decreases photosynthetic C gain. In this study, 10 diverse inbred parents were crossed in a half‐diallel design to create 45 F1 hybrids, which were tested for ozone response in the field using free air concentration enrichment (FACE). Ozone stress increased the heritability of pho‐ tosynthetic traits and altered genetic correlations among traits. Hybrids from par‐ ents Hp301 and NC338 showed greater sensitivity to ozone stress, and disrupted relationships among photosynthetic traits. The physiological responses underlying sensitivity to ozone differed in hybrids from the two parents, suggesting multiple mechanisms of response to oxidative stress. FACE technology was essential to this evaluation because genetic variation in photosynthesis under elevated ozone was not predictable based on performance at ambient ozone. These findings suggest that selection under elevated ozone is needed to identify deleterious alleles in the world's largest commodity crop

Changes in reflectin protein phosphorylation are associated with dynamic iridescence in squid

Author Posting. © The Author(s), 2009. This is the author's version of the work. It is posted here by permission of The Royal Society for personal use, not for redistribution. The definitive version was published in Journal of The Royal Society Interface 6 (2010): 549-560, doi:10.1098/rsif.2009.0299.Many cephalopods exhibit remarkable dermal iridescence, a component of their complex, dynamic camouflage and communication. In the species Euprymna scolopes, the light-organ iridescence is static and is due to reflectin protein-based platelets assembled into lamellar thin-film reflectors called iridosomes, contained within iridescent cells called iridocytes. Squid in the family Loliginidae appear to be unique in that the dermis possesses a dynamic iridescent component, with reflective, colored structures that are assembled and disassembled under the control of the muscarinic cholinergic system and the associated neurotransmitter acetylcholine (Mathger et al. 2004). Here we present the sequences and characterization of three new members of the reflectin family associated with the dynamically changeable iridescence in Loligo and not found in static Euprymna iridophores. In addition, we show that application of genistein, a protein tyrosine kinase inhibitor, suppresses acetylcholine- and calcium-induced iridescence in Loligo. We further demonstrate that two of these novel reflectins are extensively phosphorylated in concert with the activation of iridescence by exogenous acetylcholine. This phosphorylation and the correlated iridescence can be blocked with genistein. Our results suggest that tyrosine phosphorylation of reflectin proteins is involved in the regulation of dynamic iridescence in Loligo.We gratefully acknowledge support from Anteon contract F33615-03-D-5408 to the Marine Biological Laboratory, Woods Hole, MA and grant # W911NF-06-1-0285 from the Army Research Office to D.E.M

Financial considerations in the conduct of multi-centre randomised controlled trials: evidence from a qualitative study.

National Coordinating Centre for Research Methodology; Medical Research Council, UK Department of Health; Chief Scientist OfficeNot peer reviewedPublisher PD

Constraints on the χ_(c1) versus χ_(c2) polarizations in proton-proton collisions at √s = 8 TeV

The polarizations of promptly produced χ_(c1) and χ_(c2) mesons are studied using data collected by the CMS experiment at the LHC, in proton-proton collisions at √s=8  TeV. The χ_c states are reconstructed via their radiative decays χ_c → J/ψγ, with the photons being measured through conversions to e⁺e⁻, which allows the two states to be well resolved. The polarizations are measured in the helicity frame, through the analysis of the χ_(c2) to χ_(c1) yield ratio as a function of the polar or azimuthal angle of the positive muon emitted in the J/ψ → μ⁺μ⁻ decay, in three bins of J/ψ transverse momentum. While no differences are seen between the two states in terms of azimuthal decay angle distributions, they are observed to have significantly different polar anisotropies. The measurement favors a scenario where at least one of the two states is strongly polarized along the helicity quantization axis, in agreement with nonrelativistic quantum chromodynamics predictions. This is the first measurement of significantly polarized quarkonia produced at high transverse momentum

Effects of antiplatelet therapy on stroke risk by brain imaging features of intracerebral haemorrhage and cerebral small vessel diseases: subgroup analyses of the RESTART randomised, open-label trial

Background Findings from the RESTART trial suggest that starting antiplatelet therapy might reduce the risk of recurrent symptomatic intracerebral haemorrhage compared with avoiding antiplatelet therapy. Brain imaging features of intracerebral haemorrhage and cerebral small vessel diseases (such as cerebral microbleeds) are associated with greater risks of recurrent intracerebral haemorrhage. We did subgroup analyses of the RESTART trial to explore whether these brain imaging features modify the effects of antiplatelet therapy

Nurses' perceived barriers to the implementation of a Fall Prevention Clinical Practice Guideline in Singapore hospitals

BackgroundTheories of behavior change indicate that an analysis of barriers to change is helpful when trying to influence professional practice. The aim of this study was to assess the perceived barriers to practice change by eliciting nurses\u27 opinions with regard to barriers to, and facilitators of, implementation of a Fall Prevention clinical practice guideline in five acute care hospitals in Singapore.MethodsNurses were surveyed to identify their perceptions regarding barriers to implementation of clinical practice guidelines in their practice setting. The validated questionnaire, \u27Barriers and facilitators assessment instrument\u27, was administered to nurses (n = 1830) working in the medical, surgical, geriatric units, at five acute care hospitals in Singapore.ResultsAn 80.2% response rate was achieved. The greatest barriers to implementation of clinical practice guidelines reported included: knowledge and motivation, availability of support staff, access to facilities, health status of patients, and, education of staff and patients.ConclusionNumerous barriers to the use of the Fall Prevention Clinical Practice Guideline have been identified. This study has laid the foundation for further research into implementation of clinical practice guidelines in Singapore by identifying barriers to change in acute care settings.<br /

Recruitment of regulatory T cells is correlated with hypoxia-induced CXCR4 expression, and is associated with poor prognosis in basal-like breast cancers

Introduction: Basal-like breast cancers behave more aggressively despite the presence of a dense lymphoid infiltrate. We hypothesised that immune suppression in this subtype may be due to T regulatory cells (Treg) recruitment driven by hypoxia-induced up-regulation of CXCR4 in Treg.Methods: Immunoperoxidase staining for FOXP3 and CXCL12 was performed on tissue microarrays from 491 breast cancers. The hypoxia-associated marker carbonic anhydrase IX (CA9) and double FOXP3/CXCR4 staining were performed on sections from a subset of these cancers including 10 basal-like and 11 luminal cancers matched for tumour grade.Results: High Treg infiltration correlated with tumour CXCL12 positivity (OR 1.89, 95% CI 1.22 to 2.94, P = 0.004) and basal phenotype (OR 3.14, 95% CI 1.08 to 9.17, P = 0.004) in univariate and multivariate analyses. CXCL12 positivity correlated with improved survival (P = 0.005), whereas high Treg correlated with shorter survival for all breast cancers (P = 0.001), luminal cancers (P &lt; 0.001) and basal-like cancers (P = 0.040) that were confirmed in a multivariate analysis (OR 1.61, 95% CI 1.02 to 2.53, P = 0.042). In patients treated with hormone therapy, high Treg were associated with a shorter survival in a multivariate analysis (OR 1.78, 95% CI 1.01 to 3.15, P = 0.040). There was a tendency for luminal cancers to show CXCL12 expression (102/138, 74%) compared to basal-like cancers (16/27, 59%), which verged on statistical significance (P = 0.050). Up-regulation of CXCR4 in Treg correlated with the basal-like phenotype (P = 0.029) and tumour hypoxia, as indicated by CA9 expression (P = 0.049).Conclusions: Our data show that in the setting of hypoxia and CXCR4 up-regulation in Treg, CXCL12 expression may have the negative consequence of enhancing Treg recruitment and suppressing the anti-tumour immune response. © 2011 Yan et al.; licensee BioMed Central Ltd

Evolution of Genome Size and Complexity in Pinus

BACKGROUND: Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood. METHODOLOGY/PRINCIPAL FINDINGS: Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA. CONCLUSIONS/SIGNIFICANCE: Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes

Measurement of the Y(1S) pair production cross section and search for resonances decaying to Y(1S)μ⁺μ⁻ in proton-proton collisions at √s = 13 TeV

The fiducial cross section for Y(1S) pair production in proton-proton collisions at a center-of-mass energy of 13 TeV in the region where both Y(1S) mesons have an absolute rapidity below 2.0 is measured to be 79±11(stat)±6(syst)±3(B) pb assuming the mesons are produced unpolarized. The last uncertainty corresponds to the uncertainty in the Y(1S) meson dimuon branching fraction. The measurement is performed in the final state with four muons using proton-proton collision data collected in 2016 by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9fb⁻¹. This process serves as a standard model reference in a search for narrow resonances decaying to Y(1S)μ⁺μ⁻ in the same final state. Such a resonance could indicate the existence of a tetraquark that is a bound state of two b quarks and two b antiquarks. The tetraquark search is performed for masses in the vicinity of four times the bottom quark mass, between 17.5 and 19.5 GeV, while a generic search for other resonances is performed for masses between 16.5 and 27 GeV. No significant excess of events compatible with a narrow resonance is observed in the data. Limits on the production cross section times branching fraction to four muons via an intermediate Y(1S) resonance are set as a function of the resonance mass