A Case Study into Microbial Genome Assembly Gap Sequences and Finishing Strategies

A Genomics-Based Discovery of Secondary Metabolite Biosynthetic Gene Clusters in the Potential Novel Strain Streptomyces sp. 21So2-11 Isolated from Antarctic Soil

Streptomyces species are attractive sources of secondary metabolites that serve as major sources of antibiotics and other drugs. In this study, genome mining was used to determine the biosynthetic potential of Streptomyces sp. 21So2-11 isolated from Antarctic soil. 16S rRNA gene sequencing revealed that this strain is most closely related to Streptomyces drozdowiczii NBRC 101007T, with a similarity of 98.02%. Genome comparisons based on average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) showed that strain 21So2-11 represents a novel species of the genus Streptomyces. In addition to a large number of genes related to environmental adaptation and ecological function, a total of 28 putative biosynthetic gene clusters (BGCs) responsible for the biosynthesis of known and/or novel secondary metabolites, including terpenes, lantipeptides, polyketides, nonribosomal peptides, RiPPs and siderophores, were detected in the genome of strain 21So2-11. In addition, a total of 1456 BGCs were predicted to contribute to the biosynthesis of more than 300 secondary metabolites based on the genomes of 47 Streptomyces strains originating from polar regions. The results indicate the potential of Streptomyces sp. 21So2-11 for bioactive secondary metabolite production and are helpful for understanding bacterial adaptability and ecological function in cold terrestrial environments.

Genotype diversity of brucellosis agents isolated from humans and animals in Greece based on whole-genome sequencing

BACKGROUND

Brucellosis is a zoonotic disease whose causative agent, Brucella spp., is endemic in many countries of the Mediterranean basin, including Greece. Although the occurrence of brucellosis must be reported to the authorities, it is believed that the disease is under-reported in Greece, and knowledge about the genomic diversity of brucellae is lacking.

METHODS

Thus, 44 Brucella isolates, primarily B. melitensis, collected between 1999 and 2009 from humans and small ruminants in Greece were subjected to whole genome sequencing using short-read technology. The raw reads and assembled genomes were used for in silico genotyping based on single nucleotide substitutions and alleles. Further, specific genomic regions encoding putative virulence genes were screened for characteristic nucleotide changes, which arose in different genotype lineages.

RESULTS

In silico genotyping revealed that the isolates belonged to three of the known sublineages of the East Mediterranean genotype. In addition, a novel subgenotype was identified that was basal to the other East Mediterranean sublineages, comprising two Greek strains. The majority of the isolates can be assumed to be of endemic origin, as they were clustered with strains from the Western Balkans or Turkey, whereas one strain of human origin could be associated with travel to another endemic region, e.g. Portugal. Further, nucleotide substitutions in the housekeeping gene rpoB and virulence-associated genes were detected, which were characteristic of the different subgenotypes. One of the isolates originating from an aborted bovine foetus was identified as B. abortus vaccine strain RB51.

CONCLUSION

The results demonstrate the existence of several distinct persistent Brucella sp. foci in Greece. To detect these and for tracing infection chains, extensive sampling initiatives are required.

The role of plasmids in carbapenem resistant E. coli in Alameda County, California

BACKGROUND

Antimicrobial resistant infections continue to be a leading global public health crisis. Mobile genetic elements, such as plasmids, have been shown to play a major role in the dissemination of antimicrobial resistance (AMR) genes. Despite its ongoing threat to human health, surveillance of AMR in the United States is often limited to phenotypic resistance. Genomic analyses are important to better understand the underlying resistance mechanisms, assess risk, and implement appropriate prevention strategies. This study aimed to investigate the extent of plasmid mediated antimicrobial resistance that can be inferred from short read sequences of carbapenem resistant E. coli (CR-Ec) in Alameda County, California. E. coli isolates from healthcare locations in Alameda County were sequenced using an Illumina MiSeq and assembled with Unicycler. Genomes were categorized according to predefined multilocus sequence typing (MLST) and core genome multilocus sequence typing (cgMLST) schemes. Resistance genes were identified and corresponding contigs were predicted to be plasmid-borne or chromosome-borne using two bioinformatic tools (MOB-suite and mlplasmids).

RESULTS

