Assembly Of Short Reads Research Articles

Rapid, robust plasmid verification by de novo assembly of short sequencing reads.

Plasmids are a foundational tool for basic and applied research across all subfields of biology. Increasingly, researchers in synthetic biology are relying on and developing massive libraries of plasmids as vectors for directed evolution, combinatorial gene circuit tests, and for CRISPR multiplexing. Verification of plasmid sequences following synthesis is a crucial quality control step that creates a bottleneck in plasmid fabrication workflows. Crucially, researchers often elect to forego the cumbersome verification step, potentially leading to reproducibility and—depending on the application—security issues. In order to facilitate plasmid verification to improve the quality and reproducibility of life science research, we developed a fast, simple, and open source pipeline for assembly and verification of plasmid sequences from Illumina reads. We demonstrate that our pipeline, which relies on de novo assembly, can also be used to detect contaminating sequences in plasmid samples. In addition to presenting our pipeline, we discuss the role for verification and quality control in the increasingly complex life science workflows ushered in by synthetic biology.

Symbiosis islands of Loteae-nodulating Mesorhizobium comprise three radiating lineages with concordant nod gene complements and nodulation host-range groupings.

Mesorhizobium is a genus of soil bacteria, some isolates of which form an endosymbiotic relationship with diverse legumes of the Loteae tribe. The symbiotic genes of these mesorhizobia are generally carried on integrative and conjugative elements termed symbiosis islands (ICESyms). Mesorhizobium strains that nodulate Lotus spp. have been divided into host-range groupings. Group I (GI) strains nodulate L. corniculatus and L. japonicus ecotype Gifu, while group II (GII) strains have a broader host range, which includes L. pedunculatus. To identify the basis of this extended host range, and better understand Mesorhizobium and ICESym genomics, the genomes of eight Mesorhizobium strains were completed using hybrid long- and short-read assembly. Bioinformatic comparison with previously sequenced mesorhizobia genomes indicated host range was not predicted by Mesorhizobium genospecies but rather by the evolutionary relationship between ICESym symbiotic regions. Three radiating lineages of Loteae ICESyms were identified on this basis, which correlate with Lotus spp. host-range grouping and have lineage-specific nod gene complements. Pangenomic analysis of the completed GI and GII ICESyms identified 155 core genes (on average 30.1 % of a given ICESym). Individual GI or GII ICESyms carried diverse accessory genes with an average of 34.6 % of genes unique to a given ICESym. Identification and comparative analysis of NodD symbiotic regulatory motifs – nod boxes – identified 21 branches across the NodD regulons. Four of these branches were associated with seven genes unique to the five GII ICESyms. The nod boxes preceding the host-range gene nodZ in GI and GII ICESyms were disparate, suggesting regulation of nodZ may differ between GI and GII ICESyms. The broad host-range determinant(s) of GII ICESyms that confer nodulation of L. pedunculatus are likely present amongst the 53 GII-unique genes identified.

Identification of plasmids by PCR based replicon typing in bacteremic Klebsiella pneumoniae

