A Robust and Powerful Set-Valued Approach to Rare Variant Association Analyses of Secondary Traits in Case-Control Sequencing Studies

In many case-control designs of genome-wide association (GWAS) or next generation sequencing (NGS) studies, extensive data on secondary traits that may correlate and share the common genetic variants with the primary disease are available. Investigating these secondary traits can provide critical insights into the disease etiology or pathology, and enhance the GWAS or NGS results. Methods based on logistic regression (LG) were developed for this purpose. However, for the identification of rare variants (RVs), certain inadequacies in the LG models and algorithmic instability can cause severely inflated type I error, and significant loss of power, when the two traits are correlated and the RV is associated with the disease, especially at stringent significance levels. To address this issue, we propose a novel set-valued (SV) method that models a binary trait by dichotomization of an underlying continuous variable, and incorporate this into the genetic association model as a critical component. Extensive simulations and an analysis of seven secondary traits in a GWAS of benign ethnic neutropenia show that the SV method consistently controls type I error well at stringent significance levels, has larger power than the LG-based methods, and is robust in performance to effect pattern of the genetic variant (risk or protective), rare or common variants, rare or common diseases, and trait distributions. Because of the SV method’s striking and profound advantage, we strongly recommend the SV method be employed instead of the LG-based methods for secondary traits analyses in case-control sequencing studies

Multi-organ Mapping of Cancer Risk.

Cancers are distributed unevenly across the body, but the importance of cell intrinsic factors such as stem cell function in determining organ cancer risk is unknown. Therefore, we used Cre-recombination of conditional lineage tracing, oncogene, and tumor suppressor alleles to define populations of stem and non-stem cells in mouse organs and test their life-long susceptibility to tumorigenesis. We show that tumor incidence is determined by the life-long generative capacity of mutated cells. This relationship held true in the presence of multiple genotypes and regardless of developmental stage, strongly supporting the notion that stem cells dictate organ cancer risk. Using the liver as a model system, we further show that damage-induced activation of stem cell function markedly increases cancer risk. Therefore, we propose that a combination of stem cell mutagenesis and extrinsic factors that enhance the proliferation of these cell populations, creates a "perfect storm" that ultimately determines organ cancer risk. VIDEO ABSTRACT.National Institutes of Health (Grant IDs: P01CA96832, R01, P30CA021765); the American Lebanese Syrian Associated Charities; Cancer Research UKThis is the author accepted manuscript. The final version is available from Elsevier (Cell Press) via http://dx.doi.org/10.1016/j.cell.2016.07.04

A novel retinoblastoma therapy from genomic and epigenetic analyses.

Retinoblastoma is an aggressive childhood cancer of the developing retina that is initiated by the biallelic loss of RB1. Tumours progress very quickly following RB1 inactivation but the underlying mechanism is not known. Here we show that the retinoblastoma genome is stable, but that multiple cancer pathways can be epigenetically deregulated. To identify the mutations that cooperate with RB1 loss, we performed whole-genome sequencing of retinoblastomas. The overall mutational rate was very low; RB1 was the only known cancer gene mutated. We then evaluated the role of RB1 in genome stability and considered non-genetic mechanisms of cancer pathway deregulation. For example, the proto-oncogene SYK is upregulated in retinoblastoma and is required for tumour cell survival. Targeting SYK with a small-molecule inhibitor induced retinoblastoma tumour cell death in vitro and in vivo. Thus, retinoblastomas may develop quickly as a result of the epigenetic deregulation of key cancer pathways as a direct or indirect result of RB1 loss

Cellular Metabolomics Profiles Associated With Drug Chemosensitivity in AML

