Whole-genome Sequencing Reads Research Articles

Abstract B059: Assessment of ctDNA analysis methods using high purity ctDNA from the aqueous humor of retinoblastoma eyes

Abstract Background: Whole-genome sequencing (WGS) of cell-free DNA (cfDNA) paired with methods like ichorCNA and DELFI has become a valuable approach for detecting circulating tumor DNA (ctDNA). However, the mixture of ctDNA with non-tumor cfDNA in plasma complicates accurate validation of ctDNA detection methods. In contrast, aqueous humor (AH) from treatment-free retinoblastoma (RB) eyes is uniquely enriched in ctDNA, providing an optimal source for generating reference samples. This study leverages AH-derived biological ctDNA to evaluate two established ctDNA analysis methods. Methods: Biologically derived cfDNA mixtures were created by combining WGS reads from two AH samples with pooled reads from non-tumor blood samples (N=20). One AH sample (RB_023_OS) showed >99% tumor fraction (TFx) with ichorCNA, while the other (RB_025_OD) had confirmed presence of ctDNA based on >99% variant allele frequency (VAF) for an RB1 variant, despite being somatic copy number alteration (SCNA)-negative. Mixtures were prepared across a range of TFx values (0.01-0.5) and unique sequencing coverages (0.1x, 0.5x, 1x, and 2x). TFx was calculated using ichorCNA with standard parameters and optimized parameters for low TFx levels. We also tested an enrichment strategy targeting shorter ctDNA fragments (≤150 base pairs) and generated genome-wide fragmentation profiles using the DELFI method to assess short-to-long (S:L) fragment ratios across 5-Mb bins. Results: We evaluated the impact of sequencing coverage, parameter selection, and fragment size selection on TFx estimation by ichorCNA. Standard ichorCNA parameters demonstrated limited sensitivity, detecting ctDNA only in samples with TFx of 5-6% or higher. Increasing sequencing coverage did not substantially improve the limit of detection (LoD) or the accuracy of TFx quantification. Optimized parameters reduced the LoD to 3% TFx, but only in samples with ≥0.5x coverage, indicating a dependence on sequencing depth for accurate ctDNA detection and quantification. In RB_025_OD, where no SCNAs were present, ichorCNA failed to estimate the TFx despite an RB1 VAF exceeding 99%, highlighting the platform's reliance on SCNAs for accurate TFx detection. Short fragment enrichment (≤150 base pairs) significantly increased TFx values across all samples but could be mapped back to the original TFx values accurately using polynomial curve fitting (R-squared=0.97). Genome- wide fragmentation profiles showed increasing deviation from non-tumor controls as TFx increased, with SCNA-positive samples exhibiting greater deviations than SCNA-negative samples. Conclusions: These findings demonstrate the dependence on sequencing coverage and the presence of SCNAs by current ctDNA analysis methods, highlighting their limitations for accurately detecting ctDNA in challenging scenarios. Citation Format: Nikki Higa, Peter Kuhn, Jesse L Berry, Liya Xu. Assessment of ctDNA analysis methods using high purity ctDNA from the aqueous humor of retinoblastoma eyes [abstract]. In: Proceedings of the AACR Special Conference: Liquid Biopsy: From Discovery to Clinical Implementation; 2024 Nov 13-16; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(21_Suppl):Abstract nr B059.

Whole-genome sequencing resolves biochemical misidentification of Neisseria species from urogenital specimens.

