Intron Analyses Reveal Multiple Calmodulin Copies in Littorina

Spatial-temporal expression analysis of lineage-restricted shell matrix proteins reveals shell field regionalization and distinct cell populations in the slipper snail Crepidula atrasolea

Molluscs are one of the most morphologically diverse clades of metazoans, exhibiting an immense diversification of calcium carbonate structures, such as the shell. Biomineralization of the calcified shell is dependent on shell matrix proteins (SMPs). While SMP diversity is hypothesized to drive molluscan shell diversity, we are just starting to unravel SMP evolutionary history and biology. Here we leveraged two complementary model mollusc systems, Crepidula fornicata and Crepidula atrasolea , to determine the lineage-specificity of 185 Crepidula SMPs. We found that 95% of the adult C. fornicata shell proteome belongs to conserved metazoan and molluscan orthogroups, with molluscan-restricted orthogroups containing half of all SMPs in the shell proteome. The low number of C. fornicata -restricted SMPs contradicts the generally-held notion that an animal’s biomineralization toolkit is dominated by mostly novel genes. Next, we selected a subset of lineage-restricted SMPs for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval stages in C. atrasolea . We found that 12 out of 18 SMPs analyzed are expressed in the shell field. Notably, these genes are present in 5 expression patterns, which define at least three distinct cell populations within the shell field. These results represent the most comprehensive analysis of gastropod SMP evolutionary age and shell field expression patterns to date. Collectively, these data lay the foundation for future work to interrogate the molecular mechanisms and cell fate decisions underlying molluscan mantle specification and diversification.

Nuclear genome of Bulinus truncatus, an intermediate host of the carcinogenic human blood fluke Schistosoma haematobium

Some snails act as intermediate hosts (vectors) for parasitic flatworms (flukes) that cause neglected tropical diseases, such as schistosomiases. Schistosoma haematobium is a blood fluke that causes urogenital schistosomiasis and induces bladder cancer and increased risk of HIV infection. Understanding the molecular biology of the snail and its relationship with the parasite could guide development of an intervention approach that interrupts transmission. Here, we define the genome for a key intermediate host of S. haematobium-called Bulinus truncatus-and explore protein groups inferred to play an integral role in the snail's biology and its relationship with the schistosome parasite. Bu. truncatus shared many orthologous protein groups with Biomphalaria glabrata-the key snail vector for S. mansoni which causes hepatointestinal schistosomiasis in people. Conspicuous were expansions in signalling and membrane trafficking proteins, peptidases and their inhibitors as well as gene families linked to immune response regulation, such as a large repertoire of lectin-like molecules. This work provides a sound basis for further studies of snail-parasite interactions in the search for targets to block schistosomiasis transmission.

Bulinus truncatus transcriptome – a resource to enable molecular studies of snail and schistosome biology

Despite advances in high-throughput sequencing and bioinformatics, molecular investigations of snail intermediate hosts that transmit parasitic trematodes are scant. Here, we report the first transcriptome for Bulinus truncatus – a key intermediate host of Schistosoma haematobium – a blood fluke that causes urogenital schistosomiasis in humans. We assembled this transcriptome from short- and long-read RNA-sequence data. From this transcriptome, we predicted 12,998 proteins, 58% of which had orthologs in Biomphalaria glabrata – an intermediate host of Schistosoma mansoni – a blood fluke that causes hepato-intestinal schistosomiasis. We predicted that select protein groups are involved in signal transduction, cell growth and death, the immune system, environmental adaptation and/or the excretory/secretory system, suggesting roles in immune responses, pathogen defence and/or parasite-host interactions. The transcriptome of Bu. truncatus provides a useful resource to underpin future molecular investigations of this and related snail species, and its interactions with pathogens including S. haematobium. The present resource should enable comparative investigations of other molluscan hosts of socioeconomically important parasites in the future.

•

First transcriptome to represent Bulinus truncatus – a snail intermediate host of Schistosoma haematobium.

•

Select protein groups of Bu. truncatus are inferred to associate with innate immune responses against pathogens.

•

Transcriptome provides a resource for future studies of parasite-host interactions and snail-host resistance to pathogens.

Multiple calmodulin genes of the Pacific abalone, Haliotis discus hannai (Mollusca: Vetigastropoda: Haliotidae)

In this study, we identified four canonical calmodulin genes in the Pacific abalone, Haliotis discus hannai. Their full-length cDNAs were variable in the 5' and 3' untranslated regions, but highly similar (91-97%) in the coding region. Each of the genes encoded 149 amino acids, with 93-97% similarity among themselves and 94-98% similarity with human CAM I. There were 54 substitutions distributed unevenly throughout the coding regions, found mostly in the third codon position. Gene structure analysis revealed that each of the calmodulin genes comprised five exons and four introns. The intron positions and phases were identical and there were no introns in the fourth exon. The corresponding introns differed in their sequences and sizes. Expression profiles of nine tissues from abalone revealed that the calmodulin genes were transcribed in common in gill and mantle tissue, but differentially in the other tissues. A phylogenetic analysis based on the amino acid sequences revealed that calmodulin C was the most common isoform in Gastropoda and calmodulin was the most diverged isoform. An in silico analysis of the calmodulin genes identified paralogous genes in other Haliotis species, indicating that gene duplication might have occurred in the last common ancestor of Haliotis. Abbreviations: ORF: open reading frame; RACE: random amplification of cDNA end; TSA: transcriptome shotgun assembly; UTR: untranslated region.

