High levels of genetic variation in natural populations of marine lower invertebrates

High germline mutation rates, but not extreme population outbreaks, influence genetic diversity in a keystone coral predator

Lewontin's paradox, the observation that levels of genetic diversity (π) do not scale linearly with census population size (Nc) variation, is an evolutionary conundrum. The most extreme mismatches between π and Nc are found for highly abundant marine invertebrates. Yet, the influences of new mutations on π relative to extrinsic processes such as Nc fluctuations are unknown. Here, we provide the first germline mutation rate (μ) estimate for a marine invertebrate in corallivorous crown-of-thorns sea stars (Acanthaster cf. solaris). We use high-coverage whole-genome sequencing of 14 parent-offspring trios alongside empirical estimates of Nc in Australia's Great Barrier Reef to jointly examine the determinants of π in populations undergoing extreme Nc fluctuations. The A. cf. solaris mean μ was 9.13 x 10-09 mutations per-site per-generation (95% CI: 6.51 x 10-09 to 1.18 x 10-08), exceeding estimates for other invertebrates and showing greater concordance with vertebrate mutation rates. Lower-than-expected Ne (~70,000-180,000) and low Ne/Nc values (0.0047-0.048) indicated weak influences of population outbreaks on long-term π. Our findings are consistent with elevated μ evolving in response to reduced Ne and generation time length, with important implications for explaining high mutational loads and the determinants of genetic diversity in marine invertebrate taxa.

Genomic diversity of wild and cultured Yesso scallop Mizuhopecten yessoensis from Japan and Canada

The Yesso scallop Mizuhopecten yessoensis is an important aquaculture species that was introduced to Western Canada from Japan to establish an economically viable scallop farming industry. This highly fecund species has been propagated in Canadian aquaculture hatcheries for the past 40 years, raising questions about genetic diversity and genetic differences among hatchery stocks. In this study, we compare cultured Canadian and wild Japanese populations of Yesso scallop using double-digest restriction site-associated DNA (ddRAD) sequencing to genotype 21,048 variants in 71 wild-caught scallops from Japan, 65 scallops from the Vancouver Island University breeding population, and 37 scallops obtained from a commercial farm off Vancouver Island, British Columbia. The wild scallops are largely comprised of equally unrelated individuals, whereas cultured scallops are comprised of multiple families of related individuals. The polymorphism rate estimated in wild scallops was 1.7%, whereas in the cultured strains, it ranged between 1.35 and 1.07%. Interestingly, heterozygosity rates were highest in the cultured populations, which is likely due to shellfish hatchery practices of crossing divergent strains to gain benefits of heterosis and to avoid inbreeding. Evidence of founder effects and drift was observed in the cultured strains, including high genetic differentiation between cultured populations and between cultured populations and the wild population. Cultured populations had effective population sizes ranging from 9 to 26 individuals whereas the wild population was estimated at 25,048-56,291 individuals. Further, a depletion of low-frequency variants was observed in the cultured populations. These results indicate significant genetic diversity losses in cultured scallops in Canadian breeding programs.

The Iron-Responsive Genome of the Chiton Acanthopleura granulata

Molluscs biomineralize structures that vary in composition, form, and function, prompting questions about the genetic mechanisms responsible for their production and the evolution of these mechanisms. Chitons (Mollusca, Polyplacophora) are a promising system for studies of biomineralization because they build a range of calcified structures including shell plates and spine- or scale-like sclerites. Chitons also harden the calcified teeth of their rasp-like radula with a coat of iron (as magnetite). Here we present the genome of the West Indian fuzzy chiton Acanthopleura granulata, the first from any aculiferan mollusc. The A. granulata genome contains homologs of many genes associated with biomineralization in conchiferan molluscs. We expected chitons to lack genes previously identified from pathways conchiferans use to make biominerals like calcite and nacre because chitons do not use these materials in their shells. Surprisingly, the A. granulata genome has homologs of many of these genes, suggesting that the ancestral mollusc may have had a more diverse biomineralization toolkit than expected. The A. granulata genome has features that may be specialized for iron biomineralization, including a higher proportion of genes regulated directly by iron than other molluscs. A. granulata also produces two isoforms of soma-like ferritin: one is regulated by iron and similar in sequence to the soma-like ferritins of other molluscs, and the other is constitutively translated and is not found in other molluscs. The A. granulata genome is a resource for future studies of molluscan evolution and biomineralization.

The gene-rich genome of the scallop Pecten maximus

BACKGROUND

The king scallop, Pecten maximus, is distributed in shallow waters along the Atlantic coast of Europe. It forms the basis of a valuable commercial fishery and plays a key role in coastal ecosystems and food webs. Like other filter feeding bivalves it can accumulate potent phytotoxins, to which it has evolved some immunity. The molecular origins of this immunity are of interest to evolutionary biologists, pharmaceutical companies, and fisheries management.