BackgroundKlebsiella pneumoniae is a notorious pathogen with plasmid mediated resistance to all classes of antibiotics. It is important to determine the plasmid profile coding for resistance genes. Plasmid profile varies among geographical regions and tracking the types helps in determining the MDR and XDR K. pneumoniae spread especially in hospital setting. Aim of the present study was to determine the plasmid profile and types among bacteraemic K. pneumoniae. Materials and methodsNinety consecutive K. pneumoniae collected over a period of three months from blood cultures were characterised by PCR for plasmid profile. Inc plasmid types were determined by PCR based replicon typing (PBRT) and carbapenemases were determined by multiplex PCR. For a subset of isolates hybrid assemblies were developed by sequencing with Ion Torrent and MinIon. ResultsOverall, PBRT showed 29% of isolates carried four plasmids including IncHI1B, IncFIA, IncFII(K) and IncR. The most common type of plasmid was IncHI1B (93%) followed by IncFIIK (89%) and IncR (82%). IncFIA was predominant among carbapenem resistant isolates. Almost all plasmids identified in K. pneumoniae were AMR plasmids, except two isolates which had virulence plasmids. IncX3 plasmid observed in this study was previously reported to be self-disseminating. Furthermore, the hybrid genome sequencing revealed complete structural arrangements of plasmids, which would be missed in short-read sequencing. NDM and OXA48-like were co-produced in 59% of the carbapenem resistant isolates. BlaOXA-232 was present on ColKP3; aac(6’)-lb3 and rmtF on IncFIB. ConclusionDiverse plasmid profile among the successive K. pneumoniae isolates indicates the transfer of resistance genes through different types of plasmids. IncHI1B, IncFIA, IncFIIK and IncR were the prevalent plasmid types. Hybrid assembly revealed blaOXA-232 was present on ColKP3 unlike global reports of IncL/M. Hybrid assemblies provide better plasmid structure that long and short read assemblies. There was no significant association of β-lactamases with specific Inc groups in this study.

Whole-Genome Sequence of an Orf Virus Isolate Derived from a Cell Culture Infected with Contagious Ecthyma Vaccine.

This is a draft genome of an orf virus (ORFV) vaccine strain assembled via long- and short-read hybrid assembly. ORFV is a zoonotic pathogen that affects sheep and goats. The genome of the virus contained in the vaccine was found to have high similarity (98%) to those of other published strains.

A high-contiguity Brassica nigra genome localizes active centromeres and defines the ancestral Brassica genome

It is only recently, with the advent of long-read sequencing technologies, that we are beginning to uncover previously uncharted regions of complex and inherently recursive plant genomes. To comprehensively study and exploit the genome of the neglected oilseed Brassica nigra, we generated two high-quality nanopore de novo genome assemblies. The N50 contig lengths for the two assemblies were 17.1 Mb (12 contigs), one of the best among 324 sequenced plant genomes, and 0.29 Mb (424 contigs), respectively, reflecting recent improvements in the technology. Comparison with a de novo short-read assembly corroborated genome integrity and quantified sequence-related error rates (0.2%). The contiguity and coverage allowed unprecedented access to low-complexity regions of the genome. Pericentromeric regions and coincidence of hypomethylation enabled localization of active centromeres and identified centromere-associated ALE family retro-elements that appear to have proliferated through relatively recent nested transposition events (<1 Ma). Genomic distances calculated based on synteny relationships were used to define a post-triplication Brassica-specific ancestral genome, and to calculate the extensive rearrangements that define the evolutionary distance separating B. nigra from its diploid relatives.

Completion of draft bacterial genomes by long-read sequencing of synthetic genomic pools

BackgroundIllumina technology currently dominates bacterial genomics due to its high read accuracy and low sequencing cost. However, the incompleteness of draft genomes generated by Illumina reads limits their application in comprehensive genomics analyses. Alternatively, hybrid assembly using both Illumina short reads and long reads generated by single molecule sequencing technologies can enable assembly of complete bacterial genomes, yet the high per-genome cost of long-read sequencing limits the widespread use of this approach in bacterial genomics. Here we developed a protocol for hybrid assembly of complete bacterial genomes using miniaturized multiplexed Illumina sequencing and non-barcoded PacBio sequencing of a synthetic genomic pool (SGP), thus significantly decreasing the overall per-genome cost of sequencing.ResultsWe evaluated the performance of SGP hybrid assembly on the genomes of 20 bacterial isolates with different genome sizes, a wide range of GC contents, and varying levels of phylogenetic relatedness. By improving the contiguity of Illumina assemblies, SGP hybrid assembly generated 17 complete and 3 nearly complete bacterial genomes. Increased contiguity of SGP hybrid assemblies resulted in considerable improvement in gene prediction and annotation. In addition, SGP hybrid assembly was able to resolve repeat elements and identify intragenomic heterogeneities, e.g. different copies of 16S rRNA genes, that would otherwise go undetected by short-read-only assembly. Comprehensive comparison of SGP hybrid assemblies with those generated using multiplexed PacBio long reads (long-read-only assembly) also revealed the relative advantage of SGP hybrid assembly in terms of assembly quality. In particular, we observed that SGP hybrid assemblies were completely devoid of both small (i.e. single base substitutions) and large assembly errors. Finally, we show the ability of SGP hybrid assembly to differentiate genomes of closely related bacterial isolates, suggesting its potential application in comparative genomics and pangenome analysis.ConclusionOur results indicate the superiority of SGP hybrid assembly over both short-read and long-read assemblies with respect to completeness, contiguity, accuracy, and recovery of small replicons. By lowering the per-genome cost of sequencing, our parallel sequencing and hybrid assembly pipeline could serve as a cost effective and high throughput approach for completing high-quality bacterial genomes.

