Search

Environmental transitions, such as the salinity divide separating marine and fresh waters, shape biodiversity over both shallow and deep timescales, opening up new niches and creating opportunities for accelerated speciation and adaptive radiation. Understanding the genetics of environmental adaptation is central to understanding how organisms colonise and subsequently diversify in new habitats. We used time-resolved transcriptomics to contrast the hyposalinity stress responses of two diatoms. Skeletonema marinoi has deep marine ancestry but has recently invaded brackish waters. Cyclotella cryptica has deep freshwater ancestry and can withstand a much broader salinity range. Skeletonema marinoi is less adept at mitigating even mild salinity stress compared to Cyclotella cryptica, which has distinct mechanisms for rapid mitigation of hyposaline stress and long-term growth in low salinity. We show that the cellular mechanisms underlying low salinity tolerance, which has allowed diversification across freshwater habitats worldwide, includes elements that are both conserved and variable across the diatom lineage. The balance between ancestral and lineage-specific environmental responses in phytoplankton have shaped marine-freshwater transitions on evolutionary timescales and, on contemporary timescales, will affect which lineages survive and adapt to changing ocean conditions., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

The Cretaceous period is the time of the first appearance of the diatoms in the fossil record. These fossils give us direct evidence of the age and early evolution of the diatom lineage. The fossil record, however, is incomplete and therefore often extrapolated through time-calibrated phylogenies. These two approaches offer different perspectives on the early evolution of diatoms, which is still poorly understood. We compiled the first comprehensive Cretaceous Diatom Database, a tool to investigate the taxonomy, diversity, and occurrence of the earliest known diatom lineages. To further aid the integration and use of the oldest diatom fossils in molecular clock analyses, we present a set of well-documented Cretaceous fossils that can be placed onto molecular phylogenetic trees of extant and extinct species, making them ideal candidates for the calibration of molecular clocks. The analysis of the fossil record and the Cretaceous Diatom Database revealed Cretaceous diversity is substantially greater than previously thought, yet considerable taxonomic work is still needed. The Cretaceous Diatom Database and the list of Cretaceous fossils for calibrating molecular clocks represent valuable resources for future evolutionary and taxonomic studies of modern and fossil diatoms., (© 2024 The Author(s). Journal of Phycology published by Wiley Periodicals LLC on behalf of Phycological Society of America.)

A principal goal in ecology is to identify the determinants of species abundances in nature. Body size has emerged as a fundamental and repeatable predictor of abundance, with smaller organisms occurring in greater numbers than larger ones. A biogeographic component, known as Bergmann's rule, describes the preponderance, across taxonomic groups, of larger-bodied organisms in colder areas. Although undeniably important, the extent to which body size is the key trait underlying these patterns is unclear. We explored these questions in diatoms, unicellular algae of global importance for their roles in carbon fixation and energy flow through marine food webs. Using a phylogenomic dataset from a single lineage with worldwide distribution, we found that body size (cell volume) was strongly correlated with genome size, which varied by 50-fold across species and was driven by differences in the amount of repetitive DNA. However, directional models identified temperature and genome size, not cell size, as having the greatest influence on maximum population growth rate. A global metabarcoding dataset further identified genome size as a strong predictor of species abundance in the ocean, but only in colder regions at high and low latitudes where diatoms with large genomes dominated, a pattern consistent with Bergmann's rule. Although species abundances are shaped by myriad interacting abiotic and biotic factors, genome size alone was a remarkably strong predictor of abundance. Taken together, these results highlight the cascading cellular and ecological consequences of macroevolutionary changes in an emergent trait, genome size, one of the most fundamental and irreducible properties of an organism., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Roberts et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Numerous factors shape the evolution of protein-coding genes, including shifts in the strength or type of selection following gene duplications or changes in the environment. Diatoms and other silicifying organisms use a family of silicon transporters (SITs) to import dissolved silicon from the environment. Freshwaters contain higher silicon levels than oceans, and marine diatoms have more efficient uptake kinetics and less silicon in their cell walls, making them better competitors for a scarce resource. We compiled SITs from 37 diatom genomes to characterize shifts in selection following gene duplications and marine-freshwater transitions. A deep gene duplication, which coincided with a whole-genome duplication, gave rise to two gene lineages. One of them (SIT1-2) is present in multiple copies in most species and is known to actively import silicon. These SITs have evolved under strong purifying selection that was relaxed in freshwater taxa. Episodic diversifying selection was detected but not associated with gene duplications or habitat shifts. In contrast, genes in the second SIT lineage (SIT3) were present in just half the species, the result of multiple losses. Despite conservation of SIT3 in some lineages for the past 90-100 million years, repeated losses, relaxed selection, and low expression highlighted the dispensability of SIT3, consistent with a model of deterioration and eventual loss due to relaxed selection on SIT3 expression. The extensive but relatively balanced history of duplications and losses, together with paralog-specific expression patterns, suggest diatoms continuously balance gene dosage and expression dynamics to optimize silicon transport across major environmental gradients., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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 data set 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 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 Author(s) 2023. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Diatoms are ancestrally photosynthetic microalgae. However, some underwent a major evolutionary transition, losing photosynthesis to become obligate heterotrophs. The molecular and physiological basis for this transition is unclear. Here, we isolate and characterize new strains of non-photosynthetic diatoms from the coastal waters of Singapore. These diatoms occupy diverse ecological niches and display glucose-mediated catabolite repression, a classical feature of bacterial and fungal heterotrophs. Live-cell imaging reveals deposition of secreted extracellular polymeric substance (EPS). Diatoms moving on pre-existing EPS trails (runners) move faster than those laying new trails (blazers). This leads to cell-to-cell coupling where runners can push blazers to make them move faster. Calibrated micropipettes measure substantial single-cell pushing forces, which are consistent with high-order myosin motor cooperativity. Collisions that impede forward motion induce reversal, revealing navigation-related force sensing. Together, these data identify aspects of metabolism and motility that are likely to promote and underpin diatom heterotrophy.

