Compacting and correcting Trinity and Oases RNA-Seq de novo assemblies

Phylogenomics reveals the evolutionary origins of lichenization in chlorophyte algae

Mutualistic symbioses have contributed to major transitions in the evolution of life. Here, we investigate the evolutionary history and the molecular innovations at the origin of lichens, which are a symbiosis established between fungi and green algae or cyanobacteria. We de novo sequence the genomes or transcriptomes of 12 lichen algal symbiont (LAS) and closely related non-symbiotic algae (NSA) to improve the genomic coverage of Chlorophyte algae. We then perform ancestral state reconstruction and comparative phylogenomics. We identify at least three independent gains of the ability to engage in the lichen symbiosis, one in Trebouxiophyceae and two in Ulvophyceae, confirming the convergent evolution of the lichen symbioses. A carbohydrate-active enzyme from the glycoside hydrolase 8 (GH8) family was identified as a top candidate for the molecular-mechanism underlying lichen symbiosis in Trebouxiophyceae. This GH8 was acquired in lichenizing Trebouxiophyceae by horizontal gene transfer, concomitantly with the ability to associate with lichens fungal symbionts (LFS) and is able to degrade polysaccharides found in the cell wall of LFS. These findings indicate that a combination of gene family expansion and horizontal gene transfer provided the basis for lichenization to evolve in chlorophyte algae.

Macrophages modulate fibrosis during newt lens regeneration

BACKGROUND

Previous studies have suggested that macrophages are present during lens regeneration in newts, but their role in the process is yet to be elucidated.

METHODS

Here we generated a transgenic reporter line using the newt, Pleurodeles waltl, that traces macrophages during lens regeneration. Furthermore, we assessed early changes in gene expression during lens regeneration using two newt species, Notophthalmus viridescens and Pleurodeles waltl. Finally, we used clodronate liposomes to deplete macrophages during lens regeneration in both species and tested the effect of a subsequent secondary injury after macrophage recovery.

RESULTS

Macrophage depletion abrogated lens regeneration, induced the formation of scar-like tissue, led to inflammation, decreased iris pigment epithelial cell (iPEC) proliferation, and increased rates of apoptosis in the eye. Some of these phenotypes persisted throughout the last observation period of 100 days and could be attenuated by exogenous FGF2 administration. A distinct transcript profile encoding acute inflammatory effectors was established for the dorsal iris. Reinjury of the newt eye alleviated the effects of macrophage depletion, including the resolution of scar-like tissue, and re-initiated the regeneration process.

CONCLUSIONS

Together, our findings highlight the importance of macrophages for facilitating a pro-regenerative environment in the newt eye by regulating fibrotic responses, modulating the overall inflammatory landscape, and maintaining the proper balance of early proliferation and late apoptosis of the iPECs.

High fatty acid accumulation and coloration molecular mechanism of the elm mushroom (Pleurotus citrinopileatus)

Pleurotus citrinopileatus is a low-cholesterol, protein-rich, and high-nutrient food. The molecular mechanisms of the compounds and coloration have not been reported. Metabolome and transcriptome were used to clarify the molecular mechanisms of key compounds biosynthesis. K-means analysis identified 19 compounds in P. citrinopileatus, mainly lipids and alkaloids in class 8. In addition, 84 lipids were higher and that the different compounds were mainly enriched in linoleic acid metabolism. A total of 14 compounds detected in the linoleic acid metabolism pathway were significantly up-regulated, while 3 sterol regulatory element binding protein (SREBP) transcription factors were screened. Tryptophan metabolism and riboflavin biosynthesis pathway analysis indicated that 3 Unigenes had tryptophan decarboxylase similar elements, which belonged to tyrosine decarboxylase 1. Moreover, CL15618.Contig5_All had high homology with MFS. In conclusion, the expression of 3 SREBP, the synthesis of isobavachalcone D, and the regulation of riboflavin transport by MCH5 were the reasons for fatty acid accumulation and yellow cap formation in the P. citrinopileatus.

Genomic Resources and Annotations for a Colonial Ascidian, the Light-Bulb Sea Squirt Clavelina lepadiformis

Ascidian embryos have been studied since the birth of experimental embryology at the end of the 19th century. They represent textbook examples of mosaic development characterized by a fast development with very few cells and invariant cleavage patterns and lineages. Ascidians belong to tunicates, the vertebrate sister group, and their study is essential to shed light on the emergence of vertebrates. Importantly, deciphering developmental gene regulatory networks has been carried out mostly in two of the three ascidian orders, Phlebobranchia and Stolidobranchia. To infer ancestral developmental programs in ascidians, it is thus essential to carry out molecular embryology in the third ascidian order, the Aplousobranchia. Here, we present genomic resources for the colonial aplousobranch Clavelina lepadiformis: a transcriptome produced from various embryonic stages, and an annotated genome. The assembly consists of 184 contigs making a total of 233.6 Mb with a N50 of 8.5 Mb and a L50 of 11. The 32,318 predicted genes capture 96.3% of BUSCO orthologs. We further show that these resources are suitable to study developmental gene expression and regulation in a comparative framework within ascidians. Additionally, they will prove valuable for evolutionary and ecological studies.

Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Roast: a tool for reference-free optimization of supertranscriptome assemblies

BACKGROUND

Transcriptomic studies involving organisms for which reference genomes are not available typically start by generating de novo transcriptome or supertranscriptome assembly from the raw RNA-seq reads. Assembling a supertranscriptome is, however, a challenging task due to significantly varying abundance of mRNA transcripts, alternative splicing, and sequencing errors. As a result, popular de novo supertranscriptome assembly tools generate assemblies containing contigs that are partially-assembled, fragmented, false chimeras or have local mis-assemblies leading to decreased assembly accuracy. Commonly available tools for assembly improvement rely primarily on running BLAST using closely related species making their accuracy and reliability conditioned on the availability of the data for closely related organisms.

RESULTS

We present ROAST, a tool for optimization of supertranscriptome assemblies that uses paired-end RNA-seq data from Illumina sequencing platform to iteratively identify and fix assembly errors solely using the error signatures generated by RNA-seq alignment tools including soft-clips, unexpected expression coverage, and reads with mates unmapped or mapped on a different contig to identify and fix various supertranscriptome assembly errors without performing BLAST searches against other organisms. Evaluation results using simulated as well as real datasets show that ROAST significantly improves assembly quality by identifying and fixing various assembly errors.

CONCLUSION

ROAST provides a reference-free approach to optimizing supertranscriptome assemblies highlighting its utility in refining de novo supertranscriptome assemblies of non-model organisms.

Macrophages modulate fibrosis during newt lens regeneration

Background

Previous studies indicated that macrophages play a role during lens regeneration in newts, but their function has not been tested experimentally.

Methods

Here we generated a transgenic newt reporter line in which macrophages can be visualized in vivo. Using this new tool, we analyzed the location of macrophages during lens regeneration. We uncovered early gene expression changes using bulk RNAseq in two newt species, Notophthalmus viridescens and Pleurodeles waltl. Next, we used clodronate liposomes to deplete macrophages, which inhibited lens regeneration in both newt species.

Results

Macrophage depletion induced the formation of scar-like tissue, an increased and sustained inflammatory response, an early decrease in iris pigment epithelial cell (iPEC) proliferation and a late increase in apoptosis. Some of these phenotypes persisted for at least 100 days and could be rescued by exogenous FGF2. Re-injury alleviated the effects of macrophage depletion and re-started the regeneration process.

Conclusions

Together, our findings highlight the importance of macrophages in facilitating a pro-regenerative environment in the newt eye, helping to resolve fibrosis, modulating the overall inflammatory landscape and maintaining the proper balance of early proliferation and late apoptosis.

De novo assembly and annotation of the Amblyomma hebraeum tick midgut transcriptome response to Ehrlichia ruminantium infection