AbGRI4, a novel antibiotic resistance island in multiply antibiotic-resistant Acinetobacter baumannii clinical isolates.

ObjectivesTo investigate the genomic context of a novel resistance island (RI) in multiply antibiotic-resistant Acinetobacter baumannii clinical isolates and global isolates.MethodsUsing a combination of long and short reads generated from the Oxford Nanopore and Illumina platforms, contiguous chromosomes and plasmid sequences were determined. BLAST-based analysis was used to identify the RI insertion target.ResultsGenomes of four multiply antibiotic-resistant A. baumannii clinical strains, from a US hospital system, belonging to prevalent MLST ST2 (Pasteur scheme) and ST281 (Oxford scheme) clade F isolates were sequenced to completion. A class 1 integron carrying aadB (tobramycin resistance) and aadA2 (streptomycin/spectinomycin resistance) was identified. The class 1 integron was 6.8 kb, bounded by IS26 at both ends, and embedded in a new target location between an α/β-hydrolase and a reductase. Due to its novel insertion site and unique RI composition, we suggest naming this novel RI AbGRI4. Molecular analysis of global A. baumannii isolates identified multiple AbGRI4 RI variants in non-ST2 clonal lineages, including variations in the resistance gene cassettes, integron backbone and insertion breakpoints at the hydrolase gene.ConclusionsA novel RI insertion target harbouring a class 1 integron was identified in a subgroup of ST2/ST281 clinical isolates. Variants of the RI suggested evolution and horizontal transfer of the RI across clonal lineages. Long- and short-read hybrid assembly technology completely resolved the genomic context of IS-bounded RIs, which was not possible using short reads alone.

Draft Genome Sequence Resource for Blumeriella jaapii, the Cherry Leaf Spot Pathogen.

Blumeriella jaapii is the causal agent of cherry leaf spot (CLS), the most important disease of tart cherry in the Midwestern United States. Infection of leaves by B. jaapii leads to premature defoliation, which places trees at heightened risk of winter injury and death. Current management of CLS relies primarily on the application of three important fungicide classes, quinone outside inhibitors, sterol demethylation inhibitors, and succinate dehydrogenase inhibitors. Here, we present the first high-quality genome of B. jaapii through a hybrid assembly of PacBio long reads and Illumina short reads. The assembled draft genome of B. jaapii is 47.4 Mb and consists of 95 contigs with a N50 value of 1.5 Mb. The genomic information of B. jaapii, representing the most complete sequenced genome of the family Dermateaceae (Ascomycota) to date, provides a valuable resource for identifying fungicide resistance mechanisms of this pathogen and expands our knowledge of the phytopathogenic fungi in this family.

The Salmonella enterica Plasmidome as a Reservoir of Antibiotic Resistance.