Among 82 of CR-Ec identified between 2017 and 2019, twenty-five sequence types (STs) were detected. ST131 was the most prominent (n = 17) followed closely by ST405 (n = 12). bla_CTX-M were the most common ESBL genes and just over half (18/30) of these genes were predicted to be plasmid-borne by both MOB-suite and mlplasmids. Three genetically related groups of E. coli isolates were identified with cgMLST. One of the groups contained an isolate with a chromosome-borne bla_CTX-M-15 gene and an isolate with a plasmid-borne bla_CTX-M-15 gene.

CONCLUSIONS

This study provides insights into the dominant clonal groups driving carbapenem resistant E. coli infections in Alameda County, CA, USA clinical sites and highlights the relevance of whole-genome sequencing in routine local genomic surveillance. The finding of multi-drug resistant plasmids harboring high-risk resistance genes is of concern as it indicates a risk of dissemination to previously susceptible clonal groups, potentially complicating clinical and public health intervention.

Chromosome-level reference genome of Tetrastigma hemsleyanum (Vitaceae) provides insights into genomic evolution and the biosynthesis of phenylpropanoids and flavonoids

Here, we present a high-quality chromosome-scale genome assembly (2.19 Gb) and annotation of Tetrastigma hemsleyanum, a perennial herbaceous liana native to subtropical China with diverse medicinal applications. Approximately 73% of the genome was comprised of transposable elements (TEs), of which long terminal repeat retrotransposons (LTR-RTs) were a predominant group (69% of the genome). The genome size increase of T. hemsleyanum (relative to Vitis species) was mostly due to the proliferation of LTR-RTs. Of the different modes of gene duplication identified, transposed duplication (TRD) and dispersed duplication (DSD) were the predominant ones. Genes, particularly those involved in the phenylpropanoid-flavonoid (PF) pathway and those associated with therapeutic properties and environmental stress resistance, were significantly amplified through recent tandem duplications. We dated the divergence of two intraspecific lineages in Southwest (SW) versus Central-South-East (CSE) China to the late Miocene (approximately 5.2 million years ago). Of those, the former showed more upregulated genes and metabolites. Based on resequencing data of 38 individuals representing both lineages, we identified various candidate genes related to 'response to stimulus' and 'biosynthetic process', including ThFLS11, which is putatively involved in flavonoid accumulation. Overall, this study provides abundant genomic resources for future evolutionary, ecological, and functional genomics studies in T. hemsleyanum and related species.

Advances in experimental and computational methodologies for the study of microbial-surface interactions at different omics levels

The study of the biological response of microbial cells interacting with natural and synthetic interfaces has acquired a new dimension with the development and constant progress of advanced omics technologies. New methods allow the isolation and analysis of nucleic acids, proteins and metabolites from complex samples, of interest in diverse research areas, such as materials sciences, biomedical sciences, forensic sciences, biotechnology and archeology, among others. The study of the bacterial recognition and response to surface contact or the diagnosis and evolution of ancient pathogens contained in archeological tissues require, in many cases, the availability of specialized methods and tools. The current review describes advances in in vitro and in silico approaches to tackle existing challenges (e.g., low-quality sample, low amount, presence of inhibitors, chelators, etc.) in the isolation of high-quality samples and in the analysis of microbial cells at genomic, transcriptomic, proteomic and metabolomic levels, when present in complex interfaces. From the experimental point of view, tailored manual and automatized methodologies, commercial and in-house developed protocols, are described. The computational level focuses on the discussion of novel tools and approaches designed to solve associated issues, such as sample contamination, low quality reads, low coverage, etc. Finally, approaches to obtain a systems level understanding of these complex interactions by integrating multi omics datasets are presented.

Chromosome-level assembly and analysis of the Thymus genome provide insights into glandular secretory trichome formation and monoterpenoid biosynthesis in thyme

