Lower Statistical Support with Larger Data Sets: Insights from the Ochrophyta Radiation.

It is commonly assumed that increasing the number of characters has the potential to resolve evolutionary radiations. Here, we studied photosynthetic stramenopiles (Ochrophyta) using alignments of heterogeneous origin mitochondrion, plastid, and nucleus. Surprisingly while statistical support for the relationships between the six major Ochrophyta lineages increases when comparing the mitochondrion (6,762 sites) and plastid (21,692 sites) trees, it decreases in the nuclear (209,105 sites) tree. Statistical support is not simply related to the data set size but also to the quantity of phylogenetic signal available at each position and our ability to extract it. Here, we show that this ability for current phylogenetic methods is limited, because conflicting results were obtained when varying taxon sampling. Even though the use of a better fitting model improved signal extraction and reduced the observed conflicts, the plastid data set provided higher statistical support for the ochrophyte radiation than the larger nucleus data set. We propose that the higher support observed in the plastid tree is due to an acceleration of the evolutionary rate in one short deep internal branch, implying that more phylogenetic signal per position is available to resolve the Ochrophyta radiation in the plastid than in the nuclear data set. Our work therefore suggests that, in order to resolve radiations, beyond the obvious use of data sets with more positions, we need to continue developing models of sequence evolution that better extract the phylogenetic signal and design methods to search for genes/characters that contain more signal specifically for short internal branches.

Presence of male mitochondria in somatic tissues and their functional importance at the whole animal level in the marine bivalve Arctica islandica.

Metazoans normally possess a single lineage of mitochondria inherited from the mother (♀-type mitochondria) while paternal mitochondria are absent or eliminated in fertilized eggs. In doubly uniparental inheritance (DUI), which is specific to the bivalve clade including the ocean quahog, Arctica islandica, ♂-type mitochondria are retained in male gonads and, in a few species, small proportions of ♂-type mitochondria co-exist with ♀-type in somatic tissues. To the best of our knowledge, we report, for the first time in metazoan, the natural occurrence of male and female individuals with exclusively ♂-type mitochondria in somatic tissues of the bivalve A. islandica. Mitochondrial genomes differ by ~5.5% at DNA sequence level. Exclusive presence of ♂-type mitochondria affects mitochondrial complexes partially encoded by mitochondrial genes and leads to a sharp drop in respiratory capacity. Through a combination of whole mitochondrial genome sequencing and molecular assays (gene presence and expression), we demonstrate that 1) 11% of individuals of an Icelandic population appear homoplasmic for ♂-type mitochondria in somatic tissues, 2) ♂-type mitochondrial genes are transcribed and 3) individuals with ♂-type mitochondria in somatic cells lose 30% of their wild-type respiratory capacity. This mitochondrial pattern in A. islandica is a special case of DUI, highlighted in individuals from both sexes with functional consequences at cellular and conceivably whole animal level.

RNAseq and quantitative proteomic analysis of Dictyostelium knock-out cells lacking the core autophagy proteins ATG9 and/or ATG16.

BACKGROUND: Autophagy is an evolutionary ancient mechanism that sequesters substrates for degradation within autolysosomes. The process is driven by many autophagy-related (ATG) proteins, including the core members ATG9 and ATG16. However, the functions of these two core ATG proteins still need further elucidation. Here, we applied RNAseq and tandem mass tag (TMT) proteomic approaches to identify differentially expressed genes (DEGs) and proteins (DEPs) in Dictyostelium discoideum ATG9‾, ATG16‾ and ATG9‾/16‾ strains in comparison to AX2 wild-type cells. RESULT: In total, we identified 332 (279 up and 53 down), 639 (487 up and 152 down) and 260 (114 up and 146 down) DEGs and 124 (83 up and 41 down), 431 (238 up and 193 down) and 677 (347 up and 330 down) DEPs in ATG9‾, ATG16‾ and ATG9‾/16‾ strains, respectively. Thus, in the single knock-out strains, the number of DEGs was higher than the number of DEPs while in the double knock-out strain the number of DEPs was higher. Comparison of RNAseq and proteomic data further revealed, that only a small proportion of the transcriptional changes were reflected on the protein level. Gene ontology (GO) analysis revealed an enrichment of DEPs involved in lipid metabolism and oxidative phosphorylation. Furthermore, we found increased expression of the anti-oxidant enzymes glutathione reductase (gsr) and catalase A (catA) in ATG16‾ and ATG9‾/16‾ cells, respectively, indicating adaptation to excess reactive oxygen species (ROS). CONCLUSIONS: Our study provides the first combined transcriptome and proteome analysis of ATG9‾, ATG16‾ and ATG9‾/16‾ cells. Our results suggest, that most changes in protein abundance were not caused by transcriptional changes, but were rather due to changes in protein homeostasis. In particular, knock-out of atg9 and/or atg16 appears to cause dysregulation of lipid metabolism and oxidative phosphorylation.