Neisseria meningitidis (Nm) and Neisseria gonorrhoeae (Ng) are human pathogens that sometimes occupy the same anatomical niche. Ng, the causative agent of gonorrhea, infects 87 million individuals annually worldwide and is an urgent threat due to increasing drug resistance. Ng is a pathogen of the urogenital tract and may infect the oropharyngeal or rectal site, often asymptomatically. Conversely, Nm is an opportunistic pathogen. While often a commensal in the oropharyngeal tract, it is also the leading cause of bacterial meningitis with 1.2 million cases globally, causing significant morbidity and mortality. Horizontal gene transfer (HGT) is likely to occur between Ng and Nm due to their shared anatomical niches and genetic similarity, which poses challenges for accurate detection and treatment. Routine surveillance through the Gonococcal Isolate Surveillance Project and Strengthening the U.S. Response to Resistant Gonorrhea detected six concerning urogenital Neisseria isolates with contradicting species identification in Milwaukee (MIL). While all six isolates were positive for Ng using nucleic acid amplification testing (NAAT) and matrix-assisted laser desorption/ionization time of flight identified the isolates as Ng, two biochemical tests, Gonochek-II and API NH, classified them as Nm. To address this discrepancy, we performed whole-genome sequencing (WGS) using Illumina MiSeq on all isolates and employed various bioinformatics tools. Species detection analysis using BMScan, which uses WGS data, identified all isolates as Ng. Furthermore, Kraken revealed over 98% of WGS reads mapped to the Ng genome and <1% to Nm. Recombination analysis identified putative HGT in all MIL isolates within the γ-glutamyl transpeptidase (ggt) gene, a key component in the biochemical tests used to differentiate between Nm and Ng. Further analysis identified Nm as the source of HGT event. Specifically, the active Nm ggt gene replaced the Ng pseudogenes, ggt1 and ggt2. Together, this study demonstrates that closely related Neisseria species sharing a niche underwent HGT, which led to the misidentification of species following biochemical testing. Importantly, NAAT accurately detected Ng. The misidentification highlights the importance of using WGS to continually evaluate diagnostic or bacterial identification tests.

Metagenomic Study of the MESA: Detection of Gemella Morbillorum and Association With Coronary Heart Disease.

Inflammation is a feature of coronary heart disease (CHD), but the role of proinflammatory microbial infection in CHD remains understudied. CHD was defined in the MESA (Multi-Ethnic Study of Atherosclerosis) as myocardial infarction (251 participants), resuscitated arrest (2 participants), and CHD death (80 participants). We analyzed sequencing reads from 4421 MESA participants in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program using the PathSeq workflow of the Genome Analysis Tool Kit and a 65-gigabase microbial reference. Paired reads aligning to 840 microbes were detected in >1% of participants. The association of the presence of microbe reads with incident CHD (follow-up, ~18 years) was examined. First, important variables were ascertained using a single regularized Cox proportional hazard model, examining change of risk as a function of presence of microbe with age, sex, education level, Life's Simple 7, and inflammation. For variables of importance, the hazard ratio (HR) was estimated in separate (unregularized) Cox proportional hazard models including the same covariates (significance threshold Bonferroni corrected P<6×10-5, 0.05/840). Reads from 2 microbes were significantly associated with CHD: Gemella morbillorum (HR, 3.14 [95% CI, 1.92-5.12]; P=4.86×10-6) and Pseudomonas species NFACC19-2 (HR, 3.22 [95% CI, 2.03-5.41]; P=1.58×10-6). Metagenomics of whole-genome sequence reads opens a possible frontier for detection of pathogens for chronic diseases. The association of G morbillorum and Pseudomonas species reads with CHD raises the possibilities that microbes may drive atherosclerotic inflammation and that treatments for specific pathogens may provide clinical utility for CHD reduction.

VIM-1-producing Enterobacter asburiae with mobile colistin resistance genes from wastewaters

BackgroundWastewaters are considered as important players in the spread of antimicrobial resistance, thus affecting the health of humans and animals. Here, we focused on wastewaters as a possible source of carbapenemase-producing Enterobacterales for the environment.MethodsA total of 180 presumptive coliforms from hospital and municipal wastewaters, and a river in the Czech Republic were obtained by selective cultivation on meropenem-supplemented media and tested for presence of carbapenemase-encoding genes by PCR. Strains carrying genes of interest were characterized by testing antimicrobial susceptibility, carbapenemase production and combination of short- and long- read whole-genome sequencing. The phylogenetic tree including publicly available genomes of Enterobacter asburiae was conducted using Prokka, Roary and RAxML.ResultsThree VIM-producing Enterobacter asburiae isolates, members of the Enterobacter cloacae complex, were detected from hospital and municipal wastewaters, and the river. The blaVIM−1 gene was located within a class 1 integron that was carried by different F-type plasmids and one non-typeable plasmid. Furthermore, one of the isolates carried plasmid-borne colistin-resistance gene mcr-10, while in another isolate chromosomally located mcr-9 without colistin resistance phenotype was detected. In addition, the analysis of 685 publicly available E. asburiae genomes showed they frequently carry carbapenemase genes, highlighting the importance of this species in the emergence of resistance to last-line antibiotics.ConclusionOur findings pointed out the important contribution of hospital and community wastewaters in transmission of multi-drug resistant pathogens.

Easy analysis of bacterial whole-genome sequencing data for clinical microbiologists using open-source Galaxy platform: Characterization of ESBL-producing Enterobacterales from bloodstream infections