Thyme has medicinal and aromatic value because of its potent antimicrobial and antioxidant properties. However, the absence of a fully sequenced thyme genome limits functional genomic studies of Chinese native thymes. Thymus quinquecostatus Čelak., which contains large amounts of bioactive monoterpenes such as thymol and carvacrol, is an important wild medicinal and aromatic plant in China. Monoterpenoids are abundant in glandular secretory trichomes. Here, high-fidelity and chromatin conformation capture technologies were used to assemble and annotate the T. quinquecostatus genome at the chromosome level. The 13 chromosomes of T. quinquecostatus had a total length of 528.66 Mb, a contig N50 of 8.06 Mb, and a BUSCO score of 97.34%. We found that T. quinquecostatus had experienced two whole-genome duplications, with the most recent event occurring ∼4.34 million years ago. Deep analyses of the genome, in conjunction with comparative genomic, phylogenetic, transcriptomic, and metabonomic studies, uncovered many regulatory factors and genes related to monoterpenoids and glandular secretory trichome development. Genes encoding terpene synthase (TPS), cytochrome P450 monooxygenases (CYPs), short-chain dehydrogenase/reductase (SDR), R2R3-MYB, and homeodomain-leucine zipper (HD-ZIP) IV were among those present in the T. quinquecostatus genome. Notably, Tq02G002290.1 (TqTPS1) was shown to encode the terpene synthase responsible for catalyzing production of the main monoterpene product γ-terpinene from geranyl diphosphate (GPP). Our study provides significant insight into the mechanisms of glandular secretory trichome formation and monoterpenoid biosynthesis in thyme. This work will facilitate the development of molecular breeding tools to enhance the production of bioactive secondary metabolites in Lamiaceae.

Clostridium autoethanogenum isopropanol production via native plasmid pCA replicon

Clostridium autoethanogenum is a model gas-fermenting acetogen for commercial ethanol production. It is also a platform organism being developed for the carbon-negative production of acetone and isopropanol by gas fermentation. We have assembled a 5.5 kb pCA plasmid for type strain DSM10061 (JA1-1) using three genome sequence datasets. pCA is predicted to encode seven open-reading frames and estimated to be a low-copy number plasmid present at approximately 12 copies per chromosome. RNA-seq analyses indicate that pCA genes are transcribed at low levels and two proteins, CAETHG_05090 (putative replication protein) and CAETHG_05115 (hypothetical, a possible Mob protein), were detected at low levels during batch gas fermentations. Thiolase (thlA), CoA-transferase (ctfAB), and acetoacetate decarboxylase (adc) genes were introduced into a vector for isopropanol production in C. autoethanogenum using the native plasmid origin of replication. The availability of the pCA sequence will facilitate studies into its physiological role and could form the basis for genetic tool optimization.

Genomic Surveillance of Clinical Pseudomonas aeruginosa Isolates Reveals an Additive Effect of Carbapenemase Production on Carbapenem Resistance

Carbapenem resistance in Pseudomonas aeruginosa is increasing globally, and surveillance to define the mechanisms of such resistance in low- and middle-income countries is limited. This study establishes the genotypic mechanisms of β-lactam resistance by whole-genome sequencing (WGS) in 142 P. aeruginosa clinical isolates recovered from three hospitals in Islamabad and Rawalpindi, Pakistan between 2016 and 2017. Isolates were subjected to antimicrobial susceptibility testing (AST) by Kirby-Bauer disk diffusion, and their genomes were assembled from Illumina sequencing data. β-lactam resistance was high, with 46% of isolates resistant to piperacillin-tazobactam, 42% to cefepime, 48% to ceftolozane-tazobactam, and 65% to at least one carbapenem. Twenty-two percent of isolates were resistant to all β-lactams tested. WGS revealed that carbapenem resistance was associated with the acquisition of metallo-β-lactamases (MBLs) or extended-spectrum β-lactamases (ESBLs) in the blaGES, blaVIM, and blaNDM families, and mutations in the porin gene oprD. These resistance determinants were found in globally distributed lineages, including ST235 and ST664, as well as multiple novel STs which have been described in a separate investigation. Analysis of AST results revealed that acquisition of MBLs/ESBLs on top of porin mutations had an additive effect on imipenem resistance, suggesting that there is a selective benefit for clinical isolates to encode multiple resistance determinants to the same drugs. The strong association of these resistance determinants with phylogenetic background displays the utility of WGS for monitoring carbapenem resistance in P. aeruginosa, while the presence of these determinants throughout the phylogenetic tree shows that knowledge of the local epidemiology is crucial for guiding potential treatment of multidrug-resistant P. aeruginosa infections. IMPORTANCE Pseudomonas aeruginosa is associated with serious infections, and treatment can be challenging. Because of this, carbapenems and β-lactam/β-lactamase inhibitor combinations have become critical tools in treating multidrug-resistant (MDR) P. aeruginosa infections, but increasing resistance threatens their efficacy. Here, we used WGS to study the genotypic and phylogenomic patterns of 142 P. aeruginosa isolates from the Potohar region of Pakistan. We sequenced both MDR and antimicrobial susceptible isolates and found that while genotypic and phenotypic patterns of antibiotic resistance correlated with phylogenomic background, populations of MDR P. aeruginosa were found in all major phylogroups. We also found that isolates possessing multiple resistance mechanisms had significantly higher levels of imipenem resistance compared to the isolates with a single resistance mechanism. This study demonstrates the utility of WGS for monitoring patterns of antibiotic resistance in P. aeruginosa and potentially guiding treatment choices based on the local spread of β-lactamase genes.

