Proteomic comparison of the organic matrices from parietal and base plates of the acorn barnacle Amphibalanus amphitrite

Acorn barnacles are efficient colonizers on a wide variety of marine surfaces. As they proliferate on critical infrastructure, their settlement and growth have deleterious effects on performance. To address acorn barnacle biofouling, research has focused on the settlement and adhesion processes with the goal of informing the development of novel coatings. This effort has resulted in the discovery and characterization of several proteins found at the adhesive substrate interface, i.e. cement proteins, and a deepened understanding of the function and composition of the biomaterials within this region. While the adhesive properties at the interface are affected by the interaction between the proteins, substrate and mechanics of the calcified base plate, little attention has been given to the interaction between the proteins and the cuticular material present at the substrate interface. Here, the proteome of the organic matrix isolated from the base plate of the acorn barnacle Amphibalanus amphitrite is compared with the chitinous and proteinaceous matrix embedded within A. amphitrite parietal plates. The objective was to gain an understanding of how the basal organic matrix may be specialized for adhesion via an in-depth comparative proteome analysis. In general, the majority of proteins identified in the parietal matrix were also found in the basal organic matrix, including nearly all those grouped in classes of cement proteins, enzymes and pheromones. However, the parietal organic matrix was enriched with cuticle-associated proteins, of which ca 30% of those identified were unique to the parietal region. In contrast, ca 30-40% of the protease inhibitors, enzymes and pheromones identified in the basal organic matrix were unique to this region. Not unexpectedly, nearly 50% of the cement proteins identified in the basal region were significantly distinct from those found in the parietal region. The wider variety of identified proteins in the basal organic matrix indicates a greater diversity of biological function in the vicinity of the substrate interface where several processes related to adhesion, cuticle formation and expansion of the base synchronize to play a key role in organism survival.

Comparative genomics reveals the dynamic evolutionary history of cement protein genes of barnacles from intertidal to deep-sea hydrothermal vents

Thoracican barnacles are a diverse group of marine organisms for which the availability of genome assemblies is currently limited. In this study, we sequenced the genomes of two neolepadoid species (Ashinkailepas kermadecensis, Imbricaverruca yamaguchii) from hydrothermal vents, in addition to two intertidal species. Genome sizes ranged from 481 to 1054 Mb, with repetitive sequence contents of 21.2% to 50.7%. Concordance rates of orthologs and heterozygosity rates were between 82.4% and 91.7% and between 1.0% and 2.1%, respectively, indicating high genetic diversity and heterozygosity. Based on phylogenomic analyses, we revised the nomenclature of cement genes encoding cement proteins that are not homologous to any known proteins. The major cement gene, CP100A, was found in all thoracican species, including vent-associated neolepadoids, and was hypothesised to be essential for thoracican settlement. Duplicated genes, CP100B and CP100C, were found only in balanids, suggesting potential functional redundancy or acquisition of new functions associated with the calcareous base. An ancestor of CP52 genes was duplicated dynamically among lepadids, pollicipedids with multiple copies on a single scaffold, and balanids with multiple sequential repeats of the conserved regions, but no CP52 genes were found in neolepadoids, providing insights into cement gene evolution among thoracican lineages. This study enhances our understanding of the adhesion mechanisms of thoracicans in underwater environments. The newly sequenced genomes provide opportunities for studying their evolution and ecology, shedding light on their adaptation to diverse marine environments, and contributing to our knowledge of barnacle biology with valuable genomic resources for further studies in this field.

An ancient whole-genome duplication in barnacles contributes to their diversification and intertidal sessile life adaptation

INTRODUCTION

Whole-genome duplication (WGD) is one of the most sudden and dramatic events rarely reported in invertebrates, but its occurrence can lead to physiological, morphological, and behavioral diversification. WGD has also never been reported in barnacles, which is one of the most unique groups of crustaceans with extremely speciallized morphology (calcareous shells) and habits (intertidal sessile lifestyle).

OBJECTIVES

To investigate whether WGD has occurred in barnacles and examine its potential role in driving the adaptive evolution and diversification of barnacles.

METHODS

Based on a newly sequenced and assembled chromosome-level barnacle genome, a novel WGD event has been identified in barnacles through a comprehensive analysis of interchromosomal synteny, the Hox gene cluster, and synonymous substitution distribution.

RESULTS

