CRISPR-Cas Genome Editing in Ex Vivo Human Lungs to Rewire the Translational Path of Genome-Targeting Therapeutics

The ongoing advancements in CRISPR-Cas technologies can significantly accelerate the preclinical development of both in vivo and ex vivo organ genome-editing therapeutics. One of the promising applications is to genetically modify donor organs prior to implantation. The implantation of optimized donor organs with long-lasting immunomodulatory capacity holds promise for reducing the need for lifelong potent whole-body immunosuppression in recipients. However, assessing genome-targeting interventions in a clinically relevant manner prior to clinical trials remains a major challenge owing to the limited modalities available. This study introduces a novel platform for testing genome editing in human lungs ex vivo, effectively simulating preimplantation genetic engineering of donor organs. We identified gene regulatory elements whose disruption via Cas nucleases led to the upregulation of the immunomodulatory gene interleukin 10 (IL-10). We combined this approach with adenoviral vector-mediated IL-10 delivery to create favorable kinetics for early (immediate postimplantation) graft immunomodulation. Using ex vivo organ machine perfusion and precision-cut tissue slice technology, we demonstrated the feasibility of evaluating CRISPR genome editing in human lungs. To overcome the assessment limitations in ex vivo perfused human organs, we conducted an in vivo rodent study and demonstrated both early gene induction and sustained editing of the lung. Collectively, our findings lay the groundwork for a first-in-human-organ study to overcome the current translational barriers of genome-targeting therapeutics.

Immunomodulation of the donor lung with CRISPR-mediated activation of IL-10 expression

BACKGROUND

Inflammatory injury in the donor lung remains a persistent challenge in lung transplantation that limits donor organ usage and post-transplant outcomes. Inducing immunomodulatory capacity in donor organs could address this unsolved clinical problem. We sought to apply clustered regularly interspaced short palindromic repeats (CRISPR)-associated (Cas) technologies to the donor lung to fine-tune immunomodulatory gene expression, exploring for the first time the therapeutic use of CRISPR-mediated transcriptional activation in the whole donor lung.

METHODS

We explored the feasibility of CRISPR-mediated transcriptional upregulation of interleukin 10 (IL-10), a key immunomodulatory cytokine, in vitro and in vivo. We first evaluated the potency, titratability, and multiplexibility of the gene activation in rat and human cell lines. Next, in vivo CRISPR-mediated IL-10 activation was characterized in rat lungs. Finally, the IL-10-activated donor lungs were transplanted into recipient rats to assess the feasibility in a transplant setting.

RESULTS

The targeted transcriptional activation induced robust and titrable IL-10 upregulation in vitro. The combination of guide RNAs also facilitated multiplex gene modulation, that is, simultaneous activation of IL-10 and IL1 receptor antagonist. In vivo profiling demonstrated that adenoviral delivery of Cas9-based activators to the lung was feasible with the use of immunosuppression, which is routinely applied to organ transplant recipients. The transcriptionally modulated donor lungs retained IL-10 upregulation in isogeneic and allogeneic recipients.

CONCLUSIONS

Our findings highlight the potential of CRISPR epigenome editing to improve lung transplant outcomes by creating a more favorable immunomodulatory environment in the donor organ, a paradigm that may be extendable to other organ transplants.

Odderon Exchange from Elastic Scattering Differences between pp and pp[over ¯] Data at 1.96 TeV and from pp Forward Scattering Measurements

We describe an analysis comparing the pp[over ¯] elastic cross section as measured by the D0 Collaboration at a center-of-mass energy of 1.96 TeV to that in pp collisions as measured by the TOTEM Collaboration at 2.76, 7, 8, and 13 TeV using a model-independent approach. The TOTEM cross sections, extrapolated to a center-of-mass energy of sqrt[s]=1.96  TeV, are compared with the D0 measurement in the region of the diffractive minimum and the second maximum of the pp cross section. The two data sets disagree at the 3.4σ level and thus provide evidence for the t-channel exchange of a colorless, C-odd gluonic compound, also known as the odderon. We combine these results with a TOTEM analysis of the same C-odd exchange based on the total cross section and the ratio of the real to imaginary parts of the forward elastic strong interaction scattering amplitude in pp scattering for which the significance is between 3.4σ and 4.6σ. The combined significance is larger than 5σ and is interpreted as the first observation of the exchange of a colorless, C-odd gluonic compound.

scSNV: accurate dscRNA-seq SNV co-expression analysis using duplicate tag collapsing

Identifying single nucleotide variants has become common practice for droplet-based single-cell RNA-seq experiments; however, presently, a pipeline does not exist to maximize variant calling accuracy. Furthermore, molecular duplicates generated in these experiments have not been utilized to optimally detect variant co-expression. Herein, we introduce scSNV designed from the ground up to "collapse" molecular duplicates and accurately identify variants and their co-expression. We demonstrate that scSNV is fast, with a reduced false-positive variant call rate, and enables the co-detection of genetic variants and A>G RNA edits across twenty-two samples.

Author Correction: Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Morphological classification of pancreatic ductal adenocarcinoma that predicts molecular subtypes and correlates with clinical outcome

INTRODUCTION