ObjectivesClinical microbiologists require easy-to-use open access tools with graphical interfaces to perform bacterial whole-genome sequencing (WGS) in routine practice. This study aimed to build a bioinformatics pipeline on the open-source Galaxy platform, facilitating comprehensive and reproducible analysis of bacterial WGS data in a few steps. We then used it to characterize our local epidemiology of ESBL-producing Enterobacterales isolated from patients with bacteremia. MethodsWe built a bioinformatics pipeline consisting of the following sequential tools: Fastp (input data trimming); FastQC (read quality control); SPAdes (genome assembly); Quast (quality control of genome assembly); Prokka (gene annotation); Staramr (ResFinder database) and ABRicate (CARD database) for antimicrobial resistance (AMR) gene screening and molecular strain typing. Paired-end short read WGS data from all ESBL-producing Enterobacterales strains isolated from patients with bacteremia over one year were analysed. ResultsThe Galaxy platform does not require command line tools. The bioinformatics pipeline was constructed within one hour. It only required uploading fastq files and facilitated systematization of de novo assembly of genomes, MLST typing, and AMR gene screening in one step. Among the 66 ESBL-producing strains analysed, the two most frequent ESBL genes were blaCTX−M-15 (62.1%) and blaCTX−M-27 (13.6%). ConclusionsThe open-access Galaxy platform provides a graphical interface and easy-to-use tools suitable for routine use in clinical microbiology laboratories without bioinformatics specialists. We believe that this platform will facilitate fast and low-cost analysis of bacterial WGS data, especially in resource-limited settings.

Metagenome comparison (MC): A new framework for detecting unique/enriched OMUs (operational metagenomic units) derived from whole-genome sequencing reads

BackgroundCurrent methods for comparing metagenomes, derived from whole-genome sequencing reads, include top-down metrics or parametric models such as metagenome-diversity, and bottom-up, non-parametric, model-free machine learning approaches like Naïve Bayes for k-mer-profiling. However, both types are limited in their ability to effectively and comprehensively identify and catalogue unique or enriched metagenomic genes, a critical task in comparative metagenomics. This challenge is significant and complex due to its NP-hard nature, which means computational time grows exponentially, or even faster, with the problem size, rendering it impractical for even the fastest supercomputers without heuristic approximation algorithms. MethodIn this study, we introduce a new framework, MC (Metagenome-Comparison), designed to exhaustively detect and catalogue unique or enriched metagenomic genes (MGs) and their derivatives, including metagenome functional gene clusters (MFGC), or more generally, the operational metagenomic unit (OMU) that can be considered the counterpart of the OTU (operational taxonomic unit) from amplicon sequencing reads. The MC is essentially a heuristic search algorithm guided by pairs of new metrics (termed MG-specificity or OMU-specificity, MG-specificity diversity or OMU-specificity diversity). It is further constrained by statistical significance (P-value) implemented as a pair of statistical tests. ResultsWe evaluated the MC using large metagenomic datasets related to obesity, diabetes, and IBD, and found that the proportions of unique and enriched metagenomic genes ranged from 0.001% to 0.08 % and 0.08%–0.82 % respectively, and less than 10 % for the MFGC. ConclusionThe MC provides a robust method for comparing metagenomes at various scales, from baseline MGs to various function/pathway clusters of metagenomes, collectively termed OMUs.

Identification and mapping of sex-linked markers in giant abalone (Haliotis gigantea Gmelin)

Giant abalone (Haliotis gigantea Gmelin) is an economically important species within Haliotis genus, yet the genetic basis of sex in this species is still poorly understood. Therefore, we screened and mapped sex-linked SNPs and InDels, and subsequently developed a PCR-based method for genetic sex identification in giant abalone. Based on whole-genome sequencing reads, we obtained a total 6,159,941 high-quality SNPs and 2,623,263 high-quality InDels. Among them, 1,737 SNPs and 179 InDels were detected as sexually dimorphic, assuming a male heterogametic sex determination system. The majority of these variants, 99.1 % for SNPs and 99.4 % for InDels, were clustered on Chr.4 (0.02 Mb to 23.22 Mb). Furthermore, the estimation of genetic divergence between the sexes using FST showed that the genomic windows with the top 1 % FST values were concentrated on Chr.4 (0.5 Mb to 19.0 Mb) and Chr.5 (20.0 Mb to 35.0 Mb). Based on the sexually dimorphic variants on Chr.4, three PCR-based markers were successfully developed to determine the genetic sex of giant abalone. PCR validation using a tested panel of 67 females and 67 males from a cultivated population demonstrated compliance rates for these three markers ranging from 97.7 % to 99.3 %. These results suggest that the sex determination system is male heterogametic in giant abalone. Thus, this study provides valuable insights and tools for further research on the mechanisms of sex determination and differentiation, as well as for monitoring sex ratios in wild and cultivated populations in giant abalone.

