Comparative description of ten transcriptomes of newly sequenced invertebrates and efficiency estimation of genomic sampling in non-model taxa

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

Phylogenomic analyses shed light on the relationships of chiton superfamilies and shell-eye evolution

Mollusca is the second-largest animal phylum with over 100,000 extant species representing eight classes. Across 1000 extant species in the class Polyplacophora, chitons have a relatively constrained morphology but with some notable deviations. Several genera possess "shell eyes", i.e., true eyes with a lens and retina that are embedded within the dorsal shells. The phylogeny of the major chiton clades is mostly well established, in a set of superfamily-level and higher level taxa supported by various approaches, including morphological studies, multiple gene markers, mitogenome-phylogeny, and phylotranscriptomic approaches. However, one critical lineage has remained unclear, namely Schizochiton which was controversially suggested as being the potential independent origin of chiton shell eyes. Here, with the draft genome sequencing of Schizochiton incisus (superfamily Schizochitonoidea) plus assemblies of transcriptome data from other polyplacophorans, we present phylogenetic reconstructions using both mitochondrial genomes and phylogenomic approaches with multiple methods. We found that phylogenetic trees from mitogenomic data are inconsistent, reflecting larger scale confounding factors in molluscan mitogenomes. However, a consistent and robust topology was generated with protein-coding genes using different models and methods. Our results support Schizochitonoidea as the sister group to other Chitonoidea in Chitonina, in agreement with the established classification. Combined with evidence from fossils, our phylogenetic results suggest that the earliest origin of shell eyes is in Schizochitonoidea, and that these structures were also gained secondarily in other genera in Chitonoidea. Our results have generated a holistic review of the internal relationship within Polyplacophora, and a better understanding of the evolution of Polyplacophora.

Radula and Shell Microstructure Variations are Congruent with a Molecular Estimate of Shallow-Water Japanese Chitons

Variations of the radula and shell microstructures in 33 species of Japanese chiton were investigated along with molecular phylogenetic trees. The molecular phylogenetic trees indicated that Chitonida was composed of four clades, of which two clades formed Acanthochitonina and corresponded to Mopalioidea and Cryptoplacoidea, respectively, and the other clades formed Chitonina. In the radula, the shapes of the central and centro-lateral teeth and the petaloid process varied greatly among species or genera and were useful for the identification of particular species or genera. The presence of accessory and petaloid processes and the cusp shape were relatively conserved and useful for recognizing particular genera or even suborders. In the valves, four to six shell layers were found at the section, but the ventral mesostracum was not observed in Acanthochitonina. The shell microstructures in the ventral sublayer of the tegmentum varied at suborder, but those in the other layers were almost constant. The megalaesthete chamber type varied at superfamily and was helpful to identify particular families or superfamilies. The characteristics of the shell layers and shell microstructures appear to be a synapomorphy shared by the members of Acanthochitonina. The classification within Chitonina needs to be reexamined because the variations of the cusp shape and megalaesthete chamber type were relatively large and did not correspond to the current classification. Callochiton formed a sister group with Chitonida and would be equally closely related to Chitonina and Acanthochitonina because of possessing a mosaic of characteristics from both.

Third-Generation Sequencing Reveals the Adaptive Role of the Epigenome in Three Deep-Sea Polychaetes

The roles of DNA methylation in invertebrates are poorly characterized, and critical data are missing for the phylum Annelida. We fill this knowledge gap by conducting the first genome-wide survey of DNA methylation in the deep-sea polychaetes dominant in deep-sea vents and seeps: Paraescarpia echinospica, Ridgeia piscesae, and Paralvinella palmiformis. DNA methylation calls were inferred from Oxford Nanopore sequencing after assembling high-quality genomes of these animals. The genomes of these worms encode all the key enzymes of the DNA methylation metabolism and possess a mosaic methylome similar to that of other invertebrates. Transcriptomic data of these polychaetes support the hypotheses that gene body methylation strengthens the expression of housekeeping genes and that promoter methylation acts as a silencing mechanism but not the hypothesis that DNA methylation suppresses the activity of transposable elements. The conserved epigenetic profiles of genes responsible for maintaining homeostasis under extreme hydrostatic pressure suggest DNA methylation plays an important adaptive role in these worms.

Sperm Transcriptome Analysis Accurately Reveals Male Fertility Potential in Livestock

Nowadays, selection of superior male candidates in livestock as a source of frozen semen based on sperm quality at the cellular level is not considered accurate enough for predicting the potential of male fertility. Sperm transcriptome analysis approaches, such as messenger RNA levels, have been shown to correlate with fertility rates. Using this technology in livestock growth has become the principal method, which can be widely applied to predict male fertility potential in the livestock industry through the analysis of the sperm transcriptome. It provides the gene expression to validate the function of sperm in spermatogenesis, fertilization, and embryo development, as the parameters of male fertility. This review proposes a transcriptomic analysis approach as a high-throughput method to predict the fertility potential of livestock more accurately in the future.

Transcriptomic and Histological Analysis of the Greentail Prawn (Metapenaeus bennettae) Following Light Crude Oil Exposure

Oil spills pose a significant threat to marine biodiversity. Crude oil can partition into sediments where it may be persistent, placing benthic species such as decapods at particular risk of exposure. Transcriptomic and histological tools are often used to investigate the effects of hydrocarbon exposure on marine organisms following oil spill events, allowing for the identification of metabolic pathways impacted by oil exposure. However, there is limited information available for decapod crustaceans, many of which carry significant economic value. In the present study, we assess the sublethal impacts of crude oil exposure in the commercially important Australian greentail prawn (Metapenaeus bennettae) using transcriptomic and histological analyses. Prawns exposed to light, unweathered crude oil "spiked" sediments for 90 h were transferred to clean sediments for a further 72 h to assess recovery. Chemical analyses indicated that polycyclic aromatic hydrocarbons increased by approximately 65% and 91% in prawn muscle following 24 and 90 h of exposure, respectively, and significantly decreased during 24- and 72-h recovery periods. Transcriptomic responses followed an exposure and recovery pattern with innate immunity and nutrient metabolism transcripts significantly lowered in abundance after 24 h of exposure and were higher in abundance after 72 h of recovery. In addition, transcription/translation, cellular responses, and DNA repair pathways were significantly impacted after 24 h of exposure and recovered after 72 h of recovery. However, histological alterations such as tubule atrophy indicated an increase in severity after 24 and 72 h of recovery. The present study provides new insights into the sublethal impacts of crude oil exposure in greentail prawns and identifies molecular pathways altered by exposure. We expect these findings to inform future management associated with oil extraction activity and spills. Environ Toxicol Chem 2022;41:2162-2180. © 2022 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Transcriptome Analysis to Study the Molecular Response in the Gill and Hepatopancreas Tissues of Macrobrachium nipponense to Salinity Acclimation

