Footprints of natural selection at the mannose-6-phosphate isomerase locus in barnacles

New genes helped acorn barnacles adapt to a sessile lifestyle

Barnacles are the only sessile lineages among crustaceans, and their sessile life begins with the settlement of swimming larvae (cyprids) and the formation of protective shells. These processes are crucial for adaptation to a sessile lifestyle, but the underlying molecular mechanisms remain poorly understood. While investigating these mechanisms in the acorn barnacle, Amphibalanus amphitrite, we discovered a new gene, bcs-6, which is involved in the energy metabolism of cyprid settlement and originated from a transposon by acquiring the promoter and cis-regulatory element. Unlike mollusks, the barnacle shell comprises alternate layers of chitin and calcite and requires another new gene, bsf, which generates silk-like fibers that efficiently bind chitin and aggregate calcite in the aquatic environment. Our findings highlight the importance of exploring new genes in unique adaptative scenarios, and the results will provide important insights into gene origin and material development.

Genetic Diversity, Population Structure, and Environmental Adaptation Signatures of Chinese Coastal Hard-Shell Mussel Mytilus coruscus Revealed by Whole-Genome Sequencing

The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (FST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

Parasite manipulation of host phenotypes inferred from transcriptional analyses in a trematode-amphipod system

Manipulation of host phenotypes by parasites is hypothesized to be an adaptive strategy enhancing parasite transmission across hosts and generations. Characterizing the molecular mechanisms of manipulation is important to advance our understanding of host-parasite coevolution. The trematode (Levinseniella byrdi) is known to alter the colour and behaviour of its amphipod host (Orchestia grillus) presumably increasing predation of amphipods which enhances trematode transmission through its life cycle. We sampled 24 infected and 24 uninfected amphipods from a salt marsh in Massachusetts to perform differential gene expression analysis. In addition, we constructed novel genomic tools for O. grillus including a de novo genome and transcriptome. We discovered that trematode infection results in upregulation of amphipod transcripts associated with pigmentation and detection of external stimuli, and downregulation of multiple amphipod transcripts implicated in invertebrate immune responses, such as vacuolar ATPase genes. We hypothesize that suppression of immune genes and the altered expression of genes associated with coloration and behaviour may allow the trematode to persist in the amphipod and engage in further biochemical manipulation that promotes transmission. The genomic tools and transcriptomic analyses reported provide new opportunities to discover how parasites alter diverse pathways underlying host phenotypic changes in natural populations.

The molecular background of the aspartate aminotransferase polymorphism in Littorina snails maintained by strong selection on small spatial scales

Allozymes present several classical examples of divergent selection, including the variation in the cytosolic aspartate aminotransferase (AAT) in the intertidal snails Littorina saxatilis. AAT is a part of the asparate-malate shuttle, in the interidal molluscs involved in the anaerobic respiration during desiccation. Previous allozyme studies reported the sharp gradient in the frequencies of the AAT¹⁰⁰and the AAT¹²⁰ alleles between the low and high shores in the Northern Europe and the differences in their enzymatic activity, supporting the role of AAT in adaptation to desiccation. However, the populations in the Iberian Peninsula showed the opposite allele cline. Using the mRNA sequencing and the genome pool-seq analyses we characterize DNA sequences of the different AAT alleles, report the amino acid replacements behind the allozyme variation and show that same allozyme alleles in Northern and Southern populations have different protein sequences. Gene phylogeny reveals that the AAT¹⁰⁰ and the northern AAT¹²⁰ alleles represent the old polymorphism, shared among the closely related species of Littorina, while the southern AAT¹²⁰ allele is more recently derived from AAT¹⁰⁰. Further, we show that the Aat gene is expressed at constitutive level in different genotypes and conditions, supporting the role of structural variation in regulation of enzyme activity. Finally, we report the location and the structure of the gene in the L. saxatilis genome and the presence of two additional non-functional gene copies. Altogether, we provide a missing link between the classical allozyme studies and the genome scans and bring together the results produced over decades of the genetic research.

An allozyme polymorphism is associated with a large chromosomal inversion in the marine snail Littorina fabalis

Understanding the genetic targets of natural selection is one of the most challenging goals of population genetics. Some of the earliest candidate genes were identified from associations between allozyme allele frequencies and environmental variation. One such example is the clinal polymorphism in the arginine kinase (Ak) gene in the marine snail Littorina fabalis. While other enzyme loci do not show differences in allozyme frequencies among populations, the Ak alleles are near differential fixation across repeated wave exposure gradients in Europe. Here, we use this case to illustrate how a new sequencing toolbox can be employed to characterize the genomic architecture associated with historical candidate genes. We found that the Ak alleles differ by nine nonsynonymous substitutions, which perfectly explain the different migration patterns of the allozymes during electrophoresis. Moreover, by exploring the genomic context of the Ak gene, we found that the three main Ak alleles are located on different arrangements of a putative chromosomal inversion that reaches near fixation at the opposing ends of two transects covering a wave exposure gradient. This shows Ak is part of a large (3/4 of the chromosome) genomic block of differentiation, in which Ak is unlikely to be the only target of divergent selection. Nevertheless, the nonsynonymous substitutions among Ak alleles and the complete association of one allele with one inversion arrangement suggest that the Ak gene is a strong candidate to contribute to the adaptive significance of the inversion.

Lack of a genetic cline and temporal genetic stability in an introduced barnacle along the Pacific coast of Japan

Background

Large numbers of exotic marine species have been introduced worldwide. Monitoring of introduced species is important to reveal mechanisms underlying their establishment and expansion. Balanus glandula is a common intertidal barnacle native to the northeastern Pacific. However, this species has been introduced to Japan, South America, South Africa, and Europe. While a latitudinal genetic cline is well known in its native range, it is unclear whether such a genetic cline occurs in introduced areas. Twenty years have passed since it was first identified in Japan and its distribution now ranges from temperate to subarctic regions.

Methods

In the present study, we examined genotypes of cytochrome oxidase subunit I (COI) of mitochondrial (mt)-DNA and elongation factor 1a (EF1) across the distribution of B. glandula in Japan at high and mid intertidal zones.

Results

At all sampling sites, native northern genotypes are abundant and I did not detect significant effects of latitude, tide levels, or their interaction on genotypic frequencies. Further, I did not detect any change of genotype composition between data collected during a study in 2004 and samples in the present study collected in 2019. Data from the present study offer an important baseline for future monitoring of this species and supply valuable insights into the mechanisms of establishment and expansion of introduced marine species generally.

From tides to nucleotides: Genomic signatures of adaptation to environmental heterogeneity in barnacles

The northern acorn barnacle (Semibalanus balanoides) is a robust system to study the genetic basis of adaptations to highly heterogeneous environments. Adult barnacles may be exposed to highly dissimilar levels of thermal stress depending on where they settle in the intertidal (i.e., closer to the upper or lower tidal boundary). For instance, barnacles near the upper tidal limit experience episodic summer temperatures above recorded heat coma levels. This differential stress at the microhabitat level is also dependent on the aspect of sun exposure. In the present study, we used pool-seq approaches to conduct a genome wide screen for loci responding to intertidal zonation across the North Atlantic basin (Maine, Rhode Island, and Norway). Our analysis discovered 382 genomic regions containing SNPs which are consistently zonated (i.e., SNPs whose frequencies vary depending on their position in the rocky intertidal) across all surveyed habitats. Notably, most zonated SNPs are young and private to the North Atlantic. These regions show high levels of genetic differentiation across ecologically extreme microhabitats concomitant with elevated levels of genetic variation and Tajima's D, suggesting the action of non-neutral processes. Overall, these findings support the hypothesis that spatially heterogeneous selection is a general and repeatable feature for this species, and that natural selection can maintain functional genetic variation in heterogeneous environments.

Ecological Load and Balancing Selection in Circumboreal Barnacles

Acorn barnacle adults experience environmental heterogeneity at various spatial scales of their circumboreal habitat, raising the question of how adaptation to high environmental variability is maintained in the face of strong juvenile dispersal and mortality. Here, we show that 4% of genes in the barnacle genome experience balancing selection across the entire range of the species. Many of these genes harbor mutations maintained across 2 My of evolution between the Pacific and Atlantic oceans. These genes are involved in ion regulation, pain reception, and heat tolerance, functions which are essential in highly variable ecosystems. The data also reveal complex population structure within and between basins, driven by the trans-Arctic interchange and the last glaciation. Divergence between Atlantic and Pacific populations is high, foreshadowing the onset of allopatric speciation, and suggesting that balancing selection is strong enough to maintain functional variation for millions of years in the face of complex demography.

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.

Rapid Adaptation to Temperature via a Potential Genomic Island of Divergence in the Invasive Green Crab, Carcinus maenas

Widespread species often adapt easily to novel conditions – both those found in new habitats and those generated by climate change. However, rapid adaptation may be hindered in the marine realm, where long-distance dispersal and consequently high gene flow are predicted to limit potential for local adaptation. Here, we use a highly dispersive invasive marine crab to test the nature and speed of adaptation to temperature in the sea. Using single nucleotide polymorphisms (SNPs) generated from cardiac transcriptome sequencing, we characterized six populations of the European green crab (Carcinus maenas) located across parallel thermal gradients in their native and invasive ranges. We compared SNP frequencies with local temperatures and previously generated data on cardiac heat and cold tolerance to identify candidate markers associated with population-level differences in thermal physiology. Of 10,790 SNPs, 104 were identified as frequency outliers, a signal that was strongly driven by association with temperature and/or cold tolerance. Seventy-two of these outlier markers, representing 28 different genes, were in a cluster of SNPs identified as a potential inversion polymorphism using linkage disequilibrium network analysis. This SNP cluster was unique in the data set, which was otherwise characterized by low levels of linkage disequilibrium, and markers in this cluster showed a significant enrichment of coding substitutions relative to the full SNP set. These 72 outlier SNPs appear to be transmitted as a unit, and represent a putative genomic island of divergence which varied in frequency with organismal cold tolerance. This relationship was strikingly similar across both native and invasive populations, all of which showed a very strong correlation with cold tolerance (R2 = 0.96 over all six populations). Notably, three of these populations have diverged recently (<100 years) and show little to no neutral divergence, suggesting that this genomic region may be responding to temperature on a relatively short time scale. This relationship indicates adaptation to temperature based on the action of a putative genomic island of divergence, perhaps partially explaining the extraordinary invasive ability of this species.

Characterizing the cirri and gut microbiomes of the intertidal barnacle Semibalanus balanoides

BACKGROUND

Natural populations inhabiting the rocky intertidal experience multiple ecological stressors and provide an opportunity to investigate how environmental differences influence microbiomes over small geographical scales. However, very few microbiome studies focus on animals that inhabit the intertidal. In this study, we investigate the microbiome of the intertidal barnacle Semibalanus balanoides. We first describe the microbiome of two body tissues: the feeding appendages, or cirri, and the gut. Next, we examine whether there are differences between the microbiome of each body tissue of barnacles collected from the thermally extreme microhabitats of the rocky shores' upper and lower tidal zones.

RESULTS

Overall, the microbiome of S. balanoides consisted of 18 phyla from 408 genera. Our results showed that although cirri and gut microbiomes shared a portion of their amplicon sequence variants (ASVs), the microbiome of each body tissue was distinct. Over 80% of the ASVs found in the cirri were also found in the gut, and 44% of the ASVs found in the gut were also found in the cirri. Notably, the gut microbiome was not a subset of the cirri microbiome. Additionally, we identified that the cirri microbiome was responsive to microhabitat differences.

CONCLUSION

Results from this study indicate that S. balanoides maintains distinct microbiomes in its cirri and gut tissues, and that the gut microbiome is more stable than the cirri microbiome between the extremes of the intertidal.

Footprints of natural selection at the mannose-6-phosphate isomerase locus in barnacles

The rocky intertidal is a natural laboratory to study how natural selection acts on the genes and proteins responsible for organismal survival and reproduction. Alternative forms of enzymes that differ across the intertidal have been known for decades and have provided examples of selection, but the genetic basis of such enzyme variation is known in only a few cases. In this paper, we present molecular evidence of natural selection at the Mpi gene, a key enzyme in energy metabolism that alters survival of barnacles living across the stress gradient imposed by the intertidal. Our study demonstrates how natural selection can facilitate survival in highly heterogeneous environments through the maintenance of multiple molecular solutions to ecological stresses.

The mannose-6-phosphate isomerase (Mpi) locus in Semibalanus balanoides has been studied as a candidate gene for balancing selection for more than two decades. Previous work has shown that Mpi allozyme genotypes (fast and slow) have different frequencies across Atlantic intertidal zones due to selection on postsettlement survival (i.e., allele zonation). We present the complete gene sequence of the Mpi locus and quantify nucleotide polymorphism in S. balanoides, as well as divergence to its sister taxon Semibalanus cariosus. We show that the slow allozyme contains a derived charge-altering amino acid polymorphism, and both allozyme classes correspond to two haplogroups with multiple internal haplotypes. The locus shows several footprints of balancing selection around the fast/slow site: an enrichment of positive Tajima’s D for nonsynonymous mutations, an excess of polymorphism, and a spike in the levels of silent polymorphism relative to silent divergence, as well as a site frequency spectrum enriched for midfrequency mutations. We observe other departures from neutrality across the locus in both coding and noncoding regions. These include a nonsynonymous trans-species polymorphism and a recent mutation under selection within the fast haplogroup. The latter suggests ongoing allelic replacement of functionally relevant amino acid variants. Moreover, predicted models of Mpi protein structure provide insight into the functional significance of the putatively selected amino acid polymorphisms. While footprints of selection are widespread across the range of S. balanoides, our data show that intertidal zonation patterns are variable across both spatial and temporal scales. These data provide further evidence for heterogeneous selection on Mpi.