The South African bont tick Amblyomma hebraeum is a hematophagous vector for the heartwater disease pathogen Ehrlichia ruminantium in southern Africa. During feeding, the tick's enterocytes express proteins that perform vital functions in blood digestion, including proteins that may be involved in E. ruminantium acquisition, colonization or immunity. To delineate the molecular mechanism of midgut response to E. ruminantium infection, we performed comparative analyses of midgut transcriptomes of E. ruminantium infected engorged A. hebraeum nymphs, and infected adult male and female ticks with their corresponding matched uninfected controls, before and during feeding. A total of 102,036 unigenes were annotated in public databases and their expression levels analyzed for engorged nymphs as well as unfed and partly-fed adult ticks. There were 2,025 differentially expressed genes (DEGs) in midguts, of which 1,225 unigenes were up-regulated and 800 unigenes were down-regulated in the midguts of infected ticks. Annotation of DEGs revealed an increase in metabolic and cellular processes among E. ruminantium infected ticks. Notably, among the infected ticks, there was up-regulation in the expression of genes involved in tick immunity, histone proteins and oxidative stress responses. We also observed up-regulation of glycoproteins that E. ruminantium could potentially use as docking sites for host cell entry. Insights uncovered in this study offer a platform for further investigations into the molecular interaction between E. ruminantium and A. hebraeum.

Brain transcriptome of gobies inhabiting natural CO2 seeps reveal acclimation strategies to long-term acidification

Ocean acidification (OA) is known to affect the physiology, survival, behaviour and fitness of various fish species with repercussions at the population, community and ecosystem levels. Some fish species, however, seem to acclimate rapidly to OA conditions and even thrive in acidified environments. The molecular mechanisms that enable species to successfully inhabit high CO2 environments have not been fully elucidated especially in wild fish populations. Here, we used the natural CO2 seep in Vulcano Island, Italy to study the effects of elevated CO2 exposure on the brain transcriptome of the anemone goby, a species with high population density in the CO2 seep and investigate their potential for acclimation. Compared to fish from environments with ambient CO2, gobies living in the CO2 seep showed differences in the expression of transcripts involved in ion transport and pH homeostasis, cellular stress, immune response, circadian rhythm and metabolism. We also found evidence of potential adaptive mechanisms to restore the functioning of GABAergic pathways, whose activity can be affected by exposure to elevated CO2 levels. Our findings indicate that gobies living in the CO2 seep may be capable of mitigating CO2-induced oxidative stress and maintaining physiological pH while meeting the consequent increased energetic costs. The conspicuous difference in the expression of core circadian rhythm transcripts could provide an adaptive advantage by increasing the flexibility of physiological processes in elevated CO2 conditions thereby facilitating acclimation. Our results show potential molecular processes of acclimation to elevated CO2 in gobies enabling them to thrive in the acidified waters of Vulcano Island.

Comparative phylotranscriptomics reveals ancestral and derived root nodule symbiosis programmes

Symbiotic interactions such as the nitrogen-fixing root nodule symbiosis (RNS) have structured ecosystems during the evolution of life. Here we aimed at reconstructing ancestral and intermediate steps that shaped RNS observed in extant flowering plants. We compared the symbiotic transcriptomic responses of nine host plants, including the mimosoid legume Mimosa pudica for which we assembled a chromosome-level genome. We reconstructed the ancestral RNS transcriptome composed of most known symbiotic genes together with hundreds of novel candidates. Cross-referencing with transcriptomic data in response to experimentally evolved bacterial strains with gradual symbiotic proficiencies, we found the response to bacterial signals, nodule infection, nodule organogenesis and nitrogen fixation to be ancestral. By contrast, the release of symbiosomes was associated with recently evolved genes encoding small proteins in each lineage. We demonstrate that the symbiotic response was mostly in place in the most recent common ancestor of the RNS-forming species more than 90 million years ago.

Macrophages modulate fibrosis during newt lens regeneration

Previous studies indicated that macrophages play a role during lens regeneration in newts, but their function has not been tested experimentally. Here we generated a transgenic newt reporter line in which macrophages can be visualized in vivo. Using this new tool, we analyzed the location of macrophages during lens regeneration. We uncovered early gene expression changes using bulk RNAseq in two newt species, Notophthalmus viridescens and Pleurodeles waltl. Next, we used clodronate liposomes to deplete macrophages, which inhibited lens regeneration in both newt species. Macrophage depletion induced the formation of scar-like tissue, an increased and sustained inflammatory response, an early decrease in iris pigment epithelial cell (iPEC) proliferation and a late increase in apoptosis. Some of these phenotypes persisted for at least 100 days and could be rescued by exogenous FGF2. Re-injury alleviated the effects of macrophage depletion and re-started the regeneration process. Together, our findings highlight the importance of macrophages in facilitating a pro-regenerative environment in the newt eye, helping to resolve fibrosis, modulating the overall inflammatory landscape and maintaining the proper balance of early proliferation and late apoptosis.

A secreted protease-like protein in Zymoseptoria tritici is responsible for avirulence on Stb9 resistance gene in wheat

Zymoseptoria tritici is the fungal pathogen responsible for Septoria tritici blotch on wheat. Disease outcome in this pathosystem is partly determined by isolate-specific resistance, where wheat resistance genes recognize specific fungal factors triggering an immune response. Despite the large number of known wheat resistance genes, fungal molecular determinants involved in such cultivar-specific resistance remain largely unknown. We identified the avirulence factor AvrStb9 using association mapping and functional validation approaches. Pathotyping AvrStb9 transgenic strains on Stb9 cultivars, near isogenic lines and wheat mapping populations, showed that AvrStb9 interacts with Stb9 resistance gene, triggering an immune response. AvrStb9 encodes an unusually large avirulence gene with a predicted secretion signal and a protease domain. It belongs to a S41 protease family conserved across different filamentous fungi in the Ascomycota class and may constitute a core effector. AvrStb9 is also conserved among a global Z. tritici population and carries multiple amino acid substitutions caused by strong positive diversifying selection. These results demonstrate the contribution of an 'atypical' conserved effector protein to fungal avirulence and the role of sequence diversification in the escape of host recognition, adding to our understanding of host-pathogen interactions and the evolutionary processes underlying pathogen adaptation.

De Novo Transcriptome Assembly and EST-SSR Marker Development and Application in Chrysosplenium macrophyllum

Chrysosplenium macrophyllum Oliv., belonging to the family Saxifragaceae, is a traditional and unique Chinese herbal medicine. However, the lack of adequate molecular markers has hampered the progress regarding population genetics and evolution within this species. In this research, we used the DNBSEQ-T7 Sequencer (MGI) sequencing assay to analyze the transcriptome profiles of C. macrophyllum. SSR markers were developed on the basis of transcriptomic sequences and further validated on C. macrophyllum and other Chrysosplenium species. The genetic diversity and structure of the 12 populations were analyzed by using polymorphic expressed sequence tag simple sequence repeat (EST-SSR) markers. A potential pool of 3127 non-redundant EST-SSR markers were identified for C. macrophyllum in this study. The developed EST-SSR markers had high amplification rates and cross-species transferability in Chrysosplenium. Our results also showed that the natural populations of C. macrophyllum had a high level of genetic diversity. Genetic distance, principal component analysis, and popular structure analysis revealed that all 60 samples clustered into two major groups that were consistent with their geographical origins. This study provided a batch of highly polymorphic EST-SSR molecular markers that were developed via transcriptome sequencing. These markers will be of great significance for the study of the genetic diversity and evolutionary history of C. macrophyllum and other Chrysosplenium species.

Venomics survey of six myrmicine ants provides insights into the molecular and structural diversity of their peptide toxins

