Are Prorocentrum hoffmannianum and Prorocentrum belizeanum (Dinophyceae, Prorocentrales), the same species? An integration of morphological and molecular data.

The taxonomic assignment of Prorocentrum species is based on morphological characteristics; however, morphological variability has been found for several taxa isolated from different geographical regions. In this study, we evaluated species boundaries of Prorocentrum hoffmannianum and Prorocentrum belizeanum based on morphological and molecular data. A detailed morphological analysis was done, concentrating on the periflagellar architecture. Molecular analyses were performed on partial Small Sub-Unit (SSU) rDNA, partial Large Sub-Unit (LSU) rDNA, complete Internal Transcribed Spacer Regions (ITS1-5.8S-ITS2), and partial cytochrome b (cob) sequences. We concatenated the SSU-ITS-LSU fragments and constructed a phylogenetic tree using Bayesian Inference (BI) and maximum likelihood (ML) methods. Morphological analyses indicated that the main characters, such as cell size and number of depressions per valve, normally used to distinguish P. hoffmannianum from P. belizeanum, overlapped. No clear differences were found in the periflagellar area architecture. Prorocentrum hoffmannianum and P. belizeanum were a highly supported monophyletic clade separated into three subclades, which broadly corresponded to the sample collection regions. Subtle morphological overlaps found in cell shape, size, and ornamentation lead us to conclude that P. hoffmanianum and P. belizeanum might be considered conspecific. The molecular data analyses did not separate P. hoffmannianum and P. belizeanum into two morphospecies, and thus, we considered them to be the P. hoffmannianum species complex because their clades are separated by their geographic origin. These geographic and genetically distinct clades could be referred to as ribotypes: (A) Belize, (B) Florida-Cuba, (C1) India, and (C2) Australia.

A timescale for diatom evolution based on four molecular markers: reassessment of ghost lineages and major steps defining diatom evolution.

ABSTRACT. – Phylogenies and molecular clocks of the diatoms have largely been inferred from SSU rDNA sequences. A new phylogeny of diatoms was estimated using four gene markers SSU and LSU rDNA rbcL and psbA (total 4352 bp) with 42 diatom species. The four gene trees analysed with a maximum likelihood (ML) and Baysian (BI) analysis recovered a monophyletic origin of the new diatom classes with high bootstrap support, which has been controversial with single gene markers using single outgroups and alignments that do not take secondary structure of the SSU gene into account. The divergence time of the classes were calculated from a ML tree in the MultliDiv Time program using a Bayesian estimation allowing for simultaneous constraints from the fossil record and varying rates of molecular evolution of different branches in the phylogenetic tree. These divergence times are generally in agreement with those proposed by other clocks using single genes with the exception that the pennates appear much earlier and suggest a longer Cretaceous fossil record that has yet to be sampled. Ghost lineages (i.e. the discrepancy between first appearance (FA) and molecular clock age of origin from an extant taxon) were revealed in the pennate lineage, whereas those ghost lineages in the centric lineages previously reported by others are reviewed and referred to earlier literature.

A timescale for diatom evolution based on four molecular markers: reassessment of ghost lineages and major steps defining diatom evolution.

ABSTRACT. – Phylogenies and molecular clocks of the diatoms have largely been inferred from SSU rDNA sequences. A new phylogeny of diatoms was estimated using four gene markers SSU and LSU rDNA rbcL and psbA (total 4352 bp) with 42 diatom species. The four gene trees analysed with a maximum likelihood (ML) and Baysian (BI) analysis recovered a monophyletic origin of the new diatom classes with high bootstrap support, which has been controversial with single gene markers using single outgroups and alignments that do not take secondary structure of the SSU gene into account. The divergence time of the classes were calculated from a ML tree in the MultliDiv Time program using a Bayesian estimation allowing for simultaneous constraints from the fossil record and varying rates of molecular evolution of different branches in the phylogenetic tree. These divergence times are generally in agreement with those proposed by other clocks using single genes with the exception that the pennates appear much earlier and suggest a longer Cretaceous fossil record that has yet to be sampled. Ghost lineages (i.e. the discrepancy between first appearance (FA) and molecular clock age of origin from an extant taxon) were revealed in the pennate lineage, whereas those ghost lineages in the centric lineages previously reported by others are reviewed and referred to earlier literature.

