WGS Data Research Articles

GeneMates: an R package for detecting horizontal gene co-transfer between bacteria using gene-gene associations controlled for population structure

BackgroundHorizontal gene transfer contributes to bacterial evolution through mobilising genes across various taxonomical boundaries. It is frequently mediated by mobile genetic elements (MGEs), which may capture, maintain, and rearrange mobile genes and co-mobilise them between bacteria, causing horizontal gene co-transfer (HGcoT). This physical linkage between mobile genes poses a great threat to public health as it facilitates dissemination and co-selection of clinically important genes amongst bacteria. Although rapid accumulation of bacterial whole-genome sequencing data since the 2000s enables study of HGcoT at the population level, results based on genetic co-occurrence counts and simple association tests are usually confounded by bacterial population structure when sampled bacteria belong to the same species, leading to spurious conclusions.ResultsWe have developed a network approach to explore WGS data for evidence of intraspecies HGcoT and have implemented it in R package GeneMates (github.com/wanyuac/GeneMates). The package takes as input an allelic presence-absence matrix of interested genes and a matrix of core-genome single-nucleotide polymorphisms, performs association tests with linear mixed models controlled for population structure, produces a network of significantly associated alleles, and identifies clusters within the network as plausible co-transferred alleles. GeneMates users may choose to score consistency of allelic physical distances measured in genome assemblies using a novel approach we have developed and overlay scores to the network for further evidence of HGcoT. Validation studies of GeneMates on known acquired antimicrobial resistance genes in Escherichia coli and Salmonella Typhimurium show advantages of our network approach over simple association analysis: (1) distinguishing between allelic co-occurrence driven by HGcoT and that driven by clonal reproduction, (2) evaluating effects of population structure on allelic co-occurrence, and (3) direct links between allele clusters in the network and MGEs when physical distances are incorporated.ConclusionGeneMates offers an effective approach to detection of intraspecies HGcoT using WGS data.

Listeria monocytogenes clones in Italian food products: virulence and environmental adaptation

Abstract Background Food is the main source of Listeria monocytogenes (Lm) infection. Lm is a highly heterogeneous species composed of hypervirulent and hypovirulent clones. Understanding the distribution of Lm clonal complexes (CCs) in different food categories has strong implications for risk assessment. The aim of this work was to analyse collection of Lm strains of National Reference Laboratory (NRL Lm) in order to assess link between genetic profile and matrices and the level of pathogenicity of circulating strains based on CCs. Methods NRL Lm database actually consists of 906 sequenced strains isolated in 10 years from 5 food compartments (meat, fish, dairy, vegetables and composite dishes). Epidata were analysed to remove redundant strains based on the same epidemiological description. After that, WGS data from 465 Lm strains were investigated. In silico MLST was defined and Roary 3.12.0 was used to obtain a pan-genome profile. Genes were later uploaded to Pasteur Institute platform for characterization. Results In silico MLST identified 36 CCs and 6 singleton. CC9 (23.0%), CC8 (15.3%) and CC121 (13.3%) were the prevalent CCs. In particular, CC9 was present in 35.2% of meat samples and CC8 in 25.8% of fish samples. Pan genome profile revealed high prevalence (>98%) of genes related to biofilm formation and resistance to environmental stress in CC9 strains and genes involved in tolerance to quaternary ammonium compounds in CC121 strains. Conclusions Results, in particular for meat products, confirmed in Italy, the prevalence of hypovirulent Lm strains previously observed at European Union level. The high presence of stress resistance and disinfectant tolerance genes in these strains could make them able to persist in food-production environment and should be taken into account evaluating the health hazards. In fish product is also relevant the prevalence of CC8 strains which are potentially highly pathogenic and have been responsible of recent European multi country outbreak. Key messages Pangenome of Listeria monocytogenes isolated from Italian food revealed high presence of disinfectant tolerance and stress resistance genes in meat products and virulence genes in fish products. Listeria monocytogenes CC9 and CC121 prevalence in Italian meat product confirms occurrence of hypovirulent strains detected at European Union level.

Multiplex real-time PCR assays for the prediction of cephalosporin, ciprofloxacin and azithromycin antimicrobial susceptibility of positive Neisseria gonorrhoeae nucleic acid amplification test samples.

