De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide.

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

Feminising Wolbachia disrupt Armadillidium vulgare insulin-like signalling pathway.

The endosymbiont Wolbachia feminises male isopods by making them refractory to the insulin-like masculinising hormone, which shunts the autocrine development of the androgenic glands. It was, therefore, proposed that Wolbachia silences the IR receptors, either by preventing their expression or by inactivating them. We describe here the two IR paralogs of Armadillidium vulgare. They displayed a conventional structure and belonged to a family widespread among isopods. Av-IR1 displayed an ubiquist expression, whereas the expression of Av-IR2 was restricted to the gonads. Both were constitutively expressed in males and females and throughout development. However, upon silencing, altered gland physiology and gene expression therein suggested antagonistic roles for Av-IR1 (androinhibiting) and Av-IR2 (androstimulating). They may function in tandem with regulating neurohormones, as a conditional platform that conveys insulin signalling. Wolbachia infection did not alter their expression patterns: leaving the IRs unscathed, the bacteria would suppress the secretion of the neurohormones, thus inducing body-wide IR deactivation and feminisation. Adult males injected with Wolbachia acquired an intersexed physiology. Their phenotypes and gene expressions mirrored the silencing of Av-IR1 only, suggesting that imperfect feminisation stems from a flawed invasion of the androstimulating centre, whereas in fully feminised males invasion would be complete in early juveniles. TAKE AWAY: Two antagonistic Insulin Receptors were characterised in Armadillidium vulgare. The IRs were involved in androstimulating and androinhibiting functions. Wolbachia-induced feminisation did not prevent the expression of the IRs. Imperfectly feminised intersexes phenocopied the silencing of Av-IR1 only. Wolbachia would deactivate the IRs by suppressing neurosecretory co-factors.

IGFBP-rP1, a strongly conserved member of the androgenic hormone signalling pathway in Isopoda.

The first protein which has been described to interact with the malacostracan Androgenic Gland Hormone (AGH) is a binding protein called IGFBP-rP1. It has been identified and studied in several species of decapods, in which its interaction with the masculinizing hormone and its expression patterns have been established in several ways. However, this protein remains uncharacterised to date in the other malacostracan orders, like Amphipoda and Isopoda, although they were historically the first ones in which the androgenic gland and the corresponding hormone were respectively described. In this article, we identified the IGFBP-rP1 of isopods and established its implication in the pathway of the AGH with a silencing approach in the model species Armadillidium vulgare. We also showed that this gene is expressed in all the tissues of males and females, with a similar pattern in animals infected with Wolbachia, a feminizing endosymbiont of several isopod species. The expression pattern did not differ during the development of uninfected and infected animals either. We finally studied the evolution of the IGFBP-rP1 in 68 isopod species, looking for conserved motifs and evidence of natural selection. Altogether, our results showed that this gene is constitutively expressed and strongly conserved in isopods, in which it likely constitutes a key element of the insulin/IGF signalling pathway. However, we also illustrated that IGFBP-rP1 is not sufficient on its own to explain the different developmental paths taken by the males and the females or feminized genetic males.

Identification and validation of reference genes for qPCR in the terrestrial gastropod Cepaea nemoralis.