The salinity gradient separating marine and freshwater environments is a major ecological divide, and the mechanisms by which most organisms adapt to new salinity environments are poorly understood. Diatoms are a lineage of ancestrally marine microalgae that have repeatedly colonized and diversified in freshwaters. Cyclotella cryptica is a euryhaline diatom found in salinities ranging from fully freshwater to fully marine, thus providing a powerful system for understanding the genomic mechanisms for mitigating and acclimating to low salinity. To understand how diatoms mitigate acute hypo-osmotic stress, we abruptly shifted C. cryptica from seawater to freshwater and performed transcriptional profiling during the first 10 h. Freshwater shock dramatically remodelled the transcriptome, with ~50% of the genome differentially expressed in at least one time point. The peak response occurred within 1 h, with strong repression of genes involved in cell growth and osmolyte production, and strong induction of specific stress defence genes. Transcripts largely returned to baseline levels within 4-10 h, with growth resuming shortly thereafter, suggesting that gene expression dynamics may be useful for predicting acclimation. Moreover, comparison to a transcriptomics study of C. cryptica following months-long acclimation to freshwater revealed little overlap between the genes and processes differentially expressed in cells exposed to acute stress versus fully acclimated conditions. Altogether, this study highlights the power of time-resolved transcriptomics to reveal fundamental insights into how cells dynamically respond to an acute environmental shift and provides new insights into how diatoms mitigate natural salinity fluctuations and have successfully diversified across freshwater habitats worldwide., (© 2022 John Wiley & Sons Ltd.)

Molecular clocks estimate that diatom microalgae, one of Earth's foremost primary producers, originated near the Triassic-Jurassic boundary (200 Ma), which is close in age to the earliest, generally accepted diatom fossils of the genus Pyxidicula. During an extensive search for Jurassic diatoms from twenty-five sites worldwide, three sites yielded microfossils initially recognized as diatoms. After applying stringent safeguards and evaluation criteria, however, the fossils found at each of the three sites were rejected as new diatom records. This led us to systematically reexamine published evidence in support of Lower- and Middle-Jurassic Pyxidicula fossils. Although Pyxidicula resembles some extant radial centric diatoms and has character states that may have been similar to those of ancestral diatoms, we describe numerous sources of uncertainty regarding the reliability of these records. We conclude that the Lower Jurassic Pyxidicula fossils were most likely calcareous nannofossils, whereas the Middle Jurassic Pyxidicula species has been reassigned to the Lower Cretaceous and is likely a testate amoeba, not a diatom. Excluding the Pyxidicula fossils widens the gap between the estimated time of origin and the oldest abundant fossil diatom record to 75 million years. This study underscores the difficulties in discovering and validating ancient microfossils., (© 2023. The Author(s).)

In contrast to surveys based on a few genes that often provide limited taxonomic resolution, transcriptomes provide a wealth of genomic loci that can resolve relationships among taxonomically challenging lineages. Diatoms are a diverse group of aquatic microalgae that includes important bioindicator species and many such lineages. One example is Nitzschia palea, a widespread species complex with several morphologically defined taxonomic varieties, some of which are critical pollution indicators. Morphological differences among the varieties are subtle and phylogenetic studies based on a few genes fail to resolve their evolutionary relationships. We conducted morphometric and transcriptome analyses of 10 Nitzschia palea strains to resolve the relationships among strains and taxonomic varieties. Nitzschia palea was resolved into three clades, one of which corresponds to a group of strains with narrow linear-lanceolate valves. The other morphological group recovered in the shape outline analysis was not monophyletic and consisted of two clades. Gene-tree concordance analyses and phylogenetic network estimations revealed patterns of incomplete lineage sorting and gene flow between intraspecific lineages. We detected reticulated evolutionary patterns among lineages with different morphologies, resulting in a putative recent hybrid. Our study shows that phylogenomic analyses of unlinked nuclear loci, complemented with morphometrics, can resolve complex evolutionary histories of recently diverged species complexes., (© 2022 The Authors. Journal of Phycology published by Wiley Periodicals LLC on behalf of Phycological Society of America.)

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 acclimated S. marinoi grown in a common garden. Our experiment included eight strains from source populations spanning the Baltic Sea salinity cline. Gene expression analysis revealed that low salinities induced changes in the cellular metabolism of S. marinoi, including upregulation of photosynthesis and storage compound biosynthesis, increased nutrient demand, and a complex response to oxidative stress. However, the strain effect overshadowed the salinity effect, as strains differed significantly in their response, both regarding the strength and the strategy (direction of gene expression) of their response. 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 to environmental change., (© 2022. The Author(s).)

