Robust single-cell Hi-C clustering by convolution- and random-walk-based imputation.

Three-dimensional genome structure plays a pivotal role in gene regulation and cellular function. Single-cell analysis of genome architecture has been achieved using imaging and chromatin conformation capture methods such as Hi-C. To study variation in chromosome structure between different cell types, computational approaches are needed that can utilize sparse and heterogeneous single-cell Hi-C data. However, few methods exist that are able to accurately and efficiently cluster such data into constituent cell types. Here, we describe scHiCluster, a single-cell clustering algorithm for Hi-C contact matrices that is based on imputations using linear convolution and random walk. Using both simulated and real single-cell Hi-C data as benchmarks, scHiCluster significantly improves clustering accuracy when applied to low coverage datasets compared with existing methods. After imputation by scHiCluster, topologically associating domain (TAD)-like structures (TLSs) can be identified within single cells, and their consensus boundaries were enriched at the TAD boundaries observed in bulk cell Hi-C samples. In summary, scHiCluster facilitates visualization and comparison of single-cell 3D genomes.

Expanding the Diversity of Bacterioplankton Isolates and Modeling Isolation Efficacy with Large-Scale Dilution-to-Extinction Cultivation.

Cultivated bacterioplankton representatives from diverse lineages and locations are essential for microbiology, but the large majority of taxa either remain uncultivated or lack isolates from diverse geographic locales. We paired large-scale dilution-to-extinction (DTE) cultivation with microbial community analysis and modeling to expand the phylogenetic and geographic diversity of cultivated bacterioplankton and to evaluate DTE cultivation success. Here, we report results from 17 DTE experiments totaling 7,820 individual incubations over 3 years, yielding 328 repeatably transferable isolates. Comparison of isolates to microbial community data for source waters indicated that we successfully isolated 5% of the observed bacterioplankton community throughout the study; 43% and 26% of our isolates matched operational taxonomic units and amplicon single-nucleotide variants, respectively, within the top 50 most abundant taxa. Isolates included those from previously uncultivated clades such as SAR11 LD12 and Actinobacteria acIV, as well as geographically novel members from other ecologically important groups like SAR11 subclade IIIa, SAR116, and others, providing isolates in eight putatively new genera and seven putatively new species. Using a newly developed DTE cultivation model, we evaluated taxon viability by comparing relative abundance with cultivation success. The model (i) revealed the minimum attempts required for successful isolation of taxa amenable to growth on our media and (ii) identified possible subpopulation viability variation in abundant taxa such as SAR11 that likely impacts cultivation success. By incorporating viability in experimental design, we can now statistically constrain the effort necessary for successful cultivation of specific taxa on a defined medium.IMPORTANCE Even before the coining of the term "great plate count anomaly" in the 1980s, scientists had noted the discrepancy between the number of microorganisms observed under the microscope and the number of colonies that grew on traditional agar media. New cultivation approaches have reduced this disparity, resulting in the isolation of some of the "most wanted" bacterial lineages. Nevertheless, the vast majority of microorganisms remain uncultured, hampering progress toward answering fundamental biological questions about many important microorganisms. Furthermore, few studies have evaluated the underlying factors influencing cultivation success, limiting our ability to improve cultivation efficacy. Our work details the use of dilution-to-extinction (DTE) cultivation to expand the phylogenetic and geographic diversity of available axenic cultures. We also provide a new model of the DTE approach that uses cultivation results and natural abundance information to predict taxon-specific viability and iteratively constrain DTE experimental design to improve cultivation success.

Laboratory evolution reveals a two-dimensional rate-yield tradeoff in microbial metabolism.

Growth rate and yield are fundamental features of microbial growth. However, we lack a mechanistic and quantitative understanding of the rate-yield relationship. Studies pairing computational predictions with experiments have shown the importance of maintenance energy and proteome allocation in explaining rate-yield tradeoffs and overflow metabolism. Recently, adaptive evolution experiments of Escherichia coli reveal a phenotypic diversity beyond what has been explained using simple models of growth rate versus yield. Here, we identify a two-dimensional rate-yield tradeoff in adapted E. coli strains where the dimensions are (A) a tradeoff between growth rate and yield and (B) a tradeoff between substrate (glucose) uptake rate and growth yield. We employ a multi-scale modeling approach, combining a previously reported coarse-grained small-scale proteome allocation model with a fine-grained genome-scale model of metabolism and gene expression (ME-model), to develop a quantitative description of the full rate-yield relationship for E. coli K-12 MG1655. The multi-scale analysis resolves the complexity of ME-model which hindered its practical use in proteome complexity analysis, and provides a mechanistic explanation of the two-dimensional tradeoff. Further, the analysis identifies modifications to the P/O ratio and the flux allocation between glycolysis and pentose phosphate pathway (PPP) as potential mechanisms that enable the tradeoff between glucose uptake rate and growth yield. Thus, the rate-yield tradeoffs that govern microbial adaptation to new environments are more complex than previously reported, and they can be understood in mechanistic detail using a multi-scale modeling approach.