Abundantly expressed class of noncoding RNAs conserved through the multicellular evolution of dictyostelid social amoebas.

Aggregative multicellularity has evolved multiple times in diverse groups of eukaryotes, exemplified by the well-studied development of dictyostelid social amoebas, for example, Dictyostelium discoideum However, it is still poorly understood why multicellularity emerged in these amoebas while the majority of other members of Amoebozoa are unicellular. Previously, a novel type of noncoding RNA, Class I RNAs, was identified in D. discoideum and shown to be important for normal multicellular development. Here, we investigated Class I RNA evolution and its connection to multicellular development. We identified a large number of new Class I RNA genes by constructing a covariance model combined with a scoring system based on conserved upstream sequences. Multiple genes were predicted in representatives of each major group of Dictyostelia and expression analysis confirmed that our search approach identifies expressed Class I RNA genes with high accuracy and sensitivity and that the RNAs are developmentally regulated. Further studies showed that Class I RNAs are ubiquitous in Dictyostelia and share highly conserved structure and sequence motifs. In addition, Class I RNA genes appear to be unique to dictyostelid social amoebas because they could not be identified in outgroup genomes, including their closest known relatives. Our results show that Class I RNA is an ancient class of ncRNAs, likely to have been present in the last common ancestor of Dictyostelia dating back at least 600 million years. Based on previous functional analyses and the presented evolutionary investigation, we hypothesize that Class I RNAs were involved in evolution of multicellularity in Dictyostelia.

Analysis of organellar genomes in brown algae reveals an independent introduction of similar foreign sequences into the mitochondrial genome.

Kelp species (Laminariales, Phaeophyceae) are globally widespread along temperate to Polar rocky coastal lines. Here we analyse the mitochondrial and chloroplast genomes of Laminaria rodriguezii, in comparison to the organellar genomes of other kelp species. We also provide the complete mitochondrial genome sequence of another endemic kelp species from a Polar habitat, the Arctic Laminaria solidungula. We compare phylogenetic trees derived from twenty complete mitochondrial and seven complete chloroplast kelp genomes. Interestingly, we found a stretch of more than 700 bp in the mitochondrial genome of L.rodriguezii, which is not present in any other yet sequenced member of the Phaeophyceae. This stretch matches a protein coding region in the mitochondrial genome from Desmarestia viridis, another brown seaweed. Their high similarity suggests that these sequences originated through independent introduction into the two species. Their origin could have been by infection by yet unknown similar mitoviruses, currently only known from fungi and plants.

Adaptive divergence across Southern Ocean gradients in the pelagic diatom Fragilariopsis kerguelensis.

The Southern Ocean is characterized by longitudinal water circulations crossed by strong latitudinal gradients. How this oceanographic background shapes planktonic populations is largely unknown, despite the significance of this region for global biogeochemical cycles. Here, we show, based on genomic, morphometric, ecophysiological and mating compatibility data, an example of ecotypic differentiation and speciation within an endemic pelagic inhabitant, the diatom Fragilariopsis kerguelensis. We discovered three genotypic variants, one present throughout the latitudinal transect sampled, the others restricted to the north and south, respectively. The latter two showed reciprocal monophyly across all three genomes and significant ecophysiological differences consistent with local adaptation, but produced viable offspring in laboratory crosses. The third group was also reproductively isolated from the latter two. We hypothesize that this pattern originated by an adaptive expansion accompanied by ecotypic divergence, followed by sympatric speciation.

Responses of the kelp Saccharina latissima (Phaeophyceae) to the warming Arctic: from physiology to transcriptomics.