Reproducibility of next-generation-sequencing-based analysis of a CRISPR/Cas9 genome edited oil seed rape

Next-generation-sequencing (NGS) becomes increasingly important for laboratories tasked with the detection of genetically modified organisms (GMOs) in food, feed and seeds. Its implementation into standardized workflows demands reliable intra- and inter-laboratory reproducibility. Here, we analyze the reproducibility of short- and long-read targeted NGS and long-read whole genome sequencing (WGS) data between three independent laboratories. Replicate samples were submitted for sequencing and comparatively analyzed. The targeted-NGS-samples consisted of oil seed rape (OSR) sampled from a commodity shipment spiked with a genome edited (GE) OSR and the WGS-samples consisted of leaf material from the GMOs’ parental line. All laboratories delivered highly reproducible high-quality targeted NGS data with little variation. The detection of GMO-related sequences works well regardless of the facility, while the mapping to the complex genome is superior using long read data. Long read WGS is currently not suitable for routine use in enforcement laboratories, due to a large inter-laboratory variation.

A comparative investigation of single nucleotide variant calling for a personal non-Caucasian sequencing sample.

Dropping cost and increasing clinical application of whole genome sequencing (WGS) lead a necessity of efficient (accurate and rapid) variant calling procedures from a personal WGS data (n = 1). A number of variant calling pipelines have been introduced utilizing the human genome reference GRCh38 as a reference and a benchmark dataset called 'NA12878', which are both 'standard' but limited ethnic origin. Considering the nature of variant calling algorithms and recent updates in sequencing protocol, however, it is necessary to revisit the efficiency of the current best pipelines for a personal WGS data from diverse ethnicity. We discuss the most efficient practices for variant calling of a personal WGS reads, with a particular emphasis on whether (1) ethnic match or mismatch between the reference genome and a WGS data produces a distinct result and more importantly (2) there is an ethnic-specific optimal workflow. Here, we generate an appropriate WGS data, DNA array, and sufficient number of Sanger validated variants from a single Korean subject to perform such a comprehensive comparison. We applied this WGS reads and the 'NA12878' reads to 8 different variant calling pipelines with 2 different reference genomes (GRCh38 and KOREF, a Korean reference genome) to which the WGS reads from different ethnic origins are aligned. We evaluated the performance of the pipelines with the matched array genotype data and Sanger sequencing validation and demonstrated that: regardless to the ethnic match/mismatch (1) Novoalign-GATK4 showed the most efficient performance with the exceptional calls in MHC region; (2) the overall performance was better with GRCh38, while a significant difference in recall was observed. In addition, we found it is largely reduced computing cost maintaining performance to remove 'markduplication' step with PCR-free WGS data. For variant calling of a personal PCR-free WGS data, regardless of ethnicity consideration, we recommend the use of the Novoalign + GATK4 with GRCh38 and without 'markduplication'.

The Good, The Bad, and The Misleading: How to Improve the Quality of 'Genome Announcements'?

When comparing the requirements of diverse journals to publish microbial 'Genome Reports' we noticed that some mostly focus on Benchmarking Universal Single-Copy Orthologs (BUSCO) scores as a quality measure while the exclusion of possible contaminating sequences from genomic resources and the possible mis-identification of the target microbes receive less attention. To deal with these quality issues, we suggest that DNA barcodes that are widely accepted for the identification of the target microbe species should be extracted from newly reported genome resources and included in phylogenetic analyses to confirm the identity of the sequenced microorganisms before 'Genome Reports' are published. This approach, applied, for example, by the journal IMA Fungus, largely prevents the mis-identification of the microbes that are targeted for whole genome sequencing (WGS). In addition, contig similarity values, including GC-content, remapping coverage of WGS reads, and BLASTN searches against the NCBI nt database, will also reveal contamination issues. The values of these two recommendations to improve the publication criteria for microbial 'Genome Reports' in diverse journals are demonstrated here through analyses of a draft genome published in Molecular Plant-Microbe Interactions, then retracted due to contaminations.

