Menu
April 21, 2020

Accurate circular consensus long-read sequencing improves variant detection and assembly of a human genome.

The DNA sequencing technologies in use today produce either highly accurate short reads or less-accurate long reads. We report the optimization of circular consensus sequencing (CCS) to improve the accuracy of single-molecule real-time (SMRT) sequencing (PacBio) and generate highly accurate (99.8%) long high-fidelity (HiFi) reads with an average length of 13.5?kilobases (kb). We applied our approach to sequence the well-characterized human HG002/NA24385 genome and obtained precision and recall rates of at least 99.91% for single-nucleotide variants (SNVs), 95.98% for insertions and deletions <50 bp (indels) and 95.99% for structural variants. Our CCS method matches or exceeds the ability of short-read sequencing to detect small variants and structural variants. We estimate that 2,434 discordances are correctable mistakes in the 'genome in a bottle' (GIAB) benchmark set. Nearly all (99.64%) variants can be phased into haplotypes, further improving variant detection. De novo genome assembly using CCS reads alone produced a contiguous and accurate genome with a contig N50 of >15?megabases (Mb) and concordance of 99.997%, substantially outperforming assembly with less-accurate long reads.


April 21, 2020

Petunia-and Arabidopsis-Specific Root Microbiota Responses to Phosphate Supplementation

Phosphorus (P) is a limiting element for plant growth. Several root microbes, including arbuscular mycorrhizal fungi (AMF), have the capacity to improve plant nutrition and their abundance is known to depend on P fertility. However, how complex root-associated bacterial and fungal communities respond to various levels of P supplementation remains ill-defined. Here we investigated the responses of the root-associated bacteria and fungi to varying levels of P supply using 16S rRNA gene and internal transcribed spacer amplicon sequencing. We grew Petunia, which forms symbiosis with AMF, and the nonmycorrhizal model species Arabidopsis as a control in a soil that is limiting in plant-available P and we then supplemented the plants with complete fertilizer solutions that varied only in their phosphate concentrations. We searched for microbes, whose abundances varied by P fertilization, tested whether a core microbiota responding to the P treatments could be identified and asked whether bacterial and fungal co-occurrence patterns change in response to the varying P levels. Root microbiota composition varied substantially in response to the varying P application. A core microbiota was not identified as different bacterial and fungal groups responded to low-P conditions in Arabidopsis and Petunia. Microbes with P-dependent abundance patterns included Mortierellomycotina in Arabidopsis, while in Petunia, they included AMF and their symbiotic endobacteria. Of note, the P-dependent root colonization by AMF was reliably quantified by sequencing. The fact that the root microbiotas of the two plant species responded differently to low-P conditions suggests that plant species specificity would need to be considered for the eventual development of microbial products that improve plant P nutrition.


April 21, 2020

Real time monitoring of Aeromonas salmonicida evolution in response to successive antibiotic therapies in a commercial fish farm.

Our ability to predict evolutionary trajectories of pathogens in response to antibiotic pressure is one of the promising leverage to fight against the present antibiotic resistance worldwide crisis. Yet, few studies tackled this question in situ at the outbreak level, due to the difficulty to link a given pathogenic clone evolution with its precise antibiotic exposure over time. In this study, we monitored the real-time evolution of an Aeromonas salmonicida clone in response to successive antibiotic and vaccine therapies in a commercial fish farm. The clone was responsible for a four-year outbreak of furunculosis within a Recirculating Aquaculture System Salmo salar farm in China, and we reconstructed the precise tempo of mobile genetic elements (MGEs) acquisition events during this period. The resistance profile provided by the acquired MGEs closely mirrored the antibiotics used to treat the outbreak, and we evidenced that two subclonal groups developed similar resistances although unrelated MGE acquisitions. Finally, we also demonstrated the efficiency of vaccination in outbreak management and its positive effect on antibiotic resistance prevalence. Our study provides unprecedented knowledge critical to understand evolutionary trajectories of resistant pathogens outside the laboratory. © 2019 Society for Applied Microbiology and John Wiley & Sons Ltd.