Transcriptional analyses have identified several distinct molecular subtypes in pancreatic ductal adenocarcinoma (PDAC) that have prognostic and potential therapeutic significance. However, to date, an indepth, clinicomorphological correlation of these molecular subtypes has not been performed. We sought to identify specific morphological patterns to compare with known molecular subtypes, interrogate their biological significance, and furthermore reappraise the current grading system in PDAC.

DESIGN

We first assessed 86 primary, chemotherapy-naive PDAC resection specimens with matched RNA-Seq data for specific, reproducible morphological patterns. Differential expression was applied to the gene expression data using the morphological features. We next compared the differentially expressed gene signatures with previously published molecular subtypes. Overall survival (OS) was correlated with the morphological and molecular subtypes.

RESULTS

We identified four morphological patterns that segregated into two components ('gland forming' and 'non-gland forming') based on the presence/absence of well-formed glands. A morphological cut-off (≥40% 'non-gland forming') was established using RNA-Seq data, which identified two groups (A and B) with gene signatures that correlated with known molecular subtypes. There was a significant difference in OS between the groups. The morphological groups remained significantly prognostic within cancers that were moderately differentiated and classified as 'classical' using RNA-Seq.

CONCLUSION

Our study has demonstrated that PDACs can be morphologically classified into distinct and biologically relevant categories which predict known molecular subtypes. These results provide the basis for an improved taxonomy of PDAC, which may lend itself to future treatment strategies and the development of deep learning models.

Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution

Pancreatic adenocarcinoma presents as a spectrum of a highly aggressive disease in patients. The basis of this disease heterogeneity has proved difficult to resolve due to poor tumor cellularity and extensive genomic instability. To address this, a dataset of whole genomes and transcriptomes was generated from purified epithelium of primary and metastatic tumors. Transcriptome analysis demonstrated that molecular subtypes are a product of a gene expression continuum driven by a mixture of intratumoral subpopulations, which was confirmed by single-cell analysis. Integrated whole-genome analysis uncovered that molecular subtypes are linked to specific copy number aberrations in genes such as mutant KRAS and GATA6. By mapping tumor genetic histories, tetraploidization emerged as a key mutational process behind these events. Taken together, these data support the premise that the constellation of genomic aberrations in the tumor gives rise to the molecular subtype, and that disease heterogeneity is due to ongoing genomic instability during progression.

Measurement of the Effective Weak Mixing Angle in pp[over ¯]→Z/γ^{*}→ℓ^{+}ℓ^{-} Events

We present a measurement of the effective weak mixing angle parameter sin^{2}θ_{eff}^{ℓ} in pp[over ¯]→Z/γ^{*}→μ^{+}μ^{-} events at a center-of-mass energy of 1.96 TeV, collected by the D0 detector at the Fermilab Tevatron Collider and corresponding to 8.6  fb^{-1} of integrated luminosity. The measured value of sin^{2}θ_{eff}^{ℓ}[μμ]=0.23016±0.00064 is further combined with the result from the D0 measurement in pp[over ¯]→Z/γ^{*}→e^{+}e^{-} events, resulting in sin^{2}θ_{eff}^{ℓ}[comb]=0.23095±0.00040. This combined result is the most precise measurement from a single experiment at a hadron collider and is the most precise determination using the coupling of the Z/γ^{*} to light quarks.

Combined Forward-Backward Asymmetry Measurements in Top-Antitop Quark Production at the Tevatron

The CDF and D0 experiments at the Fermilab Tevatron have measured the asymmetry between yields of forward- and backward-produced top and antitop quarks based on their rapidity difference and the asymmetry between their decay leptons. These measurements use the full data sets collected in proton-antiproton collisions at a center-of-mass energy of sqrt[s]=1.96  TeV. We report the results of combinations of the inclusive asymmetries and their differential dependencies on relevant kinematic quantities. The combined inclusive asymmetry is A_{FB}^{tt[over ¯]}=0.128±0.025. The combined inclusive and differential asymmetries are consistent with recent standard model predictions.

Millimeter-Wave Polarimeters Using Kinetic Inductance Detectors for TolTEC and Beyond

Microwave Kinetic Inductance Detectors (MKIDs) provide a compelling path forward to the large-format polarimeter, imaging, and spectrometer arrays needed for next-generation experiments in millimeter-wave cosmology and astronomy. We describe the development of feedhorn-coupled MKID detectors for the TolTEC millimeter-wave imaging polarimeter being constructed for the 50-meter Large Millimeter Telescope (LMT). Observations with TolTEC are planned to begin in early 2019. TolTEC will comprise ∼7,000 polarization sensitive MKIDs and will represent the first MKID arrays fabricated and deployed on monolithic 150 mm diameter silicon wafers - a critical step towards future large-scale experiments with over 10⁵ detectors. TolTEC will operate in observational bands at 1.1, 1.4, and 2.0 mm and will use dichroic filters to define a physically independent focal plane for each passband, thus allowing the polarimeters to use simple, direct-absorption inductive structures that are impedance matched to incident radiation. This work is part of a larger program at NIST-Boulder to develop MKID-based detector technologies for use over a wide range of photon energies spanning millimeter-waves to X-rays. We present the detailed pixel layout and describe the methods, tools, and flexible design parameters that allow this solution to be optimized for use anywhere in the millimeter and sub-millimeter bands. We also present measurements of prototype devices operating in the 1.1 mm band and compare the observed optical performance to that predicted from models and simulations.