Evolution of the Diatoms: VIII. Re-Examination of the SSU-Rrna Gene Using Multiple Outgroups and a Cladistic Analysis of Valve Features.

The resolution of the SSU rRNA gene for phylogenetic analysis in the diatoms has been evaluated by Theriot et al. who claimed that the SSU rRNA gene could not be used to resolve the monophyly of the three diatoms classes described by Medlin and Kaczmarska. Although they used both only bolidomonads and heterokonts as outgroups, they did not explore outgroups further away than the heterokonts. In this study, the use of the multiple outgroups inside and outside the heterokonts with the rRNA gene for recovering the three monophyletic clades at the class level is evaluated. Trees with multiple outgroups ranging from only bolidophytes to Bacteria and Archea were analyzed with Bayesian and Maximum Likelihood analyses and two data sets were recovered with the classes being monophyletic. Other data sets were analyzed with non-weighted and weighted maximum parsimony. The latter reduced the number of clades and lengthened branch lengths between the clades. One data set using a weighted analysis recovered the three classes as monophyletic. Taking only bolidophytes as the only outgroup never produced monophyletic clades. Multiple outgroups including many heterokonts and certain members of the crown group radiation recovered monophyletic clades. The three classes can be defined by clear morphological differences primarily based on auxospore ontogeny and envelope structure, the presence or absence of a structure (tube process or sternum) associated with the annulus and the location of the cribrum in those genera with loculate areolae. A cladistic analysis of some of these features is presented and recovers the three classes.

Unexpected prevalence of parasite 18S rDNA sequences in winter among Antarctic marine protists

Parasites are not typically considered to be important components of polar marine ecosystems. It was therefore surprising when 18S rDNA surveys of protists in the West Antarctic Peninsula in winter revealed high abundances of parasite sequences. Parasite sequences made up, on average, over half (52%) of sequence reads in samples from deep water in winter. Winter surface water and sediment samples contained relatively fewer, but still strikingly high, parasite sequence reads (13 and 9%, respectively), while surface water samples in summer contained fewer parasite sequences (1.8%). A total of 1028 distinct parasite Operational Taxonomic Units were observed in winter, with the largest abundances and diversities within Syndiniales groups I and II, including Amoebophrya. Less abundant parasite sequence groups included Apicomplexa, Blastodinium, Chytriodinium, Cryptocaryon, Paradinium, Perkinsidae, Pirsonia and Ichthyophonae. Parasite sequence distributions suggested interactions with known hosts, such as diatom parasites which were mainly in the sediments, where resting spores of Chaetoceros spp. diatoms were abundant. Syndiniales sequences were correlated with radiolarian sequences, suggesting parasite–host interactions. The abundant proportions of parasite sequences indicate a potentially important role for parasites in the Antarctic marine ecosystem, with implications for plankton population dynamics, the role of the microbial loop, carbon flows and ecosystem responses to ongoing anthropogenic climate change.

Unexpected prevalence of parasite 18S rDNA sequences in winter among Antarctic marine protists

Parasites are not typically considered to be important components of polar marine ecosystems. It was therefore surprising when 18S rDNA surveys of protists in the West Antarctic Peninsula in winter revealed high abundances of parasite sequences. Parasite sequences made up, on average, over half (52%) of sequence reads in samples from deep water in winter. Winter surface water and sediment samples contained relatively fewer, but still strikingly high, parasite sequence reads (13 and 9%, respectively), while surface water samples in summer contained fewer parasite sequences (1.8%). A total of 1028 distinct parasite Operational Taxonomic Units were observed in winter, with the largest abundances and diversities within Syndiniales groups I and II, including Amoebophrya. Less abundant parasite sequence groups included Apicomplexa, Blastodinium, Chytriodinium, Cryptocaryon, Paradinium, Perkinsidae, Pirsonia and Ichthyophonae. Parasite sequence distributions suggested interactions with known hosts, such as diatom parasites which were mainly in the sediments, where resting spores of Chaetoceros spp. diatoms were abundant. Syndiniales sequences were correlated with radiolarian sequences, suggesting parasite–host interactions. The abundant proportions of parasite sequences indicate a potentially important role for parasites in the Antarctic marine ecosystem, with implications for plankton population dynamics, the role of the microbial loop, carbon flows and ecosystem responses to ongoing anthropogenic climate change.