The Arctic region is currently facing substantial environmental changes due to global warming. Melting glaciers cause reduced salinity environments in coastal Arctic habitats, which may be stressful for kelp beds. To investigate the responses of the kelp Saccharina latissima to the warming Arctic, we studied the transcriptomic changes of S. latissima from Kongsfjorden (Svalbard, Norway) over a 24-hour exposure to two salinities (Absolute Salinity [SA ] 20 and 30) after a 7-day pre-acclimation at three temperatures (0, 8 and 15°C). In addition, corresponding physiological data were assessed during an 11-days salinity/temperature experiment. Growth and maximal quantum yield for photosystem II fluorescence were positively affected by increased temperature during acclimation, whereas hyposalinity caused negative effects at the last day of treatment. In contrast, hyposalinity induced marked changes on the transcriptomic level. Compared to the control (8°C - SA 30), the 8°C - SA 20 exhibited the highest number of differentially expressed genes (DEGs), followed by the 0°C - SA 20. Comparisons indicate that S. latissima tends to convert its energy from primary metabolism (e.g. photosynthesis) to antioxidant activity under hyposaline stress. The increase in physiological performance at 15°C shows that S. latissima in the Arctic region can adjust and might even benefit from increased temperatures. However, in Arctic fjord environments its performance might become impaired by decreased salinity as a result of ice melting.

Is geographical variation driving the transcriptomic responses to multiple stressors in the kelp Saccharina latissima?

BACKGROUND: Kelps (Laminariales, Phaeophyceae) are brown macroalgae of utmost ecological, and increasingly economic, importance on temperate to polar rocky shores. Omics approaches in brown algae are still scarce and knowledge of their acclimation mechanisms to the changing conditions experienced in coastal environments can benefit from the application of RNA-sequencing. Despite evidence of ecotypic differentiation, transcriptomic responses from distinct geographical locations have, to our knowledge, never been studied in the sugar kelp Saccharina latissima so far. RESULTS: In this study we investigated gene expression responses using RNA-sequencing of S. latissima from environments with contrasting temperature and salinity conditions - Roscoff, in temperate eastern Atlantic, and Spitsbergen in the Arctic. Juvenile sporophytes derived from uniparental stock cultures from both locations were pre-cultivated at 8 °C and SA 30. Sporophytes acclimated to 0 °C, 8 °C and 15 °C were exposed to a low salinity treatment (SA 20) for 24 h. Hyposalinity had a greater impact at the transcriptomic level than the temperature alone, and its effects were modulated by temperature. Namely, photosynthesis and pigment synthesis were extensively repressed by low salinity at low temperatures. Although some responses were shared among sporophytes from the different sites, marked differences were revealed by principal component analysis, differential expression and GO enrichment. The interaction between low temperature and low salinity drove the largest changes in gene expression in sporophytes from Roscoff while specimens from Spitsbergen required more metabolic adjustment at higher temperatures. Moreover, genes related to cell wall adjustment were differentially expressed between Spitsbergen and Roscoff control samples. CONCLUSIONS: Our study reveals interactive effects of temperature and salinity on transcriptomic profiles in S. latissima. Moreover, our data suggest that under identical culture conditions sporophytes from different locations diverge in their transcriptomic responses. This is probably connected to variations in temperature and salinity in their respective environment of origin. The current transcriptomic results support the plastic response pattern in sugar kelp which is a species with several reported ecotypes. Our data provide the baseline for a better understanding of the underlying processes of physiological plasticity and may help in the future to identify strains adapted to specific environments and its genetic control.

Loss of a chloroplast encoded function could influence species range in kelp.

Kelps are important providers and constituents of marine ecological niches, the coastal kelp forests. Kelp species have differing distribution ranges, but mainly thrive in temperate and arctic regions. Although the principal factors determining biogeographic distribution ranges are known, genomics could provide additional answers to this question. We sequenced DNA from two Laminaria species with contrasting distribution ranges, Laminaria digitata and Laminaria solidungula. Laminaria digitata is found in the Northern Atlantic with a southern boundary in Brittany (France) or Massachusetts (USA) and a northern boundary in the Arctic, whereas L. solidungula is endemic to the Arctic only. From the raw reads of DNA, we reconstructed both chloroplast genomes and annotated them. A concatenated data set of all available brown algae chloroplast sequences was used for the calculation of a robust phylogeny, and sequence variations were analyzed. The two Laminaria chloroplast genomes are collinear to previously analyzed kelp chloroplast genomes with important exceptions. Rearrangements at the inverted repeat regions led to the pseudogenization of ycf37 in L. solidungula, a gene possibly required under high light conditions. This defunct gene might be one of the reasons why the habitat range of L. solidungula is restricted to lowlight sublittoral sites in the Arctic. The inheritance pattern of single nucleotide polymorphisms suggests incomplete lineage sorting of chloroplast genomes in kelp species. Our analysis of kelp chloroplast genomes shows that not only evolutionary information could be gleaned from sequence data. Concomitantly, those sequences can also tell us something about the ecological conditions which are required for species well-being.

