Whole Genome Sequencing Research Articles

Whole genome sequencing enhances molecular diagnosis of primary ciliary dyskinesia.

Primary ciliary dyskinesia (PCD) is a genetic disorder affecting motile cilia. Most cases are inherited recessively, due to variants in >50 genes that result in abnormal or absent motile cilia. This leads to chronic upper and lower airway disease, subfertility, and laterality defects. Given overlapping clinical features and genetic heterogeneity, diagnosis can be difficult and often occurs late. Of those tested an estimated 30% of genetically screened PCD patients still lack a molecular diagnosis. A molecular diagnosis allows for appropriate clinical management including prediction of phenotypic features correlated to genotype. Here, we aimed to identify how readily a genetic diagnosis could be made using whole genome sequencing (WGS) to facilitate identification of pathogenic variants in known genes as well as novel PCD candidate genes. WGS was used to screen for pathogenic variants in eight patients with PCD. 7/8 cases had homozygous or biallelic variants in DNAH5, DNAAF4 or DNAH11 classified as pathogenic or likely pathogenic. Three identified variants were deletions, ranging from 3 to 13 kb, for which WGS identified precise breakpoints, permitting confirmation by Sanger sequencing. WGS yielded identification of a de novo variant in a novel PCD gene TUBB4B. Here, WGS uplifted genetic diagnosis of PCD by identifying structural variants and novel modes of inheritance in new candidate genes. WGS could be an important component of the PCD diagnostic toolkit, increasing molecular diagnostic yield from current (70%) levels, and enhancing our understanding of fundamental biology of motile cilia and variants in the noncoding genome.

Size-selective harvesting drives genomic shifts in a harvested population.

Overfishing not only drastically reduces the number of fish in an exploited population but is also often selective for body size by removing the largest individuals from a population. Here, we study experimentally the evolutionary effects of size-selective harvesting using whole-genome sequencing on a model organism, the zebrafish (Danio rerio). We demonstrate genomic shifts in the populations exposed to size-selective harvesting for five generations and show reduced genetic diversity in all harvested lines, including the control line (non-size-selected). We also determine differences in groups of genes related to certain gene ontology annotations between size-selectively harvested lines, with enrichment in nervous system related genes in the large-selected lines. Our results illuminate the biological processes underlying fisheries-induced genetic changes and hence contribute toward the understanding of the changes potentially associated with the vulnerability of an exploited population to future stressors.

Multiplexed real-time PCR for the detection and differentiation of Klebsiella pneumoniae O-antigen serotypes.

Klebsiella pneumoniae has emerged as a global health threat due to its role in the spread of antimicrobial resistance and because it is a frequent cause of hospital-acquired infections and neonatal sepsis. Capsular and lipopolysaccharide (LPS) O-antigen polysaccharide surface antigens are major immunogens that are useful for strain classification and are candidates for vaccine development. We have developed real-time PCR reagents for molecular serotyping, subtyping, and quantitation of the most prevalent LPS O-antigen types (i.e., O1, O2, O3, and O5) of Klebsiella pneumoniae. We describe two applications for this O-typing assay: for screening culture isolates and for direct typing of Klebsiella pneumoniae present in stool samples. We find 100% concordance between the results of the O-typing assay and whole-genome sequencing of 81 culture isolates, and >90% agreement in O-typing performed directly on specimens of human stool, with disagreement arising primarily from a lack of sensitivity of the culture-based comparator method. Additionally, we find evidence for mixed O-type populations at varying levels of abundance in direct tests of stool from a hospitalized patient population. Taken together, these results demonstrate that this novel O-typing assay can be a useful tool for K. pneumoniae epidemiologic and vaccine studies.IMPORTANCEKlebsiella pneumoniae is an important opportunistic pathogen. The gastrointestinal (GI) tract is the primary reservoir of K. pneumoniae in humans, and GI carriage is believed to be a prerequisite for invasive infection. Knowledge about the dynamics and duration of GI carriage has been hampered by the lack of tools suitable for detection and strain discrimination. Real-time PCR is particularly suited to the higher-throughput workflows used in population-based studies, which are needed to improve our understanding of carriage dynamics and the factors influencing K. pneumoniae colonization.

