Whole Genome Sequencing Research Articles

Genomic Profile of a Multidrug-Resistant Klebsiella pneumoniae Strain Isolated from a Urine Specimen.

Klebsiella pneumoniaeis an opportunistic pathogen mostly found in health care-associated infections but can also be associated with community-acquired infections and is in critical need of new antimicrobial agents for strains resistant to carbapenems. The prevalence of carbapenemase-encoding genes varies among studies. Multidrug-resistant K. pneumoniae strains can harbor several antimicrobial-resistant determinants and mobile genetic elements (MGEs), along with virulence genetic determinants in community settings. We aim to determine the genetic profile of a multidrug-resistant K. pneumoniae strain isolated from a patient with community-acquired UTI. We isolated aK. pneumoniaestrain UABC-Str0120, from a urine sample of community-acquired urinary tract infection. Antimicrobial susceptibility tests and Whole-genome sequencing (WGS) were performed. The phylogenetic relationship was inferred by SNPs calling and filtering. UABC-Str0120 showed resistance toward β-lactams, combinations with β-lactamase inhibitors, and carbapenems. WGS revealed the presence of genes conferring resistance to aminoglycosides, β-lactams, carbapenems, quinolones, sulfonamides, phosphonates, phenicols, and quaternary ammonium compounds, 77 subsystems of virulence genes were identified, and an uncommon sequence type ST5889 was also determined. The sequenced strain harbors several MGEs. The UABC-Str0120 recovered from a urine sample harbors several virulence and antimicrobial resistance determinants, which assembles an endangering combination for an immunocompromised or a seemly healthy host, given its presence in a community setting.

Clinical and microbiological features of a cohort of patients with Acinetobacter baumannii bloodstream infections.

Acinetobacter baumannii is emerging as a pathogen that is a focus of global concern due to the frequent occurrence of the strains those are extensively resistant to antibiotics. This study was aimed to analyze the clinical and microbiological characteristics of a cohort of patients with A. baumannii bloodstream infections (BSIs) in western China. A retrospective study of the patients at West China Hospital of Sichuan University with A. baumannii BSIs between Jan, 2018 and May, 2023 was conducted. Antimicrobial susceptibility of A. baumannii isolates was tested by microdilution broth method. Whole-genome sequencing and genetic analysis were also performed for these isolates. Among the 117 patients included, longer intensive care unit stay, higher mortality, and more frequent invasive procedures and use of more than 3 classes of antibiotics were observed among the carbapenem-resistant A. baumannii (CRAB)-infected group (n = 76), compared to the carbapenem-susceptible A. baumannii (CSAB)-infected group (n = 41, all P ≤ 0.001). Twenty-four sequence types (STs) were determined for the117 isolates, and 98.7% (75/76) of CRAB were identified as ST2. Compared to non-ST2 isolates, ST2 isolates exhibited higher antibiotic resistance, and carried more resistance and virulence genes (P < 0.05). In addition, 80 (68.4%) isolates were CRISPR-positive, showed higher antibiotic susceptibility, and harbored less resistance and virulence genes, in comparison to CRISPR-negative ones (P < 0.05). Phylogenetic clustering based on coregenome SNPs indicated a sporadic occurrence of clonal transmission. Our findings demonstrate a high frequency of ST2 among A. baumannii causing BSIs, and high antibiotic susceptibility of non-ST2 and CRISPR-positive isolates. It is necessary to strengthen the surveillance of this pathogen.

Complete genome sequences of two tombusvirus-like viruses identified in Echinacea purpurea seeds.