Whole-Genome Sequence Datasets: A Powerful Resource for the Food Microbiology Laboratory Toolbox

Whole-genome sequencing (WGS) technologies are rapidly being adopted for routine use in food microbiology laboratories worldwide. Examples of how WGS is used to support food safety testing include gene marker discovery (e.g., virulence and anti-microbial resistance gene determination) and high-resolution typing (e.g., cg/wgMLST analysis). This has led to the establishment of large WGS databases representing the genomes of thousands of different types of food pathogenic and commensal bacteria. This information constitutes an invaluable resource that can be leveraged to develop and validate routine test methods used to support regulatory and industry food safety objectives. For example, well-curated raw and assembled genomic datasets of the key food pathogens (Salmonella enterica, Listeria monocytogenes, and Shiga-toxigenic Escherichia coli) have been used in our laboratory in studies to validate bioinformatics pipelines, as well as new molecular methods as a prelude to the laboratory phase of the “wet lab” validation process. The application of genomic information to food microbiology method development will decrease the cost of test development and lead to the generation of more robust methodologies supporting risk assessment and risk management actions.

Chromosome-level assembly of the mangrove plant Aegiceras corniculatum genome generated through Illumina, PacBio and Hi-C sequencing technologies

Aegiceras corniculatum is a major mangrove plant species adapted to waterlogging and saline conditions, grows in the coastal intertidal zone of tropical and subtropical regions. Here, we present a chromosome-level genome assembly of A. corniculatum by incorporating PacBio long-read sequencing and Hi-C technology. The results showed that the PacBio draft genome size is 906.63 Mb. Hi-C scaffolding anchored 885.06 Mb contigs (97.62% of draft assembly) onto 24 pseudochromosomes. The contig N50 and scaffold N50 were 7.1 Mb and 37.74 Mb, respectively. Out of 40,727 protein-coding genes predicted in the study, 89% have functional annotations in public databases. We also showed that of the 603.93 Mb repetitive sequences predicted in the assembled genome, long terminal repeat retrotransposons constitute 41.52%. The genome evolution analysis showed that the A. corniculatum genome experienced two whole-genome duplication events and shared the ancient γ whole-genome triplication event. A comparative genomic analysis revealed an incidence of expansion in 1,488 gene families associated with essential metabolism and biosynthetic pathways, including photosynthesis, oxidative phosphorylation, phenylalanine, glyoxylate, dicarboxylate metabolism, and DNA replication, which probably constitute adaptation traits that allow the A. corniculatum to survive in the intertidal zone. Also, the systematic characterization of genes associated with flavonoid biosynthesis pathway and the AcNHX gene family conducted in this study will provide insight into the adaptation mechanism of A. corniculatum to intertidal environments. The high-quality genome reported here can provide historical insights into genomic transformations that support the survival of A. corniculatum under harsh intertidal habitats.

What Is the Best Lens? Comparing the Resolution Power of Genome-Derived Markers and Standard Barcodes

