Life history implications of rRNA gene copy number in Escherichia coli

Illuminating ecology and distribution of the rare fungus Phellinidium pouzarii in the Bavarian Forest National Park

Due to their cryptic lifestyle, hidden diversity and a lack of ecological knowledge, conservation of wood-inhabiting fungi continues to be a niche interest. Molecular methods are able to provide deeper insights into the ecology of rare fungal species. We investigated the occurrence of the rare wood-inhabiting fungus Phellinidium pouzarii across the Bavarian Forest National Park in Germany using a fruit body survey, amplicon sequencing and qPCR. Additionally, we sequenced the genome of P. pouzarii and characterized the chemical substances responsible for its distinctive scent. Our approach gave matching results between amplicon sequencing and qPCRs, however, we found no evidence that P. pouzarii is more abundant in the National Park than we can assume based on fruit body inventories, underlining the species' critically endangered status. Genomics revealed P. pouzarii's repertoire of ligninolytic enzymes, pointing towards a white rot lifestyle. Two main components of P. pouzarii's distinct odour we identified (2-phenylethanol, methyl p-anisate) are known to act as insect attractants and/or to possess antimicrobial properties.

A sweeping view of avian mycoplasmas biology drawn from comparative genomic analyses

BACKGROUND

Avian mycoplasmas are small bacteria associated with several pathogenic conditions in many wild and poultry bird species. Extensive genomic data are available for many avian mycoplasmas, yet no comparative studies focusing on this group of mycoplasmas have been undertaken so far.

RESULTS

Here, based on the comparison of forty avian mycoplasma genomes belonging to ten different species, we provide insightful information on the phylogeny, pan/core genome, energetic metabolism, and virulence of these avian pathogens. Analyses disclosed considerable inter- and intra-species genomic variabilities, with genome sizes that can vary by twice as much. Phylogenetic analysis based on concatenated orthologous genes revealed that avian mycoplasmas fell into either Hominis or Pneumoniae groups within the Mollicutes and could split into various clusters. No host co-evolution of avian mycoplasmas can be inferred from the proposed phylogenetic scheme. With 3,237 different gene clusters, the avian mycoplasma group under study proved diverse enough to have an open pan genome. However, a set of 150 gene clusters was found to be shared between all avian mycoplasmas, which is likely encoding essential functions. Comparison of energy metabolism pathways showed that avian mycoplasmas rely on various sources of energy. Superposition between phylogenetic and energy metabolism groups revealed that the glycolytic mycoplasmas belong to two distinct phylogenetic groups (Hominis and Pneumoniae), while all the arginine-utilizing mycoplasmas belong only to Hominis group. This can stand for different evolutionary strategies followed by avian mycoplasmas and further emphasizes the diversity within this group. Virulence determinants survey showed that the involved gene arsenals vary significantly within and between species, and could even be found in species often reported apathogenic. Immunoglobulin-blocking proteins were detected in almost all avian mycoplasmas. Although these systems are not exclusive to this group, they seem to present some particular features making them unique among mycoplasmas.

CONCLUSION

This comparative genomic study uncovered the significant variable nature of avian mycoplasmas, furthering our knowledge on their biological attributes and evoking new hallmarks.

Genome and pan-genome analysis of a new exopolysaccharide-producing bacterium Pyschrobacillus sp. isolated from iron ores deposit and insights into iron uptake

Bacterial exopolysaccharides (EPS) have emerged as one of the key players in the field of heavy metal-contaminated environmental bioremediation. This study aimed to characterize and evaluate the metal biosorption potential of EPS produced by a novel Psychrobacillus strain, NEAU-3TGS, isolated from an iron ore deposit at Tamra iron mine, northern Tunisia. Genomic and pan-genomic analysis of NEAU-3TGS bacterium with nine validated published Psychrobacillus species was also performed. The results showed that the NEAU-3TGS genome (4.48 Mb) had a mean GC content of 36%, 4,243 coding sequences and 14 RNA genes. Phylogenomic analysis and calculation of nucleotide identity (ANI) values (less than 95% for new species with all strains) confirmed that NEAU-3TGS represents a potential new species. Pangenomic analysis revealed that Psychrobacillus genomic diversity represents an "open" pangenome model with 33,091 homologous genes, including 65 core, 3,738 shell, and 29,288 cloud genes. Structural EPS characterization by attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy showed uronic acid and α-1,4-glycosidic bonds as dominant components of the EPS. X-ray diffraction (XRD) analysis revealed the presence of chitin, chitosan, and calcite CaCO3 and confirmed the amorphous nature of the EPS. Heavy metal bioabsorption assessment showed that iron and lead were more adsorbed than copper and cadmium. Notably, the optimum activity was observed at 37°C, pH=7 and after 3 h contact of EPS with each metal. Genomic insights on iron acquisition and metabolism in Psychrobacillus sp. NEAU-3TGS suggested that no genes involved in siderophore biosynthesis were found, and only the gene cluster FeuABCD and trilactone hydrolase genes involved in the uptake of siderophores, iron transporter and exporter are present. Molecular modelling and docking of FeuA (protein peptidoglycan siderophore-binding protein) and siderophores ferrienterobactine [Fe+3 (ENT)]-3 and ferribacillibactine [Fe+3 (BB)]-3 ligand revealed that [Fe+3 (ENT)]-3 binds to Phe122, Lys127, Ile100, Gln314, Arg215, Arg217, and Gln252. Almost the same for [Fe+3 (ENT)]-3 in addition to Cys222 and Tyr229, but not Ile100.To the best of our knowledge, this is the first report on the characterization of EPS and the adsorption of heavy metals by Psychrobacillus species. The heavy metal removal capabilities may be advantageous for using these organisms in metal remediation.

The recovery of soil eukaryotic alpha and beta diversity after wetland restoration

Soil eukaryotes play an important role in regulating the ecological processes and ecosystem functioning. However, the recovery potential of soil eukaryotic diversity during wetland restoration is largely unknown. We compared the alpha and beta diversity of soil eukaryotes of farmlands and natural and restored wetlands to explore the underlying abiotic and biotic driving forces in the Sanjiang Plain, China. We found that there was no significant difference of the alpha diversity of soil eukaryotes, while the beta diversity of soil eukaryotes differed significantly between the three land use types, with the mean values in the restored wetlands in between those in the natural wetlands and farmlands. The composition of soil eukaryotic communities were less diverse in farmlands compared to restored and natural wetlands. Network property of soil eukaryotes community (positive: negative edges) increased from farmlands to restored wetlands to natural wetlands, indicating enhanced species positive: negative interactions during restoration. The structural equation modeling indicated that species positive: negative interactions and soil nutrients directly affected soil eukaryotic beta diversity. Soil pH and soil water content indirectly affected soil eukaryotic beta diversity by directly affecting species interactions. Our findings suggest that wetland restoration could change soil environment, strengthen microbial cooperation, and increase eukaryotic beta diversity. However, it may take a very long time to reach the original level of soil eukaryotic structure and diversity.

Lachnospiraceae are emerging industrial biocatalysts and biotherapeutics

The Lachnospiraceae is a family of anaerobic bacteria in the class Clostridia with potential to advance the bio-economy and intestinal therapeutics. Some species of Lachnospiraceae metabolize abundant, low-cost feedstocks such as lignocellulose and carbon dioxide into value-added chemicals. Others are among the dominant species of the human colon and animal rumen, where they ferment dietary fiber to promote healthy gut and immune function. Here, we summarize recent studies of the physiology, cultivation, and genetics of Lachnospiraceae, highlighting their wide substrate utilization and metabolic products with industrial applications. We examine studies of these bacteria as Live Biotherapeutic Products (LBPs), focusing on in vivo disease models and clinical studies using them to treat infection, inflammation, metabolic syndrome, and cancer. We discuss key research areas including elucidation of intra-specific diversity and genetic modification of candidate strains that will facilitate the exploitation of Lachnospiraceae in industry and medicine.

