Genomic Fragments Research Articles

Abstract 3376: Novel reference materials for the analysis of methylation in liquid biopsies

Abstract Here we describe new reference materials (RMs) for liquid biopsies that are designed to mimic the methylation found in circulating cell-free DNA (ccfDNA). Methylation is an important epigenetic modification that influences cellular differentiation and gene expression. Recently, liquid biopsies have started to incorporate analyses of epigenetic modifications for screening purposes to detect cancer-derived DNA in the blood. Additionally, epigenetic modifications are being used to also assign a tissue of origin to the cancer to direct confirmatory diagnostic procedures. Obtaining sufficient amounts of ccfDNA for assay development, validation, and proficiency testing is difficult. Furthermore, methods of analyzing the epigenetic modifications, such as bisulfite conversion, can damage a significant fraction of the input material; but, failing to carry them out to completion can result in an over-estimation of methylation. To address the challenges associated with assessing the methylation of ccfDNA, we created and characterized three types of RMs. The first type of RM is a fragmented genome of the GM24385 reference cell line that is used to create a set of RMs with different degrees of methylation (0%, ~25%, ~50%, ~75% and close to 100%). The exact degree of methylation was assigned by digital PCR using assays for heterozygous SNPs, where one allele is both a site for CpG methylation and a recognition site for a methylation-sensitive restriction enzyme while the other allele is not. The second type of RM is an amplified ccfDNA sample that has been methylated similarly but also reproduces the fragment lengths and genomic representation biases of ccfDNA. The degree of methylation was assigned by digital PCR, but only a subset of assays was useful because not all SNPs were heterozygous. The third type of RM is an amplified ccfDNA sample that has been processed to mimic the methylation found in the original ccfDNA. Analyses of the RMs by whole genome sequencing (WGS) showed that, while bisulfite conversion and/or the amplification of the resulting dU-containing DNA appeared to lead to significant biases in genomic representation, the resulting data was generally in agreement with the amount of methylation assigned by digital PCR. Additionally, it was possible to prepare an RM that generally reproduces the methylation found in ccfDNA; although, this was limited to individual molecules having either no or close to full methylation, which is unlikely to be the case for all molecules of ccfDNA in vivo. Citation Format: Jayanthi Ramprakash, Matthew G. Butler, Russell K. Garlick, Yves Konigshofer. Novel reference materials for the analysis of methylation in liquid biopsies [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3376.

CRISPR/Cas9-based simple transgenesis in Xenopus laevis

Transgenic techniques have greatly increased our understanding of the transcriptional regulation of target genes through live reporter imaging, as well as the spatiotemporal function of a gene using loss- and gain-of-function constructs. In Xenopus species, two well-established transgenic methods, restriction enzyme-mediated integration and I-SceI meganuclease-mediated transgenesis, have been used to generate transgenic animals. However, donor plasmids are randomly integrated into the Xenopus genome in both methods. Here, we established a new and simple targeted transgenesis technique based on CRISPR/Cas9 in Xenopus laevis. In this method, Cas9 ribonucleoprotein (RNP) targeting a putative harbor site (the transforming growth factor beta receptor 2-like (tgfbr2l) locus) and a preset donor plasmid DNA were co-injected into the one-cell stage embryos of X. laevis. Approximately 10% of faithful reporter expression was detected in F0 crispants in a promoter/enhancer-specific manner. Importantly, efficient germline transmission and stable transgene expression were observed in the F1 offspring. The simplicity of this method only required preparation of a donor vector containing the tgfbr2l genome fragment and Cas9 RNP targeting this site, which are common experimental procedures used in Xenopus laboratories. Our improved technique allows the simple generation of transgenic X. laevis, so is expected to become a powerful tool for reporter assay and gene function analysis.

Human parainfluenza 2 & 4: Clinical and genetic epidemiology in the UK, 2013-2017, reveals distinct disease features and co-circulating genomic subtypes.

Human Parainfluenza viruses (HPIV) comprise of four members of the genetically distinct genera of Respirovirus (HPIV1&3) and Orthorubulavirus (HPIV2&4), causing significant upper and lower respiratory tract infections worldwide, particularly in children. However, despite frequent molecular diagnosis, they are frequently considered collectively or with HPIV4 overlooked entirely. We therefore investigated clinical and viral epidemiological distinctions of the relatively less prevalent Orthorubulaviruses HPIV2&4 at a regional UK hospital across four autumn/winter epidemic seasons. A retrospective audit of clinical features of all HPIV2 or HPIV4 RT-PCR-positive patients, diagnosed between 1st September 2013 and 12th April 2017 was undertaken, alongside sequencing of viral genome fragments in a representative subset of samples. Infection was observed across all age groups, but predominantly in children under nine and adults over 40, with almost twice as many HPIV4 as HPIV2 cases. Fever, abnormal haematology, elevated C-reactive protein and hospital admission were more frequently seen in HPIV2 than HPIV4 infection. Each of the four seasonal peaks of either HPIV2, HPIV4 or both, closely matched that of RSV, occurring in November and December and preceding that of Influenza A. A subset of viruses were partially sequenced, indicating co-circulation of multiple subtypes of both HPIV2&4, but with little variation between each epidemic season or from limited global reference sequences. Despite being closest known genetic relatives, our data indicates a potential difference in associated disease between HPIV2 and HPIV4, with more hospitalisation seen in HPIV2 mono-infected individuals, but a greater overall number of HPIV4 cases.

Detection of recombinant breakpoint in the genome of human enterovirus E11 strain associated with a fatal nosocomial outbreak

BackgroundThe aim of this study was to characterize the genome of a recombinant Enterovirus associated with severe and fatal nosocomial infection; it was typed as Echovirus 11 (E-11) according to the VP1 gene. Enterovirus infection is generally asymptomatic and self-limited, but occasionally it may progress to a more severe clinical manifestation, as in the case described here. Recombination plays a crucial role in the evolution of Enteroviruses (EVs) and has been recognized as the main driving force behind the emergence of epidemic strains associated with severe infection. Therefore, it is of utmost importance to monitor the circulation of recombinant strains for surveillance purposes.MethodsEnterovirus-RNA was detected in the serum and liver biopsy of patients involved in the nosocomial cluster by commercial One-Step qRT-PCR method and the Enterovirus strains were isolated in vitro. The EVs typing was determined by analyzing the partial-length of the 5′UTR and VP1 sequences with the web-based open-access Enterovirus Genotyping Tool Version 0.1. The amplicons targeting 5′UTR, VP1 and overlapping fragments of the entire genome were sequenced with the Sanger method. Phylogenetic analysis was performed comparing the VP1 and the full-genome sequences of our strains against an appropriate reference set of Enterovirus prototypes of the Picornaviridae genera and species retrieved from the Enterovirus Genotyping Tool. Recombination analysis was performed using RDP4 software.ResultsThe Neighbor-Joining tree of the VP1 gene revealed that the 4 patients were infected with an identical molecular variant of Echovirus 11 (E-11). While the phylogenetic and the RDP4 analysis of the full-genome sequences provided evidence that it was a chimeric strain between an E-11 and a Coxsackievirus B (CV-B).ConclusionsThe chimeric structure of the E-11 genome might have contributed to the severe infection and epidemic feature of the strain, but further biological characterizations are needed. The evidence reported in this study, highlights the limit of typing techniques based on the VP1 gene, as they fail to identify the emergence of recombinant strains with potentially more pathogenic or epidemic properties, thus providing only partial information on the epidemiology and pathogenesis of Enteroviruses.

PATH-10. Nanopore sequencing reveals novel ALK fusion with interposed element in a neonate with hemispheric glioma

Abstract Here we describe the clinical course and unique molecular findings in a female neonate with infantile hemispheric glioma (IHG). The patient presented at the age of 17 days with macrocephaly and suboptimal weight gain. MRI brain revealed an 8.5cm mass over the right frontal region with significant mass effect and entrapment of ventricles. Urgent ventriculoperitoneal shunt insertion and biopsy of the highly vascular lesion was performed. Pathology was compatible with glioblastoma multiforme. Methylation profiling classified the sample as inflammatory tumor tissue (MNP v11b4), while panel RNA-sequencing revealed a fusion event between HMBOX1 and ALK which has not been described in primary CNS tumors. Long-read sequencing with the Nanopore system further revealed complex genomic rearrangement involving an interposed genomic fragment from a third chromosome with validation by Sanger sequencing. Based on an integrated diagnosis of IHG, the patient went on to receive neoadjuvant chemotherapy per the CNS-14 protocol (cyclophosphamide, carboplatin, etoposide) with significant tumor shrinkage. Near total removal of the lesion was achieved after 4 cycles of chemotherapy and adjuvant treatment with 4 additional cycles was given. The patient tolerated therapy well other than self-limiting transaminitis. At the end of treatment, the patient enjoyed intact neurology and satisfactory developmental progress. Our case report illustrates the value a multi-pronged approach to characterizing rare pediatric CNS tumors. The precise delineation of breakpoints in fusion-driven tumors might be of value in the design of cfDNA-based assays. (APYL and CTLC contributed equally to the submission)

300-OR: A Modular Massively Parallel Reporter Assay Uncovers Context-Specific Activity of Diabetes-Associated Regulatory Elements

Recent genome-wide association studies (GWAS) have identified hundreds of signals associated with type 2 diabetes (T2D) , most (∼90%) of which are found in noncoding regions of the genome enriched in tissue-specific regulatory regions. However, most GWAS signals span hundreds of tightly linked variants, and predicting each of their functional effects is difficult. Consequently, progress to translate associations to mechanisms has proceeded slowly. Here, we have designed a massively parallel reporter assay (MPRA) to systematically screen enhancer activity across a panel of genomic fragments, including (1) islet-specific TSSs defined by cap analysis of gene expression (CAGE) and (2) allelic pairs for over 10k T2D- and metabolic trait-associated variants. We cloned this library up- or downstream of a reporter gene driven by a synthetic housekeeping promoter (SCP1) or the human insulin (INS) promoter, together enabling enhancer activity measurements across different promoter and enhancer contexts. We examined expression bias within our library (FDR < 0.05) and found that 43% of CAGE-defined oligos were marked by positional bias (n = 565/1,305) , skewed towards higher expression in the upstream position (n = 539/565) . In contrast, two unique subsets of GWAS oligos displayed either positional (n = 702/11,656) or promoter (n = 698/11,656) bias. Oligos with positional bias were evenly split between the two cloning positions, while most promoter-biased oligos were more highly expressed when paired with the INS promoter (n = 512/698) . To identify sequence features associated with promoter bias, we performed Lasso regression using a set of 562 genomic annotations. This analysis revealed that presence of HNF1 motifs was associated with higher enhancer activity in the context of the INS promoter. Together, these results indicate that reporter construct design clearly influences assayed regulatory activity, particularly with use of tissue-relevant promoters. Disclosure A.Tovar: None. Funding National Institutes of Health (R01DK117960, T32DK101357)

Comparative genomics of the Western Hemisphere soft tick-borne relapsing fever borreliae highlights extensive plasmid diversity

BackgroundTick-borne relapsing fever (TBRF) is a globally prevalent, yet under-studied vector-borne disease transmitted by soft and hard bodied ticks. While soft TBRF (sTBRF) spirochetes have been described for over a century, our understanding of the molecular mechanisms facilitating vector and host adaptation is poorly understood. This is due to the complexity of their small (~ 1.5 Mb) but fragmented genomes that typically consist of a linear chromosome and both linear and circular plasmids. A majority of sTBRF spirochete genomes’ plasmid sequences are either missing or are deposited as unassembled sequences. Consequently, our goal was to generate complete, plasmid-resolved genomes for a comparative analysis of sTBRF species of the Western Hemisphere.ResultsUtilizing a Borrelia specific pipeline, genomes of sTBRF spirochetes from the Western Hemisphere were sequenced and assembled using a combination of short- and long-read sequencing technologies. Included in the analysis were the two recently isolated species from Central and South America, Borrelia puertoricensis n. sp. and Borrelia venezuelensis, respectively. Plasmid analyses identified diverse sequences that clustered plasmids into 30 families; however, only three families were conserved and syntenic across all species. We also compared two species, B. venezuelensis and Borrelia turicatae, which were isolated ~ 6,800 km apart and from different tick vector species but were previously reported to be genetically similar.ConclusionsTo truly understand the biological differences observed between species of TBRF spirochetes, complete chromosome and plasmid sequences are needed. This comparative genomic analysis highlights high chromosomal synteny across the species yet diverse plasmid composition. This was particularly true for B. turicatae and B. venezuelensis, which had high average nucleotide identity yet extensive plasmid diversity. These findings are foundational for future endeavors to evaluate the role of plasmids in vector and host adaptation.

Orienting Ordered Scaffolds: Complexity and Algorithms

Despite the recent progress in genome sequencing and assembly, many of the currently available assembled genomes come in a draft form. Such draft genomes consist of a large number of genomic fragments (scaffolds), whose order and/or orientation (i.e., strand) in the genome are unknown. There exist various scaffold assembly methods, which attempt to determine the order and orientation of scaffolds along the genome chromosomes. Some of these methods (e.g., based on FISH physical mapping, chromatin conformation capture, etc.) can infer the order of scaffolds, but not necessarily their orientation. This leads to a special case of the scaffold orientation problem (i.e., deducing the orientation of each scaffold) with a known order of the scaffolds. We address the problem of orientating ordered scaffolds as an optimization problem based on given weighted orientations of scaffolds and their pairs (e.g., coming from pair-end sequencing reads, long reads, or homologous relations). We formalize this problem using notion of a scaffold graph (i.e., a graph, where vertices correspond to the assembled contigs or scaffolds and edges represent connections between them). We prove that this problem is NP-hard, and present a polynomial-time algorithm for solving its special case, where orientation of each scaffold is imposed relatively to at most two other scaffolds. We further develop an FPT algorithm for the general case of the OOS problem.

Standards of conduct in residential long-term nursing care in view of the SARS-CoV-2 virus pandemic

Introduction: SARS - CoV - 2 virus, which causes COVID - 19 disease, is transmitted mainly by droplets. The symptoms of infection are: fever, wet cough, sore throat, shortness of breath, headache, pain in muscles and bones, general malaise accompanied by exhaustion. Anosomy (loss of smell) and dysgeusia are also common, as well as vomiting, nausea, diarrhea. The diagnosis of SARS-CoV-2 virus infection is based on molecular tests involving the detection of fragments of the viral genome in real time PCR (real time polymerase reaction) in nasopharyngeal swabs, less often throat swabs.
 Objective of the work: The aim of the study is to summarize the guidelines for long-term nursing care for a patient infected with SARS - CoV - 2 in the patient's home environment.
 Material and methods: The work was based on the method of non-systematic review of scientific literature. The following databases were searched: PubMed, EBESCO, SCOPUS, Web of Science, according to keywords in Polish and English: SARS - CoV - 2, long-term care, nursing, recommendations. The annual limits 2019 - 2020 were assumed as the search period.
 Results and conclusions:The SARS-CoV-2 virus appears in residential long-term care facilities primarily through its transmission from the external environment by the entity's employees. Proper management of long-term care facilities is the basis for reducing COVID - 19 transmission, and thus the number of confirmed cases. The knowledge about the recommendations and recommendations of authorized bodies regarding the proceedings limiting the transmission of the virus should be regularly updated.

Survey of the association between polymorphisms of CTLA-4 exon 1 49 A/G genes with rheumatoid arthritis in Iran

ABSTRACT Cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), which suppresses T cell proliferation, is a promising candidate for the susceptibility genes to rheumatic arthritis diseases (RA). This study aims to examine the association between the polymorphisms of the CTLA-4 exon 1(+ 49) genes with RA in the Qazvin city of Iran population. The polymerase chain reaction of genomic DNA-restriction fragment length polymorphism (PCR-RFLP) was applied to genotype the CTLA-4 exon 1(+ 49) polymorphisms in 105 RA patients and 90 control subjects. Laboratory diagnostic tests were also measured for RA and control groups. Our results did not demonstrate a significant difference in allele and genotype frequencies of the CTLA-4 exon 1(+ 49) between RA patients and the control group (p < .0001). There was no significant difference in age at onset, CRP, RF value in patients with RA according to the CTLA-4 polymorphisms; just anti-CCP showed a significant difference. Our data declared that polymorphisms of CTLA-4 exon 1(+ 49) genes are not correlated with RA susceptibility and its clinical and paraclinical manifestations.