Fungal species delimitation was traditionally carried out with multicopy ribosomal RNA (rRNA) genes, principally for their ease of amplification. Since the efficacy of these markers has been questioned, single-copy protein-encoding genes have been proposed alone or in combination for Multi-Locus Sequence Typing (MLST). In this context, the role of the many sequences obtained with Next-Generation Sequencing (NGS) techniques, in both genomics and metagenomics, further pushes toward an analysis of the efficacy of NGS-derived markers and of the metrics to evaluate the marker efficacy in discriminating fungal species. This paper aims at proposing MeTRe (Mean Taxonomic Resolution), a novel index that could be used both for measuring marker efficacy and for assessing the actual resolution (i.e., the level of separation) between species obtained with different markers or their combinations. In this paper, we described and then employed this index to compare the efficacy of two rRNAs and four single-copy markers obtained from public databases as both an amplicon-based approach and genome-derived sequences. Two different groups of species were used, one with a pathogenic species of Candida that was characterized by relatively well-separated taxa, whereas the other, comprising some relevant species of the sensu stricto group of the genus Saccharomyces, included close species and interspecific hybrids. The results showed the ability of MeTRe to evaluate marker efficacy in general and genome-derived markers specifically.

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.

The genome evolution and low-phosphorus adaptation in white lupin

White lupin (Lupinus albus) is a legume crop that develops cluster roots and has high phosphorus (P)-use efficiency (PUE) in low-P soils. Here, we assemble the genome of white lupin and find that it has evolved from a whole-genome triplication (WGT) event. We then decipher its diploid ancestral genome and reconstruct the three sub-genomes. Based on the results, we further reveal the sub-genome dominance and the genic expression of the different sub-genomes varying in relation to their transposable element (TE) density. The PUE genes in white lupin have been expanded through WGT as well as tandem and dispersed duplications. Furthermore, we characterize four main pathways for high PUE, which include carbon fixation, cluster root formation, soil-P remobilization, and cellular-P reuse. Among these, auxin modulation may be important for cluster root formation through involvement of potential genes LaABCG36s and LaABCG37s. These findings provide insights into the genome evolution and low-P adaptation of white lupin.

Plasmid Reconstruction from Next-Gen Data: A Detailed Protocol for the Use of PLACNETw for the Reconstruction of Plasmids from WGS Datasets

Mobile Genetic Elements (MGE) play essential roles in adaptive bacterial evolution, facilitating genetic exchange for extrachromosomal DNA, especially antibiotic resistance genes and virulence factors. For this reason, high-throughput next-generation sequencing of bacteria is of great relevance, especially for clinical pathogenic bacteria. Accurate identification of MGE from whole-genome sequencing (WGS) datasets is one of the major challenges, still hindered by methodological limitations and high sequencing costs.This chapter encompasses the protocol used for plasmid reconstruction by applying the PLACNETw methodology, from raw reads to assembled plasmids and chromosome. PLACNETw is a graphical user-friendly interface to visualize and reconstruct MGE from short-read WGS datasets. No bioinformatic background or sophisticated computational resources are required and high precision and sensitivity are achieved.

First genomic insights into carbapenem-resistant Klebsiella pneumoniae from Malaysia

OBJECTIVES

Despite the increasing reports of carbapenem-resistant Enterobacteriaceae in Malaysia, genomic resources for carbapenem-resistant clinical strains of Klebsiella pneumoniae (K. pneumoniae) remain unavailable. This study aimed to sequence the genomes of multiple carbapenem-resistant K. pneumoniae strains from Malaysia and to identify the genetic basis for their resistance.

METHODS

Illumina whole genome sequencing was performed on eight carbapenem-resistant K. pneumoniae isolated from a Malaysian hospital. Genetic diversity was inferred from the assembled genomes based on in silico multilocus sequence typing (MLST). In addition, plasmid-derived and chromosome-derived contigs were predicted using the machine learning approach. After genome annotation, genes associated with carbapenem resistance were identified based on similarity searched against the ResFinder database.

RESULTS

The eight K. pneumoniae isolates were grouped into six different sequence types, some of which were represented by a single isolate in the MLST database. Genomic potential for carbapenem-resistance was attributed to the presence of plasmid-localised bla_NDM (bla_NDM-1/bla_NDM-5) or bla_KPC (bla_KPC-2/bla_KPC-6) in these sequenced strains. The majority of these carbapenem resistance genes was flanked by repetitive (transposase or integrase) sequences, suggesting their potential mobility. This study also reported the first bla_KPC-6-harbouring plasmid contig to be assembled for K. pneumoniae, and the second for the genus Klebsiella.