Analysis of coverage of Alu repeats by aligned genomic reads

Alu repeats occupy a notable part of the human genome and greatly affect processes related to genome integrity maintenance. One of the basic methods for studying variation in a genome, including Alu repeats is genome sequencing followed by mapping the sequenced reads to a reference genome sequence. The key feature of the read alignment is the depth of reference genome region coverage by mapped reads. In this paper, a new method is proposed for analyzing the coverage of Alu repeats and their flanking regions by whole-genome sequencing reads and the distribution of mean coverage in two aforementioned region types is explored.

Adequacy of using a single nasal swab for rapid influenza diagnostic testing, PCR, and whole genome sequencing.

Rapid influenza diagnostic tests (RIDT) demonstrate varying sensitivities, often necessitating reverse transcriptase polymerase chain reaction (RT-PCR) to confirm results. The two methods generally require separate specimens. Using the same anterior nasal swab for both RIDT and molecular confirmation would reduce cost and waste and increase patient comfort. The aim of this study was to determine if RIDT residual nasal swab (rNS) specimens are adequate for RT-PCR and whole genome sequencing (WGS). We performed RT-PCR and WGS on paired rNS and nasopharyngeal or oropharyngeal (NP/OP) swab specimens that were collected from primary care patients across all ages. We randomly selected 199 and 40 paired specimens for RT-PCR and WGS, respectively, from the 962 paired surveillance specimens collected during the 2014-2015 influenza season. Sensitivity and specificity for rNS specimens were 81.3% and 96.7%, respectively, as compared to NP/OP specimens. The mean cycle threshold (Ct) value for the NP/OP specimen was significantly lower when the paired specimens were both positive than when the NP/OP swab was positive and the nasal swab was negative (25.5 vs 29.5; p<0.001). Genomic information was extracted from all 40 rNS specimens and 37 of the 40 NP/OP specimens. Complete WGS reads were available for 67.5% (14 influenza A; 13 influenza B) of the rNS specimens and 59.5% (14 influenza A; 8 influenza B) of the NP/OP specimens. It is feasible to use a single anterior nasal swab for RIDT followed by RT-PCR and/or WGS. This approach may be appropriate in situations where training and supplies are limited. Additional studies are needed to determine if residual nasal swabs from other rapid diagnostic tests produce similar results.

Transmission Dynamics of a Mycobacterium tuberculosis Complex Outbreak in an Indigenous Population in the Colombian Amazon Region.

Whole genome sequencing (WGS) has become the main tool for studying the transmission of Mycobacterium tuberculosis complex (MTBC) strains; however, the clonal expansion of one strain often limits its application in local MTBC outbreaks. The use of an alternative reference genome and the inclusion of repetitive regions in the analysis could potentially increase the resolution, but the added value has not yet been defined. Here, we leveraged short and long WGS read data of a previously reported MTBC outbreak in the Colombian Amazon Region to analyze possible transmission chains among 74 patients in the indigenous setting of Puerto Nariño (March to October 2016). In total, 90.5% (67/74) of the patients were infected with one distinct MTBC strain belonging to lineage 4.3.3. Employing a reference genome from an outbreak strain and highly confident single nucleotide polymorphisms (SNPs) in repetitive genomic regions, e.g., the proline-glutamic acid/proline-proline-glutamic-acid (PE/PPE) gene family, increased the phylogenetic resolution compared to a classical H37Rv reference mapping approach. Specifically, the number of differentiating SNPs increased from 890 to 1,094, which resulted in a more granular transmission network as judged by an increasing number of individual nodes in a maximum parsimony tree, i.e., 5 versus 9 nodes. We also found in 29.9% (20/67) of the outbreak isolates, heterogenous alleles at phylogenetically informative sites, suggesting that these patients are infected with more than one clone. In conclusion, customized SNP calling thresholds and employment of a local reference genome for a mapping approach can improve the phylogenetic resolution in highly clonal MTBC populations and help elucidate within-host MTBC diversity. IMPORTANCE The Colombian Amazon around Puerto Nariño has a high tuberculosis burden with a prevalence of 1,267/100,000 people in 2016. Recently, an outbreak of Mycobacterium tuberculosis complex (MTBC) bacteria among the indigenous populations was identified with classical MTBC genotyping methods. Here, we employed a whole-genome sequencing-based outbreak investigation in order to improve the phylogenetic resolution and gain new insights into the transmission dynamics in this remote Colombian Amazon Region. The inclusion of well-supported single nucleotide polymorphisms in repetitive regions and a de novo-assembled local reference genome provided a more granular picture of the circulating outbreak strain and revealed new transmission chains. Multiple patients from different settlements were possibly infected with at least two different clones in this high-incidence setting. Thus, our results have the potential to improve molecular surveillance studies in other high-burden settings, especially regions with few clonal multidrug-resistant (MDR) MTBC lineages/clades.