Identifying and understanding mechanisms that generate phenotypic diversity is a fundamental goal of evolutionary biology. With a diversity of pigmented shell morphotypes governed by Mendelian patterns of inheritance, the common grove snail Cepaea nemoralis (Linnaeus, 1758) has been a model for evolutionary biologists and population geneticists for decades. However, the genetic mechanisms by which C. nemoralis generates this pigmented shell diversity remain unknown. An important first step in investigating this pigmentation pattern is to establish a set of validated reference genes for differential gene expression assays. Here we have evaluated eleven candidate genes for reverse transcription quantitative polymerase chain reaction (qPCR) in C. nemoralis. Five of these were housekeeping genes traditionally employed as qPCR reference genes in other species, while six alternative genes were selected de novo from C. nemoralis transcriptome data based on the stability of their expression levels. We tested all eleven candidates for expression stability in four sub-adult tissues of C. nemoralis: pigmented mantle, unpigmented mantle, head and foot. We find that two commonly employed housekeeping genes (alpha tubulin, glyceraldehyde 3-phosphate dehydrogenase) are unsuitable for use as qPCR reference genes in C. nemoralis. The traditional housekeeping gene UBIquitin on the other hand performed very well. Additionally, an RNA-directed DNA polymerase (RNAP), a Potassium Channel Protein (KCHP) and a Prenylated Rab acceptor protein 1 (PRAP), identified de novo from transcriptomic data, were the most stably expressed genes in different tissue combinations. We also tested expression stability over two seasons and found that, although other genes are more stable within a single season, beta actin (BACT) and elongation factor 1 alpha (EF1α) were the most reliable reference genes across seasons.

In-depth proteomic analyses of Haliotis laevigata (greenlip abalone) nacre and prismatic organic shell matrix.

BACKGROUND: The shells of various Haliotis species have served as models of invertebrate biomineralization and physical shell properties for more than 20 years. A focus of this research has been the nacreous inner layer of the shell with its conspicuous arrangement of aragonite platelets, resembling in cross-section a brick-and-mortar wall. In comparison, the outer, less stable, calcitic prismatic layer has received much less attention. One of the first molluscan shell proteins to be characterized at the molecular level was Lustrin A, a component of the nacreous organic matrix of Haliotis rufescens. This was soon followed by the C-type lectin perlucin and the growth factor-binding perlustrin, both isolated from H. laevigata nacre, and the crystal growth-modulating AP7 and AP24, isolated from H. rufescens nacre. Mass spectrometry-based proteomics was subsequently applied to to Haliotis biomineralization research with the analysis of the H. asinina shell matrix and yielded 14 different shell-associated proteins. That study was the most comprehensive for a Haliotis species to date. METHODS: The shell proteomes of nacre and prismatic layer of the marine gastropod Haliotis laevigata were analyzed combining mass spectrometry-based proteomics and next generation sequencing. RESULTS: We identified 297 proteins from the nacreous shell layer and 350 proteins from the prismatic shell layer from the green lip abalone H. laevigata. Considering the overlap between the two sets we identified a total of 448 proteins. Fifty-one nacre proteins and 43 prismatic layer proteins were defined as major proteins based on their abundance at more than 0.2% of the total. The remaining proteins occurred at low abundance and may not play any significant role in shell fabrication. The overlap of major proteins between the two shell layers was 17, amounting to a total of 77 major proteins. CONCLUSIONS: The H. laevigata shell proteome shares moderate sequence similarity at the protein level with other gastropod, bivalve and more distantly related invertebrate biomineralising proteomes. Features conserved in H. laevigata and other molluscan shell proteomes include short repetitive sequences of low complexity predicted to lack intrinsic three-dimensional structure, and domains such as tyrosinase, chitin-binding, and carbonic anhydrase. This catalogue of H. laevigata shell proteins represents the most comprehensive for a haliotid and should support future efforts to elucidate the molecular mechanisms of shell assembly.

Variation in Orthologous Shell-Forming Proteins Contribute to Molluscan Shell Diversity.

