A Model of Genome Size Evolution for Prokaryotes in Stable and Fluctuating Environments

The cost of adaptability: resource availability constrains functional stability under pulsed disturbances

Global change exposes ecosystems to changes in the frequency, magnitude, and concomitancy of disturbances, which impact the composition and functioning of these systems. Here, we experimentally evaluate the effects of salinity disturbances and eutrophication on bacterial communities from coastal ecosystems. The functional stability of these communities is critically important for maintaining water quality, productivity, and ecosystem services, such as fishery yields. Microbial functional stability can be maintained via resistance and resilience, which are reflected in genomic traits such as genome size and codon usage bias and may be linked to metabolic costs. However, little is known about the mechanisms that select these traits under varying nutrient regimes. To study the impact of pulsed disturbances on community assembly and functioning depending on metabolic costs, we performed a 41-day pulse disturbance experiment across two levels of resource availability. Our setup triggered stochastic community re-assembly processes in all treatments. In contrast, we observed consistent and resource availability-dependent patterns of superordinate community functioning and structural patterns, such as functional resistance in response to disturbances, genomic trait distributions, and species diversity. Predicted genomic traits reflected the selection for taxa possessing resistant- and resilience-related traits, particularly under high nutrient availability. Our findings are a step toward unraveling the compositional and genomic underpinnings of functional resistance in microbial communities after exposure to consecutive pulse disturbances. Our work demonstrates how resource availability alleviates metabolic constraints on resistance and resilience, and this has important consequences for predicting water quality and ecosystem productivity of environments exposed to global change.

IMPORTANCE

Understanding the communities' responses to disturbances is a prerequisite to predicting ecosystem dynamics and, thus, highly relevant considering global change. Microbial communities play key roles in numerous ecosystem functions and services, and the large diversity, rapid growth, and phenotypic plasticity of microorganisms are thought to allow high resistance and resilience. While potential metabolic costs associated with adaptations to fluctuating environments have been debated, little evidence supports trade-offs between resource availability, resistance, and resilience. Here, we experimentally assessed the compositional and functional responses of an aquatic microbial model community to disturbances and systematically manipulated resource availability. Our results demonstrate that the capacity to tolerate environmental fluctuations is constrained by resource availability and reflected in the selection of genomic traits.

r/K selection of GC content in prokaryotes

The guanine/cytosine (GC) content of prokaryotic genomes is species-specific, taking values from 16% to 77%. This diversity of selection for GC content remains contentious. We analyse the correlations between GC content and a range of phenotypic and genotypic data in thousands of prokaryotes. GC content integrates well with these traits into r/K selection theory when phenotypic plasticity is considered. High GC-content prokaryotes are r-strategists with cheaper descendants thanks to a lower average amino acid metabolic cost, colonize unstable environments thanks to flagella and a bacillus form and are generalists in terms of resource opportunism and their defence mechanisms. Low GC content prokaryotes are K-strategists specialized for stable environments that maintain homeostasis via a high-cost outer cell membrane and endospore formation as a response to nutrient deprivation, and attain a higher nutrient-to-biomass yield. The lower proteome cost of high GC content prokaryotes is driven by the association between GC-rich codons and cheaper amino acids in the genetic code, while the correlation between GC content and genome size may be partly due to functional diversity driven by r/K selection. In all, molecular diversity in the GC content of prokaryotes may be a consequence of ecological r/K selection.

Bacterial genome size and gene functional diversity negatively correlate with taxonomic diversity along a pH gradient

Bacterial gene repertoires reflect adaptive strategies, contribute to ecosystem functioning and are limited by genome size. However, gene functional diversity does not necessarily correlate with taxonomic diversity because average genome size may vary by community. Here, we analyse gene functional diversity (by shotgun metagenomics) and taxonomic diversity (by 16S rRNA gene amplicon sequencing) to investigate soil bacterial communities along a natural pH gradient in 12 tropical, subtropical, and temperate forests. We find that bacterial average genome size and gene functional diversity decrease, whereas taxonomic diversity increases, as soil pH rises from acid to neutral; as a result, bacterial taxonomic and functional diversity are negatively correlated. The gene repertoire of acid-adapted oligotrophs is enriched in functions of signal transduction, cell motility, secretion system, and degradation of complex compounds, while that of neutral pH-adapted copiotrophs is enriched in functions of energy metabolism and membrane transport. Our results indicate that a mismatch between taxonomic and functional diversity can arise when environmental factors (such as pH) select for adaptive strategies that affect genome size distributions.

Genomics of soil depth niche partitioning in the Thaumarchaeota family Gagatemarchaeaceae

Knowledge of deeply-rooted non-ammonia oxidising Thaumarchaeota lineages from terrestrial environments is scarce, despite their abundance in acidic soils. Here, 15 new deeply-rooted thaumarchaeotal genomes were assembled from acidic topsoils (0-15 cm) and subsoils (30-60 cm), corresponding to two genera of terrestrially prevalent Gagatemarchaeaceae (previously known as thaumarchaeotal Group I.1c) and to a novel genus of heterotrophic terrestrial Thaumarchaeota. Unlike previous predictions, metabolic annotations suggest Gagatemarchaeaceae perform aerobic respiration and use various organic carbon sources. Evolutionary divergence between topsoil and subsoil lineages happened early in Gagatemarchaeaceae history, with significant metabolic and genomic trait differences. Reconstruction of the evolutionary mechanisms showed that the genome expansion in topsoil Gagatemarchaeaceae resulted from extensive early lateral gene acquisition, followed by progressive gene duplication throughout evolutionary history. Ancestral trait reconstruction using the expanded genomic diversity also did not support the previous hypothesis of a thermophilic last common ancestor of the ammonia-oxidising archaea. Ultimately, this study provides a good model for studying mechanisms driving niche partitioning between spatially related ecosystems.

