Accurate determination of node and arc multiplicities in de Bruijn graphs using conditional random fields

Background: De Bruijn graphs are key data structures for the analysis of next-generation sequencing data. They efficiently represent the overlap between reads and hence, also the underlying genome sequence. However, sequencing errors and repeated subsequences render the identification of the true underlying sequence difficult. A key step in this process is the inference of the multiplicities of nodes and arcs in the graph. These multiplicities correspond to the number of times eachk-mer (resp.k+1-mer) implied by a node (resp. arc) is present in the genomic sequence. Determining multiplicities thus reveals the repeat structure and presence of sequencing errors. Multiplicities of nodes/arcs in the de Bruijn graph are reflected in their coverage, however, coverage variability and coverage biases render their determination ambiguous. Current methods to determine node/arc multiplicities base their decisions solely on the information in nodes and arcs individually, under-utilising the information present in the sequencing data. Results: To improve the accuracy with which node and arc multiplicities in a de Bruijn graph are inferred, we developed a conditional random field (CRF) model to efficiently combine the coverage information within each node/arc individually with the information of surrounding nodes and arcs. Multiplicities are thus collectively assigned in a more consistent manner. Conclusions: We demonstrate that the CRF model yields significant improvements in accuracy and a more robust expectation-maximisation parameter estimation. Truek-mers can be distinguished from erroneousk-mers with a higher F(1)score than existing methods. A C++11 implementation is available atunder the GNU AGPL v3.0 license

Accurate determination of node and arc multiplicities in de bruijn graphs using conditional random fields

Background: De Bruijn graphs are key data structures for the analysis of next-generation sequencing data. They efficiently represent the overlap between reads and hence, also the underlying genome sequence. However, sequencing errors and repeated subsequences render the identification of the true underlying sequence difficult. A key step in this process is the inference of the multiplicities of nodes and arcs in the graph. These multiplicities correspond to the number of times eachk-mer (resp.k+1-mer) implied by a node (resp. arc) is present in the genomic sequence. Determining multiplicities thus reveals the repeat structure and presence of sequencing errors. Multiplicities of nodes/arcs in the de Bruijn graph are reflected in their coverage, however, coverage variability and coverage biases render their determination ambiguous. Current methods to determine node/arc multiplicities base their decisions solely on the information in nodes and arcs individually, under-utilising the information present in the sequencing data. Results: To improve the accuracy with which node and arc multiplicities in a de Bruijn graph are inferred, we developed a conditional random field (CRF) model to efficiently combine the coverage information within each node/arc individually with the information of surrounding nodes and arcs. Multiplicities are thus collectively assigned in a more consistent manner. Conclusions: We demonstrate that the CRF model yields significant improvements in accuracy and a more robust expectation-maximisation parameter estimation. Truek-mers can be distinguished from erroneousk-mers with a higher F(1)score than existing methods. A C++11 implementation is available atunder the GNU AGPL v3.0 license

Vacuum biopsy of the breast

Background. Mammotome is a percutaneous biopsy technique with a large needle, assisted by vacuum aspiration. The single insertion technique allows to collect multiple large tissue samples, which gives a more reliable histological result. The technique can be used under stereotactic, MRl or ultrasound guidance. The technique is used under local anesthesia, on an outpatient basis. Conclusions. Altough it is essentially a diagnostic technique, it has a therapeutic potential in some cases. Recently, the technique has been FDA approved for the removal of benign nodular lesions. Worldwide, over 2 million procedures have already been performed. It has the potential to replace surgical open biopsy, with a similar accuracy and at a lower cost, and it is less invasive

Influence of the π–π interaction on the hydrogen bonding capacity of stacked DNA/RNA bases

The interplay between aromatic stacking and hydrogen bonding in nucleobases has been investigated via high-level quantum chemical calculations. The experimentally observed stacking arrangement between consecutive bases in DNA and RNA/DNA double helices is shown to enhance their hydrogen bonding ability as opposed to gas phase optimized complexes. This phenomenon results from more repulsive electrostatic interactions as is demonstrated in a model system of cytosine stacked offset-parallel with substituted benzenes. Therefore, the H-bonding capacity of the N3 and O2 atoms of cytosine increases linearly with the electrostatic repulsion between the stacked rings. The local hardness, a density functional theory-based reactivity descriptor, appears to be a key index associated with the molecular electrostatic potential (MEP) minima around H-bond accepting atoms, and is inversely proportional to the electrostatic interaction between stacked molecules. Finally, the MEP minima on surfaces around the bases in experimental structures of DNA and RNA–DNA double helices show that their hydrogen bonding capacity increases when taking more neighboring (intra-strand) stacking partners into account