Functional decomposition of metabolism allows a system-level quantification of fluxes and protein allocation towards specific metabolic functions.

Quantifying the contribution of individual molecular components to complex cellular processes is a grand challenge in systems biology. Here we establish a general theoretical framework (Functional Decomposition of Metabolism, FDM) to quantify the contribution of every metabolic reaction to metabolic functions, e.g. the synthesis of biomass building blocks. FDM allowed for a detailed quantification of the energy and biosynthesis budget for growing Escherichia coli cells. Surprisingly, the ATP generated during the biosynthesis of building blocks from glucose almost balances the demand from protein synthesis, the largest energy expenditure known for growing cells. This leaves the bulk of the energy generated by fermentation and respiration unaccounted for, thus challenging the common notion that energy is a key growth-limiting resource. Moreover, FDM together with proteomics enables the quantification of enzymes contributing towards each metabolic function, allowing for a first-principle formulation of a coarse-grained model of global protein allocation based on the structure of the metabolic network.

Ecophysiology and genomics of the brackish water adapted SAR11 subclade IIIa.

The Order Pelagibacterales (SAR11) is the most abundant group of heterotrophic bacterioplankton in global oceans and comprises multiple subclades with unique spatiotemporal distributions. Subclade IIIa is the primary SAR11 group in brackish waters and shares a common ancestor with the dominant freshwater IIIb (LD12) subclade. Despite its dominance in brackish environments, subclade IIIa lacks systematic genomic or ecological studies. Here, we combine closed genomes from new IIIa isolates, new IIIa MAGS from San Francisco Bay (SFB), and 460 highly complete publicly available SAR11 genomes for the most comprehensive pangenomic study of subclade IIIa to date. Subclade IIIa represents a taxonomic family containing three genera (denoted as subgroups IIIa.1, IIIa.2, and IIIa.3) that had distinct ecological distributions related to salinity. The expansion of taxon selection within subclade IIIa also established previously noted metabolic differentiation in subclade IIIa compared to other SAR11 subclades such as glycine/serine prototrophy, mosaic glyoxylate shunt presence, and polyhydroxyalkanoate synthesis potential. Our analysis further shows metabolic flexibility among subgroups within IIIa. Additionally, we find that subclade IIIa.3 bridges the marine and freshwater clades based on its potential for compatible solute transport, iron utilization, and bicarbonate management potential. Pure culture experimentation validated differential salinity ranges in IIIa.1 and IIIa.3 and provided detailed IIIa cell size and volume data. This study is an important step forward for understanding the genomic, ecological, and physiological differentiation of subclade IIIa and the overall evolutionary history of SAR11.

The AEGEAN-169 clade of bacterioplankton is synonymous with SAR11 subclade V (HIMB59) and metabolically distinct.

Bacterioplankton of the SAR11 clade are the most abundant marine microorganisms and consist of numerous subclades spanning order-level divergence (Pelagibacterales). The assignment of the earliest diverging subclade V (a.k.a. HIMB59) to the Pelagibacterales is highly controversial, with multiple recent phylogenetic studies placing them completely separate from SAR11. Other than through phylogenomics, subclade V has not received detailed examination due to limited genomes from this group. Here, we assessed the ecogenomic characteristics of subclade V to better understand the role of this group in comparison to the Pelagibacterales. We used a new isolate genome, recently released single-amplified genomes and metagenome-assembled genomes, and previously established SAR11 genomes to perform a comprehensive comparative genomics analysis. We paired this analysis with the recruitment of metagenomes spanning the open ocean, coastal, and brackish systems. Phylogenomics, average amino acid identity, and 16S rRNA gene phylogeny indicate that SAR11 subclade V is synonymous with the ubiquitous AEGEAN-169 clade and support the contention that this group represents a taxonomic family. AEGEAN-169 shared many bulk genome qualities with SAR11, such as streamlining and low GC content, but genomes were generally larger. AEGEAN-169 had overlapping distributions with SAR11 but was metabolically distinct from SAR11 in its potential to transport and utilize a broader range of sugars as well as in the transport of trace metals and thiamin. Thus, regardless of the ultimate phylogenetic placement of AEGEAN-169, these organisms have distinct metabolic capacities that likely allow them to differentiate their niche from canonical SAR11 taxa. IMPORTANCE One goal of marine microbiologists is to uncover the roles various microorganisms are playing in biogeochemical cycles. Success in this endeavor relies on differentiating groups of microbes and circumscribing their relationships. An early-diverging group (subclade V) of the most abundant bacterioplankton, SAR11, has recently been proposed as a separate lineage that does not share a most recent common ancestor. But beyond phylogenetics, little has been done to evaluate how these organisms compare with SAR11. Our work leverages dozens of new genomes to demonstrate the similarities and differences between subclade V and SAR11. In our analysis, we also establish that subclade V is synonymous with a group of bacteria established from 16S rRNA gene sequences, AEGEAN-169. Subclade V/AEGEAN-169 has clear metabolic distinctions from SAR11 and their shared traits point to remarkable convergent evolution if they do not share a most recent common ancestor.