HIV-Quasipore: A Suite of HIV-1-Specific Nanopore Basecallers Designed to Enhance Viral Quasispecies Detection

Accounting for genetic variation is an essential consideration during human immunodeficiency virus type 1 (HIV-1) investigation. Nanopore sequencing preserves proviral integrity by passing long genomic fragments through ionic channels, allowing reads that span the entire genome of different viral quasispecies (vQS). However, this sequencing method has suffered from high error rates, limiting its utility. This was the inspiration behind HIV-Quasipore: an HIV-1-specific Nanopore basecaller suite designed to overcome these error rates through training with gold-standard data. It comprises three deep learning-based R9.4.1 basecallers: fast, high accuracy (HAC), super accuracy (SUP), and two R10.3 deep learning-based basecallers: HAC and SUP. This was accomplished by sequencing the HIV-1 J-Lat 10.6 cell line using Nanopore and high-quality Sanger techniques. Training significantly reduced basecaller error rates across all models (Student’s one-sided t-test; p = 0.0) where median error rates were 0.0189, 0.0018, 0.0008, for R9.4.1 HIV-Quasipore-fast, HAC, SUP, and 0.0007, 0.0011 for R10.3 HIV-Quasipore-HAC, and SUP, respectively. This improved quality reduces the resolution needed to accurately detect a vQS from 22.4 to 2.6% of total positional coverage for R9.4.1 HIV-Quasipore-fast, 6.9 to 0.5% for R9.4.1 HIV-Quasipore-HAC, 4.5 to 0.3% for R9.4.1 HIV-Quasipore-SUP, 8.0 to 0.3% for R10.3 HIV-Quasipore-HAC, and 5.4 to 0.3% for R10.3 HIV-Quasipore-SUP. This was consistently observed across the entire J-Lat 10.6 genome and maintained across longer reads. Reads with greater than 8,000 nucleotides display a median nucleotide identity of 0.9819, 0.9982, and 0.9991, for R9.4.1 HIV-Quasipore-fast, HAC, SUP, and 0.9993, 0.9988 for R10.3 HIV-Quasipore-HAC, and SUP, respectively. To evaluate the robustness of this tool against unseen data, HIV-Quasipore and their corresponding pretrained basecallers were used to sequence the J-Lat 9.2 cell line and a clinical isolate acquired from the Drexel Medicine CARES cohort. When sample reads were compared against their corresponding consensus sequence, all HIV-Quasipore basecallers displayed higher median alignment accuracies than their pretrained counterparts for both the J-Lat 9.2 cell line and clinical isolate. Using Nanopore sequencing can allow investigators to explore topics, such as vQS profile detection, HIV-1 integration site analysis, whole genome amplification, gene coevolution, and CRISPR-induced indel detection, among others. HIV-Quasipore basecallers can be acquired here: https://github.com/DamLabResources/HIV-Quasipore-basecallers.