Genomic and in-vitro characteristics of a novel strain Lacticaseibacillus chiayiensis AACE3 isolated from fermented blueberry

Numerous different species of LAB are used in different fields due to their unique characteristics. However, Lacticaseibacillus chiayiensis, a newly established species in 2018, has limited microorganism resources, and lacks comprehensive evaluations of its properties. In this study, L. chiayiensis AACE3, isolated from fermented blueberry, was evaluated by genomic analysis and in vitro assays of the properties. The genome identified genes associated with biofilm formation (luxS, ccpA, brpA), resistance to oxidative stress (tpx, trxA, trxB, hslO), tolerance to acidic conditions (dltA, dltC), resistance to unfavorable osmotic pressure (opuBB, gbuA, gbuB, gbuC), and adhesion (luxS, dltA, dltC). The AACE3 showed 112 unique genes, relative to the other three L. chiayiensis strains. Among them, the presence of genes such as clpP, pepO, and feoA suggests a possible advantage of AACE3 over other L. chiayiensis in terms of environmental adaptation. In vitro evaluation of the properties revealed that AACE3 had robust antibacterial activity against eight common pathogens: Streptococcus agalactiae, Staphylococcus aureus, Escherichia coli, Salmonella enteritidis, Salmonella choleraesuis, Shigella flexneri, Pseudomonas aeruginosa, and Klebsiella pneumoniae. In addition, AACE3 showed more than 80% survival rate in all tests simulating gastrointestinal fluid, and it exhibited high antioxidant capacity. Interestingly, the cell culture supernatant was superior to intact organisms and ultrasonically crushed bacterial extracts in all tests of antioxidant capacity. These results suggested that the antioxidant capacity may originate from certain metabolites and extracellular enzymes produced by AACE3. Moreover, AACE3 was a moderate biofilm producer due to the self-agglomeration effect. Taken together, L. chiayiensis AACE3 appears to be a candidate strain for combating the growing incidence of pathogen infections and antioxidant production.

A social niche breadth score reveals niche range strategies of generalists and specialists

Generalists can survive in many environments, whereas specialists are restricted to a single environment. Although a classical concept in ecology, niche breadth has remained challenging to quantify for microorganisms because it depends on an objective definition of the environment. Here, by defining the environment of a microorganism as the community it resides in, we integrated information from over 22,000 environmental sequencing samples to derive a quantitative measure of the niche, which we call social niche breadth. At the level of genera, we explored niche range strategies throughout the prokaryotic tree of life. We found that social generalists include opportunists that stochastically dominate local communities, whereas social specialists are stable but low in abundance. Social generalists have a more diverse and open pan-genome than social specialists, but we found no global correlation between social niche breadth and genome size. Instead, we observed two distinct evolutionary strategies, whereby specialists have relatively small genomes in habitats with low local diversity, but relatively large genomes in habitats with high local diversity. Together, our analysis shines data-driven light on microbial niche range strategies.

Genome size variation within Crithmum maritimum: Clues on the colonization of insular environments

Angiosperms present an astonishing diversity of genome sizes that can vary intra- or interspecifically. The remarkable new cytogenomic data shed some light on our understanding of evolution, but few studies were performed with insular and mainland populations to test possible correlations with dispersal, speciation, and adaptations to insular environments. Here, patterns of cytogenomic diversity were assessed among geographic samples (ca. 114) of Crithmum maritimum (Apiaceae), collected across the Azores and Madeira archipelagos, as well as in adjacent continental areas of Portugal. Using flow cytometry, the results indicated a significant intraspecific genome size variation, spanning from reduced sizes in the insular populations to larger ones in the mainland populations. Moreover, there was a tendency for an increase in genome size along the mainland populations, associated with lower temperatures, higher precipitation, and lower precipitation seasonality. However, this gradient might be the result of historic phylogeographical events associated with previous dispersal and extinction of local populations. Overall, our findings provided evidence that smaller genome sizes might play a critical role in the colonization of islands, corroborating other studies that argue that organisms with smaller genomes use fewer resources, having a selective advantage under insular environments. Although further studies are needed to improve our understanding of the mechanisms underlying genome size evolution on islands, conservation strategies must be promoted to protect the rich cytogenomic diversity found among C. maritimum populations, which occur in coastal areas that are particularly threatened by human activity, pollution, invasive species, and climate changes.

Abiotic selection of microbial genome size in the global ocean

Strong purifying selection is considered a major evolutionary force behind small microbial genomes in the resource-poor photic ocean. However, very little is currently known about how the size of prokaryotic genomes evolves in the global ocean and whether patterns reflect shifts in resource availability in the epipelagic and relatively stable deep-sea environmental conditions. Using 364 marine microbial metagenomes, we investigate how the average genome size of uncultured planktonic prokaryotes varies across the tropical and polar oceans to the hadal realm. We find that genome size is highest in the perennially cold polar ocean, reflecting elongation of coding genes and gene dosage effects due to duplications in the interior ocean microbiome. Moreover, the rate of change in genome size due to temperature is 16-fold higher than with depth up to 200 m. Our results demonstrate how environmental factors can influence marine microbial genome size selection and ecological strategies of the microbiome.

