Identifying local associations in biological time series: algorithms, statistical significance, and applications.

Local associations refer to spatial-temporal correlations that emerge from the biological realm, such as time-dependent gene co-expression or seasonal interactions between microbes. One can reveal the intricate dynamics and inherent interactions of biological systems by examining the biological time series data for these associations. To accomplish this goal, local similarity analysis algorithms and statistical methods that facilitate the local alignment of time series and assess the significance of the resulting alignments have been developed. Although these algorithms were initially devised for gene expression analysis from microarrays, they have been adapted and accelerated for multi-omics next generation sequencing datasets, achieving high scientific impact. In this review, we present an overview of the historical developments and recent advances for local similarity analysis algorithms, their statistical properties, and real applications in analyzing biological time series data. The benchmark data and analysis scripts used in this review are freely available at http://github.com/labxscut/lsareview.

Ammonium-Enhanced Arsenic Mobilization from Aquifer Sediments.

Ammonium-related pathways are important for groundwater arsenic (As) enrichment, especially via microbial Fe(III) reduction coupled with anaerobic ammonium oxidation; however, the key pathways (and microorganisms) underpinning ammonium-induced Fe(III) reduction and their contributions to As mobilization in groundwater are still unknown. To address this gap, aquifer sediments hosting high As groundwater from the western Hetao Basin were incubated with 15N-labeled ammonium and external organic carbon sources (including glucose, lactate, and lactate/acetate). Decreases in ammonium concentrations were positively correlated with increases in the total produced Fe(II) (Fe(II)tot) and released As. The molar ratios of Fe(II)tot to oxidized ammonium ranged from 3.1 to 3.7 for all incubations, and the δ15N values of N2 from the headspace increased in 15N-labeled ammonium-treated series, suggesting N2 as the key end product of ammonium oxidation. The addition of ammonium increased the As release by 16.1% to 49.6%, which was more pronounced when copresented with organic electron donors. Genome-resolved metagenomic analyses (326 good-quality MAGs) suggested that ammonium-induced Fe(III) reduction in this system required syntrophic metabolic interactions between bacterial Fe(III) reduction and archaeal ammonium oxidation. The current results highlight the significance of syntrophic ammonium-stimulated Fe(III) reduction in driving As mobilization, which is underestimated in high As groundwater.

Estimating maximal microbial growth rates from cultures, metagenomes, and single cells via codon usage patterns.

Maximal growth rate is a basic parameter of microbial lifestyle that varies over several orders of magnitude, with doubling times ranging from a matter of minutes to multiple days. Growth rates are typically measured using laboratory culture experiments. Yet, we lack sufficient understanding of the physiology of most microbes to design appropriate culture conditions for them, severely limiting our ability to assess the global diversity of microbial growth rates. Genomic estimators of maximal growth rate provide a practical solution to survey the distribution of microbial growth potential, regardless of cultivation status. We developed an improved maximal growth rate estimator and predicted maximal growth rates from over 200,000 genomes, metagenome-assembled genomes, and single-cell amplified genomes to survey growth potential across the range of prokaryotic diversity; extensions allow estimates from 16S rRNA sequences alone as well as weighted community estimates from metagenomes. We compared the growth rates of cultivated and uncultivated organisms to illustrate how culture collections are strongly biased toward organisms capable of rapid growth. Finally, we found that organisms naturally group into two growth classes and observed a bias in growth predictions for extremely slow-growing organisms. These observations ultimately led us to suggest evolutionary definitions of oligotrophy and copiotrophy based on the selective regime an organism occupies. We found that these growth classes are associated with distinct selective regimes and genomic functional potentials.

Phage-bacterial contig association prediction with a convolutional neural network.