Novel Multiplex PCR Method and Genome Sequence-Based Analog for High-Resolution Subclonal Assignment and Characterization of Escherichia coli Sequence Type 131 Isolates.

ABSTRACTEscherichia coli sequence type 131 (ST131) is a pandemic, multidrug-resistant extraintestinal pathogen. The multiple distinctive ST131 subclones differ for rfb and fliC alleles (O and H antigens), fimH allele (type-1 fimbriae adhesin), resistance phenotype and genotype, clinical correlates, and host predilection. Current PCR assays for detecting ST131 and its main subclones offer limited sub-ST characterization. Here we combined 22 novel and 14 published primers for a multiplex PCR assay to detect and extensively characterize ST131 isolates. The primers target mdh36, gyrB47, trpA72, sbmA, plsB, nupC, rmuC, kefC, ybbW, the O16 and O25b rfb variants, five fimH alleles (fimH22, fimH27, fimH30, fimH35, and fimH41), two fliC alleles (H4 and H5), a (subclone-specific) fluoroquinolone resistance-associated parC allele, and a (subclone-specific) prophage marker. The resulting amplicons resolve 15 molecular subsets within ST131, including 3 within clade A (H41 subclone), 5 within clade B (H22 subclone), and 7 within clade C (H30 subclone), which includes subclones C0 (H30S: 2 subsets), C1 and C1-M27 (H30R1: 2 subsets), and C2 (H30Rx: 3 subsets). Validation in three laboratories showed that this assay provides a rapid, accurate, and portable method for rapidly detecting and characterizing E. coli ST131 and its key subsets. Additionally, for users with whole genome sequencing (WGS) capability, we developed a command-line executable called ST131Typer, an in silico version of the extended multiplex PCR assay. Its accuracy was 87.8%, with most issues due to incomplete or fragmented input genome assemblies. These two novel assays should facilitate detailed ST131 subtyping using either endpoint PCR or WGS.IMPORTANCE These novel assays provide greater subclonal resolution and characterization of E. coli ST131 isolates than do the available comparable PCR assays, plus offer a novel sequence-based alternative to PCR. They may prove useful for molecular epidemiological studies, surveillance, and, potentially, clinical management.

