TIGRFAMs: a protein family resource for the functional identification of proteins

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

Darobactin Substrate Engineering and Computation Show Radical Stability Governs Ether versus C-C Bond Formation

The Gram-negative selective antibiotic darobactin A has attracted interest owing to its intriguing fused bicyclic structure and unique targeting of the outer membrane protein BamA. Darobactin, a ribosomally synthesized and post-translationally modified peptide (RiPP), is produced by a radical S-adenosyl methionine (rSAM)-dependent enzyme (DarE) and contains one ether and one C-C cross-link. Herein, we analyze the substrate tolerance of DarE and describe an underlying catalytic principle of the enzyme. These efforts produced 51 enzymatically modified darobactin variants, revealing that DarE can install the ether and C-C cross-links independently and in different locations on the substrate. Notable variants with fused bicyclic structures were characterized, including darobactin W3Y, with a non-Trp residue at the twice-modified central position, and darobactin K5F, which displays a fused diether ring pattern. While lacking antibiotic activity, quantum mechanical modeling of darobactins W3Y and K5F aided in the elucidation of the requisite features for high-affinity BamA engagement. We also provide experimental evidence for β-oxo modification, which adds support for a proposed DarE mechanism. Based on these results, ether and C-C cross-link formation was investigated computationally, and it was determined that more stable and longer-lived aromatic Cβ radicals correlated with ether formation. Further, molecular docking and transition state structures based on high-level quantum mechanical calculations support the different indole connectivity observed for ether (Trp-C7) and C-C (Trp-C6) cross-links. Finally, mutational analysis and protein structural predictions identified substrate residues that govern engagement to DarE. Our work informs on darobactin scaffold engineering and further unveils the underlying principles of rSAM catalysis.

Proteomic profiling of antibiotic-resistant Escherichia coli GW-AmxH19 isolated from hospital wastewater treated with physical plasma

Microorganisms which are resistant to antibiotics are a global threat to the health of humans and animals. Wastewater treatment plants are known hotspots for the dissemination of antibiotic resistances. Therefore, novel methods for the inactivation of pathogens, and in particular antibiotic-resistant microorganisms (ARM), are of increasing interest. An especially promising method could be a water treatment by physical plasma which provides charged particles, electric fields, UV-radiation, and reactive species. The latter are foremost responsible for the antimicrobial properties of plasma. Thus, with plasma it might be possible to reduce the amount of ARM and to establish this technology as additional treatment stage for wastewater remediation. However, the impact of plasma on microorganisms beyond a mere inactivation was analyzed in more detail by a proteomic approach. Therefore, Escherichia coli GW-AmxH19, isolated from hospital wastewater in Germany, was used. The bacterial solution was treated by a plasma discharge ignited between each of four pins and the liquid surface. The growth of E. coli and the pH-value decreased during plasma treatment in comparison with the untreated control. Proteome and antibiotic resistance profile were analyzed. Concentrations of nitrite and nitrate were determined as long-lived indicative products of a transient chemistry associated with reactive nitrogen species (RNS). Conversely, hydrogen peroxide served as indicator for reactive oxygen species (ROS). Proteome analyses revealed an oxidative stress response as a result of plasma-generated RNS and ROS as well as a pH-balancing reaction as key responses to plasma treatment. Both, the generation of reactive species and a decreased pH-value is characteristic for plasma-treated solutions. The plasma-mediated changes of the proteome are discussed also in comparison with the Gram-positive bacterium Bacillus subtilis. Furthermore, no effect of the plasma treatment, on the antibiotic resistance of E. coli, was determined under the chosen conditions. The knowledge about the physiological changes of ARM in response to plasma is of fundamental interest to understand the molecular basis for the inactivation. This will be important for the further development and implementation of plasma in wastewater remediation.

Marine particle microbiomes during a spring diatom bloom contain active sulfate-reducing bacteria

Phytoplankton blooms fuel marine food webs with labile dissolved carbon and also lead to the formation of particulate organic matter composed of living and dead algal cells. These particles contribute to carbon sequestration and are sites of intense algal-bacterial interactions, providing diverse niches for microbes to thrive. We analyzed 16S and 18S ribosomal RNA gene amplicon sequences obtained from 51 time points and metaproteomes from 3 time points during a spring phytoplankton bloom in a shallow location (6-10 m depth) in the North Sea. Particulate fractions larger than 10 µm diameter were collected at near daily intervals between early March and late May in 2018. Network analysis identified two major modules representing bacteria co-occurring with diatoms and with dinoflagellates, respectively. The diatom network module included known sulfate-reducing Desulfobacterota as well as potentially sulfur-oxidizing Ectothiorhodospiraceae. Metaproteome analyses confirmed presence of key enzymes involved in dissimilatory sulfate reduction, a process known to occur in sinking particles at greater depths and in sediments. Our results indicate the presence of sufficiently anoxic niches in the particle fraction of an active phytoplankton bloom to sustain sulfate reduction, and an important role of benthic-pelagic coupling for microbiomes in shallow environments. Our findings may have implications for the understanding of algal-bacterial interactions and carbon export during blooms in shallow-water coastal areas.

Annotating microbial functions with ProkFunFind

Analyzing microbial genomes has become an essential part of microbiology research, giving valuable insights into the functions and evolution of microbial species. Identifying genes of interest and assigning putative annotations to those genes is a central task in genome analysis, and a plethora of tools and approaches have been developed for this task. The ProkFunFind tool was developed to bridge the gap between these various annotation approaches, providing a flexible and customizable search approach to annotate microbial functions. ProkFunFind is designed around hierarchical definitions of biological functions, where individual genes can be identified using heterogeneous search terms consisting of sequences, profile hidden Markov models, protein domains, and orthology groups. This flexible and customizable search approach allows for searches to be tailored to specific biological functions, and the search results are output in multiple formats to facilitate downstream analyses. The utility of the ProkFunFind search tool was demonstrated through its application in searching for bacterial flagella, which are complex organelles composed of multiple genes. Overall, ProkFunFind provides an accessible and flexible way to integrate multiple types of annotation and sequence data while annotating biological functions in microbial genomes.IMPORTANCEGenome sequencing and analysis are increasingly important parts of microbiology, providing a way to predict metabolic functions, identify virulence factors, and understand the evolution of microbes. The expanded use of genome sequencing has also brought an abundance of search and annotation methods, but integrating the information from these different methods can be challenging and is often done through ad hoc approaches. To bridge the gap between different types of annotations, we developed ProkFunFind, a flexible and customizable search tool incorporating multiple search approaches and annotation types to annotate microbial functions. We demonstrated the utility of ProkFunFind by searching for gene clusters encoding flagellar genes using a combination of different annotation types and searches. Overall, ProkFunFind provides a reproducible and flexible way to identify gene clusters of interest, facilitating the meaningful analysis of new and existing microbial genomes.

Genome assembly of Medicago truncatula accession SA27063 provides insight into spring black stem and leaf spot disease resistance

Medicago truncatula, model legume and alfalfa relative, has served as an essential resource for advancing our understanding of legume physiology, functional genetics, and crop improvement traits. Necrotrophic fungus, Ascochyta medicaginicola, the causal agent of spring black stem (SBS) and leaf spot is a devasting foliar disease of alfalfa affecting stand survival, yield, and forage quality. Host resistance to SBS disease is poorly understood, and control methods rely on cultural practices. Resistance has been observed in M. truncatula accession SA27063 (HM078) with two recessively inherited quantitative-trait loci (QTL), rnpm1 and rnpm2, previously reported. To shed light on host resistance, we carried out a de novo genome assembly of HM078. The genome, referred to as MtHM078 v1.0, is comprised of 23 contigs totaling 481.19 Mbp. Notably, this assembly contains a substantial amount of novel centromere-related repeat sequences due to deep long-read sequencing. Genome annotation resulted in 98.4% of BUSCO fabales proteins being complete. The assembly enabled sequence-level analysis of rnpm1 and rnpm2 for gene content, synteny, and structural variation between SBS-resistant accession SA27063 (HM078) and SBS-susceptible accession A17 (HM101). Fourteen candidate genes were identified, and some have been implicated in resistance to necrotrophic fungi. Especially interesting candidates include loss-of-function events in HM078 because they fit the inverse gene-for-gene model, where resistance is recessively inherited. In rnpm1, these include a loss-of-function in a disease resistance gene due to a premature stop codon, and a 10.85 kbp retrotransposon-like insertion disrupting a ubiquitin conjugating E2. In rnpm2, we identified a frameshift mutation causing a loss-of-function in a glycosidase, as well as a missense and frameshift mutation altering an F-box family protein. This study generated a high-quality genome of HM078 and has identified promising candidates, that once validated, could be further studied in alfalfa to enhance disease resistance.

