skip to main content


Title: Comparative Analysis of rRNA Removal Methods for RNA-Seq Differential Expression in Halophilic Archaea
Despite intense recent research interest in archaea, the scientific community has experienced a bottleneck in the study of genome-scale gene expression experiments by RNA-seq due to the lack of commercial and specifically designed rRNA depletion kits. The high rRNA:mRNA ratio (80–90%: ~10%) in prokaryotes hampers global transcriptomic analysis. Insufficient ribodepletion results in low sequence coverage of mRNA, and therefore, requires a substantially higher number of replicate samples and/or sequencing reads to achieve statistically reliable conclusions regarding the significance of differential gene expression between case and control samples. Here, we show that after the discontinuation of the previous version of RiboZero (Illumina, San Diego, CA, USA) that was useful in partially or completely depleting rRNA from archaea, archaeal transcriptomics studies have experienced a slowdown. To overcome this limitation, here, we analyze the efficiency for four different hybridization-based kits from three different commercial suppliers, each with two sets of sequence-specific probes to remove rRNA from four different species of halophilic archaea. We conclude that the key for transcriptomic success with the currently available tools is the probe-specificity for the rRNA sequence hybridization. With this paper, we provide insights into the archaeal community for selecting certain reagents and strategies over others depending on the archaeal species of interest. These methods yield improved RNA-seq sensitivity and enhanced detection of low abundance transcripts.  more » « less
Award ID(s):
1651117 1936024
NSF-PAR ID:
10329346
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Biomolecules
Volume:
12
Issue:
5
ISSN:
2218-273X
Page Range / eLocation ID:
682
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Background

    RNA sequencing is a powerful approach to quantify the genome-wide distribution of mRNA molecules in a population to gain deeper understanding of cellular functions and phenotypes. However, unlike eukaryotic cells, mRNA sequencing of bacterial samples is more challenging due to the absence of a poly-A tail that typically enables efficient capture and enrichment of mRNA from the abundant rRNA molecules in a cell. Moreover, bacterial cells frequently contain 100-fold lower quantities of RNA compared to mammalian cells, which further complicates mRNA sequencing from non-cultivable and non-model bacterial species. To overcome these limitations, we report EMBR-seq (Enrichment of mRNA by Blocked rRNA), a method that efficiently depletes 5S, 16S and 23S rRNA using blocking primers to prevent their amplification.

    Results

    EMBR-seq results in 90% of the sequenced RNA molecules from anE. coliculture deriving from mRNA. We demonstrate that this increased efficiency provides a deeper view of the transcriptome without introducing technical amplification-induced biases. Moreover, compared to recent methods that employ a large array of oligonucleotides to deplete rRNA, EMBR-seq uses a single or a few oligonucleotides per rRNA, thereby making this new technology significantly more cost-effective, especially when applied to varied bacterial species. Finally, compared to existing commercial kits for bacterial rRNA depletion, we show that EMBR-seq can be used to successfully quantify the transcriptome from more than 500-fold lower starting total RNA.

    Conclusions

    EMBR-seq provides an efficient and cost-effective approach to quantify global gene expression profiles from low input bacterial samples.

     
    more » « less
  2. Komeili, Arash (Ed.)
    ABSTRACT Histone proteins are found across diverse lineages of Archaea , many of which package DNA and form chromatin. However, previous research has led to the hypothesis that the histone-like proteins of high-salt-adapted archaea, or halophiles, function differently. The sole histone protein encoded by the model halophilic species Halobacterium salinarum , HpyA, is nonessential and expressed at levels too low to enable genome-wide DNA packaging. Instead, HpyA mediates the transcriptional response to salt stress. Here we compare the features of genome-wide binding of HpyA to those of HstA, the sole histone of another model halophile, Haloferax volcanii . hstA , like hpyA , is a nonessential gene. To better understand HpyA and HstA functions, protein-DNA binding data (chromatin immunoprecipitation sequencing [ChIP-seq]) of these halophilic histones are compared to publicly available ChIP-seq data from DNA binding proteins across all domains of life, including transcription factors (TFs), nucleoid-associated proteins (NAPs), and histones. These analyses demonstrate that HpyA and HstA bind the genome infrequently in discrete regions, which is similar to TFs but unlike NAPs, which bind a much larger genomic fraction. However, unlike TFs that typically bind in intergenic regions, HpyA and HstA binding sites are located in both coding and intergenic regions. The genome-wide dinucleotide periodicity known to facilitate histone binding was undetectable in the genomes of both species. Instead, TF-like and histone-like binding sequence preferences were detected for HstA and HpyA, respectively. Taken together, these data suggest that halophilic archaeal histones are unlikely to facilitate genome-wide chromatin formation and that their function defies categorization as a TF, NAP, or histone. IMPORTANCE Most cells in eukaryotic species—from yeast to humans—possess histone proteins that pack and unpack DNA in response to environmental cues. These essential proteins regulate genes necessary for important cellular processes, including development and stress protection. Although the histone fold domain originated in the domain of life Archaea , the function of archaeal histone-like proteins is not well understood relative to those of eukaryotes. We recently discovered that, unlike histones of eukaryotes, histones in hypersaline-adapted archaeal species do not package DNA and can act as transcription factors (TFs) to regulate stress response gene expression. However, the function of histones across species of hypersaline-adapted archaea still remains unclear. Here, we compare hypersaline histone function to a variety of DNA binding proteins across the tree of life, revealing histone-like behavior in some respects and specific transcriptional regulatory function in others. 
    more » « less
  3. null (Ed.)
    Corals from the northern Red Sea and Gulf of Aqaba exhibit extreme thermal tolerance. To examine the underlying gene expression dynamics, we exposed Stylophora pistillata from the Gulf of Aqaba to short-term (hours) and long-term (weeks) heat stress with peak seawater temperatures ranging from their maximum monthly mean of 27 °C (baseline) to 29.5 °C, 32 °C, and 34.5 °C. Corals were sampled at the end of the heat stress as well as after a recovery period at baseline temperature. Changes in coral host and symbiotic algal gene expression were determined via RNA-sequencing (RNA-Seq). Shifts in coral microbiome composition were detected by complementary DNA (cDNA)-based 16S ribosomal RNA (rRNA) gene sequencing. In all experiments up to 32 °C, RNA-Seq revealed fast and pervasive changes in gene expression, primarily in the coral host, followed by a return to baseline gene expression for the majority of coral (>94%) and algal (>71%) genes during recovery. At 34.5 °C, large differences in gene expression were observed with minimal recovery, high coral mortality, and a microbiome dominated by opportunistic bacteria (including Vibrio species), indicating that a lethal temperature threshold had been crossed. Our results show that the S. pistillata holobiont can mount a rapid and pervasive gene expression response contingent on the amplitude and duration of the thermal stress. We propose that the transcriptomic resilience and transcriptomic acclimation observed are key to the extraordinary thermal tolerance of this holobiont and, by inference, of other northern Red Sea coral holobionts, up to seawater temperatures of at least 32 °C, that is, 5 °C above their current maximum monthly mean. 
    more » « less
  4. Abstract Background

    Sexually dimorphic mating behaviors differ between sexes and involve gonadal hormones and possibly sexually dimorphic gene expression in the brain. However, the associations among the brain, gonad, and sexual behavior in teleosts are still unclear. Here, we utilized germ cells-freetdrd12knockout (KO) zebrafish, and steroid synthesis enzymecyp17a1-deficient zebrafish to investigate the differences and interplays in the brain–gonad–behavior axis, and the molecular control of brain dimorphism and male mating behaviors.

    Methods

    Tdrd12+/−;cyp17a1+/−double heterozygous parents were crossed to obtaintdrd12−/−;cyp17a1+/+(tdrd12 KO),tdrd12+/+;cyp17a1−/−(cyp17a1 KO), andtdrd12−/−;cyp17a1−/−(double KO) homozygous progenies. Comparative analysis of mating behaviors were evaluated using Viewpoint zebrafish tracking software and sexual traits were thoroughly characterized based on anatomical and histological experiments in these KOs and wild types. The steroid hormone levels (testosterone, 11-ketotestosterone and 17β-estradiol) in the brains, gonads, and serum were measured using ELISA kits. To achieve a higher resolution view of the differences in region-specific expression patterns of the brain, the brains of these KOs, and control male and female fish were dissected into three regions: the forebrain, midbrain, and hindbrain for transcriptomic analysis.

    Results

    Qualitative analysis of mating behaviors demonstrated thattdrd12−/−fish behaved in the same manner as wild-type males to trigger oviposition behavior, whilecyp17a1−/−and double knockout (KO) fish did not exhibit these behaviors. Based on the observation of sex characteristics, mating behaviors and hormone levels in these mutants, we found that the maintenance of secondary sex characteristics and male mating behavior did not depend on the presence of germ cells; rather, they depended mainly on the 11-ketotestosterone and testosterone levels secreted into the brain–gonad regulatory axis. RNA-seq analysis of different brain regions revealed that the brain transcript profile oftdrd12−/−fish was similar to that of wild-type males, especially in the forebrain and midbrain. However, the brain transcript profiles ofcyp17a1−/−and double KO fish were distinct from those of wild-type males and were partially biased towards the expression pattern of the female brain. Our results revealed important candidate genes and signaling pathways, such as synaptic signaling/neurotransmission, MAPK signaling, and steroid hormone pathways, that shape brain dimorphism and modulate male mating behavior in zebrafish.

    Conclusions

    Our results provide comprehensive analyses and new insights regarding the endogenous interactions in the brain–gonad–behavior axis. Moreover, this study revealed the crucial candidate genes and neural signaling pathways of different brain regions that are involved in modulating brain dimorphism and male mating behavior in zebrafish, which would significantly light up the understanding the neuroendocrine and molecular mechanisms modulating brain dimorphism and male mating behavior in zebrafish and other teleost fish.

    Graphical Abstract 
    more » « less
  5. Hom, Erik F. (Ed.)
    ABSTRACT <p>Microbiomes have gained significant attention in ecological research, owing to their diverse interactions and essential roles within different organismal ecosystems. Microorganisms, such as bacteria, archaea, and viruses, have profound impact on host health, influencing digestion, metabolism, immune function, tissue development, and behavior. This study investigates the microbiome diversity and function of Kellet’s whelk (<italic>Kelletia kelletii</italic>) perivitelline fluid (PVF), which sustains thousands of developing<italic>K. kelletii</italic>embryos within a polysaccharide and protein matrix. Our core microbiome analysis reveals a diverse range of bacteria, with the<italic>Roseobacter</italic>genus being the most abundant. Additionally, genes related to host-microbe interactions, symbiosis, and quorum sensing were detected, indicating a potential symbiotic relationship between the microbiome and Kellet’s whelk embryos. Furthermore, the microbiome exhibits gene expression related to antibiotic biosynthesis, suggesting a defensive role against pathogenic bacteria and potential discovery of novel antibiotics. Overall, this study sheds light on the microbiome’s role in Kellet’s whelk development, emphasizing the significance of host-microbe interactions in vulnerable life history stages. To our knowledge, ours is the first study to use 16S sequencing coupled with RNA sequencing (RNA-seq) to profile the microbiome of an invertebrate PVF.</p><sec><title>IMPORTANCE

    This study provides novel insight to an encapsulated system with strong evidence of symbiosis between the microbial inhabitants and developing host embryos. The Kellet’s whelk perivitelline fluid (PVF) contains microbial organisms of interest that may be providing symbiotic functions and potential antimicrobial properties during this vulnerable life history stage. This study, the first to utilize a comprehensive approach to investigating Kellet’s whelk PVF microbiome, couples 16S rRNA gene long-read sequencing with RNA-seq. This research contributes to and expands our knowledge on the roles of beneficial host-associated microbes.

     
    more » « less