Molecular identification of micro- and macroorganisms based on nuclear markers has revolutionized our understanding of their taxonomy, phylogeny and ecology. Today, research on the diversity of eukaryotes in global ecosystems heavily relies on nuclear ribosomal RNA (rRNA) markers. Here, we present the research community-curated reference database EUKARYOME for nuclear ribosomal 18S rRNA, internal transcribed spacer (ITS) and 28S rRNA markers for all eukaryotes, including metazoans (animals), protists, fungi and plants. It is particularly useful for the identification of arbuscular mycorrhizal fungi as it bridges the four commonly used molecular markers-ITS1, ITS2, 18S V4-V5 and 28S D1-D2 subregions. The key benefits of this database over other annotated reference sequence databases are that it is not restricted to certain taxonomic groups and it includes all rRNA markers. EUKARYOME also offers a number of reference long-read sequences that are derived from (meta)genomic and (meta)barcoding-a unique feature that can be used for taxonomic identification and chimera control of third-generation, long-read, high-throughput sequencing data. Taxonomic assignments of rRNA genes in the database are verified based on phylogenetic approaches. The reference datasets are available in multiple formats from the project homepage, http://www.eukaryome.org.

The tumor suppressor p14/19ARF regulates ribosomal RNA (rRNA) synthesis by controlling the nucleolar localization of Transcription Termination Factor 1 (TTF1). However, the role played by TTF1 in regulating the rRNA genes and in potentially controlling growth has remained unclear. We now show that TTF1 expression regulates cell growth by determining the cellular complement of ribosomes. Unexpectedly, it achieves this by acting as a \"roadblock\" to synthesis of the noncoding LncRNA and pRNA that we show are generated from the \"Spacer Promoter\" duplications present upstream of the 47S pre-rRNA promoter on the mouse and human ribosomal RNA genes. Unexpectedly, the endogenous generation of these noncoding RNAs does not induce CpG methylation or gene silencing. Rather, it acts in cis to suppress 47S preinitiation complex formation and hence de novo pre-rRNA synthesis by a mechanism reminiscent of promoter interference or occlusion. Taken together, our data delineate a pathway from p19ARF to cell growth suppression via the regulation of ribosome biogenesis by noncoding RNAs and validate a key cellular growth law in mammalian cells.

Ciliates play a key role in most ecosystems. Their abundance in natural samples is crucial for answering many ecological questions. Traditional methods of quantifying individual species, which rely on microscopy, are often labour-intensive, time-consuming and can be highly biassed. As a result, we investigated the potential of digital polymerase chain reaction (dPCR) for quantifying ciliates. A significant challenge in this process is the high variation in the copy number of the taxonomic marker gene (ribosomal RNA [rRNA]). We first quantified the rRNA gene copy numbers (GCN) of the model ciliate, Paramecium tetraurelia, during different stages of the cell cycle and growth phases. The per-cell rRNA GCN varied between approximately 11,000 and 130,000, averaging around 50,000 copies per cell. Despite these variations in per-cell rRNA GCN, we found a highly significant correlation between GCN and cell numbers. This is likely due to the coexistence of different cellular stages in an uncontrolled (environmental) ciliate population. Thanks to the high sensitivity of dPCR, we were able to detect the target gene in a sample that contained only a single cell. The dPCR approach presented here is a valuable addition to the molecular toolbox in protistan ecology. It may guide future studies in quantifying and monitoring the abundance of targeted (even rare) ciliates in natural samples.

Taxonomic classification using metabarcoding is a commonly used method in microbiological studies of environmental samples and during monitoring of biotechnological processes. However, it is difficult to compare results from different laboratories, due to the variety of bioinformatics tools that have been developed and used for data analysis. This problem is compounded by different choices regarding which variable region of the 16S rRNA gene and which database is used for taxonomic identification. Therefore, this study employed the DADA2 algorithm to optimize the preprocessing of raw data obtained from the sequencing of activated sludge samples, using simultaneous analysis of three frequently used regions of 16S rRNA (V1-V3, V3-V4, V4-V5). Additionally, the study evaluated which variable region and which of the frequently used microbial databases for taxonomic classification (Greengenes2, Silva, RefSeq) more accurately classify OTUs into taxa. Adjusting the values of selected parameters of the DADA2 algorithm, we obtained the highest possible numbers of OTUs for each region. Regarding biodiversity within regions, the V3-V4 region had the highest Simpson and Shannon indexes, and the Chao1 index was similar to that of the V1-V3 region. Beta-biodiversity analysis revealed statistically significant differences between regions. When comparing databases for each of the regions studied, the highest numbers of taxonomic groups were obtained using the SILVA database. These results suggest that standardization of metabarcoding of short amplicons may be possible.

