Modeling the evolution space of breakage fusion bridge cycles with a stochastic folding process

Breakage-Fusion-Bridge cycles in cancer arise when a broken segment of DNA is duplicated and an end from each copy joined together. This structure then 'unfolds' into a new piece of palindromic DNA. This is one mechanism responsible for the localised amplicons observed in cancer genome data. The process has parallels with paper folding sequences that arise when a piece of paper is folded several times and then unfolded. Here we adapt such methods to study the breakage-fusion-bridge structures in detail. We firstly consider discrete representations of this space with 2-d trees to demonstrate that there are 2^(n(n-1)/2) qualitatively distinct evolutions involving n breakage-fusion-bridge cycles. Secondly we consider the stochastic nature of the fold positions, to determine evolution likelihoods, and also describe how amplicons become localised. Finally we highlight these methods by inferring the evolution of breakage-fusion-bridge cycles with data from primary tissue cancer samples

High-throughput, quantitative analyses of genetic interactions in E. coli.

Large-scale genetic interaction studies provide the basis for defining gene function and pathway architecture. Recent advances in the ability to generate double mutants en masse in Saccharomyces cerevisiae have dramatically accelerated the acquisition of genetic interaction information and the biological inferences that follow. Here we describe a method based on F factor-driven conjugation, which allows for high-throughput generation of double mutants in Escherichia coli. This method, termed genetic interaction analysis technology for E. coli (GIANT-coli), permits us to systematically generate and array double-mutant cells on solid media in high-density arrays. We show that colony size provides a robust and quantitative output of cellular fitness and that GIANT-coli can recapitulate known synthetic interactions and identify previously unidentified negative (synthetic sickness or lethality) and positive (suppressive or epistatic) relationships. Finally, we describe a complementary strategy for genome-wide suppressor-mutant identification. Together, these methods permit rapid, large-scale genetic interaction studies in E. coli

Modelled glacier response to centennial temperature and precipitation trends on the Antarctic Peninsula

The northern Antarctic Peninsula is currently undergoing rapid atmospheric warming. Increased glacier-surface melt during the twentieth century has contributed to ice-shelf collapse and the widespread acceleration, thinning and recession of glaciers. Therefore, glaciers peripheral to the Antarctic Ice Sheet currently make a large contribution to eustatic sea-level rise, but future melting may be offset by increased precipitation. Here we assess glacier-climate relationships both during the past and into the future, using ice-core and geological data and glacier and climate numerical model simulations. Focusing on Glacier IJR45 on James Ross Island, northeast Antarctic Peninsula, our modelling experiments show that this representative glacier is most sensitive to temperature change, not precipitation change. We determine that its most recent expansion occurred during the late Holocene a Little Ice Age' and not during the warmer mid-Holocene, as previously proposed. Simulations using a range of future Intergovernmental Panel on Climate Change climate scenarios indicate that future increases in precipitation are unlikely to offset atmospheric-warming-induced melt of peripheral Antarctic Peninsula glaciers

Influence of the incremental step size in work rate on exercise response and gas exchange in patients with pulmonary hypertension

<p>Abstract</p> <p>Background</p> <p>Cardiopulmonary exercise testing (CPET) has become increasingly important as a routine procedure in daily clinical work. So far, it is generally accepted that an individualized exercise protocol with exercise duration of 6 to 12 minutes is preferable to assess maximal exercise performance. The aim of this study was to compare an individualized NYHA adapted exercise protocol with a fixed standard protocol in patients with severe pulmonary arterial hypertension.</p> <p>Methods</p> <p>Twenty-two patients (17 female, 5 male; mean age 49 ± 14 yrs) underwent symptom limited CPET on a bicycle. On two consecutive days each subject performed a stepwise CPET according to a modified Jones protocol (16 Watt per minute stages) as well as an individualized NYHA adapted protocol with 5 or 10 Watt/min stages in a randomized order. Oxygen uptake at peak exercise (peakVO₂) and anaerobic threshold (VO₂AT), maximal ventilation (VE), breathing reserve (VE/MVV), ventilatory efficiency (VE vs. VCO₂slope), exercise time, maximal power and work rate were assessed and compared between both protocols.</p> <p>Results</p> <p>Comparing both, adapted NYHA protocol and standardized Jones protocol, we found significant differences in maximal power (56.7 ± 19 W vs. 74 ± 18 W; p < 0.001) and exercise time (332 ± 107 sec. vs. 248 ± 72 sec.; p < 0.001). In contrast, no significant differences were obvious comparing both protocols concerning work rate, VE, VE/MVV, peakVO₂, VO₂AT and VE vs. VCO₂slope.</p> <p>Conclusion</p> <p>Variations of incremental step size during CPET significantly affect exercise time and maximal power, whereas relevant parameters for clinical judgement and prognosis such as oxygen uptake, ventilation and ventilatory efficiency remain unchanged. These findings have practical implications for the exercise evaluation of patients with pulmonary hypertension. To reach maximal results for ventilation, oxygen uptake and gas exchange an individualization of incremental step size appears not to be mandatory.</p

Homozygous microdeletion of exon 5 in ZNF277 in a girl with specific language impairment

Peer reviewedPublisher PD

Reconstructing cancer genomes from paired-end sequencing data

