Reconstructing cancer genomes from paired-end sequencing data

<p>Abstract</p> <p>Background</p> <p>A cancer genome is derived from the germline genome through a series of somatic mutations. Somatic structural variants - including duplications, deletions, inversions, translocations, and other rearrangements - result in a cancer genome that is a scrambling of intervals, or "blocks" of the germline genome sequence. We present an efficient algorithm for reconstructing the block organization of a cancer genome from paired-end DNA sequencing data.</p> <p>Results</p> <p>By aligning paired reads from a cancer genome - and a matched germline genome, if available - to the human reference genome, we derive: (i) a partition of the reference genome into intervals; (ii) adjacencies between these intervals in the cancer genome; (iii) an estimated copy number for each interval. We formulate the Copy Number and Adjacency Genome Reconstruction Problem of determining the cancer genome as a sequence of the derived intervals that is consistent with the measured adjacencies and copy numbers. We design an efficient algorithm, called Paired-end Reconstruction of Genome Organization (PREGO), to solve this problem by reducing it to an optimization problem on an interval-adjacency graph constructed from the data. The solution to the optimization problem results in an Eulerian graph, containing an alternating Eulerian tour that corresponds to a cancer genome that is consistent with the sequencing data. We apply our algorithm to five ovarian cancer genomes that were sequenced as part of The Cancer Genome Atlas. We identify numerous rearrangements, or structural variants, in these genomes, analyze reciprocal vs. non-reciprocal rearrangements, and identify rearrangements consistent with known mechanisms of duplication such as tandem duplications and breakage/fusion/bridge (B/F/B) cycles.</p> <p>Conclusions</p> <p>We demonstrate that PREGO efficiently identifies complex and biologically relevant rearrangements in cancer genome sequencing data. An implementation of the PREGO algorithm is available at <url>http://compbio.cs.brown.edu/software/</url>.</p

Nasal Bone Shape Is under Complex Epistatic Genetic Control in Mouse Interspecific Recombinant Congenic Strains

Genetic determinism of cranial morphology in the mouse is still largely unknown, despite the localization of putative QTLs and the identification of genes associated with Mendelian skull malformations. To approach the dissection of this multigenic control, we have used a set of interspecific recombinant congenic strains (IRCS) produced between C57BL/6 and mice of the distant species Mus spretus (SEG/Pas). Each strain has inherited 1.3% of its genome from SEG/Pas under the form of few, small-sized, chromosomal segments.The shape of the nasal bone was studied using outline analysis combined with Fourier descriptors, and differential features were identified between IRCS BcG-66H and C57BL/6. An F2 cross between BcG-66H and C57BL/6 revealed that, out of the three SEG/Pas-derived chromosomal regions present in BcG-66H, two were involved. Segments on chromosomes 1 (∼32 Mb) and 18 (∼13 Mb) showed additive effect on nasal bone shape. The three chromosomal regions present in BcG-66H were isolated in congenic strains to study their individual effect. Epistatic interactions were assessed in bicongenic strains.Our results show that, besides a strong individual effect, the QTL on chromosome 1 interacts with genes on chromosomes 13 and 18. This study demonstrates that nasal bone shape is under complex genetic control but can be efficiently dissected in the mouse using appropriate genetic tools and shape descriptors

Gene Dosage, Expression, and Ontology Analysis Identifies Driver Genes in the Carcinogenesis and Chemoradioresistance of Cervical Cancer

Integrative analysis of gene dosage, expression, and ontology (GO) data was performed to discover driver genes in the carcinogenesis and chemoradioresistance of cervical cancers. Gene dosage and expression profiles of 102 locally advanced cervical cancers were generated by microarray techniques. Fifty-two of these patients were also analyzed with the Illumina expression method to confirm the gene expression results. An independent cohort of 41 patients was used for validation of gene expressions associated with clinical outcome. Statistical analysis identified 29 recurrent gains and losses and 3 losses (on 3p, 13q, 21q) associated with poor outcome after chemoradiotherapy. The intratumor heterogeneity, assessed from the gene dosage profiles, was low for these alterations, showing that they had emerged prior to many other alterations and probably were early events in carcinogenesis. Integration of the alterations with gene expression and GO data identified genes that were regulated by the alterations and revealed five biological processes that were significantly overrepresented among the affected genes: apoptosis, metabolism, macromolecule localization, translation, and transcription. Four genes on 3p (RYBP, GBE1) and 13q (FAM48A, MED4) correlated with outcome at both the gene dosage and expression level and were satisfactorily validated in the independent cohort. These integrated analyses yielded 57 candidate drivers of 24 genetic events, including novel loci responsible for chemoradioresistance. Further mapping of the connections among genetic events, drivers, and biological processes suggested that each individual event stimulates specific processes in carcinogenesis through the coordinated control of multiple genes. The present results may provide novel therapeutic opportunities of both early and advanced stage cervical cancers