Phylogenetic analysis of the Lancinae (Gastropoda, Lymnaeidae) with a description of the U.S. federally endangered Banbury Springs lanx

We examined the patelliform snails of the subfamily Lancinae, endemic to northwestern North America, to test whether morphological variation correlated with genetic and anatomical differences. Molecular analyses using cox1, 16S, calmodulin intron, and 28S rDNA partial sequences and anatomical data supported recognition of four species in three genera. The relationships of lancines within Lymnaeidae are not yet well-resolved. The federally endangered Banbury Springs lanx is described as a new genus and species, Idaholanxfresti, confirming its distinctiveness and narrow endemicity.

Conserved properties of individual Ca2+-binding sites in calmodulin

Calmodulin (CaM) is a Ca(2+)-sensing protein that is highly conserved and ubiquitous in eukaryotes. In humans it is a locus of life-threatening cardiomyopathies. The primary function of CaM is to transduce Ca(2+) concentration into cellular signals by binding to a wide range of target proteins in a Ca(2+)-dependent manner. We do not fully understand how CaM performs its role as a high-fidelity signal transducer for more than 300 target proteins, but diversity among its four Ca(2+)-binding sites, called EF-hands, may contribute to CaM's functional versatility. We therefore looked at the conservation of CaM sequences over deep evolutionary time, focusing primarily on the four EF-hand motifs. Expanding on previous work, we found that CaM evolves slowly but that its evolutionary rate is substantially faster in fungi. We also found that the four EF-hands have distinguishing biophysical and structural properties that span eukaryotes. These results suggest that all eukaryotes require CaM to decode Ca(2+) signals using four specialized EF-hands, each with specific, conserved traits. In addition, we provide an extensive map of sites associated with target proteins and with human disease and correlate these with evolutionary sequence diversity. Our comprehensive evolutionary analysis provides a basis for understanding the sequence space associated with CaM function and should help guide future work on the relationship between structure, function, and disease.

Characterization of cDNAs for calmodulin and calmodulin-like protein in the freshwater mussel Hyriopsis cumingii: differential expression in response to environmental Ca(2+) and calcium binding of recombinant proteins

Calmodulin and calmodulin-like protein are two crucial calcium regulators in bivalves. However, molecular characteristics and expression patterns of these genes in the freshwater mussel are poorly understood. In this study, two cDNAs encoding novel calmodulin and calmodulin-like protein (HcCaM and HcCaLP) were cloned and characterized from the freshwater pearl mussel Hyriopsis cumingii. The full-length cDNA of HcCaM was 726 bp, including a 118-bp 5'-untranslated region (UTR), a 447-bp open reading frame (ORF), and a 161-bp 3'-UTR. The 1217-bp HcCaLP cDNA comprised of a 51-bp 5'-UTR, a 447-bp ORF, and a 716-bp 3'-UTR. The potential phosphorylation sites of, Arg(80) and Phe(100) in deduced HcCaM were mutated to Thr(80) and Tyr(100) in HcCaLP. Tissue-specific expression analysis revealed that HcCaM mRNA was prominently expressed in the gill, mantle center, and foot. In contrast, HcCaLP mRNA was mainly expressed in the mantle edge. The recombinant HcCaM and HcCaLP proteins expressed in Escherichia coli showed the typical Ca(2+) dependent electrophoretic shift characterization as CaM and differed in the calcium binding affinity. The calcium stimulation test that lasted 5 weeks implied that HcCaM and HcCaLP had differential expression patterns in response to various environmental Ca(2+) concentrations (0.25-1.25 mM). The expression of HcCaM mRNA was up-regulated by low Ca(2+) concentration (0.25 mM), and the highest expression of HcCaLP mRNA occurred under Ca(2+) concentration of 1 mM.

DNA barcoding methods for invertebrates

Invertebrates comprise approximately 34 phyla, while vertebrates represent one subphylum and insects a (very large) class. Thus, the clades excepting vertebrates and insects encompass almost all of animal diversity. Consequently, the barcoding challenge in invertebrates is that of barcoding animals in general. While standard extraction, cleaning, PCR methods, and universal primers work for many taxa, taxon-specific challenges arise because of the shear genetic and biochemical diversity present across the kingdom, and because problems arising as a result of this diversity, and solutions to them, are still poorly characterized for many metazoan clades. The objective of this chapter is to emphasize general approaches, and give practical advice for overcoming the diverse challenges that may be encountered across animal taxa, but we stop short of providing an exhaustive inventory. Rather, we encourage researchers, especially those working on poorly studied taxa, to carefully consider methodological issues presented below, when standard approaches perform poorly.

Specific versus non-specific immune responses in an invertebrate species evidenced by a comparative de novo sequencing study