Land use alters bacterial growth dynamics in soil

Microbial growth and mortality are major determinants of soil carbon cycling. We measured in situ growth dynamics of individual bacterial taxa in cropped and successional soils in response to a resource pulse. We hypothesized that land use imposes selection pressures on growth characteristics. We estimated growth and death for 453 and 73 taxa, respectively. The average generation time was 5.04 ± 6.28 (SD; range 0.7-63.5) days. Lag times were shorter in cultivated than successional soils and resource amendment decreased lag times. Taxa exhibiting the greatest growth response also exhibited the greatest mortality, indicative of boom-and-bust dynamics. We observed a bimodal growth rate distribution, representing fast- and slow-growing clusters. Both clusters grew more rapidly in successional soils, which had more organic matter, than cultivated soils. Resource amendment increased the growth rate of the slower growing but not the faster-growing cluster via a mixture of increased growth rates and species turnover, indicating that competitive dynamics constrain growth rates in situ. These two clusters show that copiotrophic bacteria in soils may be subdivided into different life history groups and that these subgroups respond independently to land use and resource availability.

Cellular RNA levels define heterotrophic substrate-uptake rate sub-guilds in activated sludge microbial communities

A heterotrophic-specialist model was proposed previously to divide wastewater treatment plant (WWTP) heterotrophs into sub-guilds of consumers of readily or slowly degradable substrates (RDS or SDS, respectively). The substrate degradation rate model coupled to metabolic considerations predicted that RNA and polyhydroxyalkanoate (PHA) levels would be positively correlated in the activated sludge communities with high RNA and PHA occurring in RDS-consumers, and low RNA with no PHA accumulation occurring in SDS-consumers because their external substrates are always present. This prediction was verified in previous studies and in the current one. Thus, RNA and PHA levels were used as biomarkers of the RDS- and SDS-consumer sub-guilds for cell sorting using flow cytometry of samples from three WWTPs. Subsequently, 16S rRNA gene amplicon sequencing revealed that the sorted groups were highly similar over time and among WWTPs, and demonstrated a clear segregation by RNA levels. Predicted ecophysiological traits based on 16S rRNA phylogeny suggested that the high-RNA population showed RDS-consumer traits such as higher rrn copy numbers per genome. Using a mass-flow immigration model, it appeared that the high-RNA populations exhibited high immigration rates more frequently than low-RNA populations, but the differences in frequencies were less with increasing solids residence times.

Warmer temperatures favor slower-growing bacteria in natural marine communities

Earth's life-sustaining oceans harbor diverse bacterial communities that display varying composition across time and space. While particular patterns of variation have been linked to a range of factors, unifying rules are lacking, preventing the prediction of future changes. Here, analyzing the distribution of fast- and slow-growing bacteria in ocean datasets spanning seasons, latitude, and depth, we show that higher seawater temperatures universally favor slower-growing taxa, in agreement with theoretical predictions of how temperature-dependent growth rates differentially modulate the impact of mortality on species abundances. Changes in bacterial community structure promoted by temperature are independent of variations in nutrients along spatial and temporal gradients. Our results help explain why slow growers dominate at the ocean surface, during summer, and near the tropics and provide a framework to understand how bacterial communities will change in a warmer world.

The layered costs and benefits of translational redundancy

The rate and accuracy of translation hinges upon multiple components - including transfer RNA (tRNA) pools, tRNA modifying enzymes, and rRNA molecules - many of which are redundant in terms of gene copy number or function. It has been hypothesized that the redundancy evolves under selection, driven by its impacts on growth rate. However, we lack empirical measurements of the fitness costs and benefits of redundancy, and we have poor a understanding of how this redundancy is organized across components. We manipulated redundancy in multiple translation components of Escherichia coli by deleting 28 tRNA genes, 3 tRNA modifying systems, and 4 rRNA operons in various combinations. We find that redundancy in tRNA pools is beneficial when nutrients are plentiful and costly under nutrient limitation. This nutrient-dependent cost of redundant tRNA genes stems from upper limits to translation capacity and growth rate, and therefore varies as a function of the maximum growth rate attainable in a given nutrient niche. The loss of redundancy in rRNA genes and tRNA modifying enzymes had similar nutrient-dependent fitness consequences. Importantly, these effects are also contingent upon interactions across translation components, indicating a layered hierarchy from copy number of tRNA and rRNA genes to their expression and downstream processing. Overall, our results indicate both positive and negative selection on redundancy in translation components, depending on a species' evolutionary history with feasts and famines.

Ribosomal RNA operons define a central functional compartment in the Streptomyces chromosome

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. These bacteria possess a large linear chromosome genetically compartmentalized: core genes are grouped in the central part, while terminal regions are populated by poorly conserved genes. In exponentially growing cells, chromosome conformation capture unveiled sharp boundaries formed by ribosomal RNA (rrn) operons that segment the chromosome into multiple domains. Here we further explore the link between the genetic distribution of rrn operons and Streptomyces genetic compartmentalization. A large panel of genomes of species representative of the genus diversity revealed that rrn operons and core genes form a central skeleton, the former being identifiable from their core gene environment. We implemented a new nomenclature for Streptomyces genomes and trace their rrn-based evolutionary history. Remarkably, rrn operons are close to pericentric inversions. Moreover, the central compartment delimited by rrn operons has a very dense, nearly invariant core gene content. Finally, this compartment harbors genes with the highest expression levels, regardless of gene persistence and distance to the origin of replication. Our results highlight that rrn operons are structural boundaries of a central functional compartment prone to transcription in Streptomyces.

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.

Organelle 16S rRNA amplicon sequencing enables profiling of active gut microbiota in murine model

High-throughput sequencing of ribosomal RNA (rRNA) amplicons has served as a cornerstone in microbiome studies. Despite crucial implication of organelle 16S rRNA measurements to host gut microbial activities, genomic DNA (gDNA) was overwhelmingly targeted for amplicon sequencings. Although gDNA could be a reliable resource for gene existing validation, little information is revealed in regard to the activity of microorganisms owing to the limited changes gDNA undertaken in inactive, dormant, and dead bacteria. We applied both rRNA- and gDNA-derived sequencings on mouse cecal contents. Respective experimental designs were verified to be suitable for nucleic acid (NA) purification. Via benchmarking, mainstream 16S rRNA hypervariable region targets and reference databases were proven adequate for respective amplicon sequencing study. In phylogenetic studies, significant microbial composition differences were observed between two methods. Desulfovibrio spp. (an important group of anaerobic gut microorganisms that has caused analytical difficulties), Pediococcus spp., and Proteobacteria were drastically lower as represented by gDNA-derived compositions, while microbes like Firmicutes were higher as represented by gDNA-derived microbiome compositions. Also, using PICRUSt2 as an example, we illustrated that rRNA-derived sequencing might be more suitable for microbiome function predictions since pathways like sugar metabolism were lower as represented by rRNA-derived results. The findings of this study demonstrated that rRNA-derived amplicon sequencing could improve identification capability of specific gut microorganisms and might be more suitable for in silico microbiome function predictions. Therefore, rRNA-derived amplicon sequencings, preferably coupled with gDNA-derived ones, could be used as a capable tool to unveil active microbial components in host gut. KEY POINTS: • Conventional pipelines were adequate for the respective amplicon sequencing study • Groups, such as Desulfovibrio spp., were differently represented by two methods • Comparative amplicon sequencings could be useful in host active microbiota studies.

Comparative Genomic Analysis of the Human Pathogen Wohlfahrtiimonas Chitiniclastica Provides Insight Into the Identification of Antimicrobial Resistance Genotypes and Potential Virulence Traits