Transition to Whole Genome Sequencing Surveillance: The Impact on National Outbreak Detection and Response for Listeria monocytogenes, Salmonella, Shiga Toxin-Producing Escherichia coli, and Shigella Clusters in Canada, 2015-2021.

Between 2017 and 2019, pulsed-field gel electrophoresis was replaced by whole genome sequencing (WGS) for identifying enteric disease clusters in Canada. The number and characteristics of all clusters of Listeria monocytogenes, Salmonella, Shiga toxin-producing Escherichia coli (STEC), and Shigella spp. between 2015 and 2021 were analyzed. Following the transition to WGS, an increase in the number of Salmonella, STEC, and Shigella clusters was noted, whereas the number of clusters of L. monocytogenes decreased. Unlike previous subtyping methods, WGS provided increased resolution to identify discrete clusters of Salmonella Enteritidis. This led to the identification of a number of outbreaks linked to frozen raw breaded chicken products and ultimately a change in food safety policy to reduce the number of illnesses associated with these products. Other pathogens did not experience a similar increase in the number of outbreaks detected. Although WGS did provide increased confidence in the genetic relatedness of cases and isolates, challenges remained in collecting epidemiological data to link these illnesses to a common source.

The genomic landscape of 2,023 colorectal cancers.

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

Burkholderia cenocepacia ST-250 in cystic fibrosis patients in Switzerland: Genomic investigation of transmission routes

This report describes the characterization of Burkholderia cenocepacia isolates belonging to sequence type (ST)-250, detected in eight patients with cystic fibrosis (CF) in Switzerland. We retrospectively analyzed 18 isolates of B. cenocepacia ST-250 isolated between 2003 and 2015 by whole-genome sequencing and evaluated clinical and epidemiological data. Single nucleotide polymorphism analysis of the B.°cenocepacia ST-250 lineage showed that the isolates from all patients cluster tightly, suggesting that this cluster has a recent common ancestor. Epidemiological investigations showed that six out of eight patients acquired B.°cenocepacia ST-250 in the years 2001-2006, where participation in CF summer camps was common. Two patients were siblings. Genomic relatedness of the B. cenocepacia ST-250 isolates supported transmission by close contact, however, a common source or nosocomial routes cannot be excluded. With respect to the fatal outcome in six patients, our study shows the importance of infection control measurements in CF patients with B.°cenocepacia.

Identification of differentially expressed tumour-related genes regulated by UHRF1-driven DNA methylation

Ubiquitin-like with PHD and RING finger domains 1 (UHRF1) is an epigenetic regulator that plays critical roles in tumours. However, the DNA methylation alteration patterns driven by UHRF1 and the related differentially expressed tumour-related genes remain unclear. In this study, a UHRF1-shRNA MCF-7 cell line was constructed, and whole-genome bisulfite sequencing and RNA sequencing were performed. The DNA methylation alteration landscape was elucidated, and DNA methylation-altered regions (DMRs) were found to be distributed in both gene bodies and adjacent regions. The DMRs were annotated and categorized into 488 hypermethylated/1696 hypomethylated promoters and 1149 hypermethylated/5501 hypomethylated gene bodies. Through an integrated analysis with the RNA sequencing data, 217 methylation-regulated upregulated genes and 288 downregulated genes were identified, and these genes were primarily enriched in nervous system development and cancer signalling pathways. Further analysis revealed 21 downregulated oncogenes and 15 upregulated TSGs. We also showed that UHRF1 silencing inhibited cell proliferation and migration and suppressed tumour growth in vivo. Our study suggested that UHRF1 and the oncogenes or TSGs it regulates might serve as biomarkers and targets for breast cancer treatment.

Whole-genome sequencing for antimicrobial surveillance: species-specific quality thresholds and data evaluation from the network of the European Union Reference Laboratory for Antimicrobial Resistance genomic proficiency tests of 2021 and 2022.