Among ants, Myrmicinae represents the most speciose subfamily. The venom composition previously described for these social insects is extremely variable, with alkaloids predominant in some genera while, conversely, proteomics studies have revealed that some myrmicine ant venoms are peptide-rich. Using integrated transcriptomic and proteomic approaches, we characterized the venom peptidomes of six ants belonging to the different tribes of Myrmicinae. We identified a total of 79 myrmicitoxins precursors which can be classified into 38 peptide families according to their mature sequences. Myrmicine ant venom peptidomes showed heterogeneous compositions, with linear and disulfide-bonded monomers as well as dimeric toxins. Several peptide families were exclusive to a single venom whereas some were retrieved in multiple species. A hierarchical clustering analysis of precursor signal sequences led us to divide the myrmicitoxins precursors into eight families, including some that have already been described in other aculeate hymenoptera such as secapin-like peptides and voltage-gated sodium channel (Na_V) toxins. Evolutionary and structural analyses of two representatives of these families highlighted variation and conserved patterns that might be crucial to explain myrmicine venom peptide functional adaptations to biological targets.

Rounding up the annual ryegrass genome: High-quality reference genome of Lolium rigidum

The genome of the major agricultural weed species, annual ryegrass (Lolium rigidum) was assembled, annotated and analysed. Annual ryegrass is a major weed in grain cropping, and has the remarkable capacity to evolve resistance to herbicides with various modes of action. The chromosome-level assembly was achieved using short- and long-read sequencing in combination with Hi-C mapping. The assembly size is 2.44 Gb with N₅₀ = 361.79 Mb across 1,764 scaffolds where the seven longest sequences correspond to the seven chromosomes. Genome completeness assessed through BUSCO returned a 99.8% score for complete (unique and duplicated) and fragmented genes using the Viridiplantae set. We found evidence for the expansion of herbicide resistance-related gene families including detoxification genes. The reference genome of L. rigidum is a critical asset for leveraging genetic information for the management of this highly problematic weed species.

Comparative Transcriptomics for Genes Related to Berberine and Berbamine Biosynthesis in Berberidaceae

Berberine and berbamine are bioactive compounds of benzylisoquinoline alkaloids (BIAs) present in Berberis species. The contents of berbamine are 20 times higher than berberine in leaf tissues in three closely related species: Berberis koreana, B. thunbergii and B. amurensis. This is the first report on the quantification of berberine compared to the berbamine in the Berberis species. Comparative transcriptome analyses were carried out with mRNAs from the leaf tissues of the three-species. The comparison of the transcriptomes of B. thunbergii and B. amurensis to those of B. koreana, B. thunbergii showed a consistently higher number of differentially expressed genes than B. amurensis in KEGG and DEG analyses. All genes encoding enzymes involved in berberine synthesis were identified and their expressions were variable among the three species. There was a single copy of CYP80A/berbamunine synthase in B. koreana. Methyltransferases and cytochrome P450 mono-oxidases (CYPs) are key enzymes for BIA biosynthesis. The current report contains the copy numbers and other genomic characteristics of the methyltransferases and CYPs in Berberis species. Thus, the contents of the current research are valuable for molecular characterization for the medicinal utilization of the Berberis species.

A CLAVATA3-like Gene Acts as a Gynoecium Suppression Function in White Campion

How do separate sexes originate and evolve? Plants provide many opportunities to address this question as they have diverse mating systems and separate sexes (dioecy) that evolved many times independently. The classic "two-factor" model for evolution of separate sexes proposes that males and females can evolve from hermaphrodites via the spread of male and female sterility mutations that turn hermaphrodites into females and males, respectively. This widely accepted model was inspired by early genetic work in dioecious white campion (Silene latifolia) that revealed the presence of two sex-determining factors on the Y-chromosome, though the actual genes remained unknown. Here, we report identification and functional analysis of the putative sex-determining gene in S. latifolia, corresponding to the gynoecium suppression factor (GSF). We demonstrate that GSF likely corresponds to a Y-linked CLV3-like gene that is specifically expressed in early male flower buds and encodes the protein that suppresses gynoecium development in S. latifolia. Interestingly, GSFY has a dysfunctional X-linked homolog (GSFX) and their synonymous divergence (dS = 17.9%) is consistent with the age of sex chromosomes in this species. We propose that female development in S. latifolia is controlled via the WUSCHEL-CLAVATA feedback loop, with the X-linked WUSCHEL-like and Y-linked CLV3-like genes, respectively. Evolution of dioecy in the S. latifolia ancestor likely involved inclusion of ancestral GSFY into the nonrecombining region on the nascent Y-chromosome and GSFX loss of function, which resulted in disbalance of the WUSCHEL-CLAVATA feedback loop between the sexes and ensured gynoecium suppression in males.

De Novo Transcriptome Assembly and Analysis of the Flat Oyster Pathogenic Protozoa Bonamia Ostreae

The flat oyster Ostrea edulis is an oyster species native to Europe. It has declined to functional extinction in many areas of the NE Atlantic for several decades. Factors explaining this decline include over-exploitation of natural populations and diseases like bonamiosis, regulated across both the EU and the wider world and caused by the intracellular protozoan parasite Bonamia ostreae. To date, very limited sequence data are available for this Haplosporidian species. We present here the first transcriptome of B. ostreae. As this protozoan is not yet culturable, it remains extremely challenging to obtain high-quality -omic data. Thanks to a specific parasite isolation protocol and a dedicated bioinformatic pipeline, we were able to obtain a high-quality transcriptome for an intracellular marine micro-eukaryote, which will be very helpful to better understand its biology and to consider the development of new relevant diagnostic tools.

De novo transcriptome assembly and annotation of the third stage larvae of the zoonotic parasite Anisakis pegreffii

OBJECTIVES

Anisakis pegreffii is a zoonotic parasite requiring marine organisms to complete its life-history. Human infection (anisakiasis) occurs when the third stage larvae (L3) are accidentally ingested with raw or undercooked infected fish or squids. A new de novo transcriptome of A. pegreffii was here generated aiming to provide a robust bulk of data to be used for a comprehensive "ready-to-use" resource for detecting functional studies on genes and gene products of A. pegreffii involved in the molecular mechanisms of parasite-host interaction.

DATA DESCRIPTION

A RNA-seq library of A. pegreffii L3 was here newly generated by using Illumina TruSeq platform. It was combined with other five RNA-seq datasets previously gathered from L3 of the same species stored in SRA of NCBI. The final dataset was analyzed by launching three assembler programs and two validation tools. The use of a robust pipeline produced a high-confidence protein-coding transcriptome of A. pegreffii. These data represent a more robust and complete transcriptome of this species with respect to the actually existing resources. This is of importance for understanding the involved adaptive and immunomodulatory genes implicated in the "cross talk" between the parasite and its hosts, including the accidental one (humans).

Transcriptome of the synganglion in the tick Ixodes ricinus and evolution of the cys-loop ligand-gated ion channel family in ticks

BACKGROUND

Ticks represent a major health issue for humans and domesticated animals. Exploring the expression landscape of the tick's central nervous system (CNS), known as the synganglion, would be an important step in understanding tick physiology and in managing tick-borne diseases, but studies on that topic are still relatively scarce. Neuron-specific genes like the cys-loop ligand-gated ion channels (cys-loop LGICs, or cysLGICs) are important pharmacological targets of acaricides. To date their sequence have not been well catalogued for ticks, and their phylogeny has not been fully studied.

RESULTS