FINDINGS

Here we report the genome assembly of this species, conducted as part of the Wellcome Sanger 25 Genomes Project. This genome was assembled from PacBio reads and scaffolded with 10X Chromium and Hi-C data. Its 3,983 scaffolds have an N50 of 44.8 Mb (longest scaffold 60.1 Mb), with 92% of the assembly sequence contained in 19 scaffolds, corresponding to the 19 chromosomes found in this species. The total assembly spans 918.3 Mb and is the best-scaffolded marine bivalve genome published to date, exhibiting 95.5% recovery of the metazoan BUSCO set. Gene annotation resulted in 67,741 gene models. Analysis of gene content revealed large numbers of gene duplicates, as previously seen in bivalves, with little gene loss, in comparison with the sequenced genomes of other marine bivalve species.

CONCLUSIONS

The genome assembly of P. maximus and its annotated gene set provide a high-quality platform for studies on such disparate topics as shell biomineralization, pigmentation, vision, and resistance to algal toxins. As a result of our findings we highlight the sodium channel gene Nav1, known to confer resistance to saxitoxin and tetrodotoxin, as a candidate for further studies investigating immunity to domoic acid.

Mitochondrial DNA is unsuitable to test for isolation by distance

Tests for isolation by distance (IBD) are the most commonly used method of assessing spatial genetic structure. Many studies have exclusively used mitochondrial DNA (mtDNA) sequences to test for IBD, but this marker is often in conflict with multilocus markers. Here, we report a review of the literature on IBD, with the aims of determining (a) whether significant IBD is primarily a result of lumping spatially discrete populations, and (b) whether microsatellite datasets are more likely to detect IBD when mtDNA does not. We also provide empirical data from four species in which mtDNA failed to detect IBD by comparing these with microsatellite and SNP data. Our results confirm that IBD is mostly found when distinct regional populations are pooled, and this trend disappears when each is analysed separately. Discrepancies between markers were found in almost half of the studies reviewed, and microsatellites were more likely to detect IBD when mtDNA did not. Our empirical data rejected the lack of IBD in the four species studied, and support for IBD was particularly strong for the SNP data. We conclude that mtDNA sequence data are often not suitable to test for IBD, and can be misleading about species' true dispersal potential. The observed failure of mtDNA to reliably detect IBD, in addition to being a single-locus marker, is likely a result of a selection-driven reduction in genetic diversity obscuring spatial genetic differentiation.

DOES COSMOPOLITANISM RESULT FROM OVERCONSERVATIVE SYSTEMATICS? A CASE STUDY USING THE MARINE SPONGE CHONDRILLA NUCULA

The sponge species Chondrilla nucula has a simple morphology and a very wide geographical distribution. To verify whether the latter might be an artifact of the former, samples of this species were collected across 10,000 km of its range, in the Mediterranean, the Caribbean, and the southwestern atlantic. The classical (spicule morphology) and molecular (allozymes) systematic approaches were compared, to try to define the geographic limits between populations and detect possible cryptic species. We found five distinct genetic forms within C. nucula that sometimes showed morphological homogeneity and other times plasticity. The difference in size of spicules could not be related to the clear-cut genetic differences, suggesting that the use of spicule sizes for sponge systematics should be reappraised. The population of one of the genetic forms along 3000 km of the Brazilian coast was highly structured (F_ST = 0.21; N_e m = 0.96). Our results reject the null hypothesis of cosmopolitanism of C. nucula and indicate that the putative worldwide distribution of some marine sponges, and possibly many other benthic invertebrates, may be the result of overly conservative systematics. Cryptic species appear to be particularly prevalent when genera are well defined but species are characterized by only a few morphological characters.

Genetic variability in the natural populations of Lasioderma serricorne (F.) (Coleoptera: Anobiidae), detected by RAPD markers and by esterase isozymes

Lasioderma serricorne (F.) is a small cosmopolitan beetle regarded as a destructive pest of several stored products such as grains, flour, spices, dried fruit and tobacco. Chemical insecticides are one of the measures used against the pest. However, intensive insecticide use has resulted in the appearance of resistant insect populations. Therefore, for the elaboration of more effective control programs, it is necessary to know the biological aspects of L. serricorne. Among these aspects, the genetic variability knowledge is very important and may help in the development of new control methods. The objective of this study was to evaluate the genetic variability of 11 natural populations of L. serricorne collected respectively in three and four towns in the states of Paraná and São Paulo, Brazil, using 20 primers random amplified polymorphic DNA (RAPD) and polymorphisms of esterases. These primers produced 352 polymorphic bands. Electrophoretic analysis of esterases allowed the identification of four polymorphic loci (Est-2, Est-4, Est-5 and Est-6) and 18 alleles. Results show that populations are genetically differentiated and there is a high level of genetic variability within populations. The high degree of genetic differentiation is not directly correlated to geographical distance. Thus, our data indicate that movement of infested commodities may contribute to the dissemination of L. serricorne, facilitating gene flow.

