Single Genome Research Articles

COG database update 2024.

The Clusters of Orthologous Genes (COG) database, originally created in 1997, has been updated to reflect the constantly growing collection of completely sequenced prokaryotic genomes. This update increased the genome coverage from 1309 to 2296 species, including 2103 bacteria and 193 archaea, in most cases, with a single representative genome per genus. This set covers all genera of bacteria and archaea that included organisms with 'complete genomes' as per NCBI databases in November 2023. The number of COGs has been expanded from 4877 to 4981, primarily by including protein families involved in bacterial protein secretion. Accordingly, COG pathways and functional groups now include secretion systems of types II through X, as well as Flp/Tad and type IV pili. These groupings allow straightforward identification and examination of the prokaryotic lineages that encompass-or lack-a particular secretion system. Other developments include improved annotations for the rRNA and tRNA modification proteins, multi-domain signal transduction proteins, and some previously uncharacterized protein families. The new version of COGs is available at https://www.ncbi.nlm.nih.gov/research/COG, as well as on the NCBI FTP site https://ftp.ncbi.nlm.nih.gov/pub/COG/, which also provides archived data from previous COG releases.

Pangenome and pantranscriptome as the new reference for gene family characterisation – a case study of basic helix-loop-helix (bHLH) genes in barley

Genome-wide identification and comparative gene family analyses have been commonly performed to investigate species-specific evolution linked to various traits and molecular pathways. However, most previous studies were limited to gene screening in a single reference genome, failing to account for the gene presence/absence variations (gPAVs) in a species. Here, we propose an innovative pangenome-based approach of gene family analyses based on orthologous gene groups (OGGs). Using the basic helix-loop-helix (bHLH) transcription factor family in barley as an example, we identified 161 ∼ 176 bHLHs in 20 barley genomes, which could be classified into 201 OGGs. These 201 OGGs were further classified into 140 core, 12 soft-core, 29 shell, and 20 line-specific/cloud bHLHs, revealing a complete profile of bHLH in barley. Using a genome-scan approach, we overcome the genome annotation bias and identified on average 1.5 un-annotated core bHLHs per barley genome. We found that all core bHLHs belong to whole genome/segmental duplicates whilst dispensable bHLHs were more likely to result from small scale duplication events. Interestingly, we noticed that the dispensable bHLHs tended to enrich in specific subfamilies SF13, SF27, and SF28, implying the potential biased expansion of specific bHLHs in barley. We found that 50% of the bHLHs contain at least one intact transposon element within the 2kb upstream-to-downstream region. bHLHs with CNV have 1.48 TEs on average, significantly higher than 1.36 for core bHLH without CNV, supporting TEs’ potential role in bHLH expansion. Selection pressure analyses showed that dispensable bHLHs had experienced clear relaxed selection compared to core bHLHs, consistent with their conservation patterns. We further integrate pangenome with recently available barley pantranscriptome data in 5 tissues and discovered apparent transcriptional divergence within and across bHLH subfamilies. We conclude that pangenome-based gene family analyses can better describe the genuine evolution status of bHLHs untapped before and provided novel insights into bHLH evolution in barley. We expect this study will inspire similar analyses in many other gene families and species.

Pan-Draft: automated reconstruction of species-representative metabolic models from multiple genomes

The accurate reconstruction of genome-scale metabolic models (GEMs) for unculturable species poses challenges due to the incomplete and fragmented genetic information typical of metagenome-assembled genomes (MAGs). While existing tools leverage sequence homology from single genomes, this study introduces pan-Draft, a pan-reactome-based approach exploiting recurrent genetic evidence to determine the solid core structure of species-level GEMs. By comparing MAGs clustered at the species-level, pan-Draft addresses the issues due to the incompleteness and contamination of individual genomes, providing high-quality draft models and an accessory reactions catalog supporting the gapfilling step. This approach will improve our comprehension of metabolic functions of uncultured species.

Identification and characterization of chemotherapy resistant high-risk neuroblastoma persister cells.

