Search

Although most of the tens of thousands of diatom species are obligate photoautotrophs, many mixotrophic species can also use extracellular organic carbon for growth, and a small number of obligate heterotrophs have lost photosynthesis entirely. We sequenced the genome of a nonphotosynthetic diatom, Nitzschia sp. strain Nitz4, to determine how carbon metabolism was altered in the wake of this rare and radical trophic shift in diatoms. Like other groups that have lost photosynthesis, the genomic consequences were most evident in the plastid genome, which is exceptionally AT-rich and missing photosynthesis-related genes. The relatively small (27 Mb) nuclear genome did not differ dramatically from photosynthetic diatoms in gene or intron density. Genome-based models suggest that central carbon metabolism, including a central role for the plastid, remains relatively intact in the absence of photosynthesis. All diatom plastids lack an oxidative pentose phosphate pathway (PPP), leaving photosynthesis as the main source of plastid NADPH. Consequently, nonphotosynthetic diatoms lack the primary source of NADPH required for essential biosynthetic pathways that remain in the plastid. Genomic models highlighted similarities between nonphotosynthetic diatoms and apicomplexan parasites for provisioning NADPH in their plastids. The ancestral absence of a plastid PPP might constrain loss of photosynthesis in diatoms compared to Archaeplastida, whose plastid PPP continues to produce reducing cofactors following loss of photosynthesis. Finally, Nitzschia possesses a complete β-ketoadipate pathway. Previously known only from fungi and bacteria, this pathway may allow mixotrophic and heterotrophic diatoms to obtain energy through the degradation of abundant plant-derived aromatic compounds.

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic dataset of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a small set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

The large population sizes and high dispersal potential of microbes suggests that a given microbial species should be found in all suitable habitats worldwide. Consequently, microbes should not exhibit the kinds of biogeographic patterns seen in macroorganisms. This paradigm is challenged by a growing list of exotic microbes with biogeographic disjunctions that instead promotes microbial dispersal as inherently limited. We sampled water bodies in the United States and compiled records from the literature and public databases to characterize the distribution of the freshwater planktonic diatom, Discostella asterocostata (Xie, Lin, and Cai) Houk and Klee. Discostella asterocostata was thought to be restricted to the Far East, but we report its presence in ecologically similar water bodies across the eastern United States. Populations from the U.S. and China are indistinguishable morphometrically, suggesting they may be recently separated—a hypothesis supported by paleolimnological data, which support an introduction of D. asterocostata into the U.S. as recently as the mid-1980s. The overlapping distributions of D. asterocostata and invasive carp species, in both their native and nonnative ranges, highlighted Asian carp as a possible vector for introduction of the diatom in the U.S. The existence of exotic diatoms underscores natural constraints on microbial dispersal, resulting in biogeographic distributions that can be upended through human activity.

The diatom, Cyclotella cryptica, is a well-established model species for physiological studies and biotechnology applications of diatoms. To further facilitate its use as a model diatom, we report an improved reference genome assembly and annotation for C. cryptica strain CCMP332. We used a combination of long- and short-read sequencing to assemble a high-quality and contaminant-free genome. The genome is 171 Mb in size and consists of 662 scaffolds with a scaffold N50 of 494 kb. This represents a 176-fold decrease in scaffold number and 41-fold increase in scaffold N50 compared to the previous assembly. The genome contains 21,250 predicted genes, 75% of which were assigned putative functions. Repetitive DNA comprises 59% of the genome, and an improved classification of repetitive elements indicated that a historically steady accumulation of transposable elements has contributed to the relatively large size of the C. cryptica genome. The high-quality C. cryptica genome will serve as a valuable reference for ecological, genetic, and biotechnology studies of diatoms.

The salinity gradient separating marine and freshwater environments represents a major ecological divide for microbiota, yet the mechanisms by which marine microbes have adapted to and ultimately diversified in freshwater environments are poorly understood. Here, we take advantage of a natural evolutionary experiment: the colonization of the brackish Baltic Sea by the ancestrally marine diatom Skeletonema marinoi. To understand how diatoms respond to low salinity, we characterized transcriptomic responses of S. marinoi grown in a common garden. Our experiment included eight genotypes from source populations spanning the Baltic Sea salinity cline. Changes in gene expression revealed a shared response to salinity across genotypes, where low salinities induced profound changes in cellular metabolism, including upregulation of carbon fixation and storage compound biosynthesis, and increased nutrient demand and oxidative stress. Nevertheless, the genotype effect overshadowed the salinity effect, as genotypes differed significantly in their response, both in the magnitude and direction of gene expression. Intraspecific differences included regulation of transcription and translation, nitrogen metabolism, cell signaling, and aerobic respiration. The high degree of intraspecific variation in gene expression observed here highlights an important but often overlooked source of biological variation associated with how diatoms respond and adapt to environmental change.