Phylogenetic analysis of BKV genetic variations, based on the whole sequence of the genome and different genomic sections.

BK polyomavirus (BKV) primarily infects humans in their early life stages, and in later life stages, immunosuppressed patients may develop asymptomatic infections. The nucleotides 1744-1812 in the VP1 gene are traditionally used to determine this virus's genotype. The complete genome of the BKVsamples from patients referred to Masih Daneshvari Hospital's Virology Research Center was amplified by previously known primer sets. The phylogenetic diversity of the whole genome, different genomic sections, and the noncoding control region of BKVsamples were investigated. Using software Mega X and references, the samples' genotype was determined in separate genomic fragments and the whole genome. The samples were classified into two genotypes (I and IV) and five subtypes (Ia, Ib-1, Ib-2, IVc-1, and IVc-2), but none of the isolates belonged to genotypes II, III, V, or VI. The large T antigen-based phylogenetic tree provided 100% bootstrap values for these divisions, which were superior to those (96%-100%) used in the VP1 sequence. Among the genomic segments, large tumor antigen and VP1 had the most mutations. The noncoding control area contained mutations at the O41 position in the granulocyte/macrophage stimulus gene and the P31 position in the NF-1 gene. The validity of the phylogenetic analysis was supported by sequence analysis, which found single-nucleotide polymorphisms(SNPs) that could be useful for subclassifying isolates. More research with a large number of samples and in the wider geographical areas is needed to understand the genetic diversity of the BKV in Iran and also to determine these SNPs' clinical significance in terms of patient outcome and viral load dynamics.