Recent studies suggest that Wohlfahrtiimonas chitiniclastica may be the cause of several diseases in humans including sepsis and bacteremia making the bacterium as a previously underappreciated human pathogen. However, very little is known about the pathogenicity and genetic potential of W. chitiniclastica; therefore, it is necessary to conduct systematic studies to gain a deeper understanding of its virulence characteristics and treatment options. In this study, the entire genetic repertoire of all publicly available W. chitiniclastica genomes was examined including in silico characterization of bacteriophage content, antibiotic resistome, and putative virulence profile. The pan-genome of W. chitiniclastica comprises 3819 genes with 1622 core genes (43%) indicating a putative metabolic conserved species. Furthermore, in silico analysis indicated presumed resistome expansion as defined by the presence of genome-encoded transposons and bacteriophages. While macrolide resistance genes macA and macB are located within the core genome, additional antimicrobial resistance genotypes for tetracycline (tetH, tetB, and tetD), aminoglycosides (ant(2'')-Ia, aac(6')-Ia,aph(3'')-Ib, aph(3')-Ia, and aph(6)-Id)), sulfonamide (sul2), streptomycin (strA), chloramphenicol (cat3), and beta-lactamase (blaVEB) are distributed among the accessory genome. Notably, our data indicate that the type strain DSM 18708T does not encode any additional clinically relevant antibiotic resistance genes, whereas drug resistance is increasing within the W. chitiniclastica clade. This trend should be monitored with caution. To the best of our knowledge, this is the first comprehensive genome analysis of this species, providing new insights into the genome of this opportunistic human pathogen.

Coastal Microbial Communities Disrupted During the 2018 Hurricane Season in Outer Banks, North Carolina

Hurricane frequencies and intensities are expected to increase under warming climate scenarios, increasing potential to disrupt microbial communities from steady-state conditions and alter ecosystem function. This study shows the impact of hurricane season on microbial community dynamics within the barrier island system of Outer Banks, North Carolina. We found that the passage of two sequential energetic hurricanes in 2018 (Florence and Michael) were correlated with shifts in total and active (DNA and RNA) portions of bacterial communities but not in archaeal communities, and within surface waters but not within the sediment. These microbial community shifts were distinct from non-hurricane season conditions, suggesting significant implications for nutrient cycling in nearshore and offshore environments. Hurricane-influenced marine sites in the coastal North Atlantic region had lower microbial community evenness and Shannon diversity, in addition to increased relative abundance of copiotrophic microbes compared to non-hurricane conditions. The abundance of functional genes associated with carbon and nitrogen cycling pathways were also correlated with the storm season, potentially shifting microbial communities at offshore sites from autotroph-dominated to heterotroph-dominated and leading to impacts on local carbon budgets. Understanding the geographic- and system-dependent responses of coastal microbial communities to extreme storm disturbances is critical for predicting impacts to nutrient cycling and ecosystem stability in current and future climate scenarios.

New insights into the r/K selection theory achieved in methanogenic systems through continuous-flow and sequencing batch operational modes

Anaerobic digestion is achieved through cooperation among various types of microorganisms, and the regulation of microbial communities is key to achieving stable system operation. In this study, the r/K selection theory was adopted to examine the system performance and microbial characteristics in anaerobic reactors with different operating modes (continuous-flow reactors, CFRs; sequencing batch reactors, SBRs) and sludge retention times (25 and 10 days). Four lab-scale reactors (CFR_25d, CFR_10d, SBR_25d, and SBR_10d) were operated. In the cycle reaction, CFR_25d achieved the highest methane yield (678.0 mL/L) and methane production rate (140.8 mL/(L·h)); while those in CFR_10d were the lowest, which could have been due to an accumulation of volatile fatty acids. CFR could wash out r-strategists efficiently, such as Methanosarcina. CFR_25d and CFR_10d significantly enriched the K-strategist Geobacter, with the relative abundances of 34.0% and 72.6%, respectively. In addition, the hydrogenotrophic methanogens of Methanolinea and Methanospirillum (K-strategists) dominated in CFR_25d and CFR_10d. Methanobacterium adapted to the diverse operational conditions, but the slow grower Methanosaeta only accounted for 0.9% in CFR_10d. Failure to enrich propionate oxidizers resulted in a functional absence of propionate degradation in the CFRs.

Considerations on the Identity and Diversity of Organisms Affiliated with Sphingobacterium multivorum-Proposal for a New Species, Sphingobacterium paramultivorum

Plant biomass offers great potential as a sustainable resource, and microbial consortia are primordial in its bioconversion. The wheat-straw-biodegradative bacterial strain w15 has drawn much attention as a result of its biodegradative potential and superior degradation performance in bacterial-fungal consortia. Strain w15 was originally assigned to the species Sphingobacterium multivorum based on its 16S ribosomal RNA (rRNA) gene sequence. A closer examination of this taxonomic placement revealed that the sequence used has 98.9% identity with the ‘divergent’ 16S rRNA gene sequence of S. multivorum NCTC 11343^T, yet lower relatedness with the canonical 16S rRNA sequence. A specific region of the gene, located between positions 186 and 210, was found to be highly variable and determinative for the divergence. To solve the identity of strain w15, genome metrics and analyses of ecophysiological niches were undertaken on a selection of strains assigned to S. multivorum and related species. These analyses separated all strains into three clusters, with strain w15, together with strain BIGb0170, constituting a separate radiation, next to S. multivorum and S. siyangense. Moreover, the strains denoted FDAARGOS 1141 and 1142 were placed inside S. siyangense. We propose the renaming of strains w15 and BIGb0170 as members of the novel species, coined Sphingobacterium paramultivorum.

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits-mode of genomic inheritance-influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S-5.8S-28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA-histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction-with potential functional consequences related to the loss of sexual reproduction.

Non-Conventional Yeasts as Alternatives in Modern Baking for Improved Performance and Aroma Enhancement

Saccharomyces cerevisiae remains the baker’s yeast of choice in the baking industry. However, its ability to ferment cereal flour sugars and accumulate CO2 as a principal role of yeast in baking is not as unique as previously thought decades ago. The widely conserved fermentative lifestyle among the Saccharomycotina has increased our interest in the search for non-conventional yeast strains to either augment conventional baker’s yeast or develop robust strains to cater for the now diverse consumer-driven markets. A decade of research on alternative baker’s yeasts has shown that non-conventional yeasts are increasingly becoming important due to their wide carbon fermentation ranges, their novel aromatic flavour generation, and their robust stress tolerance. This review presents the credentials of non-conventional yeasts as attractive yeasts for modern baking. The evolution of the fermentative trait and tolerance to baking-associated stresses as two important attributes of baker’s yeast are discussed besides their contribution to aroma enhancement. The review further discusses the approaches to obtain new strains suitable for baking applications.

Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community

Temperature sensitivity (Q₁₀ ) of soil organic matter (SOM) decomposition is a crucial parameter to predict the fate of soil carbon (C) under global warming. Nonetheless, the response pattern of Q₁₀ to continuous warming and the underlying mechanisms are still under debate, especially considering the complex interactions between Q₁₀ , SOM quality, and soil microorganisms. We examined the Q₁₀ of SOM decomposition across a mean annual temperature (MAT) gradient from -1.9 to 5.1°C in temperate mixed forest ecosystems in parallel with SOM quality and bioavailability, microbial taxonomic composition, and functional genes responsible for organic carbon decomposition. Within this temperature gradient of 7.0°C, the Q₁₀ values increased with MAT, but decreased with SOM bioavailability. The Q₁₀ values increased with the prevalence of K-strategy of soil microbial community, which was characterized by: (i) high ratios of oligotrophic to copiotrophic taxa, (ii) ectomycorrhizal to saprotrophic fungi, (iii) functional genes responsible for degradation of recalcitrant to that of labile C, and (iv) low average 16S rRNA operon copy number. Because the recalcitrant organic matter was mainly utilized by the K-strategists, these findings independently support the carbon quality-temperature theory from the perspective of microbial taxonomic composition and functions. A year-long incubation experiment was performed to determine the response of labile and recalcitrant C pools to warming based on the two-pool model. The decomposition of recalcitrant SOM was more sensitive to increased temperature in southern warm regions, which might attribute to the dominance of K-selected microbial communities. It implies that climate warming would mobilize the larger recalcitrant pools in warm regions, exacerbating the positive feedback between increased MAT and CO₂ efflux. This is the first attempt to link temperature sensitivity of SOM decomposition with microbial eco-strategies by incorporating the genetic information and disentangling the complex relationship between Q₁₀ and soil microorganisms.