We carried out the sequencing of transcriptomes of the I. ricinus synganglion, for adult ticks in different conditions (unfed males, unfed females, and partially-fed females). The de novo assembly of these transcriptomes allowed us to obtain a large collection of cys-loop LGICs sequences. A reference meta-transcriptome based on synganglion and whole body transcriptomes was then produced, showing high completeness and allowing differential expression analyses between synganglion and whole body. Many of the genes upregulated in the synganglion were associated with neurotransmission and/or localized in neurons or the synaptic membrane. As the first step of a functional study of cysLGICs, we cloned the predicted sequence of the resistance to dieldrin (RDL) subunit homolog, and functionally reconstituted the first GABA-gated receptor of Ixodes ricinus. A phylogenetic study was performed for the nicotinic acetylcholine receptors (nAChRs) and other cys-loop LGICs respectively, revealing tick-specific expansions of some types of receptors (especially for Histamine-like subunits and GluCls).

CONCLUSIONS

We established a large catalogue of genes preferentially expressed in the tick CNS, including the cysLGICs. We discovered tick-specific gene family expansion of some types of cysLGIC receptors, and a case of intragenic duplication, suggesting a complex pattern of gene expression among different copies or different alternative transcripts of tick neuro-receptors.

Improving the Annotation of the Venom Gland Transcriptome of Pamphobeteus verdolaga, Prospecting Novel Bioactive Peptides

Spider venoms constitute a trove of novel peptides with biotechnological interest. Paucity of next-generation-sequencing (NGS) data generation has led to a description of less than 1% of these peptides. Increasing evidence supports the underestimation of the assembled genes a single transcriptome assembler can predict. Here, the transcriptome of the venom gland of the spider Pamphobeteus verdolaga was re-assembled, using three free access algorithms, Trinity, SOAPdenovo-Trans, and SPAdes, to obtain a more complete annotation. Assembler’s performance was evaluated by contig number, N50, read representation on the assembly, and BUSCO’s terms retrieval against the arthropod dataset. Out of all the assembled sequences with all software, 39.26% were common between the three assemblers, and 27.88% were uniquely assembled by Trinity, while 27.65% were uniquely assembled by SPAdes. The non-redundant merging of all three assemblies’ output permitted the annotation of 9232 sequences, which was 23% more when compared to each software and 28% more when compared to the previous P. verdolaga annotation; moreover, the description of 65 novel theraphotoxins was possible. In the generation of data for non-model organisms, as well as in the search for novel peptides with biotechnological interest, it is highly recommended to employ at least two different transcriptome assemblers.

Rapid evolution fuels transcriptional plasticity to ocean acidification

Ocean acidification (OA) is postulated to affect the physiology, behavior, and life-history of marine species, but potential for acclimation or adaptation to elevated pCO₂ in wild populations remains largely untested. We measured brain transcriptomes of six coral reef fish species at a natural volcanic CO₂  seep and an adjacent control reef in Papua New Guinea. We show that elevated pCO₂ induced common molecular responses related to circadian rhythm and immune system but different magnitudes of molecular response across the six species. Notably, elevated transcriptional plasticity was associated with core circadian genes affecting the regulation of intracellular pH and neural activity in Acanthochromis polyacanthus. Gene expression patterns were reversible in this species as evidenced upon reduction of CO₂ following a natural storm-event. Compared with other species, Ac. polyacanthus has a more rapid evolutionary rate and more positively selected genes in key functions under the influence of elevated CO₂ , thus fueling increased transcriptional plasticity. Our study reveals the basis to variable gene expression changes across species, with some species possessing evolved molecular toolkits to cope with future OA.

Transcriptome Analysis Reveals the Molecular Response to Salinity Challenge in Larvae of the Giant Freshwater Prawn Macrobrachium rosenbergii

Salinity is a crucial factor influencing the growth, development, immunity, and reproduction of aquatic organisms; however, little is known about the molecular mechanism of the response to salinity challenge in larvae of the giant freshwater prawn Macrobrachium rosenbergii. Herein, larvae cultured in three treatment groups with salinities of 10, 13, and 16‰ (S10, S13, and S16) were collected, and then transcriptome analysis was conducted by RNA-seq. A total of 6,473, 3,830 and 3,584 differentially expressed genes (DEGs) were identified in the S10 vs. S13 comparison, S10 vs. S16 comparison and S13 vs. S16 comparison, respectively. These genes are involved in osmoregulation, energy metabolism, molting, and the immune response. qPCR analysis was used to detect the expression patterns of 16 DEGs to verify the accuracy of the transcriptome data. Protein–protein interaction (PPI) analysis for DEGs and microsatellite marker screening were also conducted to reveal the molecular mechanism of salinity regulation. Together, our results will provide insight into the molecular genetic basis of adaptation to salinity challenge for larvae of M. rosenbergii.

Functional R2R3-MYB transcription factor NsMYB1, regulating anthocyanin biosynthesis, was relative to the fruit color differentiation in Nitraria sibirica Pall

BACKGROUND

Nitraria sibirica Pall. is an economic plant with two kinds of fruit color, widely spreads in the Qinghai Tibet Plateau. The chemical analysis and pharmacological evaluation had been carried out for several tens of years, the mechanism behind the fruit color differentiation is still unclear.

RESULTS

In this manuscript, the chemical analysis of the extractions showed that the chemical composition of fruit color was anthocyanin, and two kind of Nitraria sibirica Pall. were caused by the content differentiation with the same anthocyanin kinds. Cyanidin-3-[2"-(6'"-coumaroyl)-glucosyl]-glucoside (C3G) was the major anthocyanin. Transcriptome analysis and the qRT-PCR revealed that the structural genes relative to anthocyanin biosynthesis except CHS, F3'5'H and ANS were up-regulated in the peels of BF (Black fruit) compared with the peels of RF (Red fruit), which indicated that transcript factor should be the reason for the expression difference of the structure genes. In the unigenes of the transcript factor MYB and bHLH, relative to anthocyanin, only NsMYB1 (Cluster 8422.10600), was high-expression and up-expression in the peels of BF. NsMYB1 encoded the same length protein with four amino acid differences in the RF and BF, and both contained the intact DNA, HTH-MYB and SANT domains. NsMYB1 was close to the AtMYB114, AtMYB113 and AtPAP1, regulating anthocyanin biosynthesis, in phylogenetic relationship. Both NsMYB1r and NsMYB1b could promote the transcript of the structural genes, and induced the anthocyanin accumulation in all tissues of transgenic tobacco. The insertion of 'TATA' in the promoter of NsMYB1r gave one more promoter region, and was the reason for higher transcripts in black fruit possibly.

CONCLUSIONS

