Publications

4358

Abstract Background The RCA (Roseobacter clade affiliated) cluster belongs to the family Roseobacteracea and represents a major Roseobacter lineage in temperate to polar oceans. Despite its prevalence and abundance, only a few genomes and one described species, Planktomarina temperata, exist. To gain more insights into our limited understanding of this cluster and its taxonomic and functional diversity and biogeography, we screened metagenomic datasets from the global oceans and reconstructed metagenome-assembled genomes (MAG) affiliated to this cluster. Results The total of 82 MAGs, plus five genomes of isolates, reveal an unexpected diversity and novel insights into the genomic features, the functional diversity, and greatly refined biogeographic patterns of the RCA cluster. This cluster is subdivided into three genera: Planktomarina, Pseudoplanktomarina, and the most deeply branching Candidatus Paraplanktomarina. Six of the eight Planktomarina species have larger genome sizes (2.44–3.12 Mbp) and higher G + C contents (46.36–53.70%) than the four Pseudoplanktomarina species (2.26–2.72 Mbp, 42.22–43.72 G + C%). Cand. Paraplanktomarina is represented only by one species with a genome size of 2.40 Mbp and a G + C content of 45.85%. Three novel species of the genera Planktomarina and Pseudoplanktomarina are validly described according to the SeqCode nomenclature for prokaryotic genomes. Aerobic anoxygenic photosynthesis (AAP) is encoded in three Planktomarina species. Unexpectedly, proteorhodopsin (PR) is encoded in the other Planktomarina and all Pseudoplanktomarina species, suggesting that this light-driven proton pump is the most important mode of acquiring complementary energy of the RCA cluster. The Pseudoplanktomarina species exhibit differences in functional traits compared to Planktomarina species and adaptations to more resource-limited conditions. An assessment of the global biogeography of the different species greatly expands the range of occurrence and shows that the different species exhibit distinct biogeographic patterns. They partially reflect the genomic features of the species. Conclusions Our detailed MAG-based analyses shed new light on the diversification, environmental adaptation, and global biogeography of a major lineage of pelagic bacteria. The taxonomic delineation and validation by the SeqCode nomenclature of prominent genera and species of the RCA cluster may be a promising way for a refined taxonomic identification of major prokaryotic lineages and sublineages in marine and other prokaryotic communities assessed by metagenomics approaches.

Huanglongbing (HLB) is a devastating citrus disease caused by Candidatus Liberibacter asiaticus (CLas). Since its initial outbreak in Guangdong Province, China, it has spread to 10 provinces and caused significant economic losses. Hence, assessing CLas genetic diversity and demographic history is crucial for HLB epidemic prevention and control. In this study, we collected 500 leaf samples of CLas-infected plants from 10 provinces. We performed multi-loci sequence analysis on four gene fragments (omp, DnaA, GroEL, and SDE1) to explore the genetic differentiation and diversity of CLas in China. Our results indicated low nucleotide diversity (0.00005 ± 0.00001) in CLas, with the absence of significant systematic geographic structure in its distribution. Molecular variance analysis revealed predominant (81.7%) genetic variations within the population, with a minor variation (18.3%) occurring between populations as well as Yunnan provinces. In the Fujian population, significant gene exchange occurred with the other nine populations. Significant negative values in Tajima’s D and Fu’s FS neutrality tests indicated historical population expansions. The nucleotide mismatch distribution curve exhibits a single peak pattern, further supporting the expansion events. Our findings hold potential for advancing epidemiological research and providing suggestions for effective strategies to mitigate the spread of CLas and control HLB.

AbstractThe roles of Asgard archaea in eukaryogenesis and marine biogeochemical cycles are well studied, yet their contributions in soil ecosystems are unknown. Of particular interest are Asgard archaeal contributions to methane cycling in wetland soils. To investigate this, we reconstructed two complete genomes for soil-associated Atabeyarchaeia, a new Asgard lineage, and the first complete genome of Freyarchaeia, and defined their metabolismin situ. Metatranscriptomics highlights high expression of [NiFe]-hydrogenases, pyruvate oxidation and carbon fixation via the Wood-Ljungdahl pathway genes. Also highly expressed are genes encoding enzymes for amino acid metabolism, anaerobic aldehyde oxidation, hydrogen peroxide detoxification and glycerol and carbohydrate breakdown to acetate and formate. Overall, soil-associated Asgard archaea are predicted to be non-methanogenic acetogens, likely impacting reservoirs of substrates for methane production in terrestrial ecosystems.One-Sentence SummaryComplete genomes of Asgard archaea, coupled with metatranscriptomic data, indicate roles in production and consumption of carbon compounds that are known to serve as substrates for methane production in wetlands.