Robust AAV Genotyping Based on Genetic Distances in Rep Gene That Are Maintained by Ubiquitous Recombination.

Adeno-associated viruses (AAVs) are a convenient tool for gene therapy delivery. According to the current classification, they are divided into the species AAV A and AAV B within the genus Dependoparvovirus. Historically AAVs were also subdivided on the intraspecies level into 13 serotypes, which differ in tissue tropism and targeted gene delivery capacity. Serotype, however, is not a universal taxonomic category, and their assignment is not always robust. Cross-reactivity has been shown, indicating that classification could not rely on the results of serological tests alone. Moreover, since the isolation of AAV4, all subsequent AAVs were subdivided into serotypes based primarily on genetic differences and phylogenetic reconstructions. An increased interest in the use of AAV as a gene delivery tool justifies the need to improve the existing classification. Here, we suggest genotype-based AAV classification below the species level based on the rep gene. A robust threshold was established as 10% nt differences within the 1248 nt genome fragment, with 4 distinct AAV genotypes identified. This distinct sub-species structure is maintained by ubiquitous recombination within, but not between, rep genes of the suggested genotypes.

How clear is our current view on microbial dark matter? (Re-)assessing public MAG & SAG datasets with MDMcleaner.

As of today, the majority of environmental microorganisms remain uncultured and is therefore referred to as ‘microbial dark matter’ (MDM). Hence, genomic insights into these organisms are limited to cultivation-independent approaches such as single-cell- and metagenomics. However, without access to cultured representatives for verifying correct taxon-assignments, MDM genomes may cause potentially misleading conclusions based on misclassified or contaminant contigs, thereby obfuscating our view on the uncultured microbial majority. Moreover, gradual database contaminations by past genome submissions can cause error propagations which affect present as well as future comparative genome analyses. Consequently, strict contamination detection and filtering need to be applied, especially in the case of uncultured MDM genomes. Current genome reporting standards, however, emphasize completeness over purity and the de facto gold standard genome assessment tool, checkM, discriminates against uncultured taxa and fragmented genomes. To tackle these issues, we present a novel contig classification, screening, and filtering workflow and corresponding open-source python implementation called MDMcleaner, which was tested and compared to other tools on mock and real datasets. MDMcleaner revealed substantial contaminations overlooked by current screening approaches and sensitively detects misattributed contigs in both novel genomes and the underlying reference databases, thereby greatly improving our view on ‘microbial dark matter’.