Relapse rates in high-risk neuroblastoma remain exceedingly high. The malignant cells that are responsible for relapse have not been identified, and mechanisms of therapy resistance remain poorly understood. Here, we used single nucleus RNA sequencing and bulk whole genome sequencing to identify and characterize the residual malignant persister cells that survive chemotherapy from a cohort of 20 matched diagnosis and definitive surgery tumor samples from patients treated with high-risk neuroblastoma induction chemotherapy. We show that persister cells share common mechanisms of chemotherapy escape including suppression of MYCN activity and activation of NF-κB signaling, the latter is further enhanced by cell-cell communication between the malignant cells and the tumor microenvironment. Overall, our work dissects the transcriptional landscape of cellular persistence in high-risk neuroblastoma and paves the way to the development of new therapeutic strategies to prevent disease relapse.

Effects of aminolevulinic acid photodynamic therapy combined with antibiotics on Mycobacterium abscessus skin infections: an in vitro and in vivo study

BackgroundMycobacterium abscessus skin infections have emerged as a major medical issue. Traditional antibiotic treatments are challenging, prolonged, and often lead to recurrence, creating an urgent need for new therapies. This study investigates the effectiveness of aminolevulinic acid photodynamic therapy (ALA-PDT) combined with antibiotics in treatmenting M. abscessus, using both in vitro and in vivo methods. MethodsWe treated eight patients with M. abscessus skin infections following cosmetic surgery, using ALA-PDT (ALA concentration: 20%; red light: 80J/cm2) combined with oral or intravenous antibiotics,including clarithromycin, moxifloxacin and amikacin, to treat 8 patients with M. abscessus skin infection after medical cosmetic surgery, and assessed the treatment outcomes. Additionally, four bacterial strains (MAB-A1, MAB-A2, MAB-B1, and MAB-B2) isolated from patients were tested in vitro for ALA-PDT efficacy to determine the optimal ALA-PDT dosage. Furthermore, the strains’ single colony morphology, biofilm formation, and genome characteristics of were analyzed to explore the factors influencing ALA-PDT's bactericidal effects. Finally, a combined ALA-PDT and antibiotics sterilization experiment was conducted in vitro. ResultsClinically, ALA-PDT combined with antibiotics showed strong efficacy in treating M. abscessus skin infections, with no recurrence observed during follow-up. In vitro, ALA-PDT effectively killed M. abscessus, although MAB-B1 and MAB-B2 required a higher ALA-PDT dose compared with MAB-A1 and MAB-A2. Compared to MAB-A1 and MAB-A2, MAB-B1 and MAB-B2 exhibited stronger biofilm formation capabilities and bacterial virulence as well as genome mutations primarily affecting fatty acid synthesis and metabolism, potentially explaining their increased ALA-PDT dosage requirement. Notably, the combination of ALA-PDT and antibiotics exerted markedly higher bactericidal effects in vitro compared with antibiotics alone. ConclusionsALA-PDT combined with antibiotics emerged as an effective treatment for M. abscessus skin infections. However, optimal dosage and antibiotic combinations should be tailored to the characteristics of specific clinical strains.

Comprehensive Analysis of Complex Copy Number Variations in the Plasmodium falciparum Genome via Long-Read Sequencing: Implications for Antimalarial Resistance and Broad Genomic Insights

Abstract Introduction/Objective The genome of Plasmodium falciparum, characterized by a high AT content and a propensity for copy number variations (CNVs), presents a model for understanding how genomes adapt to stress. Our research focused on the prevalence and distribution of CNVs highlights long-read sequencing as a tool for detecting complex structural variants in other models. Methods/Case Report P. falciparum parasites were cultured and treated with DSM1 and aphidicolin that induce replication stress. Long-read sequencing was conducted using the Oxford Nanopore Technology Minion, with an emphasis on an optimizing DNA extraction method for preserving intact chromosomes. Visualization of CNVs, including tandem and inverted duplications, was achieved through a custom R Shiny application that allowed both quantitative and qualitative assessments. Results (if a Case Study enter NA) Long read sequencing identified an enrichment of multiple CNV types in treated samples, with a significant enrichment of inverted duplications as well as other types of structural variants. Aphidicolin and DSM1-treated samples demonstrated significant accumulation of these variations compared to controls (2.0 and 5.6-fold, respectively with a p value of < 0.0001 in both treatment groups). These occurred mainly in core genomic regions that do not contain variable gene copies. An increase in the number of inverted duplications may be related to stress response to the treatment and may be a potential adaptive mechanism of the parasite. Our manual single-read visualization technique made it possible to visualize the structure of various CNV types, which provides unique insight on CNVs as they arise in single genomes. Conclusion We hypothesize that inverted duplications represent stalled replication forks, and we are currently assessing their frequency in replicating and non-replicating cells. This study highlights the role of long-read sequencing in detecting complex genomic structures in P. falciparum, which may influence the parasite’s resistance to drug treatment. These findings not only advance our understanding of malaria pathology but also aid our comprehension of how cancer cells might react to drug-induced stress and the mechanisms behind CNV formation in neoplastic cells.

