Molluscan phylogenomics requires strategically selected genomes

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

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

Mitogenomes do not substantially improve phylogenetic resolution in a young non-model adaptive radiation of freshwater gastropods

BACKGROUND

Species flocks in ancient lakes, and particularly those arising from adaptive radiation, make up the bulk of overall taxonomic and morphological diversity in these insular ecosystems. For these mostly young species assemblages, classical mitochondrial barcoding markers have so far been key to disentangle interspecific relationships. However, with the rise and further development of next-generation sequencing (NGS) methods and mapping tools, genome-wide data have become an increasingly important source of information even for non-model groups.

RESULTS

Here, we provide, for the first time, a comprehensive mitogenome dataset of freshwater gastropods endemic to Sulawesi and thus of an ancient lake invertebrate species flock in general. We applied low-coverage whole-genome sequencing for a total of 78 individuals including 27 out of the 28 Tylomelania morphospecies from the Malili lake system as well as selected representatives from Lake Poso and adjacent catchments. Our aim was to assess whether mitogenomes considerably contribute to the phylogenetic resolution within this young species flock. Interestingly, we identified a high number of variable and parsimony-informative sites across the other 'non-traditional' mitochondrial loci. However, although the overall support was very high, the topology obtained was largely congruent with previously published single-locus phylogenies. Several clades remained unresolved and a large number of species was recovered polyphyletic, indicative of both rapid diversification and mitochondrial introgression.

CONCLUSIONS

This once again illustrates that, despite the higher number of characters available, mitogenomes behave like a single locus and thus can only make a limited contribution to resolving species boundaries, particularly when introgression events are involved.

Life goes on: Spatial heterogeneity promotes biodiversity in an urbanized coastal marine ecosystem

Both human populations and marine biodiversity are concentrated along coastlines, with growing conservation interest in how these ecosystems can survive intense anthropogenic impacts. Tropical urban centres provide valuable research opportunities because these megacities are often adjacent to mega-diverse coral reef systems. The Pearl River Delta is a prime exemplar, as it encompasses one of the most densely populated and impacted regions in the world and is located just northwest of the Coral Triangle. However, the spatial and taxonomic complexity of this biodiversity, most of which is small, cryptic in habitat and poorly known, make comparative analyses challenging. We deployed standardized settlement structures at seven sites differing in the intensity of human impacts and used COI metabarcoding to characterize benthic biodiversity, with a focus on metazoans. We found a total of 7184 OTUs, with an average of 665 OTUs per sampling unit; these numbers exceed those observed in many previous studies using comparable methods, despite the location of our study in an urbanized environment. Beta diversity was also high, with 52% of the OTUs found at just one site. As expected, we found that the sites close to point sources of pollution had substantially lower diversity (44% less) relative to sites bathed in less polluted oceanic waters. However, the polluted sites contributed substantially to the total animal diversity of the region, with 25% of all OTUs occurring only within polluted sites. Further analysis of Arthropoda, Annelida and Mollusca showed that phylogenetic clustering within a site was common, suggesting that environmental filtering reduced biodiversity to a subset of lineages present within the region, a pattern that was most pronounced in polluted sites and for the Arthropoda. The water quality gradients surrounding the PRD highlight the unique role of in situ studies for understanding the impacts of complex urbanization pressures on biodiversity.

Diverse pathways to speciation revealed by marine snails

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

The Genomics in Emerging Marine Systems Checklist for Clear and Reproducible Genomics in Emerging, Marine Systems

Genome sequencing becomes more accessible and powerful every year, but there is a lack of consensus on what information should be provided in publications that include genomic data. The result is a flood of sequencing data without a framework to evaluate its quality and completeness, hindering reproducibility. In non-model taxa in marine systems, a lack of detail in methods sections often hinders future researchers from adopting improved techniques, leaving them to repeat costly protocols and take up computational (wall) time with programs that are already known to fail. Here, I present a set of guidelines tailored for marine taxa (emerging model organisms) to promote consistency between publications, increase transparency of sequencing projects, and preserve the value of sequence data as sequencing technologies advance. Included is a checklist to (1) guide authors toward including more detailed information in their manuscripts, (2) expand data availability, and (3) assist reviewers to thoroughly vet methods and results of future 'omic publications. This set of guidelines will support the usefulness of 'omic data in future analyses by providing a framework to document and evaluate these data, leading to transparent and reproducible genomics research on emerging marine systems.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

Chromosome-level genome assembly of the European flat oyster (Ostrea edulis) provides insights into its evolution and adaptation

The European flat oyster (Ostrea edulis) is an endangered and economically important marine bivalve species that plays a critical role in the coastal ecosystem. Here, we report a high-quality chromosome-level genome assembly of O. edulis, generated using PacBio HiFi-CCS long reads and annotated with Nanopore full-length transcriptome. The O. edulis genome covers 946.06 Mb (scaffold N50 94.82 Mb) containing 34,495 protein-coding genes and a high proportion of repeat sequences (58.49 %). The reconstructed demographic histories show that O. edulis population might be shaped by breeding habit (embryo brooding) and historical climatic change. Comparative genomic analysis indicates that transposable elements may drive lineage-specific evolution in oysters. Notably, the O. edulis genome has a Hox gene cluster rearrangement that has never been reported in bivalves, making this species valuable for evolutionary studies of molluscan diversification. Moreover, genome expansion of O. edulis is probably central to its adaptation to filter-feeding and sessile lifestyles, as well as embryo brooding and pathogen resistance, in coastal ecosystems. This chromosome-level genome assembly provides new insights into the genome feature of oysters, and presents an important resource for genetic research, evolutionary studies, and biological conservation of O. edulis.