We provide ample evidences for WGD in the barnacle genomes. Comparative genomic analysis indicates that this WGD event predates the divergence of Thoracicalcarea, occurring more than 247 million years ago. The retained ohnologs from the WGD are primarily enriched in various pathways related to environmental information processing, shedding light on the adaptive evolution and diversification of intertidal sessile lifestyle. In addition, transcriptomic analyses show that most of these ohnologs were differentially expressed following the ebb of tide. And the cytochrome P450 ohnologs with differential expression patterns are subject to subfunctionalization and/or neofunctionalization for intertidal adaptation. Besides WGD, parallel evolution underlying intertidal adaptation has also occurred in barnacles.

CONCLUSION

This study revealed an ancient WGD event in the barnacle genomes, which is potentially associated with the origin and diversification of thoracican barnacles, and may have contributed to the adaptive evolution of their intertidal sessile lifestyle.

Gene co-option, duplication and divergence of cement proteins underpin the evolution of bioadhesives across barnacle life histories

Acquisition of new genes often results in the emergence of novel functions and is a key step in lineage-specific adaptation. As a group of sessile crustaceans, barnacles establish permanent attachment through initial cement secretion at the larval phase followed by continuous cement secretion in juveniles and adults. However, the origins and evolution of barnacle larval and adult cement proteins remain poorly understood. By performing microdissection of larval cement glands, transcriptome and shotgun proteomics and immunohistochemistry validation, we identified 30 larval and 27 adult cement proteins of the epibiotic turtle barnacle Chelonibia testudinaria, of which the majority are stage- and barnacle-specific. While only two proteins, SIPC and CP100K, were expressed in both larvae and adults, detection of protease inhibitors and the cross-linking enzyme lysyl oxidase paralogs in larvae and adult cement. Other barnacle-specific cement proteins such as CP100k and CP52k likely share a common origin dating back at least to the divergence of Rhizocephala and Thoracica. Different CP52k paralogues could be detected in larval and adult cement, suggesting stage-specific cement proteins may arise from duplication followed by changes in expression timing of the duplicates. Interestingly, the biochemical properties of larval- and adult-specific CP52k paralogues exhibited remarkable differences. We conclude that barnacle larval and adult cement systems evolved independently, and both emerged from co-option of existing genes and de novo formation, duplication and functional divergence of lineage-specific cement protein genes. Our findings provide important insights into the evolutionary mechanisms of bioadhesives in sessile marine invertebrates.

Abundance and Diversification of Repetitive Elements in Decapoda Genomes

Repetitive elements are a major component of DNA sequences due to their ability to propagate through the genome. Characterization of Metazoan repetitive profiles is improving; however, current pipelines fail to identify a significant proportion of divergent repeats in non-model organisms. The Decapoda order, for which repeat content analyses are largely lacking, is characterized by extremely variable genome sizes that suggest an important presence of repetitive elements. Here, we developed a new standardized pipeline to annotate repetitive elements in non-model organisms, which we applied to twenty Decapoda and six other Crustacea genomes. Using this new tool, we identified 10% more repetitive elements than standard pipelines. Repetitive elements were more abundant in Decapoda species than in other Crustacea, with a very large number of highly repeated satellite DNA families. Moreover, we demonstrated a high correlation between assembly size and transposable elements and different repeat dynamics between Dendrobranchiata and Reptantia. The patterns of repetitive elements largely reflect the phylogenetic relationships of Decapoda and the distinct evolutionary trajectories within Crustacea. In summary, our results highlight the impact of repetitive elements on genome evolution in Decapoda and the value of our novel annotation pipeline, which will provide a baseline for future comparative analyses.

Chromosome-scale genome assembly of the freshwater cladoceran crustacean Chydorus sphaericus: A resource for discovery of genes responsive to ecological challenges

The genus Chydorus Leach 1816 (family Chydoridae) is a diverse and ecologically important taxon within freshwater ecosystems. Despite having been widely used in ecological, evolutionary and eco-toxicological studies, no high-quality genomic resource is available for any member of the genus. Here, we present a high-quality chromosome-level assembly of the C. sphaericus genome by combining 7.40 Gb (∼ 50 × coverage) PacBio reads, 19.28 Gb (∼ 135 × coverage) Illumina paired-end reads, and 34.04 Gb Hi-C reads. Our genome assembly is approximately 151 Mb, with contig and scaffold N50 lengths of 1.09 Mb and 13.70 Mb, respectively. The assembly captured 94.9% of the complete eukaryotic BUSCO. Repetitive elements accounted for 17.6% of the genome, and 13,549 protein-coding genes were predicted (based on transcriptome sequencing data, ab-initio or homology-based prediction), of which 96.4% have been functionally annotated in the NCBI-NR database. We identified 303 gene families specific to C. sphaericus, mainly families enriched in functions related to immune response, visual senses and detoxification. Interestingly, we also found 53 significantly expanded gene families in C. sphaericus, mostly with functions related to detoxification. This high-quality assembly genome will act as a reference genome for C. sphaericus and benefit studies on functional and comparative genomics of Chydorus and other crustaceans.