MOTIVATION: Phage-host associations play important roles in microbial communities. But in natural communities, as opposed to culture-based lab studies where phages are discovered and characterized metagenomically, their hosts are generally not known. Several programs have been developed for predicting which phage infects which host based on various sequence similarity measures or machine learning approaches. These are often based on whole viral and host genomes, but in metagenomics-based studies, we rarely have whole genomes but rather must rely on contigs that are sometimes as short as hundreds of bp long. Therefore, we need programs that predict hosts of phage contigs on the basis of these short contigs. Although most existing programs can be applied to metagenomic datasets for these predictions, their accuracies are generally low. Here, we develop ContigNet, a convolutional neural network-based model capable of predicting phage-host matches based on relatively short contigs, and compare it to previously published VirHostMatcher (VHM) and WIsH. RESULTS: On the validation set, ContigNet achieves 72-85% area under the receiver operating characteristic curve (AUROC) scores, compared to the maximum of 68% by VHM or WIsH for contigs of lengths between 200 bps to 50 kbps. We also apply the model to the Metagenomic Gut Virus (MGV) catalogue, a dataset containing a wide range of draft genomes from metagenomic samples and achieve 60-70% AUROC scores compared to that of VHM and WIsH of 52%. Surprisingly, ContigNet can also be used to predict plasmid-host contig associations with high accuracy, indicating a similar genetic exchange between mobile genetic elements and their hosts. AVAILABILITY AND IMPLEMENTATION: The source code of ContigNet and related datasets can be downloaded from https://github.com/tianqitang1/ContigNet.

The power of cooperation: Experimental and computational approaches in the functional characterization of bacterial sRNAs.

Trans-acting small regulatory RNAs (sRNAs) are key players in the regulation of gene expression in bacteria. There are hundreds of different sRNAs in a typical bacterium, which in contrast to eukaryotic microRNAs are more heterogeneous in length, sequence composition, and secondary structure. The vast majority of sRNAs function post-transcriptionally by binding to other RNAs (mRNAs, sRNAs) through rather short regions of imperfect sequence complementarity. Besides, every single sRNA may interact with dozens of different target RNAs and impact gene expression either negatively or positively. These facts contributed to the view that the entirety of the regulatory targets of a given sRNA, its targetome, is challenging to identify. However, recent developments show that a more comprehensive sRNAs targetome can be achieved through the combination of experimental and computational approaches. Here, we give a short introduction into these methods followed by a description of two sRNAs, RyhB, and RsaA, to illustrate the particular strengths and weaknesses of these approaches in more details. RyhB is an sRNA involved in iron homeostasis in Enterobacteriaceae, while RsaA is a modulator of virulence in Staphylococcus aureus. Using such a combined strategy, a better appreciation of the sRNA-dependent regulatory networks is now attainable.

Genomic and transcriptomic insights into the habitat adaptation of the diazotrophic paddy-field cyanobacterium Nostoc sphaeroides.

Nitrogen-fixing cyanobacteria are common in paddy fields, one of the most productive wetland ecosystems. Here, we present the complete genome of Nostoc sphaeroides, a paddy-field diazotroph used for food and medicine for more than 1700 years and deciphered the transcriptional regulation during the developmental transition from hormogonia to vegetative filaments with heterocysts. The genome of N. sphaeroides consists of one circular chromosome (6.48 Mb), one of the largest ever reported megaplasmids (2.34 Mb), and seven plasmids. Multiple gene families involved in the adaption to high solar radiation and water fluctuation conditions were found expanded, while genes involved in anoxic adaptation and phosphonate utilization are located on the megaplasmid, suggesting its indispensable role in environmental adaptation. Distinct gene expression patterns were observed during the light-intensity-regulated transition from hormogonia to vegetative filaments, specifically, genes encoding proteins involved in photosynthetic light reaction, carbon fixation, nitrogen metabolism and heterocyst differentiation were significantly upregulated, whereas genes related to cell motility were down-regulated. Our results provide genomic and transcriptomic insights into the adaptation of a filamentous nitrogen-fixing cyanobacterium to the highly dynamic paddy-field habitat, suggesting N. sphaeroides as an excellent system to understand the transition from aquatic to terrestrial habitats and to support sustainable rice production.

Reduced vacuolar ß-1,3-glucan synthesis affects carbohydrate metabolism as well as plastid homeostasis and structure in Phaeodactylum tricornutum.