An aerobic eukaryotic parasite with functional mitochondria that likely lacks a mitochondrial genome.

Dinoflagellates are microbial eukaryotes that have exceptionally large nuclear genomes; however, their organelle genomes are small and fragmented and contain fewer genes than those of other eukaryotes. The genus Amoebophrya (Syndiniales) comprises endoparasites with high genetic diversity that can infect other dinoflagellates, such as those forming harmful algal blooms (e.g., Alexandrium). We sequenced the genome (~100 Mb) of Amoebophrya ceratii to investigate the early evolution of genomic characters in dinoflagellates. The A. ceratii genome encodes almost all essential biosynthetic pathways for self-sustaining cellular metabolism, suggesting a limited dependency on its host. Although dinoflagellates are thought to have descended from a photosynthetic ancestor, A. ceratii appears to have completely lost its plastid and nearly all genes of plastid origin. Functional mitochondria persist in all life stages of A. ceratii, but we found no evidence for the presence of a mitochondrial genome. Instead, all mitochondrial proteins appear to be lost or encoded in the A. ceratii nucleus.

Depletion of Nesprin-2 is associated with an embryonic lethal phenotype in mice.

Nesprin-2 is a nuclear envelope component and provides a link between cytoskeletal components of the cytoplasm and the nucleoplasm. Several isoforms are generated from its gene Syne2. Loss of the largest isoform Nesprin-2 Giant in mice is associated with a skin phenotype and altered wound healing, loss of C-terminal isoforms in mice leads to cardiomyopathies and neurological defects. Here we attempted to establish mice with an inducible knockout of all Nesprin-2 isoforms by inserting shRNA encoding sequences targeting the N- and C-terminus into the ROSA26 locus of mice. This caused early embryonic death of the animals harboring the mutant allele, which was presumably due to leaky expression of the shRNAs. Mutant embryos were only observed before E13. They had an altered appearance and were smaller in size than their wild type littermates. From this we conclude that the Nesprin-2 gene function is crucial during embryonic growth, differentiation and organogenesis.

The Chara Genome: Secondary Complexity and Implications for Plant Terrestrialization.

Land plants evolved from charophytic algae, among which Charophyceae possess the most complex body plans. We present the genome of Chara braunii; comparison of the genome to those of land plants identified evolutionary novelties for plant terrestrialization and land plant heritage genes. C. braunii employs unique xylan synthases for cell wall biosynthesis, a phragmoplast (cell separation) mechanism similar to that of land plants, and many phytohormones. C. braunii plastids are controlled via land-plant-like retrograde signaling, and transcriptional regulation is more elaborate than in other algae. The morphological complexity of this organism may result from expanded gene families, with three cases of particular note: genes effecting tolerance to reactive oxygen species (ROS), LysM receptor-like kinases, and transcription factors (TFs). Transcriptomic analysis of sexual reproductive structures reveals intricate control by TFs, activity of the ROS gene network, and the ancestral use of plant-like storage and stress protection proteins in the zygote.

Multiple Roots of Fruiting Body Formation in Amoebozoa.