Macrobrachium nipponense is an economically important prawn species and common in Chinese inland capture fisheries. During aquaculture, M. nipponense can survive under freshwater and low salinity conditions. The molecular mechanism underlying the response to salinity acclimation remains unclear in this species; thus, in this study, we used the Illumina RNA sequencing platform for transcriptome analyses of the gill and hepatopancreas tissues of M. nipponense exposed to salinity stress [0.4‰ (S0, control group), 6‰ (S6, low salinity group), and 12‰ (S12, high salinity group)]. Differentially expressed genes were identified, and several important salinity adaptation-related terms and signaling pathways were found to be enriched, such as "ion transport," "oxidative phosphorylation," and "glycometabolism." Quantitative real-time PCR demonstrated the participation of 12 key genes in osmotic pressure regulation in M. nipponense under acute salinity stress. Further, the role of carbonic anhydrase in response to salinity acclimation was investigated by subjecting the gill tissues of M. nipponense to in situ hybridization. Collectively, the results reported herein enhance our understanding of the mechanisms via which M. nipponense adapts to changes in salinity.

Transcriptomic and metabolomic analyses reveal the altitude adaptability and evolution of different-colored flowers in alpine Rhododendron species

Understanding the molecular mechanisms and evolutionary process of plant adaptation to the heterogeneous environment caused by altitude gradients in plateau mountain ecosystems can provide novel insight into species' responses to global changes. Flower color is the most conspicuous and highly diverse trait in nature. Herein, the gene expression patterns, evolutionary adaptation and metabolites changes of different-colored flowers of alpine Rhododendron L. species along altitude gradients were investigated based on a combined analysis of transcriptomics and metabolomics. Differentially expressed genes were found to be related to the biosynthesis of carbohydrates, fatty acids, amino acids and flavonoids, suggesting their important roles in the altitude adaptability of Rhododendron species. The evolution rate of high-altitude species was faster than that of low-altitude species. Genes related to DNA repair, mitogen-activated protein kinase and ABA signal transduction, and lipoic acid and propanoate metabolism were positively selected in the flowers of high-altitude Rhododendron species and those associated with carotenoid biosynthesis pathway, ABA signal transduction and ethylene signal transduction were positively selected in low-altitude species. These results indicated that the genes with differentiated expressions or functions exhibit varying evolution during the adaptive divergence of heterogeneous environment caused by altitude gradients. Flower-color variation might be attributed to the significant differences in gene expression or metabolites related to sucrose, flavonoids and carotenoids at the transcription or metabolism levels of Rhododendron species. This work suggests that Rhododendron species have multiple molecular mechanisms in their adaptation to changing environments caused by altitude gradients.

Understanding Diversity, Evolution, and Structure of Small Heat Shock Proteins in Annelida Through in Silico Analyses

Small heat shock proteins (sHsps) are oligomeric stress proteins characterized by an α-crystallin domain (ACD). These proteins are localized in different subcellular compartments and play critical roles in the stress physiology of tissues, organs, and whole multicellular eukaryotes. They are ubiquitous proteins found in all living organisms, from bacteria to mammals, but they have never been studied in annelids. Here, a data set of 23 species spanning the annelid tree of life, including mostly transcriptomes but also two genomes, was interrogated and 228 novel putative sHsps were identified and manually curated. The analysis revealed very high protein diversity and showed that a significant number of sHsps have a particular dimeric architecture consisting of two tandemly repeated ACDs. The phylogenetic analysis distinguished three main clusters, two of them containing both monomeric sHsps, and ACDs located downstream in the dimeric sHsps, and the other one comprising the upstream ACDs from those dimeric forms. Our results support an evolutionary history of these proteins based on duplication events prior to the Spiralia split. Monomeric sHsps 76) were further divided into five subclusters. Physicochemical properties, subcellular location predictions, and sequence conservation analyses provided insights into the differentiating elements of these putative functional groups. Strikingly, three of those subclusters included sHsps with features typical of metazoans, while the other two presented characteristics resembling non-metazoan proteins. This study provides a solid background for further research on the diversity, evolution, and function in the family of the sHsps. The characterized annelid sHsps are disclosed as essential for improving our understanding of this important family of proteins and their pleotropic functions. The features and the great diversity of annelid sHsps position them as potential powerful molecular biomarkers of environmental stress for acting as prognostic tool in a diverse range of environments.

Genetic divergence and range expansion in a western North Pacific coral

Coral poleward range expansions have recently been observed in response to warming oceans. Range expansion can lead to reduced genetic diversity and increased frequency of deleterious mutations that were rare in core populations, potentially limiting the ability for adaptation and persistence in novel environments. Successful expansions that overcome these founder effects and colonize new habitat have been attributed to multiple introductions from different sources, hybridization with native populations, or rapid adaptive evolution. Here, we investigate population genomic patterns of the reef-building coral Acropora hyacinthus along a latitudinal cline that includes a well-established range expansion front in Japan using 2b-RAD sequencing. A total of 184 coral samples were collected across seven sites spanning from ~24°N to near its northern range front at ~33°N. We uncover the presence of three cryptic lineages of A. hyacinthus, which occupy discrete reefs within this region. Only one lineage is present along the expansion front and we find evidence for its historical occupation of marginal habitats. Within this lineage we also find evidence of bottleneck pressures associated with expansion events including higher clonality, increased linkage disequilibrium, and lower genetic diversity in range edge populations compared to core populations. Asymmetric migration between populations was also detected with lower migration from edge sites. Lastly, we describe genomic signatures of local adaptation potentially attributed to lower winter temperatures experienced at the more recently expanded northern populations. Together these data illuminate the genomic consequences of range expansion in a coral and highlight how adaptation to discrete environments along expansion fronts may facilitate further range expansion in this temperate coral lineage.

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution

Background

With the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility.

Results

We present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities.

Conclusions

This study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00252-0.

Metatranscriptomic Analysis of Bacterial Communities on Laundered Textiles: A Pilot Case Study

Microbially contaminated washing machines and mild laundering conditions facilitate the survival and growth of microorganisms on laundry, promoting undesired side effects such as malodor formation. Clearly, a deeper understanding of the functionality and hygienic relevance of the laundry microbiota necessitates the analysis of the microbial gene expression on textiles after washing, which—to the best of our knowledge—has not been performed before. In this pilot case study, we used single-end RNA sequencing to generate de novo transcriptomes of the bacterial communities remaining on polyester and cotton fabrics washed in a domestic washing machine in mild conditions and subsequently incubated under moist conditions for 72 h. Two common de novo transcriptome assemblers were used. The final assemblies included 22,321 Trinity isoforms and 12,600 Spades isoforms. A large part of these isoforms could be assigned to the SwissProt database, and was further categorized into “molecular function”, “biological process” and “cellular component” using Gene Ontology (GO) terms. In addition, differential gene expression was used to show the difference in the pairwise comparison of the two tissue types. When comparing the assemblies generated with the two assemblers, the annotation results were relatively similar. However, there were clear differences between the de novo assemblies regarding differential gene expression.

Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

Genes with evidence of positive selection as potentially related to coloniality and the evolution of morphological features among the lophophorates and entoprocts