CONCLUSION

This study reported the first genomic resources for carbapenem-resistant K. pneumoniae from Malaysia. The high diversity of carbapenem resistance genes and sequence types uncovered from eight isolates from the same hospital is worrying and indicates an urgent need to improve the genomic surveillance of clinical K. pneumoniae in Malaysia.

Pushing the limits of de novo genome assembly for complex prokaryotic genomes harboring very long, near identical repeats

Generating a complete, de novo genome assembly for prokaryotes is often considered a solved problem. However, we here show that Pseudomonas koreensis P19E3 harbors multiple, near identical repeat pairs up to 70 kilobase pairs in length, which contained several genes that may confer fitness advantages to the strain. Its complex genome, which also included a variable shufflon region, could not be de novo assembled with long reads produced by Pacific Biosciences’ technology, but required very long reads from Oxford Nanopore Technologies. Importantly, a repeat analysis, whose results we release for over 9600 prokaryotes, indicated that very complex bacterial genomes represent a general phenomenon beyond Pseudomonas. Roughly 10% of 9331 complete bacterial and a handful of 293 complete archaeal genomes represented this ‘dark matter’ for de novo genome assembly of prokaryotes. Several of these ‘dark matter’ genome assemblies contained repeats far beyond the resolution of the sequencing technology employed and likely contain errors, other genomes were closed employing labor-intense steps like cosmid libraries, primer walking or optical mapping. Using very long sequencing reads in combination with assembly algorithms capable of resolving long, near identical repeats will bring most prokaryotic genomes within reach of fast and complete de novo genome assembly.

Comparative genomic analysis of monosporidial and monoteliosporic cultures for unraveling the complexity of molecular pathogenesis of Tilletia indica pathogen of wheat

Tilletia indica (Ti) - a quarantined fungal pathogen of wheat and its pathogenesis is chiefly governed by pathogen effectors secreted inside the host plant. The de novo genome sequencing of several field isolates and stages available could be used for understanding the molecular pathogenesis. The presence of gaps and low coverage of assembled genomes poses a problem in accurate functional annotation of such functions. In the present study attempts were made to improve the Ti draft genome through reconciliation of globally available datasets of three highly virulent monoteliospore cultures of Ti field isolates. It has sequence depth of 107x and N50 scaffold size of 80,772 (more than 26 times as large as achieved in the draft assembly) with highest sequence contiguity, more accurate and nearly complete. Functional annotation revealed that Ti genome contains 9209 genes evolved with many expanded gene families and arranged mostly in a cluster. About 79% of Ti genes were orthologous to other basidiomycetes fungi, Around 7.93% proteins were having secretary signals and 6.66% were identified as highly virulent pathogenicity genes. Using improved Ti genome as a reference, the genomic variation was assessed with respect to repeats, SNPs/InDel, gene families and correct set of virulence associated genes during its life cycle. The comparative intra-species, inter-stage and inter-species genomic variation will have broader implications to understand the gene regulatory networks involved in growth, mating and virulence behaviour of Tilletia f. spp. and also for better appreciation of fungal biology and disease management.

Graph analysis of fragmented long-read bacterial genome assemblies

Motivation

Long-read genome assembly tools are expected to reconstruct bacterial genomes nearly perfectly; however, they still produce fragmented assemblies in some cases. It would be beneficial to understand whether these cases are intrinsically impossible to resolve, or if assemblers are at fault, implying that genomes could be refined or even finished with little to no additional experimental cost.

Results

We propose a set of computational techniques to assist inspection of fragmented bacterial genome assemblies, through careful analysis of assembly graphs. By finding paths of overlapping raw reads between pairs of contigs, we recover potential short-range connections between contigs that were lost during the assembly process. We show that our procedure recovers 45% of missing contig adjacencies in fragmented Canu assemblies, on samples from the NCTC bacterial sequencing project. We also observe that a simple procedure based on enumerating weighted Hamiltonian cycles can suggest likely contig orderings. In our tests, the correct contig order is ranked first in half of the cases and within the top-three predictions in nearly all evaluated cases, providing a direction for finishing fragmented long-read assemblies.

Availability and implementation

https://gitlab.inria.fr/pmarijon/knot .

Supplementary information

