Brian D. O'Connor

U87MG Decoded: The Genomic Sequence of a Cytogenetically Aberrant Human Cancer Cell Line

U87MG is a commonly studied grade IV glioma cell line that has been analyzed in at least 1,700 publications over four decades. In order to comprehensively characterize the genome of this cell line and to serve as a model of broad cancer genome sequencing, we have generated greater than 30x genomic sequence coverage using a novel 50-base mate paired strategy with a 1.4kb mean insert library. A total of 1,014,984,286 mate-end and 120,691,623 single-end two-base encoded reads were generated from five slides. All data were aligned using a custom designed tool called BFAST, allowing optimal color space read alignment and accurate identification of DNA variants. The aligned sequence reads and mate-pair information identified 35 interchromosomal translocation events, 1,315 structural variations (>100 bp), 191,743 small (<21 bp) insertions and deletions (indels), and 2,384,470 single nucleotide variations (SNVs). Among these observations, the known homozygous mutation in PTEN was robustly identified, and genes involved in cell adhesion were overrepresented in the mutated gene list. Data were compared to 219,187 heterozygous single nucleotide polymorphisms assayed by Illumina 1M Duo genotyping array to assess accuracy: 93.83% of all SNPs were reliably detected at filtering thresholds that yield greater than 99.99% sequence accuracy. Protein coding sequences were disrupted predominantly in this cancer cell line due to small indels, large deletions, and translocations. In total, 512 genes were homozygously mutated, including 154 by SNVs, 178 by small indels, 145 by large microdeletions, and 35 by interchromosomal translocations to reveal a highly mutated cell line genome. Of the small homozygously mutated variants, 8 SNVs and 99 indels were novel events not present in dbSNP. These data demonstrate that routine generation of broad cancer genome sequence is possible outside of genome centers. The sequence analysis of U87MG provides an unparalleled level of mutational resolution compared to any cell line to date.

Mutations in PYCR1 cause cutis laxa with progeroid features.

Autosomal recessive cutis laxa (ARCL) describes a group of syndromal disorders that are often associated with a progeroid appearance, lax and wrinkled skin, osteopenia and mental retardation. Homozygosity mapping in several kindreds with ARCL identified a candidate region on chromosome 17q25. By high-throughput sequencing of the entire candidate region, we detected disease-causing mutations in the gene PYCR1. We found that the gene product, an enzyme involved in proline metabolism, localizes to mitochondria. Altered mitochondrial morphology, membrane potential and increased apoptosis rate upon oxidative stress were evident in fibroblasts from affected individuals. Knockdown of the orthologous genes in Xenopus and zebrafish led to epidermal hypoplasia and blistering that was accompanied by a massive increase of apoptosis. Our findings link mutations in PYCR1 to altered mitochondrial function and progeroid changes in connective tissues.

The Genomics of Disulfide Bonding and Protein Stabilization in Thermophiles

Thermophilic organisms flourish in varied high-temperature environmental niches that are deadly to other organisms. Recently, genomic evidence has implicated a critical role for disulfide bonds in the structural stabilization of intracellular proteins from certain of these organisms, contrary to the conventional view that structural disulfide bonds are exclusively extracellular. Here both computational and structural data are presented to explore the occurrence of disulfide bonds as a protein-stabilization method across many thermophilic prokaryotes. Based on computational studies, disulfide-bond richness is found to be widespread, with thermophiles containing the highest levels. Interestingly, only a distinct subset of thermophiles exhibit this property. A computational search for proteins matching this target phylogenetic profile singles out a specific protein, known as protein disulfide oxidoreductase, as a potential key player in thermophilic intracellular disulfide-bond formation. Finally, biochemical support in the form of a new crystal structure of a thermophilic protein with three disulfide bonds is presented together with a survey of known structures from the literature. Together, the results provide insight into biochemical specialization and the diversity of methods employed by organisms to stabilize their proteins in exotic environments. The findings also motivate continued efforts to sequence genomes from divergent organisms.

Multiple self-healing squamous epithelioma is caused by a disease-specific spectrum of mutations in TGFBR1

Multiple self-healing squamous epithelioma (MSSE), also known as Ferguson-Smith disease (FSD), is an autosomal-dominant skin cancer condition characterized by multiple squamous-carcinoma–like locally invasive skin tumors that grow rapidly for a few weeks before spontaneously regressing, leaving scars. High-throughput genomic sequencing of a conservative estimate (24.2 Mb) of the disease locus on chromosome 9 using exon array capture identified independent mutations in TGFBR1 in three unrelated families. Subsequent dideoxy sequencing of TGFBR1 identified 11 distinct monoallelic mutations in 18 affected families, firmly establishing TGFBR1 as the causative gene. The nature of the sequence variants, which include mutations in the extracellular ligand-binding domain and a series of truncating mutations in the kinase domain, indicates a clear genotype-phenotype correlation between loss-of-function TGFBR1 mutations and MSSE. This distinguishes MSSE from the Marfan syndrome–related disorders in which missense mutations in TGFBR1 lead to developmental defects with vascular involvement but no reported predisposition to cancer.