Prevalence of bacterial genes in the phage fraction of food viromes

Antibiotic resistance genes (ARGs) have been identified in viral DNA isolated from different kinds of food, but little is known about their origin. In this study, twenty-one viromes were analyzed from samples of food previously reported to carry ARGs, including meat (poultry, veal, and pork), fish (Mediterranean, Atlantic, frozen, farmed and shellfish) and vegetables (lettuce, cucumber, and spinach). Classification of the contigs by Kraken revealed a large percentage of unclassified contigs (43.7–98.2%) in all the viromes. Only 0.05–7.1% of the contigs were identified as viral and of these, more than 91% belonged to different bacteriophage families, Podophages and Siphophages being the most prevalent. According to VirSorter, the largest number of viral contigs were derived from viromes of shellfish, followed by spinach. Spinach viromes also included the largest number of phage sequences identified by PHASTER. The abundant presence of bacterial genes in the viromes, including 16S rRNA genes, was attributed to the phage packaging of the bacterial genome fragments, as no bacterial DNA was found outside the viral capsids. The detection of 16S rRNA genes in the different viromes allowed diverse phage bacterial hosts to be identified. The three major functional groups of genes determined were related to metabolism, detoxification/resistance, and above all, biosynthesis. Various ARGs were quantified in the viromes by qPCR, the most prevalent being β-lactamases, particularly blaTEM. Analysis of ARG diversity in the viromes by Prokka and CARD revealed various resistance-related genes, whereas a more restrictive search by ResFinder identified blaTEM in all the food viromes, blaOXA in Atlantic fish-1 and spinach-2, oqxB in lettuce-1, and dfr in spinach-2. The presence of ARGs in the food viromes points to bacterial DNA mobilization by transduction mechanisms. Transduction of resistances by phage particles may therefore contribute to the emergence of resistant strains along the food chain and should be monitored.