Our present understanding of the functioning and evolutionary history of invertebrate innate immunity derives mostly from studies on a few model species belonging to ecdysozoa. In particular, the characterization of signaling pathways dedicated to specific responses towards fungi and Gram-positive or Gram-negative bacteria in Drosophila melanogaster challenged our original view of a non-specific immunity in invertebrates. However, much remains to be elucidated from lophotrochozoan species. To investigate the global specificity of the immune response in the fresh-water snail Biomphalaria glabrata, we used massive Illumina sequencing of 5′-end cDNAs to compare expression profiles after challenge by Gram-positive or Gram-negative bacteria or after a yeast challenge. 5′-end cDNA sequencing of the libraries yielded over 12 millions high quality reads. To link these short reads to expressed genes, we prepared a reference transcriptomic database through automatic assembly and annotation of the 758,510 redundant sequences (ESTs, mRNAs) of B. glabrata available in public databases. Computational analysis of Illumina reads followed by multivariate analyses allowed identification of 1685 candidate transcripts differentially expressed after an immune challenge, with a two fold ratio between transcripts showing a challenge-specific expression versus a lower or non-specific differential expression. Differential expression has been validated using quantitative PCR for a subset of randomly selected candidates. Predicted functions of annotated candidates (approx. 700 unisequences) belonged to a large extend to similar functional categories or protein types. This work significantly expands upon previous gene discovery and expression studies on B. glabrata and suggests that responses to various pathogens may involve similar immune processes or signaling pathways but different genes belonging to multigenic families. These results raise the question of the importance of gene duplication and acquisition of paralog functional diversity in the evolution of specific invertebrate immune responses.

Exon-primed, intron-crossing (EPIC) loci for five nuclear genes in deep-sea protobranch bivalves: primer design, PCR protocols and locus utility

We describe PCR primers and amplification protocols developed to obtain introns from conserved nuclear genes in deep-sea protobranch bivalves. Because almost no sequence data for protobranchs are publically available, mollusk and other protostome sequences from GenBank were used to design degenerate primers, making these loci potentially useful in other invertebrate taxa. Amplification and sequencing success varied across the test group of 30 species, and we present five loci spanning this range of outcomes. Intron presence in the targeted regions also varied across genes and species, often within single genera; for instance, the calmodulin and β-tubulin loci contained introns with high frequency, whereas the triose phosphate isomerase locus never contained an intron. In introns for which we were able to obtain preliminary estimates of polymorphism levels in single species, polymorphism was greater than traditional mitochondrial loci. These markers will greatly increase the ability to assess population structure in the ecologically important protobranchs, and may prove useful in other taxa as well.

Analysis of Acropora muricata calmodulin (CaM) indicates that scleractinian corals possess the ancestral exon/intron organization of the eumetazoan CaM gene

Calmodulin (CaM), belonging to the tropinin C (TnC) superfamily, is one of the calcium-binding proteins that are highly conserved in their protein and gene structure. Based on the structure comparison among published vertebrate and invertebrate CaM, it is proposed that the ancestral form of eumetazoan CaM genes should have five exons and four introns (four-intron hypothesis). In this study, we determined the gene structure of CaM in the coral Acropora muricata, an anthozoan cnidarian representing the basal position in animal evolution. A CaM clone was isolated from a cDNA library constructed from the spawned eggs of A. muricata. This clone was composed of 908 nucleotides, including 162 base pairs (bp) of 5'-untranslated region (UTR), 296 bp of 3'-UTR, and an open reading frame 450 bp in length. The deduced amino acid indicated that the Acropora CaM protein is identical to that of the actiniarian, Metridinium senile, and has four putative calcium-binding domains highly similar to those of other vertebrate or invertebrate CaMs. Southern blot analysis revealed that Acropora CaM is a putative single-copy gene in the nuclear genome. Genomic sequencing showed that Acropora CaM was composed of five exons and four introns, with intron II not corresponding to any region in the actiniarian CaM gene, which possesses only four exons and three introns. Our results highlight that the coral CaM gene isolated from A. muricata has four introns at the predicted positions of the early metazoan CaM gene organization, providing the first evidence from the basal eumetazoan phylum to support the four-intron hypothesis.

Patterns of genetic structure among Hawaiian corals of the genus Pocillopora yield clusters of individuals that are compatible with morphology

Six variable sequence markers are developed and analyzed to find out species boundaries in Hawaiian corals of the genus Pocillopora: the putative mitochondrial control region; a recently discovered, hypervariable mitochondrial open reading frame; the internal transcribed spacer 2 (ITS2), located in the nuclear ribosomal DNA; three nuclear introns of calmodulin, elongation factor-1alpha and the ATP synthase beta subunit. Using the first two markers, we identify five distinct mitochondrial lineages and these lineages are compatible with morphology. The situation is more complex with nuclear markers since more than two haplotypes are observed in some individuals. To detect clusters of individuals, haplotype networks are constructed with additional connections drawn between co-occurring haplotypes to delineate potential fields for recombination: few clusters of nuclear haplotypes are found to correspond to clusters of individuals, but those that are detected (mostly in the ITS2 dataset) are also compatible with morphology.