The incidence of antimicrobial-resistant Neisseria gonorrhoeae (GC) is rising in Canada; however, antimicrobial resistance (AMR) surveillance data are unavailable for infections diagnosed directly from clinical specimens by nucleic acid amplification tests (NAATs), representing over 80% of diagnoses. We developed a set of 10 improved molecular assays for surveillance of GC-AMR and prediction of susceptibilities in NAAT specimens. Multiplex real-time PCR (RT-PCR) assays were developed to detect SNPs associated with cephalosporin (ponA, porB, mtrR -35delA, penA A311V, penA A501, N513Y, G545S), ciprofloxacin (gyrA S91, parC D86/S87/S88) and azithromycin [23S (A2059G, C2611T), mtrR meningitidis-like promoter] resistance. The assays were validated on 127 gonococcal isolates, 51 non-gonococcal isolates and 50 NAATs with matched culture isolates. SNPs determined from the assay were compared with SNPs determined from in silico analysis of WGS data. MICs were determined for culture isolates using the agar dilution method. SNP analysis of the 50 NAAT specimens had 96% agreement with the matched culture RT-PCR analysis. When compared with MICs, presence of penA A311V or penA A501 and two or more other SNPs correlated with decreased susceptibility and presence of three or more other SNPs correlated with intermediate susceptibility to cephalosporins; presence of any associated SNP correlated with ciprofloxacin or azithromycin resistance. NAAT-AMR predictions correlated with matched-culture cephalosporin, ciprofloxacin and azithromycin MICs at 94%, 100% and 98%, respectively. We expanded molecular tests for N. gonorrhoeae AMR prediction by adding new loci and multiplexing reactions to improve surveillance where culture isolates are unavailable.

A method for scoring the cell type-specific impacts of noncoding variants in personal genomes

A person's genome typically contains millions of variants which represent the differences between this personal genome and the reference human genome. The interpretation of these variants, i.e., the assessment of their potential impact on a person's phenotype, is currently of great interest in human genetics and medicine. We have developed a prioritization tool called OpenCausal which takes as inputs 1) a personal genome and 2) a reference context-specific TF expression profile and returns a list of noncoding variants prioritized according to their impact on chromatin accessibility for any given genomic region of interest. We applied OpenCausal to 6,430 samples across 18 tissues derived from the GTEx project and found that the variants prioritized by OpenCausal are highly enriched for eQTLs and caQTLs. We further propose a strategy to integrate the predicted open scores with genome-wide association studies (GWAS) data to prioritize putative causal variants and regulatory elements for a given risk locus (i.e., fine-mapping analysis). As an initial example, we applied this method to a GWAS dataset of human height and found that the prioritized putative variants and elements are correlated with the phenotype (i.e., heights of individuals) better than others.

Abstract 1067: Proteogenomic discovery of novel tumor proteins as neoantigens for personalized T cell immunotherapy in pediatric medulloblastoma

Abstract Medulloblastoma is a common childhood brain tumor for which current therapies fail to cure high-risk subgroups, including metastatic and recurrent disease. Toxicity limits further intensification of conventional chemo-radiation therapy and most survivors endure significant lifelong sequelae from treatment. There is an urgent need to develop new tools to combat this cancer without increasing the late effects burden. Neoantigens exclusively expressed on tumor cells, are one of the main targets of an effective antitumor T-cell response. The ideal target antigen is abundantly expressed by tumor cells but not by normal tissues, in order to limit off-target effects. Tumors translate a host of unique transcripts that are the result of aberrations in tumor DNA and the unmasking of alternative or novel exons. We used a novel proteogenomic approach to identify tumor-restricted peptides and used them to select and expand T cells capable of mounting a tumor-specific cytotoxic immune response, with no off-target effects for patients with MB. Using RNA-seq and WGS data, we created personalized custom searchable databases containing predicted novel proteins from somatic mutations, novel isoforms and fusion proteins from 46 medulloblastoma tumors. By searching this database with raw mass spectrometry data from the paired medulloblastoma tumor, we have identified tens of peptides arising from the translation of tumor-specific transcripts; these peptides are potential neoantigens. Our data indicates that in cases of tumors with low mutation rates, such as pediatric brain tumors, novel isoforms are the main source of neoantigens. In a pilot study, we tested these peptides for their ability to select and expand polyclonal populations of T cells from a patient whose tumor was the source of the peptides. The immunogenicity of each individual peptide was then determined. Flow cytometry cellular characterization reveals populations of both CD4+ and CD8+ cells with an activation profile marked by IFN-γ and TNF-α production. Inhibitory co-receptor profiles were also characterized for these cells. Using cytotoxicity assays, we demonstrated that tumor specific T cells can eliminate neoantigen bearing tumor cells. These findings demonstrate an initial proof of principle that proteogenomics can be used to identify immunogenic tumor specific peptides and lay the groundwork for a personalized, targeted T cell therapy for children with high risk medulloblastoma. Citation Format: Samuel Rivero-Hinojosa, Melanie Grant, Aswini Panigrahi, Huizhen Zhang, Veronika Caisova, Catherine M. Bollard, Brian R. Rood. Proteogenomic discovery of novel tumor proteins as neoantigens for personalized T cell immunotherapy in pediatric medulloblastoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1067.

