Bacterial mRNA sequencing is inefficient due to the abundance of ribosomal RNA that is challenging to deplete. While commercial kits target rRNA from common bacterial species, they are frequently inefficient when applied to divergent species, including those from environmental isolates. Similarly, other methods typically employ large probe sets that tile the entire length of rRNAs; however, such approaches are infeasible when applied to many species. Therefore, we present EMBR-seq+, which requires fewer than ten oligonucleotides per rRNA by combining rRNA blocking primers with RNase H-mediated depletion to achieve rRNA removal efficiencies of up to 99% in diverse bacterial species. Further, in more complex microbial co-cultures between F. succinogenes strain UWB7 and anerobic fungi, EMBR-seq+ depleted both bacterial and fungal rRNA, with a 4-fold improvement in bacterial rRNA depletion compared to a commercial kit, thereby demonstrating that the method can be applied to non-model microbial mixtures. Notably, for microbes with unknown rRNA sequences, EMBR-seq+ enables rapid iterations in probe design without requiring to start experiments from total RNA. Finally, efficient depletion of rRNA enabled systematic quantification of the reprogramming of the bacterial transcriptome when cultured in the presence of the anerobic fungi Anaeromyces robustus or Caecomyces churrovis. We observed that F. succinogenes strain UWB7 downregulated several lignocellulose-degrading carbohydrate-active enzymes in the presence of anerobic gut fungi, suggesting close interactions between two cellulolytic species that specialize in different aspects of biomass breakdown. Thus, EMBR-seq+ enables efficient, cost-effective and rapid quantification of the transcriptome to gain insights into non-model microbial systems.

Anaerobic fungi produce biomass-degrading enzymes and natural products that are important to harness for several biotechnology applications. Although progress has been made in the development of methods for extracting nucleic acids for genomic and transcriptomic sequencing of these fungi, most studies are limited in that they do not sample multiple fungal growth phases in batch culture. In this study, we establish a method to harvest RNA from fungal monocultures and fungal–methanogen co-cultures, and also determine an optimal time frame for high-quality RNA extraction from anaerobic fungi. Based on RNA quality and quantity targets, the optimal time frame in which to harvest anaerobic fungal monocultures and fungal-methanogen co-cultures for RNA extraction was 2-5 days of growth post-inoculation. When grown on cellulose, the fungal strain Anaeromyces robustus cocultivated with the methanogen Methanobacterium bryantii upregulated genes encoding fungal carbohydrate-active enzymes and other cellulosome components relative to fungal monocultures during this time frame, but expression patterns changed at 24-hour intervals throughout the fungal growth phase. These results demonstrate the importance of establishing methods to extract high-quality RNA from anaerobic fungi at multiple time points during batch cultivation.