The global diversity of sea pens (Cnidaria: Octocorallia: Pennatulacea)

Recent advances in deep-sea exploration technology coupled with an increase in worldwide biotic surveys, biological research, and underwater photography in shallow water marine regions such as coral reefs, has allowed for a relatively rapid expansion of our knowledge in the global diversity of many groups of marine organisms. This paper is part of the PLoS ONE review collection of WoRMS (the Worldwide Register of Marine Species), on the global diversity of marine species, and treats the pennatulacean octocorals, a group of cnidarians commonly referred to as sea pens or sea feathers. This also includes sea pansies, some sea whips, and various vermiform taxa. Pennatulaceans are a morphologically diverse group with an estimated 200 or more valid species, displaying worldwide geographic and bathymetric distributions from polar seas to the equatorial tropics and from intertidal flats to over 6100 m in depth. The paper treats new discoveries and taxa new to science, and provides greater resolution in geographic and bathymetric distributions data than was previously known, as well as descriptions of life appearances in life and in situ observations at diverse depth.

Analysis of distributions of single-locus heterozygosity as a test of neutral theory

Three tests of neutral theory were carried out using a large dataset of vertebrate allozyme studies. The first test considered the relationship between the mean and variance of single locus heterozygosity across a sample of enzymes and non-enzymatic proteins. The second test compared the distributions of heterozygosity between paired proteins in balanced datasets in which each protein is scored for the same sample of species. The third test compared the observed distribution of single locus heterozygosity with theoretical distributions predicted by neutral theory. The results show an excellent quantitative fit with the predictions of neutral theory, though some small deviations from neutrality were observed which are consistent with the action of natural selection.

A restoration genetics guide for coral reef conservation

Worldwide degradation of coral reef communities has prompted a surge in restoration efforts. They proceed largely without considering genetic factors because traditionally, coral populations have been regarded as open over large areas with little potential for local adaptation. Since, biophysical and molecular studies indicated that most populations are closed over shorter time and smaller spatial scales. Thus, it is justified to re-examine the potential for site adaptation in corals. There is ample evidence for differentiated populations, inbreeding, asexual reproduction and the occurrence of ecotypes, factors that may facilitate local adaptation. Discovery of widespread local adaptation would influence coral restoration projects mainly with regard to the physical and evolutionary distance from the source wild and/or captive bred propagules may be moved without causing a loss of fitness in the restored population. Proposed causes for loss of fitness as a result of (plant) restoration efforts include founder effects, genetic swamping, inbreeding and/or outbreeding depression. Direct evidence for any of these processes is scarce in reef corals due to a lack of model species that allow for testing over multiple generations and the separation of the relative contributions of algal symbionts and their coral hosts to the overall performance of the coral colony. This gap in our knowledge may be closed by employing novel population genetic and genomics approaches. The use of molecular tools may aid managers in the selection of appropriate propagule sources, guide spatial arrangement of transplants, and help in assessing the success of coral restoration projects by tracking the performance of transplants, thereby generating important data for future coral reef conservation and restoration projects.

Single Nucleotide polymorphisms and their relationship to codon usage bias in the Pacific oyster Crassostrea gigas

DNA sequence polymorphism and codon usage bias were investigated in a set of 41 nuclear loci in the Pacific oyster Crassostrea gigas. Our results revealed a very high level of DNA polymorphism in oysters, in the order of magnitude of the highest levels reported in animals to date. A total of 290 single nucleotide polymorphisms (SNPs) were detected, 76 of which being localised in exons and 214 in non-coding regions. Average density of SNPs was estimated to be one SNP every 60 bp in coding regions and one every 40 bp in non-coding regions. Non-synonymous substitutions contributed substantially to the polymorphism observed in coding regions. The non-synonymous to silent diversity ratio was 0.16 on average, which is fairly higher to the ratio reported in other invertebrate species recognised to display large population sizes. Therefore, purifying selection does not appear to be as strong as it could have been expected for a species with a large effective population size. The level of non-synonymous diversity varied greatly from one gene to another, in accordance with varying selective constraints. We examined codon usage bias and its relationship with DNA polymorphism. The table of optimal codons was deduced from the analysis of an EST dataset, using EST counts as a rough assessment of gene expression. As recently observed in some other taxa, we found a strong and significant negative relationship between codon bias and non-synonymous diversity suggesting correlated selective constraints on synonymous and non-synonymous substitutions. Codon bias as measured by the frequency of optimal codons for expression might therefore provide a useful indicator of the level of constraint upon proteins in the oyster genome.

Population genetics of Collisella subrugosa (Patellogastropoda: Acmaeidae): evidence of two scales of population structure