Genomic Adaptations to an Endoparasitic Lifestyle in the Morphologically Atypical Crustacean Sacculina carcini (Cirripedia: Rhizocephala)

The endoparasitic crustacean Sacculina carcini (Cirripedia: Rhizocephala) has a much simpler morphology than conventional filter-feeding barnacles, reflecting its parasitic lifestyle. To investigate the molecular basis of its refined developmental program, we produced a draft genome sequence for comparison with the genomes of nonparasitic barnacles and characterized the transcriptomes of internal and external tissues. The comparison of clusters of orthologous genes revealed the depletion of multiple gene families but also several unanticipated expansions compared to non-parasitic crustaceans. Transcriptomic analyses comparing interna and externa tissues revealed an unexpected variation of gene expression between rootlets sampled around host midgut and thoracic ganglia. Genes associated with lipid uptake were strongly expressed by the internal tissues. We identified candidate genes probably involved in host manipulation (suppression of ecdysis and gonad development) including those encoding crustacean neurohormones and the juvenile hormone binding protein. The evolution of Rhizocephala therefore appears to have involved a rapid turnover of genes (losses and expansions) as well as the fine tuning of gene expression.

Convergent evolution of barnacles and molluscs sheds lights in origin and diversification of calcareous shell and sessile lifestyle

The calcareous shell and sessile lifestyle are the representative phenotypes of many molluscs, which happen to be present in barnacles, a group of unique crustaceans. The origin of these phenotypes is unclear, but it may be embodied in the convergent genetics of such distant groups (interphylum). Herein, we perform comprehensive comparative genomics analysis in barnacles and molluscs, and reveal a genome-wide strong convergent molecular evolution between them, including coexpansion of biomineralization and organic matrix genes for shell formation, and origination of lineage-specific orphan genes for settlement. Notably, the expanded biomineralization gene encoding alkaline phosphatase evolves a novel, highly conserved motif that may trigger the origin of barnacle shell formation. Unlike molluscs, barnacles adopt novel organic matrices and cement proteins for shell formation and settlement, respectively, and their calcareous shells have potentially originated from the cuticle system of crustaceans. Therefore, our study corroborates the idea that selection pressures driving convergent evolution may strongly act in organisms inhabiting similar environments regardless of phylogenetic distance. The convergence signatures shed light on the origin of the shell and sessile lifestyle of barnacles and molluscs. In addition, notable non-convergence signatures are also present and may contribute to morphological and functional specificities.

Chromosome-level genome assembly, annotation, and phylogenomics of the gooseneck barnacle Pollicipes pollicipes

BACKGROUND

The barnacles are a group of >2,000 species that have fascinated biologists, including Darwin, for centuries. Their lifestyles are extremely diverse, from free-swimming larvae to sessile adults, and even root-like endoparasites. Barnacles also cause hundreds of millions of dollars of losses annually due to biofouling. However, genomic resources for crustaceans, and barnacles in particular, are lacking.

RESULTS

Using 62× Pacific Biosciences coverage, 189× Illumina whole-genome sequencing coverage, 203× HiC coverage, and 69× CHi-C coverage, we produced a chromosome-level genome assembly of the gooseneck barnacle Pollicipes pollicipes. The P. pollicipes genome is 770 Mb long and its assembly is one of the most contiguous and complete crustacean genomes available, with a scaffold N50 of 47 Mb and 90.5% of the BUSCO Arthropoda gene set. Using the genome annotation produced here along with transcriptomes of 13 other barnacle species, we completed phylogenomic analyses on a nearly 2 million amino acid alignment. Contrary to previous studies, our phylogenies suggest that the Pollicipedomorpha is monophyletic and sister to the Balanomorpha, which alters our understanding of barnacle larval evolution and suggests homoplasy in a number of naupliar characters. We also compared transcriptomes of P. pollicipes nauplius larvae and adults and found that nearly one-half of the genes in the genome are differentially expressed, highlighting the vastly different transcriptomes of larvae and adult gooseneck barnacles. Annotation of the genes with KEGG and GO terms reveals that these stages exhibit many differences including cuticle binding, chitin binding, microtubule motor activity, and membrane adhesion.

CONCLUSION