Association of Distinct Mutational Signatures With Correlates of Increased Immune Activity in Pancreatic Ductal Adenocarcinoma

IMPORTANCE

Outcomes for patients with pancreatic ductal adenocarcinoma (PDAC) remain poor. Advances in next-generation sequencing provide a route to therapeutic approaches, and integrating DNA and RNA analysis with clinicopathologic data may be a crucial step toward personalized treatment strategies for this disease.

OBJECTIVE

To classify PDAC according to distinct mutational processes, and explore their clinical significance.

DESIGN, SETTING, AND PARTICIPANTS

We performed a retrospective cohort study of resected PDAC, using cases collected between 2008 and 2015 as part of the International Cancer Genome Consortium. The discovery cohort comprised 160 PDAC cases from 154 patients (148 primary; 12 metastases) that underwent tumor enrichment prior to whole-genome and RNA sequencing. The replication cohort comprised 95 primary PDAC cases that underwent whole-genome sequencing and expression microarray on bulk biospecimens.

MAIN OUTCOMES AND MEASURES

Somatic mutations accumulate from sequence-specific processes creating signatures detectable by DNA sequencing. Using nonnegative matrix factorization, we measured the contribution of each signature to carcinogenesis, and used hierarchical clustering to subtype each cohort. We examined expression of antitumor immunity genes across subtypes to uncover biomarkers predictive of response to systemic therapies.

RESULTS

The discovery cohort was 53% male (n = 79) and had a median age of 67 (interquartile range, 58-74) years. The replication cohort was 50% male (n = 48) and had a median age of 68 (interquartile range, 60-75) years. Five predominant mutational subtypes were identified that clustered PDAC into 4 major subtypes: age related, double-strand break repair, mismatch repair, and 1 with unknown etiology (signature 8). These were replicated and validated. Signatures were faithfully propagated from primaries to matched metastases, implying their stability during carcinogenesis. Twelve of 27 (45%) double-strand break repair cases lacked germline or somatic events in canonical homologous recombination genes-BRCA1, BRCA2, or PALB2. Double-strand break repair and mismatch repair subtypes were associated with increased expression of antitumor immunity, including activation of CD8-positive T lymphocytes (GZMA and PRF1) and overexpression of regulatory molecules (cytotoxic T-lymphocyte antigen 4, programmed cell death 1, and indolamine 2,3-dioxygenase 1), corresponding to higher frequency of somatic mutations and tumor-specific neoantigens.

CONCLUSIONS AND RELEVANCE

Signature-based subtyping may guide personalized therapy of PDAC in the context of biomarker-driven prospective trials.

A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns

Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.

Evidence for a B_{s}^{0}π^{±} State

We report evidence for a narrow structure, X(5568), in the decay sequence X(5568)→B_{s}^{0}π^{±}, B_{s}^{0}→J/ψϕ, J/ψ→μ^{+}μ^{-}, ϕ→K^{+}K^{-}. This is evidence for the first instance of a hadronic state with valence quarks of four different flavors. The mass and natural width of this state are measured to be m=5567.8±2.9(stat)_{-1.9}^{+0.9}(syst)  MeV/c^{2} and Γ=21.9±6.4(stat)_{-2.5}^{+5.0}(syst)  MeV/c^{2}. If the decay is X(5568)→B_{s}^{*}π^{±}→B_{s}^{0}γπ^{±} with an unseen γ, m(X(5568)) will be shifted up by m(B_{s}^{*})-m(B_{s}^{0})∼49  MeV/c^{2}. This measurement is based on 10.4  fb^{-1} of pp[over ¯] collision data at sqrt[s]=1.96  TeV collected by the D0 experiment at the Fermilab Tevatron collider.

Evidence for Simultaneous Production of J/ψ and ϒ Mesons

We report evidence for the simultaneous production of J/ψ and ϒ mesons in 8.1  fb^{-1} of data collected at sqrt[s]=1.96  TeV by the D0 experiment at the Fermilab pp[over ¯] Tevatron Collider. Events with these characteristics are expected to be produced predominantly by gluon-gluon interactions. In this analysis, we extract the effective cross section characterizing the initial parton spatial distribution, σ_{eff}=2.2±0.7(stat)±0.9(syst)  mb.

Inclusive Production of the X(4140) State in pp[over ¯] Collisions at D0

We present a study of the inclusive production of the X(4140) state with the decay to the J/ψϕ final state in hadronic collisions. Based on 10.4  fb^{-1} of pp[over ¯] collision data collected by the D0 experiment at the Fermilab Tevatron collider, we report the first evidence for the prompt production of an X(4140) state and find the fraction of X(4140) events originating from b hadrons to be f_{b}=0.39±0.07(stat)±0.10(syst). The ratio of the nonprompt X(4140) production rate to the B_{s}^{0} yield in the same channel is R=0.19±0.05(stat)±0.07(syst). The values of the mass M=4152.5±1.7(stat)_{-5.4}^{+6.2}(syst)  MeV and width Γ=16.3±5.6(stat)±11.4(syst)  MeV are consistent with previous measurements.

Search for Violation of CPT and Lorentz Invariance in Bs(0) Meson Oscillations

