Dual RNA-seq in filarial nematodes and Wolbachia endosymbionts using RNase H based ribosomal RNA depletion.

Lymphatic filariasis is caused by parasitic nematodes and is a leading cause of disability worldwide. Many filarial worms contain the bacterium Wolbachia as an obligate endosymbiont. RNA sequencing is a common technique used to study their molecular relationships and to identify potential drug targets against the nematode and bacteria. Ribosomal RNA (rRNA) is the most abundant RNA species, accounting for 80-90% of the RNA in a sample. To reduce sequencing costs, it is necessary to remove ribosomal reads through poly-A enrichment or ribosomal depletion. Bacterial RNA does not contain a poly-A tail, making it difficult to sequence both the nematode and Wolbachia from the same library preparation using standard poly-A selection. Ribosomal depletion can utilize species-specific oligonucleotide probes to remove rRNA through pull-down or degradation methods. While species-specific probes are commercially available for many commonly studied model organisms, there are currently limited depletion options for filarial parasites. Here, we performed total RNA sequencing from Brugia malayi containing the Wolbachia symbiont (wBm) and designed ssDNA depletion probes against their rRNA sequences. We compared the total RNA library to poly-A enriched, Terminator 5'-Phosphate-Dependent Exonuclease treated, NEBNext Human/Bacteria rRNA depleted and our custom nematode probe depleted libraries. The custom nematode depletion library had the lowest percentage of ribosomal reads across all methods, with a 300-fold decrease in rRNA when compared to the total RNA library. The nematode depletion libraries also contained the highest percentage of Wolbachia mRNA reads, resulting in a 16-1,000-fold increase in bacterial reads compared to the other enrichment and depletion methods. Finally, we found that the Brugia malayi depletion probes can remove rRNA from the filarial worm Dirofilaria immitis and the majority of rRNA from the more distantly related free living nematode Caenorhabditis elegans. These custom filarial probes will allow for future dual RNA-seq experiments between nematodes and their bacterial symbionts from a single sequencing library.

Distribution and abundance of tetraether lipid cyclization genes in terrestrial hot springs reflect pH.

Many Archaea produce membrane-spanning lipids that enable life in extreme environments. These isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) may contain up to eight cyclopentyl and one cyclohexyl ring, where higher degrees of cyclization are associated with more acidic, hotter or energy-limited conditions. Recently, the genes encoding GDGT ring synthases, grsAB, were identified in two Sulfolobaceae; however, the distribution and abundance of grs homologs across environments inhabited by these and related organisms remain a mystery. To address this, we examined the distribution of grs homologs in relation to environmental temperature and pH, from thermal springs across Earth, where sequences derive from metagenomes, metatranscriptomes, single-cell and cultivar genomes. The abundance of grs homologs shows a strong negative correlation to pH, but a weak positive correlation to temperature. Archaeal genomes and metagenome-assembled genomes (MAGs) that carry two or more grs copies are more abundant in low pH springs. We also find grs in 12 archaeal classes, with the most representatives in Thermoproteia, followed by MAGs of the uncultured Korarchaeia, Bathyarchaeia and Hadarchaeia, while several Nitrososphaeria encodes >3 copies. Our findings highlight the key role of grs-catalysed lipid cyclization in archaeal diversification across hot and acidic environments.