DNA variation in the phenotypically-diverse brown alga Saccharina japonica

Recombinant Mitochondrial Genomes Reveal Recent Interspecific Hybridization between Invasive Salangid Fishes

The interspecific recombination of the mitochondrial (mt) genome, if not an experimental artifact, may result from interbreeding of species with broken reproductive barriers, which, in turn, is a frequent consequence of human activities including species translocations, habitat modifications, and climate change. This issue, however, has not been addressed for Protosalanx chinensis and other commercially important and, simultaneously, invasive salangid fishes that were the product of successful aquaculture in China. To assess the probability of interspecific hybridization, we analyzed the patterns of diversity and recombination in the complete mitochondrial (mt) genomes of these fishes using the GenBank resources. A sliding window analysis revealed a non-uniform distribution of the intraspecific differences in P. chinensis with four highly pronounced peaks of divergence centered at the COI, ND4L-ND4, and ND5 genes, and also at the control region. The corresponding divergent regions in P. chinensis show a high sequence similarity (99−100%) to the related salangid fishes, Neosalanx tangkahkeii and N. anderssoni. This observation suggests that the divergent regions of P. chinensis may represent a recombinant mitochondrial DNA (mtDNA) containing mt genome fragments belonging to different salangid species. Indeed, four, highly significant (pairwise homoplasy index test, P < 0.00001) signals of recombination have been revealed at coordinates closely corresponding to the divergent regions. The recombinant fragments are, however, not fixed, and different mt genomes of P. chinensis are mosaic, containing different numbers of recombinant events. These facts, along with the high similarity or full identity of the recombinant fragments between the donor and the recipient sequences, indicate a recent interspecific hybridization between P. chinensis and two Neosalanx species. Alternative hypotheses, including taxonomical misidentifications, sequence misalignments, DNA contamination, and/or artificial PCR recombinants, are not supported by the data. The recombinant fragments revealed in our study represent diagnostic genetic markers for the identification and distinguishing of hybrids, which can be used to control the invasive dynamics of hybrid salangid fishes.

Genetic Evidence for a Mixed Composition of the Genus Myoxocephalus (Cottoidei: Cottidae) Necessitates Generic Realignment

Sculpin fishes belonging to the family Cottidae represent a large and complex group, inhabiting a wide range of freshwater, brackish-water, and marine environments. Numerous studies based on analysis of their morphology and genetic makeup frequently provided controversial results. In the present work, we sequenced complete mitochondrial (mt) genomes and fragments of nuclear ribosomal DNA (rDNA) of the fourhorn sculpin Myoxocephalus quadricornis and some related cottids to increase the power of phylogenetic and taxonomic analyses of this complex fish group. A comparison of the My. quadricornis mt genomes obtained by us with other complete mt genomes available in GenBank has revealed a surprisingly low divergence (3.06 ± 0.12%) with Megalocottus platycephalus and, at the same time, a significantly higher divergence (7.89 ± 0.16%) with the species of the genus Myoxocephalus. Correspondingly, phylogenetic analyses have shown that My. quadricornis is clustered with Me. platycephalus but not with the Myoxocephalus species. Completely consistent patterns of divergence and tree topologies have been obtained based on nuclear rDNA. Thus, the multi-gene data in the present work indicates obvious contradictions in the relationships between the Myoxocephalus and Megalocottus species studied. An extensive phylogenetic analysis has provided evidence for a closer affinity of My. quadricornis with the species of the genus Megalocottus than with the species of the genus Myoxocephalus. A recombination analysis, along with the additional GenBank data, excludes introgression and/or incorrect taxonomic identification as the possible causative factors responsible for the observed closer affinity between the two species from different genera. The above facts necessitate realignment of the genera Myoxocephalus and Megalocottus. The genetic data supports the two recognized genera, Myoxocephalus and Megalocottus, but suggests changing their compositions through transferring My. quadricornis to the genus Megalocottus. The results of the present study resolve the relationships within a complex group of sculpin fishes and show a promising approach to phylogenetic systematics (as a key organizing principle in biodiversity research) for a better understanding of the taxonomy and evolution of fishes and for supplying relevant information to address various fish biodiversity conservation and management issues.