Pontiella agarivorans sp. nov., a novel marine anaerobic bacterium capable of degrading macroalgal polysaccharides and fixing nitrogen

Marine macroalgae produce abundant and diverse polysaccharides, which contribute substantially to the organic matter exported to the deep ocean. Microbial degradation of these polysaccharides plays an important role in the turnover of macroalgal biomass. Various members of the Planctomycetes-Verrucomicrobia-Chlamydia (PVC) superphylum are degraders of polysaccharides in widespread anoxic environments. In this study, we isolated a novel anaerobic bacterial strain NLcol2T from microbial mats on the surface of marine sediments offshore Santa Barbara, CA, USA. Based on 16S ribosomal RNA (rRNA) gene and phylogenomic analyses, strain NLcol2T represents a novel species within the Pontiella genus in the Kiritimatiellota phylum (within the PVC superphylum). Strain NLcol2T is able to utilize various monosaccharides, disaccharides, and macroalgal polysaccharides such as agar and ɩ-carrageenan. A near-complete genome also revealed an extensive metabolic capacity for anaerobic degradation of sulfated polysaccharides, as evidenced by 202 carbohydrate-active enzymes (CAZymes) and 165 sulfatases. Additionally, its ability of nitrogen fixation was confirmed by nitrogenase activity detected during growth on nitrogen-free medium, and the presence of nitrogenases (nifDKH) encoded in the genome. Based on the physiological and genomic analyses, this strain represents a new species of bacteria that may play an important role in the degradation of macroalgal polysaccharides and with relevance to the biogeochemical cycling of carbon, sulfur, and nitrogen in marine environments. Strain NLcol2T (= DSM 113125T = MCCC 1K08672T) is proposed to be the type strain of a novel species in the Pontiella genus, and the name Pontiella agarivorans sp. nov. is proposed.IMPORTANCEGrowth and intentional burial of marine macroalgae is being considered as a carbon dioxide reduction strategy but elicits concerns as to the fate and impacts of this macroalgal carbon in the ocean. Diverse heterotrophic microbial communities in the ocean specialize in these complex polymers such as carrageenan and fucoidan, for example, members of the Kiritimatiellota phylum. However, only four type strains within the phylum have been cultivated and characterized to date, and there is limited knowledge about the metabolic capabilities and functional roles of related organisms in the environment. The new isolate strain NLcol2T expands the known substrate range of this phylum and further reveals the ability to fix nitrogen during anaerobic growth on macroalgal polysaccharides, thereby informing the issue of macroalgal carbon disposal.

Structural and Functional Disparities within the Human Gut Virome in Terms of Genome Topology and Representative Genome Selection

Circularity confers protection to viral genomes where linearity falls short, thereby fulfilling the form follows function aphorism. However, a shift away from morphology-based classification toward the molecular and ecological classification of viruses is currently underway within the field of virology. Recent years have seen drastic changes in the International Committee on Taxonomy of Viruses' operational definitions of viruses, particularly for the tailed phages that inhabit the human gut. After the abolition of the order Caudovirales, these tailed phages are best defined as members of the class Caudoviricetes. To determine the epistemological value of genome topology in the context of the human gut virome, we designed a set of seven experiments to assay the impact of genome topology and representative viral selection on biological interpretation. Using Oxford Nanopore long reads for viral genome assembly coupled with Illumina short-read polishing, we showed that circular and linear virus genomes differ remarkably in terms of genome quality, GC skew, transfer RNA gene frequency, structural variant frequency, cross-reference functional annotation (COG, KEGG, Pfam, and TIGRfam), state-of-the-art marker-based classification, and phage-host interaction. Furthermore, the disparity profile changes during dereplication. In particular, our phage-host interaction results demonstrated that proportional abundances cannot be meaningfully compared without due regard for genome topology and dereplication threshold, which necessitates the need for standardized reporting. As a best practice guideline, we recommend that comparative studies of the human gut virome always report the ratio of circular to linear viral genomes along with the dereplication threshold so that structural and functional metrics can be placed into context when assessing biologically relevant metagenomic properties such as proportional abundance.

Anaerobic hexadecane degradation by a thermophilic Hadarchaeon from Guaymas Basin

Hadarchaeota inhabit subsurface and hydrothermally heated environments, but previous to this study, they had not been cultured. Based on metagenome-assembled genomes, most Hadarchaeota are heterotrophs that grow on sugars and amino acids, or oxidize carbon monoxide or reduce nitrite to ammonium. A few other metagenome-assembled genomes encode alkyl-coenzyme M reductases (Acrs), β-oxidation, and Wood-Ljungdahl pathways, pointing toward multicarbon alkane metabolism. To identify the organisms involved in thermophilic oil degradation, we established anaerobic sulfate-reducing hexadecane-degrading cultures from hydrothermally heated sediments of the Guaymas Basin. Cultures at 70°C were enriched in one Hadarchaeon that we propose as Candidatus Cerberiarchaeum oleivorans. Genomic and chemical analyses indicate that Ca. C. oleivorans uses an Acr to activate hexadecane to hexadecyl-coenzyme M. A β-oxidation pathway and a tetrahydromethanopterin methyl branch Wood-Ljungdahl (mWL) pathway allow the complete oxidation of hexadecane to CO2. Our results suggest a syntrophic lifestyle with sulfate reducers, as Ca. C. oleivorans lacks a sulfate respiration pathway. Comparative genomics show that Acr, mWL, and β-oxidation are restricted to one family of Hadarchaeota, which we propose as Ca. Cerberiarchaeaceae. Phylogenetic analyses further indicate that the mWL pathway is basal to all Hadarchaeota. By contrast, the carbon monoxide dehydrogenase/acetyl-coenzyme A synthase complex in Ca. Cerberiarchaeaceae was horizontally acquired from Bathyarchaeia. The Acr and β-oxidation genes of Ca. Cerberiarchaeaceae are highly similar to those of other alkane-oxidizing archaea such as Ca. Methanoliparia and Ca. Helarchaeales. Our results support the use of Acrs in the degradation of petroleum alkanes and suggest a role of Hadarchaeota in oil-rich environments.

Flexible genomic island conservation across freshwater and marine Methylophilaceae

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Unraveling the maternal and paternal origins of allotetraploid Vigna reflexo-pilosa

The genomic structures of Vigna hirtella Ridl. and Vigna trinervia (B.Heyne ex Wight & Arn.) Tateishi & Maxted, key ancestral species of the allotetraploid Vigna reflexo-pilosa var. glabra (Roxb.) N.Tomooka & Maxted, remain poorly understood. This study presents a comprehensive genomic comparison of these species to deepen our knowledge of their evolutionary trajectories. By comparing the genomic profiles of V. hirtella and V. trinervia with those of V. reflexo-pilosa, we investigate the complex genomic mechanisms underlying allopolyploid evolution within the genus Vigna. Comparison of the chloroplast genome revealed that V. trinervia is closely related to V. reflexo-pilosa. De novo assembly of the whole genome, followed by synteny analysis and Ks value calculations, confirms that V. trinervia is closely related to the A genome of V. reflexo-pilosa, and V. hirtella to its B genome. Furthermore, the comparative analyses reveal that V. reflexo-pilosa retains residual signatures of a previous polyploidization event, particularly evident in higher gene family copy numbers. Our research provides genomic evidence for polyploidization within the genus Vigna and identifies potential donor species of allotetraploid species using de novo assembly techniques. Given the Southeast Asian distribution of both V. hirtella and V. trinervia, natural hybridization between these species, with V. trinervia as the maternal ancestor and V. hirtella as the paternal donor, seems plausible.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59

IMPORTANCE

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

Benzylic Radical Stabilization Permits Ether Formation During Darobactin Biosynthesis

