Comparative Transcriptomics of Alternative Developmental Phenotypes in a Marine Gastropod: TRANSCRIPTOMICS OF NUTRITIVE EMBRYOS

Deep conservation and co-option of programmed cell death facilitates evolution of alternative phenotypes at multiple biological levels

Alternative phenotypes, such as polyphenisms and sexual dimorphisms, are widespread in nature and appear at all levels of biological organization, from genes and cells to morphology and behavior. Yet, our understanding of the mechanisms through which alternative phenotypes develop and how they evolve remains understudied. In this review, we explore the association between alternative phenotypes and programmed cell death, a mechanism responsible for the elimination of superfluous cells during development. We discuss the ancient origins and deep conservation of programmed cell death (its function, forms and underlying core regulatory gene networks), and propose that it was co-opted repeatedly to generate alternative phenotypes at the level of cells, tissues, organs, external morphology, and even individuals. We review several examples from across the tree of life to explore the conditions under which programmed cell death is likely to facilitate the evolution of alternative phenotypes.

Slipper snail tales: How Crepidula fornicata and Crepidula atrasolea became model molluscs

Despite the great abundance and diversity of molluscs, only a few have attained "model research organism" status. One of those species is the slipper snail Crepidula fornicata. Its embryos were first used for classical lineage tracing studies in the late 19th century, and over a 100 years later they were "re-discovered" by our labs and used for modern fate mapping, gene perturbation, in vivo imaging, transcriptomics, and the first application of CRISPR/Cas9-mediated genome editing among the Spiralia/Lophotrochozoa. Simultaneously, other labs made extensive examinations of taxonomy, phylogeny, ecology, life-history, mode of development, larval feeding behavior, and responses to the environment in members of the family Calyptraeidae, which includes the genus Crepidula. Recently, we developed tools, resources, and husbandry protocols for another, direct-developing species, Crepidula atrasolea. This species is an ideal "lab rat" among molluscs. Together these species will be valuable for probing the cellular and molecular mechanisms underlying molluscan biology and evolution.

Population transcriptomics reveals the effect of gene flow on the evolution of range limits

One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.

The transcriptome analysis of the whole-body of the gastropod mollusk Limax flavus and screening of putative antimicrobial peptide and protein genes

The gastropod mollusk Limax flavus, one of the most widespread pests in China, is used to treat infectious diseases in traditional Chinese medicine. However, little genomic information is available for this non-model species. In this study, the whole-body transcriptome of L. flavus was sequenced using next generation sequencing technology. A total of 6.81 Gb clean reads were obtained, which were assembled into 150,766 transcripts with 132,206 annotated unigenes. Functionally classification assigned 30,542 unigenes to 56 Gene Ontology terms, 16,745 unigenes were divided into 26 euKaryotic Ortholog Groups of proteins categories, and 13,854 unigenes were assigned to 230 Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, we identified 17,251 simple sequence repeats and several kinds of antimicrobial peptide and protein (AMPs) genes. The transcriptome data of L. flavus will provide a valuable genomic resource for further studies on this species, and the AMPs identified in L. flavus will support its medical potential.

De novo transcriptome of the whole-body of the gastropod mollusk Philomycus bilineatus, a pest with medical potential in China

Philomycus bilineatus is a highly common gastropod mollusk pest in China and is also utilized to treat infectious diseases. However, no genomic resources are available for this non-model species. In the present study, the transcriptomic analysis of P. bilineatus was completed. After sequencing using the next generation sequencing technology, 9.11 Gb of clean reads were obtained, which led to the assembly and annotation of 145,523 transcripts and 125,690 unigenes. Unigenes were functionally classified using Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG). A total of 27,554 unigenes were assigned into 55 GO terms, 13,989 unigenes were differentiated into 26 KOG categories, and 16,368 unigenes were assigned to 229 KEGG pathways. Furthermore, 16,614 simple sequence repeats (SSRs), 38 olfactory genes, and 40 antimicrobial peptide/protein genes were identified. The transcriptome profile of P. bilineatus will provide a valuable genomic resource for further study, will promote the development of new pest management strategies through interference of chemosensory communication, and will support potential medicinal uses of this species.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Mollusc genomes reveal variability in patterns of LTR-retrotransposons dynamics