The ß-1,3-glucan chrysolaminarin is the main storage polysaccharide of diatoms. In contrast to plants and green algae, diatoms and most other algal groups do not accumulate storage polysaccharides in their plastids. The diatom Phaeodactylum tricornutum possesses only a single gene encoding a putative ß-1,3-glucan synthase (PtBGS). Here, we characterize this enzyme by expressing GFP fusion proteins in P. tricornutum and by creating and investigating corresponding gene silencing mutants. We demonstrate that PtBGS is a vacuolar protein located in the tonoplast. Metabolite analyses of two mutant strains with reduced amounts of PtBGS reveal a reduction in their chrysolaminarin content and an increase of soluble sugars and lipids. This indicates that carbohydrates are shunted into alternative pathways when chrysolaminarin production is impaired. The mutant strains show reduced growth and lower photosynthetic capacities, while possessing higher photoprotective abilities than WT cells. Interestingly, a strong reduction in PtBGS expression also results in aberrations of the usually very regular thylakoid membrane patterns, including increased thylakoid thickness, reduced numbers of thylakoids per plastid, and increased numbers of lamellae per thylakoid stack. Our data demonstrate the complex intertwinement of carbohydrate storage in the vacuoles with carbohydrate metabolism, photosynthetic homeostasis, and plastid morphology.

Genomic and transcriptomic insights into the survival of the subaerial cyanobacterium Nostoc flagelliforme in arid and exposed habitats.

The cyanobacterium Nostoc flagelliforme is an extremophile that thrives under extraordinary desiccation and ultraviolet (UV) radiation conditions. To investigate its survival strategies, we performed whole-genome sequencing of N. flagelliforme CCNUN1 and transcriptional profiling of its field populations upon rehydration in BG11 medium. The genome of N. flagelliforme is 10.23 Mb in size and contains 10 825 predicted protein-encoding genes, making it one of the largest complete genomes of cyanobacteria reported to date. Comparative genomics analysis among 20 cyanobacterial strains revealed that genes related to DNA replication, recombination and repair had disproportionately high contributions to the genome expansion. The ability of N. flagelliforme to thrive under extreme abiotic stresses is supported by the acquisition of genes involved in the protection of photosynthetic apparatus, the formation of monounsaturated fatty acids, responses to UV radiation, and a peculiar role of ornithine metabolism. Transcriptome analysis revealed a distinct acclimation strategy to rehydration, including the strong constitutive expression of genes encoding photosystem I assembly factors and the involvement of post-transcriptional control mechanisms of photosynthetic resuscitation. Our results provide insights into the adaptive mechanisms of subaerial cyanobacteria in their harsh habitats and have important implications to understand the evolutionary transition of cyanobacteria from aquatic environments to terrestrial ecosystems.

CRISPR-Cas systems in multicellular cyanobacteria.

Novel CRISPR-Cas systems possess substantial potential for genome editing and manipulation of gene expression. The types and numbers of CRISPR-Cas systems vary substantially between different organisms. Some filamentous cyanobacteria harbor > 40 different putative CRISPR repeat-spacer cassettes, while the number of cas gene instances is much lower. Here we addressed the types and diversity of CRISPR-Cas systems and of CRISPR-like repeat-spacer arrays in 171 publicly available genomes of multicellular cyanobacteria. The number of 1328 repeat-spacer arrays exceeded the total of 391 encoded Cas1 proteins suggesting a tendency for fragmentation or the involvement of alternative adaptation factors. The model cyanobacterium Anabaena sp. PCC 7120 contains only three cas1 genes but hosts three Class 1, possibly one Class 2 and five orphan repeat-spacer arrays, all of which exhibit crRNA-typical expression patterns suggesting active transcription, maturation and incorporation into CRISPR complexes. The CRISPR-Cas system within the element interrupting the Anabaena sp. PCC 7120 fdxN gene, as well as analogous arrangements in other strains, occupy the genetic elements that become excised during the differentiation-related programmed site-specific recombination. This fact indicates the propensity of these elements for the integration of CRISPR-cas systems and points to a previously not recognized connection. The gene all3613 resembling a possible Class 2 effector protein is linked to a short repeat-spacer array and a single tRNA gene, similar to its homologs in other cyanobacteria. The diversity and presence of numerous CRISPR-Cas systems in DNA elements that are programmed for homologous recombination make filamentous cyanobacteria a prolific resource for their study. Abbreviations: Cas: CRISPR associated sequences; CRISPR: Clustered Regularly Interspaced Short Palindromic Repeats; C2c: Class 2 candidate; SDR: small dispersed repeat; TSS: transcriptional start site; UTR: untranslated region.