We present the first search for CPT-violating effects in the mixing of Bs(0) mesons using the full Run II data set with an integrated luminosity of 10.4  fb(-1) of proton-antiproton collisions collected using the D0 detector at the Fermilab Tevatron Collider. We measure the CPT-violating asymmetry in the decay Bs(0)→μ(±)Ds(±) as a function of celestial direction and sidereal phase. We find no evidence for CPT-violating effects and place limits on the direction and magnitude of flavor-dependent CPT- and Lorentz-invariance violating coupling coefficients. We find 95% confidence intervals of Δa⊥<1.2×10(-12)  GeV and (-0.8<ΔaT-0.396ΔaZ<3.9)×10(-13)  GeV.

Tevatron Combination of Single-Top-Quark Cross Sections and Determination of the Magnitude of the Cabibbo-Kobayashi-Maskawa Matrix Element V_{tb}

We present the final combination of CDF and D0 measurements of cross sections for single-top-quark production in proton-antiproton collisions at a center-of-mass energy of 1.96 TeV. The data correspond to total integrated luminosities of up to 9.7 fb^{-1} per experiment. The t-channel cross section is measured to be σ_{t}=2.25_{-0.31}^{+0.29} pb. We also present the combinations of the two-dimensional measurements of the s- vs t-channel cross section. In addition, we give the combination of the s+t channel cross section measurement resulting in σ_{s+t}=3.30_{-0.40}^{+0.52} pb, without assuming the standard model value for the ratio σ_{s}/σ_{t}. The resulting value of the magnitude of the top-to-bottom quark coupling is |V_{tb}|=1.02_{-0.05}^{+0.06}, corresponding to |V_{tb}|>0.92 at the 95% C.L.

Measurement of the Effective Weak Mixing Angle in pp[over ¯]→Z/γ^{*}→e^{+}e^{-} Events

We present a measurement of the fundamental parameter of the standard model, the weak mixing angle sin^{2}θ_{eff}^{ℓ} which determines the relative strength of weak and electromagnetic interactions, in pp[over ¯]→Z/γ^{*}→e^{+}e^{-} events at a center of mass energy of 1.96 TeV, using data corresponding to 9.7 fb^{-1} of integrated luminosity collected by the D0 detector at the Fermilab Tevatron. The effective weak mixing angle is extracted from the forward-backward charge asymmetry as a function of the invariant mass around the Z boson pole. The measured value of sin^{2}θ_{eff}^{ℓ}=0.23147±0.00047 is the most precise measurement from light quark interactions to date, with a precision close to the best LEP and SLD results.

RNASequel: accurate and repeat tolerant realignment of RNA-seq reads

RNA-seq is a key technology for understanding the biology of the cell because of its ability to profile transcriptional and post-transcriptional regulation at single nucleotide resolutions. Compared to DNA sequencing alignment algorithms, RNA-seq alignment algorithms have a diminished ability to accurately detect and map base pair substitutions, gaps, discordant pairs and repetitive regions. These shortcomings adversely affect experiments that require a high degree of accuracy, notably the ability to detect RNA editing. We have developed RNASequel, a software package that runs as a post-processing step in conjunction with an RNA-seq aligner and systematically corrects common alignment artifacts. Its key innovations are a two-pass splice junction alignment system that includes de novo splice junctions and the use of an empirically determined estimate of the fragment size distribution when resolving read pairs. We demonstrate that RNASequel produces improved alignments when used in conjunction with STAR or Tophat2 using two simulated datasets. We then show that RNASequel improves the identification of adenosine to inosine RNA editing sites on biological datasets. This software will be useful in applications requiring the accurate identification of variants in RNA sequencing data, the discovery of RNA editing sites and the analysis of alternative splicing.

Tevatron constraints on models of the Higgs boson with exotic spin and parity using decays to bottom-antibottom quark pairs

Combined constraints from the CDF and D0 Collaborations on models of the Higgs boson with exotic spin J and parity P are presented and compared with results obtained assuming the standard model value JP=0+. Both collaborations analyzed approximately 10  fb(-) of proton-antiproton collisions with a center-of-mass energy of 1.96 TeV collected at the Fermilab Tevatron. Two models predicting exotic Higgs bosons with JP=0- and JP=2+ are tested. The kinematic properties of exotic Higgs boson production in association with a vector boson differ from those predicted for the standard model Higgs boson. Upper limits at the 95% credibility level on the production rates of the exotic Higgs bosons, expressed as fractions of the standard model Higgs boson production rate, are set at 0.36 for both the JP=0- hypothesis and the JP=2+ hypothesis. If the production rate times the branching ratio to a bottom-antibottom pair is the same as that predicted for the standard model Higgs boson, then the exotic bosons are excluded with significances of 5.0 standard deviations and 4.9 standard deviations for the JP=0- and JP=2+ hypotheses, respectively.

Activation of an endogenous retrovirus-associated long non-coding RNA in human adenocarcinoma

Background

Long non-coding RNAs (lncRNAs) are emerging as molecules that significantly impact many cellular processes and have been associated with almost every human cancer. Compared to protein-coding genes, lncRNA genes are often associated with transposable elements, particularly with endogenous retroviral elements (ERVs). ERVs can have potentially deleterious effects on genome structure and function, so these elements are typically silenced in normal somatic tissues, albeit with varying efficiency. The aberrant regulation of ERVs associated with lncRNAs (ERV-lncRNAs), coupled with the diverse range of lncRNA functions, creates significant potential for ERV-lncRNAs to impact cancer biology.