Intricate networks of single-celled eukaryotes (protists) dominate carbon flow in the ocean. Their growth, demise, and interactions with other microorganisms drive the fluxes of biogeochemical elements through marine ecosystems. Mixotrophic protists are capable of both photosynthesis and ingestion of prey and are dominant components of open-ocean planktonic communities. Yet the role of mixotrophs in elemental cycling is obscured by their capacity to act as primary producers or heterotrophic consumers depending on factors that remain largely uncharacterized. Here, we develop and apply a machine learning model that predicts the in situ trophic mode of aquatic protists based on their patterns of gene expression. This approach leverages a public collection of protist transcriptomes as a training set to identify a subset of gene families whose transcriptional profiles predict trophic mode. We applied our model to nearly 100 metatranscriptomes obtained during two oceanographic cruises in the North Pacific and found community-level and population-specific evidence that abundant open-ocean mixotrophic populations shift their predominant mode of nutrient and carbon acquisition in response to natural gradients in nutrient supply and sea surface temperature. Metatranscriptomic data from ship-board incubation experiments revealed that abundant mixotrophic prymnesiophytes from the oligotrophic North Pacific subtropical gyre rapidly remodeled their transcriptome to enhance photosynthesis when supplied with limiting nutrients. Coupling this approach with experiments designed to reveal the mechanisms driving mixotroph physiology provides an avenue toward understanding the ecology of mixotrophy in the natural environment., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

Although most of the tens of thousands of diatom species are photoautotrophs, a small number of heterotrophic species no longer photosynthesize. We sequenced the genome of a nonphotosynthetic diatom, Nitzschia Nitz4, to determine how carbon metabolism was altered in the wake of this trophic shift. Nitzschia Nitz4 has retained its plastid and plastid genome, but changes associated with the transition to heterotrophy were cellular-wide and included losses of photosynthesis-related genes from the nuclear and plastid genomes, elimination of isoprenoid biosynthesis in the plastid, and remodeling of mitochondrial glycolysis to maximize adenosine triphosphte (ATP) yield. The genome contains a β-ketoadipate pathway that may allow Nitzschia Nitz4 to metabolize lignin-derived compounds. Diatom plastids lack an oxidative pentose phosphate pathway (oPPP), leaving photosynthesis as the primary source of NADPH to support essential biosynthetic pathways in the plastid and, by extension, limiting available sources of NADPH in nonphotosynthetic plastids. The genome revealed similarities between nonphotosynthetic diatoms and apicomplexan parasites for provisioning NADPH in their plastids and highlighted the ancestral absence of a plastid oPPP as a potentially important constraint on loss of photosynthesis, a hypothesis supported by the higher frequency of transitions to parasitism or heterotrophy in lineages that have a plastid oPPP., (© 2021 The Authors. New Phytologist © 2021 New Phytologist Foundation.)

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., (© 2020 Phycological Society of America.)

We report draft genomes of two bacterial strains in the genera Hyphobacterium and Reichenbachiella , which are associated with the diatom Cyclotella cryptica strain CCMP332. Genomes from these strains were 2,691,501 bp and 3,325,829 bp, respectively, and will be useful for understanding interactions between diatoms and bacteria., (Copyright © 2020 Gargas et al.)

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., (Copyright © 2020 Roberts et al.)

Many clades that span the marine-freshwater boundary are disproportionately more diverse in the younger, shorter lived, and scarcer freshwater environments than they are in the marine realm. This disparity is thought to be related to differences in diversification rates between marine and freshwater lineages. However, marine and freshwaters are not ecologically homogeneous, so the study of diversification across the salinity divide should also account for other potentially interacting variables. In diatoms, freshwater and substrate-associated (benthic) lineages are several-fold more diverse than their marine and suspended (planktonic) counterparts. These imbalances provide an excellent system to understand whether these variables interact with diversification. Using multistate hidden-state speciation and extinction models, we found that freshwater lineages diversify faster than marine lineages regardless of whether they inhabit the plankton or the benthos. Freshwater lineages also had higher turnover rates (speciation + extinction), suggesting that habitat transitions impact speciation and extinction rates jointly. The plankton-benthos contrast was also consistent with state-dependent diversification, but with modest differences in diversification and turnover rates. Asymmetric and bidirectional transitions rejected hypotheses about the plankton and freshwaters as absorbing, inescapable habitats. Our results further suggest that the high turnover rate of freshwater diatoms is related to high turnover of freshwater systems themselves., (© 2019 The Author(s). Evolution © 2019 The Society for the Study of Evolution.)

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., (© 2019 Botanical Society of America.)

Metabarcoding has offered unprecedented insights into microbial diversity. In many studies, short DNA sequences are binned into consecutively lower Linnaean ranks, and ranked groups (e.g., genera) are the units of biodiversity analyses. These analyses assume that Linnaean ranks are biologically meaningful and that identically ranked groups are comparable. We used a metabarcode dataset for marine planktonic diatoms to illustrate the limits of this approach. We found that the 20 most abundant marine planktonic diatom genera ranged in age from 4 to 134 million years, indicating the non-equivalence of genera because some have had more time to diversify than others. However, species richness was largely independent of genus age, suggesting that disparities in species richness among genera were better explained by variation in rates of speciation and extinction. Taxonomic classifications often do not reflect phylogeny, so genus-level analyses can include phylogenetically nested genera, further confounding rank-based analyses. These results underscore the indispensable role of phylogeny in understanding patterns of microbial diversity.