A molecular phylogeny of the Persian Gulf and the Gulf of Oman oyster species

The taxonomy of oysters along the northern coasts of the Persian Gulf and the Gulf of Oman is not well recognized. We present a phylogenetic analysis of oyster species in these regions. We combined morphological and molecular techniques to obtain the identity of oysters to the lowest taxonomic levels. Analysis of partial nucleotide sequences from mitochondrial cytochrome c oxidase subunit I (COI) was used for the phylogenetic evaluation. Based on our findings, Iranian samples nested within the genus Saccostrea and belonged to Saccostrea mordax and Saccostrea palmula clades. The shell morphology of the studied samples was variable, as in other rock oyster species. The examination of morphological features was in line with the molecular outcomes, but despite some similarities, Iranian S. palmula had well-developed and elongated chomata. The results also showed that S. mordax and S. palmula possessed significant relative abundance as dominant oysters in the Persian Gulf and the Gulf of Oman, respectively. Phylogenetic analysis revealed that Iranian samples of S. palmula formed a separate subclade from the Gulf of California and Panama samples, with large genetic distances (6–7%). Iranian specimens differed morphologically and genetically, suggesting that they could be a new species, although more research is needed.

Status in molluscan cell line development in last one decade (2010–2020): impediments and way forward

Despite the attempts that have started since the 1960s, not even a single cell line of marine molluscs is available. Considering the vast contribution of marine bivalve aquaculture to the world economy, the prevailing viral threats, and the dismaying lack of advancements in molluscan virology, the requirement of a marine molluscan cell line is indispensable. This synthetic review discusses the obstacles in developing a marine molluscan cell line concerning the choice of species, the selection of tissue and decontamination, and cell culture media, with emphasis given on the current decade 2010-2020. Detailed accounts on the experiments on the virus cultivation in vitro and molluscan cell immortalization, with a brief note on the history and applications of the molluscan cell culture, are elucidated to give a holistic picture of the current status and future trends in molluscan cell line development.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10616-022-00539-x.

Evidence of multiple genome duplication events in Mytilus evolution

Background

Molluscs remain one significantly under-represented taxa amongst available genomic resources, despite being the second-largest animal phylum and the recent advances in genomes sequencing technologies and genome assembly techniques. With the present work, we want to contribute to the growing efforts by filling this gap, presenting a new high-quality reference genome for Mytilus edulis and investigating the evolutionary history within the Mytilidae family, in relation to other species in the class Bivalvia.

Results

Here we present, for the first time, the discovery of multiple whole genome duplication events in the Mytilidae family and, more generally, in the class Bivalvia. In addition, the calculation of evolution rates for three species of the Mytilinae subfamily sheds new light onto the taxa evolution and highlights key orthologs of interest for the study of Mytilus species divergences.

Conclusions

The reference genome presented here will enable the correct identification of molecular markers for evolutionary, population genetics, and conservation studies. Mytilidae have the capability to become a model shellfish for climate change adaptation using genome-enabled systems biology and multi-disciplinary studies of interactions between abiotic stressors, pathogen attacks, and aquaculture practises.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08575-9.

How many single-copy orthologous genes from whole genomes reveal deep gastropod relationships?

The Gastropoda contains 80% of existing mollusks and is the most diverse animal class second only to the Insecta. However, the deep phylogeny of gastropods has been controversial for a long time. Especially the position of Patellogastropoda is a major uncertainty. Morphology and some mitochondria studies concluded that Patellogastropoda is likely to be sister to all other gastropods (Orthogastropoda hypothesis), while transcriptomic and other mitogenomic studies indicated that Patellogastropoda and Vetigastropoda are sister taxa (Psilogastropoda). With the release of high-quality genomes, orthologous genes can be better identified and serve as powerful candidates for phylogenetic analysis. The question is, given the current limitations on the taxon sampling side, how many markers are needed to provide robust results. Here, we identified single-copy orthologous genes (SOGs) from 14 gastropods species with whole genomes available which cover five main gastropod subclasses. We generated different datasets from 395 to 1610 SOGs by allowing species missing in different levels. We constructed gene trees of each SOG, and inferred species trees from different collections of gene trees. We found as the number of SOGs increased, the inferred topology changed from Patellogastropoda being sister to all other gastropods to Patellogastropoda being sister to Vetigastropoda + Neomphalina (Psilogastropoda s.l.), with considerable support. Our study thus rejects the Orthogastropoda concept showing that the selection of the representative species and use of sufficient informative sites greatly influence the analysis of deep gastropod phylogeny.

Mobilizing molluscan models and genomes in biology

Molluscs are among the most ancient, diverse, and important of all animal taxa. Even so, no individual mollusc species has emerged as a broadly applied model system in biology. We here make the case that both perceptual and methodological barriers have played a role in the relative neglect of molluscs as research organisms. We then summarize the current application and potential of molluscs and their genomes to address important questions in animal biology, and the state of the field when it comes to the availability of resources such as genome assemblies, cell lines, and other key elements necessary to mobilising the development of molluscan model systems. We conclude by contending that a cohesive research community that works together to elevate multiple molluscan systems to 'model' status will create new opportunities in addressing basic and applied biological problems, including general features of animal evolution. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.