Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive ...numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3-5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.
The continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) ...tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents.
We performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T isoniazid; gidB 130 bp deletion streptomycin; 1957 95% highest posterior density (HPD): 1937-1971), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P rifampicin; pncA 1 bp insertion pyrazinamide; 1984 95% HPD: 1974-1992) and 10 y (rpoB D435G rifampicin; rrs 1400 kanamycin; gyrA A90V ofloxacin; 1995 95% HPD: 1988-1999) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe.
In the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic era. Subsequent accumulation of stepwise resistance mutations, occurring over decades and prior to the explosion of HIV in this region, yielded MDR and XDR, permitting the emergence of compensatory mutations. Our results suggest that drug-resistant strains circulating today reflect not only vulnerabilities of current TB control efforts but also those that date back 50 y. In drug-resistant TB, isoniazid resistance was overwhelmingly the initial resistance mutation to be acquired, which would not be detected by current rapid molecular diagnostics employed in South Africa that assess only rifampicin resistance.
Cerebral β-amyloidosis and associated pathologies can be exogenously induced by the intracerebral injection of small amounts of pathogenic Aβ-containing brain extract into young β-amyloid precursor ...protein (APP) transgenic mice. The probable β-amyloid-inducing factor in the brain extract has been identified as a species of aggregated Aβ that is generated in its most effective conformation or composition in vivo. Here we report that Aβ in the brain extract is more proteinase K (PK) resistant than is synthetic fibrillar Aβ, and that this PK-resistant fraction of the brain extract retains the capacity to induce β-amyloid deposition upon intracerebral injection in young, pre-depositing APP23 transgenic mice. After ultracentrifugation of the brain extract, <0.05% of the Aβ remained in the supernatant fraction, and these soluble Aβ species were largely PK sensitive. However, upon intracerebral injection, this soluble fraction accounted for up to 30% of the β-amyloid induction observed with the unfractionated extract. Fragmentation of the Aβ seeds by extended sonication increased the seeding capacity of the brain extract. In summary, these results suggest that multiple Aβ assemblies, with various PK sensitivities, are capable of inducing β-amyloid aggregation in vivo. The finding that small and soluble Aβ seeds are potent inducers of cerebral β-amyloidosis raises the possibility that such seeds may mediate the spread of β-amyloidosis in the brain. If they can be identified in vivo, soluble Aβ seeds in bodily fluids also could serve as early biomarkers for cerebral β-amyloidogenesis and eventually Alzheimer's disease.
Exceptionally accurate genome reference sequences have proven to be of great value to microbial researchers. Thus, to date, about 1800 bacterial genome assemblies have been "finished" at great ...expense with the aid of manual laboratory and computational processes that typically iterate over a period of months or even years. By applying a new laboratory design and new assembly algorithm to 16 samples, we demonstrate that assemblies exceeding finished quality can be obtained from whole-genome shotgun data and automated computation. Cost and time requirements are thus dramatically reduced.
Cellular immunity is critical for controlling intracellular pathogens, but individual cellular dynamics and cell-cell cooperativity in evolving human immune responses remain poorly understood. ...Single-cell RNA-sequencing (scRNA-seq) represents a powerful tool for dissecting complex multicellular behaviors in health and disease
and nominating testable therapeutic targets
. Its application to longitudinal samples could afford an opportunity to uncover cellular factors associated with the evolution of disease progression without potentially confounding inter-individual variability
. Here, we present an experimental and computational methodology that uses scRNA-seq to characterize dynamic cellular programs and their molecular drivers, and apply it to HIV infection. By performing scRNA-seq on peripheral blood mononuclear cells from four untreated individuals before and longitudinally during acute infection
, we were powered within each to discover gene response modules that vary by time and cell subset. Beyond previously unappreciated individual- and cell-type-specific interferon-stimulated gene upregulation, we describe temporally aligned gene expression responses obscured in bulk analyses, including those involved in proinflammatory T cell differentiation, prolonged monocyte major histocompatibility complex II upregulation and persistent natural killer (NK) cell cytolytic killing. We further identify response features arising in the first weeks of infection, for example proliferating natural killer cells, which potentially may associate with future viral control. Overall, our approach provides a unified framework for characterizing multiple dynamic cellular responses and their coordination.
Biomphalaria snails are instrumental in transmission of the human blood fluke Schistosoma mansoni. With the World Health Organization's goal to eliminate schistosomiasis as a global health problem by ...2025, there is now renewed emphasis on snail control. Here, we characterize the genome of Biomphalaria glabrata, a lophotrochozoan protostome, and provide timely and important information on snail biology. We describe aspects of phero-perception, stress responses, immune function and regulation of gene expression that support the persistence of B. glabrata in the field and may define this species as a suitable snail host for S. mansoni. We identify several potential targets for developing novel control measures aimed at reducing snail-mediated transmission of schistosomiasis.
Smoldering myeloma (SMM) is associated with a high-risk of progression to myeloma (MM). We report the results of a study of 82 patients with both targeted sequencing that included a capture of the ...immunoglobulin and MYC regions. By comparing these results to newly diagnosed myeloma (MM) we show fewer NRAS and FAM46C mutations together with fewer adverse translocations, del(1p), del(14q), del(16q), and del(17p) in SMM consistent with their role as drivers of the transition to MM. KRAS mutations are associated with a shorter time to progression (HR 3.5 (1.5-8.1), p = 0.001). In an analysis of change in clonal structure over time we studied 53 samples from nine patients at multiple time points. Branching evolutionary patterns, novel mutations, biallelic hits in crucial tumour suppressor genes, and segmental copy number changes are key mechanisms underlying the transition to MM, which can precede progression and be used to guide early intervention strategies.
This study reviews contemporary management and follow-up of pediatric ovarian torsion.
This is a retrospective series of patients from birth to 19 years undergoing operative management of ovarian ...torsion from 2012 to 2016.
We studied 43 girls who underwent 51 operations for ovarian torsion. The median age was 8.3 years. Ultrasound was utilized for diagnosis in 24/29 patients (83%) evaluated in a children's hospital. In contrast, computed tomography was used initially in 7 cases (50%) in children imaged at non–children's hospitals before transfer. Initial operation for ovarian torsion was completed laparoscopically in 38 (88%). Overall, ovarian preservation was performed in 37 (86%) patients, while 6 (13%) underwent oophorectomy. Indications for oophorectomy included 5 infants with in utero torsion and an 18-year-old with a suspected malignancy. In girls with acute ovarian torsion, the oophorectomy rate was reduced to 2%. Postoperatively, 1 patient developed a small bowel obstruction requiring operation after laparoscopic ovarian detorsion. Recurrent torsion occurred in 3 patients (7%). In total, 34 patients underwent postoperative ovarian imaging. A total of 25 (74%) had follicles visualized in the previously torsed ovary.
Ovarian-sparing operations for acute torsion are safe and result in ovarian salvage and preservation of follicular development in more than 70% of children and adolescents.
T cell activation, a key early event in the adaptive immune response, is subject to elaborate transcriptional control. In the present study, we examined how the activities of eight major ...transcription factor (TF) families are integrated to shape the epigenome of naive and activated CD4 and CD8 T cells. By leveraging extensive polymorphisms in evolutionarily divergent mice, we identified the 'heavy lifters' positively influencing chromatin accessibility. Members of Ets, Runx and TCF/Lef TF families occupied the vast majority of accessible chromatin regions, acting as 'housekeepers', 'universal amplifiers' and 'placeholders', respectively, at sites that maintained or gained accessibility upon T cell activation. In addition, a small subset of strongly induced immune response genes displayed a noncanonical TF recruitment pattern. Our study provides a key resource and foundation for the understanding of transcriptional and epigenetic regulation in T cells and offers a new perspective on the hierarchical interactions between critical TFs.
Deep sequencing technologies have the potential to transform the study of highly variable viral pathogens by providing a rapid and cost-effective approach to sensitively characterize rapidly evolving ...viral quasispecies. Here, we report on a high-throughput whole HIV-1 genome deep sequencing platform that combines 454 pyrosequencing with novel assembly and variant detection algorithms. In one subject we combined these genetic data with detailed immunological analyses to comprehensively evaluate viral evolution and immune escape during the acute phase of HIV-1 infection. The majority of early, low frequency mutations represented viral adaptation to host CD8+ T cell responses, evidence of strong immune selection pressure occurring during the early decline from peak viremia. CD8+ T cell responses capable of recognizing these low frequency escape variants coincided with the selection and evolution of more effective secondary HLA-anchor escape mutations. Frequent, and in some cases rapid, reversion of transmitted mutations was also observed across the viral genome. When located within restricted CD8 epitopes these low frequency reverting mutations were sufficient to prime de novo responses to these epitopes, again illustrating the capacity of the immune response to recognize and respond to low frequency variants. More importantly, rapid viral escape from the most immunodominant CD8+ T cell responses coincided with plateauing of the initial viral load decline in this subject, suggestive of a potential link between maintenance of effective, dominant CD8 responses and the degree of early viremia reduction. We conclude that the early control of HIV-1 replication by immunodominant CD8+ T cell responses may be substantially influenced by rapid, low frequency viral adaptations not detected by conventional sequencing approaches, which warrants further investigation. These data support the critical need for vaccine-induced CD8+ T cell responses to target more highly constrained regions of the virus in order to ensure the maintenance of immunodominant CD8 responses and the sustained decline of early viremia.