Patterns of species richness are commonly linked to life history strategies. In diatoms, an exceptionally diverse lineage of photosynthetic heterokonts important for global photosynthesis and burial of atmospheric carbon, lineages with different locomotory and reproductive traits differ dramatically in species richness, but any potential association between life history strategy and diversification has not been tested in a phylogenetic framework. We constructed a time-calibrated, 11-gene, 1151-taxon phylogeny of diatoms - the most inclusive diatom species tree to date. We used this phylogeny, together with a comprehensive inventory of first-last occurrences of Cenozoic fossil diatoms, to estimate ranges of expected species richness, diversification and its variation through time and across lineages. Diversification rates varied with life history traits. Although anisogamous lineages diversified faster than oogamous ones, this increase was restricted to a nested clade with active motility in the vegetative cells. We propose that the evolution of motility in vegetative cells, following an earlier transition from oogamy to anisogamy, facilitated outcrossing and improved utilization of habitat complexity, ultimately leading to enhanced opportunity for adaptive divergence across a variety of novel habitats. Together, these contributed to a species radiation that gave rise to the majority of present-day diatom diversity., (© 2018 The Authors. New Phytologist © 2018 New Phytologist Trust.)

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.

Premise of the Study: Diatoms are one of the most species-rich lineages of microbial eukaryotes. Similarities in clade age, species richness, and primary productivity motivate comparisons to angiosperms, whose genomes have been inordinately shaped by whole-genome duplication (WGD). WGDs have been linked to speciation, increased rates of lineage diversification, and identified as a principal driver of angiosperm evolution. We synthesized a large but scattered body of evidence that suggests polyploidy may be common in diatoms as well., Methods: We used gene counts, gene trees, and distributions of synonymous divergence to carry out a phylogenomic analysis of WGD across a diverse set of 37 diatom species., Key Results: Several methods identified WGDs of varying age across diatoms. Determining the occurrence, exact number, and placement of events was greatly impacted by uncertainty in gene trees. WGDs inferred from synonymous divergence of paralogs varied depending on how redundancy in transcriptomes was assessed, gene families were assembled, and synonymous distances (Ks) were calculated. Our results highlighted a need for systematic evaluation of key methodological aspects of Ks-based approaches to WGD inference. Gene tree reconciliations supported allopolyploidy as the predominant mode of polyploid formation, with strong evidence for ancient allopolyploid events in the thalassiosiroid and pennate diatom clades., Conclusions: Our results suggest that WGD has played a major role in the evolution of diatom genomes. We outline challenges in reconstructing paleopolyploid events in diatoms that, together with these results, offer a framework for understanding the impact of genome duplication in a group that likely harbors substantial genomic diversity., (© 2018 The Authors. American Journal of Botany is published by Wiley Periodicals, Inc. on behalf of the Botanical Society of America.)

Diatoms (Bacillariophyta) are a species-rich group of eukaryotic microbes diverse in morphology, ecology, and metabolism. Previous reconstructions of the diatom phylogeny based on one or a few genes have resulted in inconsistent resolution or low support for critical nodes. We applied phylogenetic paralog pruning techniques to a data set of 94 diatom genomes and transcriptomes to infer perennially difficult species relationships, using concatenation and summary-coalescent methods to reconstruct species trees from data sets spanning a wide range of thresholds for taxon and column occupancy in gene alignments. Conflicts between gene and species trees decreased with both increasing taxon occupancy and bootstrap cutoffs applied to gene trees. Concordance between gene and species trees was lowest for short internodes and increased logarithmically with increasing edge length, suggesting that incomplete lineage sorting disproportionately affects species tree inference at short internodes, which are a common feature of the diatom phylogeny. Although species tree topologies were largely consistent across many data treatments, concatenation methods appeared to outperform summary-coalescent methods for sparse alignments. Our results underscore that approaches to species-tree inference based on few loci are likely to be misled by unrepresentative sampling of gene histories, particularly in lineages that may have diversified rapidly. In addition, phylogenomic studies of diatoms, and potentially other hyperdiverse groups, should maximize the number of gene trees with high taxon occupancy, though there is clearly a limit to how many of these genes will be available., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., (Copyright © 2016 Elsevier Inc. All rights reserved.)

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., (© 2014 Phycological Society of America.)

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.

We report the complete mitochondrial genome sequence of the flowering plant Amborella trichopoda. This enormous, 3.9-megabase genome contains six genome equivalents of foreign mitochondrial DNA, acquired from green algae, mosses, and other angiosperms. Many of these horizontal transfers were large, including acquisition of entire mitochondrial genomes from three green algae and one moss. We propose a fusion-compatibility model to explain these findings, with Amborella capturing whole mitochondria from diverse eukaryotes, followed by mitochondrial fusion (limited mechanistically to green plant mitochondria) and then genome recombination. Amborella's epiphyte load, propensity to produce suckers from wounds, and low rate of mitochondrial DNA loss probably all contribute to the high level of foreign DNA in its mitochondrial genome.

The hybrid pigeonpea (Cajanus cajan) breeding technology based on cytoplasmic male sterility (CMS) is currently unique among legumes and displays major potential for yield increase. CMS is defined as a condition in which a plant is unable to produce functional pollen grains. The novel chimeric open reading frames (ORFs) produced as a results of mitochondrial genome rearrangements are considered to be the main cause of CMS. To identify these CMS-related ORFs in pigeonpea, we sequenced the mitochondrial genomes of three C. cajan lines (the male-sterile line ICPA 2039, the maintainer line ICPB 2039, and the hybrid line ICPH 2433) and of the wild relative (Cajanus cajanifolius ICPW 29). A single, circular-mapping molecule of length 545.7 kb was assembled and annotated for the ICPA 2039 line. Sequence annotation predicted 51 genes, including 34 protein-coding and 17 RNA genes. Comparison of the mitochondrial genomes from different Cajanus genotypes identified 31 ORFs, which differ between lines within which CMS is present or absent. Among these chimeric ORFs, 13 were identified by comparison of the related male-sterile and maintainer lines. These ORFs display features that are known to trigger CMS in other plant species and to represent the most promising candidates for CMS-related mitochondrial rearrangements in pigeonpea.