Rates of gene conversions between Escherichia coli ribosomal operons

Due to their universal presence and high sequence conservation, ribosomal RNA (rRNA) sequences are used widely in phylogenetics for inferring evolutionary relationships between microbes and in metagenomics for analyzing the composition of microbial communities. Most microbial genomes encode multiple copies of rRNA genes to supply cells with sufficient capacity for protein synthesis. These copies typically undergo concerted evolution that keeps their sequences identical, or nearly so, due to gene conversion, a type of intragenomic recombination that changes one copy of a homologous sequence to exactly match another. Widely varying rates of rRNA gene conversion have previously been estimated by comparative genomics methods and using genetic reporter assays. To more directly measure rates of rRNA intragenomic recombination, we sequenced the seven Escherichia coli rRNA operons in 15 lineages that were evolved for ∼13,750 generations with frequent single-cell bottlenecks that reduce the effects of selection. We identified 38 gene conversion events and estimated an overall rate of intragenomic recombination within the 16S and 23S genes between rRNA copies of 3.6 × 10⁻⁴ per genome per generation or 8.6 × 10⁻⁶ per rRNA operon per homologous donor operon per generation. This rate varied only slightly from random expectations at different sites within the rRNA genes and between rRNA operons located at different positions in the genome. Our accurate estimate of the rate of rRNA gene conversions fills a gap in our quantitative understanding of how ribosomal sequences and other multicopy elements diversify and homogenize during microbial genome evolution.

Plasmid localization of sole rrn operon in genomes of Oecophyllibacter saccharovorans (Acetobacteraceae)

The bacterium Oecophyllibacter saccharovorans of family Acetobacteraceae is a symbiont of weaver ant Oecophylla smaragdina. In our previous study, we published the finding of novel O. saccharovorans strains Ha5^T, Ta1 and Jb2 (Chua et al. 2020) but their plasmid sequences have not been reported before. Here, we demonstrate for the first time that the sole rrn operon of their genomes was detected on a 6.6 kb circular replicon. This replicon occurred in high copy number, much smaller size and lower G + C content than the main chromosome. Based on these features, the 6.6 kb circular replicon was regarded as rrn operon-containing plasmid. Further restriction analysis on the plasmids confirmed their circular conformation. A Southern hybridization analysis also corroborated the presence of 16S rRNA gene and thus the rrn operon on a single locus in the genome of the O. saccharovorans strains. However, similar genome architecture was not observed in other closely related bacterial strains. Additional survey also detected no plasmid-borne rrn operon in available genomes of validly described taxa of family Acetobacteraceae. To date, plasmid localization of rrn operon is rarely documented. This study reports the occurrence of rrn operon on the smallest bacterial plasmid in three O. saccharovorans strains and discusses its possible importance in enhancing their competitive fitness as bacterial symbiont of O. smaragdina.

Symbiotic and Nonsymbiotic Members of the Genus Ensifer (syn. Sinorhizobium) Are Separated into Two Clades Based on Comparative Genomics and High-Throughput Phenotyping

Rhizobium–legume symbioses serve as paradigmatic examples for the study of mutualism evolution. The genus Ensifer (syn. Sinorhizobium) contains diverse plant-associated bacteria, a subset of which can fix nitrogen in symbiosis with legumes. To gain insights into the evolution of symbiotic nitrogen fixation (SNF), and interkingdom mutualisms more generally, we performed extensive phenotypic, genomic, and phylogenetic analyses of the genus Ensifer. The data suggest that SNF likely emerged several times within the genus Ensifer through independent horizontal gene transfer events. Yet, the majority (105 of 106) of the Ensifer strains with the nodABC and nifHDK nodulation and nitrogen fixation genes were found within a single, monophyletic clade. Comparative genomics highlighted several differences between the “symbiotic” and “nonsymbiotic” clades, including divergences in their pangenome content. Additionally, strains of the symbiotic clade carried 325 fewer genes, on average, and appeared to have fewer rRNA operons than strains of the nonsymbiotic clade. Initial characterization of a subset of ten Ensifer strains identified several putative phenotypic differences between the clades. Tested strains of the nonsymbiotic clade could catabolize 25% more carbon sources, on average, than strains of the symbiotic clade, and they were better able to grow in LB medium and tolerate alkaline conditions. On the other hand, the tested strains of the symbiotic clade were better able to tolerate heat stress and acidic conditions. We suggest that these data support the division of the genus Ensifer into two main subgroups, as well as the hypothesis that pre-existing genetic features are required to facilitate the evolution of SNF in bacteria.

Does Intraspecific Variation in rDNA Copy Number Affect Analysis of Microbial Communities?

Amplicon sequencing of partial regions of the ribosomal RNA loci (rDNA) is widely used to profile microbial communities. However, the rDNA is dynamic and can exhibit substantial interspecific and intraspecific variation in copy number in prokaryotes and, especially, in microbial eukaryotes. As change in rDNA copy number is a common response to environmental change, rDNA copy number is not necessarily a property of a species. Variation in rDNA copy number, especially the capacity for large intraspecific changes driven by external cues, complicates analyses of rDNA amplicon sequence data. We highlight the need to (i) interpret amplicon sequence data in light of possible interspecific and intraspecific variation, and (ii) examine the potential plasticity in rDNA copy number as an important ecological factor to better understand how microbial communities are structured in heterogeneous environments.

Long-Term Nutrient Enrichment of an Oligotroph-Dominated Wetland Increases Bacterial Diversity in Bulk Soils and Plant Rhizospheres

In nutrient-limited conditions, plants rely on rhizosphere microbial members to facilitate nutrient acquisition, and in return, plants provide carbon resources to these root-associated microorganisms. However, atmospheric nutrient deposition can affect plant-microbe relationships by changing soil bacterial composition and by reducing cooperation between microbial taxa and plants. To examine how long-term nutrient addition shapes rhizosphere community composition, we compared traits associated with bacterial (fast-growing copiotrophs, slow-growing oligotrophs) and plant (C3 forb, C4 grass) communities residing in a nutrient-poor wetland ecosystem. Results revealed that oligotrophic taxa dominated soil bacterial communities and that fertilization increased the presence of oligotrophs in bulk and rhizosphere communities. Additionally, bacterial species diversity was greatest in fertilized soils, particularly in bulk soils. Nutrient enrichment (fertilized versus unfertilized) and plant association (bulk versus rhizosphere) determined bacterial community composition; bacterial community structure associated with plant functional group (grass versus forb) was similar within treatments but differed between fertilization treatments. The core forb microbiome consisted of 602 unique taxa, and the core grass microbiome consisted of 372 unique taxa. Forb rhizospheres were enriched in potentially disease-suppressive bacterial taxa, and grass rhizospheres were enriched in bacterial taxa associated with complex carbon decomposition. Results from this study demonstrate that fertilization serves as a strong environmental filter on the soil microbiome, which leads to distinct rhizosphere communities and can shift plant effects on the rhizosphere microbiome. These taxonomic shifts within plant rhizospheres could have implications for plant health and ecosystem functions associated with carbon and nitrogen cycling.IMPORTANCE Over the last century, humans have substantially altered nitrogen and phosphorus cycling. Use of synthetic fertilizer and burning of fossil fuels and biomass have increased nitrogen and phosphorus deposition, which results in unintended fertilization of historically low-nutrient ecosystems. With increased nutrient availability, plant biodiversity is expected to decline, and the abundance of copiotrophic taxa is anticipated to increase in bacterial communities. Here, we address how bacterial communities associated with different plant functional types (forb, grass) shift due to long-term nutrient enrichment. Unlike other studies, results revealed an increase in bacterial diversity, particularly of oligotrophic bacteria in fertilized plots. We observed that nutrient addition strongly determines forb and grass rhizosphere composition, which could indicate different metabolic preferences in the bacterial communities. This study highlights how long-term fertilization of oligotroph-dominated wetlands could alter diversity and metabolism of rhizosphere bacterial communities in unexpected ways.