Ancestral Chromosomes for Family Peronosporaceae Inferred from a Telomere-to-Telomere Genome Assembly of Peronospora effusa.

Downy mildew disease of spinach, caused by the oomycete Peronospora effusa, causes major losses to spinach production. In this study, the 17 chromosomes of P. effusa were assembled telomere-to-telomere, using Pacific Biosciences high-fidelity reads. Of these, 16 chromosomes are complete and gapless; chromosome 15 contains one gap bridging the nucleolus organizer region. This is the first telomere-to-telomere genome assembly for an oomycete. Putative centromeric regions were identified on all chromosomes. This new assembly enables a reevaluation of the genomic composition of Peronospora spp.; the assembly was almost double the size and contained more repeat sequences than previously reported for any Peronospora species. Genome fragments consistently underrepresented in six previously reported assemblies of P. effusa typically encoded repeats. Some genes annotated as encoding effectors were organized into multigene clusters on several chromosomes. Putative effectors were annotated on 16 of the 17 chromosomes. The intergenic distances between annotated genes were consistent with compartmentalization of the genome into gene-dense and gene-sparse regions. Genes encoding putative effectors were enriched in gene-sparse regions. The near-gapless assembly revealed apparent horizontal gene transfer from Ascomycete fungi. Gene order was highly conserved between P. effusa and the genetically oriented assembly of the oomycete Bremia lactucae; high levels of synteny were also detected with Phytophthora sojae. Extensive synteny between phylogenetically distant species suggests that many other oomycete species may have similar chromosome organization. Therefore, this assembly provides the foundation for genomic analyses of diverse oomycetes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Host-Associated Phages Disperse across the Extraterrestrial Analogue Antarctica.