Trans-Species Polymorphism in Mitochondrial Genome of Camarodont Sea Urchins

Mitochondrial (mt) genomes of the sea urchins Strongylocentrotus intermedius and Mesocentrotus nudus demonstrate the identical patterns of intraspecific length variability of the ND6 gene, consisting of 489 bp (S variant) and 498 bp (L variant), respectively. For both species, the ND6 length difference is due to the 488A>G substitution, which changes the stop codon TAG in S variant for a tryptophan codon TGG in L variant and elongates the corresponding ND6 protein by three additional amino acids, Trp-Leu-Trp. The phylogenetic analysis based on mt genomes of sea urchins and related echinoderm groups from GenBank has shown the S and L ND6 variants as shared among the camarodont sea urchins; the rest of the echinoderms demonstrate the S variant only. The data suggest that the ND6 488A>G substitution can be the first example of the trans-species polymorphism in sea urchins, persisting at least since the time of the Odontophora diversification at the Eocene/Oligocene boundary (approximately 34 million years ago), which was characterized by an abrupt climate change and significant global ocean cooling. Alternative hypotheses, including the convergent RNA editing and/or codon reassignment, are not supported by direct comparisons of the ND6 gene sequences with the corresponding transcripts using the basic local alignment search tool (BLAST) of full sea urchin transcriptomes.

Phylogeographic diversification and postglacial range dynamics shed light on the conservation of the kelp Saccharina japonica

Studies of postglacial range shifts could enhance our understanding of seaweed species' responses to climate change and hence facilitate the conservation of natural resources. However, the distribution dynamics and phylogeographic diversification of the commercially and ecologically important kelp Saccharina japonica in the Northwest Pacific (NWP) are still poorly surveyed. In this study, we analyzed the evolutionary history of S. japonica using two mitochondrial markers and 24 nuclear microsatellites. A STRUCTURE analysis revealed two partially isolated lineages: lineage H, which is scattered along the coast of Japan; and lineage P, which occurs along the west coast of the Japan Sea. Ecological niche modeling projections to the Last Glacial Maximum (LGM) revealed that the southern coasts of the Japan Sea and the Pacific side of the Oshima and Honshu Peninsulas provided the most suitable habitats for S. japonica, implying that these regions served as ancient refugia during the LGM. Ancient isolation in different refugia may explain the observed divergence between lineages P and H. An approximate Bayesian computation analysis indicated that the two lineages experienced post-LGM range expansion and that postglacial secondary contact occurred in Sakhalin. Model projections into the year 2,100 predicted that S. japonica will shift northwards and lose its genetic diversity center on the Oshima Peninsula in Hokkaido and Shimokita Peninsula in Honshu. The range shifts and evolutionary history of S. japonica improve our understanding of how climate change impacted the distribution range and diversity of this species and provide useful information for the conservation of natural resources under ongoing environmental change in the NWP.

Strong population structure but no equilibrium yet: Genetic connectivity and phylogeography in the kelp Saccharina latissima (Laminariales, Phaeophyta)