Abstract 6333: Genomic, transcriptomic, and epigenetic profiling of triple-negative breast cancer cells after Navitoclax treatment

Abstract Cancer cells possess various anti-apoptotic defense mechanisms that can protect them from drug induced cell death. We investigated the effects of treatment with the bcl-2/bcl-XL inhibitor Navitoclax on transcriptome, methylome, and genome of MDA-MB-231 cells. Two biological replicates were treated simultaneously and profiled using 10x Genomics single-cell RNA sequencing (scRNAseq; gene expression), bulk RNA sequencing (bRNAseq; gene expression), ATAC sequencing (ATAC; chromatin accessibility), bisulphite targeted sequencing (DNAm; DNA methylation) and shallow whole-genome sequencing (WGS; copy number variants) at baseline before treatment, after 3 days of 10µM navitoclax treatment, and after 10 days of cell recovery from treatment. Variance between biological replicates on transcriptome level was minimal before the treatment (R>0.99), but we noticed divergence in gene expression after recovery from the treatment (both on bulk and single-cell level), mostly driven by cell cycle genes. On other platforms, biological replicates were similar in all 3 phases. After treatment, we observed more genes with decreased expression (n=151) in comparison to baseline, but after 10 days of recovery there were more up-regulated genes (n=655) in comparison to samples after treatment. Chromatin accessibility for most genomic regions increased after treatment in comparison to baseline, suggesting an acute response to the stress caused by drug treatment, and then returned to baseline level after recovery period. DNA methylation patterns showed more regions with decreased methylation after treatment (n=529) in comparison to baseline which remained detectable even after the recovery period. In WGS data, we found 752 genes with deletions and amplifications only in baseline samples representing genomic background of drug sensitive cells, mostly enriched in epithelial mesenchymal transition pathway. On single-cell level, we identified a subset of cells that were resistant to treatment and discovered 2,324 genes significantly up-regulated in these cells, that could be potential markers of resistance. In gene set analysis, these markers were enriched in MYC and E2F target gene sets, and were involved in angiogenesis and JAK-STAT pathway. We measured the expression of 16 top up-regulated markers of resistance in 4 different TNBC cell lines using qPCR and found that 5 were significantly enriched after treatment in 3 cell lines, 5 in 2 cell lines and 5 in a single cell line. Summarizing, we thoroughly described molecular effects of Navitoclax treatment and showed that most cells return to the basal transcriptional state after the drug recovery period, but bulk genome and methylome are permanently changed. Finally, we provided a list of new markers of resistance that may be useful in the studies of combinational therapies with other drugs. This work was funded by Breast Cancer Research Foundation. Citation Format: Michal Marczyk, Vignesh Gunasekharan, Jun Zhao, Rihao Qu, Xiaotong Li, Gauri A. Patwardhan, Vikram B. Wali, Abhishek K. Gupta, Manoj M. Pillai, Yuval Kluger, Christos Hatzis, Lajos Pusztai. Genomic, transcriptomic, and epigenetic profiling of triple-negative breast cancer cells after Navitoclax treatment [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 6333.

Strain-Level Metagenomic Data Analysis of Enriched In Vitro and In Silico Spiked Food Samples: Paving the Way towards a Culture-Free Foodborne Outbreak Investigation Using STEC as a Case Study.

Culture-independent diagnostics, such as metagenomic shotgun sequencing of food samples, could not only reduce the turnaround time of samples in an outbreak investigation, but also allow the detection of multi-species and multi-strain outbreaks. For successful foodborne outbreak investigation using a metagenomic approach, it is, however, necessary to bioinformatically separate the genomes of individual strains, including strains belonging to the same species, present in a microbial community, which has up until now not been demonstrated for this application. The current work shows the feasibility of strain-level metagenomics of enriched food matrix samples making use of data analysis tools that classify reads against a sequence database. It includes a brief comparison of two database-based read classification tools, Sigma and Sparse, using a mock community obtained by in vitro spiking minced meat with a Shiga toxin-producing Escherichia coli (STEC) isolate originating from a described outbreak. The more optimal tool Sigma was further evaluated using in silico simulated metagenomic data to explore the possibilities and limitations of this data analysis approach. The performed analysis allowed us to link the pathogenic strains from food samples to human isolates previously collected during the same outbreak, demonstrating that the metagenomic approach could be applied for the rapid source tracking of foodborne outbreaks. To our knowledge, this is the first study demonstrating a data analysis approach for detailed characterization and phylogenetic placement of multiple bacterial strains of one species from shotgun metagenomic WGS data of an enriched food sample.

Benchmarking variant callers in next-generation and third-generation sequencing analysis.

DNA variants represent an important source of genetic variations among individuals. Next- generation sequencing (NGS) is the most popular technology for genome-wide variant calling. Third-generation sequencing (TGS) has also recently been used in genetic studies. Although many variant callers are available, no single caller can call both types of variants on NGS or TGS data with high sensitivity and specificity. In this study, we systematically evaluated 11 variant callers on 12 NGS and TGS datasets. For germline variant calling, we tested DNAseq and DNAscope modes from Sentieon, HaplotypeCaller mode from GATK and WGS mode from DeepVariant. All the four callers had comparable performance on NGS data and 30× coverage of WGS data was recommended. For germline variant calling on TGS data, we tested DNAseq mode from Sentieon, HaplotypeCaller mode from GATK and PACBIO mode from DeepVariant. All the three callers had similar performance in SNP calling, while DeepVariant outperformed the others in InDel calling. TGS detected more variants than NGS, particularly in complex and repetitive regions. For somatic variant calling on NGS, we tested TNscope and TNseq modes from Sentieon, MuTect2 mode from GATK, NeuSomatic, VarScan2, and Strelka2. TNscope and Mutect2 outperformed the other callers. A higher proportion of tumor sample purity (from 10 to 20%) significantly increased the recall value of calling. Finally, computational costs of the callers were compared and Sentieon required the least computational cost. These results suggest that careful selection of a tool and parameters is needed for accurate SNP or InDel calling under different scenarios.

A 63-bp insertion in exon 2 of the porcine KIF21A gene is associated with arthrogryposis multiplex congenita.

A recessive form of arthrogryposis multiplex congenita (AMC) was detected 20years ago in the Swiss Large White (SLW) pig population. A diagnostic marker test enabled the identification of carrier animals, but the underlying causal mutation remains unknown. To identify the mutation underlying AMC, we collected SNP chip genotyping data for 11 affected piglets and 23 healthy pigs. Association testing using 47829 SNPs confirmed that AMC maps to SSC5 (P=9.4×10-13 ). Subsequent autozygosity mapping revealed a common 6.06Mb region (from 66 757970 to 72815151bp) of extended homozygosity in 11 piglets affected by AMC. Using WGS data, we detected a 63-bp insertion compatible with the recessive inheritance of AMC in the second exon of KIF21A gene encoding Kinesin Family Member 21A. The 63-bp insertion is predicted to introduce a premature stop codon in KIF21A gene (p.Val41_Phe42insTer) that truncates 1614 amino acids (~97%) from the protein. We found that this deleterious allele still segregates at a frequency of 0.1% in the SLW pig population. Carrier animals can now be detected unambiguously and excluded from breeding.

Core Genome Multilocus Sequence Typing for Food Animal Source Attribution of Human Campylobacter jejuni Infections.

Campylobacter jejuni is a major foodborne pathogen and common cause of bacterial enteritis worldwide. A total of 622 C. jejuni isolates recovered from food animals and retail meats in the United States through the National Antimicrobial Resistance Monitoring System between 2013 and 2017 were sequenced using an Illumina MiSeq. Sequences were combined with WGS data of 222 human isolates downloaded from NCBI and analyzed by core genome multilocus sequence typing (cgMLST) and traditional MLST. cgMLST allelic difference (AD) thresholds of 0, 5, 10, 25, 50, 100 and 200 identified 828, 734, 652, 543, 422, 298 and 197 cgMLST types among the 844 isolates, respectively, and traditional MLST identified 174 ST. The cgMLST scheme allowing an AD of 200 (cgMLST200) revealed strong correlation with MLST. cgMLST200 showed 40.5% retail chicken isolates, 56.5% swine, 77.4% dairy cattle and 78.9% beef cattle isolates shared cgMLST sequence type with human isolates. All ST-8 had the same cgMLST200 type (cgMLST200-12) and 74.3% of ST-8 and 75% cgMLST200-12 were confirmed as sheep abortion virulence clones by PorA analysis. Twenty-nine acquired resistance genes, including 21 alleles of blaOXA, tetO, aph(3′)-IIIa, ant(6)-Ia, aadE, aad9, aph(2′)-Ig, aph(2′)-Ih, sat4 plus mutations in gyrA, 23SrRNA and L22 were identified. Resistance genotypes were strongly linked with cgMLST200 type for certain groups including 12/12 cgMLST200-510 with the A103V substitution in L22 and 10/11 cgMLST200-608 with the T86I GyrA substitution associated with macrolide and quinolone resistance, respectively. In summary, the cgMLST200 threshold scheme combined with resistance genotype information could provide an excellent subtyping scheme for source attribution of human C. jejuni infections.

1748-P: One Novel 7.2kb Deletion in Exon 8 of INSR Gene in Diabetic Qatari Female with Type A Insulin Resistance Syndrome: The 1000 Qatar-Omics Study Cohort

Background: A 46-year-old Qatari female with a BMI of 25kg/m2 presented with hyperglycemia (blood glucose 13.5 mmol/l), hyperinsulinemia (21µU/mL) with elevated levels of HbA1c (10.3mmol/l). Biochemical investigations showed raised serum C-peptide (2.02ng/ml) and triglycerides (2.07mmol/L) and uric acid (104mg/dL), with normal BP. The objective of this case study is to validate this variant found through WGS in insulin receptor (INSR) gene that play a major cause of insulin resistant syndrome with different spectrums. Methods: Structural Variants discovered using custom multi-algorithm pipeline developed and optimized to run on HPC for short-read WGS data. In brief, 10 software packages applied (Breakdancer, Breakseq2, CNVnator, Delly, ERDS, Genomestrip, Manta, Speedseq, Svaba, WHAM), to existing short-read data. SVs annotated using AnnotSV and provides us functionally, regulatory and clinically relevant information. SVs analyzed based on their properties like pli-score, overlaps with known databases and OMIM annotations. This deletion was validated by ddPCR and the molecular docking performed by the autodock suit in order to understand the effect of deletion in the interaction between this receptor and insulin. Results: The patient was identified with 7.2k heterozygous deletion in exon 8 of INSR gene which play an essential role in the formation of extracellular alpha subunits in insulin receptor. molecular docking proved the decrease of catalytic activity of INSR and, hence signaling. Considering this and the phenotypic data from Qatar biobank, the patient developed features of type A insulin resistant syndrome as consequence of persistence hyperinsulinemia. The patient also exhibits severe Deterioration of glucose tolerant even with medication and unfortunately developed nephropathy. Conclusion: We highlighted the power of NGS in precise diagnosis and developing a potential and promising targeted therapeutics. Disclosure A.S. Akil: None. E.A.M. Yassin: None. S. Subash Padmajeya: None. L.A. Jerman: None. A. Fadda: None. W.H.O. Aamer: None. E.E. Aliyev: None. K. Fakhro: None.

Challenges in interpretation of WGS and epidemiological data to investigate nosocomial transmission of vancomycin-resistant Enterococcus faecium in an endemic region: incorporation of patient movement network and admission screening.

ObjectivesVRE are listed, by the WHO, among the leading resistant pathogens causing greatest public concern; hence the spread and transmission of VRE, especially in hospitalized patients, need to be monitored. Despite the advancements in typing methods since the implementation of WGS for outbreak investigations, data interpretation, especially for vancomycin-resistant Enterococcus faecium (VREfm) in an endemic setting, remains challenging. In this study we explored the potential added benefit of incorporating patient movement data and admission screening to accurately estimate the magnitude of an outbreak.MethodsWe sequenced 73 VREfm isolates from patients with bacteraemia (n = 43) and rectal colonization (n = 30/32). Genetic relatedness was determined by SNP distance (≤10) between isolates. Patient movements were visualized in a movement network, along with contact intensity and rectal colonization status prior to infection onset.ResultsST117, ST80 and ST203 were the predominant STs in our study population. Forty-four percent (18/41) of VREfm bacteraemia cases were of endogenous origin. SNP analysis of infection and colonization isolates revealed nine clonal groups. Eighty-six percent (37/43) of the patients were visualized in a transmission network due to spatiotemporal overlap. Nineteen out of 43 (44%) belonged to five transmission clusters. Incorporation of prior colonization status revealed that transmission was very likely in only 63% (12/19) of patients in these transmission clusters.DiscussionAlthough interpretation of WGS data is challenging, incorporation of patient movement data and colonization status by admission screening of high-risk patients may provide additional resolution when interpreting the magnitude of an outbreak in an endemic setting.

Sudden death in epilepsy and ectopic neurohypophysis in Joubert syndrome 23 diagnosed using SNVs/indels and structural variants pipelines on WGS data: a case report

BackgroundJoubert syndrome (JBTS) is a genetically heterogeneous group of neurodevelopmental syndromes caused by primary cilia dysfunction. Usually the neurological presentation starts with abnormal neonatal breathing followed by muscular hypotonia, psychomotor delay, and cerebellar ataxia. Cerebral MRI shows mid- and hindbrain anomalies including the molar tooth sign. We report a male patient with atypical presentation of Joubert syndrome type 23, thus expanding the phenotype.Case presentationClinical features were consistent with JBTS already from infancy, yet the syndrome was not suspected before cerebral MRI later in childhood showed the characteristic molar tooth sign and ectopic neurohypophysis. From age 11 years seizures developed and after few years became increasingly difficult to treat, also related to inadequate compliance to therapy. He died at 23 years of sudden unexpected death in epilepsy (SUDEP). The genetic diagnosis remained elusive for many years, despite extensive genetic testing. We reached the genetic diagnosis by performing whole genome sequencing of the family trio and analyzing the data with the combination of one analysis pipeline for single nucleotide variants (SNVs)/indels and one for structural variants (SVs). This lead to the identification of the most common variant detected in patients with JBTS23 (OMIM# 616490), rs534542684, in compound heterozygosity with a 8.3 kb deletion in KIAA0586, not previously reported.ConclusionsWe describe for the first time ectopic neurohypophysis and SUDEP in JBTS23, expanding the phenotype of this condition and raising the attention on the possible severity of the epilepsy in this disease. We also highlight the diagnostic power of WGS, which efficiently detects SNVs/indels and in addition allows the identification of SVs.

Global clonal spread of mcr-3-carrying MDR ST34 Salmonella enterica serotype Typhimurium and monophasic 1,4,[5],12:i:- variants from clinical isolates.

To investigate the prevalence and transmission of mcr-3 among Salmonella enterica serotype Typhimurium and 1,4,[5],12:i:-. A total of 4724 clinical Salmonella isolates were screened for the presence of mcr-3 in China during 2014-19. The clonal relationship of the mcr-3-positive isolates and their plasmid contents and complete sequence were also characterized based on WGS data from the Illumina and MinION platforms. We identified 10 mcr-3-positive isolates, and all were MDR, mostly resistant to colistin, cefotaxime, ciprofloxacin, doxycycline and florfenicol. mcr-3 was co-present with blaCTX-M-55-qnrS1 on hybrid ST3-IncC-FII conjugatable plasmids (n = 6) and an ST3-IncC non-conjugatable plasmid (n = 1) and embedded into a pCHL5009T-like IncFII plasmid on the Salmonella chromosome (n = 3). Four distinctive genetic contexts surrounded mcr-3 and all but one were closely related to each other and to the corresponding region of IncFII plasmid pCHL5009T. IS15DI was most likely the vehicle for integration of mcr-3-carrying IncFII plasmids into ST3-IncC plasmids and the chromosome and for shaping the MDR regions. In addition, a phylogenetic tree based on the core genome revealed a unique Salmonella lineage (≤665 SNPs) that contained these 10 mcr-3-positive isolates and another 38 (33 from patients) mcr-3-positive Salmonella from five countries. In particular, most of the 51 mcr-3-positive isolates belonged to ST34 and harboured diverse antibiotic resistance genes (ARGs), including mcr-3-blaCTX-M-55-qnrS1, and possessed similar ARG profiles. Our findings revealed global clonal spread of MDR ST34 Salmonella from clinical isolates co-harbouring mcr-3 with blaCTX-M-55 and qnrS1 and a flexibility of mcr-3 co-transmittance with other ARGs mediated by mobile genetic elements.

CNV-BAC: Copy number Variation Detection in Bacterial Circular Genome.

Whole-genome sequencing (WGS) is widely used for copy number variation (CNV) detection. However, for most bacteria, their circular genome structure and high replication rate make reads more enriched near the replication origin. CNV detection based on read depth could be seriously influenced by such replication bias. We show that the replication bias is widespread using ∼200 bacterial WGS data. We develop CNV-BAC (CNV-Bacteria) that can properly normalize the replication bias and other known biases in bacterial WGS data and can accurately detect CNVs. Simulation and real data analysis show that CNV-BAC achieves the best performance in CNV detection compared with available algorithms. CNV-BAC is available at https://github.com/XiDsLab/CNV-BAC. Supplementary data are available at Bioinformatics online.

Integration of the blaNDM-1 carbapenemase gene into a novel SXT/R391 integrative and conjugative element in Proteus vulgaris.

To characterize the genetic environment of the carbapenem resistance determinant in Proteus vulgaris of swine origin. The carbapenem-resistant P. vulgaris strain BC22 was isolated from a faecal swab from a diseased pig with diarrhoea in Sichuan Province of China in 2018. The presence of carbapenemase genes was screened by PCR. WGS and bioinformatics analysis were performed to analyse the genetic environment of the carbapenem resistance determinant. P. vulgaris strain BC22 was found to harbour the carbapenemase gene blaNDM-1. WGS data revealed that blaNDM-1 was located in a truncated ISAba125 composite transposon. The carbapenem resistance gene blaNDM-1 and 20 other resistance genes, including the multiresistance gene cfr and the bifunctional aminoglycoside/quinolone resistance gene aac(6')-lb-cr, were located in a novel SXT/R391 integrative and conjugative element (ICE). This new SXT/R391 ICE of 148.7 kb was chromosomally located, and could be transferred to Escherichia coli. Here, we report a carbapenemase gene, blaNDM-1, integrated into an SXT/R391 ICE. Our study highlights that this SXT/R391 ICE may facilitate the dissemination of clinically important resistance genes such as blaNDM-1, cfr and aac(6')-lb-cr.

Prolonged carriage of (livestock-associated) MRSA in individuals without professional livestock contact.

To investigate prolonged carriage of MRSA in adults from the general population living in a livestock-dense area, using WGS. A cross-sectional study during 2014-15 among 2492 adults without professional livestock contact identified 14 (0.6%) nasal MRSA carriers, 10 of which carried livestock-associated (LA)-MRSA of multiple-locus variable-number tandem repeat analysis (MLVA) complex (MC) 398. Two years later, 12 MRSA-positive and 88 MRSA-negative participants provided a second nasal swab and filled in a short questionnaire. Isolates from persons who were MRSA positive at both timepoints were compared using MLVA and isolates with the same MLVA type were sequenced. The WGS data were used for core-genome MLST (cgMLST) and resistome analysis, including sequenced isolates from the national MRSA surveillance. All MRSA-negative persons tested negative again, while 6 of the 12 initially MRSA-positive persons tested positive again. MLVA revealed that isolate pairs from five individuals had the same MLVA type, of which three were LA-MRSA. cgMLST showed that the distance between these isolate pairs ranged between 3 and 13 genes, while the minimum distance to unrelated isolates from the national MRSA surveillance was 38 genes. Moreover, the resistome present in the five isolate pairs was identical within each pair. None of the prolonged carriers was hospitalized during the 3 months before the sampling moment and none of them with LA-MRSA had contact with livestock in this period. Prolonged carriage of MRSA, including LA-MRSA, can be demonstrated after more than 30 months in persons without professional livestock contact.

Evolution of Listeria monocytogenes in a Food Processing Plant Involves Limited Single-Nucleotide Substitutions but Considerable Diversification by Gain and Loss of Prophages

Whole-genome sequencing (WGS) is becoming the standard method for subtyping Listeria monocytogenes Interpretation of WGS data for isolates from foods and associated environments is, however, challenging due to a lack of detailed data on Listeria evolution in processing facilities. Here, we used previously collected WGS data for 40 L. monocytogenes isolates obtained from a cold-smoked salmon processing facility between 1998 and 2015 to probe the L. monocytogenes molecular evolution in this facility, combined with phenotypic assessment of selected isolates. Isolates represented three clusters (1, 2, and 3); cluster 3 isolates (n = 32) were obtained over 18 years. The average mutation rate for cluster 3 was estimated as 1.15 × 10-7 changes per nucleotide per year (∼0.35 changes per genome per year); the most recent common ancestors (MRCAs) of subclusters 3a and 3b were estimated to have occurred around 1958 and 1974, respectively, within the age of the facility, suggesting long-term persistence in this facility. Extensive prophage diversity was observed within subclusters 3a and 3b, which have one shared and six unique prophage profiles for each subcluster (with 16 prophage profiles found among all 40 isolates). The plasmid-borne sanitizer tolerance operon bcrABC was found in all cluster 2 and 3 isolates, while the transposon-borne sanitizer tolerance gene qacH was found in one cluster 1 isolate; presence of these genes was correlated with the ability to survive increased concentrations of sanitizers. Selected isolates showed significant variation in the ability to attach to surfaces, with persistent isolates attaching better than transient isolates at 21°C.IMPORTANCE Knowledge about the genetic evolution of L. monocytogenes in food processing facilities over multiple years is generally lacking. This information is critical to interpret WGS findings involving food or food-associated isolates. This study suggests that L. monocytogenes that persists in processing facilities may evolve with a low single-nucleotide mutation rate mostly driven by negative (i.e., purifying) selection but with rapid diversification of prophages. Hence, isolation of L. monocytogenes with few single-nucleotide polymorphism (SNP) differences in different locations (e.g., supplier plants and receiving plants) is possible, highlighting the importance of epidemiological and detailed isolate metadata for interpreting WGS data in traceback investigation. Our study also shows how advanced WGS data analyses can be used to support root cause analysis efforts and may, for example, pinpoint the time when a persistence event started (which then potentially could be linked to facility changes, introduction of new equipment, etc.).

Abstract A23: Proteogenomic discovery of novel tumor proteins as neoantigens for personalized T-cell immunotherapy in pediatric medulloblastoma

Abstract Medulloblastoma is a common childhood brain tumor for which current therapies fail to cure high-risk subgroups, including metastatic and recurrent disease. Toxicity limits further intensification of conventional chemo-radiation therapy, and most survivors endure significant lifelong sequelae from treatment. There is an urgent need to develop new tools to combat this cancer without increasing the late effects burden. Neoantigens, which are exclusively expressed on tumor cells, are one of the main targets of an effective antitumor T-cell response. The ideal target antigen is abundantly expressed by tumor cells but not by normal tissues, in order to limit off-target effects. Tumors translate a host of highly novel transcripts that are the result of aberrations in tumor DNA and the unmasking of alternative or novel exons. We used a novel proteogenomic approach to identify tumor-restricted peptides and used them to select and expand T cells capable of mounting a tumor-specific cytotoxic immune response, with no off-target effects for patients with MB. Using RNA-seq and WGS data, we created personalized custom searchable databases containing predicted novel proteins from somatic mutations, novel isoforms, and fusion proteins from 56 medulloblastoma tumors. By searching this database with raw mass spectrometry data from the paired medulloblastoma tumor, we have identified tens of peptides arising from the translation of tumor-specific transcripts; these peptides are potential neoantigens. Our data indicate that in cases of tumors with low mutation rates, such as pediatric brain tumors, novel isoforms are the main source of neoantigens. In a pilot study, we have tested these peptides for their ability to select and expand polyclonal populations of T cells from the same patient whose tumor was the source of the peptides. The immunogenicity of each individual peptide was then determined. Flow cytometry cellular characterization reveals populations of both CD4+ and CD8+ cells with an activation profile marked by IFN-γ and TNF-α production. Immunosuppressive marker profiles have also been characterized for these cells. Using cytotoxicity assays, we demonstrated that tumor-specific T cells can eliminate neoantigen-bearing tumor cells. These findings demonstrate an initial proof of principle that proteogenomics can be used to identify immunogenic tumor specific peptides and lay the groundwork for a personalized, targeted T-cell therapy for children with high-risk medulloblastoma. Citation Format: Samuel Rivero-Hinojosa, Melanie Grant, Aswini Panigrahi, Huizhen Zhang, Veronika Caisova, Catherine Bollard, Brian Rood. Proteogenomic discovery of novel tumor proteins as neoantigens for personalized T-cell immunotherapy in pediatric medulloblastoma [abstract]. In: Proceedings of the AACR Special Conference on Tumor Immunology and Immunotherapy; 2019 Nov 17-20; Boston, MA. Philadelphia (PA): AACR; Cancer Immunol Res 2020;8(3 Suppl):Abstract nr A23.

Detection of Molecular Signatures of Homologous Recombination Deficiency in Prostate Cancer with or without BRCA1/2 Mutations

Prostate cancers with mutations in genes involved in homologous recombination (HR), most commonly BRCA2, respond favorably to PARP inhibition and platinum-based chemotherapy. We investigated whether other prostate tumors that do not harbor deleterious mutations in these particular genes can similarly be deficient in HR, likely rendering those sensitive to HR-directed therapies. Homologous recombination deficiency (HRD) levels can be estimated using various mutational signatures derived from next-generation sequencing data. We used this approach on whole-genome sequencing (WGS; n = 311) and whole-exome sequencing (WES) data (n = 498) of both primary and metastatic prostate adenocarcinomas to determine whether prostate cancer cases display clear signs of HRD in somatic tumor biopsies. Known BRCA-deficient samples showed all previously described HRD-associated mutational signatures in the WGS data. HRD-associated mutational signatures were also detected in a subset of patients who did not harbor germline or somatic mutations in BRCA1/2 or other HR-related genes. Similar results, albeit with lower sensitivity and accuracy, were also obtained from WES data. These findings may expand the number of cases likely to respond to PARP inhibitor treatment. On the basis of the HR-associated mutational signatures, 5% to 8% of localized prostate cancer cases may be good candidates for PARP-inhibitor treatment (including those with BRCA1/2 mutations).