Establishment of multicellularity represents a major transition in eukaryote evolution. A subgroup of Amoebozoa, the dictyosteliids, has evolved a relatively simple aggregative multicellular stage resulting in a fruiting body supported by a stalk. Protosteloid amoeba, which are scattered throughout the amoebozoan tree, differ by producing only one or few single stalked spores. Thus, one obvious difference in the developmental cycle of protosteliids and dictyosteliids seems to be the establishment of multicellularity. To separate spore development from multicellular interactions, we compared the genome and transcriptome of a Protostelium species (Protostelium aurantium var. fungivorum) with those of social and solitary members of the Amoebozoa. During fruiting body formation nearly 4,000 genes, corresponding to specific pathways required for differentiation processes, are upregulated. A comparison with genes involved in the development of dictyosteliids revealed conservation of >500 genes, but most of them are also present in Acanthamoeba castellanii for which fruiting bodies have not been documented. Moreover, expression regulation of those genes differs between P. aurantium and Dictyostelium discoideum. Within Amoebozoa differentiation to fruiting bodies is common, but our current genome analysis suggests that protosteliids and dictyosteliids used different routes to achieve this. Most remarkable is both the large repertoire and diversity between species in genes that mediate environmental sensing and signal processing. This likely reflects an immense adaptability of the single cell stage to varying environmental conditions. We surmise that this signaling repertoire provided sufficient building blocks to accommodate the relatively simple demands for cell-cell communication in the early multicellular forms.

Transcriptome reprogramming during developmental switching in Physarum polycephalum involves extensive remodeling of intracellular signaling networks.

Activation of a phytochrome photoreceptor triggers a program of Physarum polycephalum plasmodial cell differentiation through which a mitotic multinucleate protoplasmic mass synchronously develops into haploid spores formed by meiosis and rearrangement of cellular components. We have performed a transcriptome-wide RNAseq study of cellular reprogramming and developmental switching. RNAseq analysis revealed extensive remodeling of intracellular signaling and regulation in switching the expression of sets of genes encoding transcription factors, kinases, phosphatases, signal transduction proteins, RNA-binding proteins, ubiquitin ligases, regulators of the mitotic and meiotic cell cycle etc. in conjunction with the regulation of genes encoding metabolic enzymes and cytoskeletal proteins. About 15% of the differentially expressed genes shared similarity with members of the evolutionary conserved set of core developmental genes of social amoebae. Differential expression of genes encoding regulators that act at the transcriptional, translational, and post-translational level indicates the establishment of a new state of cellular function and reveals evolutionary deeply conserved molecular changes involved in cellular reprogramming and differentiation in a prototypical eukaryote.

The function of the inner nuclear envelope protein SUN1 in mRNA export is regulated by phosphorylation.

SUN1, a component of the LINC (Linker of Nucleoskeleton and Cytoskeleton) complex, functions in mammalian mRNA export through the NXF1-dependent pathway. It associates with mRNP complexes by direct interaction with NXF1. It also binds to the NPC through association with the nuclear pore component Nup153, which is involved in mRNA export. The SUN1-NXF1 association is at least partly regulated by a protein kinase C (PKC) which phosphorylates serine 113 (S113) in the N-terminal domain leading to reduced interaction. The phosphorylation appears to be important for the SUN1 function in nuclear mRNA export since GFP-SUN1 carrying a S113A mutation was less efficient in restoring mRNA export after SUN1 knockdown as compared to the wild type protein. By contrast, GFP-SUN1-S113D resembling the phosphorylated state allowed very efficient export of poly(A)+RNA. Furthermore, probing a possible role of the LINC complex component Nesprin-2 in this process we observed impaired mRNA export in Nesprin-2 knockdown cells. This effect might be independent of SUN1 as expression of a GFP tagged SUN-domain deficient SUN1, which no longer can interact with Nesprin-2, did not affect mRNA export.

A tripeptidyl peptidase 1 is a binding partner of the Golgi pH regulator (GPHR) in Dictyostelium.