Echinacea is an herbaceous plant originating from North America that is cultivated for gardening and landscaping because of its showy flowers. Using high-throughput sequencing, we identified two viral contigs from echinacea seeds that were related to the family Tombusviridae. These two viruses were similar to oat chlorotic stunt virus (OCSV) and other unassigned tombusviruses; therefore, we tentatively named them Echinacea-associated tombusviruses 1 and 2 (EaTV1 and EaTV2, respectively). The EaTVs represent putative readthrough sites and have no poly(A) tails, aligning with the common features of family Tombusviridae. The EaTVs are included in a monophyletic group of OCSV and several unassigned tombusviruses. Because OCSV is the only member of Avenavirus to date, EaTVs are tentative members of Avenavirus, or they are close sister species to OCSV with several unassigned tombusviruses. RNA-dependent RNA polymerases and coat proteins were well conserved among EaTVs and unassigned tombusviruses; however, their similarities were not correlated, implying divergent and complex evolution.

Complete genome sequence of a Pseudomonas fluorescens bacteriophage UNO-G1W1 isolated from freshwater ice in Nebraska.

We provide the complete genome sequence for a novel Pseudomonas fluorescens bacteriophage named UNO-G1W1. This phage was isolated from a single ice cover sampling. The genome was sequenced on the Nanopore MinION, generated with the direct terminal repeat-phage-pipeline and polished with Illumina short reads. Sequence identity classifies the phage as an otagovirus.

Complete genome sequence of Sphingomonas sp. strain R1, a bacterium isolated from tomato (Solanum lycopersicum).

Sphingomonas sp. strain R1 was isolated from the stem of a tomato plant and exhibited antagonism with Agrobacterium. The complete genome sequence of this bacterium consists of one 3,874,532 bp circular chromosome and two plasmids.

Complete genome sequences of Lactococcus petauri strains 473AN and 473GN isolated from the blood culture of a Japanese patient with infective endocarditis.

Here, we report the complete genome sequences of Lactococcus petauri strains 473AN and 473GN, isolated from the blood culture of a Japanese patient with infective endocarditis. The complete genomes of 473AN and 473GN consist of single chromosomes of 2,065,772 and 2,094,461 bp, respectively.

Improved complete genome sequence of Bdellovibrio bacteriovorus 109J, a widely studied laboratory strain of predatory bacteria.

The complete genome sequence of Bdellovibrio bacteriovorus 109J, a well-studied laboratory strain of predatory bacteria, first determined in 2014. Here we report an improved complete genome sequence of B. bacteriovorus 109J, incorporating 16 assembly and 87 nucleotide corrections. This revised genome will be helpful to studies on the predatory bacteria.

Association study of GBA1 variants with MSA based on comprehensive sequence analysis -Pitfalls in short-read sequence analysis depending on the human reference genome.

Multiple system atrophy (MSA) is a neurodegenerative disorder characterized by various combinations of autonomic failure, parkinsonism, and cerebellar ataxia. To elucidate variants associated with MSA, we have been conducting short-read-based whole-genome sequence analysis. In the process of the association studies, we initially focused on GBA1, a previously proposed susceptibility gene for MSA, to evaluate whether GBA1 variants can be efficiently identified despite its extraordinarily high homology with its pseudogene, GBA1LP. To accomplish this, we conducted a short-read whole-genome sequence analysis with alignment to GRCh38 as well as Sanger sequence analysis and compared the results. We identified five variants with inconsistencies between the two pipelines, of which three variants (p.L483P, p.A495P-p.V499V, p.L483_M489delinsW) were the results of misalignment due to minor alleles in GBA1P1 registered in GRCh38. The miscalling events in these variants were resolved by alignment to GRCh37 as the reference genome, where the major alleles are registered. In addition, a structural variant was not properly identified either by short-read or by Sanger sequence analyses. Having accomplished correct variant calling, we identified three variants pathogenic for Gaucher disease (p.S310G, p.L483P, and p.L483_M489delinsW). Of these variants, the allele frequency of p.L483P (0.003) in the MSA cases was higher than that (0.0011) in controls. The meta-analysis incorporating a previous report demonstrated a significant association of p.L483P with MSA with an odds ratio of 2.92 (95% CI; 1.08 - 7.90, p = 0.0353).

Congregibacter variabilis sp. nov. and Congregibacter brevis sp. nov. Within the OM60/NOR5 Clade, Isolated from Seawater, and Emended Description of the Genus Congregibacter.