Depth-dependent microbial metagenomes sampled in the northeastern Indian Ocean.

The northeastern Indian Ocean exhibits distinct hydrographic characteristics influenced by various local and remote forces. Variations in these driving factors may alter the physiochemical properties of seawater, such as dissolved oxygen levels, and affect the diversity and function of microbial communities. How the microbial communities change across water depths spanning a dissolved oxygen gradient has not been well understood. Here we employed both 16S rDNA amplicon and metagenomic sequencing approaches to study the microbial communities collected from different water depths along the E87 transect in the northeastern Indian Ocean. Samples were collected from the surface, Deep Chlorophyll Maximum (DCM), Oxygen Minimum Zone (OMZ), and bathypelagic layers. Proteobacteria were prevalent throughout the water columns, while Thermoproteota were found to be abundant in the aphotic layers. A total of 675 non-redundant metagenome-assembled genomes (MAGs) were constructed, spanning 21 bacterial and 5 archaeal phyla. The community structure and genomic information provided by this dataset offer valuable resources for the analysis of microbial biogeography and metabolism in the northeastern Indian Ocean.

DeepMicroClass sorts metagenomic contigs into prokaryotes, eukaryotes and viruses.

Sequence classification facilitates a fundamental understanding of the structure of microbial communities. Binary metagenomic sequence classifiers are insufficient because environmental metagenomes are typically derived from multiple sequence sources. Here we introduce a deep-learning based sequence classifier, DeepMicroClass, that classifies metagenomic contigs into five sequence classes, i.e. viruses infecting prokaryotic or eukaryotic hosts, eukaryotic or prokaryotic chromosomes, and prokaryotic plasmids. DeepMicroClass achieved high performance for all sequence classes at various tested sequence lengths ranging from 500 bp to 100 kbps. By benchmarking on a synthetic dataset with variable sequence class composition, we showed that DeepMicroClass obtained better performance for eukaryotic, plasmid and viral contig classification than other state-of-the-art predictors. DeepMicroClass achieved comparable performance on viral sequence classification with geNomad and VirSorter2 when benchmarked on the CAMI II marine dataset. Using a coastal daily time-series metagenomic dataset as a case study, we showed that microbial eukaryotes and prokaryotic viruses are integral to microbial communities. By analyzing monthly metagenomes collected at HOT and BATS, we found relatively higher viral read proportions in the subsurface layer in late summer, consistent with the seasonal viral infection patterns prevalent in these areas. We expect DeepMicroClass will promote metagenomic studies of under-appreciated sequence types.

Phagocytosis in Marine Coccolithophore Gephyrocapsa huxleyi: Comparison between Calcified and Non-Calcified Strains.

Coccolithophores play a significant role in marine calcium carbonate production and carbon cycles, attributing to their unique feature of producing calcareous plates, coccoliths. Coccolithophores also possess a haplo-diplontic life cycle, presenting distinct morphology types and calcification states. However, differences in nutrient acquisition strategies and mixotrophic behaviors of the two life phases remain unclear. In this study, we conducted a series of phagocytosis experiments of calcified diploid and non-calcified haploid strains of coccolithophore Gephyrocapsa huxleyi under light and dark conditions. The phagocytosis capability of each strain was examined based on characteristic fluorescent signals from ingested beads using flow cytometry and fluorescence microscopy. The results show a significantly higher phagocytosis percentage on fluorescent beads in the bacterial prey surrogates of the non-calcified haploid Gephyrocapsa huxleyi strain, than the calcified diploid strain with or without light. In addition, the non-calcified diploid cells seemingly to presented a much higher phagocytosis percentage in darkness than under light. The differential phagocytosis capacities between the calcified diploid and non-calcified haploid Gephyrocapsa huxleyi strains indicate potential distinct nutritional strategies at different coccolithophore life and calcifying stages, which may further shed light on the potential strategies that coccolithophore possesses in unfavorable environments such as twilight zones and the expanding coccolithophore niches in the natural marine environment under the climate change scenario.