Despite the evolutionary success and ancient heritage of the molluscan shell, little is known about the molecular details of its formation, evolutionary origins, or the interactions between the material properties of the shell and its organic constituents. In contrast to this dearth of information, a growing collection of molluscan shell-forming proteomes and transcriptomes suggest they are comprised of both deeply conserved, and lineage specific elements. Analyses of these sequence data sets have suggested that mechanisms such as exon shuffling, gene co-option, and gene family expansion facilitated the rapid evolution of shell-forming proteomes and supported the diversification of this phylum specific structure. In order to further investigate and test these ideas we have examined the molecular features and spatial expression patterns of two shell-forming genes (Lustrin and ML1A2) and coupled these observations with materials properties measurements of shells from a group of closely related gastropods (abalone). We find that the prominent "GS" domain of Lustrin, a domain believed to confer elastomeric properties to the shell, varies significantly in length between the species we investigated. Furthermore, the spatial expression patterns of Lustrin and ML1A2 also vary significantly between species, suggesting that both protein architecture, and the regulation of spatial gene expression patterns, are important drivers of molluscan shell evolution. Variation in these molecular features might relate to certain materials properties of the shells of these species. These insights reveal an important and underappreciated source of variation within shell-forming proteomes that must contribute to the diversity of molluscan shell phenotypes.

Combining independent de novo assemblies optimizes the coding transcriptome for nonconventional model eukaryotic organisms.

BACKGROUND: Next-generation sequencing (NGS) technologies are arguably the most revolutionary technical development to join the list of tools available to molecular biologists since PCR. For researchers working with nonconventional model organisms one major problem with the currently dominant NGS platform (Illumina) stems from the obligatory fragmentation of nucleic acid material that occurs prior to sequencing during library preparation. This step creates a significant bioinformatic challenge for accurate de novo assembly of novel transcriptome data. This challenge becomes apparent when a variety of modern assembly tools (of which there is no shortage) are applied to the same raw NGS dataset. With the same assembly parameters these tools can generate markedly different assembly outputs. RESULTS: In this study we present an approach that generates an optimized consensus de novo assembly of eukaryotic coding transcriptomes. This approach does not represent a new assembler, rather it combines the outputs of a variety of established assembly packages, and removes redundancy via a series of clustering steps. We test and validate our approach using Illumina datasets from six phylogenetically diverse eukaryotes (three metazoans, two plants and a yeast) and two simulated datasets derived from metazoan reference genome annotations. All of these datasets were assembled using three currently popular assembly packages (CLC, Trinity and IDBA-tran). In addition, we experimentally demonstrate that transcripts unique to one particular assembly package are likely to be bioinformatic artefacts. For all eight datasets our pipeline generates more concise transcriptomes that in fact possess more unique annotatable protein domains than any of the three individual assemblers we employed. Another measure of assembly completeness (using the purpose built BUSCO databases) also confirmed that our approach yields more information. CONCLUSIONS: Our approach yields coding transcriptome assemblies that are more likely to be closer to biological reality than any of the three individual assembly packages we investigated. This approach (freely available as a simple perl script) will be of use to researchers working with species for which there is little or no reference data against which the assembly of a transcriptome can be performed.

Genomic context drives transcription of insertion sequences in the bacterial endosymbiont Wolbachia wVulC.

Transposable elements (TEs) are DNA pieces that are present in almost all the living world at variable genomic density. Due to their mobility and density, TEs are involved in a large array of genomic modifications. In eukaryotes, TE expression has been studied in detail in several species. In prokaryotes, studies of IS expression are generally linked to particular copies that induce a modification of neighboring gene expression. Here we investigated global patterns of IS transcription in the Alphaproteobacterial endosymbiont Wolbachia wVulC, using both RT-PCR and bioinformatic analyses. We detected several transcriptional promoters in all IS groups. Nevertheless, only one of the potentially functional IS groups possesses a promoter located upstream of the transposase gene, that could lead up to the production of a functional protein. We found that the majority of IS groups are expressed whatever their functional status. RT-PCR analyses indicate that the transcription of two IS groups lacking internal promoters upstream of the transposase start codon may be driven by the genomic environment. We confirmed this observation with the transcription analysis of individual copies of one IS group. These results suggest that the genomic environment is important for IS expression and it could explain, at least partly, copy number variability of the various IS groups present in the wVulC genome and, more generally, in bacterial genomes.

Molecular evolution of the androgenic hormone in terrestrial isopods.