April 21, 2020

Mutation of a bHLH transcription factor allowed almond domestication.

Wild almond species accumulate the bitter and toxic cyanogenic diglucoside amygdalin. Almond domestication was enabled by the selection of genotypes harboring sweet kernels. We report the completion of the almond reference genome. Map-based cloning using an F1 population segregating for kernel taste led to the identification of a 46-kilobase gene cluster encoding five basic helix-loop-helix transcription factors, bHLH1 to bHLH5. Functional characterization demonstrated that bHLH2 controls transcription of the P450 monooxygenase-encoding genes PdCYP79D16 and PdCYP71AN24, which are involved in the amygdalin biosynthetic pathway. A nonsynonymous point mutation (Leu to Phe) in the dimerization domain of bHLH2 prevents transcription of the two cytochrome P450 genes, resulting in the sweet kernel trait. Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.


April 21, 2020

Combining orthogonal CRISPR and CRISPRi systems for genome engineering and metabolic pathway modulation in Escherichia coli.

CRISPR utilizing Cas9 from Streptococcus pyogenes (SpCas9) and CRISPR interference (CRISPRi) employing catalytically inactive SpCas9 (SpdCas9) have gained popularity for Escherichia coli engineering. To integrate the SpdCas9-based CRISPRi module using CRISPR while avoiding mutual interference between SpCas9/SpdCas9 and their cognate single-guide RNA (sgRNA), this study aimed at exploring an alternative Cas nuclease orthogonal to SpCas9. We compared several Cas9 variants from different microorganisms such as Staphylococcus aureus (SaCas9) and Streptococcus thermophilius CRISPR1 (St1Cas9) as well as Cas12a derived from Francisella novicida (FnCas12a). At the commonly used E. coli model genes  LacZ, we found that SaCas9 and St1Cas9 induced DNA cleavage more effectively than FnCas12a. Both St1Cas9 and SaCas9 were orthogonal to SpCas9 and the induced DNA cleavage promoted the integration of heterologous DNA of up to 10?kb, at which size St1Cas9 was superior to SaCas9 in recombination frequency/accuracy. We harnessed the St1Cas9 system to integrate SpdCas9 and sgRNA arrays for constitutive knockdown of three genes, knock-in pyc and knockout adhE, without compromising the CRISPRi knockdown efficiency. The combination of orthogonal CRISPR/CRISPRi for metabolic engineering enhanced succinate production while inhibiting byproduct formation and may pave a new avenue to E. coli engineering. © 2019 Wiley Periodicals, Inc.


April 21, 2020

Insights into the evolution and drug susceptibility of Babesia duncani from the sequence of its mitochondrial and apicoplast genomes.

Babesia microti and Babesia duncani are the main causative agents of human babesiosis in the United States. While significant knowledge about B. microti has been gained over the past few years, nothing is known about B. duncani biology, pathogenesis, mode of transmission or sensitivity to currently recommended therapies. Studies in immunocompetent wild type mice and hamsters have shown that unlike B. microti, infection with B. duncani results in severe pathology and ultimately death. The parasite factors involved in B. duncani virulence remain unknown. Here we report the first known completed sequence and annotation of the apicoplast and mitochondrial genomes of B. duncani. We found that the apicoplast genome of this parasite consists of a 34?kb monocistronic circular molecule encoding functions that are important for apicoplast gene transcription as well as translation and maturation of the organelle’s proteins. The mitochondrial genome of B. duncani consists of a 5.9?kb monocistronic linear molecule with two inverted repeats of 48?bp at both ends. Using the conserved cytochrome b (Cytb) and cytochrome c oxidase subunit I (coxI) proteins encoded by the mitochondrial genome, phylogenetic analysis revealed that B. duncani defines a new lineage among apicomplexan parasites distinct from B. microti, Babesia bovis, Theileria spp. and Plasmodium spp. Annotation of the apicoplast and mitochondrial genomes of B. duncani identified targets for development of effective therapies. Our studies set the stage for evaluation of the efficacy of these drugs alone or in combination against B. duncani in culture as well as in animal models.Copyright © 2018 Australian Society for Parasitology. Published by Elsevier Ltd. All rights reserved.