Mixotrophy broadens the ecological niche range of the iron oxidizer Sideroxydans sp. CL21 isolated from an iron-rich peatland

Sideroxydans sp. CL21 is a microaerobic, acid-tolerant Fe(II)-oxidizer, isolated from the Schlöppnerbrunnen fen. Since the genome size of Sideroxydans sp. CL21 is 21% larger than that of the neutrophilic Sideroxydans lithotrophicus ES-1, we hypothesized that strain CL21 contains additional metabolic traits to thrive in the fen. The common genomic content of both strains contains homologs of the putative Fe(II) oxidation genes, mtoAB and cyc2. A large part of the accessory genome in strain CL21 contains genes linked to utilization of alternative electron donors, including NiFe uptake hydrogenases, and genes encoding lactate uptake and utilization proteins, motility and biofilm formation, transposable elements, and pH homeostasis mechanisms. Next, we incubated the strain in different combinations of electron donors and characterized the fen microbial communities. Sideroxydans spp. comprised 3.33% and 3.94% of the total relative abundance in the peatland soil and peatland water, respectively. Incubation results indicate Sideroxydans sp. CL21 uses H2 and thiosulfate, while lactate only enhances growth when combined with Fe, H2, or thiosulfate. Rates of H2 utilization were highest in combination with other substrates. Thus, Sideroxydans sp. CL21 is a mixotroph, growing best by simultaneously using substrate combinations, which helps to thrive in dynamic and complex habitats.

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.

Trait-trait relationships and tradeoffs vary with genome size in prokaryotes

We report genomic traits that have been associated with the life history of prokaryotes and highlight conflicting findings concerning earlier observed trait correlations and tradeoffs. In order to address possible explanations for these contradictions we examined trait–trait variations of 11 genomic traits from ~18,000 sequenced genomes. The studied trait–trait variations suggested: (i) the predominance of two resistance and resilience-related orthogonal axes and (ii) at least in free living species with large effective population sizes whose evolution is little affected by genetic drift an overlap between a resilience axis and an oligotrophic-copiotrophic axis. These findings imply that resistance associated traits of prokaryotes are globally decoupled from resilience related traits and in the case of free-living communities also from traits associated with resource availability. However, further inspection of pairwise scatterplots showed that resistance and resilience traits tended to be positively related for genomes up to roughly five million base pairs and negatively for larger genomes. Genome size distributions differ across habitats and our findings therefore point to habitat dependent tradeoffs between resistance and resilience. This in turn may preclude a globally consistent assignment of prokaryote genomic traits to the competitor - stress-tolerator - ruderal (CSR) schema that sorts species depending on their location along disturbance and productivity gradients into three ecological strategies and may serve as an explanation for conflicting findings from earlier studies. All reviewed genomic traits featured significant phylogenetic signals and we propose that our trait table can be applied to extrapolate genomic traits from taxonomic marker genes. This will enable to empirically evaluate the assembly of these genomic traits in prokaryotic communities from different habitats and under different productivity and disturbance scenarios as predicted via the resistance-resilience framework formulated here.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Complete Genome Sequencing and Comparative Genomics of Three Potential Probiotic Strains, Lacticaseibacillus casei FBL6, Lacticaseibacillus chiayiensis FBL7, and Lacticaseibacillus zeae FBL8

Lacticaseibacillus casei, Lacticaseibacillus chiayiensis, and Lacticaseibacillus zeae are very closely related Lacticaseibacillus species. L. casei has long been proposed as a probiotic, whereas studies on functional characterization for L. chiayiensis and L. zeae are some compared to L. casei. In this study, L. casei FBL6, L. chiayiensis FBL7, and L. zeae FBL8 were isolated from raw milk, and their probiotic properties were investigated. Genomic analysis demonstrated the role of L. chiayiensis and L. zeae as probiotic candidates. The three strains were tolerant to acid and bile salt, with inhibitory action against pathogenic bacterial strains and capacity of antioxidants. Complete genome sequences of the three strains were analyzed to highlight the probiotic properties at the genetic level, which results in the discovery of genes corresponding to phenotypic characterization. Moreover, genes known to confer probiotic characteristics were identified, including genes related to biosynthesis, defense machinery, adhesion, and stress adaptation. The comparative genomic analysis with other available genomes revealed 256, 214, and 32 unique genes for FBL6, FBL7, and FBL8, respectively. These genomes contained individual genes encoding proteins that are putatively involved in carbohydrate transport and metabolism, prokaryotic immune system for antiviral defense, and physiological control processes. In particular, L. casei FBL6 had a bacteriocin gene cluster that was not present in other genomes of L. casei, resulting in this strain may exhibit a wide range of antimicrobial activity compared to other L. casei strains. Our data can help us understand the probiotic functionalities of the three strains and suggest that L. chiayiensis and L. zeae species, which are closely related to L. casei, can also be considered as novel potential probiotic candidate strains.

The genus Caulobacter and its role in plant microbiomes