Two Gram-stain-negative, aerobic, motile by means of flagella, short rod-shaped bacterial strains, designated IMCC43200T and IMCC45268T, were isolated from coastal seawater samples collected from the South Sea of Korea. Strains IMCC43200T and IMCC45268T shared 98.6% 16S rRNA gene sequence similarity and were closely related to Congregibacter litoralis KT71T (98.8% and 98.7%, respectively). Complete whole-genome sequences of IMCC43200T and IMCC45268T were 3.93 and 3.86 Mb in size with DNA G + C contents of 54.8% and 54.2%, respectively. Average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values between the two strains were 74.5% and 23.4%, respectively, revealing that they are independent species. The two strains showed ANI values of ≤ 75.8% and dDDH values of ≤ 23.0% to the type and only species of the genus Congregibacter (C. litoralis), indicating that each strain represents a novel species. Both strains contained summed feature 3 (comprising C16:1 ω6c and/or C16:1 ω7c) and summed feature 8 (comprising C18:1 ω6c and/or C18:1 ω7c) as major fatty acid constituents. The predominant isoprenoid quinone detected in both strains was ubiquinone-8 (Q-8). The major polar lipids of the two strains were phosphatidylethanolamine, phosphatidylglycerol, phospholipids, and aminolipids. Based on the phylogenetic, genomic, and phenotypic characterization, strains IMCC43200T and IMCC45268T were considered to represent two novel species within the genus Congregibacter, for which the names Congregibacter variabilis sp. nov. and Congregibacter brevis sp. nov. are proposed with IMCC43200T (= KCTC 8133T = NBRC 116295T = CCTCC AB 2023139T) and IMCC45268T (= KCTC 92921T = NBRC 116135T) as the type strains, respectively.

Whole-genome sequencing of Klebsiella pneumoniae isolated from clinical specimens in Chattogram, Bangladesh.

The emergence of multidrug-resistant Klebsiella pneumoniae (Kpn) is a global concern due to the increasing rate of mortality and hospital cost burden in the affected population. This study reports the whole-genome sequences of nine multidrug-resistant Kpn from a hospital in Chattogram city of Bangladesh.

Whole-genome sequencing of multidrug-resistant Escherichia coli causing urinary tract infection in an immunocompromised patient: a case report

BackgroundEscherichia coli is a major human pathogen responsible for a broad range of clinical illnesses. It has been linked to endemic and epidemic nosocomial diseases caused by multidrug-resistant pathogens in Sudan as well as throughout the globe.Case presentationA 76-year-old African woman arrived at Saad Rashwan Medical Centre complaining of backaches and discomfort during urination. Throughout the preceding 5 years, the patient had recurrent urinary tract infections. Following overnight incubation at 37 °C, Escherichia coli was found in her midstream urine specimen on cysteine lactose electrolyte deficient agar media. Minimum inhibitory concentration (colorimetric/turbidimetric method) was employed to test a wide range of antimicrobial drugs against this bacterial strain, and the results revealed significant multidrug resistance. QIAamp® DNA Mini Kit was used to obtain DNA Template from the purified Escherichia coli (Qiagen, Hilden, Germany). The bacterial whole-genome sequence was done by Novogene company (Hong Kong) using Illumina HiSeq 2500 (Illumina, San Diego, CA, USA), followed by whole genome reconstructions, and identification of antibiotic-resistant genes. Phylogenetic analysis revealed that our strain was related to the Escherichia coli DSM30083 ( genome sequence ID: CP033092.2) from the USA. Our strain possessed the following antimicrobial-resistant genes: aminoglycoside (kdpE, baeR, cpxA, aadA5), nitroimidazole (msbA), phosphonic acid (mdtG), tetracycline (emrY), macrolide, penam, tetracycline, (evgA, TolC, H-NS), fluoroquinolone, cephalosporin, glycylcycline, penam, tetracycline, rifamycin, phenicol antibiotic, disinfecting agents and antiseptics (acrB; marA), sulfonamide (sul1), macrolide (Mrx), cephalosporin, penam (CTX-M-15), carbapenem, cephalosporin, and penam (OXA-1).ConclusionThis study found that the isolated Escherichia coli strain had varied antimicrobial resistance genes on the basis of whole-genome sequencing and phenotypic resistance analyses. Whole-genome sequencing is critical for control and preventative methods to battle the growing threat of antimicrobial resistance. A larger investigation is recommended for improved generalization of results.