Kelp aquaculture is globally developing steadily as human food source, along with other applications. One of the newer crop species is Saccharina latissima, a northern hemisphere kelp inhabiting temperate to arctic rocky shores. To protect and document its natural genetic variation at the onset of this novel aquaculture, as well as increase knowledge on its taxonomy and phylogeography, we collected new genetic data, both nuclear and mitochondrial, and combined it with previous knowledge to estimate genetic connectivity and infer colonization history. Isolation‐with‐migration coalescent analyses demonstrate that gene flow among the sampled locations is virtually nonexistent. An updated scenario for the origin and colonization history of S. latissima is developed as follows: We propose that the species (or species complex) originated in the northwest Pacific, crossed to the northeast Pacific in the Miocene, and then crossed the Bering Strait after its opening ~5.5 Ma into the Arctic and northeast Atlantic. It subsequently crossed the Atlantic from east to west. During the Pleistocene, it was compressed in the south with evidence for northern refugia in Europe. Postglacial recolonization led to secondary contact in the Canadian Arctic. Saccharina cichorioides is shown to probably belong to the S. latissima species complex and to derive from ancestral populations in the Asian North Pacific. Our novel approach of comparing inferred gene flow based on coalescent analysis versus Wright's island model suggests that equilibrium levels of differentiation have not yet been reached in Europe and, hence, that genetic differentiation is expected to increase further if populations are left undisturbed.

Glacial vicariance drives phylogeographic diversification in the amphi-boreal kelp Saccharina latissima

Glacial vicariance is regarded as one of the most prevalent drivers of phylogeographic structure and speciation among high-latitude organisms, but direct links between ice advances and range fragmentation have been more difficult to establish in marine than in terrestrial systems. Here we investigate the evolution of largely disjunct (and potentially reproductively isolated) phylogeographic lineages within the amphi-boreal kelp Saccharina latissima s. l. Using molecular data (COI, microsatellites) we confirm that S. latissima comprises also the NE Pacific S. cichorioides complex and is composed of divergent lineages with limited range overlap and genetic admixture. Only a few genetic hybrids were detected throughout a Canadian Arctic/NW Greenland contact zone. The degree of genetic differentiation and sympatric isolation of phylogroups suggest that S. latissima s. l. represents a complex of incipient species. Phylogroup distributions compared with paleo-environmental reconstructions of the cryosphere further suggest that diversification within S. latissima results from chronic glacial isolation in disjunct persistence areas intercalated with ephemeral interglacial poleward expansions and admixture at high-latitude (Arctic) contact zones. This study thus supports a role for glaciations not just in redistributing pre-existing marine lineages but also as a speciation pump across multi-glacial cycles for marine organisms otherwise exhibiting cosmopolite amphi-boreal distributions.

Complete Mitochondrial Genomes of the Cherskii's Sculpin Cottus czerskii and Siberian Taimen Hucho taimen Reveal GenBank Entry Errors: Incorrect Species Identification and Recombinant Mitochondrial Genome

The complete mitochondrial (mt) genome is sequenced in 2 individuals of the Cherskii's sculpin Cottus czerskii. A surprisingly high level of sequence divergence (10.3%) has been detected between the 2 genomes of C czerskii studied here and the GenBank mt genome of C czerskii (KJ956027). At the same time, a surprisingly low level of divergence (1.4%) has been detected between the GenBank C czerskii (KJ956027) and the Amur sculpin Cottus szanaga (KX762049, KX762050). We argue that the observed discrepancies are due to incorrect taxonomic identification so that the GenBank accession number KJ956027 represents actually the mt genome of C szanaga erroneously identified as C czerskii. Our results are of consequence concerning the GenBank database quality, highlighting the potential negative consequences of entry errors, which once they are introduced tend to be propagated among databases and subsequent publications. We illustrate the premise with the data on recombinant mt genome of the Siberian taimen Hucho taimen (NCBI Reference Sequence Database NC_016426.1; GenBank accession number HQ897271.1), bearing 2 introgressed fragments (≈0.9 kb [kilobase]) from 2 lenok subspecies, Brachymystax lenok and Brachymystax lenok tsinlingensis, submitted to GenBank on June 12, 2011. Since the time of submission, the H taimen recombinant mt genome leading to incorrect phylogenetic inferences was propagated in multiple subsequent publications despite the fact that nonrecombinant H taimen genomes were also available (submitted to GenBank on August 2, 2014; KJ711549, KJ711550). Other examples of recombinant sequences persisting in GenBank are also considered. A GenBank Entry Error Depositary is urgently needed to monitor and avoid a progressive accumulation of wrong biological information.