The emergence of multidrug-resistant bacterial strains worldwide has become a serious problem for public health over recent decades. The increase in antimicrobial resistance has been expanding via plasmids as mobile genetic elements encoding antimicrobial resistance (AMR) genes that are transferred vertically and horizontally. This study focuses on Salmonella enterica, one of the leading foodborne pathogens in industrialized countries. S. enterica is known to carry several plasmids involved not only in virulence but also in AMR. In the current paper, we present an integrated strategy to detect plasmid scaffolds in whole genome sequencing (WGS) assemblies. We developed a two-step procedure to predict plasmids based on i) the presence of essential elements for plasmid replication and mobility, as well as ii) sequence similarity to a reference plasmid. Next, to confirm the accuracy of the prediction in 1750 S. enterica short-read sequencing data, we combined Oxford Nanopore MinION long-read sequencing with Illumina MiSeq short-read sequencing in hybrid assemblies for 84 isolates to evaluate the proportion of plasmid that has been detected. At least one scaffold with an origin of replication (ORI) was predicted in 61.3% of the Salmonella isolates tested. The results indicated that IncFII and IncI1 ORIs were distributed in many S. enterica serotypes and were the most prevalent AMR genes carrier, whereas IncHI2A/IncHI2 and IncA/C2 were more serotype restricted but bore several AMR genes. Comparison between hybrid and short-read assemblies revealed that 81.1% of plasmids were found in the short-read sequencing using our pipeline. Through this process, we established that plasmids are prevalent in S. enterica and we also substantially expand the AMR genes in the resistome of this species.

Platon: identification and characterization of bacterial plasmid contigs in short-read draft assemblies exploiting protein sequence-based replicon distribution scores.

Plasmids are extrachromosomal genetic elements that replicate independently of the chromosome and play a vital role in the environmental adaptation of bacteria. Due to potential mobilization or conjugation capabilities, plasmids are important genetic vehicles for antimicrobial resistance genes and virulence factors with huge and increasing clinical implications. They are therefore subject to large genomic studies within the scientific community worldwide. As a result of rapidly improving next-generation sequencing methods, the quantity of sequenced bacterial genomes is constantly increasing, in turn raising the need for specialized tools to (i) extract plasmid sequences from draft assemblies, (ii) derive their origin and distribution, and (iii) further investigate their genetic repertoire. Recently, several bioinformatic methods and tools have emerged to tackle this issue; however, a combination of high sensitivity and specificity in plasmid sequence identification is rarely achieved in a taxon-independent manner. In addition, many software tools are not appropriate for large high-throughput analyses or cannot be included in existing software pipelines due to their technical design or software implementation. In this study, we investigated differences in the replicon distributions of protein-coding genes on a large scale as a new approach to distinguish plasmid-borne from chromosome-borne contigs. We defined and computed statistical discrimination thresholds for a new metric: the replicon distribution score (RDS), which achieved an accuracy of 96.6 %. The final performance was further improved by the combination of the RDS metric with heuristics exploiting several plasmid-specific higher-level contig characterizations. We implemented this workflow in a new high-throughput taxon-independent bioinformatics software tool called Platon for the recruitment and characterization of plasmid-borne contigs from short-read draft assemblies. Compared to PlasFlow, Platon achieved a higher accuracy (97.5 %) and more balanced predictions (F1=82.6 %) tested on a broad range of bacterial taxa and better or equal performance against the targeted tools PlasmidFinder and PlaScope on sequenced Escherichia coli isolates. Platon is available at: http://platon.computational.bio/.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala).

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Platanus_B: an accurate de novo assembler for bacterial genomes using an iterative error-removal process.

De novo assembly of short DNA reads remains an essential technology, especially for large-scale projects and high-resolution variant analyses in epidemiology. However, the existing tools often lack sufficient accuracy required to compare closely related strains. To facilitate such studies on bacterial genomes, we developed Platanus_B, a de novo assembler that employs iterations of multiple error-removal algorithms. The benchmarks demonstrated the superior accuracy and high contiguity of Platanus_B, in addition to its ability to enhance the hybrid assembly of both short and nanopore long reads. Although the hybrid strategies for short and long reads were effective in achieving near full-length genomes, we found that short-read-only assemblies generated with Platanus_B were sufficient to obtain ≥90% of exact coding sequences in most cases. In addition, while nanopore long-read-only assemblies lacked fine-scale accuracies, inclusion of short reads was effective in improving the accuracies. Platanus_B can, therefore, be used for comprehensive genomic surveillances of bacterial pathogens and high-resolution phylogenomic analyses of a wide range of bacteria.