sparse-growth-curve: a Computational Pipeline for Parsing Cellular Growth Curves with Low Temporal Resolution.

Here, we introduce a Python-based repository, sparse-growth-curve, a software package designed for parsing cellular growth curves with low temporal resolution. The repository uses cell density and time data as the input, automatically separates different growth phases, calculates the exponential growth rates, and produces multiple graphs to aid in interpretation.

Ecophysiology of the Cosmopolitan OM252 Bacterioplankton (Gammaproteobacteria).

Among the thousands of species that comprise marine bacterioplankton communities, most remain functionally obscure. One key cosmopolitan group in this understudied majority is the OM252 clade of Gammaproteobacteria. Although frequently found in sequence data and even previously cultured, the diversity, metabolic potential, physiology, and distribution of this clade has not been thoroughly investigated. Here, we examined these features of OM252 bacterioplankton using a newly isolated strain and genomes from publicly available databases. We demonstrated that this group constitutes a globally distributed novel genus ("Candidatus Halomarinus"), sister to Litoricola, comprising two subclades and multiple distinct species. OM252 organisms have small genomes (median, 2.21 Mbp) and are predicted obligate aerobes capable of alternating between chemoorganoheterotrophic and chemolithotrophic growth using reduced sulfur compounds as electron donors. Subclade I genomes encode genes for the Calvin-Benson-Bassham cycle for carbon fixation. One representative strain of subclade I, LSUCC0096, had extensive halotolerance and a mesophilic temperature range for growth, with a maximum rate of 0.36 doublings/h at 35°C. Cells were curved rod/spirillum-shaped, ∼1.5 by 0.2 µm. Growth yield on thiosulfate as the sole electron donor under autotrophic conditions was roughly one-third that of heterotrophic growth, even though calculations indicated similar Gibbs energies for both catabolisms. These phenotypic data show that some "Ca. Halomarinus" organisms can switch between serving as carbon sources or sinks and indicate the likely anabolic cost of lithoautotrophic growth. Our results thus provide new hypotheses about the roles of these organisms in global biogeochemical cycling of carbon and sulfur. IMPORTANCE Marine microbial communities are teeming with understudied taxa due to the sheer numbers of species in any given sample of seawater. One group, the OM252 clade of Gammaproteobacteria, has been identified in gene surveys from myriad locations, and one isolated organism has even been genome sequenced (HIMB30). However, further study of these organisms has not occurred. Using another isolated representative (strain LSUCC0096) and publicly available genome sequences from metagenomic and single-cell genomic data sets, we examined the diversity within the OM252 clade and the distribution of these taxa in the world's oceans, reconstructed the predicted metabolism of the group, and quantified growth dynamics in LSUCC0096. Our results generate new knowledge about the previously enigmatic OM252 clade and point toward the importance of facultative chemolithoautotrophy for supporting some clades of ostensibly "heterotrophic" taxa.

Metagenome-Assembled Genome Sequence of Rhodobacteraceae Bacterium Strain Clear-D3, Assembled from a Dolichospermum Consortium from Clear Lake, California.

Here, we report the metagenome-assembled genome sequence of a Rhodobacteraceae bacterium strain, Clear-D3, that was reconstructed from a cyanobacterial enrichment from a eutrophic lake. The draft genome sequence shows evidence of an anoxygenic photoautotrophic lifestyle. Other potential capabilities include aerobic heterotrophy, flagellar motility, chemotaxis, and utilization of complex C-P compounds.