Plantagora: modeling whole genome sequencing and assembly of plant genomes

Background

Genomics studies are being revolutionized by the next generation sequencing technologies, which have made whole genome sequencing much more accessible to the average researcher. Whole genome sequencing with the new technologies is a developing art that, despite the large volumes of data that can be produced, may still fail to provide a clear and thorough map of a genome. The Plantagora project was conceived to address specifically the gap between having the technical tools for genome sequencing and knowing precisely the best way to use them.

Methodology/Principal Findings

For Plantagora, a platform was created for generating simulated reads from several different plant genomes of different sizes. The resulting read files mimicked either 454 or Illumina reads, with varying paired end spacing. Thousands of datasets of reads were created, most derived from our primary model genome, rice chromosome one. All reads were assembled with different software assemblers, including Newbler, Abyss, and SOAPdenovo, and the resulting assemblies were evaluated by an extensive battery of metrics chosen for these studies. The metrics included both statistics of the assembly sequences and fidelity-related measures derived by alignment of the assemblies to the original genome source for the reads. The results were presented in a website, which includes a data graphing tool, all created to help the user compare rapidly the feasibility and effectiveness of different sequencing and assembly strategies prior to testing an approach in the lab. Some of our own conclusions regarding the different strategies were also recorded on the website.

Conclusions/Significance

Plantagora provides a substantial body of information for comparing different approaches to sequencing a plant genome, and some conclusions regarding some of the specific approaches. Plantagora also provides a platform of metrics and tools for studying the process of sequencing and assembly further.     

A comprehensive genome‐scale reconstruction of Escherichia coli metabolism—2011

The initial genome‐scale reconstruction of the metabolic network of Escherichia coli K‐12 MG1655 was assembled in 2000. It has been updated and periodically released since then based on new and curated genomic and biochemical knowledge. An update has now been built, named iJO1366, which accounts for 1366 genes, 2251 metabolic reactions, and 1136 unique metabolites. iJO1366 was (1) updated in part using a new experimental screen of 1075 gene knockout strains, illuminating cases where alternative pathways and isozymes are yet to be discovered, (2) compared with its predecessor and to experimental data sets to confirm that it continues to make accurate phenotypic predictions of growth on different substrates and for gene knockout strains, and (3) mapped to the genomes of all available sequenced E. coli strains, including pathogens, leading to the identification of hundreds of unannotated genes in these organisms. Like its predecessors, the iJO1366 reconstruction is expected to be widely deployed for studying the systems biology of E. coli and for metabolic engineering applications.     

A hypomethylating variant of MTHFR, 677C>T, blunts the neural response to errors in patients with schizophrenia and healthy individuals

Background

Responding to errors is a critical first step in learning from mistakes, a process that is abnormal in schizophrenia. To gain insight into the neural and molecular mechanisms of error processing, we used functional MRI to examine effects of a genetic variant in methylenetetrahydrofolate reductase (MTHFR 677C>T, rs1801133) that increases risk for schizophrenia and that has been specifically associated with increased perseverative errors among patients. MTHFR is a key regulator of the intracellular one-carbon milieu, including DNA methylation, and each copy of the 677T allele reduces MTHFR activity by 35%.

Methodology/Principal Findings

Using an antisaccade paradigm, we found that the 677T allele induces a dose-dependent blunting of dorsal anterior cingulate cortex (dACC) activation in response to errors, a pattern that was identical in healthy individuals and patients with schizophrenia. Further, the normal relationship between dACC activation and error rate was disrupted among carriers of the 677T allele.

Conclusions/Significance

These findings implicate an epigenetic mechanism in the neural response to errors, and provide insight into normal cognitive variation through a schizophrenia risk gene.     

Lowe Syndrome protein OCRL1 supports maturation of polarized epithelial cells

