Publications

3973

AbstractGenomes are fundamental to understanding microbial ecology and evolution. The emergence of high-throughput, long-read DNA sequencing has enabled recovery of microbial genomes from environmental samples at scale. However, expanding the microbial genome catalogue of soils and sediments has been challenging due to the enormous complexity of these environments. Here, we performed deep, long-read Nanopore sequencing of 154 soil and sediment samples collected across Denmark and through an optimised bioinformatics pipeline, we recovered genomes of 15,314 novel microbial species, including 4,757 high-quality genomes. The recovered microbial genomes span 1,086 novel genera and provide the first high-quality reference genomes for 612 previously known genera, expanding the phylogenetic diversity of the prokaryotic tree of life by 8 %. The long-read assemblies also enabled the recovery of thousands of complete rRNA operons, biosynthetic gene clusters and CRISPR-Cas systems, all of which were underrepresented and highly fragmented in previous terrestrial genome catalogues. Furthermore, the incorporation of the recovered MAGs into public genome databases significantly improved species-level classification rates for soil and sediment metagenomic datasets, thereby enhancing terrestrial microbiome characterization. With this study, we demonstrate that long-read sequencing and optimised bioinformatics, allows cost-effective recovery of high-quality microbial genomes from highly complex ecosystems, which remain the largest untapped source of biodiversity for expanding genome databases and filling in the gaps of the tree of life.

ABSTRACT “ Candidatus Liberibacter asiaticus” (CLas) is associated with citrus huanglongbing, a severe disease with global importance that affects citrus production in Brazil. This study reports the first complete genome of a Brazilian strain of CLas. The genomic structure comparison of strain 9PA with those of 13 complete CLas genomes revealed 9,091 mismatches and 992 gaps/insertions, highlighting eight locally colinear blocks, among which six are in the prophage region. Phylogenetic analysis categorized 13 CLas genomes into two clusters with 9PA clustered with strains from China and the United States. Whole-genomic comparison identified diverse hypervariable genomic regions (HGRs). Three HGRs in the chromosomal region and three in the prophage region were selected and investigated by polymerase chain reaction. HGRs assessed from 68 samples, from medium- to high-huanglongbing incidence areas in Sao Paulo state, were grouped into haplotypes A to P. Haplotype A, which includes strain 9PA, is the second most prevalent, representing 19.1% of the samples. Haplotype B, the most common, accounts for 42.6%. Together with haplotype C, these make up 72% of the evaluated samples. The 9PA strain has prophage P-9PA-1, both integrated and circularized, and P-9PA-3, only found in a circularized form. Prophages show high identity with SC1 (83%) and P-JXGC-3 (98%). Co-occurrence of both type 1 and 3 prophages was observed in field samples. The approach employed provides insights into the Brazilian CLas population, providing markers for population studies and highlighting the prevalence of type 1 and 3 prophages in the population. IMPORTANCE CLas is a destructive pathogen responsible for causing the severe citrus disease known as huanglongbing. Our study presents the first fully sequenced Brazilian strain of CLas, designated as 9PA, and includes an analysis of two prophages occurring in this strain. The main objective of our research was to compare the genome features of this Brazilian strain with other fully sequenced genomes and to identify its hypervariable genetic regions. These regions were subsequently used to assess genomic variability within both the chromosomal and prophage regions in Brazilian isolates of CLas. Our findings offer valuable insights into the diversified adaptation of CLas.

ABSTRACTThe Oil Refinery (OR) consortium is a model methanogenic enrichment culture for studying anaerobic benzene degradation. Over 80% of the culture’s bacterial community is comprised of two closely related strains of benzene-fermentingDesulfobacterota(designated ORM2a and ORM2b) whose mechanism of benzene degradation is unknown. Two new metagenomes, including a fully closed metagenome-assembled genome (MAG) for ORM2a, enabled a thorough investigation of this culture’s proteome. Among the proteins identified were Bam-like subunits of an ATP-independent benzoyl-CoA degradation pathway and associated downstream beta-oxidation proteins producing acetyl-CoA. The most abundant proteins identified mapped to two ORM2 gene clusters of unknown function. Syntenic gene clusters were identified in one other known benzene degrader,Pelotomaculumcandidate BPL, as well as a handful of contigs assembled from hydrothermal vent metagenomes. Extensive searches against reference sequence and structural databases indicate that the first (“Magic”) gene cluster likely catalyzes the chemically difficult benzene activation step. The second (“Nanopod”) gene cluster is predicted to code for an efflux system that pumps excess benzene out of cells, mitigating some of its toxigenic effects. Phylogenomic analyses place ORM2a and ORM2b within a novel genus of benzene-degrading specialists which we propose naming “CandidatusBenzenivorax”. We hope to engage the research community to help in confirming the roles of the proteins in the “Magic” and “Nanopod” gene clusters, and to search through their own cultures for these features.