The diatom,Cyclotella cryptica, is a well-established experimental model for physiological studies and, more recently, biotechnology applications of diatoms. To further facilitate its use as a model diatom species, we report an improved reference genome assembly and annotation forC. crypticastrain CCMP332. We used a combination of long- and short-read sequencing to assemble a high-quality and contaminant-free genome. The genome is 171 Mb in size and consists of 662 scaffolds with a scaffold N50 of 494 kb. This represents a 176-fold decrease in scaffold number and 41-fold increase in scaffold N50 compared to the previous assembly. The genome contains 21,250 predicted genes, 75% of which were assigned putative functions. Repetitive DNA comprises 59% of the genome, and an improved classification of repetitive elements indicated that a historically steady accumulation of transposable elements has contributed to the relatively large size of theC. crypticagenome. The high-qualityC. crypticagenome will serve as a valuable reference for ecological, genetic, and biotechnology studies of diatoms.Data available fromNCBI BioProjects PRJNA628076 and PRJNA589195

We sequenced mitochondrial genomes from five diverse diatoms (Toxarium undulatum, Psammoneis japonica, Eunotia naegelii, Cylindrotheca closterium, and Nitzschia sp.), chosen to fill important phylogenetic gaps and help us characterize broadscale patterns of mitochondrial genome evolution in diatoms. Although gene content was strongly conserved, intron content varied widely across species. The vast majority of introns were of group II type and were located in the cox1 or rnl genes. Although recurrent intron loss appears to be the principal underlying cause of the sporadic distributions of mitochondrial introns across diatoms, phylogenetic analyses showed that intron distributions superficially consistent with a recurrent-loss model were sometimes more complicated, implicating horizontal transfer as a likely mechanism of intron acquisition as well. It was not clear, however, whether diatoms were the donors or recipients of horizontally transferred introns, highlighting a general challenge in resolving the evolutionary histories of many diatom mitochondrial introns. Although some of these histories may become clearer as more genomes are sampled, high rates of intron loss suggest that the origins of many diatom mitochondrial introns are likely to remain unclear.

How diatoms respond to fluctuations in osmotic pressure is important from both ecological and applied perspectives. It is well known that osmotic stress affects photosynthesis and can result in the accumulation of compounds desirable in pharmaceutical and alternative fuel industries. Gene expression responses to osmotic stress have been studied in short-term trials, but it is unclear whether the same mechanisms are recruited during long-term acclimation. We used RNA-seq to study the genome-wide transcription patterns in the euryhaline diatom, Cyclotella cryptica, following long-term acclimation to salinity that spanned the natural range of fresh to oceanic water. Long-term acclimated C. cryptica exhibited induced synthesis or repressed degradation of the osmolytes glycine betaine, taurine and dimethylsulfoniopropionate (DMSP). Although changes in proline concentration is one of the main responses in short-term osmotic stress, we did not detect a transcriptional change in proline biosynthetic pathways in our long-term experiment. Expression of membrane transporters showed a general tendency for increased import of potassium and export of sodium, consistent with the electrochemical gradients and dependence on co-transported molecules. Our results show substantial between-genotype differences in growth and gene expression reaction norms and suggest that the regulation of proline synthesis important in short-term osmotic stress might not be maintained in long-term acclimation. Further examination using time-course gene expression experiments, metabolomics and genetic validation of gene functions would reinforce patterns inferred from RNA-seq data.

