The Oyster River Protocol: a multi-assembler and kmer approach for de novo transcriptome assembly

Semimytilus algosus: first known hermaphroditic mussel with doubly uniparental inheritance of mitochondrial DNA

Doubly uniparental inheritance (DUI) of mitochondrial DNA is a rare phenomenon occurring in some freshwater and marine bivalves and is usually characterized by the mitochondrial heteroplasmy of male individuals. Previous research on freshwater Unionida mussels showed that hermaphroditic species do not have DUI even if their closest gonochoristic counterparts do. No records showing DUI in a hermaphrodite have ever been reported. Here we show for the first time that the hermaphroditic mussel Semimytilus algosus (Mytilida), very likely has DUI, based on the complete sequences of both mitochondrial DNAs and the distribution of mtDNA types between male and female gonads. The two mitogenomes show considerable divergence (34.7%). The presumably paternal M type mitogenome dominated the male gonads of most studied mussels, while remaining at very low or undetectable levels in the female gonads of the same individuals. If indeed DUI can function in the context of simultaneous hermaphroditism, a change of paradigm regarding its involvement in sex determination is needed. It is apparently associated with gonadal differentiation rather than with sex determination in bivalves.

Liver transcriptome resources of four commercially exploited teleost species

The generation of omic resources is central to develop adequate management strategies for species with economic value. Here, we provide high-coverage RNA-seq datasets of liver tissue (containing between 80,2 and 88,4 million of paired-end reads) from four wildtype teleost species with high commercial value: Trachurus trachurus (TTR; Atlantic horse mackerel), Scomber scombrus (SSC; Atlantic mackerel), Trisopterus luscus (TLU; pout), and Micromesistius poutassou (MPO; blue whiting). A comprehensive assembly pipeline, using de novo single and multi-kmer assembly approaches, produced 64 single high-quality liver transcriptomes - 16 per species. The final assemblies, with N50 values ranging from 2,543-3,700 bp and BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness values between 81.8-86.5% of the Actinopterygii gene set, were subjected to open reading frame (ORF) prediction and functional annotation. Our study provides the first transcriptomic resources for these species and offers valuable tools to evaluate both neutral and selected genetic variation among populations, and to identify candidate genes for environmental adaptation assisting in the investigation of the effects of global changes in fisheries.

De novo assembly and functional annotation of the heart + hemolymph transcriptome in the Caribbean spiny lobster Panulirus argus

The spiny lobster, Panulirus argus, is an ecologically relevant species in shallow water coral reefs and a target of the most lucrative fishery in the greater Caribbean region. This study reports, for the first time, the heart + hemolymph transcriptome of the Caribbean spiny lobster Panulirus argus assembled from short Illumina 150bp PE raw reads. A total of 80,152,094 raw reads were assembled using the Oyster River Protocol pipeline. The assembly resulted in a total of 254,773 transcripts. Functional gene annotation was conducted using the software package 'dammit'. Lastly, gene enrichment analyses were conducted using the Gene Ontology (GO) and KEGG pathway (Kaas) databases. This resource will be of utmost importance in future research aiming at exploring the effect of local and regional anthropogenic disturbances, as well as global climate change on the molecular physiology of this overexploited species.

Error, noise and bias in de novo transcriptome assemblies

De novo transcriptome assembly is a powerful tool, and has been widely used over the last decade for making evolutionary inferences. However, it relies on two implicit assumptions: that the assembled transcriptome is an unbiased representation of the underlying expressed transcriptome, and that expression estimates from the assembly are good, if noisy approximations of the relative abundance of expressed transcripts. Using publicly available data for model organisms, we demonstrate that, across assembly algorithms and data sets, these assumptions are consistently violated. Bias exists at the nucleotide level, with genotyping error rates ranging from 30% to 83%. As a result, diversity is underestimated in transcriptome assemblies, with consistent underestimation of heterozygosity in all but the most inbred samples. Even at the gene level, expression estimates show wide deviations from map-to-reference estimates, and positive bias at lower expression levels. Standard filtering of transcriptome assemblies improves the robustness of gene expression estimates but leads to the loss of a meaningful number of protein-coding genes, including many that are highly expressed. We demonstrate a computational method, length-rescaled CPM, to partly alleviate noise and bias in expression estimates. Researchers should consider ways to minimize the impact of bias in transcriptome assemblies.

Comparative Genomic Characterization of the Multimammate Mouse Mastomys coucha

Mastomys are the most widespread African rodent and carriers of various diseases such as the plague or Lassa virus. In addition, mastomys have rapidly gained a large number of mammary glands. Here, we generated a genome, variome, and transcriptomes for Mastomys coucha. As mastomys diverged at similar times from mouse and rat, we demonstrate their utility as a comparative genomic tool for these commonly used animal models. Furthermore, we identified over 500 mastomys accelerated regions, often residing near important mammary developmental genes or within their exons leading to protein sequence changes. Functional characterization of a noncoding mastomys accelerated region, located in the HoxD locus, showed enhancer activity in mouse developing mammary glands. Combined, our results provide genomic resources for mastomys and highlight their potential both as a comparative genomic tool and for the identification of mammary gland number determining factors.