Methods

We used RNA-seq analysis to identify and profile the expression of a novel lncRNA in six large cohorts, including over 7,500 samples from The Cancer Genome Atlas (TCGA).

Results

We identified the tumor-specific expression of a novel lncRNA that we have named Endogenous retroViral-associated ADenocarcinoma RNA or ‘EVADR’, by analyzing RNA-seq data derived from colorectal tumors and matched normal control tissues. Subsequent analysis of TCGA RNA-seq data revealed the striking association of EVADR with adenocarcinomas, which are tumors of glandular origin. Moderate to high levels of EVADR were detected in 25 to 53% of colon, rectal, lung, pancreas and stomach adenocarcinomas (mean = 30 to 144 FPKM), and EVADR expression correlated with decreased patient survival (Cox regression; hazard ratio = 1.47, 95% confidence interval = 1.06 to 2.04, P = 0.02). In tumor sites of non-glandular origin, EVADR expression was detectable at only very low levels and in less than 10% of patients. For EVADR, a MER48 ERV element provides an active promoter to drive its transcription. Genome-wide, MER48 insertions are associated with nine lncRNAs, but none of the MER48-associated lncRNAs other than EVADR were consistently expressed in adenocarcinomas, demonstrating the specific activation of EVADR. The sequence and structure of the EVADR locus is highly conserved among Old World monkeys and apes but not New World monkeys or prosimians, where the MER48 insertion is absent. Conservation of the EVADR locus suggests a functional role for this novel lncRNA in humans and our closest primate relatives.

Conclusions

Our results describe the specific activation of a highly conserved ERV-lncRNA in numerous cancers of glandular origin, a finding with diagnostic, prognostic and therapeutic implications.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-015-0142-6) contains supplementary material, which is available to authorized users.

Measurement of the B(s)(0) lifetime in the flavor-specific decay channel B(s)(0)→D(s)(-)μ(+)νX

We present an updated measurement of the B(s)(0) lifetime using the semileptonic decays B(s)(0)→D(s)(-)μ(+)νX, with D(s)(-)→ϕπ(-) and ϕ→K(+)K(-) (and the charge conjugate process). This measurement uses the full Tevatron Run II sample of proton-antiproton collisions at √[s]=1.96  TeV, comprising an integrated luminosity of 10.4  fb(-1). We find a flavor-specific lifetime τ(fs)(B(s)(0))=1.479±0.010(stat)±0.021(syst)  ps. This technique is also used to determine the B(0) lifetime using the analogous B(0)→D(-)μ(+)νX decay with D(-)→ϕπ(-) and ϕ→K(+)K(-), yielding τ(B(0))=1.534±0.019(stat)±0.021(syst)  ps. Both measurements are consistent with the current world averages, and the B(s)(0) lifetime measurement is one of the most precise to date. Taking advantage of the cancellation of systematic uncertainties, we determine the lifetime ratio τ(fs)(B(s)(0))/τ(B(0))=0.964±0.013(stat)±0.007(syst).

Measurement of the forward-backward asymmetry in the production of B(±) mesons in pp̄ collisions at √s=1.96 TeV

We present a measurement of the forward-backward asymmetry in the production of B(±) mesons, A(FB)(B(±)), using B(±)→J/ψK(±) decays in 10.4  fb(-1) of pp̄ collisions at sqrt[s]=1.96  TeV collected by the D0 experiment during Run II of the Tevatron collider. A nonzero asymmetry would indicate a preference for a particular flavor, i.e., b quark or ̄b antiquark, to be produced in the direction of the proton beam. We extract A(FB)(B(±)) from a maximum likelihood fit to the difference between the numbers of forward- and backward-produced B(±) mesons. We measure an asymmetry consistent with zero: A(FB)(B(±))=[-0.24±0.41  (stat)±0.19  (syst)]%.

Constraints on models for the Higgs boson with exotic spin and parity in VH → Vbb final states

We present constraints on models containing non-standard-model values for the spin J and parity P of the Higgs boson H in up to 9.7 fb(-1) of pp collisions at sqrt[s] = 1.96 TeV collected with the D0 detector at the Fermilab Tevatron Collider. These are the first studies of Higgs boson J(P) with fermions in the final state. In the ZH → ℓℓbb, WH → ℓνbb, and ZH → ννbb final states, we compare the standard model (SM) Higgs boson prediction, J(P) = 0(+), with two alternative hypotheses, J(P) = 0(-) and J(P) = 2(+). We use a likelihood ratio to quantify the degree to which our data are incompatible with non-SM J(P) predictions for a range of possible production rates. Assuming that the production rate in the signal models considered is equal to the SM prediction, we reject the J(P) = 0(-) and J(P) = 2(+) hypotheses at the 97.6% CL and at the 99.0% CL, respectively. The expected exclusion sensitivity for a J(P) = 0(-) (J(P) = 2(+)) state is at the 99.86% (99.94%) CL. Under the hypothesis that our data are the result of a combination of the SM-like Higgs boson and either a J(P) = 0(-) or a J(P) = 2(+) signal, we exclude a J(P) = 0(-) fraction above 0.80 and a J(P) = 2(+) fraction above 0.67 at the 95% CL. The expected exclusion covers J(P) = 0(-) (J(P) = 2(+)) fractions above 0.54 (0.47).