A chromosome-scale assembly of the smallest Dothideomycete genome reveals a unique genome compaction mechanism in filamentous fungi

Background

The wide variation in the size of fungal genomes is well known, but the reasons for this size variation are less certain. Here, we present a chromosome-scale assembly of ectophytic Peltaster fructicola, a surface-dwelling extremophile, based on long-read DNA sequencing technology, to assess possible mechanisms associated with genome compaction.

Results

At 18.99 million bases (Mb), P. fructicola possesses one of the smallest known genomes sequence among filamentous fungi. The genome is highly compact relative to other fungi, with substantial reductions in repeat content, ribosomal DNA copies, tRNA gene quantity, and intron sizes, as well as intergenic lengths and the size of gene families. Transposons take up just 0.05% of the entire genome, and no full-length transposon was found. We concluded that reduced genome sizes in filamentous fungi such as P. fructicola, Taphrina deformans and Pneumocystis jirovecii occurred through reduction in ribosomal DNA copy number and reduced intron sizes. These dual mechanisms contrast with genome reduction in the yeast fungus Saccharomyces cerevisiae, whose small and compact genome is associated solely with intron loss.

Conclusions

Our results reveal a unique genomic compaction architecture of filamentous fungi inhabiting plant surfaces, and broaden the understanding of the mechanisms associated with compaction of fungal genomes.

Genomic adaptations in information processing underpin trophic strategy in a whole-ecosystem nutrient enrichment experiment

Several universal genomic traits affect trade-offs in the capacity, cost, and efficiency of the biochemical information processing that underpins metabolism and reproduction. We analyzed the role of these traits in mediating the responses of a planktonic microbial community to nutrient enrichment in an oligotrophic, phosphorus-deficient pond in Cuatro Ciénegas, Mexico. This is one of the first whole-ecosystem experiments to involve replicated metagenomic assessment. Mean bacterial genome size, GC content, total number of tRNA genes, total number of rRNA genes, and codon usage bias in ribosomal protein sequences were all higher in the fertilized treatment, as predicted on the basis of the assumption that oligotrophy favors lower information-processing costs whereas copiotrophy favors higher processing rates. Contrasting changes in trait variances also suggested differences between traits in mediating assembly under copiotrophic versus oligotrophic conditions. Trade-offs in information-processing traits are apparently sufficiently pronounced to play a role in community assembly because the major components of metabolism-information, energy, and nutrient requirements-are fine-tuned to an organism's growth and trophic strategy.

Carbon Use Efficiency and Its Temperature Sensitivity Covary in Soil Bacteria

The strategy that microbial decomposers take with respect to using substrate for growth versus maintenance is one essential biological determinant of the propensity of carbon to remain in soil. To quantify the environmental sensitivity of this key physiological trade-off, we characterized the carbon use efficiency (CUE) of 23 soil bacterial isolates across seven phyla at three temperatures and with up to four substrates. Temperature altered CUE in both an isolate-specific manner and a substrate-specific manner. We searched for genes correlated with the temperature sensitivity of CUE on glucose and deemed those functional genes which were similarly correlated with CUE on other substrates to be validated as markers of CUE. Ultimately, we did not identify any such robust functional gene markers of CUE or its temperature sensitivity. However, we found a positive correlation between rRNA operon copy number and CUE, opposite what was expected. We also found that inefficient taxa increased their CUE with temperature, while those with high CUE showed a decrease in CUE with temperature. Together, our results indicate that CUE is a flexible parameter within bacterial taxa and that the temperature sensitivity of CUE is better explained by observed physiology than by genomic composition across diverse taxa. We conclude that the bacterial CUE response to temperature and substrate is more variable than previously thought.IMPORTANCE Soil microbes respond to environmental change by altering how they allocate carbon to growth versus respiration-or carbon use efficiency (CUE). Ecosystem and Earth System models, used to project how global soil C stocks will continue to respond to the climate crisis, often assume that microbes respond homogeneously to changes in the environment. In this study, we quantified how CUE varies with changes in temperature and substrate quality in soil bacteria and evaluated why CUE characteristics may differ between bacterial isolates and in response to altered growth conditions. We found that bacterial taxa capable of rapid growth were more efficient than those limited to slow growth and that taxa with high CUE were more likely to become less efficient at higher temperatures than those that were less efficient to begin with. Together, our results support the idea that the CUE temperature response is constrained by both growth rate and CUE and that this partly explains how bacteria acclimate to a warming world.

Time Series Genome-Centric Analysis Unveils Bacterial Response to Operational Disturbance in Activated Sludge

Understanding ecosystem response to disturbances and identifying the most critical traits for the maintenance of ecosystem functioning are important goals for microbial community ecology. In this study, we used 16S rRNA amplicon sequencing and metagenomics to investigate the assembly of bacterial populations in a full-scale municipal activated sludge wastewater treatment plant over a period of 3 years, including a 9-month period of disturbance characterized by short-term plant shutdowns. Following the reconstruction of 173 metagenome-assembled genomes, we assessed the functional potential, the number of rRNA gene operons, and the in situ growth rate of microorganisms present throughout the time series. Operational disturbances caused a significant decrease in bacteria with a single copy of the rRNA (rrn) operon. Despite moderate differences in resource availability, replication rates were distributed uniformly throughout time, with no differences between disturbed and stable periods. We suggest that the length of the growth lag phase, rather than the growth rate, is the primary driver of selection under disturbed conditions. Thus, the system could maintain its function in the face of disturbance by recruiting bacteria with the capacity to rapidly resume growth under unsteady operating conditions.IMPORTANCE Disturbance is a key determinant of community assembly and dynamics in natural and engineered ecosystems. Microbiome response to disturbance is thought to be influenced by bacterial growth traits and life history strategies. In this time series observational study, the response to disturbance of microbial communities in a full-scale activated sludge wastewater treatment plant was assessed by computing specific cellular traits of genomes retrieved from metagenomes. It was found that the genomes observed in disturbed periods have more copies of the rRNA operon than genomes observed in stable periods, whereas the in situ mean relative growth rates of bacteria present during stable and disturbed periods were indistinguishable. From these intriguing observations, we infer that the length of the lag phase might be a growth trait that affects the microbial response to disturbance. Further exploration of this hypothesis could contribute to better understanding of the adaptive response of microbiomes to unsteady environmental conditions.

A 16S rDNA PCR-based theoretical to actual delta approach on culturable mock communities revealed severe losses of diversity information

BACKGROUND

Subunits of ribosomal RNA genes (rDNAs) characterized by PCR-based protocols have been the proxy for studies in microbial taxonomy, phylogenetics, evolution and ecology. However, relevant factors have shown to interfere in the experimental outputs in a variety of systems. In this work, a 'theoretical' to 'actual' delta approach was applied to data on culturable mock bacterial communities (MBCs) to study the levels of losses in operational taxonomic units (OTUs) detectability. Computational and lab-bench strategies based on 16S rDNA amplification by 799F and U1492R primers were employed, using a fingerprinting method with highly improved detectability of fragments as a case-study tool. MBCs were of two major types: in silico MBCs, assembled with database-retrieved sequences, and in vitro MBCs, with AluI digestions of PCR data generated from culturable endophytes isolated from cacao trees.