This study provides high-quality genomic resources for a key group of crustaceans. This is especially valuable given the roles P. pollicipes plays in European fisheries, as a sentinel species for coastal ecosystems, and as a model for studying barnacle adhesion as well as its key position in the barnacle tree of life. A combination of genomic, phylogenetic, and transcriptomic analyses here provides valuable insights into the evolution and development of barnacles.

Histology and transcriptomic analyses of barnacles with different base materials and habitats shed lights on the duplication and chemical diversification of barnacle cement proteins

Background

Barnacles are sessile crustaceans that attach to underwater surfaces using barnacle cement proteins. Barnacles have a calcareous or chitinous membranous base, and their substratum varies from biotic (e.g. corals/sponges) to abiotic surfaces. In this study, we tested the hypothesis that the cement protein (CP) composition and chemical properties of different species vary according to the attachment substrate and/or the basal structure. We examined the histological structure of cement glands and explored the variations in cement protein homologs of 12 barnacle species with different attachment habitats and base materials.

Results

Cement gland cells in the rocky shore barnacles Tetraclita japonica formosana and Amphibalanus amphitrite are eosinophilic, while others are basophilic. Transcriptome analyses recovered CP homologs from all species except the scleractinian coral barnacle Galkinia sp. A phylogenomic analysis based on sequences of CP homologs did not reflect a clear phylogenetic pattern in attachment substrates. In some species, certain CPs have a remarkable number of paralogous sequences, suggesting that major duplication events occurred in CP genes. The examined CPs across taxa show consistent bias toward particular sets of amino acid. However, the predicted isoelectric point (pI) and hydropathy are highly divergent. In some species, conserved regions are highly repetitive.

Conclusions

Instead of developing specific cement proteins for different attachment substrata, barnacles attached to different substrata rely on a highly duplicated cementation genetic toolkit to generate paralogous CP sequences with diverse chemical and biochemical properties. This general CP cocktail might be the key genetic feature enabling barnacles to adapt to a wide variety of substrata.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08049-4.

Genome and transcriptome assemblies of the kuruma shrimp, Marsupenaeus japonicus

The kuruma shrimp Marsupenaeus japonicus (order Decapoda, family Penaeidae) is an economically important crustacean that occurs in shallow, warm seas across the Indo-Pacific. Here, using a combination of Illumina and Oxford Nanopore Technologies platforms, we produced a draft genome assembly of M. japonicus (1.70 Gbp; 18,210 scaffolds; scaffold N50 = 234.9 kbp; 34.38% GC, 93.4% BUSCO completeness) and a complete mitochondrial genome sequence (15,969 bp). As with other penaeid shrimp genomes, the M. japonicus genome is extremely rich in simple repeats, which occupies 27.4% of the assembly. A total of 26,381 protein-coding gene models (94.7% BUSCO completeness) were predicted, of which 18,005 genes (68.2%) were assigned functional description by at least one method. We also produced an Illumina-based transcriptome shotgun assembly (40,991 entries; 93.0% BUSCO completeness) and a PacBio Iso-Seq transcriptome assembly (25,415 entries; 67.5% BUSCO completeness). We envision that the M. japonicus genome and transcriptome assemblies will serve as useful resources for the basic research, fisheries management, and breeding programs of M. japonicus.

Comparative analysis of stalked and acorn barnacle adhesive proteomes

Barnacles interest the scientific community for multiple reasons: their unique evolutionary trajectory, vast diversity and economic impact—as a harvested food source and also as one of the most prolific macroscopic hard biofouling organisms. A common, yet novel, trait among barnacles is adhesion, which has enabled a sessile adult existence and global colonization of the oceans. Barnacle adhesive is primarily composed of proteins, but knowledge of how the adhesive proteome varies across the tree of life is unknown due to a lack of genomic information. Here, we supplement previous mass spectrometry analyses of barnacle adhesive with recently sequenced genomes to compare the adhesive proteomes of Pollicipes pollicipes (Pedunculata) and Amphibalanus amphitrite (Sessilia). Although both species contain the same broad protein categories, we detail differences that exist between these species. The barnacle-unique cement proteins show the greatest difference between species, although these differences are diminished when amino acid composition and glycosylation potential are considered. By performing an in-depth comparison of the adhesive proteomes of these distantly related barnacle species, we show their similarities and provide a roadmap for future studies examining sequence-specific differences to identify the proteins responsible for functional differences across the barnacle tree of life.

Genomic insights into the sessile life and biofouling of barnacles (Crustacea: Cirripedia)