Evolutionary mechanisms that underlie the origins of coloniality among organisms are diverse. Some animal colonies may be comprised strictly of clonal individuals formed from asexual budding or comprised of a chimera of clonal and sexually produced individuals that fuse secondarily. This investigation focuses on select members of the lophophorates and entoprocts whose evolutionary relationships remain enigmatic even in the age of genomics. Using transcriptomic data sets, two coloniality-based hypotheses are tested in a phylogenetic context to find candidate genes showing evidence of positive selection and potentially convergent molecular signatures among solitary species and taxa-forming colonies from aggregate groups or clonal budding. Approximately 22% of the 387 orthogroups tested showed evidence of positive selection in at least one of the three branch-site tests (CODEML, BUSTED, and aBSREL). Only 12 genes could be reliably associated with a developmental function related to traits linked with coloniality, neuroanatomy, or ciliary fields. Genes testing for both positive selection and convergent molecular characters include orthologues of Radial spoke head, Elongation translation initiation factors, SEC13, and Immediate early response gene5. Maximum likelihood analyses included here resulted in tree topologies typical of other phylogenetic investigations based on wider genomic information. Further genomic and experimental evidence will be needed to resolve whether a solitary ancestor with multiciliated cells that formed aggregate groups gave rise to colonial forms in bryozoans (and perhaps the entoprocts) or that the morphological differences exhibited by phoronids and brachiopods represent trait modifications from a colonial ancestor.

Comparing control options for time-series RNA sequencing experiments in nonmodel organisms: An example from grasses

RNA sequencing (RNA-seq) is a widely used approach to investigate gene expression and increasingly is used in time-course studies to characterize transcriptomic changes over time. Two primary options are available as controls in time-course experiments: samples collected at the first sampling time are used as controls (temporal control, TC) and samples collected in parallel at each individual sampling time are used as controls (biological control, BC). While both approaches are used in experimental studies, we know of no analyses performed to date that directly compare effects of control type choices on identifying differentially expressed genes (DEGs) and subsequent functional analysis. In the current study, we compare experimental results using these different control types for time-course RNA-seq drought stress experiments in two wild grass species in the genus Paspalum. Our results showed BC assemblies gave a higher number of loci in both species. The number of DEGs increased with increasing stress and then decreased dramatically at the recovery time point using both control types. Expression levels of the same DEGs were highly correlated between control types in both species, ranging from r = .653 to r = .852. We also observed similar rank orders of shared enriched Gene Ontology term lists using the two different control types. Collectively, our findings suggest similar results in differential gene expression and functional annotation between control types. The ultimate choice of control type will rely on the experimental length and organism type, with labour time and sequencing costs as additional factors to be considered.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora)

BACKGROUND

Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages.

RESULTS

Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses.

CONCLUSIONS

In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

Transcriptomics illuminate the phylogenetic backbone of tiger beetles

Phylogenomics is progressing rapidly, allowing large strides forward into our understanding of the tree of life. In this study, we generated transcriptomes from ethanol-preserved specimens of 13 tiger beetle species (Coleoptera: Cicindelinae) and one Scaritinae outgroup. From these 14 transcriptomes and seven publicly available transcriptomes, we recovered an average of 2538 loci for phylogenetic analysis. We constructed an evolutionary tree of tiger beetles to examine deep-level relationships and examined the extent to which the composition of the dataset, missing data, gene tree inconsistency and codon position saturation impacted phylogenetic accuracy. Ethanol-preserved specimens yielded similar numbers of loci to specimens originally preserved in costly reagents, showcasing more flexibility in transcriptomics than anticipated. The number of loci and gene tree inconsistency had less impact on downstream results than third codon position saturation and missing data. Our results recovered tiger beetles as sister to Carabidae with strong support, confirming their taxonomic status as an independent family within Adephaga. Within tiger beetles, phylogenetic relationships were robust across all nodes. This new phylogenomic backbone represents a useful framework for future endeavours in tiger beetle systematics and serves as a starting point for the development of less costly target capture toolkits to expand the taxonomic breadth of the future tiger beetle tree of life.

Recent progress in human placental transcriptomics

The placenta serves as a metabolic, respiratory, excretory, and endocrine organ that provides appropriate conditions required for adequate fetal development during pregnancy. The development of particular structures and proper functioning of the placenta are under the influence of sophisticated pathways controlled by the expression of substantial genes that are additionally regulated by long non-coding ribonucleic acids (RNAs). Disruptions to adaptive changes in the placental transcriptome as a response to alterations in the feto-maternal environment may be associated with pregnancy complications and compromised fetal outcomes. The aim of the current paper was to present recent findings in transcriptomics of the human placenta. Different approaches in bioinformatic analyses of the RNA-sequencing results were presented. Novel knowledge about the genes and mechanisms that are crucial for the proper development of the placenta is essential for the understanding what stands behind both the normal and complicated pregnancy.

De novo characterization of placental transcriptome in the Eurasian beaver (Castor fiber L.)

Our pioneering data provide the first comprehensive view of placental transcriptome of the beaver during single and multiple gestation. RNA-Seq and a de novo approach allowed global pattern identification of C. fiber placental transcriptome. Non-redundant beaver transcriptome comprised 211,802,336 nt of placental transcripts, grouped into 128,459 contigs and clustered into 83,951 unigenes. An Ensembl database search revealed 14,487, 14,994, 15,004, 15,267 and 15,892 non-redundant homologs for Ictidomys tridecemlineatus, Rattus norvegicus, Mus musculus, Homo sapiens and Castor canadensis, respectively. Due to expression levels, the identified transcripts were divided into two sets: non-redundant and highly expressed (FPKM > 2 in at least three examined samples), analysed simultaneously. Among 17,009 highly expressed transcripts, 12,147 had BLASTx hits. GO annotations (175,882) were found for 4301 transcripts that were assigned to biological process (16,386), cellular component (9149) and molecular function (8338) categories; 666 unigenes were also classified into 122 KEGG pathways. Comprehensive analyses were performed for 411 and 3078 highly expressed transcripts annotated with a list of processes linked to ‘placenta’ (31 GO terms) or ‘embryo’ (324 GO terms), respectively. Among transcripts with entire CDS annotation, 281 (placenta) and 34 (embryo) alternative splicing events were identified. A total of 8499 putative SNVs (~ 6.2 SNV/transcript and 1.7 SNV/1 kb) were predicted with 0.1 minimum frequency and maximum variant quality (p value 10e⁻⁹). Our results provide a broad-based characterization of the global expression pattern of the beaver placental transcriptome. Enhancement of transcriptomic resources for C. fiber should improve understanding of crucial pathways relevant to proper placenta development and successful reproduction.

The Toxins of Nemertean Worms

Most ribbon worms (phylum: Nemertea) are found in marine environments, where they act as predators and scavengers. They are characterized by an eversible proboscis that isused to hunt for prey and thick mucus covering their skin. Both proboscis and epidermal mucus mediate toxicity to predators and preys. Research into the chemical nature of the substances that render toxicity has not been extensive, but it has nevertheless led to the identification of several compounds of potential medicinal use or for application in biotechnology. This review provides a complete account of the current status of research into nemertean toxins.

Delegating Sex: Differential Gene Expression in Stolonizing Syllids Uncovers the Hormonal Control of Reproduction