Can genomics and meteorology predict outbreaks of legionellosis in urban settings?

Legionella pneumophila is ubiquitous and sporadically infects humans causing Legionnaire's disease (LD). Globally, reported cases of LD have risen fourfold from 2000 to 2014. In 2016, Sydney, Australia was the epicenter of an outbreak caused by L. pneumophila serogroup 1 (Lpsg1). Whole-genome sequencing was instrumental in identifying the causal clone which was found in multiple locations across the city. This study examined the epidemiology of Lpsg1 in an urban environment, assessed typing schemes to classify resident clones, and investigated the association between local climate variables and LD outbreaks. Of 223 local Lpsg1 isolates, we identified dominant clones with one clone isolated from patients in high frequency during outbreak investigations. The core genome multi-locus sequence typing scheme was the most reliable in identifying this Lpsg1 clone. While an increase in humidity and rainfall was found to coincide with a rise in LD cases, the incidence of the major L. pneumophila outbreak clone did not link to weather phenomena. These findings demonstrated the role of high-resolution typing and weather context assessment in determining source attribution for LD outbreaks in urban settings, particularly when clinical isolates remain scarce.IMPORTANCEWe investigated the genomic and meteorological influences of infections caused by Legionella pneumophila in Sydney, Australia. Our study contributes to a knowledge gap of factors that drive outbreaks of legionellosis compared to sporadic infections in urban settings. In such cases, clinical isolates can be rare, and thus, other data are needed to inform decision-making around control measures. The study revealed that core genome multi-locus sequence typing is a reliable and adaptable technique when investigating Lpsg1 outbreaks. In Sydney, the genomic profile of Lpsg1 was dominated by a single clone, which was linked to numerous community cases over a period of 40 years. Interestingly, the peak in legionellosis cases during Autumn was not associated with this prevalent outbreak clone. Incorporating meteorological data with Lpsg1 genomics can support risk assessment strategies for legionellosis in urban environments, and this approach may be relevant for other densely populated regions globally.

Integrated analyses highlight interactions between the three-dimensional genome and DNA, RNA and epigenomic alterations in metastatic prostate cancer.

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

Complete genome sequences of three Prunus-infecting nepoviruses: apricot latent ringspot virus, myrobalan latent ringspot virus, and a novel virus from smooth stone peach (Prunus mira Koehne).

The complete genome sequences of two poorly studied Prunus-infecting nepoviruses, apricot latent ringspot virus (ALRSV) and myrobalan latent ringspot virus (MLRSV) were determined, confirming that they are members of subgroup C. Serological, biological, and molecular data, in particular a low level (58.8%) of amino acid sequence identity in the coat protein, suggest that ALRSV and MLRSV should be considered taxonomically distinct. In addition, data mining of public RNASeq data from wild and ornamental Prunus identified two contigs representing the nearly complete genome of a new subgroup A nepovirus from a smooth stone peach (Prunus mira) dataset (SRR8369794) from the Himalayas, for which the name "Prunus mira virus A" is proposed.

Insights from the first chromosome-level genome assembly of the alpine gentian Gentiana straminea Maxim. (Insight from genome assembly of Gentiana straminea).