Recent omics approaches have revealed the prevalent microbial taxa that constitute the microbiome of various plant species. Across global scales and environmental conditions, strains belonging to the bacterial genus Caulobacter have consistently been found in association with various plant species. Aligned with agroecological relevance and biotechnological advances, many scientific communications have demonstrated that several Caulobacter strains (spanning several Caulobacter species) harbor the potential to enhance plant biomass for various plant species ranging from Arabidopsis to Citrullus and Zea mays. In the past several years, co-occurrence data have driven mechanistically resolved communications about select Caulobacter-plant interactions. Given the long-standing history of Caulobacter as a model organism for cell cycle regulation, genetic studies, and the prevalence of Caulobacter species in various plant microbiomes, the genus Caulobacter offers researchers a unique opportunity to leverage for investigating plant-microbe interactions and realizing targeted biotechnological applications. In this review, recent developments regarding Caulobacter-plant interactions are presented in terms of model utility for future biotechnological investigations.

A Genomic Perspective Across Earth's Microbiomes Reveals That Genome Size in Archaea and Bacteria Is Linked to Ecosystem Type and Trophic Strategy

Our view of genome size in Archaea and Bacteria has remained skewed as the data has been dominated by genomes of microorganisms that have been cultivated under laboratory settings. However, the continuous effort to catalog Earth's microbiomes, specifically propelled by recent extensive work on uncultivated microorganisms, provides an opportunity to revise our perspective on genome size distribution. We present a meta-analysis that includes 26,101 representative genomes from 3 published genomic databases; metagenomic assembled genomes (MAGs) from GEMs and stratfreshDB, and isolates from GTDB. Aquatic and host-associated microbial genomes present on average the smallest estimated genome sizes (3.1 and 3.0 Mbp, respectively). These are followed by terrestrial microbial genomes (average 3.7 Mbp), and genomes from isolated microorganisms (average 4.3 Mbp). On the one hand, aquatic and host-associated ecosystems present smaller genomes sizes in genera of phyla with genome sizes above 3 Mbp. On the other hand, estimated genome size in phyla with genomes under 3 Mbp showed no difference between ecosystems. Moreover, we observed that when using 95% average nucleotide identity (ANI) as an estimator for genetic units, only 3% of MAGs cluster together with genomes from isolated microorganisms. Although there are potential methodological limitations when assembling and binning MAGs, we found that in genome clusters containing both environmental MAGs and isolate genomes, MAGs were estimated only an average 3.7% smaller than isolate genomes. Even when assembly and binning methods introduce biases, estimated genome size of MAGs and isolates are very similar. Finally, to better understand the ecological drivers of genome size, we discuss on the known and the overlooked factors that influence genome size in different ecosystems, phylogenetic groups, and trophic strategies.

Genomics landscape of 185 Streptococcus thermophilus and identification of fermentation biomarkers

Streptococcus (S.) thermophilus, an indispensable dairy starter, has been used in autochthonous as well as industrial milk fermentation. However, the genetic architecture underlying S. thermophilus traits and phenotypes is largely unknown. Here, we sequenced 185 S. thermophilus strains, isolated from natural fermented dairy products of China and Mongolia and used comparative genomic and genome wide association study to provide novel point for genetic architecture underlying its traits and phenotypes. Genome analysis of S. thermophilus showed association of phylogeny with environmental and phenotypic features and revealed clades with high acid production potential or with substantial genome decay. A few S. thermophilus isolated from areas with high chloramphenicol emissions had a chloramphenicol-resistant gene CatB8. Most importantly, we defined a growth score and identified a missense mutation G1118698T located at the gene AcnA that were both predictive of acidification capability of S. thermophilus. Our findings provide novel insight in S. thermophilus genetic traits, antibiotic resistant and predictive of acidification capability which both may had huge help in culture starter screening.

Frequency of mispackaging of Prochlorococcus DNA by cyanophage

Prochlorococcus cells are the numerically dominant phototrophs in the open ocean. Cyanophages that infect them are a notable fraction of the total viral population in the euphotic zone, and, as vehicles of horizontal gene transfer, appear to drive their evolution. Here we examine the propensity of three cyanophages-a podovirus, a siphovirus, and a myovirus-to mispackage host DNA in their capsids while infecting Prochlorococcus, the first step in phage-mediated horizontal gene transfer. We find the mispackaging frequencies are distinctly different among the three phages. Myoviruses mispackage host DNA at low and seemingly fixed frequencies, while podo- and siphoviruses vary in their mispackaging frequencies by orders of magnitude depending on growth light intensity. We link this difference to the concentration of intracellular reactive oxygen species and protein synthesis rates, both parameters increasing in response to higher light intensity. Based on our findings, we propose a model of mispackaging frequency determined by the imbalance between the production of capsids and the number of phage genome copies during infection: when protein synthesis rate increase to levels that the phage cannot regulate, they lead to an accumulation of empty capsids, in turn triggering more frequent host DNA mispackaging errors.

Computational animal welfare: towards cognitive architecture models of animal sentience, emotion and wellbeing

To understand animal wellbeing, we need to consider subjective phenomena and sentience. This is challenging, since these properties are private and cannot be observed directly. Certain motivations, emotions and related internal states can be inferred in animals through experiments that involve choice, learning, generalization and decision-making. Yet, even though there is significant progress in elucidating the neurobiology of human consciousness, animal consciousness is still a mystery. We propose that computational animal welfare science emerges at the intersection of animal behaviour, welfare and computational cognition. By using ideas from cognitive science, we develop a functional and generic definition of subjective phenomena as any process or state of the organism that exists from the first-person perspective and cannot be isolated from the animal subject. We then outline a general cognitive architecture to model simple forms of subjective processes and sentience. This includes evolutionary adaptation which contains top-down attention modulation, predictive processing and subjective simulation by re-entrant (recursive) computations. Thereafter, we show how this approach uses major characteristics of the subjective experience: elementary self-awareness, global workspace and qualia with unity and continuity. This provides a formal framework for process-based modelling of animal needs, subjective states, sentience and wellbeing.