ONT long-read WGS for variant discovery and orthogonal confirmation of short read WGS derived genetic variants in clinical genetic testing.

Technological advances in Next-Generation Sequencing dramatically increased clinical efficiency of genetic testing, allowing detection of a wide variety of variants, from single nucleotide events to large structural aberrations. Whole Genome Sequencing (WGS) has allowed exploration of areas of the genome that might not have been targeted by other approaches, such as intergenic regions. A single technique detecting all genetic variants at once is intended to expedite the diagnostic process while making it more comprehensive and efficient. Nevertheless, there are still several shortcomings that cannot be effectively addressed by short read sequencing, such as determination of the precise size of short tandem repeat (STR) expansions, phasing of potentially compound recessive variants, resolution of some structural variants and exact determination of their boundaries, etc. Therefore, in some cases variants can only be tentatively detected by short reads sequencing and require orthogonal confirmation, particularly for clinical reporting purposes. Moreover, certain regulatory authorities, for example, New York state CLIA, require orthogonal confirmation of every reportable variant. Such orthogonal confirmations often involve numerous different techniques, not necessarily available in the same laboratory and not always performed in an expedited manner, thus negating the advantages of "one-technique-for-all" approach, and making the process lengthy, prone to logistical and analytical faults, and financially inefficient. Fortunately, those weak spots of short read sequencing can be compensated by long read technology that have comparable or better detection of some types of variants while lacking the mentioned above limitations of short read sequencing. At Variantyx we have developed an integrated clinical genetic testing approach, augmenting short read WGS-based variant detection with Oxford Nanopore Technologies (ONT) long read sequencing, providing simultaneous orthogonal confirmation of all types of variants with the additional benefit of improved identification of exact size and position of the detected aberrations. The validation study of this augmented test has demonstrated that Oxford Nanopore Technologies sequencing can efficiently verify multiple types of reportable variants, thus ensuring highly reliable detection and a quick turnaround time for WGS-based clinical genetic testing.

Molecular Epidemiology of Global Carbapenemase-Producing Citrobacter spp. (2015-2017).

The emergence of carbapenem resistance is a significant public health concern. The rate of infections caused by carbapenemase-producing Citrobacter spp., particularly C. freundii, is increasing. Concomitantly, comprehensive global genomic data on carbapenemase-producing Citrobacter spp. are scarce. We used short read whole-genome sequencing to describe the molecular epidemiology and international distribution of eighty-six carbapenemase-producing Citrobacter spp. obtained from two surveillance programs (2015 to 17). The common carbapenemases were KPC-2 (26%), VIM-1 (17%), IMP-4 (14%) and NDM-1 (10%). C. freundii and C. portucalensis were the principal species. C. freundii consisted of multiple clones obtained mainly from Colombia (with KPC-2), the United States (with KPC-2, -3), and Italy (with VIM-1). Two dominant C. freundii clones were identified: ST98 was linked with blaIMP-8 from Taiwan and blaKPC-2 from the United States, and ST22 was linked with blaKPC-2 from Colombia and blaVIM-1 from Italy. C. portucalensis consisted mainly of two clones: ST493 with blaIMP-4 which was limited to Australia, and ST545 with blaVIM-31 which was limited to Turkey. Class I integron (In916) with blaVIM-1 was circulating between multiple sequence types (STs) in Italy, Poland, and Portugal. In73 with blaIMP-8 was circulating between various STs in Taiwan, while In809 with blaIMP-4 was circulating between different STs in Australia. The global carbapenemase-producing Citrobacter spp. population is dominated by diverse STs with different characteristics and varied geographical distribution and thus requires continued monitoring. Ongoing genomic surveillance should use methodologies able to distinguish between C. freundii and C. portucalensis. IMPORTANCE Citrobacter spp. are gaining recognition as important causes of hospital-acquired infections in humans. Among Citrobacter spp., carbapenemase-producing strains are cause of utmost concern to health care services globally due to their ability to resist therapy with virtually any beta-lactam antibiotic. Here, we described the molecular characteristics of a global collection of carbapenemase-producing Citrobacter spp. C. freundii and C. portucalensis were the most common species among Citrobacter spp. with carbapenemases from this survey. Importantly, C. portucalensis was misidentified as C. freundii when using Vitek 2.0/MALDI-TOF MS (matrix-assisted laser desorption/ionization-time of flight mass spectrometry) phenotypic identification, which has important implications for future surveys. Among C. freundii, we identified two dominant clones: ST98 with blaIMP-8 from Taiwan and blaKPC-2 from the United States, and ST22 with blaKPC-2 from Colombia and blaVIM-1 from Italy. As for C. portucalensis, the dominant clones consisted of ST493 with blaIMP-4 from Australia and ST545 with blaVIM-31 from Turkey.