Stolonization in syllid annelids is a unique mode of reproduction among animals. During the breeding season, a structure resembling the adult but containing only gametes, called stolon, is formed generally at the posterior end of the animal. When stolons mature, they detach from the adult and gametes are released into the water column. The process is synchronized within each species, and it has been reported to be under environmental and endogenous control, probably via endocrine regulation. To further understand reproduction in syllids and to elucidate the molecular toolkit underlying stolonization, we generated Illumina RNA-seq data from different tissues of reproductive and nonreproductive individuals of Syllis magdalena and characterized gene expression during the stolonization process. Several genes involved in gametogenesis (ovochymase, vitellogenin, testis-specific serine/threonine-kinase), immune response (complement receptor 2), neuronal development (tyrosine-protein kinase Src42A), cell proliferation (alpha-1D adrenergic receptor), and steroid metabolism (hydroxysteroid dehydrogenase 2) were found differentially expressed in the different tissues and conditions analyzed. In addition, our findings suggest that several neurohormones, such as methyl farnesoate, dopamine, and serotonin, might trigger stolon formation, the correct maturation of gametes and the detachment of stolons when gametogenesis ends. The process seems to be under circadian control, as indicated by the expression patterns of r-opsins. Overall, our results shed light into the genes that orchestrate the onset of gamete formation and improve our understanding of how some hormones, previously reported to be involved in reproduction and metamorphosis processes in other invertebrates, seem to also regulate reproduction via stolonization.

Identification of differentially expressed placental transcripts during multiple gestations in the Eurasian beaver (Castor fiber L.)

The Eurasian beaver is one of the largest rodents that, despite its high impact on the environment, is a non-model species that lacks a reference genome. Characterising genes critical for pregnancy outcome can serve as a basis for identifying mechanisms underlying effective reproduction, which is required for the success of endangered species conservation programs. In the present study, high-throughput RNA sequencing (RNA-seq) was used to analyse global changes in the Castor fiber subplacenta transcriptome during multiple pregnancy. De novo reconstruction of the C. fiber subplacenta transcriptome was used to identify genes that were differentially expressed in placentas (n=5) from two females (in advanced twin and triple pregnancy). Analyses of the expression values revealed 124 contigs with significantly different expression; of these, 55 genes were identified using MegaBLAST. Within this group of differentially expressed genes (DEGs), 18 were upregulated and 37 were downregulated in twins. Most DEGs were associated with the following gene ontology terms: cellular process, single organism process, response to stimulus, metabolic process and biological regulation. Some genes were also assigned to the developmental process, the reproductive process or reproduction. Among this group, four genes (namely keratin 19 (Krt19) and wingless-type MMTV integration site family - member 2 (Wnt2), which were downregulated in twins, and Nik-related kinase (Nrk) and gap junction protein β2 (Gjb2), which were upregulated in twins) were assigned to placental development and nine (Krt19, Wnt2 and integrin α₇ (Itga7), downregulated in twins, and Nrk, gap junction protein β6 (Gjb6), GATA binding protein 6 (Gata6), apolipoprotein A-I (ApoA1), apolipoprotein B (ApoB) and haemoglobin subunit α₁ (HbA1), upregulated in twins) were assigned to embryo development. The results of the present study indicate that the number of fetuses affects the expression profile in the C. fiber subplacental transcriptome. Enhancement of transcriptomic resources for C. fiber will improve understanding of the pathways relevant to proper placental development and successful reproduction.

Polymorphic adaptations in metazoans to establish and maintain photosymbioses

Mutualistic symbioses are common throughout the animal kingdom. Rather unusual is a form of symbiosis, photosymbiosis, where animals are symbiotic with photoautotrophic organisms. Photosymbiosis is found among sponges, cnidarians, flatworms, molluscs, ascidians and even some amphibians. Generally the animal host harbours a phototrophic partner, usually a cyanobacteria or a unicellular alga. An exception to this rule is found in some sea slugs, which only retain the chloroplasts of the algal food source and maintain them photosynthetically active in their own cytosol - a phenomenon called 'functional kleptoplasty'. Research has focused largely on the biodiversity of photosymbiotic species across a range of taxa. However, many questions with regard to the evolution of the ability to establish and maintain a photosymbiosis are still unanswered. To date, attempts to understand genome adaptations which could potentially lead to the evolution of photosymbioses have only been performed in cnidarians. This knowledge gap for other systems is mainly due to a lack of genetic information, both for non-symbiotic and symbiotic species. Considering non-photosymbiotic species is, however, important to understand the factors that make symbiotic species so unique. Herein we provide an overview of the diversity of photosymbioses across the animal kingdom and discuss potential scenarios for the evolution of this association in different lineages. We stress that the evolution of photosymbiosis is probably based on genome adaptations, which (i) lead to recognition of the symbiont to establish the symbiosis, and (ii) are needed to maintain the symbiosis. We hope to stimulate research involving sequencing the genomes of various key taxa to increase the genomic resources needed to understand the most fundamental question: how have animals evolved the ability to establish and maintain a photosymbiosis?

Transcriptomic Analysis of Pseudoscorpion Venom Reveals a Unique Cocktail Dominated by Enzymes and Protease Inhibitors

Transcriptomic and genomic analyses have illuminated the diversity of venoms in three of the four venomous arachnid orders (scorpions, spiders, and ticks). To date, no venom gland transcriptome analysis has been available for pseudoscorpions, the fourth venomous arachnid lineage. To redress this gap, we sequenced an mRNA library generated from the venom glands of the species Synsphyronus apimelus (Garypidae). High-throughput sequencing by the Illumina protocol, followed by de novo assembly, resulted in a total of 238,331 transcripts. From those, we annotated 131 transcripts, which code for putative peptides/proteins with similar sequences to previously reported venom components available from different arachnid species in protein databases. Transcripts putatively coding for enzymes showed the richest diversity, followed by other venom components such as peptidase inhibitors, cysteine-rich peptides, and thyroglobulin 1-like peptides. Only 11 transcripts were found that code for putatively low molecular mass spider toxins. This study constitutes the first report of the diversity of components within pseudoscorpion venom.

A phylotranscriptomic backbone of the orb-weaving spider family Araneidae (Arachnida, Araneae) supported by multiple methodological approaches

The orb-weaving spider family Araneidae is extremely diverse (>3100 spp.) and its members can be charismatic terrestrial arthropods, many of them recognizable by their iconic orbicular snare web, such as the common garden spiders. Despite considerable effort to better understand their backbone relationships based on multiple sources of data (morphological, behavioral and molecular), pervasive low support remains in recent studies. In addition, no overarching phylogeny of araneids is available to date, hampering further comparative work. In this study, we analyze the transcriptomes of 33 taxa, including 19 araneids - 12 of them new to this study - representing most of the core family lineages, to examine the relationships within the family using genomic-scale datasets resulting from various methodological treatments, namely ortholog selection and gene occupancy as a measure of matrix completion. Six matrices were constructed to assess these effects by varying orthology inference method and gene occupancy threshold. Orthology methods used are the benchmarking tool BUSCO and the tree-based method UPhO; three gene occupancy thresholds (45%, 65%, 85%) were used to assess the effect of missing data. Gene tree and species tree-based methods (including multi-species coalescent and concatenation approaches, as well as maximum likelihood and Bayesian inference) were used totalling 17 analytical treatments. The monophyly of Araneidae and the placement of core araneid lineages were supported, together with some previously unsound backbone divergences; these include high support for Zygiellinae as the earliest diverging subfamily (followed by Nephilinae), the placement of Gasteracanthinae as sister group to Cyclosa and close relatives, and close relationships between the Araneus + Neoscona clade and Cyrtophorinae + Argiopinae clade. Incongruences were relegated to short branches in the clade comprising Cyclosa and its close relatives. We found congruence between most of the completed analyses, with minimal topological effects from occupancy/missing data and orthology assessment. The resulting number of genes by certain combinations of orthology and occupancy thresholds being analyzed had the greatest effect on the resulting trees, with anomalous outcomes recovered from analysis of lower numbers of genes.