<p>Abstract</p> <p>Background</p> <p>A cancer genome is derived from the germline genome through a series of somatic mutations. Somatic structural variants - including duplications, deletions, inversions, translocations, and other rearrangements - result in a cancer genome that is a scrambling of intervals, or "blocks" of the germline genome sequence. We present an efficient algorithm for reconstructing the block organization of a cancer genome from paired-end DNA sequencing data.</p> <p>Results</p> <p>By aligning paired reads from a cancer genome - and a matched germline genome, if available - to the human reference genome, we derive: (i) a partition of the reference genome into intervals; (ii) adjacencies between these intervals in the cancer genome; (iii) an estimated copy number for each interval. We formulate the Copy Number and Adjacency Genome Reconstruction Problem of determining the cancer genome as a sequence of the derived intervals that is consistent with the measured adjacencies and copy numbers. We design an efficient algorithm, called Paired-end Reconstruction of Genome Organization (PREGO), to solve this problem by reducing it to an optimization problem on an interval-adjacency graph constructed from the data. The solution to the optimization problem results in an Eulerian graph, containing an alternating Eulerian tour that corresponds to a cancer genome that is consistent with the sequencing data. We apply our algorithm to five ovarian cancer genomes that were sequenced as part of The Cancer Genome Atlas. We identify numerous rearrangements, or structural variants, in these genomes, analyze reciprocal vs. non-reciprocal rearrangements, and identify rearrangements consistent with known mechanisms of duplication such as tandem duplications and breakage/fusion/bridge (B/F/B) cycles.</p> <p>Conclusions</p> <p>We demonstrate that PREGO efficiently identifies complex and biologically relevant rearrangements in cancer genome sequencing data. An implementation of the PREGO algorithm is available at <url>http://compbio.cs.brown.edu/software/</url>.</p

Identifying work related injuries: comparison of methods for interrogating text fields

Background: Work-related injuries in Australia are estimated to cost around $57.5 billion annually, however there are currently insufficient surveillance data available to support an evidence-based public health response. Emergency departments (ED) in Australia are a potential source of information on work-related injuries though most ED’s do not have an ‘Activity Code’ to identify work-related cases with information about the presenting problem recorded in a short free text field. This study compared methods for interrogating text fields for identifying work-related injuries presenting at emergency departments to inform approaches to surveillance of work-related injury.---------- Methods: Three approaches were used to interrogate an injury description text field to classify cases as work-related: keyword search, index search, and content analytic text mining. Sensitivity and specificity were examined by comparing cases flagged by each approach to cases coded with an Activity code during triage. Methods to improve the sensitivity and/or specificity of each approach were explored by adjusting the classification techniques within each broad approach.---------- Results: The basic keyword search detected 58% of cases (Specificity 0.99), an index search detected 62% of cases (Specificity 0.87), and the content analytic text mining (using adjusted probabilities) approach detected 77% of cases (Specificity 0.95).---------- Conclusions The findings of this study provide strong support for continued development of text searching methods to obtain information from routine emergency department data, to improve the capacity for comprehensive injury surveillance

Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates

Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemination of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here

A mosaic tetracycline resistance gene tet(S/M) detected in an MDR pneumococcal CC230 lineage that underwent capsular switching in South Africa

Objectives: We reported tet(S/M) in Streptococcus pneumoniae and investigated its temporal spread in relation to nationwide clinical interventions. Methods: We whole-genome sequenced 12 254 pneumococcal isolates from 29 countries on an Illumina HiSeq sequencer. Serotype, multilocus ST and antibiotic resistance were inferred from genomes. An SNP tree was built using Gubbins. Temporal spread was reconstructed using a birth–death model. Results: We identified tet(S/M) in 131 pneumococcal isolates and none carried other known tet genes. Tetracycline susceptibility testing results were available for 121 tet(S/M)-positive isolates and all were resistant. A majority (74%) of tet(S/M)-positive isolates were from South Africa and caused invasive diseases among young children (59% HIV positive, where HIV status was available). All but two tet(S/M)-positive isolates belonged to clonal complex (CC) 230. A global phylogeny of CC230 (n=389) revealed that tet(S/M)-positive isolates formed a sublineage predicted to exhibit resistance to penicillin, co-trimoxazole, erythromycin and tetracycline. The birth–death model detected an unrecognized outbreak of this sublineage in South Africa between 2000 and 2004 with expected secondary infections (effective reproductive number, R) of ∼2.5. R declined to ∼1.0 in 2005 and <1.0 in 2012. The declining epidemic could be related to improved access to ART in 2004 and introduction of pneumococcal conjugate vaccine (PCV) in 2009. Capsular switching from vaccine serotype 14 to non-vaccine serotype 23A was observed within the sublineage. Conclusions: The prevalence of tet(S/M) in pneumococci was low and its dissemination was due to an unrecognized outbreak of CC230 in South Africa. Capsular switching in this MDR sublineage highlighted its potential to continue to cause disease in the post-PCV13 era

Tracing the dynamic life story of a Bronze Age Female

YesAncient human mobility at the individual level is conventionally studied by the diverse application of suitable techniques (e.g. aDNA, radiogenic strontium isotopes, as well as oxygen and lead isotopes) to either hard and/or soft tissues. However, the limited preservation of coexisting hard and soft human tissues hampers the possibilities of investigating high-resolution diachronic mobility periods in the life of a single individual. Here, we present the results of a multidisciplinary study of an exceptionally well preserved circa 3.400-year old Danish Bronze Age female find, known as the Egtved Girl. We applied biomolecular, biochemical and geochemical analyses to reconstruct her mobility and diet. We demonstrate that she originated from a place outside present day Denmark (the island of Bornholm excluded), and that she travelled back and forth over large distances during the final months of her life, while consuming a terrestrial diet with intervals of reduced protein intake. We also provide evidence that all her garments were made of non-locally produced wool. Our study advocates the huge potential of combining biomolecular and biogeochemical provenance tracer analyses to hard and soft tissues of a single ancient individual for the reconstruction of high-resolution human mobility.The Danish National Research Foundation; The Carlsberg Foundation, L'Oreal Denmark-UNESCO; The ERC agreement no. 26944