PREMISE OF THE STUDY Loss of photosynthesis is a common and often repeated trajectory in nearly all major groups of photosynthetic eukaryotes. One small subset of "apochloritic" diatoms in the genus Nitzschia have lost their ability to photosynthesize and require extracellular carbon for growth. Similar to other secondarily nonphotosynthetic taxa, apochloritic diatoms maintain colorless plastids with highly reduced plastid genomes. Although the narrow taxonomic breadth of apochloritic Nitzschia suggests a single loss of photosynthesis in their common ancestor, previous phylogenetic analyses suggested that photosynthesis was lost multiple times. METHODS We analyzed genes from the nuclear, plastid, and mitochondrial genomes for a broad set of taxa to test whether photosynthesis was lost one or multiple times in Bacillariales. We also sequenced and characterized the plastid genome of a nonphotosynthetic Nitzschia species. KEY RESULTS Phylogenetic analyses showed that genes from all three genetic compartments either supported or failed to reject monophyly of apochloritic Nitzschia species, consistent with a single loss of photosynthesis in this group. The plastid genomes of two apochloritic Nitzschia are highly similar in all respects, indicating streamlining of the plastid genome before the split of these two species. CONCLUSIONS A better understanding of the phylogeny and ecology of apochloritic Nitzschia, together with emerging genomic resources, will help identify the factors that have driven and maintained the loss of photosynthesis in this group of diatoms. Finally, some habitats host diverse communities of co-occurring nonphotosynthetic diatoms, reflecting resource abundance or resource partitioning in ecologically favorable habitats.

Loss of photosynthesis is a common and often repeated trajectory in nearly all major groups of photosynthetic eukaryotes. One small subset of ‘apochloritic’ diatoms in the genus Nitzschia have lost their ability to photosynthesize and require extracellular carbon for growth. Similar to other secondarily nonphotosynthetic taxa, apochloritic diatoms maintain colorless plastids with highly reduced plastid genomes. Although the narrow taxonomic breadth of apochloritic diatoms suggests a single loss of photosynthesis in the common ancestor of these species, previous phylogenetic analyses suggested that photosynthesis was lost multiple times. We sequenced additional phylogenetic markers from the nuclear and mitochondrial genomes for a larger set of taxa and found that the best trees for datasets representing all three genetic compartments provided low to moderate support for monophyly of apochloritic Nitzschia, consistent with a single loss of photosynthesis in diatoms. We sequenced the plastid genome of one apochloritic species and found that it was highly similar in all respects to the plastid genome of another apochloritic Nitzschia species, indicating that streamlining of the plastid genome had completed prior to the split of these two species. Finally, it is increasingly clear that some locales host relatively large numbers apochloritic Nitzschia species that span the phylogenetic diversity of the group, indicating that these species co-exist because of resource abundance or resource partitioning in ecologically favorable habitats. A better understanding of the phylogeny and ecology of this group, together with emerging genomic resources, will help identify the factors that have driven and maintained the loss of photosynthesis in this group, a rare event in diatoms.

Premise of the studyDiatoms are one of the most species-rich lineages of microbial eukaryotes. Similarities in clade age, species richness, and contributions to primary production motivate comparisons to flowering plants, whose genomes have been inordinately shaped by whole genome duplication (WGD). These events that have been linked to speciation and increased rates of lineage diversification, identifying WGDs as a principal driver of angiosperm evolution. We synthesized a relatively large but scattered body of evidence that, taken together, suggests that polyploidy may be common in diatoms.MethodsWe used data from gene counts, gene trees, and patterns of synonymous divergence to carry out the first large-scale phylogenomic analysis of genome-scale duplication histories for a phylogenetically diverse set of 37 diatom taxa.Key resultsSeveral methods identified WGD events of varying age across diatoms, though determining the exact number and placement of events and, more broadly, inferences of WGD at all, were greatly impacted by gene-tree uncertainty. Gene-tree reconciliations supported allopolyploidy as the predominant mode of polyploid formation, with particularly strong evidence for ancient allopolyploid events in the thalassiosiroid and pennate diatom clades.ConclusionsWhole genome duplication appears to have been an important driver of genome evolution in diatoms. Denser taxon sampling will better pinpoint the timing of WGDs and likely reveal many more of them. We outline potential challenges in reconstructing paleopolyploid events in diatoms that, together with these results, offer a framework for understanding the evolutionary roles of genome duplication in a group that likely harbors substantial genomic diversity.