PCR from the CPR offers a historical perspective on marine population ecology.

The Continuous Plankton Recorder (CPR) survey has collected plankton samples from regular tracks across the world's oceans for almost 70 y. Over 299,000 spatially extensive CPR samples are archived and stored in buffered formalin. This CPR archive offers huge potential to study changes in marine communities using molecular data from a period when marine pollution, exploitation and global anthropogenic impact were much less pronounced. However, to harness the amount of data available within the CPR archive fully, it is necessary to improve techniques of larval identification, to genus and species preferably, and to obtain genetic information for historical studies of population ecology. To increase the potential of the CPR database this paper describes the first extraction, amplification by the polymerase chain reaction and utilization of a DNA sequence (mitochondrial 16S rDNA) from a CPR sample, a formalin fixed larval sandeel.

Morphological and genetic variation in the north Atlantic copepod, Centropages typicus

This study describes phenotypic and genotypic variations in the planktonic copepod, Centropages typicus (Copepoda: Calanoida) that indicate differentiation between geographical samples. We found consistent differences in the morphology of the chela of the sexually modified fifth pereiopod (P5) of male C. typicus between samples from the Mediterranean, western North Atlantic and eastern North Atlantic. A 560 base pairs (bp) region of the C. typicus mitochondrial cytochrome c oxidase subunit I (COI) and a 462 bp fragment of the nuclear rDNA internal transcribed spacer (ITS) tandem array were analysed to determine whether these morphological variations reflect population genetic differentiation. Mitochondrial haplotype diversity was found to be high with 100 unique COI haplotypes among 116 individuals. Analysis of mtCOI variation suggested differentiation between the Mediterranean and Atlantic populations but no separation was detected within the Atlantic. Intragenomic variation in the ITS array suggested genetic differentiation between samples from the western North Atlantic and those from the eastern North Atlantic and Mediterranean. Breeding experiments would be required to elucidate the extent of genetic isolation between C. typicus from the different population centres.

Validation of the detection of Pseudo-nitzschia spp. using specific RNA probes tested in a microarray format: Calibration of signal based on variability of RNA content with environmental conditions.

Harmful algal blooms (HAB) occur worldwide and cause health problems and economic damage to fisheries and tourism. Monitoring for toxic algae is therefore essential but is based primarily on light microscopy, which is time consuming and can be limited by insufficient morphological characters such that more time is needed to examine critical features with electron microscopy. Monitoring with molecular tools is done in only a few places world-wide. EU FP7 MIDTAL (Microarray Detection of Toxic Algae) used SSU and LSU rRNA genes as targets on microarrays to identify toxic species. In order to comply with current monitoring requirements to report cell numbers as the relevant threshold measurement to trigger closure of fisheries, it was necessary to calibrate our microarray to convert the hybridisation signal obtained to cell numbers. Calibration curves for two species of Pseudo-nitzschia for use with the MIDTAL microarray are presented to obtain cell numbers following hybridisation. It complements work presented by Barra et al. (2012b. Environ. Sci. Pollut. Res. doi: 10.1007/s11356-012-1330-1v) for two other Pseudo-nitzschia spp., Dittami and Edvardsen (2012a. J. Phycol. 48, 1050) for Pseudochatonella, Blanco et al. (2013. Harmful Algae 24, 80) for Heterosigma, McCoy et al. (2013. FEMS. doi: 10.1111/1574-6941.12277) for Prymnesium spp., Karlodinium veneficum, and cf. Chatonella spp. and Taylor et al. (2014. Harmful Algae, in press) for Alexandrium.

Evolution of the diatoms: major steps in their evolution and a review of the supporting molecular and morphological evidence.