GAMBIT (Genomic Approximation Method for Bacterial Identification and Tracking): A methodology to rapidly leverage whole genome sequencing of bacterial isolates for clinical identification.

Whole genome sequencing (WGS) of clinical bacterial isolates has the potential to transform the fields of diagnostics and public health. To realize this potential, bioinformatic software that reports identification results needs to be developed that meets the quality standards of a diagnostic test. We developed GAMBIT (Genomic Approximation Method for Bacterial Identification and Tracking) using k-mer based strategies for identification of bacteria based on WGS reads. GAMBIT incorporates this algorithm with a highly curated searchable database of 48,224 genomes. Herein, we describe validation of the scoring methodology, parameter robustness, establishment of confidence thresholds and the curation of the reference database. We assessed GAMBIT by way of validation studies when it was deployed as a laboratory-developed test in two public health laboratories. This method greatly reduces or eliminates false identifications which are often detrimental in a clinical setting.

A reference genome for Bluegill (Centrarchidae: Lepomis macrochirus).

North American sunfishes (Family Centrarchidae) are among the most popular sportfish throughout the United States and Canada. Despite the popularity of sunfishes, their ecological importance, and their extensive stocking and aquacultural history, few molecular studies have examined the evolutionary relationships and species boundaries among members of this group, many of which are known to hybridize. Here, we describe a chromosome-scale genome assembly representing Bluegill (Lepomis macrochirus), one of the most widespread centrarchid species. By combining long-read, Oxford Nanopore sequencing data with short-insert, whole-genome and HiC sequence reads, we produced an assembly (Lm_LA_1.1) having a total length of 889 Mb including 1,841 scaffolds and having a scaffold N50 of 36 Mb, L50 of 12, N90 of 29 Mb, and L90 of 22. We detected 99% (eukaryota_odb10) and 98% (actinopterygii_odb10) universal single-copy orthologs (BUSCOs), and ab initio gene prediction performed using this new assembly identified a set of 17,233 genes that were supported by external (OrthoDB v10) data. This new assembly provides an important addition to the growing set of assemblies already available for spiny-rayed fishes (Acanthomorpha), and it will serve as a resource for future studies that focus on the complex evolutionary history of centrarchids.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation.

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

Diagnostic Utility of Whole Genome Sequencing in Adults with B-Other Acute Lymphoblastic Leukaemia