The Gram-negative selective antibiotic darobactin A has attracted interest owing to its intriguing fused bicyclic structure and unique mode of action. Biosynthetic studies have revealed that darobactin is a ribosomally synthesized and post-translationally modified peptide (RiPP). During maturation, the darobactin precursor peptide (DarA) is modified by a radical S-adenosyl methionine (rSAM)-dependent enzyme (DarE) to contain ether and C-C crosslinks. In this work, we describe the enzymatic tolerance of DarE using a panel of DarA variants, revealing that DarE can install the ether and C-C crosslinks independently and in different locations on DarA. These efforts produced 57 darobactin variants, 50 of which were enzymatically modified. Several new variants with fused bicyclic structures were characterized, including darobactin W3Y, which replaces tryptophan with tyrosine at the twice-modified central position, and darobactin K5F, which displays a fused diether ring pattern. Three additional darobactin variants contained fused diether macrocycles, leading us to investigate the origin of ether versus C-C crosslink formation. Computational analyses found that more stable and long-lived Cβ radicals found on aromatic amino acids correlated with ether formation. Further, molecular docking and calculated transition state structures provide support for the different indole connectivity observed for ether (Trp-C7) and C-C (Trp-C6) crosslink formation. We also provide experimental evidence for a β-oxotryptophan modification, a proposed intermediate during ether crosslink formation. Finally, mutational analysis of the DarA leader region and protein structural predictions identified which residues were dispensable for processing and others that govern substrate engagement by DarE. Our work informs on darobactin scaffold engineering and sheds additional light on the underlying principles of rSAM catalysis.

Bacteria can maintain rRNA operons solely on plasmids for hundreds of millions of years

It is generally assumed that all bacteria must have at least one rRNA operon (rrn operon) on the chromosome, but some strains of the genera Aureimonas and Oecophyllibacter carry their sole rrn operon on a plasmid. However, other related strains and species have chromosomal rrn loci, suggesting that the exclusive presence of rrn operons on a plasmid is rare and unlikely to be stably maintained over long evolutionary periods. Here, we report the results of a systematic search for additional bacteria without chromosomal rrn operons. We find that at least four bacterial clades in the phyla Bacteroidota, Spirochaetota, and Pseudomonadota (Proteobacteria) lost chromosomal rrn operons independently. Remarkably, Persicobacteraceae have apparently maintained this peculiar genome organization for hundreds of millions of years. In our study, all the rrn-carrying plasmids in bacteria lacking chromosomal rrn loci possess replication initiator genes of the Rep_3 family. Furthermore, the lack of chromosomal rrn operons is associated with differences in copy numbers of rrn operons, plasmids, and chromosomal tRNA genes. Thus, our findings indicate that the absence of rrn loci in bacterial chromosomes can be stably maintained over long evolutionary periods.

Role of Relebactam in the Antibiotic Resistance Acquisition in Pseudomonas aeruginosa: In Vitro Study

BACKGROUND

Pseudomonas aeruginosa shows resistance to several antibiotics and often develops such resistance during patient treatment.

OBJECTIVE

Develop an in vitro model, using clinical isolates of P. aeruginosa, to compare the ability of the imipenem and imipenem/relebactam to generate resistant mutants to imipenem and to other antibiotics. Perform a genotypic analysis to detect how the selective pressure changes their genomes.

METHODS

The antibiotics resistance was studied by microdilution assays and e-test, and the genotypic study was performed by NGS.

RESULTS

The isolates acquired resistance to imipenem in an average of 6 days, and to imipenem/relebactam in 12 days (p value = 0.004). After 30 days of exposure, 75% of the isolates reached a MIC > 64 mg/L for imipenem and 37.5% for imipenem/relebactam (p value = 0.077). The 37.5% and the 12.5% imipenem/relebactam mutants developed resistance to piperacillin/tazobactam and ceftazidime, respectively, while the 87.5% and 37.5% of the imipenem mutants showed resistance to these drugs (p value = 0.003, p value = 0.015). The main biological processes altered by the SNPs were the glycosylation pathway, transcriptional regulation, histidine kinase response, porins, and efflux pumps.

DISCUSSION

The addition of relebactam delays the generation of resistance to imipenem and limits the cross-resistance to other beta-lactams. The clinical relevance of this phenomenon, which has the limitation that it has been performed in vitro, should be evaluated by stewardship programs in clinical practice, as it could be useful in controlling multi-drug resistance in P. aeruginosa.

Multi-heme cytochrome-mediated extracellular electron transfer by the anaerobic methanotroph 'Candidatus Methanoperedens nitroreducens'

Anaerobic methanotrophic archaea (ANME) carry out anaerobic oxidation of methane, thus playing a crucial role in the methane cycle. Previous genomic evidence indicates that multi-heme c-type cytochromes (MHCs) may facilitate the extracellular electron transfer (EET) from ANME to different electron sinks. Here, we provide experimental evidence supporting cytochrome-mediated EET for the reduction of metals and electrodes by 'Candidatus Methanoperedens nitroreducens', an ANME acclimated to nitrate reduction. Ferrous iron-targeted fluorescent assays, metatranscriptomics, and single-cell imaging suggest that 'Ca. M. nitroreducens' uses surface-localized redox-active cytochromes for metal reduction. Electrochemical and Raman spectroscopic analyses also support the involvement of c-type cytochrome-mediated EET for electrode reduction. Furthermore, several genes encoding menaquinone cytochrome type-c oxidoreductases and extracellular MHCs are differentially expressed when different electron acceptors are used.

Investigating in vivo Mycobacterium avium subsp. paratuberculosis microevolution and mixed strain infections

Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne's Disease (JD) in ruminants, which is responsible for significant economic loss to the global dairy industry. Mixed strain infection (MSI) refers to the concurrent infection of a susceptible host with genetically distinct strains of a pathogen, whereas within-host changes in an infecting strain leading to genetically distinguishable progeny is called microevolution. The two processes can influence host-pathogen dynamics, disease progression and outcomes, but not much is known about their prevalence and impact on JD. Therefore, we obtained up to 10 MAP isolates each from 14 high-shedding animals and subjected them to whole-genome sequencing. Twelve of the 14 animals examined showed evidence for the presence of MSIs and microevolution, while the genotypes of MAP isolates from the remaining two animals could be attributed solely to microevolution. All MAP isolates that were otherwise isogenic had differences in short sequence repeats (SSRs), of which SSR1 and SSR2 were the most diverse and homoplastic. Variations in SSR1 and SSR2, which are located in ORF1 and ORF2, respectively, affect the genetic reading frame, leading to protein products with altered sequences and computed structures. The ORF1 gene product is predicted to be a MAP surface protein with possible roles in host immune modulation, but nothing could be inferred regarding the function of ORF2. Both genes are conserved in Mycobacterium avium complex members, but SSR1-based modulation of ORF1 reading frames seems to only occur in MAP, which could have potential implications on the infectivity of this pathogen. IMPORTANCE Johne's disease (JD) is a major problem in dairy animals, and concerns have been raised regarding the association of Mycobacterium avium subsp. paratuberculosis (MAP) with Crohn's disease in humans. MAP is an extremely slow-growing bacterium with low genome evolutionary rates. Certain short sequence repeats (SSR1 and SSR2) in the MAP chromosome are highly variable and evolve at a faster rate than the rest of the chromosome. In the current study, multiple MAP isolates with genetic variations such as single-nucleotide polymorphisms, and more noticeably, diverse SSRs, could simultaneously infect animals. Variations in SSR1 and SSR2 affect the products of the respective genes containing them. Since multiple MAP isolates can infect the same animal and the possibility that the pathogen undergoes further changes within the host due to unstable SSRs, this could provide a compensative mechanism for an otherwise slow-evolving pathogen to increase phenotypic diversity for overcoming host responses.

Toxic antiphage defense proteins inhibited by intragenic antitoxin proteins

Recombination-promoting nuclease (Rpn) proteins are broadly distributed across bacterial phyla, yet their functions remain unclear. Here, we report that these proteins are toxin-antitoxin systems, comprised of genes-within-genes, that combat phage infection. We show the small, highly variable Rpn C-terminal domains (RpnS), which are translated separately from the full-length proteins (RpnL), directly block the activities of the toxic RpnL. The crystal structure of RpnAS revealed a dimerization interface encompassing α helix that can have four amino acid repeats whose number varies widely among strains of the same species. Consistent with strong selection for the variation, we document that plasmid-encoded RpnP2L protects Escherichia coli against certain phages. We propose that many more intragenic-encoded proteins that serve regulatory roles remain to be discovered in all organisms.