Marine invertebrate populations usually show high levels of genetic variability that has frequently been associated with spatial and temporal environmental heterogeneity. One of the most heterogeneous marine environments is the intertidal zone, the habitat of Collisella subrugosa, the most widespread and abundant Brazilian limpet. C. subrugosa has planktonic larvae that can disperse over long distances, what can promote gene flow among shores, working against interpopulational differentiation. In this study we investigated the genetic variability and populational substructure of C. subrugosa through analysis of 24 allozyme loci in 14 samples (590 individuals) collected along 2,700 km of the Brazilian coast. The genetic variability was high ([Formula: see text] and [Formula: see text]), as expected for intertidal species. Genetic differentiation among samples was low (F (ST) = 0.03) what may reflect intensive gene flow associated with larval dispersal. However, we detected an isolation-by-distance pattern of population substructure in one sampled region. High levels of heterozygote deficiency were also observed for many loci in each sample. Alternative hypothesis are discussed, and the "breeding groups" is suggested to explain these pattern, indicating the main cause as environmental heterogeneity.

High level of genetic divergence between sympatric color morphs of the littoral sea anemone Anthopleura orientalis (Anthozoa: Actiniaria)

Using enzyme electrophoresis and nematocyst analysis, the sympatrically occurring "light" and "dark" color morphs of the sea anemone Anthopleura orientalis from the Sea of Japan were shown to be two valid species. The "light" morph was identified as A. orientalis (Averintsev, 1967 Issledovaniya fauny morei: Vyp. 5 (13). Nanka, Leningrad, pp. 62-77), while the "dark" morph was designated as Anthopleura sp. The analysis of 21 isozyme loci revealed high value of Nei's genetic distance (D=1.284) between the two species, which are indistinguishable in their external morphology. The mean values of observed and expected heterozygosities for A. orientalis and Anthopleura sp. are high (H(o)=0.252+/-0.061, H(e)=0.250+/-0.061 and H(o)=0.327+/-0.052, H(e)=0.351+/-0.054, respectively). The species differ significantly in the size of spirocysts and nematocysts, among which the atrichs from acrorhagi and the basitrichs from actinopharynx contribute most to the observed difference. Strong qualitative difference is revealed between distributions of nematocysts in mesenteric filaments of the two sea anemone species studied. The possible conspecificity of Anthopleura sp. with Anthopleura artemisia (Dana, 1848) is discussed and the conclusion made that these are two separate species.

Multiple Ig-like featuring genes divergent within and among individuals of the marine sponge Geodia cydonium

The receptor tyrosine kinase of the marine sponge Geodia cydonium features two extracellular Ig-like domains in which we recently documented RT-PCR polymorphism among individuals. Genomic-PCR analysis presented here revealed 14 unique sequences from four sponges, differing predominantly in the sequence of an intron which splits the Ig-like domains. Nevertheless, analysis of putative coding regions in 19 distinct clones (156-159 aa) from seven sponges revealed 69 positions of nucleotide substitutions, 67.6% of them non-synonymous, translating into 43 positions of divergent residues. Excluding aa deletions, these 19 sequences share pairwise aa identities of 89-99%. In three sponges, four or five unique Ig-like coding regions were scored. PCR amplification across this intron revealed multiple bands, polymorphic among five of six sponges. Further substantiated by Southern and Northern blots, it is evident that the genome of G. cydonium harbors multicopies of moderately divergent Ig-like domains. Presently, this only appears paralleled by the human KIR multigene family of NK cells MHC class I-specific receptors, which consist of two or three moderately divergent extracellular Ig-like domains.

The main protein of the aggregation factor responsible for species-specific cell adhesion in the marine sponge Microciona prolifera is highly polymorphic

Species-specific cell recognition in sponges, the oldest living metazoans, is based on a proteoglycan-like aggregation factor. We have screened individual sponge cDNA libraries, identifying multiple related forms for the aggregation factor core protein (MAFp3). Northern blots show the presence in several human tissues of transcripts strongly binding a MAFp3-specific probe. The open reading frame for MAFp3 is not interrupted in the 5' direction, revealing variable protein sequences that contain numerous introns equally spaced. We have studied tissue histocompatibility within a sponge population, finding 100% correlation between rejection behavior and the individual-specific restriction fragment length polymorphism pattern using aggregation factor-related probes. PCR amplifications with specific primers showed that at least some of the MAFp3 forms are allelic and distribute in the population used. A pronounced polymorphism is also observed when analyzing purified aggregation factor in polyacrylamide gels. Protease digestion of the polymorphic glycosaminoglycan-containing bands indicates that glycans are also responsible for the variability. The data presented reveal a high polymorphism of aggregation factor components, which matches the elevated sponge alloincompatibility, suggesting an involvement of the cell adhesion system in sponge allogeneic reactions.