Mutations in the inositol polyphosphate 5-phosphatase OCRL1 cause Lowe Syndrome, leading to cataracts, mental retardation and renal failure. We noted that cell types affected in Lowe Syndrome are highly polarized, and therefore we studied OCRL1 in epithelial cells as they mature from isolated individual cells into polarized sheets and cysts with extensive communication between neighbouring cells. We show that a proportion of OCRL1 targets intercellular junctions at the early stages of their formation, co-localizing both with adherens junctional components and with tight junctional components. Correlating with this distribution, OCRL1 forms complexes with junctional components α-catenin and zonula occludens (ZO)-1/2/3. Depletion of OCRL1 in epithelial cells growing as a sheet inhibits maturation; cells remain flat, fail to polarize apical markers and also show reduced proliferation. The effect on shape is reverted by re-expressed OCRL1 and requires the 5'-phosphatase domain, indicating that down-regulation of 5-phosphorylated inositides is necessary for epithelial development. The effect of OCRL1 in epithelial maturation is seen more strongly in 3-dimensional cultures, where epithelial cells lacking OCRL1 not only fail to form a central lumen, but also do not have the correct intracellular distribution of ZO-1, suggesting that OCRL1 functions early in the maturation of intercellular junctions when cells grow as cysts. A role of OCRL1 in junctions of polarized cells may explain the pattern of organs affected in Lowe Syndrome.     

Levels and patterns of nucleotide variation in domestication QTL regions on rice chromosome 3 suggest lineage-specific selection

Oryza sativa or Asian cultivated rice is one of the major cereal grass species domesticated for human food use during the Neolithic. Domestication of this species from the wild grass Oryza rufipogon was accompanied by changes in several traits, including seed shattering, percent seed set, tillering, grain weight, and flowering time. Quantitative trait locus (QTL) mapping has identified three genomic regions in chromosome 3 that appear to be associated with these traits. We would like to study whether these regions show signatures of selection and whether the same genetic basis underlies the domestication of different rice varieties. Fragments of 88 genes spanning these three genomic regions were sequenced from multiple accessions of two major varietal groups in O. sativa--indica and tropical japonica--as well as the ancestral wild rice species O. rufipogon. In tropical japonica, the levels of nucleotide variation in these three QTL regions are significantly lower compared to genome-wide levels, and coalescent simulations based on a complex demographic model of rice domestication indicate that these patterns are consistent with selection. In contrast, there is no significant reduction in nucleotide diversity in the homologous regions in indica rice. These results suggest that there are differences in the genetic and selective basis for domestication between these two Asian rice varietal groups.     

Association of reactive oxygen species-mediated signal transduction with in vitro apoptosis sensitivity in chronic lymphocytic leukemia B cells

Background

Chronic lymphocytic leukemia (CLL) is a B cell malignancy with a variable clinical course and unpredictable response to therapeutic agents. Single cell network profiling (SCNP) utilizing flow cytometry measures alterations in signaling biology in the context of molecular changes occurring in malignancies. In this study SCNP was used to identify proteomic profiles associated with in vitro apoptotic responsiveness of CLL B cells to fludarabine, as a basis for ultimately linking these with clinical outcome.

Methodology/Principal Finding

SCNP was used to quantify modulated-signaling of B cell receptor (BCR) network proteins and in vitro F-ara-A mediated apoptosis in 23 CLL samples. Of the modulators studied the reactive oxygen species, hydrogen peroxide (H₂O₂), a known intracellular second messenger and a general tyrosine phosphatase inhibitor stratified CLL samples into two sub-groups based on the percentage of B cells in a CLL sample with increased phosphorylation of BCR network proteins. Separately, in the same patient samples, in vitro exposure to F-ara-A also identified two sub-groups with B cells showing competence or refractoriness to apoptotic induction. Statistical analysis showed that in vitro F-ara-A apoptotic proficiency was highly associated with the proficiency of CLL B cells to undergo H₂O₂-augmented signaling.

Conclusions/Significance

This linkage in CLL B cells among the mechanisms governing chemotherapy-induced apoptosis increased signaling of BCR network proteins and a likely role of phosphatase activity suggests a means of stratifying patients for their response to F-ara-A based regimens. Future studies will examine the clinical applicability of these findings and also the utility of this approach in relating mechanism to function of therapeutic agents.