As antimicrobial resistance (AMR) surveillance shifts to genomics, ensuring the quality of whole-genome sequencing (WGS) data produced across laboratories is critical. Participation in genomic proficiency tests (GPTs) not only increases individual laboratories' WGS capacity but also provides a unique opportunity to improve species-specific thresholds for WGS quality control (QC) by repeated resequencing of distinct isolates. Here, we present the results of the EU Reference Laboratory for Antimicrobial Resistance (EURL-AR) network GPTs of 2021 and 2022, which included 25 EU national reference laboratories (NLRs). A total of 392 genomes from 12 AMR-bacteria were evaluated based on WGS QC metrics. Two percent (n = 9) of the data were excluded, due to contamination, and 11% (n = 41) of the remaining genomes were identified as outliers in at least one QC metric and excluded from computation of the adjusted QC thresholds (AQT). Two QC metric correlation groups were identified through linear regression. Eight percent (n = 28) of the submitted genomes, from 11 laboratories, failed one or more of the AQTs. However, only three laboratories (12%) were identified as underperformers, failing across AQTs for uncorrelated QC metrics in at least two genomes. Finally, new species-specific thresholds for "N50" and "number of contigs > 200 bp" are presented for guidance in routine laboratory QC. The continued participation of NRLs in GPTs will reveal WGS workflow flaws and improve AMR surveillance data. GPT data will continue to contribute to the development of reliable species-specific thresholds for routine WGS QC, standardizing sequencing data QC and ensure inter- and intranational laboratory comparability.IMPORTANCEIllumina next-generation sequencing is an integral part of antimicrobial resistance (AMR) surveillance and the most widely used whole-genome sequencing (WGS) platform. The high-throughput, relative low-cost, high discriminatory power, and rapid turnaround time of WGS compared to classical biochemical methods means the technology will likely remain a fundamental tool in AMR surveillance and public health. In this study, we present the current level of WGS capacity among national reference laboratories in the EU Reference Laboratory for AMR network, summarizing applied methodology and statistically evaluating the quality of the obtained sequence data. These findings provide the basis for setting new and revised thresholds for quality metrics used in routine WGS, which have previously been arbitrarily defined. In addition, underperforming participants are identified and encouraged to evaluate their workflows to produce reliable results.

Whole-exome and whole-genome sequencing of 1064 individuals with type 1 diabetes reveals novel genes for diabetic kidney disease.

Diabetic kidney disease (DKD) is a severe diabetic complication that affects one third of individuals with type 1 diabetes. Although several genes and common variants have been shown to be associated with DKD, much of the predicted inheritance remains unexplained. Here, we performed next-generation sequencing to assess whether low-frequency variants, extending to a minor allele frequency (MAF) ≤10% (single or aggregated) contribute to the missing heritability in DKD. We performed whole-exome sequencing (WES) of 498 individuals and whole-genome sequencing (WGS) of 599 individuals with type 1 diabetes. After quality control, next-generation sequencing data were available for a total of 1064 individuals, of whom 541 had developed either severe albuminuria or end-stage kidney disease, and 523 had retained normal albumin excretion despite a long duration of type 1 diabetes. Single-variant and gene-aggregate tests for protein-altering variants (PAV) and protein-truncating variants (PTV) were performed separately for WES and WGS data and combined in a meta-analysis. We also performed genome-wide aggregate analyses on genomic windows (sliding window), promoters and enhancers using the WGS dataset. In the single-variant meta-analysis, no variant reached genome-wide significance, but a suggestively associated common THAP7 rs369250 variant (p=1.50×10-5, MAF=49%) was replicated in the FinnGen general population genome-wide association study (GWAS) data for chronic kidney disease and DKD phenotypes. The gene-aggregate meta-analysis provided suggestive evidence (p<4.0×10-4) at four genes for DKD, of which NAT16 (MAFPAV≤10%) and LTA (also known as TNFβ, MAFPAV≤5%) are replicated in the FinnGen general population GWAS data. The LTA rs2229092 Callele was associated with significantly lower TNFR1, TNFR2 and TNFR3 serum levels in a subset of FinnDiane participants. Of the intergenic regions suggestively associated with DKD, the enhancer on chromosome 18q12.3 (p=3.94×10-5, MAFvariants≤5%) showed interaction with the METTL4 gene; the lead variant was replicated, and predicted to alter binding of the MafB transcription factor. Our sequencing-based meta-analysis revealed multiple genes, variants and regulatory regions that were suggestively associated with DKD. However, as no variant or gene reached genome-wide significance, further studies are needed to validate the findings.