Differences in HIV-1 reservoir size, landscape characteristics and decay dynamics in acute and chronic treated HIV-1 Clade C infection.

Background Persisting HIV reservoir viruses in resting CD4 T cells and other cellular subsets are the main barrier to cure efforts. Antiretroviral therapy (ART) intensification by early initiation has been shown to enable post-treatment viral control in some cases but the underlying mechanisms are not fully understood. We hypothesized that ART initiated during the hyperacute phase of infection before peak will affect the size, decay dynamics and landscape characteristics of HIV-1 subtype C viral reservoirs. Methods We studied 35 women at high risk of infection from Durban, South Africa identified with hyperacute HIV infection by twice weekly testing for plasma HIV-1 RNA. Study participants included 11 who started ART at a median of 456 (297-1203) days post onset of viremia (DPOV), and 24 who started ART at a median of 1 (1-3) DPOV. We used peripheral blood mononuclear cells (PBMC) to measure total HIV-1 DNA by ddPCR and to sequence reservoir viral genomes by full length individual proviral sequencing (FLIP-seq) from onset of detection of HIV up to 1 year post treatment initiation. Results Whereas ART in hyperacute infection blunted peak viremia compared to untreated individuals (p<0.0001), there was no difference in total HIV-1 DNA measured contemporaneously (p=0.104). There was a steady decline of total HIV DNA in early treated persons over 1 year of ART (p=0.0004), with no significant change observed in the late treated group. Total HIV-1 DNA after one year of treatment was lower in the early treated compared to the late treated group (p=0.02). Generation of 697 single viral genome sequences revealed a difference in the longitudinal proviral genetic landscape over one year between untreated, late treated, and early treated infection: the relative contribution of intact genomes to the total pool of HIV-1 DNA after 1 year was higher in untreated infection (31%) compared to late treated (14%) and early treated infection (0%). Treatment initiated in both late and early infection resulted in a more rapid decay of intact (13% and 51% per month) versus defective (2% and 35% per month) viral genomes. However, intact genomes were still observed one year post chronic treatment initiation in contrast to early treatment where intact genomes were no longer detectable. Moreover, early ART reduced phylogenetic diversity of intact genomes and limited the seeding and persistence of cytotoxic T lymphocyte immune escape variants in the reservoir. Conclusions Overall, our results show that whereas ART initiated in hyperacute HIV-1 subtype C infection did not impact reservoir seeding, it was nevertheless associated with more rapid decay of intact viral genomes, decreased genetic complexity and immune escape in reservoirs, which could accelerate reservoir clearance when combined with other interventional strategies.

Revealing the Demographic History of the European Nightjar (Caprimulgus europaeus).