The diatom genus Entomoneis is known from the benthos and plankton of marine, brackish, and freshwaters. Entomoneis includes diatoms with a bilobate keel elevated above the valve surface, a sigmoid canal raphe, and numerous girdle bands. Owing mostly to the scarcity of molecular data for a diverse set of species, the phylogeny of Entomoneis has not been investigated in depth. The few previous studies that included Entomoneis were focused on broader questions and the available data were from a small number of either unidentified Entomoneis or well-known species (e.g., E. paludosa). Since the first description of new species combining both molecular and morphological characters (E. tenera), we have continued to cultivate and investigate Entomoneis in the plankton of the Adriatic Sea. Combined multigene phylogeny (SSU rDNA sequences, rbcL, and psbC genes) and morphological observations (LM, SEM and TEM) revealed six new Entomoneis species supported by phylogenetic and morphological data: E. pusilla, E. gracilis, E. vilicicii, E. infula, E. adriatica, and E. umbratica. The most important morphological features for species delineation were cell shape, the degree and mode of torsion, valve apices, the appearance and structure of the transition between keel and valve body, the ultrastructure and the shape of the girdle bands, and the arrangement and density of perforations along the valve and valvocopulae. Our results highlight the underappreciated diversity of Entomoneis and call for a more in-depth morphological and molecular investigation of this genus especially in planktonic habitats.

We reconstructed the phylogeny of representatives from nine genera and three families of the Cymbellales using two nuclear and three chloroplast genes. After rooting with Anomoeoneis, Placoneis was...

Photosynthesis by diatoms accounts for roughly one-fifth of global primary production, but despite this, relatively little is known about their plastid genomes. We report the completely sequenced plastid genomes for eight phylogenetically diverse diatoms and show them to be variable in size, gene and foreign sequence content, and gene order. The genomes contain a core set of 122 protein-coding genes, with 15 additional genes exhibiting complex patterns of 1) gene losses at varying phylogenetic scales, 2) functional transfers to the nucleus, 3) gene duplication, divergence, and differential retention of paralogs, and 4) acquisitions of putatively functional recombinase genes from resident plasmids. The newly sequenced genomes also contain several previously unreported genes, highlighting how poorly characterized diatom plastid genomes are overall. Genome size variation reflects major expansions of the inverted repeat region in some cases but, more commonly, large-scale expansions of intergenic regions, many of which contain unique open reading frames of likely foreign origin. Although many gene clusters are conserved across species, rearrangements appear to be frequent in most lineages.

Although the plastid genomes of diatoms maintain a conserved architecture and core gene set, considerable variation about this core theme exists and can be traced to several different processes. Gene duplication, pseudogenization, and loss, as well as intracellular transfer of genes to the nuclear genome, have all contributed to variation in gene content among diatom species. In addition, some noncoding sequences have highly restricted phylogenetic distributions that suggest a recent foreign origin. We sequenced the plastid genome of the marine diatom, Toxarium undulatum, and found that the genome contains three genes (chlB, chlL, and chlN) involved in light-independent chlorophyll a biosynthesis that were not previously known from diatoms. Phylogenetic and syntenic data suggest that these genes were differentially retained in this one lineage as they were repeatedly lost from most other diatoms. Unique among diatoms and other heterokont algae sequenced so far, the genome also contains a large group II intron within an otherwise intact psaA gene. Although the intron is most similar to one in the plastid-encoded psaA gene of some green algae, high sequence divergence between the diatom and green algal introns rules out recent shared ancestry. We conclude that the psaA intron was likely introduced into the plastid genome of T. undulatum, or some earlier ancestor, by horizontal transfer from an unknown donor. This genome further highlights the myriad processes driving variation in gene and intron content in the plastid genomes of diatoms, one of the world's foremost primary producers.

The Surirellales and Rhopalodiales are large, widespread, and morphologically diverse groups of raphid pennate diatoms (Bacillariphyta) whose raphe, a structure that facilitates active motility, opens internally into a siliceous canal. We collected 202 representatives of the lineage and sequenced genes from the nuclear, plastid, and mitochondrial genomes to infer phylogenetic relationships as a basis for comparative study of ecology and morphological evolution as well as reclassification. The lineage was ancestrally marine, and we report the first evidence for a ‘stepping stone’ model of marine–freshwater transitions in which freshwater invasions were preceded by adaptation to intermediate brackish habitats. Phylogenetic comparative analyses also showed that the shift from an apical (e.g., Entomoneis) to transapical major axis of development (e.g., Surirella) did not have to proceed through subcircular intermediate forms (i.e., Campylodiscus). Rather, subcircular forms evolved both within lineages with longer apical axis or longer transapical axis. We also used the inferred phylogeny as a basis for genus-level reclassification of the lineage. Campylodiscus now includes the fastuosoid members of Surirella and Campylodiscus, but excludes other marine Campylodiscus which are now classified as Coronia. Surirella includes the Surirella striatula clade, Surirella Pinnatae group, and species formerly classified as Cymatopleura. We resurrected the genus Iconella to accommodate Stenopterobia and the robustoid members of Surirella and Campylodiscus. We broadened Epithemia to include members of the paraphyletic genus Rhopalodia. Finally, we discuss the challenges of constructing a classification that best leverages available phylogenetic data, while minimizing disruption to the research community and recognizing practical considerations stemming from the slow rate of progress on systematic studies of understudied organisms.