BackgroundAcute myeloid leukemia (AML) is a hematological malignancy with a dismal prognosis. For over four decades, AML has primarily been treated by cytarabine combined with an anthracycline. Although a significant proportion of patients achieve remission with this regimen, roughly 40% of children and 70% of adults relapse. Over 90% of patients with resistant or relapsed AML die within 3 years. Thus, relapsed and resistant disease following treatment with standard therapy are the most common clinical failures that occur in treating this disease. In this study, we evaluated the relationship between AML cell line global metabolomes and variation in chemosensitivity.MethodsWe performed global metabolomics on seven AML cell lines with varying chemosensitivity to cytarabine and the anthracycline doxorubicin (MV4.11, KG-1, HL-60, Kasumi-1, AML-193, ME1, THP-1) using ultra-high performance liquid chromatography – mass spectrometry (UHPLC-MS). Univariate and multivariate analyses were performed on the metabolite peak intensity values from UHPLC-MS using MetaboAnalyst to identify cellular metabolites associated with drug chemosensitivity.ResultsA total of 1,624 metabolic features were detected across the leukemic cell lines. Of these, 187 were annotated to known metabolites. With respect to doxorubicin, we observed significantly greater abundance of a carboxylic acid (1-aminocyclopropane-1-carboxylate) and several amino acids in resistant cell lines. Pathway analysis found enrichment of several amino acid biosynthesis and metabolic pathways. For cytarabine resistance, nine annotated metabolites were significantly different in resistance vs. sensitive cell lines, including D-raffinose, guanosine, inosine, guanine, aldopentose, two xenobiotics (allopurinol and 4-hydroxy-L-phenylglycine) and glucosamine/mannosamine. Pathway analysis associated these metabolites with the purine metabolic pathway.ConclusionOverall, our results demonstrate that metabolomics differences contribute toward drug resistance. In addition, it could potentially identify predictive biomarkers for chemosensitivity to various anti-leukemic drugs. Our results provide opportunity to further explore these metabolites in patient samples for association with clinical response

ChIP-PaM: an algorithm to identify protein-DNA interaction using ChIP-Seq data

<p>Abstract</p> <p>Background</p> <p>ChIP-Seq is a powerful tool for identifying the interaction between genomic regulators and their bound DNAs, especially for locating transcription factor binding sites. However, high cost and high rate of false discovery of transcription factor binding sites identified from ChIP-Seq data significantly limit its application.</p> <p>Results</p> <p>Here we report a new algorithm, ChIP-PaM, for identifying transcription factor target regions in ChIP-Seq datasets. This algorithm makes full use of a protein-DNA binding pattern by capitalizing on three lines of evidence: 1) the tag count modelling at the peak position, 2) pattern matching of a specific tag count distribution, and 3) motif searching along the genome. A novel data-based two-step eFDR procedure is proposed to integrate the three lines of evidence to determine significantly enriched regions. Our algorithm requires no technical controls and efficiently discriminates falsely enriched regions from regions enriched by true transcription factor (TF) binding on the basis of ChIP-Seq data only. An analysis of real genomic data is presented to demonstrate our method.</p> <p>Conclusions</p> <p>In a comparison with other existing methods, we found that our algorithm provides more accurate binding site discovery while maintaining comparable statistical power.</p

SVSI: Fast and Powerful Set-Valued System Identification Approach to Identifying Rare Variants in Sequencing Studies for Ordered Categorical Traits: SVSIfor Genetic Association Studies

For genetic association studies that involve an ordered categorical phenotype, we usually either regroup multiple categories of the phenotype into two categories (“cases” and “controls”) and then apply the standard logistic regression (LG), or apply ordered logistic (oLG) or ordered probit (oPRB) regression which accounts for the ordinal nature of the phenotype. However, these approaches may lose statistical power or may not control type I error rate due to their model assumption and/or instable parameter estimation algorithm when the genetic variant is rare or sample size is limited. Here to solve this problem, we propose a set-valued (SV) system model, which assumes that an underlying continuous phenotype follows a normal distribution, to identify genetic variants associated with an ordinal categorical phenotype. We couple this model with a set-valued system identification algorithm to identify all the key system parameters. Simulations and two real data analyses show that SV and LG accurately controlled the Type I error rate even at a significance level of 10−6 but not oLG and oPRB in some cases. LG had significantly smaller power than the other three methods due to disregarding of the ordinal nature of the phenotype, and SV had similar or greater power than oLG and oPRB. For instance, in a simulation with data generated from an additive SV model with odds ratio of 7.4 for a phenotype with three categories, a single nucleotide polymorphism with minor allele frequency of 0.75% and sample size of 999 (333 per category), the power of SV, oLG and LG models were 70%, 40% and <1%, respectively, at a significance level of 10−6. Thus, SV should be employed in genetic association studies for ordered categorical phenotype

Pratos e mais pratos: louças domésticas, divisões culturais e limites sociais no Rio de Janeiro, século XIX

Reply to ten comments on a paper published in the last issue of this journal. The discussion follows along six main lines: History museums, identity, ideology and the category of nation; the need of material collections and their modalities: patrimonial, operational, virtual; theater versus laboratory; visitors and their ambiguities; Public History: the museum and the academy.Resposta aos comentários de dez especialistas que contribuíram no debate de texto publicado no último número desta revista. A discussão orientou-se segundo seis tópicos principais: museus históricos, identidade, ideologia e a categoria de nação; a necessidade de acervos materiais e suas modalidades: acervo patrimonial, operacional, virtual; teatro versus laboratório; o público e suas ambigüidades; História Pública: o museu e a Academia