De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species

Background

Salmonid fishes exhibit high levels of phenotypic and ecological variation and are thus ideal model systems for studying evolutionary processes of adaptive divergence and speciation. Furthermore, salmonids are of major interest in fisheries, aquaculture, and conservation research. Improving understanding of the genetic mechanisms underlying traits in these species would significantly progress research in these fields. Here we generate high quality de novo transcriptomes for four salmonid species: Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Arctic charr (Salvelinus alpinus), and European whitefish (Coregonus lavaretus). All species except Atlantic salmon have no reference genome publicly available and few if any genomic studies to date.

Results

We used paired-end RNA-seq on Illumina to generate high coverage sequencing of multiple individuals, yielding between 180 and 210 M reads per species. After initial assembly, strict filtering was used to remove duplicated, redundant, and low confidence transcripts. The final assemblies consisted of 36,505 protein-coding transcripts for Atlantic salmon, 35,736 for brown trout, 33,126 for Arctic charr, and 33,697 for European whitefish and are made publicly available. Assembly completeness was assessed using three approaches, all of which supported high quality of the assemblies: 1) ~78% of Actinopterygian single-copy orthologs were successfully captured in our assemblies, 2) orthogroup inference identified high overlap in the protein sequences present across all four species (40% shared across all four and 84% shared by at least two), and 3) comparison with the published Atlantic salmon genome suggests that our assemblies represent well covered (~98%) protein-coding transcriptomes. Thorough comparison of the generated assemblies found that 84-90% of transcripts in each assembly were orthologous with at least one of the other three species. We also identified 34-37% of transcripts in each assembly as paralogs. We further compare completeness and annotation statistics of our new assemblies to available related species.

Conclusion

New, high-confidence protein-coding transcriptomes were generated for four ecologically and economically important species of salmonids. This offers a high quality pipeline for such complex genomes, represents a valuable contribution to the existing genomic resources for these species and provides robust tools for future investigation of gene expression and sequence evolution in these and other salmonid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4379-x) contains supplementary material, which is available to authorized users.

The ‘Peripatos' in Eurogondwana? — Lack of evidence that southeast Asian onychophorans walked through Europe

Onychophorans, or velvet worms, are cryptic but extremely charismatic terrestrial invertebrates that have often been the subject of interesting biogeographic debate. Despite great interest, a well resolved and complete phylogeny of the group and a reliable chronogram have been elusive due to their broad geographic distribution, paucity of samples, and challenging molecular composition. Here we present a molecular phylogenetic analysis of Onychophora that includes previously unsampled and undersampled lineages and we analyse the expanded dataset using a series of nested taxon sets designed to increase the amount of information available for particular subclades. These include a dataset with outgroups, one restricted to the ingroup taxa, and three others for Peripatopsidae, Peripatidae and Neopatida (= the Neotropical Peripatidae). To explore competing biogeographic scenarios we generate a new time tree for Onychophora using the few available reliable fossils as calibration points. Comparing our results to those of Cyphophthalmi, we reconsider the hypothesis that velvet worms reached Southeast Asia via Eurogondwana, and conclude that a more likely scenario is that they reached Southeast Asia by rafting on the Sibumasu terrane. Our phylogenetic results support the reciprocal monophyly of both families as well as an early division between East and West Gondwana, also in both families, each beginning to diversify between the Permian and the Jurassic. Peripatopsidae clearly supports paraphyly of South Africa with respect to southern South America (Chile) and a sister group relationship of the Southeast Asian/New Guinean Paraperipatus to the Australian/New Zealand taxa. The latter includes a clade that divides between Western Australia and Eastern Australia and two sister clades of trans-Tasman species (one oviparous and one viviparous). This pattern clearly shows that oviparity is secondarily derived in velvet worms. Peripatidae finds a sister group relationship between the Southeast Asian Eoperipatus and the West Gondwanan clade, which divides into the African Mesoperipatus and Neopatida. The latter shows a well supported split between the Pacific Oroperipatus (although it is unclear whether they form one or two clades) and a sister clade that includes the members of the genera Peripatus, Epiperipatus, Macroperipatus and representatives of the monotypic genera Cerradopatus, Plicatoperipatus and Principapillatus. However, Peripatus, Epiperipatus and Macroperipatus are not monophyletic, and all the species from the monotypic genera are related to geographically close species. The same goes for the type species of Macroperipatus (from Trinidad, and sister group to other Trinidad and Tobago species of Epiperipatus) and Epiperipatus (from French Guiana, and related to other Guyana shield species of Epiperipatus and Peripatus). Geographic structure within Neopatida is largely obscured by an unresolved backbone, but many well supported instances of generic non-monophyly challenge the current taxonomic framework, which has often relied on anatomical characters that are untested phylogenetically.

Identification of differentially expressed genes in flower, leaf and bulb scale of Lilium oriental hybrid 'Sorbonne' and putative control network for scent genes

BACKGROUND

Lily is an economically important plant, with leaves and bulbs consisting of overlapping scales, large ornamental flowers and a very large genome. Although it is recognized that flowers and bulb scales are modified leaves, very little is known about the genetic control and biochemical differentiation underlying lily organogenesis and development. Here we examined the differentially expressed genes in flower, leaf and scale of lily, using RNA-sequencing, and identified organ-specific genes, including transcription factors, genes involved in photosynthesis in leaves, carbohydrate metabolism in bulb scales and scent and color production in flowers.

RESULTS

Over 11Gb data were obtained and 2685, 2296, and 1709 differentially expressed genes were identified in the three organs, with 581, 662 and 977 unique DEGs in F-vs-S, L-vs-S and L-vs-F comparisons. By functional enrichment analysis, genes likely to be involved in biosynthetic pathways leading to floral scent production, such as 1-deoxy-D-xylulose-5-phosphate synthase (DXS), 3-ketoacyl-CoA thiolase (KAT), hydroperoxide lyase (HPL), geranylgeranyl pyrophosphate (GGPP) 4-hydroxy-3-methylbut-2-en-1-yl diphosphate (HDS) and terpene synthase (TPS), and floral color genes, such as dihydroflavonol 4-reductase (DFR), chalcone synthase (CHS), chalcone isomerase (CHI), flavonol synthase (FLS) were identified. Distinct groups of genes that participate in starch and sucrose metabolism, such as sucrose synthase (SS), invertase (INV), sucrose phosphate synthase (SPS), starch synthase (SSS), starch branching enzyme (SBE), ADP-glucose pyrophosphorylase (AGP) andβ-amylase (BAM) and photosynthesis genes (Psa, Psb, Pet and ATP) were also identified. The expression of six floral fragrance-related DGEs showed agreement between qRT-PCR results and RPKM values, confirming the value of the data obtained by RNA-seq. We obtained the open reading frame of the terpene synthase gene from Lilium 'Sorbonne', designated LsTPS, which had 99.55% homology to transcript CL4520.Contig5_All. In addition, 54, 48 and 50 differently expressed transcription factor were identified by pairwise comparisons between the three organs and a regulatory network for monoterpene biosynthesis was constructed.