Light-harvesting antennae are critical for collecting energy from sunlight and providing it to photosynthetic reaction centers. Their abundance and composition are tightly regulated to maintain efficient photosynthesis in changing light conditions. Many cyanobacteria alter their light-harvesting antennae in response to changes in ambient light-color conditions through the process of chromatic acclimation. The control of green light induction (Cgi) pathway is a light-color-sensing system that controls the expression of photosynthetic genes during chromatic acclimation, and while some evidence suggests that it operates via transcription attenuation, the components of this pathway have not been identified. We provide evidence that translation initiation factor 3 (IF3), an essential component of the prokaryotic translation initiation machinery that binds the 30S subunit and blocks premature association with the 50S subunit, is part of the control of green light induction pathway. Light regulation of gene expression has not been previously described for any translation initiation factor. Surprisingly, deletion of the IF3-encoding gene infCa was not lethal in the filamentous cyanobacterium Fremyella diplosiphon, and its genome was found to contain a second, redundant, highly divergent infC gene which, when deleted, had no effect on photosynthetic gene expression. Either gene could complement an Escherichia coli infC mutant and thus both encode bona fide IF3s. Analysis of prokaryotic and eukaryotic genome databases established that multiple infC genes are present in the genomes of diverse groups of bacteria and land plants, most of which do not undergo chromatic acclimation. This suggests that IF3 may have repeatedly evolved important roles in the regulation of gene expression in both prokaryotes and eukaryotes.

The phenotype represents a critical interface between the genome and the environment in which organisms live and evolve. Phenotypic characters also are a rich source of biodiversity data for tree building, and they enable scientists to reconstruct the evolutionary history of organisms, including most fossil taxa, for which genetic data are unavailable. Therefore, phenotypic data are necessary for building a comprehensive Tree of Life. In contrast to recent advances in molecular sequencing, which has become faster and cheaper through recent technological advances, phenotypic data collection remains often prohibitively slow and expensive. The next-generation phenomics project is a collaborative, multidisciplinary effort to leverage advances in image analysis, crowdsourcing, and natural language processing to develop and implement novel approaches for discovering and scoring the phenome, the collection of phentotypic characters for a species. This research represents a new approach to data collection that has the potential to transform phylogenetics research and to enable rapid advances in constructing the Tree of Life. Our goal is to assemble large phenomic datasets built using new methods and to provide the public and scientific community with tools for phenomic data assembly that will enable rapid and automated study of phenotypes across the Tree of Life.

Background: The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms., Results: We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced., Conclusions: Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

Genome size and complexity vary tremendously among eukaryotic species and their organelles. Comparisons across deeply divergent eukaryotic lineages have suggested that variation in mutation rates may explain this diversity, with increased mutational burdens favoring reduced genome size and complexity. The discovery that mitochondrial mutation rates can differ by orders of magnitude among closely related angiosperm species presents a unique opportunity to test this hypothesis. We sequenced the mitochondrial genomes from two species in the angiosperm genus Silene with recent and dramatic accelerations in their mitochondrial mutation rates. Contrary to theoretical predictions, these genomes have experienced a massive proliferation of noncoding content. At 6.7 and 11.3 Mb, they are by far the largest known mitochondrial genomes, larger than most bacterial genomes and even some nuclear genomes. In contrast, two slowly evolving Silene mitochondrial genomes are smaller than average for angiosperms. Consequently, this genus captures approximately 98% of known variation in organelle genome size. The expanded genomes reveal several architectural changes, including the evolution of complex multichromosomal structures (with 59 and 128 circular-mapping chromosomes, ranging in size from 44 to 192 kb). They also exhibit a substantial reduction in recombination and gene conversion activity as measured by the relative frequency of alternative genome conformations and the level of sequence divergence between repeat copies. The evolution of mutation rate, genome size, and chromosome structure can therefore be extremely rapid and interrelated in ways not predicted by current evolutionary theories. Our results raise the hypothesis that changes in recombinational processes, including gene conversion, may be a central force driving the evolution of both mutation rate and genome structure., Competing Interests: The authors have declared that no competing interests exist.

The angiosperm genus Silene exhibits some of the most extreme and rapid divergence ever identified in mitochondrial genome architecture and nucleotide substitution rates. These patterns have been considered mitochondrial specific based on the absence of correlated changes in the small number of available nuclear and plastid gene sequences. To better assess the relationship between mitochondrial and plastid evolution, we sequenced the plastid genomes from four Silene species with fully sequenced mitochondrial genomes. We found that two species with fast-evolving mitochondrial genomes, S. noctiflora and S. conica, also exhibit accelerated rates of sequence and structural evolution in their plastid genomes. The nature of these changes, however, is markedly different from those in the mitochondrial genome. For example, in contrast to the mitochondrial pattern, which appears to be genome wide and mutationally driven, the plastid substitution rate accelerations are restricted to a subset of genes and preferentially affect nonsynonymous sites, indicating that altered selection pressures are acting on specific plastid-encoded functions in these species. Indeed, some plastid genes in S. noctiflora and S. conica show strong evidence of positive selection. In contrast, two species with more slowly evolving mitochondrial genomes, S. latifolia and S. vulgaris, have correspondingly low rates of nucleotide substitution in plastid genes as well as a plastid genome structure that has remained essentially unchanged since the origin of angiosperms. These results raise the possibility that common evolutionary forces could be shaping the extreme but distinct patterns of divergence in both organelle genomes within this genus.