Recovery of 1887 metagenome-assembled genomes from the South China Sea.

The South China Sea (SCS) is a marginal sea characterized by strong land-sea biogeochemical interactions. SCS has a distinctive landscape with a multitude of seamounts in its basin. Seamounts create "seamount effects" that influence the diversity and distribution of planktonic microorganisms in the surrounding oligotrophic waters. Although the vertical distribution and community structure of marine microorganisms have been explored in certain regions of the global ocean, there is a lack of comprehensive microbial genomic surveys for uncultured microorganisms in SCS, particularly in the seamount regions. Here, we employed a metagenomic approach to study the uncultured microbial communities sampled from the Xianbei seamount region to the North Coast waters of SCS. A total of 1887 non-redundant prokaryotic metagenome-assembled genomes (MAGs) were reconstructed, of which, 153 MAGs were classified as high-quality MAGs based on the MIMAG standards. The community structure and genomic information provided by this dataset could be used to analyze microbial distribution and metabolism in the SCS.

[The principle of dRNA-seq and its applications in prokaryotic tran-scriptome analyses].

The genomic and transcriptomic studies have been greatly accelerated by the next-generation sequencing technology. Currently, there are two main methods in transcriptomic studies, namely, Tiling microarray and RNA-seq. Comparatively, because of its overwhelming superiority in transcript coverage and resolution, as well as decrease of costs, RNA-seq has been employed in transcriptome analyses by an increasing number of research institutes. According to the mRNA enrichment or rRNA depletion methods, RNA-seq can be performed in six ways. Among them, dRNA-seq can descriminate the original transcripts from the processed RNAs based on a differential exonuclease treatment. This kind of method has been broadly applied in studies of prokaryotic transcriptional start sites, small regulatory RNAs, promoters and operons, etc. Here, we reviewed the detailed principle and technical process of dRNA-seq and its applications in prokaryotic transcriptome analyses. Besides, we also summarized the advantages and limitations of dRNA-seq at present stage, and then gave our perspective of its future development and potential applications, expecting to present useful references for civil researchers in related fields.

Complete Genome Sequence of a Cognatishimia activa Strain Assembled from the Phycosphere of the Marine Diatom Skeletonema tropicum.

Cognatishimia activa, previously known as Thalassobius activus, has been frequently isolated from marine environments. Here, we present the complete genome sequence of C. activa strain SOCE 004, assembled from the phycosphere of a long-term laboratory-maintained culture of the diatom Skeletonema tropicum. The complete genome is 3,211,994 bp long, with an average G+C content of 53.69%. The genome contains 3,195 genes, including 3,133 protein-coding genes, 50 tRNAs, and 3 copies each of 5S, 16S, and 23S rRNA genes.

Complete Genome Sequence of Stutzerimonas stutzeri Strain SOCE 002, a Marine Bacterium Isolated from the Surface Seawater of Dapeng Bay.

We report the complete genome sequence of Stutzerimonas stutzeri strain SOCE 002, obtained from Illumina and Oxford Nanopore sequencing. The genome is 4.68 Mb long, with a GC content of 63.5%, and contains 4,334 protein-coding genes, 60 tRNAs, and 12 rRNAs. We expect that this complete genome sequence will provide a reference for both genomic and metabolic analyses of S. stutzeri.

Diversity and Distribution of Harmful Algal Bloom Species from Seamount to Coastal Waters in the South China Sea.

