Ribosomal RNA-Depletion Provides an Efficient Method for Successful Dual RNA-Seq Expression Profiling of a Marine Sponge Holobiont

Investigations of host-symbiont interactions can benefit enormously from a complete and reliable holobiont gene expression profiling. The most efficient way to acquire holobiont transcriptomes is to perform RNA-Seq on both host and symbionts simultaneously. However, optimal methods for capturing both host and symbiont mRNAs are still under development, particularly when the host is a eukaryote and the symbionts are bacteria or archaea. Traditionally, poly(A)-enriched libraries have been used to capture eukaryotic mRNA, but the ability of this method to adequately capture bacterial mRNAs is unclear because of the short half-life of the bacterial transcripts. Here, we address this gap in knowledge with the aim of helping others to choose an appropriate RNA-Seq approach for analysis of animal host-bacterial symbiont transcriptomes. Specifically, we compared transcriptome bias, depth and coverage achieved by two different mRNA capture and sequencing strategies applied to the marine demosponge Amphimedon queenslandica holobiont. Annotated genomes of the sponge host and the three most abundant bacterial symbionts, which can comprise up to 95% of the adult microbiome, are available. Importantly, this allows for transcriptomes to be accurately mapped to these genomes, and thus quantitatively assessed and compared. The two strategies that we compare here are (i) poly(A) captured mRNA-Seq (Poly(A)-RNA-Seq) and (ii) ribosomal RNA depleted RNA-Seq (rRNA-depleted-RNA-Seq). For the host sponge, we find no significant difference in transcriptomes generated by the two different mRNA capture methods. However, for the symbiont transcriptomes, we confirm the expectation that the rRNA-depleted-RNA-Seq performs much better than the Poly(A)-RNA-Seq. This comparison demonstrates that RNA-Seq by ribosomal RNA depletion is an effective and reliable method to simultaneously capture gene expression in host and symbionts and thus to analyse holobiont transcriptomes.

Mitochondrial evolution in the Demospongiae (Porifera): Phylogeny, divergence time, and genome biology

The sponge class Demospongiae is the most speciose and morphologically diverse in the phylum Porifera, and the species within it are vital components of a range of ecosystems worldwide. Despite their ubiquity, a number of recalcitrant problems still remain to be solved regarding their phylogenetic inter-relationships, the timing of their appearance, and their mitochondrial biology, the latter of which is only beginning to be investigated. Here we generated 14 new demosponge mitochondrial genomes which, alongside previously published mitochondrial resources, were used to address these issues. In addition to phylogenomic analysis, we have used syntenic data and analysis of coding regions to forge a framework for understanding the inter-relationships between Demospongiae sub-classes and orders. We have also leveraged our new resources to study the mitochondrial biology of these clades in terms of codon usage, optimisation and gene expression, to understand how these vital cellular components may have contributed to the success of the Porifera. Our results strongly support a sister relationship between Keratosa and (Verongimorpha + Heteroscleromorpha), contradicting previous studies using nuclear markers. Our study includes one species of Clionaida, and show for the first time support for a grouping of Suberitida+(Clionaida+(Tethyida + Poecilosclerida). The findings of our phylogenetic analyses are supported by in-depth examination of structural and coding-level evidence from our mitochondrial data. A time-calibrated phylogeny estimated the origin of Demospongiae in the Cambrian (~529 Mya), and suggests that most demosponge order crown-groups emerged in the Mesozoic. This work therefore provides a robust basis for considering demosponge phylogenetic relationships, as well as essential mitochondrial data for understanding the biological basis for their success and diversity.

Diverse RNA interference strategies in early-branching metazoans

BACKGROUND

Micro RNAs (miRNAs) and piwi interacting RNAs (piRNAs), along with the more ancient eukaryotic endogenous small interfering RNAs (endo-siRNAs) constitute the principal components of the RNA interference (RNAi) repertoire of most animals. RNAi in non-bilaterians - sponges, ctenophores, placozoans and cnidarians - appears to be more diverse than that of bilaterians, and includes structurally variable miRNAs in sponges, an enormous number of piRNAs in cnidarians and the absence of miRNAs in ctenophores and placozoans.

RESULTS

Here we identify thousands of endo-siRNAs and piRNAs from the sponge Amphimedon queenslandica, the ctenophore Mnemiopsis leidyi and the cnidarian Nematostella vectensis using a computational approach that clusters mapped small RNA sequences and annotates each cluster based on the read length and relative abundance of the constituent reads. This approach was validated on 11 small RNA libraries in Drosophila melanogaster, demonstrating the successful annotation of RNAi-associated loci with properties consistent with previous reports. In the non-bilaterians we uncover seven new miRNAs from Amphimedon and four from Nematostella as well as sub-populations of candidate cis-natural antisense transcript (cis-NAT) endo-siRNAs. We confirmed the absence of miRNAs in Mnemiopsis but detected an abundance of endo-siRNAs in this ctenophore. Analysis of putative piRNA structure suggests that conserved localised secondary structures in primary transcripts may be important for the production of mature piRNAs in Amphimedon and Nematostella, as is also the case for endo-siRNAs.

CONCLUSION

Together, these findings suggest that the last common ancestor of extant animals did not have the entrained RNAi system that typifies bilaterians. Instead it appears that bilaterians, cnidarians, ctenophores and sponges express unique repertoires and combinations of miRNAs, piRNAs and endo-siRNAs.

The complete mitochondrial genome of sponge Halichondria okadai (Demospongiae, Suberitida, Halichondriidae) from Korea water

The mitogenome sequence of Sponge Halichondria okadai (Kadota, 1922) (Suberitida, Halichondriidae) was determined for the first time in this study. The circular genome is 20,722 bp in length, containing 14 protein coding genes (PCGs), two ribosomal RNAs (rRNAs), and 25 transfer RNAs (tRNAs). The nucleotide composition of mitogenome consists of 29.5% A, 14.2% C, 21.5% G, 34.7% T, showing a high content of A + T similar to the other Suberitid sponges. These results will be useful for inferring the phylogenetic relationships among the members of family Halichondriidae within the Suberitids.

The complete mitochondrial genome of sponge Halichondria (Halichondria) sp. (Demospongiae, Suberitida, Halichondriidae)

In this study, we found out the complete mitochondrial genome of Halichondria (Halichondria) sp., a common demosponge in China. This is the first complete mitochondrial report on the genus Halichondria. The mitochondrial genome of Halichondria (Halichondria) sp. is 20 746 bp in length, with 14 protein-coding genes, two rRNA genes and 25 tRNA genes. When compared with the complete mitochondrial genome of Hymeniacidon sinapium, our phylogenetic result suggested that Halichondria (Halichondria) sp. converged well according to morphological result.

The complete mitochondrial genome of sponge Tethya sp. (Demospongiae, Tethyida, Tethyidae)

The complete mitochondrial genome of Tethya sp. was studied. This is the second complete mitochondrial report on the family Tethyidae. The mitochondrial genome of Tethya sp. is 20,582 bp in length, containing 14 protein-coding genes and 25 tRNA genes, with 2 rRNA genes. Our phylogenetic result suggested that Tethya sp. converged well with Tethya actinia, which further verified the morphological result. We anticipate our study to shed light on future molecular studies of demosponges.