ABSTRACTMembers of the “CandidatusAccumulibacter” genus are widely studied as key polyphosphate-accumulating organisms (PAOs) in biological nutrient removal (BNR) facilities performing enhanced biological phosphorus removal (EBPR). This diverse lineage includes 18 “Ca. Accumulibacter” species, which have been proposed based on the phylogenetic divergence of the polyphosphate kinase 1 (ppk1) gene and genome-scale comparisons of metagenome-assembled genomes (MAGs). Phylogenetic classification based on the 16S rRNA genetic marker has been difficult to attain because most “Ca. Accumulibacter” MAGs are incomplete and often do not include the rRNA operon. Here, we investigate the “Ca. Accumulibacter” diversity in pilot-scale treatment trains performing BNR under low dissolved oxygen (DO) conditions using genome-resolved metagenomics. Using long-read sequencing, we recovered medium and high-quality MAGs for 5 of the 18 “Ca. Accumulibacter” species, all with rRNA operons assembled, which allowed a reassessment of the 16S rRNA-based phylogeny of this genus and an analysis of phylogeny based on the 23S rRNA gene. In addition, we recovered a cluster of MAGs that based on 16S rRNA, 23S rRNA,ppk1, and genome-scale phylogenetic analyses do not belong to any of the currently recognized “Ca. Accumulibacter” species for which we propose the new species designation “Ca. Accumulibacter jenkinsii” sp. nov. Relative abundance evaluations of the genus across all pilot plant operations revealed that regardless of the operational mode, “Ca. A. necessarius” and “Ca. A. propinquus” accounted for more than 40% of the “Ca. Accumulibacter” community, whereas the newly proposed “Ca. A. jenkinsii” accounted for about 5% of the “Ca. Accumulibacter” community.IMPORTANCEOne of the main drivers of energy use and operational costs in activated sludge processes is the amount of oxygen provided to enable biological phosphorus and nitrogen removal. Wastewater treatment facilities are increasingly considering reduced aeration to decrease energy consumption, and whereas successful BNR has been demonstrated in systems with minimal aeration, an adequate understanding of the microbial communities that facilitate nutrient removal under these conditions is still lacking. In this study, we used genome-resolved metagenomics to evaluate the diversity of the “CandidatusAccumulibacter” genus in pilot-scale plants operating with minimal aeration. We identified the “Ca.Accumulibacter” species enriched under these conditions, including one novel species for which we propose “Ca.Accumulibacter jenkinsii” sp. nov. as its designation. Furthermore, the MAGs obtained for 5 additional “Ca.Accumulibacter” species further refine the phylogeny of the “Ca.Accumulibacter” genus and provide new insight into its diversity within unconventional biological nutrient removal systems.

AbstractCandidatusAltiarchaea are widespread across aquatic subsurface ecosystems and possess a highly conserved core genome, yet adaptations of this core genome to different biotic and abiotic factors based on gene expression remain unknown. Here, we investigated the metatranscriptome of twoCa. Altiarchaeum populations that thrive in two substantially different subsurface ecosystems. In Crystal Geyser, a high-CO2groundwater system in the USA,Ca. Altiarchaeum crystalense co-occurs with the symbiontCa. Huberiarchaeum crystalense, while in the Muehlbacher sulfidic spring in Germany, an artesian spring high in sulfide concentration,Ca. A. hamiconexum is heavily infected with viruses. We here mapped metatranscriptome reads against their genomes to analyze thein situexpression profile of their core genomes. Out of 537 shared gene clusters, 331 were functionally annotated and 130 differed significantly in expression between the two sites. Main differences were related to genes involved in cell defense like CRISPR-Cas, virus defense, replication, and transcription as well as energy and carbon metabolism. Our results demonstrate that altiarchaeal populations in the subsurface are likely adapted to their environment while influenced by other biological entities that tamper with their core metabolism. We consequently posit that viruses and symbiotic interactions can be major energy sinks for organisms in the deep biosphere.(Originality-Significance StatementOrganisms of the uncultivated phylumCa. Altiarchaeota are globally widespread and fulfill essential roles in carbon cycling,e.g., carbon fixation in the continental subsurface. Here, we show that the transcriptional activity of organisms in the continental subsurface differ significantly depending on the geological and microbial setting of the ecosystem explaining many of the previously observed physiological traits of this organism group.)