Gentiana straminea Maxim. is a perennial herb and mainly distributed in the Qinghai-Tibetan Plateau. To adapt to the extreme environment, it has developed particular morphological, physiological and genetic structures. Also, rich in iridoids, it is one of the original plants of traditional Chinese herb "Qinjiao". Herein, we present its first chromosome-level genome sequence assembly, and compare it with the genomes of other Gentiana species to facilitate the analysis of genomic characteristics. The assembled genome size of G. straminea was 1.25 Gb, with a contig N50 of 7.5Mb. A total of 96.08% of the genome sequences was anchored on 13 pseudochromosomes, with a scaffold N50 of 92.70Mb. A total of 54,310 protein-coding genes were predicted, 80.25% of which were functionally annotated. Comparative genomic analyses indicated that G. straminea experienced two whole-genome duplication events after the γ whole-genome triplication with other eudicots, and it diverged from other Gentiana species at ~3.2 Mya. A total of 142 enzyme-coding genes related to iridoid biosynthesis were identified in its genome. Additionally, we identified differences in the number and expression patterns of iridoid biosynthetic pathway genes in G. straminea compared with two other Gentiana species by integrating whole-genome sequence and transcriptomic analyses.

Correction for de Castro et al., "Draft Whole-Genome Sequence of Bacillus paramycoides LB_RP2, a Putative Polyhydroxyalkanoate-Producing Bacterium Isolated from an Amazonian Blackwater River".

Correction for de Castro et al., "Draft Whole-Genome Sequence of Bacillus paramycoides LB_RP2, a Putative Polyhydroxyalkanoate-Producing Bacterium Isolated from an Amazonian Blackwater River".

Environmental memory alters the fitness effects of adaptive mutations in fluctuating environments.

Evolution in a static laboratory environment often proceeds via large-effect beneficial mutations that may become maladaptive in other environments. Conversely, natural settings require populations to endure environmental fluctuations. A sensible assumption is that the fitness of a lineage in a fluctuating environment is the time average of its fitness over the sequence of static conditions it encounters. However, transitions between conditions may pose entirely new challenges, which could cause deviations from this time average. To test this, we tracked hundreds of thousands of barcoded yeast lineages evolving in static and fluctuating conditions and subsequently isolated 900 mutants for pooled fitness assays in 15 environments. Here we find that fitness in fluctuating environments indeed often deviates from the time average, leading to fitness non-additivity. Moreover, closer examination reveals that fitness in one component of a fluctuating environment is often strongly influenced by the previous component. We show that this environmental memory is especially common for mutants with high variance in fitness across tested environments. We use a simple mathematical model and whole-genome sequencing to propose mechanisms underlying this effect, including lag time evolution and sensing mutations. Our results show that environmental fluctuations impact fitness and suggest that variance in static environments can explain these impacts.

Genomic characterisation of generic Escherichia coli from food-producing animals and products of animal origin in South Africa

Escherichia coli is found in diverse environmental niches, including meat and meat products, and is known for its significance in both food safety and public health. In South Africa, whole genomic sequencing (WGS) efforts for E. coli are overwhelmingly skewed toward human isolates with limited studies conducted on non-human isolates. Therefore, the aim of this study was to use WGS to characterise generic strains of E. coli isolated from animal specimens, meat, and meat-based products in South Africa. Based on WGS analysis, a total of 35 E. coli strains were grouped into five phylogroups (A, B1, B2, C, and E), with A (46%) being the most predominant. Virulence-associated genes identified the isolates as either extra-intestinal pathogenic E. coli (69%) or intestinal pathogenic E. coli (31%) pathotypes. Twenty-three different serotypes were identified, with O101:H37 (17%), O2:H4 (17%), O6:H5 (13%), and O64:H19 (9%) being the predominant ones. Among the 19 different sequence types (STs), ST1858, ST975, and ST10 were the most prevalent (11% each). Various virulence genes, antimicrobial resistance genes, and genetic mobile elements carrying Tn2, IS26, and Tn6196 elements were detected, with the disinfectant resistance sitABCD being the most predominant. The type 1 CRISPR system which functions by storing records of previous invasions to provide immunological memory for a rapid and robust response upon subsequent viral infections was detected in all isolates, consisting of subtypes I-E (86%), I-A (57%), and I-F (11%). The findings of this study provide an insight into the genetic diversity of generic E. coli isolates from animal species, meat, and meat-based products in South Africa.