Lobban C.S., Schefter M. and Ruck E.C. 2011. Licmophora flucticulata sp. nov. (Licmophoraceae, Bacillariophyceae), an unusual new flabellate species from Guam and Palau. Phycologia 50: 11–22. DOI: 10.2216/09-85.1 A new species of Licmophora, L. flucticulata sp. nov., forms distinctive, rippled, fan-shaped colonies on coral reef seaweeds on Western Pacific islands. It is readily observed underwater with the naked eye but is very weakly silicified, and its frustule does not survive normal acid cleaning. Its identity as a Licmophora species was demonstrated with scanning electron microscopy and small-subunit rDNA sequencing. Licmophora flucticulata has exceptionally long, narrow cells cemented into fascicles that attach to the substratum by short, multistranded mucilage stalks. Its valve and colony morphology are compared to that of L. remulus, L. flabellata and L. aurivillii.

1Texas Natural Science Center, 2400 Trinity Street, University of Texas, Austin TX 78619, U.S.A. 2Section of Integrative Biology, The University of Texas at Austin, 1 University Station (A6700), Austin, TX 78712, U.S.A. 3Plant Biology Graduate Program, The University of Texas at Austin, 1 University Station (A6700), Austin, TX 78712, U.S.A. 4Institute of Cellular and Molecular Biology, The University of Texas at Austin, 1 University Station (A6700), Austin, TX 78712, U.S.A. *Author for correspondence: etheriot@austin.utexas.edu

Previous analyses of single diatom chloroplast protein-encoded genes recovered results highly incongruent with both traditional phylogenies and phylogenies derived from the nuclear encoded small subunit (SSU) gene. Our analysis here of six individual chloroplast genes (atpB, psaA, psaB, psbA, psbC and rbcL) obtained similar anomalous results. However, phylogenetic noise in these genes did not appear to be correlated, and their concatenation appeared to effectively sum their collective signal. We empirically demonstrated the value of combining phylogenetic information profiling, partitioned Bremer support and entropy analysis in examining the utility of various partitions in phylogenetic analysis. Noise was low in the 1st and 2nd codon positions, but so was signal. Conversely, high noise levels in the 3rd codon position was accompanied by high signal. Perhaps counterintuitively, simple exclusion experiments demonstrated this was especially true at deeper nodes where the 3rd codon position contributed most to a result congruent with morphology and SSU (and the total evidence tree here). Correlated with our empirical findings, probability of correct signal (derived from information profiling) increased and the statistical significance of substitutional saturation decreased as data were aggregated. In this regard, the aggregated 3rd codon position performed as well or better than more slowly evolving sites. Simply put, direct methods of noise removal (elimination of fast-evolving sites) disproportionately removed signal. Information profiling and partitioned Bremer support suggest that addition of chloroplast data will rapidly improve our understanding of the diatom phylogeny, but conversely also illustrate that some parts of the diatom tree are likely to remain recalcitrant to addition of molecular data. The methods based on information profiling have been criticized for their numerous assumptions and parameter estimates and the fact that they are based on quartets of taxa. Our empirical results support theoretical arguments that the simplifying assumptions made in these methods are robust to “real-life” situations.

Diatoms are perhaps the most diverse lineage of eukaryotic algae, with their siliceous cell wall and diplontic life history often considered to have played important roles in their extraordinary diversification. The characteristic diminution of the diatom cell wall over the course of vegetative growth provides a reliable, intrinsic trigger for sexual reproduction, establishing a direct link between the evolution of their cell-wall and life-history features. It is unclear, however, whether the diplontic life cycle of diatoms represents an ancestral or derived trait. This uncertainty is based in part on our lack of understanding of the life cycle of the sister lineage to diatoms, which includes a mix of two free-living and separately classified forms: naked biflagellate unicells in the genus Bolidomonas and silicified forms in the order Parmales. These two forms might represent different life-history stages, although directly establishing such links can be difficult. We sequenced transcriptomes for Bolidomonas and two diatoms and found that ~0.1% of the coding regions in the two diploid diatoms are heterozygous, whereas Bolidomonas is virtually devoid of heterozygous alleles, consistent with expectations for a haploid genome. These results suggest that Bolidomonas is haploid and predict that parmaleans represent the diploid phase of a haplodiplontic life cycle. These data fill an important gap in our understanding of the origin of the diplontic life history of diatoms, which may represent an evolutionarily derived, adaptive feature.