Fragile Spectral and Temporal Auditory Processing in Adolescents with Autism Spectrum Disorder and Early Language Delay

We investigated low-level auditory spectral and temporal processing in adolescents with autism spectrum disorder (ASD) and early language delay compared to matched typically developing controls. Auditory measures were designed to target right versus left auditory cortex processing (i.e. frequency discrimination and slow amplitude modulation (AM) detection versus gap-in-noise detection and faster AM detection), and to pinpoint the task and stimulus characteristics underlying putative superior spectral processing in ASD. We observed impaired frequency discrimination in the ASD group and suggestive evidence of poorer temporal resolution as indexed by gap-in-noise detection thresholds. These findings question the evidence of enhanced spectral sensitivity in ASD and do not support the hypothesis of superior right and inferior left hemispheric auditory processing in ASD.University of Leuven. Research Council (Grant IDO/08/013

Diversity in kinetics correlated with structure in nano body-stabilized LacY.

Funder: research foundation-flandersThe structure of lactose permease, stabilized in a periplasmic open conformation by two Gly to Trp replacements (LacYww) and complexed with a nanobody directed against this conformation, provides the highest resolution structure of the symporter. The nanobody binds in a different manner than two other nanobodies made against the same mutant, which also bind to the same general region on the periplasmic side. This region of the protein may represent an immune hotspot. The CDR3 loop of the nanobody is held by hydrogen bonds in a conformation that partially blocks access to the substrate-binding site. As a result, kon and koff for galactoside binding to either LacY or the double mutant complexed with the nanobody are lower than for the other two LacY/nanobody complexes though the Kd values are similar, reflecting the fact that the nanobodies rigidify structures along the pathway. While the wild-type LacY/nanobody complex clearly stabilizes a similar 'extracellular open' conformation in solution, judged by binding kinetics, the complex with wild-type LacY did not yet crystallize, suggesting the nanobody does not bind strongly enough to shift the equilibrium to stabilize a periplasmic side-open conformation suitable for crystallization. However, the similarity of the galactoside binding kinetics for the nanobody-bound complexes with wild type LacY and with LacYWW indicates that they have similar structures, showing that the reported co-structures reliably show nanobody interactions with LacY

In vitro reconstitution of dynamically interacting integral membrane subunits of energy-coupling factor transporters

Energy-coupling factor (ECF) transporters mediate import of micronutrients in prokaryotes. They consist of an integral membrane S-component (that binds substrate) and ECF module (that powers transport by ATP hydrolysis). It has been proposed that different S-components compete for docking onto the same ECF module, but a minimal liposome-reconstituted system, required to substantiate this idea, is lacking. Here, we co-reconstituted ECF transporters for folate (ECF-FolT2) and pantothenate (ECF-PanT) into proteoliposomes, and assayed for crosstalk during active transport. The kinetics of transport showed that exchange of S-components is part of the transport mechanism. Competition experiments suggest much slower substrate association with FolT2 than with PanT. Comparison of a crystal structure of ECF-PanT with previously determined structures of ECF-FolT2 revealed larger conformational changes upon binding of folate than pantothenate, which could explain the kinetic differences. Our work shows that a minimal in vitro system with two reconstituted transporters recapitulates intricate kinetics behaviour observed in vivo

The Molecular Mechanism of Transport by the Mitochondrial ADP/ATP Carrier.

Mitochondrial ADP/ATP carriers transport ADP into the mitochondrial matrix for ATP synthesis, and ATP out to fuel the cell, by cycling between cytoplasmic-open and matrix-open states. The structure of the cytoplasmic-open state is known, but it has proved difficult to understand the transport mechanism in the absence of a structure in the matrix-open state. Here, we describe the structure of the matrix-open state locked by bongkrekic acid bound in the ADP/ATP-binding site at the bottom of the central cavity. The cytoplasmic side of the carrier is closed by conserved hydrophobic residues, and a salt bridge network, braced by tyrosines. Glycine and small amino acid residues allow close-packing of helices on the matrix side. Uniquely, the carrier switches between states by rotation of its three domains about a fulcrum provided by the substrate-binding site. Because these features are highly conserved, this mechanism is likely to apply to the whole mitochondrial carrier family. VIDEO ABSTRACT