Celsius: a community resource for Affymetrix microarray data

Celsius is a data warehousing system to aggregate Affymetrix CEL files and associated metadata. It provides mechanisms for importing, storing, querying, and exporting large volumes of primary and pre-processed microarray data. Celsius contains ten billion assay measurements and affiliated metadata. It is the largest publicly available source of Affymetrix microarray data, and through sheer volume it allows a sophisticated, broad view of transcription that has not previously been possible.

Pathogenicity of a disease-associated human IL-4 receptor allele in experimental asthma

Polymorphisms in the interleukin-4 receptor α chain (IL-4Rα) have been linked to asthma incidence and severity, but a causal relationship has remained uncertain. In particular, a glutamine to arginine substitution at position 576 (Q576R) of IL-4Rα has been associated with severe asthma, especially in African Americans. We show that mice carrying the Q576R polymorphism exhibited intense allergen-induced airway inflammation and remodeling. The Q576R polymorphism did not affect proximal signal transducer and activator of transcription (STAT) 6 activation, but synergized with STAT6 in a gene target– and tissue-specific manner to mediate heightened expression of a subset of IL-4– and IL-13–responsive genes involved in allergic inflammation. Our findings indicate that the Q576R polymorphism directly promotes asthma in carrier populations by selectively augmenting IL-4Rα–dependent signaling.

GMODWeb: a web framework for the generic model organism database

The Generic Model Organism Database (GMOD) initiative provides species-agnostic data models and software tools for representing curated model organism data. Here we describe GMODWeb, a GMOD project designed to speed the development of model organism database (MOD) websites. Sites created with GMODWeb provide integration with other GMOD tools and allow users to browse and search through a variety of data types. GMODWeb was built using the open source Turnkey web framework and is available from http://turnkey.sourceforge.net.

Utilizing logical relationships in genomic data to decipher cellular processes

The wealth of available genomic data has spawned a corresponding interest in computational methods that can impart biological meaning and context to these experiments. Traditional computational methods have drawn relationships between pairs of proteins or genes based on notions of equality or similarity between their patterns of occurrence or behavior. For example, two genes displaying similar variation in expression, over a number of experiments, may be predicted to be functionally related. We have introduced a natural extension of these approaches, instead identifying logical relationships involving triplets of proteins. Triplets provide for various discrete kinds of logic relationships, leading to detailed inferences about biological associations. For instance, a protein C might be encoded within an organism if, and only if, two other proteins A and B are also both encoded within the organism, thus suggesting that gene C is functionally related to genes A and B. The method has been applied fruitfully to both phylogenetic and microarray expression data, and has been used to associate logical combinations of protein activity with disease state phenotypes, revealing previously unknown ternary relationships among proteins, and illustrating the inherent complexities that arise in biological data.

GDAP: a web tool for genome-wide protein disulfide bond prediction

The Genomic Disulfide Analysis Program (GDAP) provides web access to computationally predicted protein disulfide bonds for over one hundred microbial genomes, including both bacterial and achaeal species. In the GDAP process, sequences of unknown structure are mapped, when possible, to known homologous Protein Data Bank (PDB) structures, after which specific distance criteria are applied to predict disulfide bonds. GDAP also accepts user-supplied protein sequences and subsequently queries the PDB sequence database for the best matches, scans for possible disulfide bonds and returns the results to the client. These predictions are useful for a variety of applications and have previously been used to show a dramatic preference in certain thermophilic archaea and bacteria for disulfide bonds within intracellular proteins. Given the central role these stabilizing, covalent bonds play in such organisms, the predictions available from GDAP provide a rich data source for designing site-directed mutants with more stable thermal profiles. The GDAP web application is a gateway to this information and can be used to understand the role disulfide bonds play in protein stability both in these unusual organisms and in sequences of interest to the individual researcher. The prediction server can be accessed at http://www.doe-mbi.ucla.edu/Services/GDAP.

Correction: U87MG Decoded: The Genomic Sequence of a Cytogenetically Aberrant Human Cancer Cell Line

[This corrects the article DOI: 10.1371/journal.pgen.1000832.].