A species' demographic history gives important context to contemporary population genetics and a possible insight into past responses to climate change; with an individual's genome providing a window into the evolutionary history of contemporary populations. Pairwise sequentially Markovian coalescent (PSMC) analysis uses information from a single genome to derive fluctuations in effective population size change over the last ~5 million years. Here, we apply PSMC analysis to two European nightjar (Caprimulgus europaeus) genomes, sampled in Northwest and Southern Europe, with the aim of revealing the demographic history of nightjar in Europe. We successfully reconstructed effective population size over the last 5 million years. Our analysis shows that in response to global climate change, the effective population size of nightjar broadly increased under stable warm periods and decreased during cooler spans and prolonged glacial periods. PSMC analysis on the pseudo-diploid combination of the two genomes revealed fluctuations in gene flow between ancestral populations over time, with gene flow ceasing by the last-glacial period. Our results are tentatively suggestive of divergence in the European nightjar population, with timings consistent with differentiation being driven by restriction to different refugia during periods of glaciation. Finally, our results suggest that migratory behaviour in nightjar likely evolved prior to the last-glacial period, with long-distance migration seemingly persisting throughout the Pleistocene. However, further genetic structure analysis of individuals from known breeding sites across the species' contemporary range is needed to understand the extent and origins of range-wide differentiation in nightjar.

Leveraging transcriptome sequence read archives for virus detection in wild and colony populations of triatomines (Hemiptera: Reduviidae: Triatominae)

Triatomines are infamous as vectors of the parasite Trypanosoma cruzi, the causative agent of Chagas disease. However, climate-driven range expansion and urbanization adaptation of triatomine populations, coupled with their highly diverse feeding strategies (vertebrate haematophagy, kleptohaematophagy, and coprophagy), and has elevated interest in triatomines as potential arboviral vectors. Information on the triatomine virome is scant, with prior records including only eight insect-specific viruses: Triatoma virus (TrV) and Rhodnius prolixus viruses 1–7. Here, we leverage publicly available transcriptome datasets to assess viral diversity in 122 wild and colony kissing bugs representing eight species from six countries. In total, six viruses were detected (including Rhodnius prolixus viruses 4–6), and TrV was detected in almost half of all screened triatomines. This is the first report of TrV in Triatoma brasiliensis and in members of the genus Mepraia (M. gajardoi, M. spinolai, and M. parapatrica), and this effort has vastly expanded the publicly available genomic resources of TrV, adding 39 genome sequences to the single genome sequence currently available in the GenBank database. Furthermore, two additional viruses—Meccus longipennis virus 1 and Drosophila melanogaster Nora virus—are herein reported for the first time from kissing bugs. Meccus longipennis virus 1 was detected in Triatoma infestans from Argentina, Brazil, Chile, and Peru, and Drosophila melanogaster Nora virus was found in T. infestans from Argentina. Our results illustrate the advantage and utility of low-cost transcriptome data mining for the discovery of known and novel arboviruses in triatomines and other potential insect vectors.

Genealogical asymmetry under the IM model and a two-taxon test for gene flow.

Methods for detecting gene flow between populations often rely on asymmetry in the average length of particular genealogical branches, with the ABBA-BABA test being a well known example. Currently, asymmetry-based methods cannot be applied to a pair of populations and such analyses are instead performed using model-based methods. Here we investigate genealogical asymmetry under a two-population Isolation with Migration model. We focus on genealogies where the first coalescence event is between lineages sampled from different populations, as the external branches of these genealogies have equal expected length as long as there is no post-divergence gene flow. We show that unidirectional gene flow breaks this symmetry and results in the recipient population having longer external branches. We derive expectations for the probability of this genealogical asymmetry and propose a simple statistic (Am) to detect it from genome sequence data. Am provides a two-taxon test for gene flow that only requires a single unphased diploid genome from each population, with no outgroup information. We use analytic expectations and simulations to explore how recombination, unequal effective population sizes, bidirectional gene flow and background selection influence Am and find that the statistic provides unambiguous evidence for gene flow under a continent-island history. We estimate Am for genome sequence data from Heliconius butterflies and Odocoileus deer, generating results consistent with previous model-based analyses. Our work highlights a signal of gene flow overlooked to date and provides a method that complements existing approaches for investigating the demographic history of recently diverged populations.

Development of targeted whole genome sequencing approaches for Crimean-Congo haemorrhagic fever virus (CCHFV)