Over the years, many reviews of different aspects of diatom biology, ecology and evolution have appeared. Since 1993 many molecular trees have been produced to infer diatom phylogeny. In 2004, Medlin & Kaczmarska revised the systematics of the diatoms based on more than 20 years of consistent recovery of two major clades of diatoms that did not correspond to a traditional concept of centrics and pennates and established three classes of diatoms: Clade 1 = Coscinodiscophyceae (radial centrics) and Clade 2 = Mediophyceae (polar centrics + radial Thalassiosirales) and Bacillariophyceae (pennates). However, under certain analytical conditions, an alternative view of diatom evolution, a grades of clades, has been recovered that suggests a gradual evolution from centric to pennate symmetry. These two schemes of diatom evolution are evaluated in terms of whether or not the criteria advocated by Medlin & Kaczmarska that should be met to recover monophyletic classes have been used. The monophyly of the three diatom classes can only be achieved if (1) a secondary structure of the small subunit (SSU) rRNA gene was used to construct the alignment and not an alignment based on primary structure and (2) multiple outgroups were used. These requirements have not been met in each study of diatom evolution; hence, the grade of clades, which is useful in reconstructing the sequence of evolution, is not useful for accepting the new classification of the diatoms. Evidence for how these two factors affect the recovery of the three monophyletic classes is reviewed here. The three classes have been defined by clear morphological differences primarily based on gametangia and auxospore ontogeny and envelope structure, the presence or absence of a structure (tube process or sternum) associated with the annulus and the location of the cribrum in those genera with loculate areolae. New evidence supporting the three clades is reviewed. Other features of the cell are examined to determine whether they can also be used to support the monophyly of the three classes.

Evolution of the diatoms: major steps in their evolution and a review of the supporting molecular and morphological evidence.

Over the years, many reviews of different aspects of diatom biology, ecology and evolution have appeared. Since 1993 many molecular trees have been produced to infer diatom phylogeny. In 2004, Medlin & Kaczmarska revised the systematics of the diatoms based on more than 20 years of consistent recovery of two major clades of diatoms that did not correspond to a traditional concept of centrics and pennates and established three classes of diatoms: Clade 1 = Coscinodiscophyceae (radial centrics) and Clade 2 = Mediophyceae (polar centrics + radial Thalassiosirales) and Bacillariophyceae (pennates). However, under certain analytical conditions, an alternative view of diatom evolution, a grades of clades, has been recovered that suggests a gradual evolution from centric to pennate symmetry. These two schemes of diatom evolution are evaluated in terms of whether or not the criteria advocated by Medlin & Kaczmarska that should be met to recover monophyletic classes have been used. The monophyly of the three diatom classes can only be achieved if (1) a secondary structure of the small subunit (SSU) rRNA gene was used to construct the alignment and not an alignment based on primary structure and (2) multiple outgroups were used. These requirements have not been met in each study of diatom evolution; hence, the grade of clades, which is useful in reconstructing the sequence of evolution, is not useful for accepting the new classification of the diatoms. Evidence for how these two factors affect the recovery of the three monophyletic classes is reviewed here. The three classes have been defined by clear morphological differences primarily based on gametangia and auxospore ontogeny and envelope structure, the presence or absence of a structure (tube process or sternum) associated with the annulus and the location of the cribrum in those genera with loculate areolae. New evidence supporting the three clades is reviewed. Other features of the cell are examined to determine whether they can also be used to support the monophyly of the three classes.

Recent Reticulate Evolution in the Ecologically Dominant Lineage of Coccolithophores