Precision measurement of the top quark mass in lepton + jets final States

We measure the mass of the top quark in lepton+jets final states using the full sample of pp collision data collected by the D0 experiment in Run II of the Fermilab Tevatron Collider at sqrt[s] = 1.96 TeV, corresponding to 9.7 fb(-1) of integrated luminosity. We use a matrix element technique that calculates the probabilities for each event to result from tt production or background. The overall jet energy scale is constrained in situ by the mass of the W boson. We measure m(t) = 174.98 ± 0.76 GeV. This constitutes the most precise single measurement of the top-quark mass.

Observation of s-channel production of single top quarks at the Tevatron

We report the first observation of single-top-quark production in the s channel through the combination of the CDF and D0 measurements of the cross section in proton-antiproton collisions at a center-of-mass energy of 1.96 TeV. The data correspond to total integrated luminosities of up to 9.7 fb(-1) per experiment. The measured cross section is σ(s) = 1.29(-0.24)(+0.26) pb. The probability of observing a statistical fluctuation of the background to a cross section of the observed size or larger is 1.8 × 10(-10), corresponding to a significance of 6.3 standard deviations for the presence of an s-channel contribution to the production of single-top quarks.

Measurement of the W boson production charge asymmetry in pp̄→W+X→eν+X events at √s=1.96 TeV

We present a measurement of the W boson production charge asymmetry in pp̄→W+X→eν+X events at a center of mass energy of 1.96 TeV, using 9.7  fb-1 of integrated luminosity collected with the D0 detector at the Fermilab Tevatron Collider. The neutrino longitudinal momentum is determined by using a neutrino weighting method, and the asymmetry is measured as a function of the W boson rapidity. The measurement extends over wider electron pseudorapidity region than previous results and is the most precise to date, allowing for precise determination of proton parton distribution functions in global fits.

Measurement of the direct CP-violating charge asymmetry in Ds ± → ϕ π ± decays

We present a measurement of the direct CP-violating charge asymmetry in Ds+ → ϕ π(±) decays where the ϕ meson decays to K(+)K(-), using the full Run II data set with an integrated luminosity of 10.4  fb-1 of proton-antiproton collisions collected using the D0 detector at the Fermilab Tevatron Collider. The normalized difference A(CP) in the yield of D(s)(+) and D(s)(-) mesons in these decays is measured by fitting the difference between their reconstructed invariant mass distributions. This results in an asymmetry of A(CP) = [-0.38 ± 0.27]%, which is the most precise measurement of this quantity to date. The result is consistent with the standard model prediction of zero CP asymmetry in this decay.

Measurement of associated production of z bosons with charm quark jets in pp collisions at √s=1.96 TeV

We present the first measurements of the ratios of cross sections σ(pp → Z+c jet)/σ(pp → Z+jet) and σ(pp → Z+c jet)/σ(pp → Z+b jet) for the associated production of a Z boson with at least one charm or bottom quark jet. Jets have transverse momentum ​pT(jet)​>20  GeV and pseudorapidity |​η(jet)​|<2.5. These cross section ratios are measured differentially as a function of jet and Z boson transverse momenta, based on 9.7  fb(-1) of pp collisions collected with the D0 detector at the Fermilab Tevatron Collider at √s=1.96  TeV. The measurements show significant deviations from perturbative QCD calculations and predictions from various event generators.

Measurement of direct CP violation parameters in B± → J/ψK± and B± → J/ψπ± decays with 10.4 fb-1 of Tevatron data

We present a measurement of the direct CP-violating charge asymmetry in B(±) mesons decaying to J/ψK(±) and J/ψπ(±) where J/ψ decays to μ(+) μ(-), using the full run II data set of 10.4 fb(-1) of proton-antiproton collisions collected using the D0 detector at the Fermilab Tevatron Collider. A difference in the yield of B(-) and B(+) mesons in these decays is found by fitting to the difference between their reconstructed invariant mass distributions resulting in asymmetries of A(J/ψK) = [0.59 ± 0.37]%, which is the most precise measurement to date, and A(J/ψπ) = [-4.2 ± 4.5]%. Both measurements are consistent with standard model predictions.

Measurement of the semileptonic charge asymmetry using B(s)(0) → D(s)μX decays

We present a measurement of the time-integrated flavor-specific semileptonic charge asymmetry in the decays of B(s)(0) mesons that have undergone flavor mixing, a(sl)(s), using B(s)(0)(B(s)(0)) → D(s)(-/+) μ(±) X decays, with D(s)(-/+) → [symbol: see text]π(-/+) and [symbol: see text]→ K(+) K(-), using 10.4 fb(-1) of proton-antiproton collisions collected by the D0 detector during Run II at the Fermilab Tevatron Collider. A fit to the difference between the time-integrated D(s)(-) and D(s)(+) mass distributions of the B(s)(0) and B(s)(0) candidates yields the flavor-specific asymmetry a(sl)(s) = [-1.12 ± 0.74(stat) ± 0.17(syst)]%, which is the most precise measurement and in agreement with the standard model prediction.

Search for the standard model Higgs boson in associated WH production in 9.7 fb(-1) of pp collisions with the D0 detector