Complete genome sequence of a novel mitovirus isolated from the phytopathogenic fungus Alternaria alternata causing apple leaf blotch.

In this study, a novel mitovirus, tentatively designated as "Alternaria alternata mitovirus 2" (AaMV2), was isolated from the fungus Alternaria alternata f. sp. mali causing apple leaf blotch disease. The complete genome of AaMV2 is 3,157 nucleotides in length, with an A+U content of 68.10%. The genome has a single large open reading frame (ORF) encoding an RNA-dependent RNA polymerase (RdRp) protein with a molecular mass of 98.10 kDa. BLAST analysis revealed that AaMV2 has the highest sequence identity to Leptosphaeria biglobosa mitovirus 6, with 79.76% and 82.86% identity at the amino acid and nucleotide level, respectively. Phylogenetic analysis suggested that AaMV2 is a new member of the genus Duamitovirus within the family Mitoviridae. This is the first report of the complete genome sequence analysis of a mitovirus in A. alternata.

Molecular epidemiology of string test-positive Klebsiella pneumoniae isolates in Huzhou, China, 2020-2023

ObjectiveThis study used whole-genome sequencing (WGS) to explore the genetic diversity, virulence factors, and antimicrobial resistance determinants of string test-positive Klebsiella pneumoniae (KP) over a 4-year surveillance period in Huzhou, China.MethodsIn total, 632 clinical isolates were collected via hospital surveillance from 2020 to 2023; 100 were positive in the string test and these 100 strains were subjected to antimicrobial susceptibility testing using an agar dilution method followed by WGS.ResultsThe resistance rates to cefotaxime (77.0%), trimethoprim-sulfamethoxazole (67.0%), and nalidixic acid (64.0%) were high. Multilocus sequence typing revealed high genetic diversity; there were 33 sequence types (STs) and 15 capsular serotypes. The most common ST was ST23 (16.0%) and the most common capsular serotype was K1 (22.5%). Virulome analysis revealed among-strain differences in virulence factors that affected bacterial adherence, efflux pump action, iron uptake, nutritional factors, metabolic regulation, the secretion system, and toxin production. The Kleborate strain-specific virulence scores of all 100 string test-positive KPs were derived: 28 strains scored 5, 28 scored 4, 21 scored 3, 12 scored 1, and 11 scored 0. All 77 strains with scores of 3 to 5 contained the iucA gene. The phylogeny based on whole-genome single nucleotide polymorphisms (wgSNPs) indicated high clonality; the string test-positive KP strains were grouped into six clades. Closely related isolates in each genetic cluster usually shared STs.ConclusionThe present study highlights the significance of the KP iucA gene in terms of hypervirulence and the diverse genotypes of string test-positive KP strains isolated in Huzhou hospitals.

A Systems Biology Analysis of Chronic Lymphocytic Leukemia

Whole-genome sequencing has revealed that TP53, NOTCH1, ATM, SF3B1, BIRC3, ABL, NXF1, BCR, and ZAP70 are often mutated in CLL, but not consistently across all CLL patients. This paper employs a statistical thermodynamics approach in combination with the systems biology of the CLL protein–protein interaction networks to identify the most significant participant proteins in the cancerous transformation. Betti number (a topology of complexity) estimates highlight a protein hierarchy, primarily in the Wnt pathway known for aberrant CLL activation. These individually identified proteins suggest a network-targeted strategy over single-target drug development. The findings advocate for a multi-target inhibition approach, limited to several key proteins to minimize side effects, thereby providing a foundation for designing therapies. This study emphasizes a shift towards a comprehensive, multi-scale analysis to enhance personalized treatment strategies for CLL, which could be experimentally validated using siRNA or small-molecule inhibitors. The result is not just the identification of these proteins but their rank-order, offering a potent signal amplification in the context of the 20,000 proteins produced by the human body, thus providing a strategic basis for therapeutic intervention in CLL, underscoring the necessity for a more holistic, cellular, chromosomal, and genome-wide study to develop tailored treatments for CLL patients.