April 21, 2020

Recompleting the Caenorhabditis elegans genome.

Caenorhabditis elegans was the first multicellular eukaryotic genome sequenced to apparent completion. Although this assembly employed a standard C. elegans strain (N2), it used sequence data from several laboratories, with DNA propagated in bacteria and yeast. Thus, the N2 assembly has many differences from any C. elegans available today. To provide a more accurate C. elegans genome, we performed long-read assembly of VC2010, a modern strain derived from N2. Our VC2010 assembly has 99.98% identity to N2 but with an additional 1.8 Mb including tandem repeat expansions and genome duplications. For 116 structural discrepancies between N2 and VC2010, 97 structures matching VC2010 (84%) were also found in two outgroup strains, implying deficiencies in N2. Over 98% of N2 genes encoded unchanged products in VC2010; moreover, we predicted =53 new genes in VC2010. The recompleted genome of C. elegans should be a valuable resource for genetics, genomics, and systems biology. © 2019 Yoshimura et al.; Published by Cold Spring Harbor Laboratory Press.


April 21, 2020

Competition between mobile genetic elements drives optimization of a phage-encoded CRISPR-Cas system: insights from a natural arms race.

CRISPR-Cas systems function as adaptive immune systems by acquiring nucleotide sequences called spacers that mediate sequence-specific defence against competitors. Uniquely, the phage ICP1 encodes a Type I-F CRISPR-Cas system that is deployed to target and overcome PLE, a mobile genetic element with anti-phage activity in Vibrio cholerae. Here, we exploit the arms race between ICP1 and PLE to examine spacer acquisition and interference under laboratory conditions to reconcile findings from wild populations. Natural ICP1 isolates encode multiple spacers directed against PLE, but we find that single spacers do not interfere equally with PLE mobilization. High-throughput sequencing to assay spacer acquisition reveals that ICP1 can also acquire spacers that target the V. cholerae chromosome. We find that targeting the V. cholerae chromosome proximal to PLE is sufficient to block PLE and is dependent on Cas2-3 helicase activity. We propose a model in which indirect chromosomal spacers are able to circumvent PLE by Cas2-3-mediated processive degradation of the V. cholerae chromosome before PLE mobilization. Generally, laboratory-acquired spacers are much more diverse than the subset of spacers maintained by ICP1 in nature, showing how evolutionary pressures can constrain CRISPR-Cas targeting in ways that are often not appreciated through in vitro analyses. This article is part of a discussion meeting issue ‘The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems’.


April 21, 2020

Plasmid analysis of Escherichia coli isolates from South Korea co-producing NDM-5 and OXA-181 carbapenemases.

Recently, Escherichia coli isolates co-producing New Delhi metallo-ß-lactamase (NDM)-5 and oxacillinase (OXA)-181 were identified in a tertiary-care hospital of South Korea. Isolate CC1702-1 was collected from urine in January 2017 and isolate CC1706-1 was recovered from a transtracheal aspirate of a hospitalized patient in May 2017. Carbapenemase genes were identified by multiplex PCR and sequencing, and whole genome sequencing was performed subsequently using the PacBio RSII system. Both E. coli isolates belonged to the same clone (ST410) and were resistant to all ß-lactams including carbapenems. We obtained whole plasmid sequences of the isolates: pCC1702-NDM-5 from CC1702-1 and pCC1706-NDM-5 and pCC1706-OXA-181 from CC1706-1. The two E. coli isolates belonged to the same clone (ST410) and they were completely resistant to all ß-lactams, as well as carbapenems. Two blaNDM-5-harboring plasmids belonged to the same incompatibility group, IncFIA/B, and consisted of 79,613?bp and 111,890?bp with 87 and 130 coding sequences, respectively. The genetic structures of the two blaNDM-5-bearing plasmids, which were distinct from the blaNDM-5-bearing plasmids from the Klebsiella pneumoniae isolates previously transmitted from the United Arab Emirates (UAE) to South Korea, differed from each other. While pCC1702-NDM-5 showed high degree of identity with the plasmid from a multidrug-resistant isolate of Citrobacter fruendii P5571 found in China, pCC1706-NDM-5 was very similar to the plasmid from a multidrug-resistant isolate of E. coli AMA1176 found in Denmark. pCC1706-OXA-181, which was a 51?kb, self-transmissible IncX3 plasmid, was identical to the E. coli plasmids pAMA1167-OXA-181 from Denmark and pOXA-181-WCHEC14828 from China. Plasmids harboring blaNDM-5 in E. coli isolates might not be transferred from K. pneumoniae isolates co-producing NDM-5 and OXA-181. They probably originated from multiple sources.Copyright © 2019 Elsevier Inc. All rights reserved.