Extreme Antarctic conditions provide one of the closest analogues of extraterrestrial environments. Since air and snow samples, especially from polar regions, yield DNA amounts in the lower picogram range, binning of prokaryotic genomes is challenging and renders studying the dispersal of biological entities across these environments difficult. Here, we hypothesized that dispersal of host-associated bacteriophages (adsorbed, replicating, or prophages) across the Antarctic continent can be tracked via their genetic signatures, aiding our understanding of virus and host dispersal across long distances. Phage genome fragments (PGFs) reconstructed from surface snow metagenomes of three Antarctic stations were assigned to four host genomes, mainly Betaproteobacteria, including Ralstonia spp. We reconstructed the complete genome of a temperate phage with nearly complete alignment to a prophage in the reference genome of Ralstonia pickettii 12D. PGFs from different stations were related to each other at the genus level and matched similar hosts. Metagenomic read mapping and nucleotide polymorphism analysis revealed a wide dispersal of highly identical PGFs, 13 of which were detected in seawater from the Western Antarctic Peninsula at a distance of 5,338 km from the snow sampling stations. Our results suggest that host-associated phages, especially of Ralstonia sp., disperse over long distances despite the harsh conditions of the Antarctic continent. Given that 14 phages associated with two R. pickettii draft genomes isolated from space equipment were identified, we conclude that Ralstonia phages are ideal mobile genetic elements to track dispersal and contamination in ecosystems relevant for astrobiology. IMPORTANCE Host-associated phages of the bacterium Ralstonia identified in snow samples can be used to track microbial dispersal over thousands of kilometers across the Antarctic continent, which functions as an extraterrestrial analogue because of its harsh environmental conditions. Due to the presence of these bacteria carrying genome-integrated prophages on space-related equipment and the potential for dispersal of host-associated phages demonstrated here, our work has implications for planetary protection, a discipline in astrobiology interested in preventing contamination of celestial bodies with alien biomolecules or forms of life.

Inactivation mechanisms of human adenovirus by e-beam irradiation in water environments.

This study aims to study the kinetics and mechanisms of human adenovirus inactivation by electron beam. Human adenovirus type 5 (HAdV-5) was inoculated in two types of aqueous substrates (phosphate-buffered saline - PBS, domestic wastewater - WW) treated by electron beam at a dose range between 3 and 21kGy. Samples were evaluated for virus infectivity, PCR amplification of fragments of HAdV-5 genome and abundance and antigenicity of the virion structural proteins. The maximum reduction in viral titre, in plaque-forming units (PFU) per millilitre, was about 7 and 5 log PFU/mL for e-beam irradiation at 20kGy in PBS and 19kGy in wastewater, respectively. Among the virion structural proteins detected, the hexon protein showed the higher radioresistance. Long (10.1 kbp) genomic DNA fragments were differently PCR amplified, denoting a substrate effect on HAdV-5 genome degradation by e-beam. The differences observed between the two substrates can be explained by the protective effect that the organic matter present in the substrate may have on viral irradiation. According to the obtained results, the decrease in viral viability/infectivity may be due to DNA damage and to protein alterations. In summary, electron beam irradiation at a dose of 13kGy is capable of reducing HAdV-5 viral titres by more than 99.99% (4 log PFU/mL) in both substrates assayed, indicating that this type of technology is effective for viral wastewater disinfection and may be used as a tertiary treatment in water treatment plants. KEY POINTS: • The substrate in which the virus is suspended has an impact on its sensitivity to e-beam treatment. • E-beam irradiation at 13 kGy is capable of reducing by 4 Log PFU/mL the HAdV-5 viral titre. • The decrease in viral viability/infectivity may be due to DNA damage and to protein alterations.

SunUp and Sunset genomes revealed impact of particle bombardment mediated transformation and domestication history in papaya.

Transgenic papaya is widely publicized for controlling papaya ringspot virus. However, the impact of particle bombardment on the genome remains unknown. The transgenic SunUp and its progenitor Sunset genomes were assembled into 351.5 and 350.3 Mb in nine chromosomes, respectively. We identified a 1.64 Mb insertion containing three transgenic insertions in SunUp chromosome 5, consisting of 52 nuclear-plastid, 21 nuclear-mitochondrial and 1 nuclear genomic fragments. A 591.9 kb fragment in chromosome 5 was translocated into the 1.64 Mb insertion. We assembled a gapless 9.8 Mb hermaphrodite-specific region of the Yh chromosome and its 6.0 Mb X counterpart. Resequencing 86 genomes revealed three distinct groups, validating their geographic origin and breeding history. We identified 147 selective sweeps and defined the essential role of zeta-carotene desaturase in carotenoid accumulation during domestication. Our findings elucidated the impact of particle bombardment and improved our understanding of sex chromosomes and domestication to expedite papaya improvement.