Isolation, biochemical characterization, and genome sequencing of two high-quality genomes of a novel chitinolytic Jeongeupia species

Chitin is the second most abundant polysaccharide worldwide as part of arthropods' exoskeletons and fungal cell walls. Low concentrations in soils and sediments indicate rapid decomposition through chitinolytic organisms in terrestrial and aquatic ecosystems. The enacting enzymes, so-called chitinases, and their products, chitooligosaccharides, exhibit promising characteristics with applications ranging from crop protection to cosmetics, medical, textile, and wastewater industries. Exploring novel chitinolytic organisms is crucial to expand the enzymatical toolkit for biotechnological chitin utilization and to deepen our understanding of diverse catalytic mechanisms. In this study, we present two long-read sequencing-based genomes of highly similar Jeongeupia species, which have been screened, isolated, and biochemically characterized from chitin-amended soil samples. Through metabolic characterization, whole-genome alignments, and phylogenetic analysis, we could demonstrate how the investigated strains differ from the taxonomically closest strain Jeongeupia naejangsanensis BIO-TAS4-2T (DSM 24253). In silico analysis and sequence alignment revealed a multitude of highly conserved chitinolytic enzymes in the investigated Jeongeupia genomes. Based on these results, we suggest that the two strains represent a novel species within the genus of Jeongeupia, which may be useful for environmentally friendly N-acetylglucosamine production from crustacean shell or fungal biomass waste or as a crop protection agent.

Small protein modules dictate prophage fates during polylysogeny

Most bacteria in the biosphere are predicted to be polylysogens harbouring multiple prophages1-5. In studied systems, prophage induction from lysogeny to lysis is near-universally driven by DNA-damaging agents6. Thus, how co-residing prophages compete for cell resources if they respond to an identical trigger is unknown. Here we discover regulatory modules that control prophage induction independently of the DNA-damage cue. The modules bear little resemblance at the sequence level but share a regulatory logic by having a transcription factor that activates the expression of a neighbouring gene that encodes a small protein. The small protein inactivates the master repressor of lysis, which leads to induction. Polylysogens that harbour two prophages exposed to DNA damage release mixed populations of phages. Single-cell analyses reveal that this blend is a consequence of discrete subsets of cells producing one, the other or both phages. By contrast, induction through the DNA-damage-independent module results in cells producing only the phage sensitive to that specific cue. Thus, in the polylysogens tested, the stimulus used to induce lysis determines phage productivity. Considering the lack of potent DNA-damaging agents in natural habitats, additional phage-encoded sensory pathways to lysis likely have fundamental roles in phage-host biology and inter-prophage competition.

VIRify: An integrated detection, annotation and taxonomic classification pipeline using virus-specific protein profile hidden Markov models

The study of viral communities has revealed the enormous diversity and impact these biological entities have on various ecosystems. These observations have sparked widespread interest in developing computational strategies that support the comprehensive characterisation of viral communities based on sequencing data. Here we introduce VIRify, a new computational pipeline designed to provide a user-friendly and accurate functional and taxonomic characterisation of viral communities. VIRify identifies viral contigs and prophages from metagenomic assemblies and annotates them using a collection of viral profile hidden Markov models (HMMs). These include our manually-curated profile HMMs, which serve as specific taxonomic markers for a wide range of prokaryotic and eukaryotic viral taxa and are thus used to reliably classify viral contigs. We tested VIRify on assemblies from two microbial mock communities, a large metagenomics study, and a collection of publicly available viral genomic sequences from the human gut. The results showed that VIRify could identify sequences from both prokaryotic and eukaryotic viruses, and provided taxonomic classifications from the genus to the family rank with an average accuracy of 86.6%. In addition, VIRify allowed the detection and taxonomic classification of a range of prokaryotic and eukaryotic viruses present in 243 marine metagenomic assemblies. Finally, the use of VIRify led to a large expansion in the number of taxonomically classified human gut viral sequences and the improvement of outdated and shallow taxonomic classifications. Overall, we demonstrate that VIRify is a novel and powerful resource that offers an enhanced capability to detect a broad range of viral contigs and taxonomically classify them.

In-silico functional analysis of hypothetical proteins from Lactiplantibacillus plantarum plasmids reveals enrichment of cell envelope proteins

Lactiplantibacillus plantarum is one of the important species of lactic acid bacterium (LAB) found in diverse environments, with many strains exhibiting probiotic properties. In our previous study, 41.6% of protein families (PFs) encoded by 395 plasmids from several L. plantarum strains were found to be hypothetical proteins with no predicted function. This study aimed at predicting the functions of these 647 hypothetical proteins using 21 different bioinformatics methods. As a result, 160 PFs could be newly annotated. A lower proportion of plasmid-specific functions was annotated as compared to the functions shared between plasmids and chromosomes. Also, hypothetical proteins were less conserved than the annotated proteins across L.plantarum plasmids. Based on the subcellular localization, cell envelope proteins represented the biggest category in the newly annotated proteins. Transporters (112 PFs) which was a part of cell envelop proteins represented the largest functional group. Additionally, 40 and 25 other PFs were predicted to contain signal peptides and transmembrane helices, respectively. We speculate that such hypothetical proteins might be involved in the transport of various chemicals and environmental interactions in L. plantarum. In the future, functional characterization of these proteins through wet-lab experimental approach can provide novel insights into their contribution to the physiology, probiotic properties, and industrial utility of these bacteria.

Proteomics analysis of urine and catheter-associated biofilms in spinal cord injury patients

After spinal cord injury (SCI), use chronic urinary catheters for bladder management is common, making these patients especially vulnerable to catheter-associated complications. Chronic catheterization is associated with bacterial colonization and frequent catheter-associated urinary tract infections (CAUTI). One determinant of infection success and treatment resistance is production of catheter-associated biofilms, composed of microorganisms and host- and microbial-derived components. To better understand the biofilm microenvironment, we performed proteomics analysis of catheter-associated biofilms and paired urine samples from four people with SCI with chronic indwelling urinary catheters. We developed a novel method for the removal of adhered cellular components on catheters that contained both human and microbial homologous proteins. Proteins from seven microbial species were identified including: Escherichia coli, Klebsiella species (spp), Enterococcus spp, Proteus mirabilis, Pseudomonas spp, Staphylococcus spp, and Candida spp. Peptides identified from catheter biofilms were assigned to 4,820 unique proteins, with 61% of proteins assigned to the biofilm-associated microorganisms, while the remainder were human-derived. Contrastingly, in urine, only 51% were assigned to biofilm-associated microorganisms and 4,554 proteins were identified as a human-derived. Of the proteins assigned to microorganisms in the biofilm and paired urine, Enterococcus, Candida spp, and P. mirabilis had greater associations with the biofilm phase, whereas E. coli and Klebsiella had greater associations with the urine phase, thus demonstrating a significant difference between the urine and adhered microbial communities. The microbial proteins that differed significantly between the biofilm and paired urine samples mapped to pathways associated with amino acid synthesis, likely related to adaptation to high urea concentrations in the urine, and growth and protein synthesis in bacteria in the biofilm. Human proteins demonstrated enrichment for immune response in the catheter-associated biofilm. Proteomic analysis of catheter-associated biofilms and paired urine samples has the potential to provide detailed information on host and bacterial responses to chronic indwelling urinary catheters and could be useful for understanding complications of chronic indwelling catheters including CAUTIs, urinary stones, and catheter blockages.

Genomic Assessment of the Contribution of the Wolbachia Endosymbiont of Eurosta solidaginis to Gall Induction