Crimean-Congo haemorrhagic fever (CCHF) is the most prevalent human tick-borne viral disease, with a reported case fatality rate of 30 % or higher. The virus contains a tri-segmented, negative-sense RNA genome consisting of the small (S), medium (M) and large (L) segments encoding respectively the nucleoprotein (NP), the glycoproteins precursor (GPC) and the viral RNA-dependent RNA polymerase (RDRP). CCHFV is one of the most genetically diverse arboviruses, with seven distinct lineages named after the region they were first reported in and based on S segment phylogenetic analysis.Due to the high genetic divergence of the virus, a single targeted tiling PCR strategy to enrich for viral nucleic acids prior to sequencing is difficult to develop, and previously we have developed and validated a tiling PCR enrichment method for the Europe 1 genetic lineage.We have developed a targeted, probe hybridisation capture method and validated its performance on clinical as well as cell-cultured material of CCHFV from different genetic lineages, including Europe 1, Europe 2, Africa 2 and Africa 3. The method produced over 95 % reference coverages with at least 10x sequencing depth. While we were only able to recover a single complete genome sequence from the tested Europe 1 clinical samples with the capture hybridisation protocol, the data provides evidence of its applicability to different CCHFV genetic lineages.CCHFV is an important tick-borne human pathogen with wide geographical distribution. Environmental as well as anthropogenic factors are causing increased CCHFV transmission. Development of strategies to recover CCHFV sequences from genetically diverse lineages of the virus is of paramount importance to monitor the presence of the virus in new areas, and in public health responses for CCHFV molecular surveillance to rapidly detect, diagnose and characterise currently circulating strains.

Linked regulation of genome integrity and senescence-associated inflammation by p53.

Genomic instability and inflammation are distinct hallmarks of aging, but the connection between them is poorly understood. Understanding their interrelationship will help unravel new mechanisms and therapeutic targets of aging and age-associated diseases. Here we report a novel mechanism directly linking genomic instability and inflammation in senescent cells through a mitochondria-regulated molecular circuit driven by p53 and cytoplasmic chromatin fragments (CCF). We show, through activation or inactivation of p53 by genetic and pharmacologic approaches, that p53 suppresses CCF accumulation and the downstream inflammatory senescence-associated secretory phenotype (SASP), without affecting cell cycle arrest. p53 activation suppressed CCF formation by promoting DNA repair, and this is reflected in maintenance of genomic integrity, particularly in subtelomeric regions, as shown by single cell genome resequencing. Activation of p53 in aged mice by pharmacological inhibition of MDM2 reversed signatures of aging, including age- and senescence-associated transcriptomic signatures of inflammation and age-associated accumulation of monocytes and macrophages in liver. Remarkably, mitochondria in senescent cells suppressed p53 activity by promoting CCF formation and thereby restricting ATM-dependent nuclear DNA damage signaling. These data provide evidence for a mitochondria-regulated p53 signaling circuit in senescent cells that controls DNA repair, genome integrity, and senescence- and age-associated inflammation. This pathway is immunomodulatory in mice and a potential target for healthy aging interventions by small molecules already shown to activate p53.

Assembly, comparative analysis, and utilization of a single haplotype reference genome for soybean.

Cultivar Williams 82 has served as the reference genome for the soybean research community since 2008, but is known to have areas of genomic heterogeneity among different sub-lines. This work provides an updated assembly (version Wm82.a6) derived from a specific sub-line known as Wm82-ISU-01 (seeds available under USDA accession PI 704477). The genome was assembled using Pacific BioSciences HiFi reads and integrated into chromosomes using HiC. The 20 soybean chromosomes assembled into a genome of 1.01Gb, consisting of 36 contigs. The genome annotation identified 48 387 gene models, named in accordance with previous assembly versions Wm82.a2 and Wm82.a4. Comparisons of Wm82.a6 with other near-gapless assemblies of Williams 82 reveal large regions of genomic heterogeneity, including regions of differential introgression from the cultivar Kingwa within approximately 30 Mb and 25 Mb segments on chromosomes 03 and 07, respectively. Additionally, our analysis revealed a previously unknown large (>20 Mb) heterogeneous region in the pericentromeric region of chromosome 12, where Wm82.a6 matches the 'Williams' haplotype while the other two near-gapless assemblies do not match the haplotype of either parent of Williams 82. In addition to the Wm82.a6 assembly, we also assembled the genome of 'Fiskeby III,' a rich resource for abiotic stress resistance genes. A genome comparison of Wm82.a6 with Fiskeby III revealed the nucleotide and structural polymorphisms between the two genomes within a QTL region for iron deficiency chlorosis resistance. The Wm82.a6 and Fiskeby III genomes described here will enhance comparative and functional genomics capacities and applications in the soybean community.