Dendritic Cells Generated From Mops condylurus, a Likely Filovirus Reservoir Host, Are Susceptible to and Activated by Zaire Ebolavirus Infection

Ebola virus infection of human dendritic cells (DCs) induces atypical adaptive immune responses and thereby exacerbates Ebola virus disease (EVD). Human DCs, infected with Ebola virus aberrantly express low levels of the DC activation markers CD80, CD86, and MHC class II. The T cell responses ensuing are commonly anergic rather than protective against EVD. We hypothesize that DCs derived from potential reservoir hosts such as bats, which do not develop disease signs in response to Ebola virus infection, would exhibit features associated with activation. In this study, we have examined Zaire ebolavirus (EBOV) infection of DCs derived from the Angolan free-tailed bat species, Mops condylurus. This species was previously identified as permissive to EBOV infection in vivo, in the absence of disease signs. M. condylurus has also been recently implicated as the reservoir host for Bombali ebolavirus, a virus species that is closely related to EBOV. Due to the absence of pre-existing M. condylurus species-specific reagents, we characterized its de novo assembled transcriptome and defined its phylogenetic similarity to other mammals, which enabled the identification of cross-reactive reagents for M. condylurus bone marrow-derived DC (bat-BMDC) differentiation and immune cell phenotyping. Our results reveal that bat-BMDCs are susceptible to EBOV infection as determined by detection of EBOV specific viral RNA (vRNA). vRNA increased significantly 72 h after EBOV-infection and was detected in both cells and in culture supernatants. Bat-BMDC infection was further confirmed by the observation of GFP expression in DC cultures infected with a recombinant GFP-EBOV. Bat-BMDCs upregulated CD80 and chemokine ligand 3 (CCL3) transcripts in response to EBOV infection, which positively correlated with the expression levels of EBOV vRNA. In contrast to the aberrant responses to EBOV infection that are typical for human-DC, our findings from bat-BMDCs provide evidence for an immunological basis of asymptomatic EBOV infection outcomes.

Evaluating the performance of targeted sequence capture, RNA-Seq, and degenerate-primer PCR cloning for sequencing the largest mammalian multigene family

Multigene families evolve from single-copy ancestral genes via duplication, and typically encode proteins critical to key biological processes. Molecular analyses of these gene families require high-confidence sequences, but the high sequence similarity of the members can create challenges for sequencing and downstream analyses. Focusing on the common vampire bat, Desmodus rotundus, we evaluated how different sequencing approaches performed in recovering the largest mammalian protein-coding multigene family: olfactory receptors (OR). Using the genome as a reference, we determined the proportion of intact protein-coding receptors recovered by: (a) amplicons from degenerate primers sequenced via Sanger technology, (b) RNA-Seq of the main olfactory epithelium, and (c) those genes captured with probes designed from transcriptomes of closely-related species. Our initial re-annotation of the high-quality vampire bat genome resulted in >400 intact OR genes, more than doubling the original estimate. Sanger-sequenced amplicons performed the poorest among the three approaches, detecting <33% of receptors in the genome. In contrast, the transcriptome reliably recovered >50% of the annotated genomic ORs, and targeted sequence capture recovered nearly 75% of annotated genes. Each sequencing approach assembled high-quality sequences, even if it did not recover all receptors in the genome. While some variation may be due to limitations of the study design (e.g., different individuals), variation among approaches was mostly caused by low coverage of some receptors rather than high rates of assembly error. Given this variability, we caution against using the counts of intact receptors per species to model the birth-death process of multigene families. Instead, our results support the use of orthologous sequences to explore and model the evolutionary processes shaping these genes.

Expressed Vomeronasal Type-1 Receptors (V1rs) in Bats Uncover Conserved Sequences Underlying Social Chemical Signaling

In mammals, social and reproductive behaviors are mediated by chemical cues encoded by hyperdiverse families of receptors expressed in the vomeronasal organ. Between species, the number of intact receptors can vary by orders of magnitude. However, the evolutionary processes behind variation in receptor number, and its link to fitness-related behaviors are not well understood. From vomeronasal transcriptomes, we discovered the first evidence of intact vomeronasal type-1 receptor (V1r) genes in bats, and we tested whether putatively functional bat receptors were orthologous to those of related taxa, or whether bats have evolved novel receptors. Instead of lineage-specific duplications, we found that bat V1rs show high levels of orthology to those of their relatives, and receptors are under comparative levels of purifying selection as non-bats. Despite widespread vomeronasal organ loss in bats, V1r copies have been retained for >65 million years. The highly conserved nature of bat V1rs challenges our current understanding of mammalian V1r function and suggests roles other than conspecific recognition or mating initiation in social behavior.