Psyllids are major vectors of plant diseases, including Candidatus Liberibacter solanacearum (CLso), the bacterial agent associated with 'zebra chip' disease in potatoes and 'carrot yellows' disease in carrot. Despite their agricultural significance, there is limited knowledge on the genome structure and genetic diversity of psyllids. In this study, we provide chromosome-level genome assemblies for three psyllid species known to transmit CLso: Dyspersa apicalis (carrot psyllid), Dyspersa pallida, and Trioza urticae (nettle psyllid). As D. apicalis is recognised as the primary vector of CLso by carrot growers in Northern Europe, we also resequenced populations of this species from Finland, Norway, and Austria. Genome assemblies were constructed using PacBio HiFi and Hi-C sequencing data, yielding genome sizes of: 594.01 Mbp for D. apicalis; 587.80 Mbp for D. pallida; and 655.58 Mbp for T. urticae. Over 90% of sequences anchored into 13 pseudo-chromosomes per species. The assemblies for D. apicalis and D. pallida exhibited high completeness, capturing over 92% of conserved Hemiptera single-copy orthologues, as assessed by Benchmarking Universal Single-Copy Orthologues (BUSCO) analysis. Furthermore, we identified sequences of the primary psyllid symbiont, Candidatus Carsonella ruddii, in all three species. Comparative genomic analyses demonstrated synteny with other psyllid species. Notably, we observed significant expansions in gene families, particularly those linked to potential insecticide detoxification, within the Dyspersa lineage. Resequencing efforts also revealed the existence of multiple subpopulations of D. apicalis across Europe. These high-quality genome resources will support future research on genome evolution, insect-plant-pest interactions, and strategies for disease management.

ABSTRACTThe Coxiellaceae bacterial family, within the order Legionellales, is defined by a collection of poorly characterized obligate intracellular bacteria. The zoonotic pathogen and causative agent of human Q fever,Coxiella burnetii, represents the best characterized member of this family. Coxiellaceae establish replicative niches within diverse host cells and rely on their host for survival, making them challenging to isolate and cultivate within a laboratory setting. Here we describe a new genus within the Coxiellaceae family that has been previously shown to infect economically significant freshwater crayfish. Using culture-independent long-read metagenomics, we reconstructed the complete genome of this novel organism and demonstrate that the previously referred to as CandidatusCoxiella cheraxirepresents a novel genus within this family, herein denoted CandidatusParacoxiella cheracis. Interestingly, we demonstrate that CandidatusP. cheracisencodes for a complete, putatively functional Dot/Icm type 4 secretion system that likely mediates the intracellular success of this pathogen. In silico analysis defined a unique repertoire of Dot/Icm effector proteins and highlighted homologues of several importantC. burnetiieffectors including a homologue of CpeB that was demonstrated to be a Dot/Icm substrate inC. burnetii.IMPORTANCEUsing long-read sequencing technology we have uncovered the full genome sequence of CandidatusParacoxiellacheracis, a pathogen of economic importance in aquaculture. Analysis of this sequence has revealed new insight into this novel member of the Coxiellaceae family, demonstrating that it represents a new genus within this poorly characterized family of intracellular organisms. Importantly, the genome sequence reveals invaluable information that will support diagnostics and potentially both preventative and treatment strategies within crayfish breeding facilities. CandidatusP. cheracisalso represents a new member of Dot/Icm pathogens that rely on this system to establish an intracellular niche. CandidatusP. cheracispossesses a unique cohort of putative Dot/Icm substrates that constitute a collection of new eukaryotic cell biology manipulating effector proteins.

Abstract  The Gemmatimonadota phylum has been widely detected in diverse natural environments, yet their specific ecological roles in many habitats remain poorly investigated. Similarly, the Candidatus ARS69 phylum has been identified only in a few habitats, and literature on their metabolic functions is relatively scarce. In the present study, we investigated the ecological significance of phyla Ca. ARS69 and Gemmatimonadota in the Arctic glacier foreland (GF) ecosystems through genome-resolved metagenomics. We have reconstructed the first high-quality metagenome-assembled genome (MAG) belonging to Ca. ARS69 and 12 other MAGs belonging to phylum Gemmatimonadota from the three different Arctic GF samples. We further elucidated these two groups phylogenetic lineage and their metabolic function through phylogenomic and pangenomic analysis. The analysis showed that all the reconstructed MAGs potentially belonged to novel species. The MAGs belonged to Ca. ARS69 consist about 8296 gene clusters, of which only about 8% of single-copy core genes (n = 980) were shared among them. The study also revealed the potential ecological role of Ca. ARS69 is associated with carbon fixation, denitrification, sulfite oxidation, and reduction biochemical processes in the GF ecosystems. Similarly, the study demonstrates the widespread distribution of different classes of Gemmatimonadota across wide ranges of ecosystems and their metabolic functions, including in the polar region. Key points • Glacier foreland ecosystems act as a natural laboratory to study microbial community structure. • We have reconstructed 13 metagenome-assembled genomes from the soil samples. • All the reconstructed MAGs belonged to novel species with different metabolic processes. • Ca. ARS69 and Gemmatimonadota MAGs were found to participate in carbon fixation and denitrification processes.

Cottony ash psyllid (CAP, Psyllopsis discrepans) is an important, invasive insect pest of ash trees in North America, where it has established populations and is the host of a newly identified strain of ‘ Candidatus Liberibacter solanacearum’. However, not much is known about the diversity of its introduced population. In this study, a CAP mitochondrial genome (mitogenome) sequence was obtained from a collection in Saskatoon, Saskatchewan, Canada. The CAP mitogenome is a circular DNA of 18,824 bp, encoding 13 protein-coding genes, 21 transfer RNA genes, and two ribosomal RNA genes. BLAST search using the CAP mitogenome as a query against the GenBank sequence database showed the mitogenome of Euphyllura phillyreae (15,202 bp) was the most similar (query coverage = 77%; percentage identity = 78.90%). The CAP mitogenome is significantly different from other known psyllid mitogenomes with the presence of a 4,357-bp control region. The mitogenome sequence will further the genomic understanding of CAP. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 “No Rights Reserved” license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2024.