Memory-Efficient Assembly Using Flye

Genomic analysis of antibiotic resistance genes and mobile genetic elements in eight strains of nontyphoid Salmonella

Nontyphoidal Salmonella (NTS) is the main etiological agent of human nontyphoidal salmonellosis. The aim of this study was to analyze the epidemiological characteristics and horizontal transfer mechanisms of antimicrobial resistance (AMR) genes from eight strains of NTS detected in Zhenjiang City, Jiangsu Province, China. Fecal samples from outpatients with food-borne diarrhea were collected in 2022. The NTS isolates were identified, and their susceptibility was tested with the Vitek 2 Compact system. The genomes of the NTS isolates were sequenced with the Illumina NovaSeq platform and Oxford Nanopore Technologies platform. The AMR genes and mobile genetic elements (MGEs) were predicted with the relevant open access resources. Eight strains of NTS were isolated from 153 specimens, and Salmonella Typhimurium ST19 was the most prevalent serotype. The AMR gene with the highest detection rate was AAC(6')-Iaa (10.5%) followed by TEM-1 (7.9%), sul2 (6.6%), and tet(A) (5.3%). Eleven MGEs carrying 34 AMR genes were identified on the chromosomes of 3 of the 8 NTS, including 3 resistance islands, 6 composite transposons (Tns), and 2 integrons. Eighteen plasmids carrying 40 AMR genes were detected in the 8 NTS strains, including 6 mobilizable plasmids, 3 conjugative plasmids, and 9 nontransferable plasmids, 7 of which carried 10 composite Tns and 3 integrons. This study provided a theoretical basis, from a genetic perspective, for the prevention and control of NTS resistance in Zhenjiang City.

IMPORTANCE

Human nontyphoidal salmonellosis is one of the common causes of bacterial food-borne illnesses, with significant social and economic impacts, especially those caused by invasive multidrug-resistant nontyphoidal Salmonella, which entails high morbidity and mortality. Antimicrobial resistance is mainly mediated by drug resistance genes, and mobile genetic elements play key roles in the capture, accumulation, and dissemination of antimicrobial resistance genes. Therefore, it is necessary to study the epidemiological characteristics and horizontal transfer mechanisms of antimicrobial resistance genes of nontyphoidal Salmonella to prevent the spread of multidrug-resistant nontyphoidal Salmonella.

Characterization and Genome Analysis of the Delftia lacustris Strain LzhVag01 Isolated from Vaginal Discharge

Delftia has been separated from freshwater, sludge, and soil and has emerged as a novel opportunistic pathogen in the female vagina. However, the genomic characteristics, pathogenicity, and biotechnological properties still need to be comprehensively investigated. In this study, a Delftia strain was isolated from the vaginal discharge of a 43-year-old female with histologically confirmed cervical intraepithelial neoplasm (CIN III), followed by whole-genome sequencing. Phylogenetic analysis and average nucleotide identity (ANI) analysis demonstrated that it belongs to Delftia lacustris, named D. lacustris strain LzhVag01. LzhVag01 was sensitive to β-lactams, macrolides, and tetracyclines but exhibited resistance to lincoamines, nitroimidazoles, aminoglycosides, and fluoroquinolones. Its genome is a single, circular chromosome of 6,740,460 bp with an average GC content of 66.59%. Whole-genome analysis identified 16 antibiotic resistance-related genes, which match the antimicrobial susceptibility profile of this strain, and 11 potential virulence genes. These pathogenic factors may contribute to its colonization in the vaginal environment and its adaptation and accelerate the progression of cervical cancer. This study sequenced and characterized the whole-genome of Delftia lacustris isolated from vaginal discharge, which provides investigators and clinicians with valuable insights into this uncommon species.

Complete genome sequence of Enterococcus faecalis phage EF_RCK

A lytic bacteriophage EF_RCK infecting Enterococcus faecalis was isolated from a water sample collected in a raw cockle storage container at Taman Ria market, Sungai Petani, Malaysia. The phage has a 57,848-bp double-stranded DNA genome harboring 107 protein-encoding genes and shares 90.9% nucleotide similarity with Enterococcus phage EFKL (Saphexavirus genus).

Klebsiella pneumoniae sequence type 147: a high-risk clone increasingly associated with plasmids carrying both resistance and virulence elements

Introduction. The first hybrid resistance/virulence plasmid, combining elements from virulence plasmids described in hypervirulent types of Klebsiella pneumoniae with those from conjugative resistance plasmids, was described in an isolate of sequence type (ST) 147 from 2016. Subsequently, this type has been increasingly associated with these plasmids.Hypothesis or gap statement. The extent of carriage of hybrid virulence/resistance plasmids in nosocomial isolates of K. pneumoniae requires further investigation.Aim. To describe the occurrence of virulence/resistance plasmids among isolates of K. pneumoniae received by the UK reference laboratory, particularly among representatives of ST147, and to compare their sequences.Methodology. Isolates received by the laboratory during 2022 and the first half of 2023 (n=1278) were screened for virulence plasmids by PCR detection of rmpA/rmpA2 and typed by variable-number tandem repeat analysis. Twenty-nine representatives of ST147 (including a single-locus variant) from seven hospital laboratories were subjected to long-read nanopore sequencing using high-accuracy q20 chemistry to provide complete assemblies.Results. rmpA/rmpA2 were detected in 110 isolates, of which 59 belonged to hypervirulent K1-ST23, K2-ST86 and K2-ST65/375. Of the remainder, representatives of ST147 formed the largest group, with 22 rmpA/rmpA2-positive representatives (out of 47 isolates). Representatives were from 19 hospital laboratories, with rmpA/rmpA2-positive isolates from 10. Nanopore sequencing of 29 representatives of ST147 divided them into those with no virulence plasmid (n=12), those with non-New Delhi metallo-β-lactamase (NDM) virulence plasmids (n=6) and those carrying bla NDM-5 (n=9) or bla NDM-1 (n=2) virulence plasmids. These plasmids were of IncFIB(pNDM-Mar)/IncHI1B(pNDM-MAR) replicon types. Most of the non-NDM virulence plasmids were highly similar to the originally described KpvST147L_NDM plasmid. Those carrying bla NDM-5 were highly similar to one another and to previously described plasmids in ST383 and carried an extensive array of resistance genes. Comparison of the fully assembled chromosomes indicated multiple introductions of ST147 in UK hospitals.Conclusion. This study highlights the high proportion of representatives of ST147 that carry IncFIB(pNDM-Mar)/IncHI1B(pNDM-MAR) hybrid resistance virulence plasmids. It is important to be aware of the high probability that representatives of this type carry these plasmids combining resistance and virulence determinants and of the consequent increased risk to patients.

Complete genome sequences of a Canadian strain of enteroaggregative Escherichia coli (EAEC) with multiple metals and antimicrobial resistance genes isolated from municipal waste-activated sludge

Enteroaggregative Escherichia coli (EAEC) is an emerging food-borne pathogen causing acute or persistent diarrhea in humans. Here, we report the complete genome sequence of a strain of EAEC with multiple metals and antimicrobial resistance genes isolated from a waste-activated sludge collected from a Canadian municipal wastewater treatment plant.

Comparative genome analysis of the freshwater fungus Filosporella fistucella indicates potential for plant-litter degradation at cold temperatures

Freshwater fungi play an important role in the decomposition of organic matter of leaf litter in rivers and streams. They also possess the necessary mechanisms to endure lower temperatures caused by habitat and weather variations. This includes the production of cold-active enzymes and antifreeze proteins. To better understand the physiological activities of freshwater fungi in their natural environment, different methods are being applied, and genome sequencing is one in the spotlight. In our study, we sequenced the first genome of the freshwater fungus Filosporella fistucella (45.7 Mb) and compared the genome with the evolutionary close-related species Tricladium varicosporioides (48.2 Mb). The genomes were annotated using the carbohydrate-active enzyme database where we then filtered for leaf-litter degradation-related enzymes (cellulase, hemicellulase, laccase, pectinase, cutinase, amylase, xylanase, and xyloglucanase). Those enzymes were analyzed for antifreeze properties using a machine-learning approach. We discovered that F. fistucella has more enzymes to participate in the breakdown of sugar, leaf, and wood than T. varicosporioides (855 and 719, respectively). Filosporella fistucella shows a larger set of enzymes capable of resisting cold temperatures than T. varicosporioides (75 and 66, respectively). Our findings indicate that in comparison with T. varicosporioides, F. fistucella has a greater capacity for aquatic growth, adaptability to freshwater environments, and resistance to low temperatures.

Genotypic Evolution of Klebsiella pneumoniae Sequence Type 512 during Ceftazidime/Avibactam, Meropenem/Vaborbactam, and Cefiderocol Treatment, Italy

In February 2022, a critically ill patient colonized with a carbapenem-resistant K. pneumoniae producing KPC-3 and VIM-1 carbapenemases was hospitalized for SARS-CoV-2 in the intensive care unit of Policlinico Umberto I hospital in Rome, Italy. During 95 days of hospitalization, ceftazidime/avibactam, meropenem/vaborbactam, and cefiderocol were administered consecutively to treat 3 respiratory tract infections sustained by different bacterial agents. Those therapies altered the resistome of K. pneumoniae sequence type 512 colonizing or infecting the patient during the hospitalization period. In vivo evolution of the K. pneumoniae sequence type 512 resistome occurred through plasmid loss, outer membrane porin alteration, and a nonsense mutation in the cirA siderophore gene, resulting in high levels of cefiderocol resistance. Cross-selection can occur between K. pneumoniae and treatments prescribed for other infective agents. K. pneumoniae can stably colonize a patient, and antimicrobial-selective pressure can promote progressive K. pneumoniae resistome evolution, indicating a substantial public health threat.

Genetic background of neomycin resistance in clinical Escherichia coli isolated from Danish pig farms

Neomycin is the first-choice antibiotic for the treatment of porcine enteritis caused by enterotoxigenic Escherichia coli. Resistance to this aminoglycoside is on the rise after the increased use of neomycin due to the ban on zinc oxide. We identified the neomycin resistance determinants and plasmid contents in a historical collection of 128 neomycin-resistant clinical E. coli isolates from Danish pig farms. All isolates were characterized by whole-genome sequencing and antimicrobial susceptibility testing, followed by conjugation experiments and long-read sequencing of eight selected representative strains. We detected 35 sequence types (STs) with ST100 being the most prevalent lineage (38.3%). Neomycin resistance was associated with two resistance genes, namely aph(3')-Ia and aph(3')-Ib, which were identified in 93% and 7% of the isolates, respectively. The aph(3')-Ia was found on different large conjugative plasmids belonging to IncI1α, which was present in 67.2% of the strains, on IncHI1, IncHI2, and IncN, as well as on a multicopy ColRNAI plasmid. All these plasmids except ColRNAI carried genes encoding resistance to other antimicrobials or heavy metals, highlighting the risk of co-selection. The aph(3')-Ib gene occurred on a 19 kb chimeric, mobilizable plasmid that contained elements tracing back its origin to distantly related genera. While aph(3')-Ia was flanked by either Tn903 or Tn4352 derivatives, no clear association was observed between aph(3')-Ib and mobile genetic elements. In conclusion, the spread of neomycin resistance in porcine clinical E. coli is driven by two resistance determinants located on distinct plasmid scaffolds circulating within a highly diverse population dominated by ST100. IMPORTANCE Neomycin is the first-choice antibiotic for the management of Escherichia coli enteritis in pigs. This work shows that aph(3')-Ia and to a lesser extent aph(3')-Ib are responsible for the spread of neomycin resistance that has been recently observed among pig clinical isolates and elucidates the mechanisms of dissemination of these two resistance determinants. The aph(3')-Ia gene is located on different conjugative plasmid scaffolds and is associated with two distinct transposable elements (Tn903 and Tn4352) that contributed to its spread. The diffusion of aph(3')-Ib is mediated by a small non-conjugative, mobilizable chimeric plasmid that likely derived from distantly related members of the Pseudomonadota phylum and was not associated with any detectable mobile genetic element. Although the spread of neomycin resistance is largely attributable to horizontal transfer, both resistance determinants have been acquired by a predominant lineage (ST100) associated with enterotoxigenic E. coli, which accounted for approximately one-third of the strains.

Multiplicity of blaKPC Genes and pKpQIL Plasmid Plasticity in the Development of Ceftazidime-Avibactam and Meropenem Coresistance in Klebsiella pneumoniae Sequence Type 307

In 2021, Klebsiella pneumoniae sequence type 307 (ST307) strains causing pulmonary and bloodstream infections identified in a hospital in Rome, Italy, reached high levels of resistance to ceftazidime-avibactam (CZA). One of these strains reached high levels of resistance to both CZA and carbapenems and carried two copies of blaKPC-3 and one copy of blaKPC-31 located on plasmid pKpQIL. The genomes and plasmids of CZA-resistant ST307 strains were analyzed to identify the molecular mechanisms leading to the evolution of resistance and compared with ST307 genomes at local and global levels. A complex pattern of multiple plasmids in rearranged configurations, coresident within the CZA-carbapenem-resistant K. pneumoniae strain, was observed. Characterization of these plasmids revealed recombination and segregation events explaining why K. pneumoniae isolates from the same patient had different antibiotic resistance profiles. This study illustrates the intense genetic plasticity occurring in ST307, one of the most worldwide-diffused K. pneumoniae high-risk clones.

The mitochondrial genome of the diploid oat Avena longiglumis

BACKGROUND

Avena longiglumis Durieu (2n = 2x = 14) is a wild relative of cultivated oat (Avena sativa, 2n = 6x = 42) with good agronomic and nutritional traits. The plant mitochondrial genome has a complex organization and carries genetic traits of value in exploiting genetic resources, not least male sterility alleles used to generate F₁ hybrid seeds. Therefore, we aim to complement the chromosomal-level nuclear and chloroplast genome assemblies of A. longiglumis with the complete assembly of the mitochondrial genome (mitogenome) based on Illumina and ONT long reads, comparing its structure with Poaceae species.

RESULTS

The complete mitochondrial genome of A. longiglumis can be represented by one master circular genome being 548,445 bp long with a GC content of 44.05%. It can be represented by linear or circular DNA molecules (isoforms or contigs), with multiple alternative configurations mediated by long (4,100-31,235 bp) and medium (144-792 bp) size repeats. Thirty-five unique protein-coding genes, three unique rRNA genes, and 11 unique tRNA genes are identified. The mitogenome is rich in duplications (up to 233 kb long) and multiple tandem or simple sequence repeats, together accounting for more than 42.5% of the total length. We identify homologous sequences between the mitochondrial, plastid and nuclear genomes, including the exchange of eight plastid-derived tRNA genes, and nuclear-derived retroelement fragments. At least 85% of the mitogenome is duplicated in the A. longiglumis nuclear genome. We identify 269 RNA editing sites in mitochondrial protein-coding genes including stop codons truncating ccmFC transcripts.

CONCLUSIONS

Comparative analysis with Poaceae species reveals the dynamic and ongoing evolutionary changes in mitochondrial genome structure and gene content. The complete mitochondrial genome of A. longiglumis completes the last link of the oat reference genome and lays the foundation for oat breeding and exploiting the biodiversity in the genus.

Whole-genome-based characterization of Campylobacter jejuni from human patients with gastroenteritis collected over an 18 year period reveals increasing prevalence of antimicrobial resistance

Campylobacteriosis is the most common cause of acute gastrointestinal bacterial infection in Europe, with most infections linked to the consumption of contaminated food. While previous studies found an increasing rate of antimicrobial resistance (AMR) in Campylobacter spp. over the past decades, the investigation of additional clinical isolates is likely to provide novel insights into the population structure and mechanisms of virulence and drug resistance of this important human pathogen. Therefore, we combined whole-genome sequencing and antimicrobial-susceptibility testing of 340 randomly selected Campylobacter jejuni isolates from humans with gastroenteritis, collected in Switzerland over an 18 year period. In our collection, the most common multilocus sequence types (STs) were ST-257 (n=44), ST-21 (n=36) and ST-50 (n=35); the most common clonal complexes (CCs) were CC-21 (n=102), CC-257 (n=49) and CC-48 (n=33). High heterogeneity was observed among STs, with the most abundant STs recurring over the entire study period, while others were observed only sporadically. Source attribution based on ST assigned more than half of the strains to the 'generalist' category (n=188), 25  % as 'poultry specialist' (n=83), and only a few to 'ruminant specialist' (n=11) or 'wild bird' origin (n=9). The isolates displayed an increased frequency of AMR from 2003 to 2020, with the highest rates of resistance observed for ciprofloxacin and nalidixic acid (49.8 %), followed by tetracycline (36.9 %). Quinolone-resistant isolates carried chromosomal gyrA mutations T86I (99.4 %) and T86A (0.6 %), whereas tetracycline-resistant isolates carried tet(O) (79.8 %) or mosaic tetO/32/O (20.2 %) genes. A novel chromosomal cassette carrying several resistance genes, including aph(3')-III, satA and aad(6), and flanked by insertion sequence elements was detected in one isolate. Collectively, our data revealed an increasing prevalence of resistance to quinolones and tetracycline in C. jejuni isolates from Swiss patients over time, linked to clonal expansion of gyrA mutants and acquisition of the tet(O) gene. Investigation of source attribution suggests that infections are most likely related to isolates from poultry or generalist backgrounds. These findings are relevant to guide future infection prevention and control strategies.

From Samples to Germline and Somatic Sequence Variation: A Focus on Next-Generation Sequencing in Melanoma Research

Next-generation sequencing (NGS) applications have flourished in the last decade, permitting the identification of cancer driver genes and profoundly expanding the possibilities of genomic studies of cancer, including melanoma. Here we aimed to present a technical review across many of the methodological approaches brought by the use of NGS applications with a focus on assessing germline and somatic sequence variation. We provide cautionary notes and discuss key technical details involved in library preparation, the most common problems with the samples, and guidance to circumvent them. We also provide an overview of the sequence-based methods for cancer genomics, exposing the pros and cons of targeted sequencing vs. exome or whole-genome sequencing (WGS), the fundamentals of the most common commercial platforms, and a comparison of throughputs and key applications. Details of the steps and the main software involved in the bioinformatics processing of the sequencing results, from preprocessing to variant prioritization and filtering, are also provided in the context of the full spectrum of genetic variation (SNVs, indels, CNVs, structural variation, and gene fusions). Finally, we put the emphasis on selected bioinformatic pipelines behind (a) short-read WGS identification of small germline and somatic variants, (b) detection of gene fusions from transcriptomes, and (c) de novo assembly of genomes from long-read WGS data. Overall, we provide comprehensive guidance across the main methodological procedures involved in obtaining sequencing results for the most common short- and long-read NGS platforms, highlighting key applications in melanoma research.

Factors Affecting the Quality of Bacterial Genomes Assemblies by Canu after Nanopore Sequencing

Long-read sequencing (LRS), like Oxford Nanopore Technologies, is usually associated with higher error rates compared to previous generations. Factors affecting the assembly quality are the integrity of DNA, the flowcell efficiency, and, not least all, the raw data processing. Among LRS-intended de novo assemblers, Canu is highly flexible, with its dozens of adjustable parameters. Different Canu parameters were compared for assembling reads of Salmonellaenterica ser. Bovismorbificans (genome size of 4.8 Mbp) from three runs on MinION (N50 651, 805, and 5573). Two of them, with low quality and highly fragmented DNA, were not usable alone for assembly, while they were successfully assembled when combining the reads from all experiments. The best results were obtained by modifying Canu parameters related to the error correction, such as corErrorRate (exclusion of overlaps above a set error rate, set up at 0.40), corMhapSensitivity (the coarse sensitivity level, set to “high”), corMinCoverage (set to 0 to correct all reads, regardless the overlaps length), and corOutCoverage (corrects the longest reads up to the imposed coverage, set to 100). This setting produced two contigs corresponding to the complete sequences of the chromosome and a plasmid. The overall results highlight the importance of a tailored bioinformatic analysis.