Single cell genome and epigenome co-profiling reveals hardwiring and plasticity in breast cancer.

Understanding the impact of genetic alterations on epigenomic phenotypes during breast cancer progression is challenging with unimodal measurements. Here, we report wellDA-seq, the first high-genomic resolution, high-throughput method that can simultaneously measure the whole genome and chromatin accessibility profiles of thousands of single cells. Using wellDA-seq, we profiled 22,123 single cells from 2 normal and 9 tumors breast tissues. By directly mapping the epigenomic phenotypes to genetic lineages across cancer subclones, we found evidence of both genetic hardwiring and epigenetic plasticity. In 6 estrogen-receptor positive breast cancers, we directly identified the ancestral cancer cells, and found that their epithelial cell-of-origin was Luminal Hormone Responsive cells. We also identified cell types with copy number aberrations (CNA) in normal breast tissues and discovered non-epithelial cell types in the microenvironment with CNAs in breast cancers. These data provide insights into the complex relationship between genetic alterations and epigenomic phenotypes during breast tumor evolution.

Metagenomic and genomic sequences from a nitrate-reducing benzene-degrading enrichment culture.

Metagenome-assembled genomes (MAGs) were recovered from metagenomic assemblies from a nitrate-reducing benzene-degrading enrichment culture. Ten MAGs of high quality or functional interest to benzene degradation are reported, seven of which are single contig genomes.

Pangenome Identification and Analysis of Terpene Synthase Gene Family Members in Gossypium.

Terpene synthases (TPSs), key gatekeepers in the biosynthesis of herbivore-induced terpenes, are pivotal in the diversity of terpene chemotypes across and within plant species. Here, we constructed a gene-based pangenome of the Gossypium genus by integrating the genomes of 17 diploid and 10 tetraploid species. Within this pangenome, 208 TPS syntelog groups (SGs) were identified, comprising 2 core SGs (TPS5 and TPS42) present in all 27 analyzed genomes, 6 softcore SGs (TPS11, TPS12, TPS13, TPS35, TPS37, and TPS47) found in 25 to 26 genomes, 131 dispensable SGs identified in 2 to 24 genomes, and 69 private SGs exclusive to a single genome. The mutational load analysis of these identified TPS genes across 216 cotton accessions revealed a great number of splicing variants and complex splicing patterns. The nonsynonymous/synonymous Ka/Ks value for all 52 analyzed TPS SGs was less than one, indicating that these genes were subject to purifying selection. Of 208 TPS SGs encompassing 1795 genes, 362 genes derived from 102 SGs were identified as atypical and truncated. The structural analysis of TPS genes revealed that gene truncation is a major mechanism contributing to the formation of atypical genes. An integrated analysis of three RNA-seq datasets from cotton plants subjected to herbivore infestation highlighted nine upregulated TPSs, which included six previously characterized TPSs in G. hirsutum (AD1_TPS10, AD1_TPS12, AD1_TPS40, AD1_TPS42, AD1_TPS89, and AD1_TPS104), two private TPSs (AD1_TPS100 and AD2_TPS125), and one atypical TPS (AD2_TPS41). Also, a TPS-associated coexpression module of eight genes involved in the terpenoid biosynthesis pathway was identified in the transcriptomic data of herbivore-infested G. hirsutum. These findings will help us understand the contributions of TPS family members to interspecific terpene chemotypes within Gossypium and offer valuable resources for breeding insect-resistant cotton cultivars.

First report of squash vein yellowing virus naturally infecting cucumber, squash and melon in Jordan.