A morphological and phylogenetic investigation into divergence among sympatric Australian southern bull kelps (Durvillaea potatorum and D. amatheiae sp. nov.)

Genetic analyses can reveal a wealth of hitherto undiscovered cryptic biodiversity. For co-occurring and morphologically similar species, the combination of molecular, ecological and morphological analyses provides an excellent opportunity for understanding some of the processes that can lead to divergence and speciation. The Australian endemic brown macroalga Durvillaea potatorum (Phaeophyceae) was examined with a combination of genetic and morphological approaches to confirm the presence of two separate species and to infer the processes that led to their divergence. A total of 331 individuals from 11 sites around coastal Tasmania were collected and measured in situ for a range of morphological and ecological characteristics. Tissue samples were also collected for each individual to allow genetic analyses using mitochondrial (COI) and nuclear (28S) markers. Genetic analyses confirmed the presence of two deeply divergent clades. The significant morphological differentiation, despite high levels of intra-lineage variability, further supported their recognition as distinct species. We describe a new species, D. amatheiae sp. nov., which is characterised by a narrower and proportionately shorter stipe, shorter total length, and higher number of stipitate lateral blades and branches than D. potatorum (sensu stricto). The occurrence of both species in sympatry along Tasmania's eastern and western coasts, as well as their contrasting patterns of haplotype diversity, supports a hypothesis of geographical isolation, allopatric speciation and subsequent secondary contact in response to sea level and ocean current change throughout the Pleistocene glaciation cycles. This research contributes to resolving the phylogenetic relationships, taxonomy and evolution of the ecologically keystone kelp genus Durvillaea.

Identification of Genes under Positive Selection Reveals Differences in Evolutionary Adaptation between Brown-Algal Species

Brown algae are an important taxonomic group in coastal ecosystems. The model brown algal species Ectocarpus siliculosus and Saccharina japonica are closely related lineages. Despite their close phylogenetic relationship, they vary greatly in morphology and physiology. To obtain further insights into the evolutionary forces driving divergence in brown algae, we analyzed 3,909 orthologs from both species to identify Genes Under Positive Selection (GUPS). About 12% of the orthologs in each species were considered to be under positive selection. Many GUPS are involved in membrane transport, regulation of homeostasis, and sexual reproduction in the small sporophyte of E. siliculosus, which is known to have a complex life cycle and to occupy a wide range of habitats. Genes involved in photosynthesis and cell division dominated the group of GUPS in the large kelp of S. japonica, which might explain why this alga has evolved the ability to grow very rapidly and to form some of the largest sporophytes. A significant number of molecular chaperones (e.g., heat-shock proteins) involved in stress responses were identified to be under positive selection in both species, potentially indicating their important roles for macroalgae to cope with the relatively variable environment of coastal ecosystems. Moreover, analysis of previously published microarray data of E. siliculosus showed that many GUPS in E. siliculosus were responsive to stress conditions, such as oxidative and hyposaline stress, whereas our RNA-seq data of S. japonica showed that GUPS in this species were most highly expressed in large sporophytes, which supports the suggestion that selection largely acts on different sets of genes in both marcoalgal species, potentially reflecting their adaptation to different ecological niches.

DNA polymorphism and selection at the bindin locus in three Strongylocentrotus sp. (Echinoidea)

BACKGROUND

The sperm gene bindin encodes a gamete recognition protein, which plays an important role in conspecific fertilization and reproductive isolation of sea urchins. Molecular evolution of the gene has been extensively investigated with the attention focused on the protein coding regions. Intron evolution has been investigated to a much lesser extent. We have studied nucleotide variability in the complete bindin locus, including two exons and one intron, in the sea urchin Strongylocentrotus intermedius represented by two morphological forms. We have also analyzed all available bindin sequences for two other sea urchin species, S. pallidus and S. droebachiensis.

RESULTS

