Explore the vast range of titles available.

The new generation of massively parallel DNA sequencers, combined with the challenge of whole human genome resequencing, result in the need for rapid and accurate alignment of billions of short DNA sequence reads to a large reference genome. Speed is obviously of great importance, but equally important is maintaining alignment accuracy of short reads, in the 25-100 base range, in the presence of errors and true biological variation. We introduce a new algorithm specifically optimized for this task, as well as a freely available implementation, BFAST, which can align data produced by any of current sequencing platforms, allows for user-customizable levels of speed and accuracy, supports paired end data, and provides for efficient parallel and multi-threaded computation on a computer cluster. The new method is based on creating flexible, efficient whole genome indexes to rapidly map reads to candidate alignment locations, with arbitrary multiple independent indexes allowed to achieve robustness against read errors and sequence variants. The final local alignment uses a Smith-Waterman method, with gaps to support the detection of small indels. We compare BFAST to a selection of large-scale alignment tools -- BLAT, MAQ, SHRiMP, and SOAP -- in terms of both speed and accuracy, using simulated and real-world datasets. We show BFAST can achieve substantially greater sensitivity of alignment in the context of errors and true variants, especially insertions and deletions, and minimize false mappings, while maintaining adequate speed compared to other current methods. We show BFAST can align the amount of data needed to fully resequence a human genome, one billion reads, with high sensitivity and accuracy, on a modest computer cluster in less than 24 hours. BFAST is available at (http://bfast.sourceforge.net).

Mapping the methylome A newly developed method of characterizing an organism's 'methylome', that is the pattern of DNA methylation in the genome, has been used to generate a map of methylated cytosines in Arabidopsis to single base-pair resolution. The procedure, termed BS-Seq, combines bisulphite treatment of genomic DNA with ultra-high-throughput DNA sequencing to achieve a more precise and comprehensive result than previously possible. DNA methylation is an important factor in regulating gene expression, and this method, which can be applied to larger genomes like the mouse as well as to Arabidopsis , could prove a significant advance in the study of this form of gene regulation. In Arabidopsis , a map of methylated cytosines is generated at single base pair resolution by combining bisulphite treatment of genomic DNA with ultra-high-throughput sequencing. Cytosine DNA methylation is important in regulating gene expression and in silencing transposons and other repetitive sequences 1 , 2 . Recent genomic studies in Arabidopsis thaliana have revealed that many endogenous genes are methylated either within their promoters or within their transcribed regions, and that gene methylation is highly correlated with transcription levels 3 , 4 , 5 . However, plants have different types of methylation controlled by different genetic pathways, and detailed information on the methylation status of each cytosine in any given genome is lacking. To this end, we generated a map at single-base-pair resolution of methylated cytosines for Arabidopsis , by combining bisulphite treatment of genomic DNA with ultra-high-throughput sequencing using the Illumina 1G Genome Analyser and Solexa sequencing technology 6 . This approach, termed BS-Seq, unlike previous microarray-based methods, allows one to sensitively measure cytosine methylation on a genome-wide scale within specific sequence contexts. Here we describe methylation on previously inaccessible components of the genome and analyse the DNA methylation sequence composition and distribution. We also describe the effect of various DNA methylation mutants on genome-wide methylation patterns, and demonstrate that our newly developed library construction and computational methods can be applied to large genomes such as that of mouse.

The study of reaction–diffusion processes is much more complicated on general curved surfaces than on standard Cartesian coordinate spaces. Here we show how to formulate and solve systems of reaction–diffusion equations on surfaces in an extremely simple way, using only the standard Cartesian form of differential operators, and a discrete unorganized point set to represent the surface. Our method decouples surface geometry from the underlying differential operators. As a consequence, it becomes possible to formulate and solve rather general reaction–diffusion equations on general surfaces without having to consider the complexities of differential geometry or sophisticated numerical analysis. To illustrate the generality of the method, computations for surface diffusion, pattern formation, excitable media, and bulk-surface coupling are provided for a variety of complex point cloud surfaces.

Schizophrenia is a devastating neurodevelopmental disorder whose genetic influences remain elusive. We hypothesize that individually rare structural variants contribute to the illness. Microdeletions and microduplications >100 kilobases were identified by microarray comparative genomic hybridization of genomic DNA from 150 individuals with schizophrenia and 268 ancestry-matched controls. All variants were validated by high-resolution platforms. Novel deletions and duplications of genes were present in 5% of controls versus 15% of cases and 20% of young-onset cases, both highly significant differences. The association was independently replicated in patients with childhood-onset schizophrenia as compared with their parents. Mutations in cases disrupted genes disproportionately from signaling networks controlling neurodevelopment, including neuregulin and glutamate pathways. These results suggest that multiple, individually rare mutations altering genes in neurodevelopmental pathways contribute to schizophrenia.

Oral-facial-digital (OFD) syndromes are a heterogeneous group of congenital disorders characterized by malformations of the face and oral cavity, and digit anomalies. Mutations within 12 cilia-related genes have been identified that cause several types of OFD, suggesting that OFDs constitute a subgroup of developmental ciliopathies. Through homozygosity mapping and exome sequencing of two families with variable OFD type 2, we identified distinct germline variants in INTS13 , a subunit of the Integrator complex. This multiprotein complex associates with RNA Polymerase II and cleaves nascent RNA to modulate gene expression. We determined that INTS13 utilizes its C-terminus to bind the Integrator cleavage module, which is disrupted by the identified germline variants p.S652L and p.K668Nfs*9. Depletion of INTS13 disrupts ciliogenesis in human cultured cells and causes dysregulation of a broad collection of ciliary genes. Accordingly, its knockdown in Xenopus embryos leads to motile cilia anomalies. Altogether, we show that mutations in INTS13 cause an autosomal recessive ciliopathy, which reveals key interactions between components of the Integrator complex. The integrator complex is required for the synthesis of protein coding and non-coding RNA and contains the protein INTS13. Here, the authors find germline mutations in INTS13 in two families with oral facial digital syndrome and show that the mutation affects the c-terminal domain of the protein and disrupts cilliogenesis.

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.

Background Since the introduction of next-generation DNA sequencers the rapid increase in sequencer throughput, and associated drop in costs, has resulted in more than a dozen human genomes being resequenced over the last few years. These efforts are merely a prelude for a future in which genome resequencing will be commonplace for both biomedical research and clinical applications. The dramatic increase in sequencer output strains all facets of computational infrastructure, especially databases and query interfaces. The advent of cloud computing, and a variety of powerful tools designed to process petascale datasets, provide a compelling solution to these ever increasing demands. Results In this work, we present the SeqWare Query Engine which has been created using modern cloud computing technologies and designed to support databasing information from thousands of genomes. Our backend implementation was built using the highly scalable, NoSQL HBase database from the Hadoop project. We also created a web-based frontend that provides both a programmatic and interactive query interface and integrates with widely used genome browsers and tools. Using the query engine, users can load and query variants (SNVs, indels, translocations, etc) with a rich level of annotations including coverage and functional consequences. As a proof of concept we loaded several whole genome datasets including the U87MG cell line. We also used a glioblastoma multiforme tumor/normal pair to both profile performance and provide an example of using the Hadoop MapReduce framework within the query engine. This software is open source and freely available from the SeqWare project ( http://seqware.sourceforge.net ). Conclusions The SeqWare Query Engine provided an easy way to make the U87MG genome accessible to programmers and non-programmers alike. This enabled a faster and more open exploration of results, quicker tuning of parameters for heuristic variant calling filters, and a common data interface to simplify development of analytical tools. The range of data types supported, the ease of querying and integrating with existing tools, and the robust scalability of the underlying cloud-based technologies make SeqWare Query Engine a nature fit for storing and searching ever-growing genome sequence datasets.

Birgitte Lane and colleagues show that Ferguson-Smith disease, an autosomal dominant skin cancer condition characterized by the development of multiple self-healing tumors, 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 1 , 2 . 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.

Daniel Cohn and colleagues identify mutations in the gene encoding the calcium-permeable cation channel TRPV4 in families with autosomal dominant brachyolmia. Functional studies show that the mutations result in gain-of-function of channel activation. The brachyolmias constitute a clinically and genetically heterogeneous group of skeletal dysplasias characterized by a short trunk, scoliosis and mild short stature 1 . Here, we identify a locus for an autosomal dominant form of brachyolmia on chromosome 12q24.1–12q24.2. Among the genes in the genetic interval, we selected TRPV4 , which encodes a calcium permeable cation channel of the transient receptor potential (TRP) vanilloid family, as a candidate gene because of its cartilage-selective gene expression pattern. In two families with the phenotype, we identified point mutations in TRPV4 that encoded R616Q and V620I substitutions, respectively. Patch clamp studies of transfected HEK cells showed that both mutations resulted in a dramatic gain of function characterized by increased constitutive activity and elevated channel activation by either mechano-stimulation or agonist stimulation by arachidonic acid or the TRPV4-specific agonist 4α-phorbol 12,13-didecanoate (4αPDD). This study thus defines a previously unknown mechanism, activation of a calcium-permeable TRP ion channel, in skeletal dysplasia pathogenesis.

Hospitalizations for asthma exacerbations are common in the United States, but there are no national estimates of outcomes in this population. It is also not known if race disparities in asthma deaths exist among hospitalized patients. To estimate outcomes of patients hospitalized for asthma in the United States and to determine if the risk of death in this population is higher among black patients compared with white patients. We used the Nationwide Inpatient Sample for 2000. Admissions for asthma exacerbations among patients > 5 yr of age were included. Mortality was the primary outcome; secondary outcomes were length of stay and total hospital charges. In-hospital asthma mortality was 0.5% (99% confidence interval [CI], 0.4-0.6), with mean hospital stay of 2.7 d (99% CI, 2.6-2.8 d) and 9,078 dollars (99% CI, 8,300-9,855 dollars) in hospital charges. Deaths in this population accounted for about one-third of all asthma deaths reported in the United States. Black patients hospitalized for asthma exacerbations were less likely to die when compared with white patients (0.3 vs. 0.6%; p < 0.001). However, in multivariable analyses, there were no significant race differences in hospital deaths. Mortality among patients hospitalized for asthma exacerbations accounts for one-third of all deaths from asthma. The higher overall risk of death from asthma in black patients compared with white patients in the United States is not explained by race differences in hospital deaths and therefore is attributable to factors preceding hospitalization.