We explored the genome of the Wolbachia strain, wEsol, symbiotic with the plant-gall-inducing fly Eurosta solidaginis with the goal of determining if wEsol contributes to gall induction by its insect host. Gall induction by insects has been hypothesized to involve the secretion of the phytohormones cytokinin and auxin and/or proteinaceous effectors to stimulate cell division and growth in the host plant. We sequenced the metagenome of E. solidaginis and wEsol and assembled and annotated the genome of wEsol. The wEsol genome has an assembled length of 1.66 Mbp and contains 1878 protein-coding genes. The wEsol genome is replete with proteins encoded by mobile genetic elements and shows evidence of seven different prophages. We also detected evidence of multiple small insertions of wEsol genes into the genome of the host insect. Our characterization of the genome of wEsol indicates that it is compromised in the synthesis of dimethylallyl pyrophosphate (DMAPP) and S-adenosyl L-methionine (SAM), which are precursors required for the synthesis of cytokinins and methylthiolated cytokinins. wEsol is also incapable of synthesizing tryptophan, and its genome contains no enzymes in any of the known pathways for the synthesis of indole-3-acetic acid (IAA) from tryptophan. wEsol must steal DMAPP and L-methionine from its host and therefore is unlikely to provide cytokinin and auxin to its insect host for use in gall induction. Furthermore, in spite of its large repertoire of predicted Type IV secreted effector proteins, these effectors are more likely to contribute to the acquisition of nutrients and the manipulation of the host's cellular environment to contribute to growth and reproduction of wEsol than to aid E. solidaginis in manipulating its host plant. Combined with earlier work that shows that wEsol is absent from the salivary glands of E. solidaginis, our results suggest that wEsol does not contribute to gall induction by its host.

Genome of tetraploid sour cherry (Prunus cerasus L.) 'Montmorency' identifies three distinct ancestral Prunus genomes