RESULTS

Interfering factors for the 16 s rDNA amplifications, such as the type of template, direct and nested PCR, proportion of chloroplast DNA from a tropical plant source (Virola officinalis), and biased-amplification by the primers resulted in altered bacterial 16S rDNA amplification, both on MBCs and V. officinalis leaf-extracted DNA. For the theoretical data, the maximum number of fragments for in silico and in vitro cuts were not significantly different from each other. Primers' preferences for certain sequences were detected, depending on the MBCs' composition prior to PCR. The results indicated overall losses from 2.3 up to 8.2 times in the number of OTUs detected from actual AluI digestions of MBCs when compared to in silico and in vitro theoretical data.

CONCLUSIONS

Due to all those effects, the final amplification profile of the bacterial community assembled was remarkably simplified when compared to the expected number of detectable fragments known to be present in the MBC. From these findings, the scope of hypotheses generation and conclusions from experiments based on PCR amplifications of bacterial communities was discussed.

Deletion of rRNA Operons of Sinorhizobium fredii Strain NGR234 and Impact on Symbiosis With Legumes

During their lifecycle, from free-living soil bacteria to endosymbiotic nitrogen-fixing bacteroids of legumes, rhizobia must colonize, and cope with environments where nutrient concentrations and compositions vary greatly. Bacterial colonization of legume rhizospheres and of root surfaces is subject to a fierce competition for plant exudates. By contrast root nodules offer to rhizobia sheltered nutrient-rich environments within which the cells that successfully propagated via infection threads can rapidly multiply. To explore the effects on symbiosis of a slower rhizobia growth and metabolism, we deleted one or two copies of the three functional rRNA operons of the promiscuous Sinorhizobium fredii strain NGR234 and examined the impact of these mutations on free-living and symbiotic lifestyles. Strains with two functional rRNA operons (NGRΔrRNA1 and NGRΔrRNA3) grew almost as rapidly as NGR234, and NGRΔrRNA1 was as proficient as the parent strain on all of the five legume species tested. By contrast, the NGRΔrRNA1,3 double mutant, which carried a single rRNA operon and grew significantly slower than NGR234, had a reduced symbiotic proficiency on Cajanus cajan, Macroptilium atropurpureum, Tephrosia vogelii, and Vigna unguiculata. In addition, while NGRΔrRNA1 and NGR234 equally competed for nodulation of V. unguiculata, strain NGRΔrRNA1,3 was clearly outcompeted by wild-type. Surprisingly, on Leucaena leucocephala, NGRΔrRNA1,3 was the most proficient strain and competed equally NGR234 for nodule occupation. Together, these results indicate that for strains with otherwise identical repertoires of symbiotic genes, a faster growth on roots and/or inside plant tissues may contribute to secure access to nodules of some hosts. By contrast, other legumes such as L. leucocephala appear as less selective and capable of providing symbiotic environments susceptible to accommodate strains with a broader spectrum of competences.

Integrated evolution of ribosomal RNAs, introns, and intron nurseries

The initial components of ribosomes first appeared more than 3.8 billion years ago during a time when many types of RNAs were evolving. While modern ribosomes are complex molecular machines consisting of rRNAs and proteins, they were assembled during early evolution by the association and joining of small functional RNA units. Introns may have provided the means to ligate many of these pieces together. All four classes of introns (group I, group II, spliceosomal, and archaeal) are present in many rRNA gene loci over a broad phylogenetic range. A survey of rRNA intron sequences across the three major life domains suggests that some of the classes of introns may have diverged from one another within rRNA gene loci. Analyses of rRNA sequences revealed self-splicing group I and group II introns are present in ancestral regions of the SSU (small subunit) and LSU (large subunit), whereas spliceosomal and archaeal introns appeared in sections of the rRNA that evolved later. Most classes of introns increased in number for approximately 1 billion years. However, their frequencies are low in the most recently evolved regions added to the SSU and LSU rRNAs. Furthermore, many of the introns appear to have been in the same locations for billions of years, suggesting an ancient origin for these sequences. In this Perspectives paper, I reviewed and analyzed rRNA intron sequences, locations, structural characteristics, and splicing mechanisms; and suggest that rRNA gene loci may have served as evolutionary nurseries for intron formation and diversification.

Translational Selection for Speed Is Not Sufficient to Explain Variation in Bacterial Codon Usage Bias

Increasing growth rate across bacteria strengthens selection for faster translation, concomitantly increasing the total number of tRNA genes and codon usage bias (CUB: enrichment of specific synonymous codons in highly expressed genes). Typically, enriched codons are translated by tRNAs with higher gene copy numbers (GCN). A model of tRNA–CUB coevolution based on fast growth-associated selection on translational speed recapitulates these patterns. A key untested implication of the coevolution model is that translational selection should favor higher tRNA GCN for more frequently used amino acids, potentially weakening the effect of growth-associated selection on CUB. Surprisingly, we find that CUB saturates with increasing growth rate across γ-proteobacteria, even as the number of tRNA genes continues to increase. As predicted, amino acid-specific tRNA GCN is positively correlated with the usage of corresponding amino acids, but there is no correlation between growth rate associated changes in CUB and amino acid usage. Instead, we find that some amino acids—cysteine and those in the NNA/G codon family—show weak CUB that does not increase with growth rate, despite large variation in the corresponding tRNA GCN. We suggest that amino acid-specific variation in CUB is not explained by tRNA GCN because GCN does not influence the difference between translation times of synonymous codons as expected. Thus, selection on translational speed alone cannot fully explain quantitative variation in overall or amino acid-specific CUB, suggesting a significant role for other functional constraints and amino acid-specific codon features.

Multiple Ribosomal RNA Operons in Bacteria; Their Concerted Evolution and Potential Consequences on the Rate of Evolution of Their 16S rRNA

Bacterial species differ greatly in the number and location of the rRNA operons which may be present in the bacterial chromosomes and plasmids. Most bacterial species contain more than one ribosomal RNA operon copy in their genomes, with some species containing up to 15 such copies. We review the number and location of the rRNA operons and discuss evolution of 16S rRNA (rrs) genes -which are considered as ultimate chronometers for phylogenetic classification- in bacteria with multiple copies of these genes. In these bacterial species, the rrs genes must evolve in concert and sequence changes generated by mutation or horizontal gene transfer must be either erased or spread to every gene copy to avoid divergence, as it occurs when they are present in different species. Analysis of polymorphic sites in intra-genomic rrs copies identifies putative conversion events and demonstrates that sequence conversion is patchy and occurs in small conversion tracts. Sequence conversion probably arises by a non-reciprocal transfer between two or more copies where one copy contributes only a small contiguous segment of DNA, whereas the other copy contributes the rest of the genome in a fairly well understood molecular process. Because concerted evolution implies that a mutation in any of the rrs copies is either eliminated or transferred to every rrs gene in the genome, this process should slow their evolution rate relative to that of single copy genes. However, available data on the rrs genes in bacterial genomes do not show a clear relationship between their evolution rates and the number of their copies in the genome.

Long-term use of cover crops and no-till shift soil microbial community life strategies in agricultural soil