We present a search for the standard model Higgs boson in final states with a charged lepton (electron or muon), missing transverse energy, and two or three jets, at least one of which is identified as a b-quark jet. The search is primarily sensitive to WH→ℓνbb production and uses data corresponding to 9.7 fb(-1) of integrated luminosity collected with the D0 detector at the Fermilab Tevatron pp Collider at √s = 1.96 TeV. We observe agreement between the data and the expected background. For a Higgs boson mass of 125 GeV, we set a 95% C.L. upper limit on the production of a standard model Higgs boson of 5.2 × σ(SM), where σ(SM) is the standard model Higgs boson production cross section, while the expected limit is 4.7 × σ(SM).

Combined search for the standard model Higgs boson decaying to bb using the D0 run II data set

We present the results of the combination of searches for the standard model Higgs boson produced in association with a W or Z boson and decaying into bb using the data sample collected with the D0 detector in pp collisions at √s = 1.96 TeV at the Fermilab Tevatron Collider. We derive 95% C.L. upper limits on the Higgs boson cross section relative to the standard model prediction in the mass range 100 GeV ≤ M(H) ≤ 150 GeV, and we exclude Higgs bosons with masses smaller than 102 GeV at the 95% C.L. In the mass range 120 GeV ≤ M(H) ≤145 GeV, the data exhibit an excess above the background prediction with a global significance of 1.5 standard deviations, consistent with the expectation in the presence of a standard model Higgs boson.

Search for the standard model Higgs boson in ZH→ℓ(+)ℓ(-)bb production with the D0 detector in 9.7 fb(-1) of pp collisions at √s = 1.96 TeV

We present a search for the standard model (SM) Higgs boson produced in association with a Z boson in 9.7 fb(-1) of pp collisions collected with the D0 detector at the Fermilab Tevatron Collider at √s = 1.96 TeV. Selected events contain one reconstructed Z→e(+)e(-) or Z→μ(+)μ(-) candidate and at least two jets, including at least one jet identified as likely to contain a b quark. To validate the search procedure, we also measure the cross section for ZZ production in the same final state. It is found to be consistent with its SM prediction. We set upper limits on the ZH production cross section times branching ratio for H→bb at the 95% C.L. for Higgs boson masses 90 ≤ M(H) ≤ 150 GeV. The observed (expected) limit for M(H) = 125 GeV is 7.1 (5.1) times the SM cross section.

Evidence for a particle produced in association with weak bosons and decaying to a bottom-antibottom quark pair in higgs boson searches at the tevatron

We combine searches by the CDF and D0 Collaborations for the associated production of a Higgs boson with a W or Z boson and subsequent decay of the Higgs boson to a bottom-antibottom quark pair. The data, originating from Fermilab Tevatron pp collisions at √s = 1.96 TeV, correspond to integrated luminosities of up to 9.7 fb(-1). The searches are conducted for a Higgs boson with mass in the range 100-150 GeV/c(2). We observe an excess of events in the data compared with the background predictions, which is most significant in the mass range between 120 and 135 GeV/c(2). The largest local significance is 3.3 standard deviations, corresponding to a global significance of 3.1 standard deviations. We interpret this as evidence for the presence of a new particle consistent with the standard model Higgs boson, which is produced in association with a weak vector boson and decays to a bottom-antibottom quark pair.

Search for violation of Lorentz invariance in top quark pair production and decay

Using data collected with the D0 detector at the Fermilab Tevatron Collider, corresponding to 5.3 fb(-1) of integrated luminosity, we search for violation of Lorentz invariance by examining the tt[over ¯] production cross section in lepton+jets final states. We quantify this violation using the standard-model extension framework, which predicts a dependence of the tt[over ¯] production cross section on sidereal time as the orientation of the detector changes with the rotation of the Earth. Within this framework, we measure components of the matrices (c(Q))(μν33) and (c(U))(μν33) containing coefficients used to parametrize violation of Lorentz invariance in the top quark sector. Within uncertainties, these coefficients are found to be consistent with zero.

Measurements of WW and WZ production in W + jets final states in pp collisions

We study WW and WZ production with ℓνqq (ℓ=e,μ) final states using data collected by the D0 detector at the Fermilab Tevatron Collider corresponding to 4.3  fb(-1) of integrated luminosity from pp collisions at sqrt[s]=1.96  TeV. Assuming the ratio between the production cross sections σ(WW) and σ(WZ) as predicted by the standard model, we measure the total WV (V=W,Z) cross section to be σ(WV)=19.6(-3.0)(+3.2)  pb and reject the background-only hypothesis at a level of 7.9 standard deviations. We also use b-jet discrimination to separate the WZ component from the dominant WW component. Simultaneously fitting WW and WZ contributions, we measure σ(WW)=15.9(-3.2)(+3.7)  pb and σ(WZ)=3.3(-3.3)(+4.1)  pb, which is consistent with the standard model predictions.

Measurement of the W boson mass with the D0 detector

We present a measurement of the W boson mass using data corresponding to 4.3  fb(-1) of integrated luminosity collected with the D0 detector during Run II at the Fermilab Tevatron pp collider. With a sample of 1,677,394 W → eν candidate events, we measure M(W) = 80.367 ± 0.026  GeV. This result is combined with an earlier D0 result determined using an independent Run II data sample, corresponding to 1  fb(-1) of integrated luminosity, to yield M(W) = 80.375 ± 0.023  GeV.