The methanogenic strain Mx-05T was isolated from the human fecal microbiome. A phylogenetic analysis based on the 16S rRNA gene and protein marker genes indicated that the strain is affiliated with the order Methanomassiliicoccales. It shares 86.9% 16S rRNA gene sequence identity with Methanomassiliicoccus luminyensis, the only member of this order previously isolated. The cells of Mx-05T were non-motile cocci, with a diameter range of 0.4–0.7 μm. They grew anaerobically and reduced methanol, monomethylamine, dimethylamine, and trimethylamine into methane, using H2 as an electron donor. H2/CO2, formate, ethanol, and acetate were not used as energy sources. The growth of Mx-05T required an unknown medium factor(s) provided by Eggerthella lenta and present in rumen fluid. Mx-05T grew between 30 °C and 40 °C (optimum 37 °C), over a pH range of 6.9–8.3 (optimum pH 7.5), and between 0.02 and 0.34 mol.L−1 NaCl (optimum 0.12 mol.L−1 NaCl). The genome is 1.67 Mbp with a G+C content of 55.5 mol%. Genome sequence annotation confirmed the absence of the methyl branch of the H4MPT Wood–Ljungdahl pathway, as described for other Methanomassiliicoccales members. Based on an average nucleotide identity analysis, we propose strain Mx-05T as being a novel representative of the order Methanomassiliicoccales, within the novel family Methanomethylophilaceae, for which the name Methanomethylophilus alvi gen. nov, sp. nov. is proposed. The type strain is Mx-05T (JCM 31474T).

Maternally inherited obligate endosymbionts codiverge with their invertebrate hosts and reflect their host’s evolutionary history. Whiteflies (Hemiptera: Aleyrodidae) harbor one obligate endosymbiont, Candidatus Portiera aleyrodidarum (hereafter Portiera). Portiera was anciently acquired by whitefly and has been coevolving with its host ever since. Uncovering the divergence of endosymbionts provides a fundamental basis for inspecting the coevolutionary processes between the bacteria and their hosts. To illustrate the divergence of Portiera lineages across different whitefly species, we sequenced the Portiera genome from Aleyrodes shizuokensis and conducted a comparative analysis on the basic features and gene evolution with bacterial genomes from five whitefly genera, namely Aleurodicus, Aleyrodes, Bemisia, Pealius, and Trialeurodes. The results indicated that Portiera from Bemisia possessed significantly larger genomes, fewer coding sequences (CDSs), and a lower coding density. Their gene arrangement differed notably from those of other genera. The phylogeny of the nine Portiera lineages resembled that of their hosts. Moreover, the lineages were classified into three distinct genetic groups based on the genetic distance, one from Aleurodicus (Aleurodicinae), one from Bemisia (Aleyrodinae), and another from Aleyrodes, Pealius, and Trialeurrodes (Aleyrodinae). Synonymous and nonsynonymous rate analyses, parity rule 2 plot analyses, neutrality plot analyses, and effective number of codons analyses supported the distinction of the three genetic groups. Our results indicated that Portiera from distant hosts exhibit distinct genomic contents, implying codivergence between hosts and their endosymbionts. This work will enhance our understanding of coevolution between hosts and their endosymbionts.