Human papillomavirus-mediated carcinogenesis and HPV-associated oral and oropharyngeal squamous cell carcinoma. Part 1: Human papillomavirus-mediated carcinogenesis

High-risk human papillomavirus (HPV) E6 and E7 oncoproteins are essential factors for HPV-induced carcinogenesis, and for the maintenance of the consequent neoplastic growth. Cellular transformation is achieved by complex interaction of these oncogenes with several cellular factors of cell cycle regulation including p53, Rb, cyclin-CDK complexes, p21 and p27. Both persistent infection with high-risk HPV genotypes and immune dysregulation are associated with increased risk of HPV-induced squamous cell carcinoma

Gastrointestinal decontamination in the acutely poisoned patient

ObjectiveTo define the role of gastrointestinal (GI) decontamination of the poisoned patient.Data sourcesA computer-based PubMed/MEDLINE search of the literature on GI decontamination in the poisoned patient with cross referencing of sources.Study selection and data extractionClinical, animal and in vitro studies were reviewed for clinical relevance to GI decontamination of the poisoned patient.Data synthesisThe literature suggests that previously, widely used, aggressive approaches including the use of ipecac syrup, gastric lavage, and cathartics are now rarely recommended. Whole bowel irrigation is still often recommended for slow-release drugs, metals, and patients who "pack" or "stuff" foreign bodies filled with drugs of abuse, but with little quality data to support it. Activated charcoal (AC), single or multiple doses, was also a previous mainstay of GI decontamination, but the utility of AC is now recognized to be limited and more time dependent than previously practiced. These recommendations have resulted in several treatment guidelines that are mostly based on retrospective analysis, animal studies or small case series, and rarely based on randomized clinical trials.ConclusionsThe current literature supports limited use of GI decontamination of the poisoned patient

Proteomic Changes Resulting from Gene Copy Number Variations in Cancer Cells

Along the transformation process, cells accumulate DNA aberrations, including mutations, translocations, amplifications, and deletions. Despite numerous studies, the overall effects of amplifications and deletions on the end point of gene expression—the level of proteins—is generally unknown. Here we use large-scale and high-resolution proteomics combined with gene copy number analysis to investigate in a global manner to what extent these genomic changes have a proteomic output and therefore the ability to affect cellular transformation. We accurately measure expression levels of 6,735 proteins and directly compare them to the gene copy number. We find that the average effect of these alterations on the protein expression is only a few percent. Nevertheless, by using a novel algorithm, we find the combined impact that many of these regional chromosomal aberrations have at the protein level. We show that proteins encoded by amplified oncogenes are often overexpressed, while adjacent amplified genes, which presumably do not promote growth and survival, are attenuated. Furthermore, regulation of biological processes and molecular complexes is independent of general copy number changes. By connecting the primary genome alteration to their proteomic consequences, this approach helps to interpret the data from large-scale cancer genomics efforts

Signatures of Selection in Fusion Transcripts Resulting From Chromosomal Translocations in Human Cancer

BACKGROUND: The recurrence and non-random distribution of translocation breakpoints in human tumors are usually attributed to local sequence features present in the vicinity of the breakpoints. However, it has also been suggested that functional constraints might contribute to delimit the position of translocation breakpoints within the genes involved, but a quantitative analysis of such contribution has been lacking. METHODOLOGY: We have analyzed two well-known signatures of functional selection, such as reading-frame compatibility and non-random combinations of protein domains, on an extensive dataset of fusion proteins resulting from chromosomal translocations in cancer. CONCLUSIONS: Our data provide strong experimental support for the concept that the position of translocation breakpoints in the genome of cancer cells is determined, to a large extent, by the need to combine certain protein domains and to keep an intact reading frame in fusion transcripts. Additionally, the information that we have assembled affords a global view of the oncogenic mechanisms and domain architectures that are used by fusion proteins. This can be used to assess the functional impact of novel chromosomal translocations and to predict the position of breakpoints in the genes involved