Diatoms are often referred to one of six structural groups. The two major groups are centrics and pennates, and each is further subdivided. Centrics are either radial centrics or (bi-)multipolar (or simply polar) centrics. The former typically are circular and lack any prominent structures which may be paired or multiply arranged so as to give some sort of visually prominent polarity to the cell. Polar centrics have such structures and often have elongate outlines. Pennates are either araphid pennates or raphid pennates, depending on whether or not they possess a raphe. These structural groups have been arranged differently through time, whether the source of data was morphology and the method of analysis noncanonical or whether the data were molecular and analyzed by more formal methods. Both congruence and conflict between these various approaches have been claimed. Diatomists have rejected traditional views because they conflict with molecular results in some instances, and yet reject molecular results because they conflict with morphologically based results in others. Such conflicts are rarely formally tested. Here, we formally test several traditional hypotheses and a recent molecular-based reclassification of diatoms against a three-gene combined molecular dataset. Centrics are strongly rejected as monophyletic. However, some relationships could not be rejected. Monophyly of araphids is not statistically worse than the best tree (in which araphids are recovered as a grade). Monophyly of radial centrics and of polar centrics cannot be rejected, nor can a competing hypothesis (in which radial centrics are a grade and the Thalassiosirales are part of that grade). This last result is congruent with complex morphological characters and is an example of the value of formally testing conflict and congruence between datasets, and of the potential value of formal phylogenetic analysis of diatom morphology.

One lineage of pennate diatoms has a slit through the siliceous cell wall, called a “raphe,” that functions in motility. Raphid pennate diatoms number in the perhaps tens of thousands of species, with the diversity of raphe forms potentially matching this number. Three lineages—the Bacillariales, Rhopalodiales, and Surirellales—possess a complex and presumably highly derived raphe that is physically separated from the cell interior, most often by a set of siliceous braces. Because the relationship among these three lineages is unclear, the number of origins of the canal raphe system and the homology of it and its constitutive parts among these lineages, is equally unclear. We reconstructed the phylogeny of raphid pennate diatoms and included, for the first time, members of all three canal raphid diatom lineages, and used the phylogeny to test specific hypotheses about the origin of the canal raphe. The canal raphe appears to have evolved twice, once in the common ancestor of Bacillariales and once in the common ancestor of Rhopalodiales and Surirellales, which form a monophyletic group in our analyses. These results recommend careful follow-up morphogenesis studies of the canal raphe in these two lineages to determine the underlying developmental basis for this remarkable case of parallel evolution.

The plastid genome of Trifolium subterraneum is 144,763 bp, about 20 kb longer than those of closely related legumes, which also lost one copy of the large inverted repeat (IR). The genome has undergone extensive genomic reconfiguration, including the loss of six genes (accD, infA, rpl22, rps16, rps18, and ycf1) and two introns (clpP and rps12) and numerous gene order changes, attributable to 14-18 inversions. All endpoints of rearranged gene clusters are flanked by repeated sequences, tRNAs, or pseudogenes. One unusual feature of the Trifolium subterraneum genome is the large number of dispersed repeats, which comprise 19.5% (ca. 28 kb) of the genome (versus about 4% for other angiosperms) and account for part of the increase in genome size. Nine genes (psbT, rbcL, clpP, rps3, rpl23, atpB, psbN, trnI-cau, and ycf3) have also been duplicated either partially or completely. rpl23 is the most highly duplicated gene, with portions of this gene duplicated six times. Comparisons of the Trifolium plastid genome with the Plant Repeat Database and searches for flanking inverted repeats suggest that the high incidence of dispersed repeats and rearrangements is not likely the result of transposition. Trifolium has 19.5 kb of unique DNA distributed among 160 fragments ranging in size from 30 to 494 bp, greatly surpassing the other five sequenced legume plastid genomes in novel DNA content. At least some of this unique DNA may represent horizontal transfer from bacterial genomes. These unusual features provide direction for the development of more complex models of plastid genome evolution.