Reducing tillage and growing cover crops, widely recommended practices for boosting soil health, have major impacts on soil communities. Surprisingly little is known about their impacts on soil microbial functional diversity, and especially so in irrigated Mediterranean ecosystems. In long-term experimental plots at the West Side Research and Extension Center in California's Central Valley, we characterized soil microbial communities in the presence or absence of physical disturbance due to tillage, in the presence or absence of cover crops, and at three depths: 0-5, 5-15 and 15-30 cm. This characterization included qPCR for bacterial and archaeal abundances, DNA sequencing of the 16S rRNA gene, and phylogenetic estimation of two ecologically important microbial traits (rRNA gene copy number and genome size). Total (bacterial + archaeal) diversity was higher in no-till than standard till; diversity increased with depth in no-till but decreased with depth in standard till. Total bacterial numbers were higher in cover cropped plots at all depths, while no-till treatments showed higher numbers in 0-5 cm but lower numbers at lower depths compared to standard tillage. Trait estimates suggested that different farming practices and depths favored distinctly different microbial life strategies. Tillage in the absence of cover crops shifted microbial communities towards fast growing competitors, while no-till shifted them toward slow growing stress tolerators. Across all treatment combinations, increasing depth resulted in a shift towards stress tolerators. Cover crops shifted the communities towards ruderals-organisms with wider metabolic capacities and moderate rates of growth. Overall, our results are consistent with decreasing nutrient availability with soil depth and under no-till treatments, bursts of nutrient availability and niche homogenization under standard tillage, and increases in C supply and variety provided by cover crops. Understanding how agricultural practices shift microbial abundance, diversity and life strategies, such as presented here, can assist with designing farming systems that can support high yields, while enhancing C sequestration and increasing resilience to climate change.

On the intrinsic constraint of bacterial growth rate: M. tuberculosis's view of the protein translation capacity

In nature, the maximal growth rates vary widely among different bacteria species. Fast-growing bacteria species such as Escherichia coli can have a shortest generation time of 20 min. Slow-growing bacteria species are perhaps best known for Mycobacterium tuberculosis, a human pathogen with a generation time being no less than 16 h. Despite of the significant progress made on understanding the pathogenesis of M. tuberculosis, we know little on the origin of its intriguingly slow growth. From a global view, the intrinsic constraint of the maximal growth rate of bacteria remains to be a fundamental question in microbiology. In this review, we analyze and discuss this issue from the angle of protein translation capacity, which is the major demand for cell growth. Based on quantitative analysis, we propose four parameters: rRNA chain elongation rate, abundance of RNA polymerase engaged in rRNA synthesis, polypeptide chain elongation rate, and active ribosome fraction, which potentially limit the maximal growth rate of bacteria. We further discuss the relation of these parameters with the growth rate for M. tuberculosis as well as other bacterial species. We highlight future comprehensive investigation of these parameters for different bacteria species to understand how bacteria set their own specific growth rates.

Diversity of Translation Initiation Mechanisms across Bacterial Species Is Driven by Environmental Conditions and Growth Demands

The Shine-Dalgarno (SD) sequence motif is frequently found upstream of protein coding genes and is thought to be the dominant mechanism of translation initiation used by bacteria. Experimental studies have shown that the SD sequence facilitates start codon recognition and enhances translation initiation by directly interacting with the highly conserved anti-SD sequence on the 30S ribosomal subunit. However, the proportion of SD-led genes within a genome varies across species and the factors governing this variation in translation initiation mechanisms remain largely unknown. Here, we conduct a phylogenetically informed analysis and find that species capable of rapid growth contain a higher proportion of SD-led genes throughout their genomes. We show that SD sequence utilization covaries with a suite of genomic features that are important for efficient translation initiation and elongation. In addition to these endogenous genomic factors, we further show that exogenous environmental factors may influence the evolution of translation initiation mechanisms by finding that thermophilic species contain significantly more SD-led genes than mesophiles. Our results demonstrate that variation in translation initiation mechanisms across bacterial species is predictable and is a consequence of differential life-history strategies related to maximum growth rate and environmental-specific constraints.

Multidomain analyses of a longitudinal human microbiome intestinal cleanout perturbation experiment

Our work focuses on the stability, resilience, and response to perturbation of the bacterial communities in the human gut. Informative flash flood-like disturbances that eliminate most gastrointestinal biomass can be induced using a clinically-relevant iso-osmotic agent. We designed and executed such a disturbance in human volunteers using a dense longitudinal sampling scheme extending before and after induced diarrhea. This experiment has enabled a careful multidomain analysis of a controlled perturbation of the human gut microbiota with a new level of resolution. These new longitudinal multidomain data were analyzed using recently developed statistical methods that demonstrate improvements over current practices. By imposing sparsity constraints we have enhanced the interpretability of the analyses and by employing a new adaptive generalized principal components analysis, incorporated modulated phylogenetic information and enhanced interpretation through scoring of the portions of the tree most influenced by the perturbation. Our analyses leverage the taxa-sample duality in the data to show how the gut microbiota recovers following this perturbation. Through a holistic approach that integrates phylogenetic, metagenomic and abundance information, we elucidate patterns of taxonomic and functional change that characterize the community recovery process across individuals. We provide complete code and illustrations of new sparse statistical methods for high-dimensional, longitudinal multidomain data that provide greater interpretability than existing methods.

Nutrient and Rainfall Additions Shift Phylogenetically Estimated Traits of Soil Microbial Communities

Microbial traits related to ecological responses and functions could provide a common currency facilitating synthesis and prediction; however, such traits are difficult to measure directly for all taxa in environmental samples. Past efforts to estimate trait values based on phylogenetic relationships have not always distinguished between traits with high and low phylogenetic conservatism, limiting reliability, especially in poorly known environments, such as soil. Using updated reference trees and phylogenetic relationships, we estimated two phylogenetically conserved traits hypothesized to be ecologically important from DNA sequences of the 16S rRNA gene from soil bacterial and archaeal communities. We sampled these communities from an environmental change experiment in California grassland applying factorial addition of late-season precipitation and soil nutrients to multiple soil types for 3 years prior to sampling. Estimated traits were rRNA gene copy number, which contributes to how rapidly a microbe can respond to an increase in resources and may be related to its maximum growth rate, and genome size, which suggests the breadth of environmental and substrate conditions in which a microbe can thrive. Nutrient addition increased community-weighted mean estimated rRNA gene copy number and marginally increased estimated genome size, whereas precipitation addition decreased these community means for both estimated traits. The effects of both treatments on both traits were associated with soil properties, such as ammonium, available phosphorus, and pH. Estimated trait responses within several phyla were opposite to the community mean response, indicating that microbial responses, although largely consistent among soil types, were not uniform across the tree of life. Our results show that phylogenetic estimation of microbial traits can provide insight into how microbial ecological strategies interact with environmental changes. The method could easily be applied to any of the thousands of existing 16S rRNA sequence data sets and offers potential to improve our understanding of how microbial communities mediate ecosystem function responses to global changes.

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.

A Numbers Game: Ribosome Densities, Bacterial Growth, and Antibiotic-Mediated Stasis and Death

We postulate that the inhibition of growth and low rates of mortality of bacteria exposed to ribosome-binding antibiotics deemed bacteriostatic can be attributed almost uniquely to these drugs reducing the number of ribosomes contributing to protein synthesis, i.e., the number of effective ribosomes. We tested this hypothesis with Escherichia coli K-12 MG1655 and constructs that had been deleted for 1 to 6 of the 7 rRNA (rrn) operons. In the absence of antibiotics, constructs with fewer rrn operons have lower maximum growth rates and longer lag phases than those with more ribosomal operons. In the presence of the ribosome-binding “bacteriostatic” antibiotics tetracycline, chloramphenicol, and azithromycin, E. coli strains with 1 and 2 rrn operons are killed at a substantially higher rate than those with more rrn operons. This increase in the susceptibility of E. coli with fewer rrn operons to killing by ribosome-targeting bacteriostatic antibiotics is not reflected in their greater sensitivity to killing by the bactericidal antibiotic ciprofloxacin, which does not target ribosomes, but also to killing by gentamicin, which does. Finally, when such strains are exposed to these ribosome-targeting bacteriostatic antibiotics, the time before these bacteria start to grow again when the drugs are removed, referred to as the post-antibiotic effect (PAE), is markedly greater for constructs with fewer rrn operons than for those with more rrn operons. We interpret the results of these other experiments reported here as support for the hypothesis that the reduction in the effective number of ribosomes due to binding to these structures provides a sufficient explanation for the action of bacteriostatic antibiotics that target these structures.