Functional Copy-Number Alterations in Cancer

Understanding the molecular basis of cancer requires characterization of its genetic defects. DNA microarray technologies can provide detailed raw data about chromosomal aberrations in tumor samples. Computational analysis is needed (1) to deduce from raw array data actual amplification or deletion events for chromosomal fragments and (2) to distinguish causal chromosomal alterations from functionally neutral ones. We present a comprehensive computational approach, RAE, designed to robustly map chromosomal alterations in tumor samples and assess their functional importance in cancer. To demonstrate the methodology, we experimentally profile copy number changes in a clinically aggressive subtype of soft-tissue sarcoma, pleomorphic liposarcoma, and computationally derive a portrait of candidate oncogenic alterations and their target genes. Many affected genes are known to be involved in sarcomagenesis; others are novel, including mediators of adipocyte differentiation, and may include valuable therapeutic targets. Taken together, we present a statistically robust methodology applicable to high-resolution genomic data to assess the extent and function of copy-number alterations in cancer

The HLH-6 Transcription Factor Regulates C. elegans Pharyngeal Gland Development and Function

The Caenorhabditis elegans pharynx (or foregut) functions as a pump that draws in food (bacteria) from the environment. While the “organ identity factor” PHA-4 is critical for formation of the C. elegans pharynx as a whole, little is known about the specification of distinct cell types within the pharynx. Here, we use a combination of bioinformatics, molecular biology, and genetics to identify a helix-loop-helix transcription factor (HLH-6) as a critical regulator of pharyngeal gland development. HLH-6 is required for expression of a number of gland-specific genes, acting through a discrete cis-regulatory element named PGM1 (Pharyngeal Gland Motif 1). hlh-6 mutants exhibit a frequent loss of a subset of glands, while the remaining glands have impaired activity, indicating a role for hlh-6 in both gland development and function. Interestingly, hlh-6 mutants are also feeding defective, ascribing a biological function for the glands. Pharyngeal pumping in hlh-6 mutants is normal, but hlh-6 mutants lack expression of a class of mucin-related proteins that are normally secreted by pharyngeal glands and line the pharyngeal cuticle. An interesting possibility is that one function of pharyngeal glands is to secrete a pharyngeal lining that ensures efficient transport of food along the pharyngeal lumen

Gene Expression and Functional Studies of the Optic Nerve Head Astrocyte Transcriptome from Normal African Americans and Caucasian Americans Donors

To determine whether optic nerve head (ONH) astrocytes, a key cellular component of glaucomatous neuropathy, exhibit differential gene expression in primary cultures of astrocytes from normal African American (AA) donors compared to astrocytes from normal Caucasian American (CA) donors.We used oligonucleotide Affymetrix microarray (HG U133A & HG U133A 2.0 chips) to compare gene expression levels in cultured ONH astrocytes from twelve CA and twelve AA normal age matched donor eyes. Chips were normalized with Robust Microarray Analysis (RMA) in R using Bioconductor. Significant differential gene expression levels were detected using mixed effects modeling and Statistical Analysis of Microarray (SAM). Functional analysis and Gene Ontology were used to classify differentially expressed genes. Differential gene expression was validated by quantitative real time RT-PCR. Protein levels were detected by Western blots and ELISA. Cell adhesion and migration assays tested physiological responses. Glutathione (GSH) assay detected levels of intracellular GSH.Multiple analyses selected 87 genes differentially expressed between normal AA and CA (P<0.01). The most relevant genes expressed in AA were categorized by function, including: signal transduction, response to stress, ECM genes, migration and cell adhesion.These data show that normal astrocytes from AA and CA normal donors display distinct expression profiles that impact astrocyte functions in the ONH. Our data suggests that differences in gene expression in ONH astrocytes may be specific to the development and/or progression of glaucoma in AA