CONCLUSIONS

Analysis of differentially expressed genes in flower, leaf and bulb scale of lily, using second generation sequencing technology, yielded detailed information on lily metabolic differentiation in three organs. Analysis of the expression of flower scent biosynthesis genes has provided a model for the regulation of the pathway and identified a candidate gene encoding an enzyme catalyzing the final step in scent production. These digital gene expression profiles provide a valuable and informative database for the further identification and analysis of structural genes and transcription factors in different lily organs and elucidation of their function.

Development and Characterization of 12 Novel Polymorphic Microsatellite Loci for the Mammal Chewing Louse Geomydoecus aurei (Insecta: Phthiraptera) and a Comparison of Next-Generation Sequencing Approaches for Use in Parasitology

Next-generation sequencing methodologies open the door for evolutionary studies of wildlife parasites. We used 2 next-generation sequencing approaches to discover microsatellite loci in the pocket gopher chewing louse Geomydoecus aurei for use in population genetic studies. In one approach, we sequenced a library enriched for microsatellite loci; in the other approach, we mined microsatellites from genomic sequences. Following microsatellite discovery, promising loci were tested for amplification and polymorphism in 390 louse individuals from 13 pocket gopher hosts. In total, 12 loci were selected for analysis (6 from each methodology), none of which exhibited evidence of null alleles or heterozygote deficiencies. These 12 loci showed adequate genetic diversity for population-level analyses, with 3-9 alleles per locus with an average H_E per locus ranging from 0.32 to 0.70. Analysis of Molecular Variance (AMOVA) indicated that genetic variation among infrapopulations accounts for a low, but significant, percentage of the overall genetic variation, and individual louse infrapopulations showed F_ST values that were significantly different from zero in the majority of pairwise infrapopulation comparisons, despite all 13 infrapopulations being taken from the same locality. Therefore, these 12 polymorphic markers will be useful at the infrapopulation and population levels for future studies involving G. aurei. This study shows that next-generation sequencing methodologies can successfully be used to efficiently obtain data for a variety of evolutionary questions.

De novo transcriptome assembly and developmental mode specific gene expression of Pygospio elegans

Species with multiple different larval developmental modes are interesting models for the study of mechanisms underlying developmental mode transitions and life history evolution. Pygospio elegans, a small, tube-dwelling polychaete worm commonly found in estuarine and marine habitats around the northern hemisphere, is one species with variable developmental modes. To provide new genomic resources for studying P. elegans and to address the differences in gene expression between individuals producing offspring with different larval developmental modes, we performed whole transcriptome Illumina RNA sequencing of adult worms from two populations and prepared a de novo assembly of the P. elegans transcriptome. The transcriptome comprises 66,233 unigenes, of which 33,807 contain predicted coding sequences, 26,448 have at least one functional annotation, and 3,076 are classified as putative long non-coding RNAs. We found more than 8,000 unigenes significantly differentially expressed between adult worms from populations producing either planktonic or benthic larvae. This comprehensive transcriptome resource for P. elegans adds to the available genomic data for annelids and can be used to uncover mechanisms allowing developmental variation in this and potentially other marine invertebrate species.

A genome-wide transcriptome map of pistachio (Pistacia vera L.) provides novel insights into salinity-related genes and marker discovery

Background

Pistachio (Pistacia vera L.) is one of the most important commercial nut crops worldwide. It is a salt-tolerant and long-lived tree, with the largest cultivation area in Iran. Climate change and subsequent increased soil salt content have adversely affected the pistachio yield in recent years. However, the lack of genomic/global transcriptomic sequences on P. vera impedes comprehensive researches at the molecular level. Hence, whole transcriptome sequencing is required to gain insight into functional genes and pathways in response to salt stress.

Results

RNA sequencing of a pooled sample representing 24 different tissues of two pistachio cultivars with contrasting salinity tolerance under control and salt treatment by Illumina Hiseq 2000 platform resulted in 368,953,262 clean 100 bp paired-ends reads (90 Gb). Following creating several assemblies and assessing their quality from multiple perspectives, we found that using the annotation-based metrics together with the length-based parameters allows an improved assessment of the transcriptome assembly quality, compared to the solely use of the length-based parameters. The generated assembly by Trinity was adopted for functional annotation and subsequent analyses. In total, 29,119 contigs annotated against all of five public databases, including NR, UniProt, TAIR10, KOG and InterProScan. Among 279 KEGG pathways supported by our assembly, we further examined the pathways involved in the plant hormone biosynthesis and signaling as well as those to be contributed to secondary metabolite biosynthesis due to their importance under salinity stress. In total, 11,337 SSRs were also identified, which the most abundant being dinucleotide repeats. Besides, 13,097 transcripts as candidate stress-responsive genes were identified. Expression of some of these genes experimentally validated through quantitative real-time PCR (qRT-PCR) that further confirmed the accuracy of the assembly. From this analysis, the contrasting expression pattern of NCED3 and SOS1 genes were observed between salt-sensitive and salt-tolerant cultivars.

Conclusion

This study, as the first report on the whole transcriptome survey of P. vera, provides important resources and paves the way for functional and comparative genomic studies on this major tree to discover the salinity tolerance-related markers and stress response mechanisms for breeding of new pistachio cultivars with more salinity tolerance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3989-7) contains supplementary material, which is available to authorized users.

Origin and Evolution of the Sponge Aggregation Factor Gene Family

Although discriminating self from nonself is a cardinal animal trait, metazoan allorecognition genes do not appear to be homologous. Here, we characterize the Aggregation Factor (AF) gene family, which encodes putative allorecognition factors in the demosponge Amphimedon queenslandica, and trace its evolution across 24 sponge (Porifera) species. The AF locus in Amphimedon is comprised of a cluster of five similar genes that encode Calx-beta and Von Willebrand domains and a newly defined Wreath domain, and are highly polymorphic. Further AF variance appears to be generated through individualistic patterns of RNA editing. The AF gene family varies between poriferans, with protein sequences and domains diagnostic of the AF family being present in Amphimedon and other demosponges, but absent from other sponge classes. Within the demosponges, AFs vary widely with no two species having the same AF repertoire or domain organization. The evolution of AFs suggests that their diversification occurs via high allelism, and the continual and rapid gain, loss and shuffling of domains over evolutionary time. Given the marked differences in metazoan allorecognition genes, we propose the rapid evolution of AFs in sponges provides a model for understanding the extensive diversification of self-nonself recognition systems in the animal kingdom.

Differential transcriptome analysis supports Rhodnius montenegrensis and Rhodnius robustus (Hemiptera, Reduviidae, Triatominae) as distinct species