Insights into the Metabolism and Evolution of the Genus Acidiphilium, a Typical Acidophile in Acid Mine Drainage

Here, we report three new Acidiphilium genomes, reclassified existing Acidiphilium species, and performed the first comparative genomic analysis on Acidiphilium in an attempt to address the metabolic potential, ecological functions, and evolutionary history of the genus Acidiphilium In the genomes of Acidiphilium, we found an abundant repertoire of horizontally transferred genes (HTGs) contributing to environmental adaption and metabolic expansion, including genes conferring photosynthesis (puf, puh), CO2 assimilation (rbc), capacity for methane metabolism (mmo, mdh, frm), nitrogen source utilization (nar, cyn, hmp), sulfur compound utilization (sox, psr, sqr), and multiple metal and osmotic stress resistance capacities (czc, cop, ect). Additionally, the predicted donors of horizontal gene transfer were present in a cooccurrence network of Acidiphilium Genome-scale positive selection analysis revealed that 15 genes contained adaptive mutations, most of which were multifunctional and played critical roles in the survival of extreme conditions. We proposed that Acidiphilium originated in mild conditions and adapted to extreme environments such as acidic mineral sites after the acquisition of many essential functions.IMPORTANCE Extremophiles, organisms that thrive in extreme environments, are key models for research on biological adaption. They can provide hints for the origin and evolution of life, as well as improve the understanding of biogeochemical cycling of elements. Extremely acidophilic bacteria such as Acidiphilium are widespread in acid mine drainage (AMD) systems, but the metabolic potential, ecological functions, and evolutionary history of this genus are still ambiguous. Here, we sequenced the genomes of three new Acidiphilium strains and performed comparative genomic analysis on this extremely acidophilic bacterial genus. We found in the genomes of Acidiphilium an abundant repertoire of horizontally transferred genes (HTGs) contributing to environmental adaption and metabolic ability expansion, as indicated by phylogenetic reconstruction and gene context comparison. This study has advanced our understanding of microbial evolution and biogeochemical cycling in extreme niches.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

Small and mighty: adaptation of superphylum Patescibacteria to groundwater environment drives their genome simplicity

BACKGROUND

The newly defined superphylum Patescibacteria such as Parcubacteria (OD1) and Microgenomates (OP11) has been found to be prevalent in groundwater, sediment, lake, and other aquifer environments. Recently increasing attention has been paid to this diverse superphylum including > 20 candidate phyla (a large part of the candidate phylum radiation, CPR) because it refreshed our view of the tree of life. However, adaptive traits contributing to its prevalence are still not well known.

RESULTS

Here, we investigated the genomic features and metabolic pathways of Patescibacteria in groundwater through genome-resolved metagenomics analysis of > 600 Gbp sequence data. We observed that, while the members of Patescibacteria have reduced genomes (~ 1 Mbp) exclusively, functions essential to growth and reproduction such as genetic information processing were retained. Surprisingly, they have sharply reduced redundant and nonessential functions, including specific metabolic activities and stress response systems. The Patescibacteria have ultra-small cells and simplified membrane structures, including flagellar assembly, transporters, and two-component systems. Despite the lack of CRISPR viral defense, the bacteria may evade predation through deletion of common membrane phage receptors and other alternative strategies, which may explain the low representation of prophage proteins in their genomes and lack of CRISPR. By establishing the linkages between bacterial features and the groundwater environmental conditions, our results provide important insights into the functions and evolution of this CPR group.

CONCLUSIONS

We found that Patescibacteria has streamlined many functions while acquiring advantages such as avoiding phage invasion, to adapt to the groundwater environment. The unique features of small genome size, ultra-small cell size, and lacking CRISPR of this large lineage are bringing new understandings on life of Bacteria. Our results provide important insights into the mechanisms for adaptation of the superphylum in the groundwater environments, and demonstrate a case where less is more, and small is mighty.

Revisiting the hypothesis of an energetic barrier to genome complexity between eukaryotes and prokaryotes

The absence of genome complexity in prokaryotes, being the evolutionary precursors to eukaryotic cells comprising all complex life (the prokaryote-eukaryote divide), is a long-standing question in evolutionary biology. A previous study hypothesized that the divide exists because prokaryotic genome size is constrained by bioenergetics (prokaryotic power per gene or genome being significantly lower than eukaryotic ones). However, this hypothesis was evaluated using a relatively small dataset due to lack of data availability at the time, and is therefore controversial. Accordingly, we constructed a larger dataset of genomes, metabolic rates, cell sizes and ploidy levels to investigate whether an energetic barrier to genome complexity exists between eukaryotes and prokaryotes while statistically controlling for the confounding effects of cell size and phylogenetic signals. Notably, we showed that the differences in bioenergetics between prokaryotes and eukaryotes were less significant than those previously reported. More importantly, we found a limited contribution of power per genome and power per gene to the prokaryote-eukaryote dichotomy. Our findings indicate that the prokaryote-eukaryote divide is hard to explain from the energetic perspective. However, our findings may not entirely discount the traditional hypothesis; in contrast, they indicate the need for more careful examination.

Genome Sequencing of Mesonia algae K4-1 Reveals Its Adaptation to the Arctic Ocean

