A Common Polymorphism in HIBCH Influences Methylmalonic Acid Concentrations in Blood Independently of Cobalamin.

Methylmalonic acid (MMA) is a by-product of propionic acid metabolism through the vitamin B12 (cobalamin)-dependent enzyme methylmalonyl CoA mutase. Elevated MMA concentrations are a hallmark of several inborn errors of metabolism and indicators of cobalamin deficiency in older persons. In a genome-wide analysis of 2,210 healthy young Irish adults (median age 22 years) we identified a strong association of plasma MMA with SNPs in 3-hydroxyisobutyryl-CoA hydrolase (HIBCH, p = 8.42 × 10(-89)) and acyl-CoA synthetase family member 3 (ACSF3, p = 3.48 × 10(-19)). These loci accounted for 12% of the variance in MMA concentration. The most strongly associated SNP (HIBCH rs291466; c:2T>C) causes a missense change of the initiator methionine codon (minor-allele frequency = 0.43) to threonine. Surprisingly, the resulting variant, p.Met1?, is associated with increased expression of HIBCH mRNA and encoded protein. These homozygotes had, on average, 46% higher MMA concentrations than methionine-encoding homozygotes in young adults with generally low MMA concentrations (0.17 [0.14-0.21] µmol/L; median [25(th)-75(th) quartile]). The association between MMA levels and HIBCH rs291466 was highly significant in a replication cohort of 1,481 older individuals (median age 79 years) with elevated plasma MMA concentrations (0.34 [0.24-0.51] µmol/L; p = 4.0 × 10(-26)). In a longitudinal study of 185 pregnant women and their newborns, the association of this SNP remained significant across the gestational trimesters and in newborns. HIBCH is unique to valine catabolism. Studies evaluating flux through the valine catabolic pathway in humans should account for these variants. Furthermore, this SNP could help resolve equivocal clinical tests where plasma MMA values have been used to diagnose cobalamin deficiency.

Accurate and comprehensive sequencing of personal genomes.

As whole-genome sequencing becomes commoditized and we begin to sequence and analyze personal genomes for clinical and diagnostic purposes, it is necessary to understand what constitutes a complete sequencing experiment for determining genotypes and detecting single-nucleotide variants. Here, we show that the current recommendation of ∼30× coverage is not adequate to produce genotype calls across a large fraction of the genome with acceptably low error rates. Our results are based on analyses of a clinical sample sequenced on two related Illumina platforms, GAII(x) and HiSeq 2000, to a very high depth (126×). We used these data to establish genotype-calling filters that dramatically increase accuracy. We also empirically determined how the callable portion of the genome varies as a function of the amount of sequence data used. These results help provide a "sequencing guide" for future whole-genome sequencing decisions and metrics by which coverage statistics should be reported.

Replication and exploratory analysis of 24 candidate risk polymorphisms for neural tube defects.

BACKGROUND: Neural tube defects (NTDs), which are among the most common congenital malformations, are influenced by environmental and genetic factors. Low maternal folate is the strongest known contributing factor, making variants in genes in the folate metabolic pathway attractive candidates for NTD risk. Multiple studies have identified nominally significant allelic associations with NTDs. We tested whether associations detected in a large Irish cohort could be replicated in an independent population. METHODS: Replication tests of 24 nominally significant NTD associations were performed in racially/ethnically matched populations. Family-based tests of fifteen nominally significant single nucleotide polymorphisms (SNPs) were repeated in a cohort of NTD trios (530 cases and their parents) from the United Kingdom, and case-control tests of nine nominally significant SNPs were repeated in a cohort (190 cases, 941 controls) from New York State (NYS). Secondary hypotheses involved evaluating the latter set of nine SNPs for NTD association using alternate case-control models and NTD groupings in white, African American and Hispanic cohorts from NYS. RESULTS: Of the 24 SNPs tested for replication, ADA rs452159 and MTR rs10925260 were significantly associated with isolated NTDs. Of the secondary tests performed, ARID1A rs11247593 was associated with NTDs in whites, and ALDH1A2 rs7169289 was associated with isolated NTDs in African Americans. CONCLUSIONS: We report a number of associations between SNP genotypes and neural tube defects. These associations were nominally significant before correction for multiple hypothesis testing. These corrections are highly conservative for association studies of untested hypotheses, and may be too conservative for replication studies. We therefore believe the true effect of these four nominally significant SNPs on NTD risk will be more definitively determined by further study in other populations, and eventual meta-analysis.