Chemotherapeutic agents, including antibiotics, have been used for more than a century; nevertheless, there are still major gaps in our understanding of how these drugs operate which limit future advances in antibacterial chemotherapy. Although the molecular mechanisms by which antibiotics bind to their target structures are largely known, fundamental questions about how these drugs actually kill and/or inhibit the replication of bacteria remain unanswered and subjects of controversy. We postulate that for the broad class of ribosome-binding bacteriostatic antibiotics, their reducing the number of active (functional) ribosomes per cell provides a sufficient explanation for the abatement of replication and the low rate of decline in densities of viable cells of bacteria exposed to these drugs. Using E. coli K-12 constructs with deletions of from one to six of the seven ribosome-RNA operons and the ribosome-binding bacteriostatic antibiotics tetracycline, chloramphenicol, and azithromycin, we tested this hypothesis. The results of our experiments are consistent with this “numbers game” hypothesis.

Transcriptome Sequencing Reveals Large-Scale Changes in Axenic Aedes aegypti Larvae

Mosquitoes host communities of microbes in their digestive tract that consist primarily of bacteria. We previously reported that Aedes aegypti larvae colonized by a native community of bacteria and gnotobiotic larvae colonized by only Escherichia coli develop very similarly into adults, whereas axenic larvae never molt and die as first instars. In this study, we extended these findings by first comparing the growth and abundance of bacteria in conventional, gnotobiotic, and axenic larvae during the first instar. Results showed that conventional and gnotobiotic larvae exhibited no differences in growth, timing of molting, or number of bacteria in their digestive tract. Axenic larvae in contrast grew minimally and never achieved the critical size associated with molting by conventional and gnotobiotic larvae. In the second part of the study we compared patterns of gene expression in conventional, gnotobiotic and axenic larvae by conducting an RNAseq analysis of gut and nongut tissues (carcass) at 22 h post-hatching. Approximately 12% of Ae. aegypti transcripts were differentially expressed in axenic versus conventional or gnotobiotic larvae. However, this profile consisted primarily of transcripts in seven categories that included the down-regulation of select peptidases in the gut and up-regulation of several genes in the gut and carcass with roles in amino acid transport, hormonal signaling, and metabolism. Overall, our results indicate that axenic larvae exhibit alterations in gene expression consistent with defects in acquisition and assimilation of nutrients required for growth.

rRNA Operon Copy Number Can Explain the Distinct Epidemiology of Hospital-Associated Methicillin-Resistant Staphylococcus aureus

The distinct epidemiology of original hospital-associated methicillin-resistant Staphylococcus aureus (HA-MRSA) and early community-associated MRSA (CA-MRSA) is largely unexplained. S. aureus carries either five or six rRNA operon copies. Evidence is provided for a scenario in which MRSA has adapted to the hospital environment by rRNA operon loss (six to five copies) due to antibiotic pressure. Early CA-MRSA, in contrast, results from wild-type methicillin-susceptible S. aureus (MSSA) that acquired mecA without loss of an rRNA operon. Of the HA-MRSA isolates (n = 77), 67.5% had five rRNA operon copies, compared to 23.2% of the CA-MRSA isolates (n = 69) and 7.7% of MSSA isolates (n = 195) (P < 0.001). In addition, 105 MSSA isolates from cystic fibrosis patients were tested, because these patients are repeatedly treated with antibiotics; 32.4% of these isolates had five rRNA operon copies. For all subsets, a correlation between resistance profile and rRNA copy number was found. Furthermore, we showed that in vitro antibiotic pressure may result in rRNA operon copy loss. We also showed that without antibiotic pressure, S. aureus isolates containing six rRNA copies are more fit than isolates with five copies. We conclude that HA-MRSA and cystic fibrosis isolates most likely have adapted to an environment with high antibiotic pressure by the loss of an rRNA operon copy. This loss has facilitated resistance development, which promoted survival in these niches. However, strain fitness decreased, which explains their lack of success in the community. In contrast, CA-MRSA isolates retained six rRNA operon copies, rendering them fitter and thereby able to survive and spread in the community.

Exploiting rRNA operon copy number to investigate bacterial reproductive strategies

The potential for rapid reproduction is a hallmark of microbial life, but microbes in nature must also survive and compete when growth is constrained by resource availability. Successful reproduction requires different strategies when resources are scarce compared to when they are abundant^1,2, but a systematic framework for predicting these reproductive strategies in bacteria has not been available. Here we show that the number of ribosomal RNA operons (rrn) in bacterial genomes predicts two important components of reproduction – growth rate and growth efficiency – which are favored under contrasting regimes of resource availability^3,4. We find that the maximum reproductive rate of bacteria doubles with a doubling of rrn copy number, while the efficiency of carbon use is inversely related to maximal growth rate and rrn copy number. We also identify a feasible explanation for these patterns: the rate and yield of protein synthesis mirror the overall pattern in maximum growth rate and growth efficiency. Furthermore, comparative analysis of genomes from 1,167 bacterial species reveals that rrn copy number predicts traits associated with resource availability, including chemotaxis and genome streamlining. Genome-wide patterns of orthologous gene content covary with rrn copy number, suggesting convergent evolution in response to resource availability. Our findings indicate that basic cellular processes adapt in contrasting ways to long-term differences in resource availability. They also establish a basis for predicting changes in bacterial community composition in response to resource perturbations using rrn copy number measurements⁵ or inferences^6,7.

Biogeochemistry and Geographical Ecology: Embracing All Twenty-Five Elements Required to Build Organisms

Biogeochemistry is a key but relatively neglected part of the abiotic template that underlies ecology. The template has a geography, one that is increasingly being rearranged in this era of global change. Justus von Liebig's law of the minimum has played a useful role in focusing attention on biogeochemical regulation of populations, but given that ∼25+ elements are required to build organisms and that these organisms use and deplete nutrients in aggregates of communities and ecosystems, we make the case that it is time to move on. We review available models that suggest the many different mechanisms that give rise to multiple elements, or colimitation. We then review recent empirical data that show that rates of decomposition and primary productivity may be limited by multiple elements. In that light, given the tropics' high species diversity and generally more weathered soils, we predict that colimitation at community and ecosystem scales is more prevalent closer to the equator. We conclude with suggestions for how to move forward with experimental studies of colimitation.

The Ecology of Acidobacteria: Moving beyond Genes and Genomes

The phylum Acidobacteria is one of the most widespread and abundant on the planet, yet remarkably our knowledge of the role of these diverse organisms in the functioning of terrestrial ecosystems remains surprisingly rudimentary. This blatant knowledge gap stems to a large degree from the difficulties associated with the cultivation of these bacteria by classical means. Given the phylogenetic breadth of the Acidobacteria, which is similar to the metabolically diverse Proteobacteria, it is clear that detailed and functional descriptions of acidobacterial assemblages are necessary. Fortunately, recent advances are providing a glimpse into the ecology of members of the phylum Acidobacteria. These include novel cultivation and enrichment strategies, genomic characterization and analyses of metagenomic DNA from environmental samples. Here, we couple the data from these complementary approaches for a better understanding of their role in the environment, thereby providing some initial insights into the ecology of this important phylum. All cultured acidobacterial type species are heterotrophic, and members of subdivisions 1, 3, and 4 appear to be more versatile in carbohydrate utilization. Genomic and metagenomic data predict a number of ecologically relevant capabilities for some acidobacteria, including the ability to: use of nitrite as N source, respond to soil macro-, micro nutrients and soil acidity, express multiple active transporters, degrade gellan gum and produce exopolysaccharide (EPS). Although these predicted properties allude to a competitive life style in soil, only very few of these prediction shave been confirmed via physiological studies. The increased availability of genomic and physiological information, coupled to distribution data in field surveys and experiments, should direct future progress in unraveling the ecology of this important but still enigmatic phylum.