The special ecological environment of the Arctic has brought about a large number of salt-tolerant and psychrotolerant microorganisms. We isolated two culturable bacterial strains of the genus Mesonia; one from the Arctic ocean, Mesonia algae K4-1, and one from the tropical sea, Mesonia sp. HuA40. Our genome analysis and phenotypic experiments indicated that Mesonia algae K4-1 is a moderately halophilic and psychrophilic bacterium. Mesonia algae K4-1 can tolerate 3–14% NaCl and grow at a wide range of temperatures from 4 to 50°C. Mesonia sp. HuA40 is a mesophilic bacterium that can only grow with 3–9% NaCl. In addition, the salt adaptation strategy of Mesonia algae K4-1 accumulates organic osmolytes in the cell. RNA helicases, glutathione and organic compatible solutes may play important roles in maintaining the metabolism and physiological function of Mesonia algae K4-1 under cold stress. Moreover, the ability of Mesonia algae K4-1 to adapt to an oligotrophic marine environment is likely due to the synthesis of a large number of extracellular polysaccharides and the secretion of various families of extracellular proteases. This study systematically analyzed the relationship between genomic differentiation and environmental factors of the Mesonia genus and revealed the possible adaptation mechanism of Mesonia algae K4-1 in the extreme Arctic marine environment at the genomic level.

Fast and accurate average genome size and 16S rRNA gene average copy number computation in metagenomic data

Background

Metagenomics caused a quantum leap in microbial ecology. However, the inherent size and complexity of metagenomic data limit its interpretation. The quantification of metagenomic traits in metagenomic analysis workflows has the potential to improve the exploitation of metagenomic data. Metagenomic traits are organisms’ characteristics linked to their performance. They are measured at the genomic level taking a random sample of individuals in a community. As such, these traits provide valuable information to uncover microorganisms’ ecological patterns. The Average Genome Size (AGS) and the 16S rRNA gene Average Copy Number (ACN) are two highly informative metagenomic traits that reflect microorganisms’ ecological strategies as well as the environmental conditions they inhabit.

Results

Here, we present the ags.sh and acn.sh tools, which analytically derive the AGS and ACN metagenomic traits. These tools represent an advance on previous approaches to compute the AGS and ACN traits. Benchmarking shows that ags.sh is up to 11 times faster than state-of-the-art tools dedicated to the estimation AGS. Both ags.sh and acn.sh show comparable or higher accuracy than existing tools used to estimate these traits. To exemplify the applicability of both tools, we analyzed the 139 prokaryotic metagenomes of TARA Oceans and revealed the ecological strategies associated with different water layers.

Conclusion

We took advantage of recent advances in gene annotation to develop the ags.sh and acn.sh tools to combine easy tool usage with fast and accurate performance. Our tools compute the AGS and ACN metagenomic traits on unassembled metagenomes and allow researchers to improve their metagenomic data analysis to gain deeper insights into microorganisms’ ecology. The ags.sh and acn.sh tools are publicly available using Docker container technology at https://github.com/pereiramemo/AGS-and-ACN-tools.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3031-y) contains supplementary material, which is available to authorized users.

Genome Size Reversely Correlates With Host Plant Range in Helicoverpa Species

In organisms with very low percentages of transposable elements (TEs), genome size may positively or negatively correlate with host range, depending on whether host adaptation or host modification is the main route to host generalism. To test if this holds true for insect herbivores with greater percentages of TEs, we conducted flow cytometry to measure the endopolyploidy levels and C-values of the host modification (salivary gland and mandibular gland in head), host adaptation (midgut), and host use-independent tissues (male gonad, hemolymph, and Malpighian tubules) of the generalist Helicoverpa armigera and the head of its older specialist sister H. assulta. Larval salivary gland displayed a consecutive chain of endopolyploidy particles from 8Cx to higher than 32Cx and larval head and midgut had endopolyploidy nuclei clusters of 16Cx and 32Cx, whereas larval male gonad, hemolymph, and Malpighian tubules possessed no endopolyploidy nuclei of higher than 8Cx. The estimated genome size of the Solanaceae plant specialist H. assulta is 430 Mb, significantly larger than that of its older generalist sister Heliothis virescens (408 Mb) and those of its two generalist descendants H. armigera (394 Mb) and H. zea (363 Mb). These data not only reveal a negative correlation between host plant range and genome size in this terminal lineage, but also imply that Helicoverpa species appear to depend more on host modification than on host adaptation to achieve polyphagy.

Novel haloarchaeon Natrinema thermophila having the highest growth temperature among haloarchaea with a large genome size

Environmental temperature is one of the most important factors for the growth and survival of microorganisms. Here we describe a novel extremely halophilic archaeon (haloarchaea) designated as strain CBA1119^T isolated from solar salt. Strain CBA1119^T had the highest maximum and optimal growth temperatures (66 °C and 55 °C, respectively) and one of the largest genome sizes among haloarchaea (5.1 Mb). It also had the largest number of strain-specific pan-genome orthologous groups and unique pathways among members of the genus Natrinema in the class Halobacteria. A dendrogram based on the presence/absence of genes and a phylogenetic tree constructed based on OrthoANI values highlighted the particularities of strain CBA1119^T as compared to other Natrinema species and other haloarchaea members. The large genome of strain CBA1119^T may provide information on genes that confer tolerance to extreme environmental conditions, which may lead to the discovery of other thermophilic strains with potential applications in industrial biotechnology.