Chagas disease is one of the main parasitic diseases found in Latin America and it is estimated that between six and seven million people are infected worldwide. Its etiologic agent, the protozoan Trypanosoma cruzi, is transmitted by triatomines, some of which from the genus Rhodnius. Twenty species are currently recognized in this genus, including some closely related species with low levels of morphological differentiation, such as Rhodnius montenegrensis and Rhodnius robustus. In order to investigate genetic differences between these two species, we generated large-scale RNA-sequencing data (consisting of four RNA-seq libraries) from the heads and salivary glands of males of R. montenegrensis and R. robustus. Transcriptome assemblies produced for each species resulted in 64,952 contigs for R. montenegrensis and 70,894 contigs for R. robustus, with N50 of approximately 2,100 for both species. SNP calling based on the more complete R. robustus assembly revealed 3,055 fixed interspecific differences and 216 transcripts with high levels of divergence which contained only fixed differences between the two species. A gene ontology enrichment analysis revealed that these highly differentiated transcripts were enriched for eight GO terms related to AP-2 adaptor complex, as well as other interesting genes that could be involved in their differentiation. The results show that R. montenegrensis and R. robustus have a substantial quantity of fixed interspecific polymorphisms, which suggests a high degree of genetic divergence between the two species and likely corroborates the species status of R. montenegrensis.

Comparative transcriptomics of Entelegyne spiders (Araneae, Entelegynae), with emphasis on molecular evolution of orphan genes

Next-generation sequencing technology is rapidly transforming the landscape of evolutionary biology, and has become a cost-effective and efficient means of collecting exome information for non-model organisms. Due to their taxonomic diversity, production of interesting venom and silk proteins, and the relative scarcity of existing genomic resources, spiders in particular are excellent targets for next-generation sequencing (NGS) methods. In this study, the transcriptomes of six entelegyne spider species from three genera (Cicurina travisae, C. vibora, Habronattus signatus, H. ustulatus, Nesticus bishopi, and N. cooperi) were sequenced and de novo assembled. Each assembly was assessed for quality and completeness and functionally annotated using gene ontology information. Approximately 100 transcripts with evidence of homology to venom proteins were discovered. After identifying more than 3,000 putatively orthologous genes across all six taxa, we used comparative analyses to identify 24 instances of positively selected genes. In addition, between ~ 550 and 1,100 unique orphan genes were found in each genus. These unique, uncharacterized genes exhibited elevated rates of amino acid substitution, potentially consistent with lineage-specific adaptive evolution. The data generated for this study represent a valuable resource for future phylogenetic and molecular evolutionary research, and our results provide new insight into the forces driving genome evolution in taxa that span the root of entelegyne spider phylogeny.

High-quality RNA extraction from the sea urchin Paracentrotus lividus embryos

The sea urchin Paracentrotus lividus (Lamarck, 1816) is a keystone herbivore in the Mediterranean Sea due to its ability to transform macroalgal-dominated communities into barren areas characterized by increased cover of bare substrates and encrusting coralline algae, reduced biodiversity and altered ecosystem functions. P. lividus is also an excellent animal model for toxicology, physiology and biology investigations having been used for more than a century as a model for embryological studies with synchronously developing embryos which are easy to manipulate and analyze for morphological aberrations. Despite its importance for the scientific community, the complete genome is still not fully annotated. To date, only a few molecular tools are available and a few Next Generation Sequencing (NGS) studies have been performed. Here we aimed at setting-up an RNA extraction method to obtain high quality and sufficient quantity of RNA for NGS from P. lividus embryos at the pluteus stage. We compared five different RNA extraction protocols from four different pools of plutei (500, 1000, 2500 and 5000 embryos): TRIzol®, and four widely-used Silica Membrane kits, GenElute™ Mammalian Total RNA Miniprep Kit, RNAqueous® Micro Kit, RNeasy® Micro Kit and Aurum™ Total RNA Mini Kit. The quantity of RNA isolated was evaluated using NanoDrop. The quality, considering the purity, was measured as A260/A280 and A260/230 ratios. The integrity was measured by RNA Integrity Number (RIN). Our results demonstrated that the most efficient procedures were GenElute, RNeasy and Aurum, producing a sufficient quantity of RNA for NGS. The Bioanalyzer profiles and RIN values revealed that the most efficient methods guaranteeing for RNA integrity were RNeasy and Aurum combined with an initial preservation in RNAlater. This research represents the first attempt to standardize a method for high-quality RNA extraction from sea urchin embryos at the pluteus stage, providing a new resource for this established model marine organism.

Implications for management and conservation of the population genetic structure of the wedge clam Donax trunculus across two biogeographic boundaries

In a resource management perspective, the understanding of the relative influence of the physical factors on species connectivity remains a major challenge and is also of great ecological and conservation biology interest. Despite the overfishing threat on the wedge clam Donax trunculus in Europe, relatively little information is known about its population genetic structure and connectivity and their consequences on conservation policies. We employed 16 microsatellite loci to characterise the genetic diversity and population structure of D. trunculus. A total of 514 samples from seven different localities along the Atlantic-Mediterranean transition, from the Atlantic (Gulf of Cádiz) to the north-western Mediterranean were genotyped. The analysis of the population genetic structure displayed a clear distinction along the Atlantic-Mediterranean transition with different clusters in the Atlantic Ocean, the Alboran Sea and the northwestern Mediterranean. Consequently, we recommend that these three areas should be considered as different management units. We showed that all populations seem to be at high long-term risk of extinction with the exception of the protected Doñana National Park population which still seems to have evolutionary potential. Therefore, our results emphasized the necessity of protection of this economic resource and the validity of molecular tools to evaluate the population dynamics.

Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks

BACKGROUND

Mollusks display a striking morphological disparity, including, among others, worm-like animals (the aplacophorans), snails and slugs, bivalves, and cephalopods. This phenotypic diversity renders them ideal for studies into animal evolution. Despite being one of the most species-rich phyla, molecular and in silico studies concerning specific key developmental gene families are still scarce, thus hampering deeper insights into the molecular machinery that governs the development and evolution of the various molluscan class-level taxa.

RESULTS

Next-generation sequencing was used to retrieve transcriptomes of representatives of seven out of the eight recent class-level taxa of mollusks. Similarity searches, phylogenetic inferences, and a detailed manual curation were used to identify and confirm the orthology of numerous molluscan Hox and ParaHox genes, which resulted in a comprehensive catalog that highlights the evolution of these genes in Mollusca and other metazoans. The identification of a specific molluscan motif in the Hox paralog group 5 and a lophotrochozoan ParaHox motif in the Gsx gene is described. Functional analyses using KEGG and GO tools enabled a detailed description of key developmental genes expressed in important pathways such as Hedgehog, Wnt, and Notch during development of the respective species. The KEGG analysis revealed Wnt8, Wnt11, and Wnt16 as Wnt genes hitherto not reported for mollusks, thereby enlarging the known Wnt complement of the phylum. In addition, novel Hedgehog (Hh)-related genes were identified in the gastropod Lottia cf. kogamogai, demonstrating a more complex gene content in this species than in other mollusks.

CONCLUSIONS

The use of de novo transcriptome assembly and well-designed in silico protocols proved to be a robust approach for surveying and mining large sequence data in a wide range of non-model mollusks. The data presented herein constitute only a small fraction of the information retrieved from the analysed molluscan transcriptomes, which can be promptly employed in the identification of novel genes and gene families, phylogenetic inferences, and other studies using molecular tools. As such, our study provides an important framework for understanding some of the underlying molecular mechanisms involved in molluscan body plan diversification and hints towards functions of key developmental genes in molluscan morphogenesis.

Elucidating and mining the Tulipa and Lilium transcriptomes