A new strategy for genome assembly using short sequence reads and reduced representation libraries.

We have developed a novel approach for using massively parallel short-read sequencing to generate fast and inexpensive de novo genomic assemblies comparable to those generated by capillary-based methods. The ultrashort (<100 base) sequences generated by this technology pose specific biological and computational challenges for de novo assembly of large genomes. To account for this, we devised a method for experimentally partitioning the genome using reduced representation (RR) libraries prior to assembly. We use two restriction enzymes independently to create a series of overlapping fragment libraries, each containing a tractable subset of the genome. Together, these libraries allow us to reassemble the entire genome without the need of a reference sequence. As proof of concept, we applied this approach to sequence and assembled the majority of the 125-Mb Drosophila melanogaster genome. We subsequently demonstrate the accuracy of our assembly method with meaningful comparisons against the current available D. melanogaster reference genome (dm3). The ease of assembly and accuracy for comparative genomics suggest that our approach will scale to future mammalian genome-sequencing efforts, saving both time and money without sacrificing quality.

Systematic comparison of three genomic enrichment methods for massively parallel DNA sequencing.

Massively parallel DNA sequencing technologies have greatly increased our ability to generate large amounts of sequencing data at a rapid pace. Several methods have been developed to enrich for genomic regions of interest for targeted sequencing. We have compared three of these methods: Molecular Inversion Probes (MIP), Solution Hybrid Selection (SHS), and Microarray-based Genomic Selection (MGS). Using HapMap DNA samples, we compared each of these methods with respect to their ability to capture an identical set of exons and evolutionarily conserved regions associated with 528 genes (2.61 Mb). For sequence analysis, we developed and used a novel Bayesian genotype-assigning algorithm, Most Probable Genotype (MPG). All three capture methods were effective, but sensitivities (percentage of targeted bases associated with high-quality genotypes) varied for an equivalent amount of pass-filtered sequence: for example, 70% (MIP), 84% (SHS), and 91% (MGS) for 400 Mb. In contrast, all methods yielded similar accuracies of >99.84% when compared to Infinium 1M SNP BeadChip-derived genotypes and >99.998% when compared to 30-fold coverage whole-genome shotgun sequencing data. We also observed a low false-positive rate with all three methods; of the heterozygous positions identified by each of the capture methods, >99.57% agreed with 1M SNP BeadChip, and >98.840% agreed with the whole-genome shotgun data. In addition, we successfully piloted the genomic enrichment of a set of 12 pooled samples via the MGS method using molecular bar codes. We find that these three genomic enrichment methods are highly accurate and practical, with sensitivities comparable to that of 30-fold coverage whole-genome shotgun data.

Wnt10b induces chemotaxis of osteosarcoma and correlates with reduced survival.

BACKGROUND: Osteosarcoma (OS) is a primary malignant tumor of the bone that typically presents in the second decade of life and has a poor prognosis, especially in metastatic cases. Wnt signaling contributes to the pathogenesis of tumors such as colon cancer and malignant melanoma. Wnt signaling controls normal bone formation during embryogenesis and homeostasis in adult organisms, thus we evaluated Wnt signaling in OS. PROCEDURE: We surveyed the expression of Wnts, their receptors, Frizzleds and LRPs, and soluble Wnt inhibitors (sFRPs) in four OS cell lines by RT-PCR. We also tested biological response of OS cell lines to exogenous Wnts by measuring beta-catenin stabilization, Dvl phosphorylation, TOPFLASH activity and chemotaxis. Human OS tumor microarrays were evaluated for expression of Wnt10b by immunohistochemistry. RESULTS: All cell lines tested showed expression of at least three Wnts and one Frizzled. Exogenous Wnt3a and Wnt10b treatment induced Dvl phosphorylation, beta-catenin stabilization and TCF4 transcriptional activity in both metastatic and non-metastatic murine OS cell lines. Metastatic OS cell lines showed better chemotaxis response to Wnts than the non-metastatic OS cell lines. Immunohistochemistry studies of 44 human OS samples demonstrated that Wnt10b expression correlated with decreased overall survival. CONCLUSIONS: These results further supports a possible autocrine or paracrine Wnt pathway in metastatic potential of OS.