An improved pig reference genome sequence to enable pig genetics and genomics research.

BackgroundThe domestic pig (Sus scrofa) is important both as a food source and as a biomedical model given its similarity in size, anatomy, physiology, metabolism, pathology, and pharmacology to humans. The draft reference genome (Sscrofa10.2) of a purebred Duroc female pig established using older clone-based sequencing methods was incomplete, and unresolved redundancies, short-range order and orientation errors, and associated misassembled genes limited its utility.ResultsWe present 2 annotated highly contiguous chromosome-level genome assemblies created with more recent long-read technologies and a whole-genome shotgun strategy, 1 for the same Duroc female (Sscrofa11.1) and 1 for an outbred, composite-breed male (USMARCv1.0). Both assemblies are of substantially higher (>90-fold) continuity and accuracy than Sscrofa10.2.ConclusionsThese highly contiguous assemblies plus annotation of a further 11 short-read assemblies provide an unprecedented view of the genetic make-up of this important agricultural and biomedical model species. We propose that the improved Duroc assembly (Sscrofa11.1) become the reference genome for genomic research in pigs.

NtJoin: Fast and lightweight assembly-guided scaffolding using minimizer graphs.

SummaryThe ability to generate high-quality genome sequences is cornerstone to modern biological research. Even with recent advancements in sequencing technologies, many genome assemblies are still not achieving reference-grade. Here, we introduce ntJoin, a tool that leverages structural synteny between a draft assembly and reference sequence(s) to contiguate and correct the former with respect to the latter. Instead of alignments, ntJoin uses a lightweight mapping approach based on a graph data structure generated from ordered minimizer sketches. The tool can be used in a variety of different applications, including improving a draft assembly with a reference-grade genome, a short-read assembly with a draft long-read assembly and a draft assembly with an assembly from a closely related species. When scaffolding a human short-read assembly using the reference human genome or a long-read assembly, ntJoin improves the NGA50 length 23- and 13-fold, respectively, in under 13 m, using <11 GB of RAM. Compared to existing reference-guided scaffolders, ntJoin generates highly contiguous assemblies faster and using less memory.Availability and implementationntJoin is written in C++ and Python and is freely available at https://github.com/bcgsc/ntjoin.Supplementary information Supplementary data are available at Bioinformatics online.

Plasmid carrying mcr-9 from an extensively drug-resistant NDM-1-producing Klebsiella quasipneumoniae subsp. quasipneumoniae clinical isolate

The extensively drug-resistant NDM-1-producing Klebsiella quasipneumoniae subsp. quasipneumoniae M17277 (Kqsq-M17277) was only susceptible to colistin and fosfomycin. Whole genome sequencing of Kqsq-M17277 was performed through Illumina (short reads) and Oxford Nanopore Technologies (long reads). Hybrid assembly of short and long reads rendered circular contigs of the Kqsq-M17277 chromosome and three plasmids of 477,340 (p17277A_477), 138,998 (p17277C_138, which harbored blaNDM-1) and 123,307 bp (p17277B_123). In silico analysis showed that p17277A_477 was not previously described, belonged to the IncHI2 incompatibility group and contained 51 complete or partial (≥300 bp) insertion sequences and Tn3 family transposons, which encompassed, as a whole, 14% of the p17277A_477 sequence. In silico search for antimicrobial resistance genes revealed that p17277A_477 harbored the recently described plasmid-encoded mcr-9 gene (colistin resistance), besides other resistance genes. This is the first worldwide report of a K. quasipneumoniae harboring mcr-9 and blaNDM-1 and the first report of mcr-9 in Latin America. This gene was embedded in a genetic structure, i.e., IS903B-like/mcr-9/wbuC/IS26, that was found in 71% of the mcr-9-harboring sequences of the NCBI Nucleotide Collection database. The qseC and qseB genes involved in mcr-9 expression were not found in Kqsq-M17277, which could explain its colistin susceptibility. Fifteen additional resistance genes were detected in p17277A_477. Ten of them [aac(6′)-IIc, aadA2, aph(3′)-Ia, blaDHA-1, dfrA12, ereA2, mphA, qnrB4, sul1 and tet(D)] were found clustered within a 44,907-bp-long array of IS26-bounded resistance genes. The plasmid-mediated quinolone resistance gene qnrB4 was located in a complex class 1 integron within this IS26-bounded resistance gene cluster. p17277A_477 was self-transferable by conjugation to Escherichia coli J53, azide-resistant. The Tra1 and Tra2 conjugative transfer regions of the IncHI2 plasmid R478 were found by in silico search. p17277A_477 is a complex plasmid that provides a large baggage of both resistance genes and genetic resources for plasmid rearrangements and additional resistance gene acquisition.