Members of the flowering plant family Cucurbitaceae harbor the largest known mitochondrial genomes. Here, we report the 1685-kb mitochondrial genome of cucumber (Cucumis sativus). We help solve a 30-year mystery about the origins of its large size by showing that it mainly reflects the proliferation of dispersed repeats, expansions of existing introns, and the acquisition of sequences from diverse sources, including the cucumber nuclear and chloroplast genomes, viruses, and bacteria. The cucumber genome has a novel structure for plant mitochondria, mapping as three entirely or largely autonomous circular chromosomes (lengths 1556, 84, and 45 kb) that vary in relative abundance over a twofold range. These properties suggest that the three chromosomes replicate independently of one another. The two smaller chromosomes are devoid of known functional genes but nonetheless contain diagnostic mitochondrial features. Paired-end sequencing conflicts reveal differences in recombination dynamics among chromosomes, for which an explanatory model is developed, as well as a large pool of low-frequency genome conformations, many of which may result from asymmetric recombination across intermediate-sized and sometimes highly divergent repeats. These findings highlight the promise of genome sequencing for elucidating the recombinational dynamics of plant mitochondrial genomes.

Background: Publication of the first diatom genome, that of Thalassiosira pseudonana, established it as a model species for experimental and genomic studies of diatoms. Virtually every ensuing study has treated T. pseudonana as a marine diatom, with genomic and experimental data valued for their insights into the ecology and evolution of diatoms in the world's oceans., Results: The natural distribution of T. pseudonana spans both marine and fresh waters, and phylogenetic analyses of morphological and molecular datasets show that, 1) T. pseudonana marks an early divergence in a major freshwater radiation by diatoms, and 2) as a species, T. pseudonana is likely ancestrally freshwater. Marine strains therefore represent recent recolonizations of higher salinity habitats. In addition, the combination of a relatively nondescript form and a convoluted taxonomic history has introduced some confusion about the identity of T. pseudonana and, by extension, its phylogeny and ecology. We resolve these issues and use phylogenetic criteria to show that T. pseudonana is more appropriately classified by its original name, Cyclotella nana. Cyclotella contains a mix of marine and freshwater species and so more accurately conveys the complexities of the phylogenetic and natural histories of T. pseudonana., Conclusions: The multitude of physical barriers that likely must be overcome for diatoms to successfully colonize freshwaters suggests that the physiological traits of T. pseudonana, and the genes underlying those traits, might differ from those of strictly marine diatoms. The freshwater ancestry of T. pseudonana might therefore confound generalizations about the physiological and metabolic properties of marine diatoms. The freshwater component of T. pseudonana's history merits careful consideration in the interpretation of experimental data collected for this important model species.

The mitochondrial genomes of seed plants are exceptionally fluid in size, structure, and sequence content, with the accumulation and activity of repetitive sequences underlying much of this variation. We report the first fully sequenced mitochondrial genome of a legume, Vigna radiata (mung bean), and show that despite its unexceptional size (401,262 nt), the genome is unusually depauperate in repetitive DNA and "promiscuous" sequences from the chloroplast and nuclear genomes. Although Vigna lacks the large, recombinationally active repeats typical of most other seed plants, a PCR survey of its modest repertoire of short (38-297 nt) repeats nevertheless revealed evidence for recombination across all of them. A set of novel control assays showed, however, that these results could instead reflect, in part or entirely, artifacts of PCR-mediated recombination. Consequently, we recommend that other methods, especially high-depth genome sequencing, be used instead of PCR to infer patterns of plant mitochondrial recombination. The average-sized but repeat- and feature-poor mitochondrial genome of Vigna makes it ever more difficult to generalize about the factors shaping the size and sequence content of plant mitochondrial genomes.

Background: Mitochondrial gene loss and functional transfer to the nucleus is an ongoing process in many lineages of plants, resulting in substantial variation across species in mitochondrial gene content. The Caryophyllaceae represents one lineage that has experienced a particularly high rate of mitochondrial gene loss relative to other angiosperms., Results: In this study, we report the first complete mitochondrial genome sequence from a member of this family, Silene latifolia. The genome can be mapped as a 253,413 bp circle, but its structure is complicated by a large repeated region that is present in 6 copies. Active recombination among these copies produces a suite of alternative genome configurations that appear to be at or near "recombinational equilibrium". The genome contains the fewest genes of any angiosperm mitochondrial genome sequenced to date, with intact copies of only 25 of the 41 protein genes inferred to be present in the common ancestor of angiosperms. As observed more broadly in angiosperms, ribosomal proteins have been especially prone to gene loss in the S. latifolia lineage. The genome has also experienced a major reduction in tRNA gene content, including loss of functional tRNAs of both native and chloroplast origin. Even assuming expanded wobble-pairing rules, the mitochondrial genome can support translation of only 17 of the 61 sense codons, which code for only 9 of the 20 amino acids. In addition, genes encoding 18S and, especially, 5S rRNA exhibit exceptional sequence divergence relative to other plants. Divergence in one region of 18S rRNA appears to be the result of a gene conversion event, in which recombination with a homologous gene of chloroplast origin led to the complete replacement of a helix in this ribosomal RNA., Conclusions: These findings suggest a markedly expanded role for nuclear gene products in the translation of mitochondrial genes in S. latifolia and raise the possibility of altered selective constraints operating on the mitochondrial translational apparatus in this lineage.