The familyHalomonadaceaeis the largest family composed of halophilic bacteria, with more than 160 species with validly published names as of July 2023. Several classifications to circumscribe this family are available in major resources, such as those provided by the List of Prokaryotic names with Standing in Nomenclature (LPSN), NCBI Taxonomy, Genome Taxonomy Database (GTDB), and Bergey’s Manual of Systematics of Archaea and Bacteria (BMSAB), with some degree of disagreement between them. Moreover, regardless of the classification adopted, the genusHalomonasis not phylogenetically consistent, likely because it has been used as a catch-all for newly described species within the familyHalomonadaceaethat could not be clearly accommodated in otherHalomonadaceaegenera. In the past decade, some taxonomic rearrangements have been conducted on theHalomonadaceaebased on ribosomal and alternative single-copy housekeeping gene sequence analysis. High-throughput technologies have enabled access to the genome sequences of many type strains belonging to the familyHalomonadaceae; however, genome-based studies specifically addressing its taxonomic status have not been performed to date. In this study, we accomplished the genome sequencing of 17 missing type strains ofHalomonadaceaespecies that, together with other publicly available genome sequences, allowed us to re-evaluate the genetic relationship, phylogeny, and taxonomy of the species and genera within this family. The approach followed included the estimate of the Overall Genome Relatedness Indexes (OGRIs) such as the average amino acid identity (AAI), phylogenomic reconstructions using amino acid substitution matrices customized for the familyHalomonadaceae, and the analysis of clade-specific signature genes. Based on our results, we conclude that the genusHalovibriois obviously out of place within the familyHalomonadaceae,and, on the other hand, we propose a division of the genusHalomonasinto seven separate genera and the transfer of seven species fromHalomonasto the genusModicisalibacter, together with the emendation of the latter. Additionally, data from this study demonstrate the existence of various synonym species names in this family.

Methane emission by terrestrial invertebrates is restricted to millipedes, termites, cockroaches, and scarab beetles. The arthropod-associated archaea known to date belong to the orders Methanobacteriales, Methanomassiliicoccales, Methanomicrobiales, and Methanosarcinales, and in a few cases also to non-methanogenic Nitrososphaerales and Bathyarchaeales. However, all major host groups are severely undersampled, and the taxonomy of existing lineages is not well developed. Full-length 16S rRNA gene sequences and genomes of arthropod-associated archaea are scarce, reference databases lack resolution, and the names of many taxa are either not validly published or under-classified and require revision. Here, we investigated the diversity of archaea in a wide range of methane-emitting arthropods, combining phylogenomic analysis of isolates and metagenome-assembled genomes (MAGs) with amplicon sequencing of full-length 16S rRNA genes. Our results allowed us to describe numerous new species in hitherto undescribed taxa among the orders Methanobacteriales (Methanacia, Methanarmilla, Methanobaculum, Methanobinarius, Methanocatella, Methanoflexus, Methanorudis, and Methanovirga, all gen. nova), Methanomicrobiales (Methanofilum and Methanorbis, both gen. nova), Methanosarcinales (Methanofrustulum and Methanolapillus, both gen. nova), Methanomassiliicoccales (Methanomethylophilaceae fam. nov., Methanarcanum, Methanogranum, Methanomethylophilus, Methanomicula, Methanoplasma, Methanoprimaticola, all gen. nova), and the new family Bathycorpusculaceae (Bathycorpusculum gen. nov.). Reclassification of amplicon libraries from this and previous studies using this new taxonomic framework revealed that arthropods harbor only CO2 and methyl-reducing hydrogenotrophic methanogens. Numerous genus-level lineages appear to be present exclusively in arthropods, suggesting long evolutionary trajectories with their termite, cockroach, and millipede hosts, and a radiation into various microhabitats and ecological niches provided by their digestive tracts (e.g., hindgut compartments, gut wall, or anaerobic protists). The distribution patterns among the different host groups are often complex, indicating a mixed mode of transmission and a parallel evolution of invertebrate and vertebrate-associated lineages.

As the name of the genus Pantoea (“of all sorts and sources”) suggests, this genus includes bacteria with a wide range of provenances, including plants, animals, soils, components of the water cycle, and humans. Some members of the genus are pathogenic to plants, and some are suspected to be opportunistic human pathogens; while others are used as microbial pesticides or show promise in biotechnological applications. During its taxonomic history, the genus and its species have seen many revisions. However, evolutionary and comparative genomics studies have started to provide a solid foundation for a more stable taxonomy. To move further toward this goal, we have built a 2,509-gene core genome tree of 437 public genome sequences representing the currently known diversity of the genus Pantoea. Clades were evaluated for being evolutionarily and ecologically significant by determining bootstrap support, gene content differences, and recent recombination events. These results were then integrated with genome metadata, published literature, descriptions of named species with standing in nomenclature, and circumscriptions of yet-unnamed species clusters, 15 of which we assigned names under the nascent SeqCode. Finally, genome-based circumscriptions and descriptions of each species and each significant genetic lineage within species were uploaded to the LINbase Web server so that newly sequenced genomes of isolates belonging to any of these groups could be precisely and accurately identified.