Variation in pigmentation gene expression is associated with distinct aposematic color morphs in the poison frog Dendrobates auratus

Background

Color and pattern phenotypes have clear implications for survival and reproduction in many species. However, the mechanisms that produce this coloration are still poorly characterized, especially at the genomic level. Here we have taken a transcriptomics-based approach to elucidate the underlying genetic mechanisms affecting color and pattern in a highly polytypic poison frog. We sequenced RNA from the skin from four different color morphs during the final stage of metamorphosis and assembled a de novo transcriptome. We then investigated differential gene expression, with an emphasis on examining candidate color genes from other taxa.

Results

Overall, we found differential expression of a suite of genes that control melanogenesis, melanocyte differentiation, and melanocyte proliferation (e.g., tyrp1, lef1, leo1, and mitf) as well as several differentially expressed genes involved in purine synthesis and iridophore development (e.g., arfgap1, arfgap2, airc, and gart).

Conclusions

Our results provide evidence that several gene networks known to affect color and pattern in vertebrates play a role in color and pattern variation in this species of poison frog.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1410-7) contains supplementary material, which is available to authorized users.

Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes

BACKGROUND

De novo transcriptome assemblies are required prior to analyzing RNA sequencing data from a species without an existing reference genome or transcriptome. Despite the prevalence of transcriptomic studies, the effects of using different workflows, or "pipelines," on the resulting assemblies are poorly understood. Here, a pipeline was programmatically automated and used to assemble and annotate raw transcriptomic short-read data collected as part of the Marine Microbial Eukaryotic Transcriptome Sequencing Project. The resulting transcriptome assemblies were evaluated and compared against assemblies that were previously generated with a different pipeline developed by the National Center for Genome Research.

RESULTS

New transcriptome assemblies contained the majority of previous contigs as well as new content. On average, 7.8% of the annotated contigs in the new assemblies were novel gene names not found in the previous assemblies. Taxonomic trends were observed in the assembly metrics. Assemblies from the Dinoflagellata showed a higher number of contigs and unique k-mers than transcriptomes from other phyla, while assemblies from Ciliophora had a lower percentage of open reading frames compared to other phyla.

CONCLUSIONS

Given current bioinformatics approaches, there is no single "best" reference transcriptome for a particular set of raw data. As the optimum transcriptome is a moving target, improving (or not) with new tools and approaches, automated and programmable pipelines are invaluable for managing the computationally intensive tasks required for re-processing large sets of samples with revised pipelines and ensuring a common evaluation workflow is applied to all samples. Thus, re-assembling existing data with new tools using automated and programmable pipelines may yield more accurate identification of taxon-specific trends across samples in addition to novel and useful products for the community.

The transcriptome of the veiled chameleon (Chamaeleo calyptratus): A resource for studying the evolution and development of vertebrates

PURPOSE

The veiled chameleon (Chamaeleo calyptratus) is an emerging model system for studying functional morphology and evolutionary developmental biology (evo-devo). Chameleons possess body plans that are highly adapted to an arboreal life style, featuring laterally compressed bodies, split hands/ft for grasping, a projectile tongue, turreted independently moving eyes, and a prehensile tail. Despite being one of the most phenotypically divergent clades of tetrapods, genomic resources for chameleons are severely lacking.

METHODS

To address this lack of resources, we used RNAseq to generate 288 million raw Illumina sequence reads from four adult tissues (male and female eyes and gonads) and whole embryos at three distinct developmental stages. We used these data to assemble a largely complete de novo transcriptome consisting of only 82 952 transcripts. In addition, a majority of assembled transcripts (67%) were successfully annotated.

RESULTS

We then demonstrated the utility of these data in the context of studying visual system evolution by examining the content of veiled chameleon opsin genes to show that chameleons possess all five ancestral tetrapod opsins.

CONCLUSION

We present this de novo, annotated, multi-tissue transcriptome assembly for the Veiled Chameleon, Chamaeleo calyptratus, as a resource to address a range of evolutionary and developmental questions. The associated raw reads and final annotated transcriptome assembly are freely available for use on NCBI and Figshare, respectively.

Keeping it light: (re)analyzing community-wide datasets without major infrastructure

DNA sequencing technology has revolutionized the field of biology, shifting biology from a data-limited to data-rich state. Central to the interpretation of sequencing data are the computational tools and approaches that convert raw data into biologically meaningful information. Both the tools and the generation of data are actively evolving, yet the practice of re-analysis of previously generated data with new tools is not commonplace. Re-analysis of existing data provides an affordable means of generating new information and will likely become more routine within biology, yet necessitates a new set of considerations for best practices and resource development. Here, we discuss several practices that we believe to be broadly applicable when re-analyzing data, especially when done by small research groups.