BACKGROUND: Sequencing variable regions of the 16S rRNA gene (≃300 bp) with Illumina technology is commonly used to study the composition of human microbiota. Unfortunately, short reads are unable to differentiate between highly similar species. Considering that species from the same genus can be associated with health or disease it is important to identify them at the lowest possible taxonomic rank. Third-generation sequencing platforms such as PacBio SMRT, increase read lengths allowing to sequence the whole gene with the maximum taxonomic resolution. Despite its potential, full length 16S rRNA gene sequencing is not widely used yet. The aim of the current study was to compare the sequencing output and taxonomic annotation performance of the two approaches (Illumina short read sequencing and PacBio long read sequencing of 16S rRNA gene) in different human microbiome samples. DNA from saliva, oral biofilms (subgingival plaque) and faeces of 9 volunteers was isolated. Regions V3-V4 and V1-V9 were amplified and sequenced by Illumina Miseq and by PacBio Sequel II sequencers, respectively.
RESULTS: With both platforms, a similar percentage of reads was assigned to the genus level (94.79% and 95.06% respectively) but with PacBio a higher proportion of reads were further assigned to the species level (55.23% vs 74.14%). Regarding overall bacterial composition, samples clustered by niche and not by sequencing platform. In addition, all genera with > 0.1% abundance were detected in both platforms for all types of samples. Although some genera such as Streptococcus tended to be observed at higher frequency in PacBio than in Illumina (20.14% vs 14.12% in saliva, 10.63% vs 6.59% in subgingival plaque biofilm samples) none of the differences were statistically significant when correcting for multiple testing.
CONCLUSIONS: The results presented in the current manuscript suggest that samples sequenced using Illumina and PacBio are mostly comparable. Considering that PacBio reads were assigned at the species level with higher accuracy than Illumina, our data support the use of PacBio technology for future microbiome studies, although a higher cost is currently required to obtain an equivalent number of reads per sample.

Research has evolved on aerobic granular sludge (AGS) process, but still there are very few studies on the treatment of excess AGS sludge, with almost none considering its aerobic digestion. Here therefore, the aerobic digestibility of typical AGS sludge was assessed. Granules were produced from acetate-based synthetic wastewater (WW) and were subjected to aerobic digestion for 64 d. The stabilization process was monitored over time through physical-chemical parameters, oxygen uptake rates (OUR) and 16S rRNA gene sequencing. The microbial analyses revealed that the cultivated granules were dominated by slow-growing bacteria, mainly ordinary heterotrophic organisms with potential for polyhydroxyalkanoates (PHA) aerobic storage (PHA-OHOs), polyphosphate and glycogen accumulating organisms (PAOs and GAOs), fermentative anaerobes and nitrifiers (AOB and NOB). Differential abundance analysis of the bacterial data (before versus after digestion) discriminated between the most vulnerable microbiome genera and those most resistant to aerobic digestion. Furthermore, modeling of the stabilization process determined that the endogenous decay rate constant (bH) for the heterotrophs present in the granules was notably low; bH = 0.05 d-1 (average), four times less than for common activated sludge (AS), which is rated at 0.2 d-1. For first time, the research reveals another important feature of AGS sludge, i.e. the slow-decaying character of its bacteria (along with their known slow-growing character). This results in slower stabilization, need of bigger digesters and reconsideration of the specific OUR limits in biosolids regulations (SOUR limit of 1.5 mg/gTSS.h), for waste AGS compared to conventional waste AS. The study suggests that aerobic digestion of waste AGS (fully-granulated) could differ from that of conventional AS. Future work is needed on aerobic digestibility of real AGS sludges from municipal and industrial WWs, compared to synthetic WWs.

Bacteria are primarily responsible for biological water treatment processes in constructed wetland systems. Gravel in constructed wetlands serves as an essential substrate onto which complex bacterial biofilms may successfully grow and evolve. To fully understand the bacterial community in these systems it is crucial to properly isolate biofilms and process DNA from such substrates. This study looked at how best to isolate bacterial biofilms from gravel substrates in terms of bacterial richness. It considered factors including the duration of agitation during extraction, extraction temperature, and enzyme usage. Further, the 16S taxonomy data subsequently produced from Illumina MiSeq reads (using the SILVA 132 ribosomal RNA (rRNA) database on the DADA2 pipeline) were compared with the 16S data produced from Oxford Nanopore Technologies (ONT) MinION reads (using the NCBI 16S database on the EPI2ME pipeline). Finally, performance was tested by comparing the taxonomy data generated from the Illumina MiSeq and ONT MinION reads using the same (SILVA 132) database. We found no significant differences in the effective number of species observed when using different bacterial biofilm detachment techniques. However, enzyme treatment enhanced the total concentration of DNA. In terms of wetland community profiles, relative abundance differences within each sample type were clearer at the genus level. For genus-level taxonomic classification, MinION sequencing with the EPI2ME pipeline (NCBI database) produced bacterial abundance information that was poorly correlated with that from the Illumina MiSeq and DADA2 pipelines (SILVA132 database). When using the same database for each sequencing technology (SILVA132), the correlation between relative abundances at genus-level improved from negligible to moderate. This study provides detailed information of value to researchers working on constructed wetlands regarding efficient biofilm detachment techniques for DNA isolation and 16 s metabarcoding platforms for sequencing and data analysis.