Chromosome-level assembly of wild Bactrian camel genome reveals organization of immune gene loci.

Camelids are characterized by their unique adaptive immune system that exhibits the generation of homodimeric heavy-chain immunoglobulins, somatic hypermutation of T-cell receptors, and low genetic diversity of major histocompatibility complex (MHC) genes. However, short-read assemblies are typically highly fragmented in these gene loci owing to their repetitive and polymorphic nature. Here, we constructed a chromosome-level assembly of wild Bactrian camel genome based on high-coverage long-read sequencing and chromatin interaction mapping. The assembly with a contig N50 of 5.37Mb and a scaffold N50 of 76.03Mb, represents the most contiguous camelid genome to date. The genomic organization of immunoglobulin heavy-chain locus was similar between the wild Bactrian camel and alpaca, and genes encoding for conventional and heavy-chain antibodies were intermixed. The organizations of two immunoglobulin light-chain loci and four T cell receptor loci were also fully deciphered using the new assembly. Additionally, the complete classical MHC region was resolved into a single contig. The high-quality assembly presented here provides an essential reference for future investigations examining the camelid immune system.

Assembly-free single-molecule sequencing recovers complete virus genomes from natural microbial communities.

Viruses are the most abundant biological entities on Earth and play key roles in host ecology, evolution, and horizontal gene transfer. Despite recent progress in viral metagenomics, the inherent genetic complexity of virus populations still poses technical difficulties for recovering complete virus genomes from natural assemblages. To address these challenges, we developed an assembly-free, single-molecule nanopore sequencing approach, enabling direct recovery of complete virus genome sequences from environmental samples. Our method yielded thousands of full-length, high-quality draft virus genome sequences that were not recovered using standard short-read assembly approaches. Additionally, our analyses discriminated between populations whose genomes had identical direct terminal repeats versus those with circularly permuted repeats at their termini, thus providing new insight into native virus reproduction and genome packaging. Novel DNA sequences were discovered, whose repeat structures, gene contents, and concatemer lengths suggest they are phage-inducible chromosomal islands, which are packaged as concatemers in phage particles, with lengths that match the size ranges of co-occurring phage genomes. Our new virus sequencing strategy can provide previously unavailable information about the genome structures, population biology, and ecology of naturally occurring viruses and viral parasites.

Genome Assembly and Analysis of the North American Mountain Goat (Oreamnos americanus) Reveals Species-Level Responses to Extreme Environments.

The North American mountain goat (Oreamnos americanus) is an iconic alpine species that faces stressors from climate change, industrial development, and recreational activities. This species’ phylogenetic position within the Caprinae lineage has not been resolved and their phylogeographic history is dynamic and controversial. Genomic data could be used to address these questions and provide valuable insights to conservation and management initiatives. We sequenced short-read genomic libraries constructed from a DNA sample of a 2.5-year-old female mountain goat at 80X coverage. We improved the short-read assembly by generating Chicago library data and scaffolding using the HiRise approach. The final assembly was 2,506 Mbp in length with an N50 of 66.6 Mbp, which is within the length range and in the upper quartile for N50 published ungulate genome assemblies. Comparative analysis identified 84 gene families unique to the mountain goat. The species demographic history in terms of effective population size generally mirrored climatic trends over the past one hundred thousand years and showed a sharp decline during the last glacial maximum. This genome assembly will provide a reference basis for future population and comparative genomic analyses.