Genome sequencing remains a challenge for species with large and complex genomes containing extensive repetitive sequences, of which the bulbous and monocotyledonous plants tulip and lily are examples. In such a case, sequencing of only the active part of the genome, represented by the transcriptome, is a good alternative to obtain information about gene content. In this study we aimed to generate a high quality transcriptome of tulip and lily and to make this data available as an open-access resource via a user-friendly web-based interface. The Illumina HiSeq 2000 platform was applied and the transcribed RNA was sequenced from a collection of different lily and tulip tissues, respectively. In order to obtain good transcriptome coverage and to facilitate effective data mining, assembly was done using different filtering parameters for clearing out contamination and noise of the RNAseq datasets. This analysis revealed limitations of commonly applied methods and parameter settings used in de novo transcriptome assembly. The final created transcriptomes are publicly available via a user friendly Transcriptome browser ( http://www.bioinformatics.nl/bulbs/db/species/index ). The usefulness of this resource has been exemplified by a search for all potential transcription factors in lily and tulip, with special focus on the TCP transcription factor family. This analysis and other quality parameters point out the quality of the transcriptomes, which can serve as a basis for further genomics studies in lily, tulip, and bulbous plants in general.

Phylogenomics and Divergence Dating of Fungus-Farming Ants (Hymenoptera: Formicidae) of the Genera Sericomyrmex and Apterostigma

Fungus-farming ("attine") ants are model systems for studies of symbiosis, coevolution, and advanced eusociality. A New World clade of nearly 300 species in 15 genera, all attine ants cultivate fungal symbionts for food. In order to better understand the evolution of ant agriculture, we sequenced, assembled, and analyzed transcriptomes of four different attine ant species in two genera: three species in the higher-attine genus Sericomyrmex and a single lower-attine ant species, Apterostigma megacephala, representing the first genomic data for either genus. These data were combined with published genomes of nine other ant species and the honey bee Apis mellifera for phylogenomic and divergence-dating analyses. The resulting phylogeny confirms relationships inferred in previous studies of fungus-farming ants. Divergence-dating analyses recovered slightly older dates than most prior analyses, estimating that attine ants originated 53.6-66.7 million of years ago, and recovered a very long branch subtending a very recent, rapid radiation of the genus Sericomyrmex. This result is further confirmed by a separate analysis of the three Sericomyrmex species, which reveals that 92.71% of orthologs have 99% - 100% pairwise-identical nucleotide sequences. We searched the transcriptomes for genes of interest, most importantly argininosuccinate synthase and argininosuccinate lyase, which are functional in other ants but which are known to have been lost in seven previously studied attine ant species. Loss of the ability to produce the amino acid arginine has been hypothesized to contribute to the obligate dependence of attine ants upon their cultivated fungi, but the point in fungus-farming ant evolution at which these losses occurred has remained unknown. We did not find these genes in any of the sequenced transcriptomes. Although expected for Sericomyrmex species, the absence of arginine anabolic genes in the lower-attine ant Apterostigma megacephala strongly suggests that the loss coincided with the origin of attine ants.

Testis-specific RNA-Seq of Calligrapha (Chrysomelidae) as a transcriptomic resource for male-biased gene inquiry in Coleoptera

We report the architecture of testis transcriptomes of four closely related species of Calligrapha (Chrysomelidae) beetles, which diverged during the last 3 million years. Five cDNA libraries were sequenced using Illumina HiSeq technology, retrieving 102 884-176 514 assembled contigs, of which ~33-45% of these longer than 499 nt were functionally annotated. Annotation and sequence similarity comparisons of these libraries revealed high homogeneity in gene composition and the presence of several functional candidates related to reproduction or reproductive processes (0.72-1.08% of annotated sequences). Stringent sequence similarity analyses of these transcriptomes against empirically demonstrated male-biased genes in Drosophila melanogaster and Tribolium castaneum allowed the identification of 77 homologues in Calligrapha, possible candidates of male-biased expression. Some of these genes - including CG9313, Tektin-A or tomboy40 - were confirmed as orthologs of these male-biased genes using phylogenetic inference and available model insect data, increasing our confidence that they represent functional homologues too. Our transcriptomes are a valuable transcriptomic resource for the analysis of male-biased genes in Calligrapha, which has the added interest of including several female-only species. But it simultaneously represents a landmark for similar studies in Coleoptera, broadening the taxonomic diversity currently represented by the model species T. castaneum, and incipient genomic data in other herbivorous lineages, including weevils, longhorn beetles and leaf beetles.

Transcriptome analysis of the copepod Eurytemora affinis upon exposure to endocrine disruptor pesticides: Focus on reproduction and development

Copepods-which include freshwater and marine species-represent the most abundant group of aquatic invertebrates. Among them, the calanoid copepod Eurytemora affinis is widely represented in the northern hemisphere estuaries and has become a species of interest in ecotoxicology. Like other non-target organisms, E. affinis may be exposed to a wide range of chemicals such as endocrine disruptors (EDs). This study investigated the gene expression variation in E. affinis after exposure to ED pesticides-chosen as model EDs-in order to (i) improve the knowledge on their effects in crustaceans, and (ii) highlight relevant transcripts for further development of potential biomarkers of ED exposure/effect. The study focused on the reproduction function in response to ED. Copepods were exposed to sublethal concentrations of pyriproxyfen (PXF) and chlordecone (CLD) separately. After 48h, males and females (400 individuals each) were sorted for RNA extraction. Their transcriptome was pyrosequenced using the Illumina(®) technology. Contigs were blasted and functionally annotated using Blast2GO(®). The differential expression analysis between ED- and acetone-exposed organisms was performed according to sexes and contaminants. Half of the 19,721 contigs provided by pyrosequencing were annotated, mostly (80%) from arthropod sequences. Overall, 2,566 different genes were differentially expressed after ED exposures in comparison with controls. As many genes were differentially expressed after PXF exposure as after CLD exposure. In contrast, more genes were differentially expressed in males than in females after both exposures. Ninety-seven genes overlapped in all conditions. Finally, 31 transcripts involved in reproduction, growth and development, and changed in both chemical exposures were selected as potential candidates for future development of biomarkers.

Tentacle Transcriptome and Venom Proteome of the Pacific Sea Nettle, Chrysaora fuscescens (Cnidaria: Scyphozoa)

Jellyfish venoms are rich sources of toxins designed to capture prey or deter predators, but they can also elicit harmful effects in humans. In this study, an integrated transcriptomic and proteomic approach was used to identify putative toxins and their potential role in the venom of the scyphozoan jellyfish Chrysaora fuscescens. A de novo tentacle transcriptome, containing more than 23,000 contigs, was constructed and used in proteomic analysis of C. fuscescens venom to identify potential toxins. From a total of 163 proteins identified in the venom proteome, 27 were classified as putative toxins and grouped into six protein families: proteinases, venom allergens, C-type lectins, pore-forming toxins, glycoside hydrolases and enzyme inhibitors. Other putative toxins identified in the transcriptome, but not the proteome, included additional proteinases as well as lipases and deoxyribonucleases. Sequence analysis also revealed the presence of ShKT domains in two putative venom proteins from the proteome and an additional 15 from the transcriptome, suggesting potential ion channel blockade or modulatory activities. Comparison of these potential toxins to those from other cnidarians provided insight into their possible roles in C. fuscescens venom and an overview of the diversity of potential toxin families in cnidarian venoms.