Micro-costing of genetic diagnostics in acute leukemia in Sweden: from standard-of-care to whole-genome sequencing

ABSTRACT Aims and background Whole-genome sequencing (WGS) is increasingly applied in clinical practice and expected to replace standard-of-care (SoC) diagnostics in hematological malignancies. This study aims to assess and compare the fully burdened cost (‘micro-costing’) per patient for Swedish laboratories using WGS and SoC, respectively, in pediatric and adult patients with acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML). Methods The resource use and cost details associated with SoC, e.g. chromosome banding analysis, fluorescent in situ hybridization, and targeted sequencing analysis, were collected via activity-based costing methods from four diagnostic laboratories. For WGS, corresponding data was collected from two of the centers. A simulation-based scenario model was developed for analyzing the WGS cost based on different annual sample throughput to evaluate economy of scale. Results The average SoC total cost per patient was €2,465 for pediatric AML and €2,201 for pediatric ALL, while in adults, the corresponding cost was €2,458 for AML and €1,207 for ALL. The average WGS cost (90x tumor/30x normal; sequenced on the Illumina NovaSeq 6000 platform) was estimated to €3,472 based on an annual throughput of 2,500 analyses, however, with an annual volume of 7,500 analyses the average cost would decrease by 23% to €2,671. Conclusion In summary, WGS is currently more costly than SoC, however the cost can be reduced by utilizing laboratories with higher throughput and by the expected decline in cost of reagents. Our data provides guidance to decision-makers for the resource allocation needed when implementing WGS in diagnostics of hematological malignancies.

Multidisciplinary Effort Leading to Effective Tuberculosis Community Outbreak Containment in Israel

Tuberculosis (TB) is the second-most prevalent cause of mortality resulting from infectious diseases worldwide. It is caused by bacteria belonging to the Mycobacterium tuberculosis complex (MTBC). In Israel, TB incidence is low, acknowledged by the WHO as being in a pre-elimination phase. Most cases occur among immigrants from high TB incidence regions like the Horn of Africa and the former Soviet Union (FSU), with occasional outbreaks. The outbreak described in this report occurred between 2018 and 2024, increasing the incidence rate of TB in the region. Control of this outbreak posed challenges due to factors including a diverse population (including Ethiopian immigrants, Israeli-born citizens, and immigrants from other countries), economic and social barriers, and hesitancy to disclose information. The unique multidisciplinary team formed to address these challenges, involving the local TB clinic, district health ministry, health maintenance organization (HMO) infectious disease consultant, neighborhood clinic, and National Mycobacterium Reference Laboratory (NMRL), achieved effective treatment and containment. Whole genome sequencing (WGS) proved pivotal in unraveling patient connections during the outbreak. It pinpointed those patients overlooked in initial field investigations, established connections between patients across different health departments, and uncovered the existence of two distinct clusters with separate transmission chains within the same neighborhood. This study underscores collaborative efforts across sectors that successfully contained a challenging outbreak.

Genetic Characterization of Antibiotic-Resistant Staphylococcus spp. and Mammaliicoccus sciuri from Healthy Humans and Poultry in Nigeria