Plant viruses are major restrictive pathogens to cucurbits production in Jordan. During field surveys conducted in September 2022 in the main cucurbit growing areas (desert area, Jordan Valley, and highlands), virus-like symptoms such as vein clearing, mosaic patterns, interveinal chlorosis, and fruit malformation were observed, in the presence of high whitefly populations (Bemisia tabaci MEAM1, Mdanat et al., 2022). A total of 80 leaf samples from different species [48 cucumber (Cucumis sativus), 11 squash (Cucurbita pepo), 14 melon (Cucumis melo) and 7 watermelon (Citrullus lanatus)], including 70 symptomatic and 10 asymptomatic samples, were collected for further investigations. Total RNA and DNA were extracted using RNeasy and DNeasy plant mini kit (QIAGEN), respectively, and molecular detection against an array of cucurbit-infecting viruses was conducted using protocols available in the literature (Suppl. Table). Squash vein yellowing virus (SqVYV; Potyviridae, Ipomovirus) was detected in combination with other cucurbit-infecting viruses, in 10 of 70 symptomatic samples with primers (SqYVV-v4762; SqYVV-c5512) targeting a portion of the cylindrical inclusion (CI) gene (Hernandez et al., 2021), including 3 cucumber, 4 squash and 3 melon samples, however, it was not detected in watermelon. Among other RNA viruses, cucumber green mottle mosaic virus, cucurbit yellow stunting disorder virus, cucurbit aphid-borne yellows virus, cucumber vein yellowing virus, and cucurbit chlorotic yellows virus were detected in 9, 34, 7, 18, and 23 samples, respectively. The DNA viruses, squash leaf curl virus and watermelon chlorotic stunt virus, were detected in 36 and 30 samples. None of the detected viruses were present in asymptomatic samples. All 80 samples were negative for watermelon mosaic virus and tomato leaf curl New Delhi virus. Cloning (pGEM T-Easy Vector; Promega), sequencing, and BLASTn analysis of 4 CI-specific cloned amplicons (~ 0.75 kb, GenBank Acc no. PP908660-PP908663) confirmed the identification of SqVYV, with highest BLASTn identity of ~91% and ~99% to isolates ESF3 from USA (MW584342) and SqVYV-IL from Israel (KT721735), respectively. To reconstruct the complete genome sequence of the SqVYV isolate from Jordan, total RNA from a pool of squash and melon plants was used to construct a cDNA library with the Illumina DNA Prep kit, which was sequenced on a NextSeq2000 instrument as paired reads (2x150 bp) at Leibniz Institute DSMZ, to generate 18,723,252 total reads. Bioinformatic analysis in Geneious (Biomatters) resulted in the assembly of a single genome sequence of 9,831 nt (GenBank Acc. no. PP916052), covered by 83,995 reads, with mean coverage of 1,197 (Geneious mapper, 10% maximum mismatches per read). The complete genome sequence and the deduced polyprotein sequence shared over 99% identity with SqVYV-IL from Israel. Mechanical inoculation of 10 cucumber (cv. Giant Global) and 10 squash (cv. Lebanese) plants with inoculum from infected cucumber, resulted in vein clearing and mosaic symptoms, after three weeks from inoculation, while no symptoms were observed in the six negative controls. Symptomatic plants were confirmed by PCR to be infected by SqVYV as described earlier. In this study, we report the occurrence of SqVYV for the first time in cucurbit crops in Jordan. SqVYV has been reported to cause large economic losses in cucurbits in USA and Israel, posing a major threat to watermelon growers (Adkins et al., 2007; 2008 a & b; 2013; Reingold et al., 2016; Webster et al., 2013). Our findings should encourage further studies on the incidence and prevalence of SqVYV in cucurbit and non-cucurbit crops, including weeds, to understand its epidemiology in Jordan, including its natural host range.

Pan-pathogen deep sequencing of nosocomial bacterial pathogens in Italy in spring 2020: a prospective cohort study