Genomic diversity and analysis of resistance determinants of &lt;i&gt;Salmonella enterica&lt;/i&gt; subspecies &lt;i&gt;enterica&lt;/i&gt; serotype Kentucky isolated in Russia

Introduction. Salmonella Kentucky sequence type ST198 is one of the epidemiologically significant non-typhoidal Salmonella clones worldwide and is characterized by the presence of highly resistant strains and the ability to adapt to different animal hosts and environmental conditions. The aim of this study was to analyze S. Kentucky strains isolated from various sources in Russia in terms of their phylogenetic position within the global diversity of the pathogen and their genetic characteristics. Materials and methods. We examined 55 strains of S. Kentucky by whole-genome sequencing, which were isolated from 2010 to 2022 from various sources (clinical strains, food, as well as from farm animals, feed and environmental samples). Whole genome sequencing was performed using Illumina platforms. Phylogenetic analysis based on nucleotide variation analysis included an additional 390 S. Kentucky strains. Results. Most of the Russian strains (n = 50) belonged to the ST198 sequence type, four strains were ST314 and one strain was ST152. Of the 50 Russian sequence-type ST198 strains, 44 belonged to the international monophyletic MDR lineage S. Kentucky ST198, and belonged to four separate sublineages, six strains occupying a basal position in relation to the MDR lineage. A total of 320 genes and mutations responsible for resistance to antimicrobial agents were identified. The most common were point mutations in the QRDR region. In most cases, Russian strains were characterized by the presence of variants of the SGI1-K genomic island. Moreover, the putative structure of SGI1 was correlated with the phylogenetic clustering of S. Kentucky sublineages. Conclusions. The results of the study made it possible to assess the population structure of Russian S. Kentucky ST198 strains on a global scale and determine the genetic determinants of antibiotic resistance, including the structure of the SGI1 genomic island.

AcornHRD: an HRD algorithm highly associated with anthracycline-based neoadjuvant chemotherapy in breast cancer in China

PurposeOur study aimed to develop and validate a homologous recombination deficiency (HRD) scoring algorithm in the Chinese breast cancer population.Methods and materialsNinety-six in-house breast cancer (BC) samples and 6 HRD-positive standard cells were analyzed by whole-genome sequencing (WGS). Besides, 122 BCs from the TCGA database were down-sampled to ~ 1X WGS. We constructed an algorithm named AcornHRD for HRD score calculated based on WGS at low coverage as input data to estimate large-scale copy number alteration (LCNA) events on the genome. A clinical cohort of 50 BCs (15 cases carrying BRCA mutation) was used to assess the association between HRD status and anthracyclines-based neoadjuvant treatment outcomes.ResultsA 100-kb window was defined as the optimal size using 41 in-house cases and the TCGA dataset. HRD score high threshold was determined as HRD score ≥ 10 using 55 in-house BCs with BRCA mutation to achieve a 95% BRCA-positive agreement rate. Furthermore, the HRD status agreement rate of AcornHRD is 100%, while the ShallowHRD is 60% in standard cells. BRCA mutation was significantly associated with a high HRD score evaluated by AcornHRD and ShallowHRD (p = 0.008 and p = 0.003, respectively) in the TCGA dataset. However, AcornHRD showed a higher positive agreement rate than did the ShallowHRD algorithm (70% vs 60%). In addition, the BRCA-positive agreement rate of AcornHRD was superior to that of ShallowHRD (87% vs 13%) in the clinical cohort. Importantly, the high HRD score assessed by AcornHRD was significantly correlated with a residual cancer burden score of 0 or 1 (RCB0/1). Besides, the HRD-positive group was more likely to respond to anthracycline-based chemotherapy than the HRD-negative group (pCR [OR = 9.5, 95% CI 1.11–81.5, p = 0.040] and RCB0/1 [OR = 10.29, 95% CI 2.02–52.36, p = 0.005]).ConclusionUsing the AcornHRD algorithm evaluation, our analysis demonstrated the high performance of the LCNA genomic signature for HRD detection in breast cancers.