Next-generation sequencing technology has driven the rapid advancement of human microbiome studies by enabling community-level sequence profiling of microbiomes. Although all microbiome sequencing methods depend on recovering the DNA from a sample as a first critical step, lysis methods can be a major determinant of microbiome profile bias. Gentle enzyme-based DNA preparation methods preserve DNA quality but can bias the results by failing to open difficult-to-lyse bacteria. Mechanical methods like bead beating can also bias DNA recovery because the mechanical energy required to break tougher cell walls may shear the DNA of the more easily lysed microbes, and shearing can vary depending on the time and intensity of beating, influencing reproducibility. We introduce a non-mechanical, non-enzymatic, novel rapid microbial DNA extraction procedure suitable for 16S rRNA gene-based microbiome profiling applications that eliminates bead beating. The simultaneous application of alkaline, heat, and detergent (\'Rapid\' protocol) to milligram quantity samples provided consistent representation across the population of difficult and easily lysed bacteria equal to or better than existing protocols, producing sufficient high-quality DNA for full-length 16S rRNA gene PCR. The novel \'Rapid\' method was evaluated using mock bacterial communities containing both difficult and easily lysed bacteria. Human fecal sample testing compared the novel Rapid method with a standard Human Microbiome Project (HMP) protocol for samples from lung cancer patients and controls. DNA recovered from both methods was analyzed using 16S rRNA gene sequencing of the V1V3 and V4 regions on the Illumina platform and the V1V9 region on the PacBio platform. Our findings indicate that the \'Rapid\' protocol consistently yielded higher levels of Firmicutes species, which reflected the profile of the bacterial community structure more accurately, which was confirmed by mock community evaluation. The novel \'Rapid\' DNA lysis protocol reduces population bias common to bead beating and enzymatic lysis methods, presenting opportunities for improved microbial community profiling, combined with the reduction in sample input to 10 milligrams or less, and it enables rapid transfer and simultaneous lysis of 96 samples in a standard plate format. This results in a 20-fold reduction in sample handling time and an overall 2-fold time advantage when compared to widely used commercial methods. We conclude that the novel \'Rapid\' DNA extraction protocol offers a reliable alternative for preparing fecal specimens for 16S rRNA gene amplicon sequencing.

The taxonomic characterization of ancient microbiomes is a key step in the rapidly growing field of paleomicrobiology. While PCR amplification of the 16S ribosomal RNA (rRNA) gene is a widely used technique in modern microbiota studies, this method has systematic biases when applied to ancient microbial DNA. Shotgun metagenomic sequencing has proven to be the most effective method in reconstructing taxonomic profiles of ancient dental calculus samples. Nevertheless, shotgun sequencing approaches come with inherent limitations that could be addressed through hybridization enrichment capture. When employed together, shotgun sequencing and hybridization capture have the potential to enhance the characterization of ancient microbial communities. Here, we develop, test, and apply a hybridization enrichment capture technique to selectively target 16S rRNA gene fragments from the libraries of ancient dental calculus samples generated with shotgun techniques. We simulated data sets generated from hybridization enrichment capture, indicating that taxonomic identification of fragmented and damaged 16S rRNA gene sequences was feasible. Applying this enrichment approach to 15 previously published ancient calculus samples, we observed a 334-fold increase of ancient 16S rRNA gene fragments in the enriched samples when compared to unenriched libraries. Our results suggest that 16S hybridization capture is less prone to the effects of background contamination than 16S rRNA amplification, yielding a higher percentage of on-target recovery. While our enrichment technique detected low abundant and rare taxa within a given sample, these assignments may not achieve the same level of specificity as those achieved by unenriched methods.

Gene sequence has been widely used in molecular ecology. For instance, the ribosomal RNA (rRNA) gene has been widely used as a biological marker to understand microbial communities. The variety of the detected rRNA gene sequences reflects the diversity of the microorganisms existing in the analyzed sample. Their biomass can also be estimated by applying quantitative sequencing with information on rRNA gene copy numbers in genomes; however, information on rRNA gene copy numbers is still limited. Especially, the copy number in microbial eukaryotes is much less understood than that of prokaryotes, possibly because of the large and complex structure of eukaryotic genomes. In this study, we report an alternative approach that is more appropriate than the existing method of quantitative sequencing and demonstrate that the copy number of eukaryotic rRNA can be measured efficiently and comprehensively. By applying this approach widely, information on the eukaryotic rRNA copy number can be determined, and their community structures can be depicted and compared more efficiently.