April 21, 2020

Development of CRISPR-Cas systems for genome editing and beyond

The development of clustered regularly interspaced short-palindromic repeat (CRISPR)-Cas systems for genome editing has transformed the way life science research is conducted and holds enormous potential for the treatment of disease as well as for many aspects of biotech- nology. Here, I provide a personal perspective on the development of CRISPR-Cas9 for genome editing within the broader context of the field and discuss our work to discover novel Cas effectors and develop them into additional molecular tools. The initial demonstra- tion of Cas9-mediated genome editing launched the development of many other technologies, enabled new lines of biological inquiry, and motivated a deeper examination of natural CRISPR-Cas systems, including the discovery of new types of CRISPR-Cas systems. These new discoveries in turn spurred further technological developments. I review these exciting discoveries and technologies as well as provide an overview of the broad array of applications of these technologies in basic research and in the improvement of human health. It is clear that we are only just beginning to unravel the potential within microbial diversity, and it is quite likely that we will continue to discover other exciting phenomena, some of which it may be possible to repurpose as molecular technologies. The transformation of mysterious natural phenomena to powerful tools, however, takes a collective effort to discover, characterize, and engineer them, and it has been a privilege to join the numerous researchers who have contributed to this transformation of CRISPR-Cas systems.


April 21, 2020

Genome and plasmid diversity of Extended-Spectrum ß-Lactamase-producing Escherichia coli ST131 – tracking phylogenetic trajectories with Bayesian inference.

Clonal lineages of ESBL (Extended-Spectrum ß-Lactamase)-producing E. coli belonging to sequence type 131 (ST131) have disseminated globally during the last 30 years, leading to an increased prevalence of resistance to fluoroquinolones and extended-spectrum cephalosporins in clinical isolates of E. coli. We aimed to study if Swedish ESBL-producing ST131 isolates originated from single or multiple introductions to the population by assessing the amount of genetic variation, on chromosomal and plasmid level, between Swedish and international E. coli ST131. Bayesian inference of Swedish E. coli ST131 isolates (n?=?29), sequenced using PacBio RSII, together with an international ST131 dataset showed that the Swedish isolates were part of the international ST131 A, C1 and C2 clades. Highly conserved plasmids were identified in three clusters although they were separated by several years, which indicates a strong co-evolution between some ST131 lineages and specific plasmids. In conclusion, the tight clonal relationship observed within the ST131 clades, together with highly conserved plasmids, challenges investigation of strain transmission events. A combination of few SNPs on a genome-wide scale and an epidemiological temporospatial link, are needed to track the spread of the ST131 subclones.


April 21, 2020

SMRT sequencing reveals differential patterns of methylation in two O111:H- STEC isolates from a hemolytic uremic syndrome outbreak in Australia.