Nosocomial infections pose a considerable risk to patients who are susceptible, and this is particularly acute in intensive care units when hospital-associated bacteria are endemic. During the first wave of the COVID-19 pandemic, the surge of patients presented a significant obstacle to the effectiveness of infection control measures. We aimed to assess the risks and extent of nosocomial pathogen transmission under a high patient burden by designing a novel bacterial pan-pathogen deep-sequencing approach that could be integrated with standard clinical surveillance and diagnostics workflows. We did a prospective cohort study in a region of northern Italy that was severely affected by the first wave of the COVID-19 pandemic. Inpatients on both ordinary and intensive care unit (ICU) wards at the San Matteo hospital, Pavia were sampled on multiple occasions to identify bacterial pathogens from respiratory, nasal, and rectal samples. Diagnostic samples collected between April 7 and May 10, 2020were cultured on six different selective media designed to enrich for Acinetobacter baumannii, Escherichia coli, Enterococcus faecium, Enterococcus faecalis, Klebsiella pneumoniae, Pseudomonas aeruginosa, Staphylococcus aureus, and Streptococcus pneumoniae, and DNA from each plate with positive growth was deep sequenced en masse. We used mSWEEP and mGEMS to bin sequencing reads by sequence cluster for each species, followed by mapping with snippy to generate high quality alignments. Antimicrobial resistance genes were detected by use of ARIBA and CARD. Estimates of hospital transmission were obtained from pairwise bacterial single nucleotide polymorphism distances, partitioned by within-patient and between-patient samples. Finally, we compared the accuracy of our binned Acinetobacter baumannii genomes with those obtained by single colony whole-genome sequencing of isolates from the same hospital. We recruited patients from March 1to May 7, 2020. The pathogen population among the patients was large and diverse, with 2148species detections overall among the 2418sequenced samples from the 256patients. In total, 55sequence clusters from key pathogen species were detected at least five times. The antimicrobial resistance gene prevalence was correspondingly high, with key carbapenemase and extended spectrum ß-lactamase genes detected in at least 50 (40%) of 125patients in ICUs. Using high-resolution mapping to infer transmission, we established that hospital transmission was likely to be a significant mode of acquisition for each of the pathogen species. Finally, comparison with single colony Acinetobacter baumannii genomes showed that the resolution offered by deep sequencing was equivalent to single-colony sequencing, with the additional benefit of detection of co-colonisation of highly similar strains. Our study shows that a culture-based deep-sequencing approach is a possible route towards improving future pathogen surveillance and infection control at hospitals. Future studies should be designed to directly compare the accuracy, cost, and feasibility of culture-based deep sequencing with single colony whole-genome sequencing on a range of bacterial species. Wellcome Trust, European Research Council, Academy of Finland Flagship program, Trond Mohn Foundation, and Research Council of Norway.

Genome sequence of Lactobacillus johnsonii XZ17, isolated from a C57BL/6J mouse lung.

Lactobacillus johnsonii XZ17 was isolated from the lung homogenate of a healthy C57BL/6J mouse. XZ17 is diminished in the lungs of syngeneic bone marrow-transplanted recipient mice. Long-read sequencing of XZ17 yielded a single genome of 1,948,140 bp, with a GC content of 34.64%.

An enhanced Eco1 retron editor enables precision genome engineering in human cells from a single-copy integrated lentivirus.

Retrons are a retroelement class found in diverse prokaryotes that can be adapted to augment CRISPR-Cas9 genome engineering technology to efficiently rewrite short stretches of genetic information in bacteria and yeast; however, efficiency in human cells has been limited by unknown factors. We identified non-coding RNA (ncRNA) instability and impaired Cas9 activity as major contributors to poor retron editor efficiency. We re-engineered the Eco1 ncRNA to incorporate an exoribonuclease-resistant RNA pseudoknot from the Zika virus 3' UTR and devised an RNA processing strategy using Csy4 ribonuclease to liberate the sgRNA and ncRNA. These modifications yielded a ncRNA with 5'- and 3'-end protection and an sgRNA with minimal 5' extension. This strategy increased steady-state ncRNA levels and rescued Cas9 activity leading to enhanced efficiency of the Eco1 retron editor in human cells. The enhanced Eco1 retron editor enabled the insertion of missense mutations in human cells from a single integrated lentivirus, thereby ensuring genotype-phenotype linkage over multiple cell divisions. This work reveals a previously unappreciated role for ncRNA stability in retron editor efficiency in human cells. Here we present an enhanced Eco1 retron editor that enables efficient introduction of missense mutations in human cells from a single heritable genome copy.