Structure-based design of potent small-molecule inhibitors of anti-apoptotic Bcl-2 proteins.

A structure-based approach was employed to design a new class of small-molecule inhibitors of Bcl-2. The most potent compound 5 (TW-37) binds to Bcl-2 with a K(i) value of 290 nM and also to Bcl-xL and Mcl-1 with high affinities. Compound 5 potently inhibits cell growth in PC-3 prostate cancer cells with an IC(50) value of 200 nM and effectively induces apoptosis in a dose-dependent manner.

Genomics of CpG methylation in developing and developed zebrafish.

DNA methylation is a dynamic process through which specific chromatin modifications can be stably transmitted from parent to daughter cells. A large body of work has suggested that DNA methylation influences gene expression by silencing gene promoters. However, these conclusions were drawn from data focused mostly on promoter regions. Regarding the entire genome, it is unclear how methylation and gene transcription patterns are related during vertebrate development. To identify the genome-wide distribution of CpG methylation, we created series of high-resolution methylome maps of Danio rerio embryos during development and in mature, differentiated tissues. We found that embryonic and terminal tissues have unique methylation signatures in CpG islands and repetitive sequences. Fully differentiated tissues have increased CpG and LTR methylation and decreased SINE methylation relative to embryonic tissues. Unsupervised clustering analyses reveal that the embryonic and terminal tissues can be classified solely by their methylation patterning. Novel analyses also identify a previously undescribed genome-wide exon methylation signature. We also compared whole genome methylation with genome-wide mRNA expression levels using publicly available RNA-seq datasets. These comparisons revealed previously unrecognized relationships between gene expression, alternative splicing, and exon methylation. Surprisingly, we found that exonic methylation is a better predictor of mRNA expression level than promoter methylation. We also found that transcriptionally skipped exons have significantly less methylation than retained exons. Our integrative analyses reveal highly complex interplay between gene expression, alternative splicing, development, and methylation patterning in zebrafish.

A transcriptomic atlas of mouse neocortical layers.

In the mammalian cortex, neurons and glia form a patterned structure across six layers whose complex cytoarchitectonic arrangement is likely to contribute to cognition. We sequenced transcriptomes from layers 1-6b of different areas (primary and secondary) of the adult (postnatal day 56) mouse somatosensory cortex to understand the transcriptional levels and functional repertoires of coding and noncoding loci for cells constituting these layers. A total of 5,835 protein-coding genes and 66 noncoding RNA loci are differentially expressed ("patterned") across the layers, on the basis of a machine-learning model (naive Bayes) approach. Layers 2-6b are each associated with specific functional and disease annotations that provide insights into their biological roles. This new resource (http://genserv.anat.ox.ac.uk/layers) greatly extends currently available resources, such as the Allen Mouse Brain Atlas and microarray data sets, by providing quantitative expression levels, by being genome-wide, by including novel loci, and by identifying candidate alternatively spliced transcripts that are differentially expressed across layers.

Local DNA topography correlates with functional noncoding regions of the human genome.

The three-dimensional molecular structure of DNA, specifically the shape of the backbone and grooves of genomic DNA, can be dramatically affected by nucleotide changes, which can cause differences in protein-binding affinity and phenotype. We developed an algorithm to measure constraint on the basis of similarity of DNA topography among multiple species, using hydroxyl radical cleavage patterns to interrogate the solvent-accessible surface area of DNA. This algorithm found that 12% of bases in the human genome are evolutionarily constrained-double the number detected by nucleotide sequence-based algorithms. Topography-informed constrained regions correlated with functional noncoding elements, including enhancers, better than did regions identified solely on the basis of nucleotide sequence. These results support the idea that the molecular shape of DNA is under selection and can identify evolutionary history.