Intestinal microbiota domination under extreme selective pressures characterized by metagenomic read cloud sequencing and assembly

BackgroundLow diversity of the gut microbiome, often progressing to the point of intestinal domination by a single species, has been linked to poor outcomes in patients undergoing hematopoietic cell transplantation (HCT). Our ability to understand how certain organisms attain intestinal domination over others has been restricted in part by current metagenomic sequencing technologies that are typically unable to reconstruct complete genomes for individual organisms present within a sequenced microbial community. We recently developed a metagenomic read cloud sequencing and assembly approach that generates improved draft genomes for individual organisms compared to conventional short-read sequencing and assembly methods. Herein, we applied metagenomic read cloud sequencing to four stool samples collected longitudinally from an HCT patient preceding treatment and over the course of heavy antibiotic exposure.ResultsCharacterization of microbiome composition by taxonomic classification of reads reveals that that upon antibiotic exposure, the subject’s gut microbiome experienced a marked decrease in diversity and became dominated by Escherichia coli. While diversity is restored at the final time point, this occurs without recovery of the original species and strain-level composition. Draft genomes for individual organisms within each sample were generated using both read cloud and conventional assembly. Read clouds were found to improve the completeness and contiguity of genome assemblies compared to conventional assembly. Moreover, read clouds enabled the placement of antibiotic resistance genes present in multiple copies both within a single draft genome and across multiple organisms. The occurrence of resistance genes associates with the timing of antibiotics administered to the patient, and comparative genomic analysis of the various intestinal E. coli strains across time points as well as the bloodstream isolate showed that the subject’s E. coli bloodstream infection likely originated from the intestine. The E. coli genome from the initial pre-transplant stool sample harbors 46 known antimicrobial resistance genes, while all other species from the pre-transplant sample each contain at most 5 genes, consistent with a model of heavy antibiotic exposure resulting in selective outgrowth of the highly antibiotic-resistant E. coli.ConclusionThis study demonstrates the application and utility of metagenomic read cloud sequencing and assembly to study the underlying strain-level genomic factors influencing gut microbiome dynamics under extreme selective pressures in the clinical context of HCT.

Capture of complete ciliate chromosomes in single sequencing reads reveals widespread chromosome isoforms

BackgroundWhole-genome shotgun sequencing, which stitches together millions of short sequencing reads into a single genome, ushered in the era of modern genomics and led to a rapid expansion of the number of genome sequences available. Nevertheless, assembly of short reads remains difficult, resulting in fragmented genome sequences. Ultimately, only a sequencing technology capable of capturing complete chromosomes in a single run could resolve all ambiguities. Even “third generation” sequencing technologies produce reads far shorter than most eukaryotic chromosomes. However, the ciliate Oxytricha trifallax has a somatic genome with thousands of chromosomes averaging only 3.2 kbp, making it an ideal candidate for exploring the benefits of sequencing whole chromosomes without assembly.ResultsWe used single-molecule real-time sequencing to capture thousands of complete chromosomes in single reads and to update the published Oxytricha trifallax JRB310 genome assembly. In this version, over 50% of the completed chromosomes with two telomeres derive from single reads. The improved assembly includes over 12,000 new chromosome isoforms, and demonstrates that somatic chromosomes derive from variable rearrangements between somatic segments encoded up to 191,000 base pairs away. However, while long reads reduce the need for assembly, a hybrid approach that supplements long-read sequencing with short reads for error correction produced the most complete and accurate assembly, overall.ConclusionsThis assembly provides the first example of complete eukaryotic chromosomes captured by single sequencing reads and demonstrates that traditional approaches to genome assembly can mask considerable structural variation.