Efficient parallel and out of core algorithms for constructing large bi-directed de Bruijn graphs

<p>Abstract</p> <p>Background</p> <p>Assembling genomic sequences from a set of overlapping reads is one of the most fundamental problems in computational biology. Algorithms addressing the assembly problem fall into two broad categories <b>- </b>based on the data structures which they employ. The first class uses an overlap/string graph and the second type uses a de Bruijn graph. However with the recent advances in short read sequencing technology, de Bruijn graph based algorithms seem to play a vital role in practice. Efficient algorithms for building these massive de Bruijn graphs are very essential in large sequencing projects based on short reads. In an earlier work, an <it>O</it>(<it>n/p</it>) time parallel algorithm has been given for this problem. Here <it>n </it>is the size of the input and <it>p </it>is the number of processors. This algorithm enumerates all possible bi-directed edges which can overlap with a node and ends up generating Θ(<it>n</it>Σ) messages (Σ being the size of the alphabet).</p> <p>Results</p> <p>In this paper we present a Θ(<it>n/p</it>) time parallel algorithm with a communication complexity that is equal to that of parallel sorting and is not sensitive to Σ. The generality of our algorithm makes it very easy to extend it even to the out-of-core model and in this case it has an optimal I/O complexity of <inline-formula><m:math name="1471-2105-11-560-i1" xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mrow><m:mo>Θ</m:mo><m:mo stretchy="false">(</m:mo><m:mfrac><m:mrow><m:mi>n</m:mi><m:mi>log</m:mi><m:mo stretchy="false">(</m:mo><m:mi>n</m:mi><m:mo>/</m:mo><m:mi>B</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mrow><m:mi>B</m:mi><m:mi>log</m:mi><m:mo stretchy="false">(</m:mo><m:mi>M</m:mi><m:mo>/</m:mo><m:mi>B</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:mfrac><m:mo stretchy="false">)</m:mo></m:mrow></m:math></inline-formula> (<it>M </it>being the main memory size and <it>B </it>being the size of the disk block). We demonstrate the scalability of our parallel algorithm on a SGI/Altix computer. A comparison of our algorithm with the previous approaches reveals that our algorithm is faster <b>- </b>both asymptotically and practically. We demonstrate the scalability of our sequential out-of-core algorithm by comparing it with the algorithm used by VELVET to build the bi-directed de Bruijn graph. Our experiments reveal that our algorithm can build the graph with a constant amount of memory, which clearly outperforms VELVET. We also provide efficient algorithms for the bi-directed chain compaction problem.</p> <p>Conclusions</p> <p>The bi-directed de Bruijn graph is a fundamental data structure for any sequence assembly program based on Eulerian approach. Our algorithms for constructing Bi-directed de Bruijn graphs are efficient in parallel and out of core settings. These algorithms can be used in building large scale bi-directed de Bruijn graphs. Furthermore, our algorithms do not employ any all-to-all communications in a parallel setting and perform better than the prior algorithms. Finally our out-of-core algorithm is extremely memory efficient and can replace the existing graph construction algorithm in VELVET.</p

Buttressing staples with cholecyst-derived extracellular matrix (CEM) reinforces staple lines in an ex vivo peristaltic inflation model

This is the author's accepted manuscript. The final published article is available from the link below. Copyright @ Springer Science + Business Media, LLC 2008Background - Staple line leakage and bleeding are the most common problems associated with the use of surgical staplers for gastrointestinal resection and anastomotic procedures. These complications can be reduced by reinforcing the staple lines with buttressing materials. The current study reports the potential use of cholecyst-derived extracellular matrix (CEM) in non-crosslinked (NCEM) and crosslinked (XCEM) forms, and compares their mechanical performance with clinically available buttress materials [small intestinal submucosa (SIS) and bovine pericardium (BP)] in an ex vivo small intestine model. Methods - Three crosslinked CEM variants (XCEM0005, XCEM001, and XCEM0033) with different degree of crosslinking were produced. An ex vivo peristaltic inflation model was established. Porcine small intestine segments were stapled on one end, using buttressed or non-buttressed surgical staplers. The opened, non-stapled ends were connected to a peristaltic pump and pressure transducer and sealed. The staple lines were then exposed to increased intraluminal pressure in a peristaltic manner. Both the leak and burst pressures of the test specimens were recorded. Results - The leak pressures observed for non-crosslinked NCEM (137.8 ± 22.3 mmHg), crosslinked XCEM0005 (109.1 ± 14.1 mmHg), XCEM001 (150.1 ± 16.0 mmHg), XCEM0033 (98.8 ± 10.5 mmHg) reinforced staple lines were significantly higher when compared to non-buttressed control (28.3 ± 10.8 mmHg) and SIS (one and four layers) (62.6 ± 11.8 and 57.6 ± 12.3 mmHg, respectively) buttressed staple lines. NCEM and XCEM were comparable to that observed for BP buttressed staple lines (138.8 ± 3.6 mmHg). Only specimens with reinforced staple lines were able to achieve high intraluminal pressures (ruptured at the intestinal mesentery), indicating that buttress reinforcements were able to withstand pressure higher than that of natural tissue (physiological failure). Conclusions - These findings suggest that the use of CEM and XCEM as buttressing materials is associated with reinforced staple lines and increased leak pressures when compared to non-buttressed staple lines. CEM and XCEM were found to perform comparably with clinically available buttress materials in this ex vivo model.Enterprise Irelan

The influence of the landscape structure within buffer zones, catchment land use and instream environmental variables on mollusc communities in a medium-sized lowland river

The world’s freshwater molluscan fauna is facing unprecedented threats from habitat loss and degradation. Declines in native populations are mostly attributed to the human impact, which results in reduced water quality. The objectives of our survey were to analyse the structure of the mollusc communities in a medium-sized lowland river and to determine the most important environmental variables at different spatial scales, including landscape structure, catchment land use and instream environmental factors that influence their structure. Our survey showed that a medium-sized river, that flows through areas included in the European Ecological Natura 2000 Network Programme of protected sites, provides diverse instream habitats and niches that support 47 mollusc species including Unio crassus, a bivalve of Community interest, whose conservation requires the designation of a special conservation area under the Habitats Directive Natura 2000. This survey showed that mollusc communities are impacted by several environmental variables that act together at multiple scales. The landscape structure within buffer zones, catchment land use and instream environmental variables were all important and influenced the structure of mollusc communities. Therefore, they should all be taken into consideration in the future restoration of the river, future management projects and programmes for the conservation of biodiversity in running waters. The results of this study may be directly applicable for the rehabilitation of river ecosystems and are recommended to stakeholders in their future decision concerning landscape planning, monitoring species and their habitats, conservation plans and management in accordance with the requirements of sustainable development

Temporal Dynamics of European Bat Lyssavirus Type 1 and Survival of Myotis myotis Bats in Natural Colonies

Many emerging RNA viruses of public health concern have recently been detected in bats. However, the dynamics of these viruses in natural bat colonies is presently unknown. Consequently, prediction of the spread of these viruses and the establishment of appropriate control measures are hindered by a lack of information. To this aim, we collected epidemiological, virological and ecological data during a twelve-year longitudinal study in two colonies of insectivorous bats (Myotis myotis) located in Spain and infected by the most common bat lyssavirus found in Europe, the European bat lyssavirus subtype 1 (EBLV-1). This active survey demonstrates that cyclic lyssavirus infections occurred with periodic oscillations in the number of susceptible, immune and infected bats. Persistence of immunity for more than one year was detected in some individuals. These data were further used to feed models to analyze the temporal dynamics of EBLV-1 and the survival rate of bats. According to these models, the infection is characterized by a predicted low basic reproductive rate (R0 = 1.706) and a short infectious period (D = 5.1 days). In contrast to observations in most non-flying animals infected with rabies, the survival model shows no variation in mortality after EBLV-1 infection of M. myotis. These findings have considerable public health implications in terms of management of colonies where lyssavirus-positive bats have been recorded and confirm the potential risk of rabies transmission to humans. A greater understanding of the dynamics of lyssavirus in bat colonies also provides a model to study how bats contribute to the maintenance and transmission of other viruses of public health concern

Stress-Induced Reinstatement of Drug Seeking: 20 Years of Progress

In human addicts, drug relapse and craving are often provoked by stress. Since 1995, this clinical scenario has been studied using a rat model of stress-induced reinstatement of drug seeking. Here, we first discuss the generality of stress-induced reinstatement to different drugs of abuse, different stressors, and different behavioral procedures. We also discuss neuropharmacological mechanisms, and brain areas and circuits controlling stress-induced reinstatement of drug seeking. We conclude by discussing results from translational human laboratory studies and clinical trials that were inspired by results from rat studies on stress-induced reinstatement. Our main conclusions are (1) The phenomenon of stress-induced reinstatement, first shown with an intermittent footshock stressor in rats trained to self-administer heroin, generalizes to other abused drugs, including cocaine, methamphetamine, nicotine, and alcohol, and is also observed in the conditioned place preference model in rats and mice. This phenomenon, however, is stressor specific and not all stressors induce reinstatement of drug seeking. (2) Neuropharmacological studies indicate the involvement of corticotropin-releasing factor (CRF), noradrenaline, dopamine, glutamate, kappa/dynorphin, and several other peptide and neurotransmitter systems in stress-induced reinstatement. Neuropharmacology and circuitry studies indicate the involvement of CRF and noradrenaline transmission in bed nucleus of stria terminalis and central amygdala, and dopamine, CRF, kappa/dynorphin, and glutamate transmission in other components of the mesocorticolimbic dopamine system (ventral tegmental area, medial prefrontal cortex, orbitofrontal cortex, and nucleus accumbens). (3) Translational human laboratory studies and a recent clinical trial study show the efficacy of alpha-2 adrenoceptor agonists in decreasing stress-induced drug craving and stress-induced initial heroin lapse

High Pro-Inflammatory Cytokine Secretion and Loss of High Avidity Cross-Reactive Cytotoxic T-Cells during the Course of Secondary Dengue Virus Infection

BACKGROUND: Dengue is one of the most important human diseases transmitted by an arthropod vector and the incidence of dengue virus infection has been increasing - over half the world's population now live in areas at risk of infection. Most infections are asymptomatic, but a subset of patients experience a potentially fatal shock syndrome characterised by plasma leakage. Severe forms of dengue are epidemiologically associated with repeated infection by more than one of the four dengue virus serotypes. Generally attributed to the phenomenon of antibody-dependent enhancement, recent observations indicate that T-cells may also influence disease phenotype. METHODS AND FINDINGS: Virus-specific cytotoxic T lymphocytes (CTL) showing high level cross reactivity between dengue serotypes could be expanded from blood samples taken during the acute phase of secondary dengue infection. These could not be detected in convalescence when only CTL populations demonstrating significant serotype specificity were identified. Dengue cross-reactive CTL clones derived from these patients were of higher avidity than serotype-specific clones and produced much higher levels of both type 1 and certain type 2 cytokines, many previously implicated in dengue pathogenesis. CONCLUSION: Dengue serotype cross-reactive CTL clones showing high avidity for antigen produce higher levels of inflammatory cytokines than serotype-specific clones. That such cells cannot be expanded from convalescent samples suggests that they may be depleted, perhaps as a consequence of activation-induced cell death. Such high avidity cross-reactive memory CTL may produce inflammatory cytokines during the course of secondary infection, contributing to the pathogenesis of vascular leak. These cells appear to be subsequently deleted leaving a more serotype-specific memory CTL pool. Further studies are needed to relate these cellular observations to disease phenotype in a large group of patients. If confirmed they have significant implications for understanding the role of virus-specific CTL in pathogenesis of dengue disease

The Warburg Effect Suppresses Oxidative Stress Induced Apoptosis in a Yeast Model for Cancer

BACKGROUND: Otto Warburg observed that cancer cells are often characterized by intense glycolysis in the presence of oxygen and a concomitant decrease in mitochondrial respiration. Research has mainly focused on a possible connection between increased glycolysis and tumor development whereas decreased respiration has largely been left unattended. Therefore, a causal relation between decreased respiration and tumorigenesis has not been demonstrated. METHODOLOGY/PRINCIPAL FINDINGS: For this purpose, colonies of Saccharomyces cerevisiae, which is suitable for manipulation of mitochondrial respiration and shows mitochondria-mediated cell death, were used as a model. Repression of respiration as well as ROS-scavenging via glutathione inhibited apoptosis and conferred a survival advantage during seeding and early development of this fast proliferating solid cell population. In contrast, enhancement of respiration triggered cell death. CONCLUSION/SIGNIFICANCE: Thus, the Warburg effect might directly contribute to the initiation of cancer formation--not only by enhanced glycolysis--but also via decreased respiration in the presence of oxygen, which suppresses apoptosis

Biology and conservation of freshwater bivalves : past, present and future perspectives

Freshwater bivalves have been highly threatened by human activities, and recently their global decline has been causing conservational and social concern. In this paper, we review the most important research events in freshwater bivalve biology calling attention to the main scientific achievements. A great bias exists in the research effort, with much more information available for bivalve species belonging to the Unionida in comparison to other groups. The same is true for the origin of these studies, since the publishing pattern does not always correspond to the hotspots of biodiversity but is concentrated in the northern hemisphere mainly in North America, Europe and Russia, with regions such as Africa and Southeast Asia being quite understudied. We also summarize information about past, present and future perspectives concerning the most important research topics that include taxonomy, systematics, anatomy, physiology, ecology and conservation of freshwater bivalves. Finally, we introduce the articles published in this Hydrobiologia special issue related with the International Meeting on Biology and Conservation of Freshwater Bivalves held in 2012 in Braganc¸a, Portugal.We would like to express our gratitude to our sponsors and institutions, especially to the Polytechnic Institute of Braganca for all the logistic support. We acknowledge all keynote speakers, authors, session chairpersons and especially to all attendees whose contributions were fundamental for the success of this meeting. We would also like to thank all referees of this special issue and to Koen Martens, Editor-in-Chief of Hydrobiologia, for all the valuable comments and suggestions. The chronogram was built with the help of the expert opinion of fellow colleagues Rafael Araujo, Arthur Bogan, Kevin Cummings, Dan Graf, Wendell Haag, Karl-Otto Nagel and David Strayer to whom we are very grateful. The authors acknowledge the support provided by Portuguese Foundation for Science and Technology (FCT) and COMPETE funds-projects CONBI (Contract: PTDC/AAC-AMB/117688/2010) and ECO-IAS (Contract: PTDC/AAC-AMB/116685/2010), and by the European Regional Development Fund (ERDF) through the COMPETE, under the project "PEst-C/MAR/LA0015/2011"