Complete Sequence of pCY-CTX, a Plasmid Carrying a Phage-Like Region and an ISEcp1-Mediated Tn2 Element from Enterobacter cloacae

A plasmid pCY-CTX carrying a phage-like backbone from an extensively drug-resistant Enterobacter cloacae strain Guangzhou-ECL001 (previously known as CY01) was identified in this study. By Illumina MiSeq 2 × 250-bp paired-end sequencing, de novo assembly, and PCR, full sequence of pCY-CTX was obtained. Plasmid pCY-CTX was a circular plasmid with a length of 116,700 bp, harboring 136 putative open reading frames with the average G + C content of 50.8%. The backbone of pCY-CTX showed high identity to previously reported phage-like plasmid pHCM2 and phage SSU5. In addition, pCY-CTX contained a distinctive ISEcp1-mediated Tn2 region with two resistance genes bla_TEM-1 and bla_CTX-M-3. Transposition unit "ISEcp1- bla_CTX-M-3- orf477" was inserted into the Tn2 structure, dividing Tn2 into two parts. This represents the first identification of a plasmid carrying a phage-like backbone and a distinctive ISEcp1-mediated Tn2 region within bla_TEM-1 and bla_CTX-M-3 in clinical E. cloacae. The finding of phage-like regions located in plasmids provides a new perspective in gene transfer associated with antimicrobial resistance.

Relating genomic characteristics to environmental preferences and ubiquity in different microbial taxa

Background

Despite the important role that microorganisms play in environmental processes, the low percentage of cultured microbes (5%) has limited, until now, our knowledge of their ecological strategies. However, the development of high-throughput sequencing has generated a huge amount of genomic and metagenomic data without the need of culturing that can be used to study ecological questions. This study aims to estimate the functional capabilities, genomic sizes and 16S copy number of different taxa in relation to their ubiquity and their environmental preferences.

Results

To achieve this goal, we compiled data regarding the presence of each prokaryotic genera in diverse environments. Then, genomic characteristics such as genome size, 16S rRNA gene copy number, and functional content of the genomes were related to their ubiquity and different environmental preferences of the corresponding taxa. The results showed clear correlations between genomic characteristics and environmental conditions.

Conclusions

Ubiquity and adaptation were linked to genome size, while 16S copy number was not directly related to ubiquity. We observed that different combinations of these two characteristics delineate the different environments. Besides, the analysis of functional classes showed some clear signatures linked to particular environments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3888-y) contains supplementary material, which is available to authorized users.

Following Rapoport's Rule: the geographic range and genome size of bacterial taxa decline at warmer latitudes

We sought to test whether stream bacterial communities conform to Rapoport's Rule, a pattern commonly observed for plants and animals whereby taxa exhibit decreased latitudinal range sizes closer to the equator. Using a DNA sequencing approach, we explored the biogeography of biofilm bacterial communities in 204 streams across a ∼1000 km latitudinal gradient. The range sizes of bacterial taxa were strongly correlated with latitude, decreasing closer to the equator, which coincided with a greater than fivefold increase in bacterial taxonomic richness. The relative richness and range size of bacteria were associated with spatially correlated variation in temperature and rainfall. These patterns were observed despite enormous variability in catchment environmental characteristics. Similar results were obtained when restricting the same analyses to native forest catchments, thereby controlling for spatial biases in land use. We analysed genomic data from ∼500 taxa detected in this study, for which data were available and found that bacterial communities at cooler latitudes also tended to possess greater potential metabolic potential. Collectively, these data provide the first evidence of latitudinal variation in the range size distributions of freshwater bacteria, a trend which may be determined, in part, by a trade-off between bacterial genome size and local variation in climatic conditions.

The multi-omics promise in context: from sequence to microbial isolate

The numbers and diversity of microbes in ecosystems within and around us is unmatched, yet most of these microorganisms remain recalcitrant to in vitro cultivation. Various high-throughput molecular techniques, collectively termed multi-omics, provide insights into the genomic structure and metabolic potential as well as activity of complex microbial communities. Nonetheless, pure or defined cultures are needed to (1) decipher microbial physiology and thus test multi-omics-based ecological hypotheses, (2) curate and improve database annotations and (3) realize novel applications in biotechnology. Cultivation thus provides context. In turn, we here argue that multi-omics information awaits integration into the development of novel cultivation strategies. This can build the foundation for a new era of omics information-guided microbial cultivation technology and reduce the inherent trial-and-error search space. This review discusses how information that can be extracted from multi-omics data can be applied for the cultivation of hitherto uncultured microorganisms. Furthermore, we summarize groundbreaking studies that successfully translated information derived from multi-omics into specific media formulations, screening techniques and selective enrichments in order to obtain novel targeted microbial isolates. By integrating these examples, we conclude with a proposed workflow to facilitate future omics-aided cultivation strategies that are inspired by the microbial complexity of the environment.

Exosomes in mammals with greater habitat variability contain more proteins and RNAs

Factors determining habitat variability are poorly understood despite possible explanations based on genome and physiology. This is because previous studies only focused on primary measures such as genome size and body size. In this study, we hypothesize that specific gene functions determine habitat variability in order to explore new factors beyond primary measures. We comprehensively evaluate the relationship between gene functions and the climate envelope while statistically controlling for potentially confounding effects by using data on the habitat range, genome, body size and metabolism of various mammals. Our analyses show that the number of proteins and RNAs contained in exosomes is predominantly associated with the climate envelope. This finding indicates the importance of exosomes to habitat range expansion of mammals and provides a new hypothesis for the relationship between the genome and habitat variability.