The results show that the bindin sequences from the two forms of S. intermedius are intermingled with no evidence of genetic divergence; however, the forms exhibit slightly different patterns in bindin variability. The level of the bindin nucleotide diversity is close for S. intermedius and S. droebachiensis, but noticeably higher for S. pallidus. The distribution of variability is non-uniform along the gene; however there are striking similarities among the species, indicating similar evolutionary trends in this gene engaged in reproductive function. The patterns of nucleotide variability and divergence are radically different in the bindin coding and intron regions. Positive selection is detected in the bindin coding region. The neutrality tests as well as the maximum likelihood approaches suggest the action of diversifying selection in the bindin intron.

CONCLUSIONS

Significant deviation from neutrality has been detected in the bindin coding region and suggested in the intron, indicating the possible functional importance of the bindin intron variability. To clarify the question concerning possible involvement of diversifying selection in the bindin intron evolution more data combining population genetic and functional approaches are necessary.

Phylogeographic data revealed shallow genetic structure in the kelp Saccharina japonica (Laminariales, Phaeophyta)

BACKGROUND

Population structure and genetic diversity of marine organisms in the Northwestern Pacific Ocean exhibited complex patterns. Saccharina japonica is a commercially and ecologically important kelp species widely distributed along the coast of Japan Sea. However, it is still poorly known about population genetics and phylogeographic patterns of wild S. japonica populations on a large geographic scale, which is an important contribution to breeding and conservation of this marine crop.

RESULTS

We collected 612 mitochondrial COI and trnW-trnL sequences. Diversity indices suggested that S. japonica populations along the coast of Hokkaido exhibited the highest genetic diversity. Bayesian Analysis of Population Structure (BAPS) revealed four clusters in the kelp species (cluster 1: Hokkaido and South Korea; cluster 2: northwestern Hokkaido; cluster 3: Far Eastern Russia; cluster 4: China). The network inferred from concatenated data exhibited two shallow genealogies corresponding to two BAPS groups (cluster 2 and cluster 3). We did not detect gene flow between the two shallow genealogies, but populations within genealogy have asymmetric gene exchange. Bayesian skyline plots and neutrality tests suggested that S. japonica experienced postglacial expansion around 10.45 ka.

CONCLUSIONS

The coast of Hokkaido might be the origin and diversification center of S. japonica. Gene exchange among S. japonica populations could be caused by anthropogenic interference and oceanographic regimes. Postglacial expansions and gene exchange apparently led to more shared haplotypes and less differentiation that in turn led to the present shallow phylogeographical patterns in S. japonica.

Symbiotic associations in the phenotypically-diverse brown alga Saccharina japonica

The brown alga Saccharina japonica (Areschoug) Lane, Mayes, Druehl et Saunders is a highly polymorphic representative of the family Laminariaceae, inhabiting the northwest Pacific region. We have obtained 16S rRNA sequence data in symbiont microorganisms of the typical form (TYP) of S. japonica and its common morphological varieties, known as "longipes" (LON) and "shallow-water" (SHA), which show contrasting bathymetric distribution and sharp morphological, life history traits, and ecological differences. Phylogenetic analysis of the 16S rRNA sequences shows that the microbial communities are significantly different in the three forms studied and consist of mosaic sets of common and form-specific bacterial lineages. The divergence in bacterial composition is substantial between the TYP and LON forms in spite of their high genetic similarity. The symbiont distribution in the S. japonica forms and in three other laminarialean species is not related to the depth or locality of the algae settlements. Combined with our previous results on symbiont associations in sea urchins and taking into account the highly specific character of bacteria-algae associations, we propose that the TYP and LON forms may represent incipient species passing through initial steps of reproductive isolation. We suggest that phenotype differences between genetically similar forms may be caused by host-symbiont interactions that may be a general feature of evolution in algae and other eukaryote organisms. Bacterial symbionts could serve as sensitive markers to distinguish genetically similar algae forms and also as possible growth-promoting inductors to increase algae productivity.