Background

The three superfamilies of Long Terminal Repeat (LTR) retrotransposons are a widespread kind of transposable element and a major factor in eukaryotic genome evolution. In metazoans, recent studies suggested that Copia LTR-retrotransposons display specific dynamic compared to the more abundant and diverse Gypsy elements. Indeed, Copia elements show a relative scarcity and the prevalence of only a few clades in specific hosts. Thus, BEL/Pao seems to be the second most abundant superfamily. However, the generality of these assumptions remains to be assessed. Therefore, we carried out the first large-scale comparative genomic analysis of LTR-retrotransposons in molluscs. The aim of this study was to analyse the diversity, copy numbers, genomic proportions and distribution of LTR-retrotransposons in a large host phylum.

Results

We compare nine genomes of molluscs and further added LTR-retrotransposons sequences detected in databases for 47 additional species. We identified 1709 families, which enabled us to define 31 clades. We show that clade richness was highly dependent on the considered superfamily. We found only three Copia clades, including GalEa and Hydra which appear to be widely distributed and highly dominant as they account for 96% of the characterised Copia elements. Among the seven BEL/Pao clades identified, Sparrow and Surcouf are characterised for the first time. We find no BEL or Pao elements, but the rare clades Dan and Flow are present in molluscs. Finally, we characterised 21 Gypsy clades, only five of which had been previously described, the C-clade being the most abundant one. Even if they are found in the same number of host species, Copia elements are clearly less abundant than BEL/Pao elements in copy number or genomic proportions, while Gypsy elements are always the most abundant ones whatever the parameter considered.

Conclusions

Our analysis confirms the contrasting dynamics of Copia and Gypsy elements in metazoans and indicates that BEL/Pao represents the second most abundant superfamily, probably reflecting an intermediate dynamic. Altogether, the data obtained in several taxa highly suggest that these patterns can be generalised for most metazoans. Finally, we highlight the importance of using database information in complement of genome analyses when analyzing transposable element diversity.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5200-1) contains supplementary material, which is available to authorized users.

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.

Decoding the architecture and origins of mechanisms for developmental polyphenism

Developmental polyphenism affords a single genotype multiple solutions to match an organism to its environment. Because polyphenism is the extreme example of how development deviates from a linear genetic blueprint, it demands a genetic explanation for how environmental cues shunt development to hypothetically alternative modules. We highlight several recent advances that have begun to illuminate genetic mechanisms for polyphenism and how this recurring developmental novelty may arise. An emerging genetic knowledge of polyphenism is providing precise targets for testing hypotheses of how switch mechanisms are built-out of olfactory, nutrient-sensing, hormone-reception, and developmental and genetic buffering systems-to accommodate plasticity. Moreover, classic and new model systems are testing the genetic basis of polyphenism's proposed causal roles in evolutionary change.

Molluscan models: Crepidula fornicata

Gastropod snails in the genus Crepidula have emerged as model systems for studying a metazoan super clade, the Spiralia. Recent work on one species in particular, Crepidula fornicata, has produced high-resolution cell lineage fate maps, details of morphogenetic events during gastrulation, key insights into the molecular underpinnings of early development, and the first demonstration of CRISPR/Cas9 genome editing in the Spiralia. Furthermore, invasive species of Crepidula are a significant ecological threat, while one of these, C. fornicata, is also being harvested for food. This review highlights progress towards developing these animals as models for evolutionary, developmental, and ecological studies. Such studies have contributed greatly to our understanding of biology in a major clade of bilaterians. This information may also help us to control and cultivate these snails.