Staphylococcus spp. poses a significant threat to human and animal health due to their capacity to cause a wide range of infections in both. In this study, resistance genes conferring antibiotic resistance in Staphylococcus spp. and Mammaliicoccus sciuri isolates from humans and poultry in Edo state, Nigeria, were investigated. In April 2017, 61 Staphylococcus spp. isolates were obtained from urine, wounds, nasal and chicken fecal samples. Species identification was carried out by matrix-assisted laser desorption ionization-time of flight mass spectrometry. Antimicrobial susceptibility testing was performed using the Kirby-Bauer method for 16 antibiotics. Whole-genome sequencing was used for characterization of the isolates. The 61 investigated isolates included Staphylococcus aureus, S. arlettae, M. sciuri, S. haemolyticus, and S. epidermidis. A total of 47 isolates (77%) belonged to human samples and 14 (23%) isolates were collected from poultry samples. All were phenotypically resistant to at least three antimicrobial(s). Multiple resistance determinants were detected in the human and poultry isolates analyzed. Phylogenetic analysis revealed close relatedness among the isolates within each species for S. arlettae, M. sciuri, and S. haemolyticus, respectively. This study delivered comprehensive genomic insights into antibiotic-resistant Staphylococcus species and M. sciuri isolates from human and poultry sources in Edo state, Nigeria, from a One Health perspective.

Plasmodium falciparum population dynamics in East Africa and genomic surveillance along the Kenya-Uganda border

East African countries accounted for ~ 10% of all malaria prevalence worldwide in 2022, with an estimated 23.8 million cases and > 53,000 deaths. Despite recent increases in malaria incidence, high-resolution genome-wide analyses of Plasmodium parasite populations are sparse in Kenya, Tanzania, and Uganda. The Kenyan-Ugandan border region is a particular concern, with Uganda confirming the emergence and spread of artemisinin resistant P. falciparum parasites. To establish genomic surveillance along the Kenyan-Ugandan border and analyse P. falciparum population dynamics within East Africa, we generated whole-genome sequencing (WGS) data for 38 parasites from Bungoma, Western Kenya. These sequences were integrated into a genomic analysis of available East African isolate data (n = 599) and revealed parasite subpopulations with distinct genetic structure and diverse ancestral origins. Ancestral admixture analysis of these subpopulations alongside isolates from across Africa (n = 365) suggested potential independent ancestral populations from other major African populations. Within isolates from Western Kenya, the prevalence of biomarkers associated with chloroquine resistance (e.g. Pfcrt K76T) were significantly reduced compared to wider East African populations and a single isolate contained the PfK13 V568I variant, potentially linked to reduced susceptibility to artemisinin. Overall, our work provides baseline WGS data and analysis for future malaria genomic surveillance in the region.

Quantitative drug susceptibility testing for Mycobacterium tuberculosis using unassembled sequencing data and machine learning.

There remains a clinical need for better approaches to rapid drug susceptibility testing in view of the increasing burden of multidrug resistant tuberculosis. Binary susceptibility phenotypes only capture changes in minimum inhibitory concentration when these cross the critical concentration, even though other changes may be clinically relevant. We developed a machine learning system to predict minimum inhibitory concentration from unassembled whole-genome sequencing data for 13 anti-tuberculosis drugs. We trained, validated and tested the system on 10,859 isolates from the CRyPTIC dataset. Essential agreement rates (predicted MIC within one doubling dilution of observed MIC) were above 92% for first-line drugs, 91% for fluoroquinolones and aminoglycosides, and 90% for new and repurposed drugs, albeit with a significant drop in performance for the very few phenotypically resistant isolates in the latter group. To further validate the model in the absence of external MIC datasets, we predicted MIC and converted values to binary for an external set of 15,239 isolates with binary phenotypes, and compare their performance against a previously validated mutation catalogue, the expected performance of existing molecular assays, and World Health Organization Target Product Profiles. The sensitivity of the model on the external dataset was greater than 90% for all drugs except ethionamide, clofazimine and linezolid. Specificity was greater than 95% for all drugs except ethambutol, ethionamide, bedaquiline, delamanid and clofazimine. The proposed system can provide quantitative susceptibility phenotyping to help guide antimicrobial therapy, although further data collection and validation are required before machine learning can be used clinically for all drugs.

Variation in the spectrum of new mutations among inbred strains of mice.