Theoretical arguments suggest that mutation rates influence the proliferation and maintenance of RNA editing. We identified RNA editing sites in five species within the angiosperm genus Silene that exhibit highly divergent mitochondrial mutation rates. We found that mutational acceleration has been associated with rapid loss of mitochondrial editing sites. In contrast, we did not find a significant difference in the frequency of editing in chloroplast genes, which lack the mutation rate variation observed in the mitochondrial genome. As found in other angiosperms, the rate of substitution at RNA editing sites in Silene greatly exceeds the rate at synonymous sites, a pattern that has previously been interpreted as evidence for selection against RNA editing. Alternatively, we suggest that editing sites may experience higher rates of C-to-T mutation than other portions of the genome. Such a pattern could be caused by gene conversion with reverse-transcribed mRNA (i.e., retroprocessing). If so, the genomic distribution of RNA editing site losses in Silene suggests that such conversions must be occurring at a local scale such that only one or two editing sites are affected at a time. Because preferential substitution at editing sites appears to occur in angiosperms regardless of the mutation rate, we conclude that mitochondrial rate accelerations within Silene have "fast-forwarded" a preexisting pattern but have not fundamentally changed the evolutionary forces acting on RNA editing sites.

The mitochondrial genomes of seed plants are unusually large and vary in size by at least an order of magnitude. Much of this variation occurs within a single family, the Cucurbitaceae, whose genomes range from an estimated 390 to 2,900 kb in size. We sequenced the mitochondrial genomes of Citrullus lanatus (watermelon: 379,236 nt) and Cucurbita pepo (zucchini: 982,833 nt)--the two smallest characterized cucurbit mitochondrial genomes--and determined their RNA editing content. The relatively compact Citrullus mitochondrial genome actually contains more and longer genes and introns, longer segmental duplications, and more discernibly nuclear-derived DNA. The large size of the Cucurbita mitochondrial genome reflects the accumulation of unprecedented amounts of both chloroplast sequences (>113 kb) and short repeated sequences (>370 kb). A low mutation rate has been hypothesized to underlie increases in both genome size and RNA editing frequency in plant mitochondria. However, despite its much larger genome, Cucurbita has a significantly higher synonymous substitution rate (and presumably mutation rate) than Citrullus but comparable levels of RNA editing. The evolution of mutation rate, genome size, and RNA editing are apparently decoupled in Cucurbitaceae, reflecting either simple stochastic variation or governance by different factors.

A recent reclassification of diatoms based on phylogenies recovered using the nuclear-encoded SSU rRNA gene contains three major classes, Coscinodiscophyceae, Mediophyceae and the Bacillariophyceae (the CMB hypothesis). We evaluated this with a sequence alignment of 1336 protist and heterokont algae SSU rRNAs, which includes 673 diatoms. Sequences were aligned to maintain structural elements conserved within this dataset. Parsimony analysis rejected the CMB hypothesis, albeit weakly. Morphological data are also incongruent with this recent CMB hypothesis of three diatom clades. We also reanalyzed a recently published dataset which purports to support the CMB hypothesis. Our reanalysis found that the original analysis had not converged on the true bipartition posterior probability distribution, and rejected the CMB hypothesis. Thus we conclude that a reclassification of the evolutionary relationships of the diatoms according to the CMB hypothesis is premature.

Horizontal gene transfer is surprisingly common among plant mitochondrial genomes. The first well-established case involves a homing group I intron in the mitochondrial cox1 gene shown to have been frequently acquired via horizontal transfer in angiosperms. Here, we report extensive additional sampling of angiosperms, including 85 newly sequenced introns from 30 families. Analysis of all available data leads us to conclude that, among the 640 angiosperms (from 212 families) whose cox1 intron status has been characterized thus far, the intron has been acquired via roughly 70 separate horizontal transfer events. We propose that the intron was originally seeded into angiosperms by a single transfer from fungi, with all subsequent inferred transfers occurring from one angiosperm to another. The pattern of angiosperm-to-angiosperm transfer is biased toward exchanges between plants belonging to the same family. Illegitimate pollination is proposed as one potential factor responsible for this pattern, given that aberrant, cross-species pollination is more likely between close relatives. Other potential factors include shared vectoring agents or common geographic locations. We report the first apparent cases of loss of the cox1 intron; losses are accompanied by retention of the exonic coconversion tract, which is located immediately downstream of the intron and which is a product of the intron's self-insertion mechanism. We discuss the many reasons why the cox1 intron is so frequently and detectably transferred, and rarely lost, and conclude that it should be regarded as the "canary in the coal mine" with respect to horizontal transfer in angiosperm mitochondria.