The effect of extrinsic mortality on genome size evolution in prokaryotes

Mortality has a significant role in prokaryotic ecology and evolution, yet the impact of variations in extrinsic mortality on prokaryotic genome evolution has received little attention. We used both mathematical and agent-based models to reveal how variations in extrinsic mortality affect prokaryotic genome evolution. Our results suggest that the genome size of bacteria increases with increased mortality. A high extrinsic mortality increases the pool of free resources and shortens life expectancy, which selects for faster reproduction, a phenotype we called 'scramblers'. This phenotype is realised by the expansion of gene families involved in nutrient acquisition and metabolism. In contrast, a low mortality rate increases an individual's life expectancy, which results in natural selection favouring tolerance to starvation when conditions are unfavourable. This leads to the evolution of small, streamlined genomes ('stayers'). Our models predict that large genomes, gene family expansion and horizontal gene transfer should be observed in prokaryotes occupying ecosystems exposed to high abiotic stress, as well as those under strong predator- and/or pathogen-mediated selection. A comparison of genome size of cyanobacteria in relatively stable marine versus more turbulent freshwater environments corroborates our predictions, although other factors between these environments could also be responsible.

Assessment of the bimodality in the distribution of bacterial genome sizes

Bacterial genome sizes have previously been shown to exhibit a bimodal distribution. This phenomenon has prompted discussion regarding the evolutionary forces driving genome size in bacteria and its ecological significance. We investigated the level of inherent redundancy in the public database and the effect it has on the shape of the apparent bimodal distribution. Our study reveals that there is a significant bias in the genome sequencing efforts towards a certain group of species, and that correcting the bias using species nomenclature and clustering of the 16S rRNA gene, results in a unimodal rather than the previously published bimodal distribution. The true genome size distribution and its wider ecological implications will soon emerge as we are currently witnessing rapid growth in the number of sequenced genomes from diverse environmental niches across a range of habitats at an unprecedented rate.

Absence of genome reduction in diverse, facultative endohyphal bacteria

Fungi interact closely with bacteria, both on the surfaces of the hyphae and within their living tissues (i.e. endohyphal bacteria, EHB). These EHB can be obligate or facultative symbionts and can mediate diverse phenotypic traits in their hosts. Although EHB have been observed in many lineages of fungi, it remains unclear how widespread and general these associations are, and whether there are unifying ecological and genomic features can be found across EHB strains as a whole. We cultured 11 bacterial strains after they emerged from the hyphae of diverse Ascomycota that were isolated as foliar endophytes of cupressaceous trees, and generated nearly complete genome sequences for all. Unlike the genomes of largely obligate EHB, the genomes of these facultative EHB resembled those of closely related strains isolated from environmental sources. Although all analysed genomes encoded structures that could be used to interact with eukaryotic hosts, pathways previously implicated in maintenance and establishment of EHB symbiosis were not universally present across all strains. Independent isolation of two nearly identical pairs of strains from different classes of fungi, coupled with recent experimental evidence, suggests horizontal transfer of EHB across endophytic hosts. Given the potential for EHB to influence fungal phenotypes, these genomes could shed light on the mechanisms of plant growth promotion or stress mitigation by fungal endophytes during the symbiotic phase, as well as degradation of plant material during the saprotrophic phase. As such, these findings contribute to the illumination of a new dimension of functional biodiversity in fungi.

Metagenomics and Bioinformatics in Microbial Ecology: Current Status and Beyond

Metagenomic approaches are now commonly used in microbial ecology to study microbial communities in more detail, including many strains that cannot be cultivated in the laboratory. Bioinformatic analyses make it possible to mine huge metagenomic datasets and discover general patterns that govern microbial ecosystems. However, the findings of typical metagenomic and bioinformatic analyses still do not completely describe the ecology and evolution of microbes in their environments. Most analyses still depend on straightforward sequence similarity searches against reference databases. We herein review the current state of metagenomics and bioinformatics in microbial ecology and discuss future directions for the field. New techniques will allow us to go beyond routine analyses and broaden our knowledge of microbial ecosystems. We need to enrich reference databases, promote platforms that enable meta- or comprehensive analyses of diverse metagenomic datasets, devise methods that utilize long-read sequence information, and develop more powerful bioinformatic methods to analyze data from diverse perspectives.

Functional profiling of cyanobacterial genomes and its role in ecological adaptations

With the availability of complete genome sequences of many cyanobacterial species, it is becoming feasible to study the broad prospective of the environmental adaptation and the overall changes at transcriptional and translational level in these organisms. In the evolutionary phase, niche-specific competitive forces have resulted in specific features of the cyanobacterial genomes. In this study, functional composition of the 84 different cyanobacterial genomes and their adaptations to different environments was examined by identifying the genomic composition for specific cellular processes, which reflect their genomic functional profile and ecological adaptation. It was identified that among cyanobacterial genomes, metabolic genes have major share over other categories and differentiation of genomic functional profile was observed for the species inhabiting different habitats. The cyanobacteria of freshwater and other habitats accumulate large number of poorly characterized genes. Strain specific functions were also reported in many cyanobacterial members, of which an important feature was the occurrence of phage-related sequences. From this study, it can be speculated that habitat is one of the major factors in giving the shape of functional composition of cyanobacterial genomes towards their ecological adaptations.