Mount Xianbei is one of the largest shallow seamounts located in the middle of the South China Sea (SCS), which might play a role in shaping the biodiversity of surrounding continental coastal waters, particularly the diversity of phytoplankton species causing frequent harmful algal blooms (HABs) in northern SCS. However, the diversity, composition, and distribution of phytoplankton species in the seamount regions of Xianbei remain largely unexplored. In this study, samples around and outside the seamount regions were collected during a late summer cruise of 2021 to test whether seamounts play a role in HAB species propagation. In total, we identified 19 HAB species across all samples using the ASV-based DNA metabarcoding approach, 6 of which had not been reported previously in the SCS, suggesting a diverse HAB species in the SCS. Specifically, 16 HAB species were found in the seamount region of Xianbei, and 5 of them were also found in the coastal waters, indicating a close connection between seamount and coastal waters. This study was the first attempt to explore HAB species' spatial diversity and vertical distribution in the seamount region of Xianbei at single-nucleotide resolution, which provides a novel explanation for the coastal HAB occurrence in the northern SCS. IMPORTANCE There are a number of seamounts under the water of the South China Sea (SCS). The seamounts might play a role in shaping the biodiversity of surrounding continental coastal waters. However, there is no direct evidence revealing the relationship of the biodiversity of phytoplankton between seamounts and coastal waters in the SCS, especially those species having the potential to form harmful algal blooms (HABs). Some HAB species might proliferate in certain geographic locations, while others may be broadly distributed across oceanic provinces. In this study, we provided a detailed analysis of phytoplankton composition and molecular detection of HAB species from seamount to coastal waters in the SCS, which suggested a strong interaction in the HAB species between the two areas. This finding provides new insights into the diversity and distribution of HABs in seamounts and their role in shaping the composition and the occurrence of HABs in coastal water.

Virus impacted community adaptation in oligotrophic groundwater environment revealed by Hi-C coupled metagenomic and viromic study.

Viruses play a crucial role in microbial mortality, diversity and biogeochemical cycles. Groundwater is the largest global freshwater and one of the most oligotrophic aquatic systems on Earth, but how microbial and viral communities are shaped in this special habitat is largely unexplored. In this study, we collected groundwater samples from 23 to 60 m aquifers at Yinchuan Plain, China. In total, 1920 non-reductant viral contigs were retrieved from metagenomes and viromes constructed by Illumina and Nanopore hybrid sequencing. Only 3% of them could be clustered with known viruses, most of which were Caudoviricetes. Coupling 1.2 Tb Hi-C sequencing with CRISPR matching and homology search, we connected 469 viruses with their hosts while some viral clusters presented a broad-host-range trait. Meanwhile, a large proportion of biosynthesis related auxiliary metabolism genes were identified. Those characteristics might benefit viruses for a better survival in this special oligotrophic environment. Additionally, the groundwater virome showed genomic features distinct from those of the open ocean and wastewater treatment facilities in GC distribution and unannotated gene compositions. This paper expands the current knowledge of the global viromic records and serves as a foundation for a more thorough understanding of viruses in groundwater.

Identification of the HetR recognition sequence upstream of hetZ in Anabaena sp. strain PCC 7120.

HetR is the master regulator of heterocyst differentiation in Anabaena sp. strain PCC 7120 and has been found to specifically bind to an inverted-repeat-containing region upstream of hetP, a heterocyst differentiation gene. However, no such inverted-repeat sequence can be found in promoters of other genes in the genome. hetZ is a gene involved in early heterocyst differentiation. As shown with the gfp reporter gene, transcription from P(hetZ) was correlated to the expression level of hetR and inhibition by RGSGR, the pentapeptide derived from the C terminus of PatS. As detected by electrophoretic mobility shift assay, a recombinant HetR showed specific binding to the region upstream of hetZ, and the binding was inhibited by RGSGR. Tests of a series of the upstream fragments delimited the HetR-binding site to a 40-bp region that shows similarity to that upstream of hetP. The introduction of substitutions of bases conserved in the two HetR-binding sites showed that at least 12 bases are required for recognition by HetR. Deletion of a 51-bp region containing the HetR-binding site completely eliminated the transcription activity of P(hetZ). Based on the HetR recognition sequence of hetZ, those upstream of hetR and patA are proposed.