In 1995 a severe haemolytic-uremic syndrome (HUS) outbreak in Adelaide occurred. A recent genomic analysis of Shiga toxigenic Escherichia coli (STEC) O111:H- strains 95JB1 and 95NR1 from this outbreak found that the more virulent isolate, 95NR1, harboured two additional copies of the Shiga toxin 2 (Stx2) genes encoded within prophage regions. The structure of the Stx2-converting prophages could not be fully resolved using short-read sequence data alone and it was not clear if there were other genomic differences between 95JB1 and 95NR1. In this study we have used Pacific Biosciences (PacBio) single molecule real-time (SMRT) sequencing to characterise the genome and methylome of 95JB1 and 95NR1. We completely resolved the structure of all prophages including two, tandemly inserted, Stx2-converting prophages in 95NR1 that were absent from 95JB1. Furthermore we defined all insertion sequences and found an additional IS1203 element in the chromosome of 95JB1. Our analysis of the methylome of 95NR1 and 95JB1 identified hemi-methylation of a novel motif (5′-CTGCm6AG-3′) in more than 4000 sites in the 95NR1 genome. These sites were entirely unmethylated in the 95JB1 genome, and included at least 177 potential promoter regions that could contribute to regulatory differences between the strains. IS1203 mediated deactivation of a novel type IIG methyltransferase in 95JB1 is the likely cause of the observed differential patterns of methylation between 95NR1 and 95JB1. This study demonstrates the capability of PacBio SMRT sequencing to resolve complex prophage regions and reveal the genetic and epigenetic heterogeneity within a clonal population of bacteria.


April 21, 2020

Extended insight into the Mycobacterium chelonae-abscessus complex through whole genome sequencing of Mycobacterium salmoniphilum outbreak and Mycobacterium salmoniphilum-like strains.

Members of the Mycobacterium chelonae-abscessus complex (MCAC) are close to the mycobacterial ancestor and includes both human, animal and fish pathogens. We present the genomes of 14 members of this complex: the complete genomes of Mycobacterium salmoniphilum and Mycobacterium chelonae type strains, seven M. salmoniphilum isolates, and five M. salmoniphilum-like strains including strains isolated during an outbreak in an animal facility at Uppsala University. Average nucleotide identity (ANI) analysis and core gene phylogeny revealed that the M. salmoniphilum-like strains are variants of the human pathogen Mycobacterium franklinii and phylogenetically close to Mycobacterium abscessus. Our data further suggested that M. salmoniphilum separates into three branches named group I, II and III with the M. salmoniphilum type strain belonging to group II. Among predicted virulence factors, the presence of phospholipase C (plcC), which is a major virulence factor that makes M. abscessus highly cytotoxic to mouse macrophages, and that M. franklinii originally was isolated from infected humans make it plausible that the outbreak in the animal facility was caused by a M. salmoniphilum-like strain. Interestingly, M. salmoniphilum-like was isolated from tap water suggesting that it can be present in the environment. Moreover, we predicted the presence of mutational hotspots in the M. salmoniphilum isolates and 26% of these hotspots overlap with genes categorized as having roles in virulence, disease and defense. We also provide data about key genes involved in transcription and translation such as sigma factor, ribosomal protein and tRNA genes.


April 21, 2020

An African Salmonella Typhimurium ST313 sublineage with extensive drug-resistance and signatures of host adaptation.

Bloodstream infections by Salmonella enterica serovar Typhimurium constitute a major health burden in sub-Saharan Africa (SSA). These invasive non-typhoidal (iNTS) infections are dominated by isolates of the antibiotic resistance-associated sequence type (ST) 313. Here, we report emergence of ST313 sublineage II.1 in the Democratic Republic of the Congo. Sublineage II.1 exhibits extensive drug resistance, involving a combination of multidrug resistance, extended spectrum ß-lactamase production and azithromycin resistance. ST313 lineage II.1 isolates harbour an IncHI2 plasmid we name pSTm-ST313-II.1, with one isolate also exhibiting decreased ciprofloxacin susceptibility. Whole genome sequencing reveals that ST313 II.1 isolates have accumulated genetic signatures potentially associated with altered pathogenicity and host adaptation, related to changes observed in biofilm formation and metabolic capacity. Sublineage II.1 emerged at the beginning of the 21st century and is involved in on-going outbreaks. Our data provide evidence of further evolution within the ST313 clade associated with iNTS in SSA.


April 21, 2020

Genome-wide mutational biases fuel transcriptional diversity in the Mycobacterium tuberculosis complex.

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.