Cyanidin-3-[2''-(6'"-coumaroyl)-glucosyl]-glucoside was the major anthocyanin in black fruit of Nitraria sibirica Pall.. NsMYB1 was a functional R2R3-MYB transcription factor, regulated the anthocyanin biosynthesis, and led to the fruit color differentiation in Nitraria sibirica Pall.

Chromosome-Level Genome Assembly Reveals Dynamic Sex Chromosomes in Neotropical Leaf-Litter Geckos (Sphaerodactylidae: Sphaerodactylus)

Sex determination is a critical element of successful vertebrate development, suggesting that sex chromosome systems might be evolutionarily stable across lineages. For example, mammals and birds have maintained conserved sex chromosome systems over long evolutionary time periods. Other vertebrates, in contrast, have undergone frequent sex chromosome transitions, which is even more amazing considering we still know comparatively little across large swaths of their respective phylogenies. One reptile group in particular, the gecko lizards (infraorder Gekkota), shows an exceptional lability with regard to sex chromosome transitions and may possess the majority of transitions within squamates (lizards and snakes). However, detailed genomic and cytogenetic information about sex chromosomes is lacking for most gecko species, leaving large gaps in our understanding of the evolutionary processes at play. To address this, we assembled a chromosome-level genome for a gecko (Sphaerodactylidae: Sphaerodactylus) and used this assembly to search for sex chromosomes among six closely related species using a variety of genomic data, including whole-genome re-sequencing, RADseq, and RNAseq. Previous work has identified XY systems in two species of Sphaerodactylus geckos. We expand upon that work to identify between two and four sex chromosome cis-transitions (XY to a new XY) within the genus. Interestingly, we confirmed two different linkage groups as XY sex chromosome systems that were previously unknown to act as sex chromosomes in tetrapods (syntenic with Gallus chromosome 3 and Gallus chromosomes 18/30/33), further highlighting a unique and fascinating trend that most linkage groups have the potential to act as sex chromosomes in squamates.

Unstable Relationship Between Braarudosphaera bigelowii (= Chrysochromulina parkeae) and Its Nitrogen-Fixing Endosymbiont

Marine phytoplankton are major primary producers, and their growth is primarily limited by nitrogen in the oligotrophic ocean environment. The haptophyte Braarudosphaera bigelowii possesses a cyanobacterial endosymbiont (UCYN-A), which plays a major role in nitrogen fixation in the ocean. However, host-symbiont interactions are poorly understood because B. bigelowii was unculturable. In this study, we sequenced the complete genome of the B. bigelowii endosymbiont and showed that it was highly reductive and closely related to UCYN-A2 (an ecotype of UCYN-A). We succeeded in establishing B. bigelowii strains and performed microscopic observations. The detailed observations showed that the cyanobacterial endosymbiont was surrounded by a single host derived membrane and divided synchronously with the host cell division. The transcriptome of B. bigelowii revealed that B. bigelowii lacked the expression of many essential genes associated with the uptake of most nitrogen compounds, except ammonia. During cultivation, some of the strains completely lost the endosymbiont. Moreover, we did not find any evidence of endosymbiotic gene transfer from the endosymbiont to the host. These findings illustrate an unstable morphological, metabolic, and genetic relationship between B. bigelowii and its endosymbiont.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Transcriptome Profiling of the Pacific Oyster Crassostrea gigas Visceral Ganglia over a Reproduction Cycle Identifies Novel Regulatory Peptides

The neuropeptides involved in the regulation of reproduction in the Pacific oyster (Crassostrea gigas) are quite diverse. To investigate this diversity, a transcriptomic survey of the visceral ganglia (VG) was carried out over an annual reproductive cycle. RNA-seq data from 26 samples corresponding to VG at different stages of reproduction were de novo assembled to generate a specific reference transcriptome of the oyster nervous system and used to identify differentially expressed transcripts. Transcriptome mining led to the identification of novel neuropeptide precursors (NPPs) related to the bilaterian Eclosion Hormone (EH), crustacean female sex hormone/Interleukin 17, Nesfatin, neuroparsin/IGFBP, prokineticins, and urotensin I; to the protostome GNQQN, pleurin, prohormones 3 and 4, prothoracotropic hormones (PTTH), and QSamide/PXXXamide; to the lophotrochozoan CCWamide, CLCCY, HFAamide, and LXRX; and to the mollusk-specific NPPs CCCGS, clionin, FYFY, GNamide, GRWRN, GSWN, GWE, IWMPxxGYxx, LXRYamide, RTLFamide, SLRFamide, and WGAGamide. Among the complete repertoire of NPPs, no sex-biased expression was observed. However, 25 NPPs displayed reproduction stage-specific expression, supporting their involvement in the control of gametogenesis or associated metabolisms.

Lipid exchanges drove the evolution of mutualism during plant terrestrialization

Symbiosis with arbuscular mycorrhizal fungi (AMF) improves plant nutrition in most land plants, and its contribution to the colonization of land by plants has been hypothesized. Here, we identify a conserved transcriptomic response to AMF among land plants, including the activation of lipid metabolism. Using gain of function, we show the transfer of lipids from the liverwort Marchantia paleacea to AMF and its direct regulation by the transcription factor WRINKLED (WRI). Arbuscules, the nutrient-exchange structures, were not formed in loss-of-function wri mutants in M. paleacea, leading to aborted mutualism. Our results show the orthology of the symbiotic transfer of lipids across land plants and demonstrate that mutualism with arbuscular mycorrhizal fungi was present in the most recent ancestor of land plants 450 million years ago.

A chromosome-level genome assembly of the parasitoid wasp, Cotesia glomerata (Hymenoptera: Braconidae)

Hymenopterans make up about 20% of all animal species, but most are poorly known and lack high-quality genomic resources. One group of important, yet under-studied hymenopterans, are parasitoid wasps in the family Braconidae. Among this under-studied group are braconid wasps in the genus Cotesia; a clade of ~1,000 species routinely used in studies of physiology, ecology, biological control, and genetics. However, our ability to understand these organisms has been hindered by a lack of genomic resources. We helped bridge this gap by generating a high-quality genome assembly for the parasitoid wasp, Cotesia glomerata (Braconidae; Microgastrinae). We generated this assembly using multiple sequencing technologies, including Oxford Nanopore, whole-genome shotgun sequencing, and 3-D chromatin contact information (Hi-C). Our assembly is one of the most contiguous, complete, and publicly available hymenopteran genomes, represented by 3,355 scaffolds with a scaffold N50 of ~28Mb and a BUSCO score of ~99%. Given the genome sizes found in closely related species, our genome assembly was ~50% larger than expected, which was apparently induced by runaway amplification of three types of repetitive elements: simple repeats, Long Terminal Repeats (LTRs), and Long Interspersed Nuclear Elements (LINEs). This assembly is another step forward for genomics across this hyper-diverse, yet understudied, order of insects. The assembled genomic data and metadata files are publicly available via Figshare (https://doi.org/10.6084/m9.figshare.13010549).

Quick and efficient approach to develop genomic resources in orphan species: Application in Lavandula angustifolia

Next-Generation Sequencing (NGS) technologies, by reducing the cost and increasing the throughput of sequencing, have opened doors to generate genomic data in a range of previously poorly studied species. In this study, we propose a method for the rapid development of a large-scale molecular resources for orphan species. We studied as an example the true lavender (Lavandula angustifolia Mill.), a perennial sub-shrub plant native from the Mediterranean region and whose essential oil have numerous applications in cosmetics, pharmaceuticals, and alternative medicines. The heterozygous clone “Maillette” was used as a reference for DNA and RNA sequencing. We first built a reference Unigene, compound of coding sequences, thanks to de novo RNA-seq assembly. Then, we reconstructed the complete genes sequences (with introns and exons) using an Unigene-guided DNA-seq assembly approach. This aimed to maximize the possibilities of finding polymorphism between genetically close individuals despite the lack of a reference genome. Finally, we used these resources for SNP mining within a collection of 16 commercial lavender clones and tested the SNP within the scope of a genetic distance analysis. We obtained a cleaned reference of 8, 030 functionally in silico annotated genes. We found 359K polymorphic sites and observed a high SNP frequency (mean of 1 SNP per 90 bp) and a high level of heterozygosity (more than 60% of heterozygous SNP per genotype). On overall, we found similar genetic distances between pairs of clones, which is probably related to the out-crossing nature of the species and the restricted area of cultivation. The proposed method is transferable to other orphan species, requires little bioinformatics resources and can be realized within a year. This is also the first reported large-scale SNP development on Lavandula angustifolia. All the genomics resources developed herein are publicly available and provide a rich pool of molecular resources to explore and exploit lavender genetic diversity in breeding programs.

Transfer and Transcriptomic Profiling in Liver and Brain of European Eels (Anguilla anguilla) After Diet-borne Exposure to Gold Nanoparticles

A nanometric revolution is underway, promising technical innovations in a wide range of applications and leading to a potential boost in environmental discharges. The propensity of nanoparticles (NPs) to be transferred throughout trophic chains and to generate toxicity was mainly assessed in primary consumers, whereas a lack of knowledge for higher trophic levels persists. The present study focused on a predatory fish, the European eel (Anguilla anguilla) exposed to gold NPs (AuNPs; 10 nm, polyethylene glycol-coated) for 21 d at 3 concentration levels in food: 0 (NP0), 1 (NP1), and 10 (NP10) mg Au kg^-1 . Transfer was assessed by Au quantification in eel tissues, and transcriptomic responses in the liver and brain were revealed by a high-throughput RNA-sequencing approach. Eels fed at NP10 presented an erratic feeding behavior, whereas Au quantification only indicated transfer to intestine and kidney of NP1-exposed eels. Sequencing of RNA was performed in NP0 and NP1 eels. A total of 258 genes and 156 genes were significantly differentially transcribed in response to AuNP trophic exposure in the liver and brain, respectively. Enrichment analysis highlighted modifications in the immune system-related processes in the liver. In addition, results pointed out a shared response of both organs regarding 13 genes, most of them being involved in immune functions. This finding may shed light on the mode of action and toxicity of AuNPs in fish. Environ Toxicol Chem 2020;39:2450-2461. © 2020 SETAC.

Siberian sturgeon multi-tissue reference transcriptome database

Motivation: Siberian sturgeon is a long lived and late maturing fish farmed for caviar production in 50 countries. Functional genomics enable to find genes of interest for fish farming. In the absence of a reference genome, a reference transcriptome is very useful for sequencing based functional studies.

Results: We present here a high-quality transcriptome assembly database built using RNA-seq reads coming from brain, pituitary, gonadal, liver, stomach, kidney, anterior kidney, heart, embryonic and pre-larval tissues. It will facilitate crucial research on topics such as puberty, reproduction, growth, food intake and immunology. This database represents a major contribution to the publicly available sturgeon transcriptome reference datasets.

Availability: The database is publicly available at http://siberiansturgeontissuedb.sigenae.org

Supplementary information:  Supplementary data are available at Database online.

Characterization and complete genome sequences of two novel variants of the family Closteroviridae from Chinese kiwifruit

Two distinct closterovirus-like genome sequences (termed AdV-1 v1 and v2) were identified in Actinidia chinensis var. deliciosa 'Miliang-1' that had no disease symptoms using high-throughput sequencing. Using overlapping reverse transcription-polymerase chain reaction and rapid amplification of cDNA ends, the genomic sequences of AdV-1 v1 and v2 were confirmed as 17,646 and 18,578 nucleotides in length, respectively. The two complete genomes contained 9 and 15 open reading frames, respectively, coding for proteins having domains typical of Closteroviridae, such as RNA-dependent RNA polymerase (RdRp), heat shock protein 70 homolog (HSP70h) and coat protein (CP). Sequence analysis showed that the amino acid sequences of RdRp, HSP70h, and CP of the two variants exhibited high similarity (> 80%), while their genomic organization was somewhat different. This suggested that the two viral genomes identified here are variants of the family Closteroviridae in a single kiwifruit host. Furthermore, phylogenetic relationship analysis revealed that the two variants had a closer relationship with the unclassified virus Persimmon virus B (PeVB) and Actinidia virus 1 (AcV-1) than with other members of the family Closteroviridae, as did their genomic organization. It is speculated that the two variants, together with PeVB and AcV-1 belong to a new subfamily of Closteroviridae.

Development of a relevant strategy using de novo transcriptome assembly method for transcriptome comparisons between Muscovy and common duck species and their reciprocal inter-specific mule and hinny hybrids fed ad libitum and overfed

Background

Common Pekin and Muscovy ducks and their intergeneric hinny and mule hybrids have different abilities for fatty liver production. RNA-Seq analyses from the liver of these different genetic types fed ad libitum or overfed would help to identify genes with different response to overfeeding between them. However RNA-seq analyses from different species and comparison is challenging. The goal of this study was develop a relevant strategy for transcriptome analysis and comparison between different species.

Results

Transcriptomes were first assembled with a reference-based approach. Important mapping biases were observed when heterologous mapping were conducted on common duck reference genome, suggesting that this reference-based strategy was not suited to compare the four different genetic types. De novo transcriptome assemblies were then performed using Trinity and Oases. Assemblies of transcriptomes were not relevant when more than a single genetic type was considered. Finally, single genetic type transcriptomes were assembled with DRAP in a mega-transcriptome. No bias was observed when reads from the different genetic types were mapped on this mega-transcriptome and differences in gene expression between the four genetic types could be identified.

Conclusions

Analyses using both reference-based and de novo transcriptome assemblies point out a good performance of the de novo approach for the analysis of gene expression in different species. It also allowed the identification of differences in responses to overfeeding between Pekin and Muscovy ducks and hinny and mule hybrids.

The development of early pioneer neurons in the annelid Malacoceros fuliginosus

BACKGROUND

Nervous system development is an interplay of many processes: the formation of individual neurons, which depends on whole-body and local patterning processes, and the coordinated growth of neurites and synapse formation. While knowledge of neural patterning in several animal groups is increasing, data on pioneer neurons that create the early axonal scaffold are scarce. Here we studied the first steps of nervous system development in the annelid Malacoceros fuliginosus.

RESULTS

We performed a dense expression profiling of a broad set of neural genes. We found that SoxB expression begins at 4 h postfertilization, and shortly later, the neuronal progenitors can be identified at the anterior and the posterior pole by the transient and dynamic expression of proneural genes. At 9 hpf, the first neuronal cells start differentiating, and we provide a detailed description of axonal outgrowth of the pioneer neurons that create the primary neuronal scaffold. Tracing back the clonal origin of the ventral nerve cord pioneer neuron revealed that it is a descendant of the blastomere 2d (2d²²¹), which after 7 cleavages starts expressing Neurogenin, Acheate-Scute and NeuroD.

CONCLUSIONS

We propose that an anterior and posterior origin of the nervous system is ancestral in annelids. We suggest that closer examination of the first pioneer neurons will be valuable in better understanding of nervous system development in spirally cleaving animals, to determine the potential role of cell-intrinsic properties in neuronal specification and to resolve the evolution of nervous systems.

Domestication modulates the expression of genes involved in neurogenesis in high-quality eggs of Sander lucioperca

Pikeperch, Sander lucioperca, is a species of high interest to the aquaculture. The expansion of its production can only be achieved by furthering domestication level. However, the mechanisms driving the domestication process in finfishes are poorly understood. Transcriptome profiling of eggs was found to be a useful tool allowing understanding of the domestication process in teleosts. In this study, using next-generation sequencing, the first pikeperch transcriptome has been generated as well as pikeperch-specific microarray comprising 35,343 unique probes. Next, we performed transcriptome profiling of eggs obtained from wild and domesticated populations. We found 710 differentially expressed genes that were linked mostly to nervous system development. These results provide new insights into processes that are directly involved in the domestication of finfishes. It can be suggested that all the identified processes were predetermined by the maternally derived set of genes contained in the unfertilized eggs. This allows us to suggest that fish behavior, along with many other processes, can be predetermined at the cellular level and may have significant implications on the adaptation of cultured fish to the natural environment. This also allows to suggest that fish behavior should be considered as a very important pikeperch aquaculture selection trait.

Analysis of Natural Selection of Immune Genes in Spinibarbus caldwelli by Transcriptome Sequencing

Spinibarbus caldwelli is an omnivorous cyprinid fish that is distributed widely in China. To investigate the adaptive evolution of S. caldwelli, the muscle transcriptome was sequenced by Illumina HiSeq 4000 platform. A total of 80,447,367 reads were generated by next-generation sequencing. Also, 211,386 unigenes were obtained by de novo assembly. Additionally, we calculated that the divergence time between S. caldwelli and Sinocyclocheilus grahami is 23.14 million years ago (Mya). And both of them diverged from Ctenopharyngodon idellus 46.95 Mya. Furthermore, 38 positive genes were identified by calculating Ka/Ks ratios from 9225 orthologs. Among them, several immune-related genes were identified as positively selected, such as POLR3B, PIK3C3, TOPORS, FASTKD3, CYPLP1A1, and UACA. Our results throw light on the nature of the natural selection of S. caldwelli and contribute to future immunological and transcriptome studies.

De novo assembly and annotation of transcriptomes from two cultivars of Cannabis sativa with different cannabinoid profiles

Cannabis has been cultivated for millennia for medicinal, industrial and recreational uses. Our long-term goal is to compare the transcriptomes of cultivars with different cannabinoid profiles for therapeutic purposes. Here we describe the de novo assembly, annotation and initial analysis of two cultivars of Cannabis, a high THC variety and a CBD plus THC variety. Cultivars were grown under different lighting conditions; flower buds were sampled over 71 days. Cannabinoid profiles were determined by ESI-LC/MS. RNA samples were sequenced using the HiSeq4000 platform. Transcriptomes were assembled using the DRAP pipeline and annotated using the BLAST2GO pipeline and other tools. Each transcriptome contained over twenty thousand protein encoding transcripts with ORFs and flanking sequence. Identification of transcripts for cannabinoid pathway and related enzymes showed full-length ORFs that align with the draft genomes of the Purple Kush and Finola cultivars. Two transcripts were found for olivetolic acid cyclase (OAC) that mapped to distinct locations on the Purple Kush genome suggesting multiple genes for OAC are expressed in some cultivars. The ability to make high quality annotated reference transcriptomes in Cannabis or other plants can promote rapid comparative analysis between cultivars and growth conditions in Cannabis and other organisms without annotated genome assemblies.

Biological rhythms in the deep-sea hydrothermal mussel Bathymodiolus azoricus

Biological rhythms are a fundamental property of life. The deep ocean covers 66% of our planet surface and is one of the largest biomes. The deep sea has long been considered as an arrhythmic environment because sunlight is totally absent below 1,000 m depth. In the present study, we have sequenced the temporal transcriptomes of a deep-sea species, the ecosystem-structuring vent mussel Bathymodiolus azoricus. We reveal that tidal cycles predominate in the transcriptome and physiology of mussels fixed directly at hydrothermal vents at 1,688 m depth at the Mid-Atlantic Ridge, whereas daily cycles prevail in mussels sampled after laboratory acclimation. We identify B. azoricus canonical circadian clock genes, and show that oscillations observed in deep-sea mussels could be either a direct response to environmental stimulus, or be driven endogenously by one or more biological clocks. This work generates in situ insights into temporal organisation in a deep-sea organism.

Little is known about gene expression of organisms in the deep sea, partially owing to constraints on sampling these organisms in situ. Here the authors circumvent this problem, fixing tissue of a deep-sea mussel at 1,688 m in depth, and later analyzing transcriptomes to reveal gene expression patterns showing tidal oscillations.

Multi-tissue Siberian sturgeon RNA sequencing data

Siberian sturgeon, Acipenser baerii, is a commercially valuable fish for flesh and caviar production and a threatened species. We produced transcriptomic data for ten tissues with relevance to puberty, reproduction, early development, growth and food intake. The data includes RNA-Seq read sets of brain, pituitary, anterior-kidney, kidney, stomach, liver, heart, embryonic, pre-larval, and immature gonad sequences. Tissues were collected from sex differentiated fish (17 to 42 months of age, 66 to 85 cm) RNA was extracted and sequenced. Our purpose is to facilitate fundamental studies of sturgeon physiology to wild and aquaculture populations management.

Venom Peptide Repertoire of the European Myrmicine Ant Manica rubida: Identification of Insecticidal Toxins

Using an integrated transcriptomic and proteomic approach, we characterized the venom peptidome of the European red ant, Manica rubida. We identified 13 "myrmicitoxins" that share sequence similarities with previously identified ant venom peptides, one of them being identified as an EGF-like toxin likely resulting from a threonine residue modified by O-fucosylation. Furthermore, we conducted insecticidal assays of reversed-phase HPLC venom fractions on the blowfly Lucilia caesar, permitting us to identify six myrmicitoxins (i.e., U₃-, U₁₀-, U₁₃-, U₂₀-MYRTX-Mri1a, U₁₀-MYRTX-Mri1b, and U₁₀-MYRTX-Mri1c) with an insecticidal activity. Chemically synthesized U₁₀-MYRTX-Mri1a, -Mri1b, -Mri1c, and U₂₀-MYRTX-Mri1a irreversibly paralyzed blowflies at the highest doses tested (30-125 nmol·g^-1). U₁₃-MYRTX-Mri1a, the most potent neurotoxic peptide at 1 h, had reversible effects after 24 h (150 nmol·g^-1). Finally, U₃-MYRTX-Mri1a has no insecticidal activity, even at up to 55 nmol·g^-1. Thus, M. rubida employs a paralytic venom rich in linear insecticidal peptides, which likely act by disrupting cell membranes.

Genome Sequence of the Euryhaline Javafish Medaka, Oryzias javanicus: A Small Aquarium Fish Model for Studies on Adaptation to Salinity

The genus Oryzias consists of 35 medaka-fish species each exhibiting various ecological, morphological and physiological peculiarities and adaptations. Beyond of being a comprehensive phylogenetic group for studying intra-genus evolution of several traits like sex determination, behavior, morphology or adaptation through comparative genomic approaches, all medaka species share many advantages of experimental model organisms including small size and short generation time, transparent embryos and genome editing tools for reverse and forward genetic studies. The Java medaka, Oryzias javanicus, is one of the two species of medaka perfectly adapted for living in brackish/sea-waters. Being an important component of the mangrove ecosystem, O. javanicus is also used as a valuable marine test-fish for ecotoxicology studies. Here, we sequenced and assembled the whole genome of O. javanicus, and anticipate this resource will be catalytic for a wide range of comparative genomic, phylogenetic and functional studies. Complementary sequencing approaches including long-read technology and data integration with a genetic map allowed the final assembly of 908 Mbp of the O. javanicus genome. Further analyses estimate that the O. javanicus genome contains 33% of repeat sequences and has a heterozygosity of 0.96%. The achieved draft assembly contains 525 scaffolds with a total length of 809.7 Mbp, a N50 of 6,3 Mbp and a L50 of 37 scaffolds. We identified 21454 predicted transcripts for a total transcriptome size of 57, 146, 583 bps. We provide here a high-quality chromosome scale draft genome assembly of the euryhaline Javafish medaka (321 scaffolds anchored on 24 chromosomes (representing 97.7% of the total bases)), and give emphasis on the evolutionary adaptation to salinity.

Detecting, Categorizing, and Correcting Coverage Anomalies of RNA-Seq Quantification

Because of incomplete reference transcriptomes, incomplete sequencing bias models, or other modeling defects, algorithms to infer isoform expression from RNA sequencing (RNA-seq) sometimes do not accurately model expression. We present a computational method to detect instances where a quantification algorithm could not completely explain the input reads. Our approach identifies regions where the read coverage significantly deviates from expectation. We call these regions "expression anomalies." We further present a method to attribute their cause to either the incompleteness of the reference transcriptome or algorithmic mistakes. We detect anomalies for 30 GEUVADIS and 16 Human Body Map samples. By correcting anomalies when possible, we reduce the number of falsely predicted instances of differential expression. Anomalies that cannot be corrected are suspected to indicate the existence of isoforms unannotated by the reference. We detected 88 common anomalies of this type and find that they tend to have a lower-than-expected coverage toward their 3' ends.

The Peptide Venom Composition of the Fierce Stinging Ant Tetraponera aethiops (Formicidae: Pseudomyrmecinae)

In the mutualisms involving certain pseudomyrmicine ants and different myrmecophytes (i.e., plants sheltering colonies of specialized “plant-ant” species in hollow structures), the ant venom contributes to the host plant biotic defenses by inducing the rapid paralysis of defoliating insects and causing intense pain to browsing mammals. Using integrated transcriptomic and proteomic approaches, we identified the venom peptidome of the plant-ant Tetraponera aethiops (Pseudomyrmecinae). The transcriptomic analysis of its venom glands revealed that 40% of the expressed contigs encoded only seven peptide precursors related to the ant venom peptides from the A-superfamily. Among the 12 peptide masses detected by liquid chromatography-mass spectrometry (LC–MS), nine mature peptide sequences were characterized and confirmed through proteomic analysis. These venom peptides, called pseudomyrmecitoxins (PSDTX), share amino acid sequence identities with myrmeciitoxins known for their dual offensive and defensive functions on both insects and mammals. Furthermore, we demonstrated through reduction/alkylation of the crude venom that four PSDTXs were homo- and heterodimeric. Thus, we provide the first insights into the defensive venom composition of the ant genus Tetraponera indicative of a streamlined peptidome.

Identification of potential chemosignals in the European water vole Arvicola terrestris

The water vole Arvicola terrestris is endemic to Europe where its outbreak generates severe economic losses for farmers. Our project aimed at characterising putative chemical signals used by this species, to develop new sustainable methods for population control that could also be used for this species protection in Great Britain. The water vole, as well as other rodents, uses specific urination sites as territorial and sex pheromone markers, still unidentified. Lateral scent glands and urine samples were collected from wild males and females caught in the field, at different periods of the year. Their volatile composition was analysed for each individual and not on pooled samples, revealing a specific profile of flank glands in October and a specific profile of urinary volatiles in July. The urinary protein content appeared more contrasted as males secrete higher levels of a lipocalin than females, whenever the trapping period. We named this protein arvicolin. Male and female liver transcript sequencing did not identify any expression of other odorant-binding protein sequence. This work demonstrates that even in absence of genome, identification of chemical signals from wild animals is possible and could be helpful in strategies of species control and protection.

Comparative Analysis of Strategies for De Novo Transcriptome Assembly in Prokaryotes: Streptomyces clavuligerus as a Case Study

The performance of software tools for de novo transcriptome assembly greatly depends on the selection of software parameters. Up to now, the development of de novo transcriptome assembly for prokaryotes has not been as remarkable as that for eukaryotes. In this contribution, Rockhopper2 was used to perform a comparative transcriptome analysis of Streptomyces clavuligerus exposed to diverse environmental conditions. The study focused on assessing the incidence of software parameters on software performance for the identification of differentially expressed genes as a final goal. For this, a statistical optimization was performed using the Transrate Assembly Score (TAS). TAS was also used for evaluating the software performance and for comparing it with related tools, e.g., Trinity. Transcriptome redundancy and completeness were also considered for this analysis. Rockhopper2 and Trinity reached a TAS value of 0.55092 and 0.58337, respectively. Trinity assembles transcriptomes with high redundancy, with 55.6% of transcripts having some duplicates. Additionally, we observed that the total number of differentially expressed genes (DEG) and their annotation greatly depends on the method used for removing redundancy and the tools used for transcript quantification. To our knowledge, this is the first work aimed at assessing de novo assembly software for prokaryotic organisms.

Environmental transcriptomes of invasive dreissena, a model species in ecotoxicology and invasion biology

Dreissenids are established model species for ecological and ecotoxicological studies, since they are sessile and filter feeder organisms and reflect in situ freshwater quality. Despite this strong interest for hydrosystem biomonitoring, omics data are still scarce. In the present study, we achieved full de novo assembly transcriptomes of digestive glands to gain insight into Dreissena polymorpha and D. rostriformis bugensis molecular knowledge. Transcriptomes were obtained by Illumina RNA sequencing of seventy-nine organisms issued from fifteen populations inhabiting sites that exhibits multiple freshwater contamination levels and different hydrosystem topographies (open or closed systems). Based on a recent de novo assembly algorithm, we carried out a complete, quality-checked and annotated transcriptomes. The power of the present study lies in the completeness of transcriptomes gathering multipopulational organisms sequencing and its full availability through an open access interface that gives a friendly and ready-to-use access to data. The use of such data for proteogenomic and targeted biological pathway investigations purpose is promising as they are first full transcriptomes for this two Dreissena species.

Grouper: graph-based clustering and annotation for improved de novo transcriptome analysis

Motivation

De novo transcriptome analysis using RNA-seq offers a promising means to study gene expression in non-model organisms. Yet, the difficulty of transcriptome assembly means that the contigs provided by the assembler often represent a fractured and incomplete view of the transcriptome, complicating downstream analysis. We introduce Grouper, a new method for clustering contigs from de novo assemblies that are likely to belong to the same transcripts and genes; these groups can subsequently be analyzed more robustly. When provided with access to the genome of a related organism, Grouper can transfer annotations to the de novo assembly, further improving the clustering.

Results

On de novo assemblies from four different species, we show that Grouper is able to accurately cluster a larger number of contigs than the existing state-of-the-art method. The Grouper pipeline is able to map greater than 10% more reads against the contigs, leading to accurate downstream differential expression analyses. The labeling module, in the presence of a closely related annotated genome, can efficiently transfer annotations to the contigs and use this information to further improve clustering. Overall, Grouper provides a complete and efficient pipeline for processing de novo transcriptomic assemblies.

Availability and implementation

The Grouper software is freely available at https://github.com/COMBINE-lab/grouper under the 2-clause BSD license.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comprehensive Influences of Overexpression of a MYB Transcriptor Regulating Anthocyanin Biosynthesis on Transcriptome and Metabolome of Tobacco Leaves

Overexpression of R2R3-MYB transcriptor can induce up-expression of anthocyanin biosynthesis structural genes, and improve the anthocyanin content in plant tissues, but it is not clear whether the MYB transcription factor overexpression does effect on other genes transcript and chemical compounds accumulation. In this manuscript, RNA-sequencing and the stepwise multiple ion monitoring-enhanced product ions (stepwise MIM-EPI) strategy were employed to evaluate the comprehensive effect of the MYB transcription factor LrAN2 in tobacco. Overexpression of LrAN2 could promote anthocyanin accumulation in a lot of tissues of tobacco cultivar Samsun. Only 185 unigenes express differently in a total of 160,965 unigenes in leaves, and 224 chemical compounds were differently accumulated. Three anthocyanins, apigeninidin chloride, pelargonidin 3-O-beta-D-glucoside and cyanidin 3,5-O-diglucoside, were detected only in transgenic lines, which could explain the phenotype of purple leaves. Except for anthocyanins, the phenylpropanoid, polyphenol (catechin), flavonoid, flavone and flavonol, belong to the same subgroups of flavonoids biosynthesis pathway with anthocyanin and were also up-accumulated. Overexpression of LrAN2 activated the bHLH (basic helix-loop-helix protein) transcription factor AN1b, relative to anthocyanin biosynthesis and the MYB transcription factor MYB3, relative to proanthocyanin biosynthesis. Then, the structural genes, relative to the phenylpropanoid pathway, were activated, which led to the up-accumulation of phenylpropanoid, polyphenol (catechin), flavonoid, flavone, flavonol and anthocyanin. The MYB transcription factor CPC, negative to anthocyanin biosynthesis, also induced up-expression in transgenic lines, which implied that a negative regulation mechanism existed in the anthocyanin biosynthesis pathway. The relative contents of all 19 differently accumulated amino and derivers were decreased in transgenic lines, which meant the phenylalanine biosynthesis pathway completed the same substrates with other amino acids. Interestingly, the acetylalkylglycerol acetylhydrolase was down-expressed in transgenic lines, which caused 19 lyso-phosphatidylcholine and derivatives of lipids to be up-accumulated, and 8 octodecane and derivatives were down-accumulated. This research will give more information about the function of MYB transcription factors on the anthocyanin biosynthesis and other chemical compounds and be of benefit to obtaining new plant cultivars with high anthocyanin content by biotechnology.