In crustaceans, the androgenic gland (AG), thanks to the synthesis of the androgenic gland hormone (AGH), controls the differentiation of the primary and secondary male sexual characters. In this study, we amplified 12 new AGH cDNAs in species belonging to five different families of the infra-order Ligiamorpha of terrestrial isopods. Putative essential amino acids for the production of a functional AGH protein exhibit signatures of negative selection and are strictly conserved including typical proteolytic cleavage motifs, a putative N-linked glycosylation motif on the A chains and the eight Cys positions. An insulin-like growth factor motif was also identified in Armadillidium AGH sequences. The phylogenetic relationships of AGH sequences allowed one to distinguish two main clades, corresponding to members of the Armadillidiidae and the Porcellionidae families which are congruent with the narrow specificity of AG heterospecific grafting. An in-depth understanding of the regulation of AGH expression would help deciphering the interaction between Wolbachia, widespread feminizing endosymbiotic bacteria in isopods, and the sex differentiation of their hosts.

Horizontal transfer of transposons between and within crustaceans and insects.

BACKGROUND: Horizontal transfer of transposable elements (HTT) is increasingly appreciated as an important source of genome and species evolution in eukaryotes. However, our understanding of HTT dynamics is still poor in eukaryotes because the diversity of species for which whole genome sequences are available is biased and does not reflect the global eukaryote diversity. RESULTS: In this study we characterized two Mariner transposable elements (TEs) in the genome of several terrestrial crustacean isopods, a group of animals particularly underrepresented in genome databases. The two elements have a patchy distribution in the arthropod tree and they are highly similar (>93% over the entire length of the element) to insect TEs (Diptera and Hymenoptera), some of which were previously described in Ceratitis rosa (Crmar2) and Drosophila biarmipes (Mariner-5_Dbi). In addition, phylogenetic analyses and comparisons of TE versus orthologous gene distances at various phylogenetic levels revealed that the taxonomic distribution of the two elements is incompatible with vertical inheritance. CONCLUSIONS: We conclude that the two Mariner TEs each underwent at least three HTT events. Both elements were transferred once between isopod crustaceans and insects and at least once between isopod crustacean species. Crmar2 was also transferred between tephritid and drosophilid flies and Mariner-5 underwent HT between hymenopterans and dipterans. We demonstrate that these various HTTs took place recently (most likely within the last 3 million years), and propose iridoviruses and/or Wolbachia endosymbionts as potential vectors of these transfers.

Short- and long-term evolutionary dynamics of bacterial insertion sequences: insights from Wolbachia endosymbionts.

Transposable elements (TE) are one of the major driving forces of genome evolution, raising the question of the long-term dynamics underlying their evolutionary success. Long-term TE evolution can readily be reconstructed in eukaryotes, thanks to many degraded copies constituting genomic fossil records of past TE proliferations. By contrast, bacterial genomes usually experience high sequence turnover and short TE retention times, thereby obscuring ancient TE evolutionary patterns. We found that Wolbachia bacterial genomes contain 52-171 insertion sequence (IS) TEs. IS account for 11% of Wolbachia wRi, which is one of the highest IS genomic coverage reported in prokaryotes to date. We show that many IS groups are currently expanding in various Wolbachia genomes and that IS horizontal transfers are frequent among strains, which can explain the apparent synchronicity of these IS proliferations. Remarkably, >70% of Wolbachia IS are nonfunctional. They constitute an unusual bacterial IS genomic fossil record providing direct empirical evidence for a long-term IS evolutionary dynamics following successive periods of intense transpositional activity. Our results show that comprehensive IS annotations have the potential to provide new insights into prokaryote TE evolution and, more generally, prokaryote genome evolution. Indeed, the identification of an important IS genomic fossil record in Wolbachia demonstrates that IS elements are not always of recent origin, contrary to the conventional view of TE evolution in prokaryote genomes. Our results also raise the question whether the abundance of IS fossils is specific to Wolbachia or it may be a general, albeit overlooked, feature of prokaryote genomes.