Salinity imposes a significant barrier to the distribution of many organisms, including diatoms. Diatoms are ancestrally marine, and the number of times they have independently colonized fresh waters and the physiological adaptations that facilitated these transitions remain outstanding questions in diatom evolution. The colonization of fresh waters by diatoms has been compared to "crossing the Rubicon," implying that successful colonization events are rare, irreversible, and lead to substantial species diversification. To test these hypotheses, we reconstructed the phylogeny of Thalassiosirales, a diatom lineage with high diversity in both marine and fresh waters. We collected approximately 5.3kb of DNA sequence data from the nuclear (SSU and partial LSU rDNA) and chloroplast genomes (psbC and rbcL) and reconstructed the phylogeny using parsimony and Bayesian methods. Alternative topology tests strongly reject all previous colonization hypotheses, including monophyly of the predominantly freshwater Stephanodiscaceae. Results showed at least three independent colonizations of fresh waters, and whereas previous accounts of freshwater-to-marine transitions have been discounted, these results provide compelling evidence for as many as three independent re-colonizations of the marine habitat, two of which led to speciation events. This study adds valuable phylogenetic context to previous debate about the nature of the salinity barrier in diatoms and provides compelling evidence that, at least for Thalassiosirales, the salinity barrier might be less formidable than previously thought.

Background: The Vitaceae (grape) is an economically important family of angiosperms whose phylogenetic placement is currently unresolved. Recent phylogenetic analyses based on one to several genes have suggested several alternative placements of this family, including sister to Caryophyllales, asterids, Saxifragales, Dilleniaceae or to rest of rosids, though support for these different results has been weak. There has been a recent interest in using complete chloroplast genome sequences for resolving phylogenetic relationships among angiosperms. These studies have clarified relationships among several major lineages but they have also emphasized the importance of taxon sampling and the effects of different phylogenetic methods for obtaining accurate phylogenies. We sequenced the complete chloroplast genome of Vitis vinifera and used these data to assess relationships among 27 angiosperms, including nine taxa of rosids., Results: The Vitis vinifera chloroplast genome is 160,928 bp in length, including a pair of inverted repeats of 26,358 bp that are separated by small and large single copy regions of 19,065 bp and 89,147 bp, respectively. The gene content and order of Vitis is identical to many other unrearranged angiosperm chloroplast genomes, including tobacco. Phylogenetic analyses using maximum parsimony and maximum likelihood were performed on DNA sequences of 61 protein-coding genes for two datasets with 28 or 29 taxa, including eight or nine taxa from four of the seven currently recognized major clades of rosids. Parsimony and likelihood phylogenies of both data sets provide strong support for the placement of Vitaceae as sister to the remaining rosids. However, the position of the Myrtales and support for the monophyly of the eurosid I clade differs between the two data sets and the two methods of analysis. In parsimony analyses, the inclusion of Gossypium is necessary to obtain trees that support the monophyly of the eurosid I clade. However, maximum likelihood analyses place Cucumis as sister to the Myrtales and therefore do not support the monophyly of the eurosid I clade., Conclusion: Phylogenies based on DNA sequences from complete chloroplast genome sequences provide strong support for the position of the Vitaceae as the earliest diverging lineage of rosids. Our phylogenetic analyses support recent assertions that inadequate taxon sampling and incorrect model specification for concatenated multi-gene data sets can mislead phylogenetic inferences when using whole chloroplast genomes for phylogeny reconstruction.

During the past decade, there has been a rapid increase in our understanding of plastid genome organization and evolution due to the availability of many new completely sequenced genomes. There are 45 complete genomes published and ongoing projects are likely to increase this sampling to nearly 200 genomes during the next 5 years. Several groups of researchers including ours have been developing new techniques for gathering and analyzing entire plastid genome sequences and details of these developments are summarized in this chapter. The most important developments that enhance our ability to generate whole chloroplast genome sequences involve the generation of pure fractions of chloroplast genomes by whole genome amplification using rolling circle amplification, cloning genomes into Fosmid or bacterial artificial chromosome (BAC) vectors, and the development of an organellar annotation program (Dual Organellar GenoMe Annotator [DOGMA]). In addition to providing details of these methods, we provide an overview of methods for analyzing complete plastid genome sequences for repeats and gene content, as well as approaches for using gene order and sequence data for phylogeny reconstruction. This explosive increase in the number of sequenced plastid genomes and improved computational tools will provide many insights into the evolution of these genomes and much new data for assessing relationships at deep nodes in plants and other photosynthetic organisms.

Much attention is being directed toward understanding and co-opting the mechanisms by which diatoms control patterning of silica for nanotechnology applications. Given the enormous diversity of cell-wall patterns among diatom taxa, a complete understanding of these mechanisms will require comparative analysis of nanopatterning strategies from a diversity of diatom taxa. A well-supported phylogenetic tree provides the best basis for such comparisons because closely related taxa will have many similar attributes (morphological, physiological, ecological, etc.) simply because they share a recent common ancestor. Phylogenetic information is exploited in as diverse a set of fields as pharmacology, epidemiology, and geology. Nanotechnologists can use the phylogenetic tree of diatoms to help select exemplar taxa and to streamline the search for alternative nanopatterning strategies. We review the progress made thus far in reconstructing the phylogeny of diatoms. Most analyses have been based on small subunit ribosomal DNA sequences, and inferences have varied substantially depending on the number and diversity of taxa included in the analyses. We review several seminal studies on diatom phylogeny in light of theoretical and empirical investigations that have emphasized the critical importance of taxonomic sampling on phylogenetic inference. One consistent result is that centric diatoms grade into araphid pennates, and araphid pennates grade the raphe-bearing pennate diatoms, which are a natural group. Including more taxa and more characters (molecular and morphological) should yield a better supported hypothesis of diatom relationships.