In recent years, massively parallel complementary DNA sequencing (RNA sequencing [RNA-Seq]) has emerged as a fast, cost-effective, and robust technology to study entire transcriptomes in various manners. In particular, for non-model organisms and in the absence of an appropriate reference genome, RNA-Seq is used to reconstruct the transcriptome de novo. Although the de novo transcriptome assembly of non-model organisms has been on the rise recently and new tools are frequently developing, there is still a knowledge gap about which assembly software should be used to build a comprehensive de novo assembly. Here, we present a large-scale comparative study in which 10 de novo assembly tools are applied to 9 RNA-Seq data sets spanning different kingdoms of life. Overall, we built >200 single assemblies and evaluated their performance on a combination of 20 biological-based and reference-free metrics. Our study is accompanied by a comprehensive and extensible Electronic Supplement that summarizes all data sets, assembly execution instructions, and evaluation results. Trinity, SPAdes, and Trans-ABySS, followed by Bridger and SOAPdenovo-Trans, generally outperformed the other tools compared. Moreover, we observed species-specific differences in the performance of each assembler. No tool delivered the best results for all data sets. We recommend a careful choice and normalization of evaluation metrics to select the best assembling results as a critical step in the reconstruction of a comprehensive de novo transcriptome assembly.

Sequence analyses of RNA virus genomes remain challenging owing to the exceptional genetic plasticity of these viruses. Because of high mutation and recombination rates, genome replication by viral RNA-dependent RNA polymerases leads to populations of closely related viruses, so-called “quasispecies.” Standard (short-read) sequencing technologies are ill-suited to reconstruct large numbers of full-length haplotypes of (1) RNA virus genomes and (2) subgenome-length (sg) RNAs composed of noncontiguous genome regions. Here, we used a full-length, direct RNA sequencing (DRS) approach based on nanopores to characterize viral RNAs produced in cells infected with a human coronavirus. By using DRS, we were able to map the longest (∼26-kb) contiguous read to the viral reference genome. By combining Illumina and Oxford Nanopore sequencing, we reconstructed a highly accurate consensus sequence of the human coronavirus (HCoV)-229E genome (27.3 kb). Furthermore, by using long reads that did not require an assembly step, we were able to identify, in infected cells, diverse and novel HCoV-229E sg RNAs that remain to be characterized. Also, the DRS approach, which circumvents reverse transcription and amplification of RNA, allowed us to detect methylation sites in viral RNAs. Our work paves the way for haplotype-based analyses of viral quasispecies by showing the feasibility of intra-sample haplotype separation. Even though several technical challenges remain to be addressed to exploit the potential of the nanopore technology fully, our work illustrates that DRS may significantly advance genomic studies of complex virus populations, including predictions on long-range interactions in individual full-length viral RNA haplotypes.

The terrestrial subsurface contains nearly all of Earth’s freshwater reserves and harbours the majority of our planet’s total prokaryotic biomass. Although genetic surveys suggest these organisms rely on in situ carbon fixation, rather than the photosynthetically derived organic carbon transported from surface environments, direct measurements of carbon fixation in the subsurface are absent. Using an ultra-low level¹⁴C-labelling technique, we estimate in situ carbon fixation rates in a carbonate aquifer. We find these rates are similar to those measured in oligotrophic marine surface waters and up to six-fold greater than those observed in the lower euphotic zone. Our empirical carbon fixation rates agree with nitrification rate data. Metagenomic analyses reveal abundant putative chemolithoautotrophic members of an uncharacterized order of Nitrospiria that may be behind the carbon fixation. On the basis of our determined carbon fixation rates, we conservatively extrapolate global primary production in carbonate groundwaters (10% of global reserves) to be 0.11 Pg carbon per year. These rates fall within the range found for oligotrophic marine surface waters, indicating a substantial contribution of in situ primary production to subsurface ecosystem processes. We further suggest that, just as phototrophs are for marine biogeochemical cycling, such subsurface carbon fixation is potentially foundational to subsurface trophic webs.

Background: The highly diverseCand. Patescibacteria are predicted to have minimal biosynthetic and metabolic pathways, which hinders understanding of how their populations differentiate in response to environmental drivers or host organisms. Their mechanisms employed to cope with oxidative stress are largely unknown. Here, we utilized genome-resolved metagenomics to investigate the adaptive genome repertoire of Patescibacteria in oxic and anoxic groundwaters, and to infer putative host ranges.Results: Within six groundwater wells,Cand. Patescibacteria was the most dominant (up to 79%) super-phylum across 32 metagenomes sequenced from DNA retained on 0.2 and 0.1 µm filters after sequential filtration. Of the reconstructed 1275 metagenome-assembled genomes (MAGs), 291 high-quality MAGs were classified asCand. Patescibacteria.Cand. Paceibacteria andCand. Microgenomates were enriched exclusively in the 0.1 µm fractions, whereas candidate division ABY1 andCand. Gracilibacteria were enriched in the 0.2 µm fractions. On average, Patescibacteria enriched in the smaller 0.1 µm filter fractions had 22% smaller genomes, 13.4% lower replication measures, higher proportion of rod-shape determining proteins, and of genomic features suggesting type IV pili mediated cell–cell attachments. Near-surface wells harbored Patescibacteria with higher replication rates than anoxic downstream wells characterized by longer water residence time. Except prevalence of superoxide dismutase genes in Patescibacteria MAGs enriched in oxic groundwaters (83%), no major metabolic or phylogenetic differences were observed. The most abundant Patescibacteria MAG in oxic groundwater encoded a nitrate transporter, nitrite reductase, and F-type ATPase, suggesting an alternative energy conservation mechanism. Patescibacteria consistently co-occurred with one another or with members of phyla Nanoarchaeota, Bacteroidota, Nitrospirota, and Omnitrophota. Among the MAGs enriched in 0.2 µm fractions,, only 8% Patescibacteria showed highly significant one-to-one correlation, mostly with Omnitrophota. Motility and transport related genes in certain Patescibacteria were highly similar to genes from other phyla (Omnitrophota, Proteobacteria and Nanoarchaeota).Conclusion: Other than genes to cope with oxidative stress, we found little genomic evidence for niche adaptation of Patescibacteria to oxic or anoxic groundwaters. Given that we could detect specific host preference only for a few MAGs, we speculate that the majority of Patescibacteria is able to attach multiple hosts just long enough to loot or exchange supplies.

Organism aging is characterized by increased inflammation and decreased stem cell function, yet the relationship between these factors remains incompletely understood. This study shows that aged hematopoietic stem and progenitor cells (HSPCs) exhibit increased ground-stage NF-κB activity, which enhances their responsiveness to undergo differentiation and loss of self-renewal in response to inflammation. The study identifies Rad21/cohesin as a critical mediator of NF-κB signaling, which increases chromatin accessibility in the vicinity of NF-κB target genes in response to inflammation. Rad21 is required for normal differentiation, but limits self-renewal of hematopoietic stem cells (HSCs) during aging and inflammation in an NF-κB–dependent manner. HSCs from aged mice fail to down-regulate Rad21/cohesin and inflammation/differentiation signals in the resolution phase of inflammation. Inhibition of cohesin/NF-κB reverts hypersensitivity of aged HSPCs to inflammation-induced differentiation and myeloid-biased HSCs with disrupted/reduced expression of Rad21/cohesin are increasingly selected during aging. Together, Rad21/cohesin-mediated NF-κB signaling limits HSPC function during aging and selects for cohesin-deficient HSCs with myeloid-skewed differentiation.