Members of the infraclass Cirripedia, commonly called barnacles, are unique among the subphylum Crustacea in that they exhibit a biphasic life cycle with a planktonic larval stage and a sessile adult stage. Understanding their unique sessile life and mechanisms of attachment are hampered by the lack of genomic resources. Here, we present a 746 Mb genome assembly of Lepas anserifera – the first sequenced stalked barnacle genome. We estimate that Cirripedia first arose ~495 million years ago (MYA) and further diversified since Mesozoic. A demographic analysis revealed remarkable population changes of the barnacle in relation to sea-level fluctuations in the last 2 MYA. Comparative genomic analyses revealed the expansion of a number of developmental related genes families in barnacle genomes, such as Br–C, PCP20 and Lola, which are potentially important for the evolution of metamorphosis, cuticle development and central nervous system. Phylogenetic analysis and tissue expression profiling showed the possible roles of gene duplication, functional diversification and co-option in shaping the genomic evolution of barnacles. Overall, our study provides not only a valuable draft genome for comparative genomic analysis of crustacean evolution, but also facilitates studies of biofouling control.

Genomic Characterization of the Barnacle Balanus improvisus Reveals Extreme Nucleotide Diversity in Coding Regions

Barnacles are key marine crustaceans in several habitats, and they constitute a common practical problem by causing biofouling on man-made marine constructions and ships. Despite causing considerable ecological and economic impacts, there is a surprising void of basic genomic knowledge, and a barnacle reference genome is lacking. We here set out to characterize the genome of the bay barnacle Balanus improvisus (= Amphibalanus improvisus) based on short-read whole-genome sequencing and experimental genome size estimation. We show both experimentally (DNA staining and flow cytometry) and computationally (k-mer analysis) that B. improvisus has a haploid genome size of ~ 740 Mbp. A pilot genome assembly rendered a total assembly size of ~ 600 Mbp and was highly fragmented with an N50 of only 2.2 kbp. Further assembly-based and assembly-free analyses revealed that the very limited assembly contiguity is due to the B. improvisus genome having an extremely high nucleotide diversity (π) in coding regions (average π ≈ 5% and average π in fourfold degenerate sites ≈ 20%), and an overall high repeat content (at least 40%). We also report on high variation in the α-octopamine receptor OctA (average π = 3.6%), which might increase the risk that barnacle populations evolve resistance toward antifouling agents. The genomic features described here can help in planning for a future high-quality reference genome, which is urgently needed to properly explore and understand proteins of interest in barnacle biology and marine biotechnology and for developing better antifouling strategies.

The Significance of Genotypic Diversity in Coral Competitive Interaction: A Transcriptomic Perspective

Competitive interactions shape coral assemblages and govern the dynamics of coral ecosystems. Although competition is an ecological concept, the outcomes of competitive interactions are ultimately determined by patterns of gene expression. These patterns are subject to genotypic variation on both sides of any interaction. Such variation is typically treated as “noise”, but it is sometimes possible to identify patterns within it that reveal important hidden factors in an experiment. To incorporate genotypic variation into the investigation of coral competitive interactions, we used RNA-sequencing to study changes in gene expression in a hard coral (Porites cylindrica) resulting from non-contact competition experiment with a soft coral (Lobophytum pauciflorum). Hard coral genotype explained the largest proportion of variation between samples; however, it was also possible to detect gene expression changes in 76 transcripts resulting from interaction with the soft coral. In addition, we found a group of 20 short secreted proteins that were expressed as a coordinated unit in three interacting Porites-Lobophytum pairs. The presence of this secretion response was idiosyncratic in that it could not be predicted based on polyp behaviour, or the genotype of hard or soft coral alone. This study illustrates the significance of individual variation as a determinant of competitive behaviour, and also provides some intriguing glimpses into the molecular mechanisms employed by hard corals competing at a distance.

Chitin Synthesis and Degradation in Crustaceans: A Genomic View and Application

Chitin is among the most important components of the crustacean cuticular exoskeleton and intestinal peritrophic matrix. With the progress of genomics and sequencing technology, a large number of gene sequences related to chitin metabolism have been deposited in the GenBank database in recent years. Here, we summarized the genes and pathways associated with the biosynthesis and degradation of chitins in crustaceans based on genomic analyses. We found that chitin biosynthesis genes typically occur in single or two copies, whereas chitin degradation genes are all multiple copies. Moreover, the chitinase genes are significantly expanded in most crustacean genomes. The gene structure and expression pattern of these genes are similar to those of insects, albeit with some specific characteristics. Additionally, the potential applications of the chitin metabolism genes in molting regulation and immune defense, as well as industrial chitin degradation and production, are also summarized in this review.