The mouse serves as a mammalian model for understanding the nature of variation from new mutations, a question that has both evolutionary and medical significance. Previous studies suggest that the rate of single nucleotide mutations (SNMs) in mice is approximately 50% of that in humans. However, information largely comes from studies involving the C57BL/6 strain, and there is little information from other mouse strains. Here, we study the mutations that accumulated in 59 mouse lines derived from four inbred strains that are commonly used in genetics and clinical research (BALB/cAnNRj, C57BL/6JRj, C3H/HeNRj, and FVB/NRj), maintained for 8-9 generations by brother-sister mating. By analysing Illumina whole-genome sequencing data, we estimate that the average rate of new SNMs in mice is approximately μ = 6.7 × 10-9. However, there is substantial variation in the spectrum of SNMs among strains, so the burden from new mutations also varies among strains. For example, the FVB strain has a spectrum that is markedly skewed towards C→A transversions, and is likely to experience a higher deleterious load than other strains, due to an increased frequency of nonsense mutations in glutamic acid codons. Finally, we observe substantial variation in the rate of new SNMs among DNA sequence contexts, CpG sites and their adjacent nucleotides playing an important role.

Inferring enterovirus D68 transmission dynamics from the genomic data of two 2022 North American outbreaks

Enterovirus D68 (EV-D68) has emerged as a significant cause of acute respiratory illness in children globally, notably following its extensive outbreak in North America in 2014. A recent outbreak of EV-D68 was observed in Ontario, Canada, from August to October 2022. Our phylogenetic analysis revealed a notable genetic similarity between the Ontario outbreak and a concurrent outbreak in Maryland, USA. Utilizing Bayesian phylodynamic modeling on whole genome sequences (WGS) from both outbreaks, we determined the median peak time-varying reproduction number (Rt) to be 2.70, 95% HPD (1.76, 4.08) in Ontario and 2.10, 95% HPD (1.41, 3.17) in Maryland. The Rt trends in Ontario closely matched those derived via EpiEstim using reported case numbers. Our study also provides new insights into the median infection duration of EV-D68, estimated at 7.94 days, 95% HPD (4.55, 12.8) in Ontario and 10.8 days, 95% HPD (5.85, 18.6) in Maryland, addressing the gap in the existing literature surrounding EV-D68’s infection period. We observed that the estimated Time since the Most Recent Common Ancestor (TMRCA) and the epidemic’s origin coincided with the easing of COVID-19 related social contact restrictions in both areas. This suggests that the relaxation of non-pharmaceutical interventions, initially implemented to control COVID-19, may have inadvertently facilitated the spread of EV-D68. These findings underscore the effectiveness of phylodynamic methods in public health, demonstrating their broad application from local to global scales and underscoring the critical role of pathogen genomic data in enhancing public health surveillance and outbreak characterization.

Occurrence of blaOXA-116 Carbapenemase in Escherichia coli ST2519 of Clinical Origin: A Report from Northeast India.

Carbapenem-resistant Escherichia coli pose a significant threat to global public health due to the dearth of available treatment options, resulting in infections with high mortality and morbidity. The study aimed to investigate the mechanism of carbapenem resistance in a carbapenem non-susceptible E. coli isolate recovered from an urinary tract infection patient admitted to a tertiary referral hospital, through whole-genome sequencing using Illumina NovaSeq 6000 platform. Carbapenemase production followed by antibiotic susceptibility testing were performed following Clinical Laboratory Standard Institute guidelines.Polymerase chain reaction targeting carbapenemase genes was performed followed by an investigation of horizontal transferability. The Center for Genomic Epidemiology database was used to analyze the sequenced data. ST2519 E. coli BJD_EC1808 with a genome size of 5.8 Mb harbored Col440I plasmid and a chromosomally located blaOXA-116 gene with an IS18 element upstream, along with multiple antibiotic resistance genes conferring clinical resistance toward beta-lactams, aminoglycosides, amphenicols, sulfonamides, tetracyclines, trimethoprim, rifampin, macrolide, and streptogramin antibiotics and antiseptics. E. coli ST2519 harboring blaOXA-116 associated with a mobile genetic element exhibiting carbapenem resistance is a public health threat due to its limiting effect on the therapeutic usage of carbapenem and their dissemination into carbapenem non-susceptible phenotypes will contribute to carbapenem resistance burden and, therefore, warrants urgent monitoring and clinical intervention.