The coccolithophore family Noëlaerhabdaceae contains a number of taxa that are very abundant in modern oceans, including the cosmopolitan bloom-forming Emiliania huxleyi. Introgressive hybridization has been suggested to account for incongruences between nuclear, mitochondrial and plastidial phylogenies of morphospecies within this lineage, but the number of species cultured to date remains rather limited. Here, we present the characterization of 5 new Noëlaerhabdaceae culture strains isolated from samples collected in the south-east Pacific Ocean. These were analyzed morphologically using scanning electron microscopy and phylogenetically by sequencing 5 marker genes (nuclear 18S and 28S rDNA, plastidial tufA, and mitochondrial cox1 and cox3 genes). Morphologically, one of these strains corresponded to Gephyrocapsa ericsonii and the four others to Reticulofenestra parvula. Ribosomal gene sequences were near identical between these new strains, but divergent from G. oceanica, G. muellerae, and E. huxleyi. In contrast to the clear distinction in ribosomal phylogenies, sequences from other genomic compartments clustered with those of E. huxleyi strains with which they share an ecological range (i.e., warm temperate to tropical waters). These data provide strong support for the hypothesis of past (and potentially ongoing) introgressive hybridization within this ecologically important lineage and for the transfer of R. parvula to Gephyrocapsa. These results have important implications for understanding the role of hybridization in speciation in vast ocean meta-populations of phytoplankton.

Recent Reticulate Evolution in the Ecologically Dominant Lineage of Coccolithophores

The coccolithophore family Noëlaerhabdaceae contains a number of taxa that are very abundant in modern oceans, including the cosmopolitan bloom-forming Emiliania huxleyi. Introgressive hybridization has been suggested to account for incongruences between nuclear, mitochondrial and plastidial phylogenies of morphospecies within this lineage, but the number of species cultured to date remains rather limited. Here, we present the characterization of 5 new Noëlaerhabdaceae culture strains isolated from samples collected in the south-east Pacific Ocean. These were analyzed morphologically using scanning electron microscopy and phylogenetically by sequencing 5 marker genes (nuclear 18S and 28S rDNA, plastidial tufA, and mitochondrial cox1 and cox3 genes). Morphologically, one of these strains corresponded to Gephyrocapsa ericsonii and the four others to Reticulofenestra parvula. Ribosomal gene sequences were near identical between these new strains, but divergent from G. oceanica, G. muellerae, and E. huxleyi. In contrast to the clear distinction in ribosomal phylogenies, sequences from other genomic compartments clustered with those of E. huxleyi strains with which they share an ecological range (i.e., warm temperate to tropical waters). These data provide strong support for the hypothesis of past (and potentially ongoing) introgressive hybridization within this ecologically important lineage and for the transfer of R. parvula to Gephyrocapsa. These results have important implications for understanding the role of hybridization in speciation in vast ocean meta-populations of phytoplankton.

Morphological and Phylogenetic Characterization of New Gephyrocapsa Isolates Suggests Introgressive Hybridization in the Emiliania/Gephyrocapsa Complex (Haptophyta)

The coccolithophore genus Gephyrocapsa contains a cosmopolitan assemblage of pelagic species, including the bloom-forming Gephyrocapsa oceanica, and is closely related to the emblematic coccolithophore Emiliania huxleyi within the Noëlaerhabdaceae. These two species have been extensively studied and are well represented in culture collections, whereas cultures of other species of this family are lacking. We report on three new strains of Gephyrocapsa isolated into culture from samples from the Chilean coastal upwelling zone using a novel flow cytometric single-cell sorting technique. The strains were characterized by morphological analysis using scanning electron microscopy and phylogenetic analysis of 6 genes (nuclear 18S and 28S rDNA, plastidial 16S and tufA, and mitochondrial cox1 and cox3 genes). Morphometric features of the coccoliths indicate that these isolates are distinct from G. oceanica and best correspond to G. muellerae. Surprisingly, both plastidial and mitochondrial gene phylogenies placed these strains within the E. huxleyi clade and well separated from G. oceanica isolates, making Emiliania appear polyphyletic. The only nuclear sequence difference, 1 bp in the 28S rDNA region, also grouped E. huxleyi with the new Gephyrocapsa isolates and apart from G. oceanica. Specifically, the G. muellerae morphotype strains clustered with the mitochondrial β clade of E. huxleyi, which, like G. muellerae, has been associated with cold (temperate and sub-polar) waters. Among putative evolutionary scenarios that could explain these results we discuss the possibility that E. huxleyi is not a valid taxonomic unit, or, alternatively the possibility of past hybridization and introgression between each E. huxleyi clade and older Gephyrocapsa clades. In either case, the results support the transfer of Emiliania to Gephyrocapsa. These results have important implications for relating morphological species concepts to ecological and evolutionary units of diversity.