Search for universal extra dimensions in pp collisions

We present a search for Kaluza-Klein (KK) particles predicted by models with universal extra dimensions (UED) using a data set corresponding to an integrated luminosity of 7.3 fb(-1), collected by the D0 detector at a pp center-of-mass energy of 1.96 TeV. The decay chain of KK particles can lead to a final state with two muons of the same charge. This signature is used to set a lower limit on the compactification scale of R(-1)>260 GeV in a minimal UED model.

Search for charged massive long-lived particles

We report on a search for charged massive long-lived particles (CMLLPs), based on 5.2  fb(-1) of integrated luminosity collected with the D0 detector at the Fermilab Tevatron pp collider. We search for events in which one or more particles are reconstructed as muons but have speed and ionization energy loss (dE/dx) inconsistent with muons produced in beam collisions. CMLLPs are predicted in several theories of physics beyond the standard model. We exclude pair-produced long-lived gauginolike charginos below 267 GeV and Higgsino-like charginos below 217 GeV at 95% C.L., as well as long-lived scalar top quarks with mass below 285 GeV.

Evidence for spin correlation in t ̅t production

We present a measurement of the ratio of events with correlated t and t ̅t spins to the total number of t ̅t events. This ratio f is evaluated using a matrix-element-based approach in 729 t ̅t candidate events with a single lepton ℓ (electron or muon) and at least four jets. The analyzed p ̅p collisions data correspond to an integrated luminosity of 5.3 fb(-1) and were collected with the D0 detector at the Fermilab Tevatron collider operating at a center-of-mass energy √s=1.96 TeV. Combining this result with a recent measurement of f in dileptonic final states, we find f in agreement with the standard model. In addition, the combination provides evidence for the presence of spin correlation in t ̅t events with a significance of more than 3 standard deviations.

Search for doubly charged Higgs boson pair production in pp¯ collisions at √s=1.96 TeV

We present a search for pair production of doubly charged Higgs bosons in the processes qq¯→H(++)H(--) decaying through H(±±)→τ(±)τ(±),μ(±)τ(±),μ(±)μ(±). The search is performed in pp¯ collisions at a center-of-mass energy of √s = 1.96 TeV using an integrated luminosity of up to 7.0 fb(-1) collected by the D0 experiment at the Fermilab Tevatron Collider. The results are used to set 95% C.L. limits on the pair production cross section of doubly charged Higgs bosons and on their mass for different H(±±) branching fractions. Models predicting different H(±±) decays are investigated. Assuming B(H(±±)→τ(±)τ(±))=1 yields an observed (expected) lower limit on the mass of a left-handed H(L)(±±) boson of 128 (116) GeV and assuming B(H(±±)→μ(±)τ(±))=1 the corresponding limits are 144 (149) GeV. In a model with B(H(±±)→τ(±)τ(±))=B(H(±±)→μ(±)τ(±))=B(H(±±)→μ(±)μ(±))=1/3, we obtain M(H(L)(±±))>130 (138) GeV.

Wγ production and limits on anomalous WWγ couplings in pp collisions at sqrt[s] = 1.96 TeV

We measure the cross section and the difference in rapidities between photons and charged leptons for inclusive W(→lν) + γ production in eγ and μγ final states. Using data corresponding to an integrated luminosity of 4.2  fb(-1) collected with the D0 detector at the Fermilab Tevatron Collider, the measured cross section times branching fraction for the process pp → Wγ + X → lνγ + X and the distribution of the charge-signed photon-lepton rapidity difference are found to be in agreement with the standard model. These results provide the most stringent limits on anomalous WWγ couplings for data from hadron colliders: -0.4<Δκ(γ)<0.4 and -0.08<λ(γ)<0.07 at the 95% C.L.

Search for a fermiophobic and standard model Higgs boson in diphoton final states

We present a search for the standard model Higgs boson and a fermiophobic Higgs boson in the diphoton final states based on 8.2  fb(-1) of pp collisions at sqrt[s]=1.96  TeV collected with the D0 detector at the Fermilab Tevatron Collider. No excess of data above background predictions is observed and upper limits at the 95% C.L. on the cross section multiplied by the branching fraction are set which are the most restrictive to date. A fermiophobic Higgs boson with a mass below 112.9 GeV is excluded at the 95% C.L.

Search for neutral minimal supersymmetric standard model Higgs bosons decaying to tau pairs produced in association with b quarks in pp collisions at √s = 1.96 TeV

We report results from a search for neutral Higgs bosons produced in association with b quarks using data recorded by the D0 experiment at the Fermilab Tevatron Collider and corresponding to an integrated luminosity of 7.3 fb(-1). This production mode can be enhanced in several extensions of the standard model (SM) such as in its minimal supersymmetric extension (MSSM) at high tanβ. We search for Higgs bosons decaying to tau pairs with one tau decaying to a muon and neutrinos and the other to hadrons. The data are found to be consistent with SM expectations, and we set upper limits on the cross section times branching ratio in the Higgs boson mass range from 90 to 320 GeV/c(2). We interpret our result in the MSSM parameter space, excluding tanβ values down to 25 for Higgs boson masses below 170 GeV/c(2).