Supplementary data are available at Bioinformatics online.

Nanopore ultra-long read sequencing technology for antimicrobial resistance detection in Mannheimia haemolytica

Disruptive innovations in long-range, cost-effective direct template nucleic acid sequencing are transforming clinical and diagnostic medicine. A multidrug resistant strain and a pan-susceptible strain of Mannheimia haemolytica, isolated from pneumonic bovine lung samples, were sequenced at 146× and 111× coverage, respectively with Oxford Nanopore Technologies MinION. De novo assembly produced a complete genome for the non-resistant strain and a nearly complete assembly for the drug resistant strain. Functional annotation using RAST (Rapid Annotations using Subsystems Technology), CARD (Comprehensive Antibiotic Resistance Database) and ResFinder databases identified genes conferring resistance to different classes of antibiotics including β-lactams, tetracyclines, lincosamides, phenicols, aminoglycosides, sulfonamides and macrolides. Resistance phenotypes of the M. haemolytica strains were determined by minimum inhibitory concentration (MIC) of the antibiotics. Sequencing with a highly portable MinION device corresponded to MIC assays with most of the antimicrobial resistant determinants being identified with as few as 5437 reads, except for the genes responsible for resistance to Fluoroquinolones. The resulting quality assemblies and AMR gene annotation highlight the efficiency of ultra-long read, whole-genome sequencing (WGS) as a valuable tool in diagnostic veterinary medicine.

Improved genome of Agrobacterium radiobacter type strain provides new taxonomic insight into Agrobacterium genomospecies 4

The reported Agrobacterium radiobacter DSM 30174^T genome is highly fragmented, hindering robust comparative genomics and genome-based taxonomic analysis. We re-sequenced the Agrobacterium radiobacter type strain, generating a dramatically improved genome with high contiguity. In addition, we sequenced the genome of Agrobacterium tumefaciens B6^T, enabling for the first time, a proper comparative genomics of these contentious Agrobacterium species. We provide concrete evidence that the previously reported Agrobacterium radiobacter type strain genome (Accession Number: ASXY01) is contaminated which explains its abnormally large genome size and fragmented assembly. We propose that Agrobacterium tumefaciens be reclassified as Agrobacterium radiobacter subsp. tumefaciens and that Agrobacterium radiobacter retains it species status with the proposed name of Agrobacterium radiobacter subsp. radiobacter. This proposal is based, first on the high pairwise genome-scale average nucleotide identity supporting the amalgamation of both Agrobacterium radiobacter and Agrobacterium tumefaciens into a single species. Second, maximum likelihood tree construction based on the concatenated alignment of shared genes (core genes) among related strains indicates that Agrobacterium radiobacter NCPPB3001 is sufficiently divergent from Agrobacterium tumefaciens to propose two independent sub-clades. Third, Agrobacterium tumefaciens demonstrates the genomic potential to synthesize the L configuration of fucose in its lipid polysaccharide, fostering its ability to colonize plant cells more effectively than Agrobacterium radiobacter.

Polyglutamine Repeats in Viruses

This review explores the presence and functions of polyglutamine (polyQ) in viral proteins. In mammals, mutations in polyQ segments (and CAG repeats at the nucleotide level) have been linked to neural disorders and ataxias. PolyQ regions in normal human proteins have documented functional roles, in transcription factors and, more recently, in regulating autophagy. Despite the high frequency of polyQ repeats in eukaryotic genomes, little attention has been given to the presence or possible role of polyQ sequences in virus genomes. A survey described here revealed that polyQ repeats occur rarely in RNA viruses, suggesting that they have detrimental effects on virus replication at the nucleotide or protein level. However, there have been sporadic reports of polyQ segments in potyviruses and in reptilian nidoviruses (among the largest RNA viruses known). Conserved polyQ segments are found in the regulatory control proteins of many DNA viruses. Variable length polyQ tracts are found in proteins that contribute to transmissibility (cowpox A-type inclusion protein (ATI)) and control of latency (herpes viruses). New longer-read sequencing methods, using original biological samples, should reveal more details on the presence and functional role of polyQ in viruses, as well as the nucleotide regions that encode them. Given the known toxic effects of polyQ repeats, the role of these segments in neurovirulent and tumorigenic viruses should be further explored.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.