Background B-other ALL represents a working definition for patients with B cell precursor (BCP) ALL without a known primary chromosomal abnormality at diagnosis. In this study we use whole genome sequencing (WGS) to characterize the molecular pathogenesis of adult B-other cases (age ≥ 25yrs) from the UKALL14 trial (NCT01085617). Methods Of 652 patients aged 25-65yrs enrolled onto UKALL14, 333 (51%) had B-other ALL. Sufficient material was available to screen 156/333 B-other cases for recurrent Ph-like fusion events (CRLF2, JAK2, ABL1, ABL2 & PDGFRB) using FISH and MLPA (kit P335). This identified 28 (18%) "Ph-like" fusion events (21 CRLF2, 5 ABL-class and 2 JAK2). Of the remaining 128 B-other cases 57 had available samples for tumor normal paired WGS (read depth 60x and 30x respectively). Diagnostic RNA was available for 33 of the 57 cases Bioinformatic analysis was performed to determine small somatic mutations; single nucleotide variants (SNVs) and insertion/deletions (INDELs) as well as copy number aberrations (CNA) and structural variants (SV). We developed WGS approaches to identify IGH enhancer hijack and DUX4 rearrangements, utilizing gGnome and GRIDSS respectively. For 33 cases RNA-seq data was analyzed to validate fusion genes and genetic classification into one of 18 ALL subtypes using the ALLSorts v2 classifier. Results 52/57 cases passed QC. In 5 (10%) cases an established WHO2016 genetic subtype was identified that was missed by standard of care including three high-risk (2 low hypodiploid and a near haploid) and 2 standard risk (high hyperdiploid and TCF3::PBX1) cases. Cytogenetic analysis for these 5 cases either failed (n=2) or produced a normal karyotype (n=3). The WGS inferred tumor purity for these 5 cases was variable with a median of 46% and a range of 33% to 96% compared to 89% (range 13%-99%) for the remaining cases (n=47) in the dataset. Among the remaining 47 B-other cases, WGS analysis identified abnormalities enabling 41 (87%) cases to be classified into a genetic ALL subtype: DUX4 rearrangements (DUX4r, n=8), PAX5alt (n=7), PAX5 P80R (n=4), ZNF384 rearrangements (ZNF384r, n=5), ZEB2/CEBP (n=4), MEF2D rearrangements (MEF2Dr, n=3), UBTF::ATXN7L3 fusions (n=2), IDH1/2 mutations (n=2), JAK-STAT abnormalities (n=2), along with single cases of IGK::BCL2, IGH::CEBPA, IGH::ID4 and IGH::MIR125B1. WGS resolves 8 cases with complex karyotype into one of 6 different genomic subtypes: 2xDUX4r, 2xMEF2Dr, ZEB2/CEBP, UBTF::ATXN7L3, IGK::BCL2 & IGH::ID4. The MEF2Dr subtype (n=3), is enriched in chromothripsis (2/3), and complex SVs targeting the MEF2D gene. Lastly, in all cases with PAX5 mutations (PAX5 P80R (n=4) and PAX5alt (n=4)) we observe bi-allelic targeting of locus by mutations, deletions or LOH. The remaining 6 (13%) B-other cases could not be classified into a recognized B-other subtype. A MYO18A::FGFR1 fusion (t(8;17)(p11;q23)), which most likely maps to the WHO2016 "Lymphoid neoplasms with FGFR1 rearrangement" classification. Two cases had clonal frameshifts in genes associated with T-ALL (PTPN2 p.A108fs*5 and WT1 p.V371fs*14 respectively). A case with a somatic ETV6 p.R399C previously reported in germline risk with a ZNF292::KDM6A fusion. Lastly a case with a clonal mutation typically associated with myeloid disease ASXL1 p.R417*. Among the 52 cases studied, WGS only failed to identify a driver or putative driver in one case, which had the lowest burden of SNVs (861) and an estimated tumor purity, 28%. For the 31 cases with both DNA and RNA, WGS and RNA analysis were concordant for 100% of cases, with RNA identifying an additional KDM6A::CGA fusion. In 3 cases WGS provides evidence that would have challenged diagnosis and/or treatment; IGK::BCL2 with BCL2 and MYC overexpression, IGH::DUX4 with hypermutation (MSH6 loss) and IGH::ID4 with homozygous loss of CD58 and concurrent LOH of HLA-B. Conclusions WGS assigned 88% (46/52) of cases called B-other to an established genetic subtype of ALL. This includes 5 cases with WHO2016 defined subtypes, 3 were high-risk ploidy subtypes that would have changed their UKALL14 risk group and post-induction treatment in the absence of other risk factors. Among the true B-other cases, 87% (41/47) were assigned to one of the newly described genetic subtypes of ALL. Taken together our findings demonstrate the clinical utility of WGS as a diagnostic assay to inform and improve the management of adult B-Other ALL patients in the future.

Gamete Binning to Achieve Haplotype-Resolved Genome Assembly.

Haplotype-resolved genome assemblies remain a challenge in practice. Here, we provide a step-by-step guide on gamete binning, a method to generate haplotype-resolved genome assemblies for diploid species. The protocol starts by phasing heterozygous variants to individual haplotypes of specific chromosomes using the genome information of individual haploid gametes of the focal individual. Using phased variants, the whole-genome sequencing reads from the diploid genome can be genotyped and assigned into groups, which represent the individual haplotypes of each of the chromosomes. Finally, haplotype-specific chromosomes can be assembled independently using standard assembly tools. First applications of gamete binning revealed a haplotyping accuracy over 99%, which outperformed sequence-only or Hi-C-based haplotype-resolved genome assemblies.Availability: github.com/schneebergerlab/GameteBinning_prac .