Sour cherry (Prunus cerasus L.) is a valuable fruit crop in the Rosaceae family and a hybrid between progenitors closely related to extant Prunus fruticosa (ground cherry) and Prunus avium (sweet cherry). Here we report a chromosome-scale genome assembly for sour cherry cultivar Montmorency, the predominant cultivar grown in the USA. We also generated a draft assembly of P. fruticosa to use alongside a published P. avium sequence for syntelog-based subgenome assignments for 'Montmorency' and provide compelling evidence P. fruticosa is also an allotetraploid. Using hierarchal k-mer clustering and phylogenomics, we show 'Montmorency' is trigenomic, containing two distinct subgenomes inherited from a P. fruticosa-like ancestor (A and A') and two copies of the same subgenome inherited from a P. avium-like ancestor (BB). The genome composition of 'Montmorency' is AA'BB and little-to-no recombination has occurred between progenitor subgenomes (A/A' and B). In Prunus, two known classes of genes are important to breeding strategies: the self-incompatibility loci (S-alleles), which determine compatible crosses, successful fertilization, and fruit set, and the Dormancy Associated MADS-box genes (DAMs), which strongly affect dormancy transitions and flowering time. The S-alleles and DAMs in 'Montmorency' and P. fruticosa were manually annotated and support subgenome assignments. Lastly, the hybridization event 'Montmorency' is descended from was estimated to have occurred less than 1.61 million years ago, making sour cherry a relatively recent allotetraploid. The 'Montmorency' genome highlights the evolutionary complexity of the genus Prunus and will inform future breeding strategies for sour cherry, comparative genomics in the Rosaceae, and questions regarding neopolyploidy.

PHERI-Phage Host ExploRation Pipeline

Antibiotic resistance is becoming a common problem in medicine, food, and industry, with multidrug-resistant bacterial strains occurring in all regions. One of the possible future solutions is the use of bacteriophages. Phages are the most abundant form of life in the biosphere, so we can highly likely purify a specific phage against each target bacterium. The identification and consistent characterization of individual phages was a common form of phage work and included determining bacteriophages' host-specificity. With the advent of new modern sequencing methods, there was a problem with the detailed characterization of phages in the environment identified by metagenome analysis. The solution to this problem may be to use a bioinformatic approach in the form of prediction software capable of determining a bacterial host based on the phage whole-genome sequence. The result of our research is the machine learning algorithm-based tool called PHERI. PHERI predicts the suitable bacterial host genus for the purification of individual viruses from different samples. In addition, it can identify and highlight protein sequences that are important for host selection.

Toxic anti-phage defense proteins inhibited by intragenic antitoxin proteins

Recombination-promoting nuclease (Rpn) proteins are broadly distributed across bacterial phyla, yet their functions remain unclear. Here we report these proteins are new toxin-antitoxin systems, comprised of genes-within-genes, that combat phage infection. We show the small, highly variable Rpn C -terminal domains (Rpn _S ), which are translated separately from the full-length proteins (Rpn _L ), directly block the activities of the toxic full-length proteins. The crystal structure of RpnA _S revealed a dimerization interface encompassing a helix that can have four amino acid repeats whose number varies widely among strains of the same species. Consistent with strong selection for the variation, we document plasmid-encoded RpnP2 _L protects Escherichia coli against certain phages. We propose many more intragenic-encoded proteins that serve regulatory roles remain to be discovered in all organisms.

Significance

Here we document the function of small genes-within-genes, showing they encode antitoxin proteins that block the functions of the toxic DNA endonuclease proteins encoded by the longer rpn genes. Intriguingly, a sequence present in both long and short protein shows extensive variation in the number of four amino acid repeats. Consistent with a strong selection for the variation, we provide evidence that the Rpn proteins represent a phage defense system.

A hydrogenotrophic Sulfurimonas is globally abundant in deep-sea oxygen-saturated hydrothermal plumes

Members of the bacterial genus Sulfurimonas (phylum Campylobacterota) dominate microbial communities in marine redoxclines and are important for sulfur and nitrogen cycling. Here we used metagenomics and metabolic analyses to characterize a Sulfurimonas from the Gakkel Ridge in the Central Arctic Ocean and Southwest Indian Ridge, showing that this species is ubiquitous in non-buoyant hydrothermal plumes at Mid Ocean Ridges across the global ocean. One Sulfurimonas species, ^USulfurimonas pluma, was found to be globally abundant and active in cold (<0-4 °C), oxygen-saturated and hydrogen-rich hydrothermal plumes. Compared with other Sulfurimonas species, ^US. pluma has a reduced genome (>17%) and genomic signatures of an aerobic chemolithotrophic metabolism using hydrogen as an energy source, including acquisition of A2-type oxidase and loss of nitrate and nitrite reductases. The dominance and unique niche of ^US. pluma in hydrothermal plumes suggest an unappreciated biogeochemical role for Sulfurimonas in the deep ocean.

Dual domestications and origin of traits in grapevine evolution

We elucidate grapevine evolution and domestication histories with 3525 cultivated and wild accessions worldwide. In the Pleistocene, harsh climate drove the separation of wild grape ecotypes caused by continuous habitat fragmentation. Then, domestication occurred concurrently about 11,000 years ago in Western Asia and the Caucasus to yield table and wine grapevines. The Western Asia domesticates dispersed into Europe with early farmers, introgressed with ancient wild western ecotypes, and subsequently diversified along human migration trails into muscat and unique western wine grape ancestries by the late Neolithic. Analyses of domestication traits also reveal new insights into selection for berry palatability, hermaphroditism, muscat flavor, and berry skin color. These data demonstrate the role of the grapevines in the early inception of agriculture across Eurasia.

Docsubty: FLALipid extract derived from newly isolated Rhodotorula toruloides LAB-07 for cosmetic applications

Rhodotorula toruloides is a non-conventional yeast with a natural carotenoid pathway. In particular, R. toruloides is an oleaginous yeast that can accumulate lipids in high content, thereby gaining interest as a promising industrial host. In this study, we isolated and taxonomically identified a new R. toruloides LAB-07 strain. De novo genome assembly using PacBio and Illumina hybrid platforms yielded 27 contigs with a 20.78 Mb genome size. Subsequent genome annotation analysis based on RNA-seq predicted 5296 protein-coding genes, including the fatty acid production pathway. We compared lipid production under different media; it was highest in the yeast extract salt medium with glycerol as a carbon source. Polyunsaturated α-linolenic acid was detected among the fatty acids, and docking phosphatidylcholine as a substrate to modeled Fad2, which annotated as Δ12-fatty acid desaturase showed bifunctional Δ12, 15-desaturation is structurally possible in that the distances between the diiron center and the carbon-carbon bond in which desaturation occurs were similar to those of structurally identified mouse stearoyl-CoA desaturase. Finally, the applicability of the extracted total lipid fraction of R. toruloides was investigated, demonstrating an increase in filaggrin expression and suppression of heat-induced MMP-1 expression when applied to keratinocytes, along with the additional antioxidant activity. This work presents a new R. toruloides LAB-07 strain with genomic and lipidomic data, which would help understand the physiology of R. toruloides. Also, the various skin-related effect of R. toruloides lipid extract indicates its potential usage as a promising cosmetic ingredient.

A novel and diverse group of Candidatus Patescibacteria from bathypelagic Lake Baikal revealed through long-read metagenomics

BACKGROUND

Lake Baikal, the world's deepest freshwater lake, contains important numbers of Candidatus Patescibacteria (formerly CPR) in its deepest reaches. However, previously obtained CPR metagenome-assembled genomes recruited very poorly indicating the potential of other groups being present. Here, we have applied for the first time a long-read (PacBio CCS) metagenomic approach to analyze in depth the Ca. Patescibacteria living in the bathypelagic water column of Lake Baikal at 1600 m.

RESULTS

The retrieval of nearly complete 16S rRNA genes before assembly has allowed us to detect the presence of a novel and a likely endemic group of Ca. Patescibacteria inhabiting bathypelagic Lake Baikal. This novel group seems to possess extremely high intra-clade diversity, precluding complete genomes' assembly. However, read binning and scaffolding indicate that these microbes are similar to other Ca. Patescibacteria (i.e. parasites or symbionts), although they seem to carry more anabolic pathways, likely reflecting the extremely oligotrophic habitat they inhabit. The novel bins have not been found anywhere, but one of the groups appears in small amounts in an oligotrophic and deep alpine Lake Thun. We propose this novel group be named Baikalibacteria.

CONCLUSION

The recovery of 16S rRNA genes via long-read metagenomics plus the use of long-read binning to uncover highly diverse "hidden" groups of prokaryotes are key strategies to move forward in ecogenomic microbiology. The novel group possesses enormous intraclade diversity akin to what happens with Ca. Patescibacteria at the interclade level, which is remarkable in an environment that has changed little in the last 25 million years.

Single-amplified genomes reveal most streamlined free-living marine bacteria

Evolutionary adaptations of prokaryotes to the environment sometimes result in genome reduction. Our knowledge of this phenomenon among free-living bacteria remains scarce. We address the dynamics and limits of genome reduction by examining one of the most abundant bacteria in the ocean, the SAR86 clade. Despite its abundance, comparative genomics has been limited by the absence of pure cultures and the poor representation in metagenome-assembled genomes. We co-assembled multiple previously available single-amplified genomes to obtain the first complete genomes from members of the four families. All families showed a convergent evolutionary trajectory with characteristic features of streamlined genomes, most pronounced in the TMED112 family. This family has a genome size of ca. 1 Mb and only 1 bp as median intergenic distance, exceeding values found in other abundant microbes such as SAR11, OM43 and Prochlorococcus. This genomic simplification led to a reduction in the biosynthesis of essential molecules, DNA repair-related genes, and the ability to sense and respond to environmental factors, which could suggest an evolutionary dependence on other co-occurring microbes for survival (Black Queen hypothesis). Therefore, these reconstructed genomes within the SAR86 clade provide new insights into the limits of genome reduction in free-living marine bacteria.

Genomic epidemiology of Mycobacterium avium subsp. paratuberculosis isolates from Canadian dairy herds provides evidence for multiple infection events

Mycobacterium avium subsp. paratuberculosis (MAP) is the pathogen responsible for paratuberculosis or Johne's Disease (JD) in ruminants, which is responsible for substantial economic losses worldwide. MAP transmission primarily occurs through the fecal-oral route, and the introduction of an MAP infected animal into a herd is an important transmission route. In the current study, we characterized MAP isolates from 67 cows identified in 20 herds from the provinces of Quebec and Ontario, Canada. Whole genome sequencing (WGS) was performed and an average genome coverage (relative to K-10) of ∼14.9 fold was achieved. The total number of SNPs present in each isolate varied from 51 to 132 and differed significantly between herds. Isolates with the highest genetic variability were generally present in herds from Quebec. The isolates were broadly separated into two main clades and this distinction was not influenced by the province from which they originated. Analysis of 8 MIRU-VNTR loci and 11 SSR loci was performed on the 67 isolates from the 20 dairy herds and publicly available references, notably major genetic lineages and six isolates from the province of Newfoundland and Labrador. All 67 field isolates were phylogenetically classified as Type II (C-type) and according to MIRU-VNTR, the predominant type was INMV 2 (76.1%) among four distinct patterns. Multilocus SSR typing identified 49 distinct INMV SSR patterns. The discriminatory index of the multilocus SSR typing was 0.9846, which was much higher than MIRU-VNTR typing (0.3740). Although multilocus SSR analysis provides good discriminatory power, the resolution was not informative enough to determine inter-herd transmission. In select cases, SNP-based analysis was the only approach able to document disease transmission between herds, further validated by animal movement data. The presence of SNPs in several virulence genes, notably for PE, PPE, mce and mmpL, is expected to explain differential antigenic or pathogenetic host responses. SNP-based studies will provide insight into how MAP genetic variation may impact host-pathogen interactions. Our study highlights the informative power of WGS which is now recommended for epidemiological studies and to document mixed genotypes infections.

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

The gut microbiota in multiple sclerosis varies with disease activity

BACKGROUND

Multiple sclerosis is a chronic immune-mediated disease of the brain and spinal cord resulting in physical and cognitive impairment in young adults. It is hypothesized that a disrupted bacterial and viral gut microbiota is a part of the pathogenesis mediating disease impact through an altered gut microbiota-brain axis. The aim of this study is to explore the characteristics of gut microbiota in multiple sclerosis and to associate it with disease variables, as the etiology of the disease remains only partially known.

METHODS

Here, in a case-control setting involving 148 Danish cases with multiple sclerosis and 148 matched healthy control subjects, we performed shotgun sequencing of fecal microbial DNA and associated bacterial and viral microbiota findings with plasma cytokines, blood cell gene expression profiles, and disease activity.

RESULTS

We found 61 bacterial species that were differentially abundant when comparing all multiple sclerosis cases with healthy controls, among which 31 species were enriched in cases. A cluster of inflammation markers composed of blood leukocytes, CRP, and blood cell gene expression of IL17A and IL6 was positively associated with a cluster of multiple sclerosis-related species. Bacterial species that were more abundant in cases with disease-active treatment-naïve multiple sclerosis were positively linked to a group of plasma cytokines including IL-22, IL-17A, IFN-β, IL-33, and TNF-α. The bacterial species richness of treatment-naïve multiple sclerosis cases was associated with number of relapses over a follow-up period of 2 years. However, in non-disease-active cases, we identified two bacterial species, Faecalibacterium prausnitzii and Gordonibacter urolithinfaciens, whose absolute abundance was enriched. These bacteria are known to produce anti-inflammatory metabolites including butyrate and urolithin. In addition, cases with multiple sclerosis had a higher viral species diversity and a higher abundance of Caudovirales bacteriophages.

CONCLUSIONS

Considerable aberrations are present in the gut microbiota of patients with multiple sclerosis that are directly associated with blood biomarkers of inflammation, and in treatment-naïve cases bacterial richness is positively associated with disease activity. Yet, the finding of two symbiotic bacterial species in non-disease-active cases that produce favorable immune-modulating compounds provides a rationale for testing these bacteria as adjunct therapeutics in future clinical trials.

Both Disease Activity and HLA-B27 Status Are Associated With Gut Microbiome Dysbiosis in Spondyloarthritis Patients

OBJECTIVE

Gut microbiome dysbiosis has previously been reported in spondyloarthritis (SpA) patients and could be critically involved in the pathogenesis of this disorder. The objectives of this study were to further characterize the microbiota structure in SpA patients and to investigate the relationship between dysbiosis and disease activity in light of the putative influence of the genetic background.

METHODS

Shotgun sequencing was performed on fecal DNA isolated from stool samples from 2 groups of adult volunteers: SpA patients (n = 102) and healthy controls (n = 63). A subset of the healthy controls comprised the age-matched siblings of patients whose HLA-B27 status was known. Changes in gut microbiota composition were assessed based on species diversity, enterotypes, and taxonomic and functional differences.

RESULTS

Dysbiosis was confirmed in SpA patients as compared to healthy controls. The restriction of microbiota diversity was detected in patients with the most active disease, and the abundance of several bacterial species was correlated with Bath Ankylosing Spondylitis Disease Activity Index score. Among healthy controls, significant differences in microbiota composition were also detected between the HLA-B27-positive and the HLA-B27-negative siblings of SpA patients. We highlighted a decreased abundance of several species of bacteria in SpA patients, especially those bacteria belonging to the Clostridiales order. Among the few species of bacteria showing increased abundance, Ruminococcus gnavus was one of the top differentiating species.

CONCLUSION

These findings reveal that genetic background and level of disease activity are likely to influence the composition of the gut microbiota of patients with SpA. It may be appropriate for further research on chronic arthritis to focus on these key parameters.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

Metagenomic Features Characterized with Microbial Iron Oxidoreduction and Mineral Interaction in Southwest Indian Ridge

The Southwest Indian Ridge (SWIR) is one of the typical representatives of deep-sea ultraslow-spreading ridges, and has increasingly become a hot spot of studying subsurface geological activities and deep-sea mining management. However, the understanding of microbial activities is still limited on active hydrothermal vent chimneys in SWIR. In this study, samples from an active black smoker and a diffuse vent located in the Longqi hydrothermal region were collected for deep metagenomic sequencing, which yielded approximately 290 GB clean data and 295 mid-to-high-quality metagenome-assembled genomes (MAGs). Sulfur oxidation conducted by a variety of Gammaproteobacteria, Alphaproteobacteria, and Campylobacterota was presumed to be the major energy source for chemosynthesis in Longqi hydrothermal vents. Diverse iron-related microorganisms were recovered, including iron-oxidizing Zetaproteobacteria, iron-reducing Deferrisoma, and magnetotactic bacterium. Twenty-two bacterial MAGs from 12 uncultured phyla harbored iron oxidase Cyc2 homologs and enzymes for organic carbon degradation, indicated novel chemolithoheterotrophic iron-oxidizing bacteria that affected iron biogeochemistry in hydrothermal vents. Meanwhile, potential interactions between microbial communities and chimney minerals were emphasized as enriched metabolic potential of siderophore transportation, and extracellular electron transfer functioned by multi-heme proteins was discovered. Composition of chimney minerals probably affected microbial iron metabolic potential, as pyrrhotite might provide more available iron for microbial communities. Collectively, this study provides novel insights into microbial activities and potential mineral-microorganism interactions in hydrothermal vents. IMPORTANCE Microbial activities and interactions with minerals and venting fluid in active hydrothermal vents remain unclear in the ultraslow-spreading SWIR (Southwest Indian Ridge). Understanding about how minerals influence microbial metabolism is currently limited given the obstacles in cultivating microorganisms with sulfur or iron oxidoreduction functions. Here, comprehensive descriptions on microbial composition and metabolic profile on 2 hydrothermal vents in SWIR were obtained based on cultivation-free metagenome sequencing. In particular, autotrophic sulfur oxidation supported by minerals was presumed, emphasizing the role of chimney minerals in supporting chemosynthesis. Presence of novel heterotrophic iron-oxidizing bacteria was also indicated, suggesting overlooked biogeochemical pathways directed by microorganisms that connected sulfide mineral dissolution and organic carbon degradation in hydrothermal vents. Our findings offer novel insights into microbial function and biotic interactions on minerals in ultraslow-spreading ridges.

Genome-Centric Dynamics Shape the Diversity of Oral Bacterial Populations

Two major viewpoints have been put forward for how microbial populations change, differing in whether adaptation is driven principally by gene-centric or genome-centric processes. Longitudinal sampling at microbially relevant timescales, i.e., days to weeks, is critical for distinguishing these mechanisms. Because of its significance for both microbial ecology and human health and its accessibility and high level of curation, we used the oral microbiota to study bacterial intrapopulation genome dynamics. Metagenomes were generated by shotgun sequencing of total community DNA from the healthy tongues of 17 volunteers at four to seven time points obtained over intervals of days to weeks. We obtained 390 high-quality metagenome-assembled genomes (MAGs) defining population genomes from 55 genera. The vast majority of genes in each MAG were tightly linked over the 2-week sampling window, indicating that the majority of the population's genomes were temporally stable at the MAG level. MAG-defined populations were composed of up to 5 strains, as determined by single-nucleotide-variant frequencies. Although most were stable over time, individual strains carrying over 100 distinct genes that rose from low abundance to dominance in a population over a period of days were detected. These results indicate a genome-wide as opposed to a gene-level process of population change. We infer that genome-wide selection of ecotypes is the dominant mode of adaptation in the oral populations over short timescales. IMPORTANCE The oral microbiome represents a microbial community of critical relevance to human health. Recent studies have documented the diversity and dynamics of different bacteria to reveal a rich, stable ecosystem characterized by strain-level dynamics. However, bacterial populations and their genomes are neither monolithic nor static; their genomes are constantly evolving to lose, gain, or alter their functional potential. To better understand how microbial genomes change in complex communities, we used culture-independent approaches to reconstruct the genomes (MAGs) for bacterial populations that approximated different species, in 17 healthy donors' mouths over a 2-week window. Our results underscored the importance of strain-level dynamics, which agrees with and expands on the conclusions of previous research. Altogether, these observations reveal patterns of genomic dynamics among strains of oral bacteria occurring over a matter of days.

Bioinformatic prediction and experimental validation of RiPP recognition elements

Ribosomally synthesized and post-translationally modified peptides (RiPPs) are a family of natural products for which discovery efforts have rapidly grown over the past decade. There are currently 38 known RiPP classes encoded by prokaryotes. Half of the prokaryotic RiPP classes include a protein domain called the RiPP Recognition Element (RRE) for successful installation of post-translational modifications on a RiPP precursor peptide. In most cases, the RRE domain binds to the N-terminal "leader" region of the precursor peptide, facilitating enzymatic modification of the C-terminal "core" region. The prevalence of the RRE domain renders it a theoretically useful bioinformatic handle for class-independent RiPP discovery; however, first-in-class RiPPs have yet to be isolated and experimentally characterized using an RRE-centric strategy. Moreover, with most known RRE domains engaging their cognate precursor peptide(s) with high specificity and nanomolar affinity, evaluation of the residue-specific interactions that govern RRE:substrate complexation is a necessary first step to leveraging the RRE domain for various bioengineering applications. This chapter details protocols for developing custom bioinformatic models to predict and annotate RRE domains in a class-specific manner. Next, we outline methods for experimental validation of precursor peptide binding using fluorescence polarization binding assays and in vitro enzyme activity assays. We anticipate the methods herein will guide and enhance future critical analyses of the RRE domain, eventually enabling its future use as a customizable tool for molecular biology.

Multispecies transcriptomes reveal core fruit development genes

During angiosperm evolution there have been repeated transitions from an ancestral dry fruit to a derived fleshy fruit, often with dramatic ecological and economic consequences. Following the transition to fleshy fruits, domestication may also dramatically alter the fruit phenotype via artificial selection. Although the morphologies of these fruits are well documented, relatively less is known about the molecular basis of these developmental and evolutionary shifts. We generated RNA-seq libraries from pericarp tissue of desert tobacco and both cultivated and wild tomato species at common developmental time points and combined this with corresponding, publicly available data from Arabidopsis and melon. With this broadly sampled dataset consisting of dry/fleshy fruits and wild/domesticated species, we applied novel bioinformatic methods to investigate conserved and divergent patterns of gene expression during fruit development and evolution. A small set of 121 orthologous "core" fruit development genes show a common pattern of expression across all five species. These include key players in developmental patterning such as orthologs of KNOLLE, PERIANTHIA, and ARGONAUTE7. GO term enrichment suggests that these genes function in basic cell division processes, cell wall biosynthesis, and developmental patterning. We furthermore uncovered a number of "accessory" genes with conserved expression patterns within but not among fruit types, and whose functional enrichment highlights the conspicuous differences between these phenotypic classes. We observe striking conservation of gene expression patterns despite large evolutionary distances, and dramatic phenotypic shifts, suggesting a conserved function for a small subset of core fruit development genes.

Draft Genome Sequence of the Tomato Stem Endophyte Bacillus safensis TS3

Tomato stem endophyte Bacillus safensis TS3 was isolated from surface-sterilized stems of greenhouse tomato plants. Here, we sequenced the complete genome of this strain to understand the molecular mechanisms underlying its beneficial activities.

α-cyanobacteria possessing form IA RuBisCO globally dominate aquatic habitats

RuBisCO (ribulose 1,5-bisphosphate carboxylase/oxygenase) is one the most abundant enzymes on Earth. Virtually all food webs depend on its activity to supply fixed carbon. In aerobic environments, RuBisCO struggles to distinguish efficiently between CO2 and O2. To compensate, organisms have evolved convergent solutions to concentrate CO2 around the active site. The genetic engineering of such inorganic carbon concentrating mechanisms (CCMs) into plants could help facilitate future global food security for humankind. In bacteria, the carboxysome represents one such CCM component, of which two independent forms exist: α and β. Cyanobacteria are important players in the planet's carbon cycle and the vast majority of the phylum possess a β-carboxysome, including most cyanobacteria used as laboratory models. The exceptions are the exclusively marine Prochlorococcus and Synechococcus that numerically dominate open ocean systems. However, the reason why marine systems favor an α-form is currently unknown. Here, we report the genomes of 58 cyanobacteria, closely related to marine Synechococcus that were isolated from freshwater lakes across the globe. We find all these isolates possess α-carboxysomes accompanied by a form 1A RuBisCO. Moreover, we demonstrate α-cyanobacteria dominate freshwater lakes worldwide. Hence, the paradigm of a separation in carboxysome type across the salinity divide does not hold true, and instead the α-form dominates all aquatic systems. We thus question the relevance of β-cyanobacteria as models for aquatic systems at large and pose a hypothesis for the reason for the success of the α-form in nature.

Phage-Plasmids Spread Antibiotic Resistance Genes through Infection and Lysogenic Conversion

Antibiotic resistance is rapidly spreading via the horizontal transfer of resistance genes in mobile genetic elements. While plasmids are key drivers of this process, few integrative phages encode antibiotic resistance genes. Here, we find that phage-plasmids, elements that are both phages and plasmids, often carry antibiotic resistance genes. We found 60 phage-plasmids with 184 antibiotic resistance genes, providing resistance for broad-spectrum-cephalosporins, carbapenems, aminoglycosides, fluoroquinolones, and colistin. These genes are in a few hot spots, seem to have been cotranslocated with transposable elements, and are often in class I integrons, which had not been previously found in phages. We tried to induce six phage-plasmids with resistance genes (including four with resistance integrons) and succeeded in five cases. Other phage-plasmids and integrative prophages were coinduced in these experiments. As a proof of concept, we focused on a P1-like element encoding an extended spectrum β-lactamase, bla_CTX-M-55. After induction, we confirmed that it is capable of infecting and converting four other E. coli strains. Its reinduction led to the further conversion of a sensitive strain, confirming that it is a fully functional phage. This study shows that phage-plasmids carry a large diversity of clinically relevant antibiotic resistance genes that they can transfer across bacteria. As plasmids, these elements seem plastic and capable of acquiring genes from other plasmids. As phages, they may provide novel paths of transfer for resistance genes because they can infect bacteria that are distant in time and space from the original host. As a matter of alarm, they may also mediate transfer to other types of phages.

Deep-branching ANME-1c archaea grow at the upper temperature limit of anaerobic oxidation of methane

In seafloor sediments, the anaerobic oxidation of methane (AOM) consumes most of the methane formed in anoxic layers, preventing this greenhouse gas from reaching the water column and finally the atmosphere. AOM is performed by syntrophic consortia of specific anaerobic methane-oxidizing archaea (ANME) and sulfate-reducing bacteria (SRB). Cultures with diverse AOM partners exist at temperatures between 12°C and 60°C. Here, from hydrothermally heated sediments of the Guaymas Basin, we cultured deep-branching ANME-1c that grow in syntrophic consortia with Thermodesulfobacteria at 70°C. Like all ANME, ANME-1c oxidize methane using the methanogenesis pathway in reverse. As an uncommon feature, ANME-1c encode a nickel-iron hydrogenase. This hydrogenase has low expression during AOM and the partner Thermodesulfobacteria lack hydrogen-consuming hydrogenases. Therefore, it is unlikely that the partners exchange hydrogen during AOM. ANME-1c also does not consume hydrogen for methane formation, disputing a recent hypothesis on facultative methanogenesis. We hypothesize that the ANME-1c hydrogenase might have been present in the common ancestor of ANME-1 but lost its central metabolic function in ANME-1c archaea. For potential direct interspecies electron transfer (DIET), both partners encode and express genes coding for extracellular appendages and multiheme cytochromes. Thermodesulfobacteria encode and express an extracellular pentaheme cytochrome with high similarity to cytochromes of other syntrophic sulfate-reducing partner bacteria. ANME-1c might associate specifically to Thermodesulfobacteria, but their co-occurrence is so far only documented for heated sediments of the Gulf of California. However, in the deep seafloor, sulfate-methane interphases appear at temperatures up to 80°C, suggesting these as potential habitats for the partnership of ANME-1c and Thermodesulfobacteria.

Elucidating the picocyanobacteria salinity divide through ecogenomics of new freshwater isolates

BACKGROUND

Cyanobacteria are the major prokaryotic primary producers occupying a range of aquatic habitats worldwide that differ in levels of salinity, making them a group of interest to study one of the major unresolved conundrums in aquatic microbiology which is what distinguishes a marine microbe from a freshwater one? We address this question using ecogenomics of a group of picocyanobacteria (cluster 5) that have recently evolved to inhabit geographically disparate salinity niches. Our analysis is made possible by the sequencing of 58 new genomes from freshwater representatives of this group that are presented here, representing a 6-fold increase in the available genomic data.

RESULTS

Overall, freshwater strains had larger genomes (≈2.9 Mb) and %GC content (≈64%) compared to brackish (2.69 Mb and 64%) and marine (2.5 Mb and 58.5%) isolates. Genomic novelties/differences across the salinity divide highlighted acidic proteomes and specific salt adaptation pathways in marine isolates (e.g., osmolytes/compatible solutes - glycine betaine/ggp/gpg/gmg clusters and glycerolipids glpK/glpA), while freshwater strains possessed distinct ion/potassium channels, permeases (aquaporin Z), fatty acid desaturases, and more neutral/basic proteomes. Sulfur, nitrogen, phosphorus, carbon (photosynthesis), or stress tolerance metabolism while showing distinct genomic footprints between habitats, e.g., different types of transporters, did not obviously translate into major functionality differences between environments. Brackish microbes show a mixture of marine (salt adaptation pathways) and freshwater features, highlighting their transitional nature.

CONCLUSIONS

The plethora of freshwater isolates provided here, in terms of trophic status preference and genetic diversity, exemplifies their ability to colonize ecologically diverse waters across the globe. Moreover, a trend towards larger and more flexible/adaptive genomes in freshwater picocyanobacteria may hint at a wider number of ecological niches in this environment compared to the relatively homogeneous marine system.

Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event

The anaerobic ammonium oxidation (anammox) bacteria can transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remain enigmatic. Here, we performed a comprehensive phylogenomic and molecular clock analysis of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which include implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based molecular dating approach, we estimated a consistent origin of anammox bacteria at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32–2.5 Ga) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress adaptation, bioenergetics, and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria.

Discovery of lignin-transforming bacteria and enzymes in thermophilic environments using stable isotope probing

Characterizing microorganisms and enzymes involved in lignin biodegradation in thermal ecosystems can identify thermostable biocatalysts. We integrated stable isotope probing (SIP), genome-resolved metagenomics, and enzyme characterization to investigate the degradation of high-molecular weight, ¹³C-ring-labeled synthetic lignin by microbial communities from moderately thermophilic hot spring sediment (52 °C) and a woody "hog fuel" pile (53 and 62 °C zones). ¹³C-Lignin degradation was monitored using IR-GCMS of ¹³CO₂, and isotopic enrichment of DNA was measured with UHLPC-MS/MS. Assembly of 42 metagenomic libraries (72 Gb) yielded 344 contig bins, from which 125 draft genomes were produced. Fourteen genomes were significantly enriched with ¹³C from lignin, including genomes of Actinomycetes (Thermoleophilaceae, Solirubrobacteraceae, Rubrobacter sp.), Firmicutes (Kyrpidia sp., Alicyclobacillus sp.) and Gammaproteobacteria (Steroidobacteraceae). We employed multiple approaches to screen genomes for genes encoding putative ligninases and pathways for aromatic compound degradation. Our analysis identified several novel laccase-like multi-copper oxidase (LMCO) genes in ¹³C-enriched genomes. One of these LMCOs was heterologously expressed and shown to oxidize lignin model compounds and minimally transformed lignin. This study elucidated bacterial lignin depolymerization and mineralization in thermal ecosystems, establishing new possibilities for the efficient valorization of lignin at elevated temperature.