Mutations in tripeptidyl peptidase 1 (TPP1) have been associated with late infantile neuronal ceroid lipofuscinosis (NCL), a neurodegenerative disorder. TPP1 is a lysosomal serine protease, which removes tripeptides from the N-terminus of proteins and is composed of an N-terminal prodomain and a catalytic domain. It is conserved in mammals, amphibians, fish and the amoeba Dictyostelium discoideum. D. discoideum harbors at least six genes encoding TPP1, tpp1A to tpp1F We identified TPP1F as binding partner of Dictyostelium GPHR (Golgi pH regulator), which is an evolutionarily highly conserved intracellular transmembrane protein. A region encompassing the DUF3735 (GPHR_N) domain of GPHR was responsible for the interaction. In TPP1F, the binding site is located in the prodomain of the protein. The tpp1F gene is transcribed throughout development and translated into a polypeptide of ∼65 kDa. TPP1 activity was demonstrated for TPP1F-GFP immunoprecipitated from D. discoideum cells. Its activity could be inhibited by addition of the recombinant DUF3735 domain of GPHR. Knockout tpp1F mutants did not display any particular phenotype, and TPP1 activity was not abrogated, presumably because tpp1B compensates as it has the highest expression level of all the TPP1 genes during growth. The GPHR interaction was not restricted to TPP1F but occurred also with TPP1B. As previous reports show that the majority of the TPP1 mutations in NCL resulted in reduction or loss of enzyme activity, we suggest that Dicyostelium could be used as a model system in which to test new reagents that could affect the activity of the protein and ameliorate the disease.

TRE5-A retrotransposition profiling reveals putative RNA polymerase III transcription complex binding sites on the Dictyostelium extrachromosomal rDNA element.

The amoeba Dictyostelium discoideum has a haploid genome in which two thirds of the DNA encodes proteins. Consequently, the space available for selfish mobile elements to expand without excess damage to the host genome is limited. The non-long terminal repeat retrotransposon TRE5-A maintains an active population in the D. discoideum genome and apparently adapted to this gene-dense environment by targeting positions ~47 bp upstream of tRNA genes that are devoid of protein-coding regions. Because only ~24% of tRNA genes are associated with a TRE5-A element in the reference genome, we evaluated whether TRE5-A retrotransposition is limited to this subset of tRNA genes. We determined that a tagged TRE5-A element (TRE5-Absr) integrated at 384 of 405 tRNA genes, suggesting that expansion of the current natural TRE5-A population is not limited by the availability of targets. We further observed that TRE5-Absr targets the ribosomal 5S gene on the multicopy extrachromosomal DNA element that carries the ribosomal RNA genes, indicating that TRE5-A integration may extend to the entire RNA polymerase III (Pol III) transcriptome. We determined that both natural TRE5-A and cloned TRE5-Absr retrotranspose to locations on the extrachromosomal rDNA element that contain tRNA gene-typical A/B box promoter motifs without displaying any other tRNA gene context. Based on previous data suggesting that TRE5-A targets tRNA genes by locating Pol III transcription complexes, we propose that A/B box loci reflect Pol III transcription complex assembly sites that possess a function in the biology of the extrachromosomal rDNA element.

Improved annotation with de novo transcriptome assembly in four social amoeba species.

BACKGROUND: Annotation of gene models and transcripts is a fundamental step in genome sequencing projects. Often this is performed with automated prediction pipelines, which can miss complex and atypical genes or transcripts. RNA sequencing (RNA-seq) data can aid the annotation with empirical data. Here we present de novo transcriptome assemblies generated from RNA-seq data in four Dictyostelid species: D. discoideum, P. pallidum, D. fasciculatum and D. lacteum. The assemblies were incorporated with existing gene models to determine corrections and improvement on a whole-genome scale. This is the first time this has been performed in these eukaryotic species. RESULTS: An initial de novo transcriptome assembly was generated by Trinity for each species and then refined with Program to Assemble Spliced Alignments (PASA). The completeness and quality were assessed with the Benchmarking Universal Single-Copy Orthologs (BUSCO) and Transrate tools at each stage of the assemblies. The final datasets of 11,315-12,849 transcripts contained 5,610-7,712 updates and corrections to >50% of existing gene models including changes to hundreds or thousands of protein products. Putative novel genes are also identified and alternative splice isoforms were observed for the first time in P. pallidum, D. lacteum and D. fasciculatum. CONCLUSIONS: In taking a whole transcriptome approach to genome annotation with empirical data we have been able to enrich the annotations of four existing genome sequencing projects. In doing so we have identified updates to the majority of the gene annotations across all four species under study and found putative novel genes and transcripts which could be worthy for follow-up. The new transcriptome data we present here will be a valuable resource for genome curators in the Dictyostelia and we propose this effective methodology for use in other genome annotation projects.

Evolutionary genomics of the cold-adapted diatom Fragilariopsis cylindrus.

The Southern Ocean houses a diverse and productive community of organisms. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-adapted diatom from the Southern Ocean, Fragilariopsis cylindrus, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO2. Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean.