Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change.

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

18S rDNA gene metabarcoding of microeukaryotes and epi-endophytes in the holobiome of seven species of large brown algae.

Brown algae (Phaeophyceae) are habitat-forming species in coastal ecosystems and include kelp forests and seaweed beds that support a wide diversity of marine life. Host-associated microbial communities are an integral part of phaeophyte biology, and whereas the bacterial microbial partners have received considerable attention, the microbial eukaryotes associated with brown algae have hardly been studied. Here, we used broadly targeted "pan-eukaryotic" primers (metabarcoding) to investigate brown algal-associated eukaryotes (the eukaryome). Using this approach, we aimed to investigate the eukaryome of seven large brown algae that are important and common species in coastal ecosystems. We also aimed to assess whether these macroalgae harbor novel eukaryotic diversity and to ascribe putative functional roles to the host-associated eukaryome based on taxonomic affiliation and phylogenetic placement. We detected a significant diversity of microeukaryotic and algal lineages associated with the brown algal species investigated. The operational taxonomic units (OTUs) were taxonomically assigned to 10 of the eukaryotic major supergroups, including taxonomic groups known to be associated with seaweeds as epibionts, endobionts, parasites, and commensals. Additionally, we revealed previously unrecorded sequence types, including novel phaeophyte OTUs, particularly in the Fucus spp. samples, that may represent fucoid genomic variants, sequencing artifacts, or undescribed epi-/endophytes. Our results provide baseline data and technical insights that will be useful for more comprehensive seaweed eukaryome studies investigating the evidently lineage-rich and functionally diverse symbionts of brown algae.

Structural evidence for MADS-box type I family expansion seen in new assemblies of Arabidopsis arenosa and A. lyrata.

Arabidopsis thaliana diverged from A. arenosa and A. lyrata at least 6 million years ago. The three species differ by genome-wide polymorphisms and morphological traits. The species are to a high degree reproductively isolated, but hybridization barriers are incomplete. A special type of hybridization barrier is based on the triploid endosperm of the seed, where embryo lethality is caused by endosperm failure to support the developing embryo. The MADS-box type I family of transcription factors is specifically expressed in the endosperm and has been proposed to play a role in endosperm-based hybridization barriers. The gene family is well known for its high evolutionary duplication rate, as well as being regulated by genomic imprinting. Here we address MADS-box type I gene family evolution and the role of type I genes in the context of hybridization. Using two de-novo assembled and annotated chromosome-level genomes of A. arenosa and A. lyrata ssp. petraea we analyzed the MADS-box type I gene family in Arabidopsis to predict orthologs, copy number, and structural genomic variation related to the type I loci. Our findings were compared to gene expression profiles sampled before and after the transition to endosperm cellularization in order to investigate the involvement of MADS-box type I loci in endosperm-based hybridization barriers. We observed substantial differences in type-I expression in the endosperm of A. arenosa and A. lyrata ssp. petraea, suggesting a genetic cause for the endosperm-based hybridization barrier between A. arenosa and A. lyrata ssp. petraea.

Disentangling temporal associations in marine microbial networks.

BACKGROUND: Microbial interactions are fundamental for Earth's ecosystem functioning and biogeochemical cycling. Nevertheless, they are challenging to identify and remain barely known. Omics-based censuses are helpful in predicting microbial interactions through the statistical inference of single (static) association networks. Yet, microbial interactions are dynamic and we have limited knowledge of how they change over time. Here, we investigate the dynamics of microbial associations in a 10-year marine time series in the Mediterranean Sea using an approach inferring a time-resolved (temporal) network from a single static network. RESULTS: A single static network including microbial eukaryotes and bacteria was built using metabarcoding data derived from 120 monthly samples. For the decade, we aimed to identify persistent, seasonal, and temporary microbial associations by determining a temporal network that captures the interactome of each individual sample. We found that the temporal network appears to follow an annual cycle, collapsing, and reassembling when transiting between colder and warmer waters. We observed higher association repeatability in colder than in warmer months. Only 16 associations could be validated using observations reported in literature, underlining our knowledge gap in marine microbial ecological interactions. CONCLUSIONS: Our results indicate that marine microbial associations follow recurrent temporal dynamics in temperate zones, which need to be accounted for to better understand the functioning of the ocean microbiome. The constructed marine temporal network may serve as a resource for testing season-specific microbial interaction hypotheses. The applied approach can be transferred to microbiome studies in other ecosystems. Video Abstract.

Lymphocyte subsets in Atlantic cod (Gadus morhua) interrogated by single-cell sequencing.

Atlantic Cod (Gadus morhua) has lost the major histocompatibility complex class II presentation pathway. We recently identified CD8-positive T cells, B cells, and plasma cells in cod, but further characterisation of lymphocyte subsets is needed to elucidate immune adaptations triggered by the absence of CD4-positive T lymphocytes. Here, we use single-cell RNA sequencing to examine the lymphocyte heterogeneity in Atlantic cod spleen. We describe five T cell subsets and eight B cell subsets and propose a B cell trajectory of differentiation. Notably, we identify a subpopulation of T cells that are CD8-negative. Most of the CD8-negative T lymphocytes highly express the homologue of monocyte chemotactic protein 1b, and another subset of CD8-negative T lymphocytes express the homologue of the scavenger receptor m130. Uncovering the multiple lymphocyte cell sub-clusters reveals the different immune states present within the B and T cell populations, building a foundation for further work.

Long-term patterns of an interconnected core marine microbiota.

BACKGROUND: Ocean microbes constitute ~ 70% of the marine biomass, are responsible for ~ 50% of the Earth's primary production and are crucial for global biogeochemical cycles. Marine microbiotas include core taxa that are usually key for ecosystem function. Despite their importance, core marine microbes are relatively unknown, which reflects the lack of consensus on how to identify them. So far, most core microbiotas have been defined based on species occurrence and abundance. Yet, species interactions are also important to identify core microbes, as communities include interacting species. Here, we investigate interconnected bacteria and small protists of the core pelagic microbiota populating a long-term marine-coastal observatory in the Mediterranean Sea over a decade. RESULTS: Core microbes were defined as those present in > 30% of the monthly samples over 10 years, with the strongest associations. The core microbiota included 259 Operational Taxonomic Units (OTUs) including 182 bacteria, 77 protists, and 1411 strong and mostly positive (~ 95%) associations. Core bacteria tended to be associated with other bacteria, while core protists tended to be associated with bacteria. The richness and abundance of core OTUs varied annually, decreasing in stratified warmers waters and increasing in colder mixed waters. Most core OTUs had a preference for one season, mostly winter, which featured subnetworks with the highest connectivity. Groups of highly associated taxa tended to include protists and bacteria with predominance in the same season, particularly winter. A group of 13 highly-connected hub-OTUs, with potentially important ecological roles dominated in winter and spring. Similarly, 18 connector OTUs with a low degree but high centrality were mostly associated with summer or autumn and may represent transitions between seasonal communities. CONCLUSIONS: We found a relatively small and dynamic interconnected core microbiota in a model temperate marine-coastal site, with potential interactions being more deterministic in winter than in other seasons. These core microbes would be essential for the functioning of this ecosystem over the year. Other non-core taxa may also carry out important functions but would be redundant and non-essential. Our work contributes to the understanding of the dynamics and potential interactions of core microbes possibly sustaining ocean ecosystem function.

Disentangling environmental effects in microbial association networks.

BACKGROUND: Ecological interactions among microorganisms are fundamental for ecosystem function, yet they are mostly unknown or poorly understood. High-throughput-omics can indicate microbial interactions through associations across time and space, which can be represented as association networks. Associations could result from either ecological interactions between microorganisms, or from environmental selection, where the association is environmentally driven. Therefore, before downstream analysis and interpretation, we need to distinguish the nature of the association, particularly if it is due to environmental selection or not. RESULTS: We present EnDED (environmentally driven edge detection), an implementation of four approaches as well as their combination to predict which links between microorganisms in an association network are environmentally driven. The four approaches are sign pattern, overlap, interaction information, and data processing inequality. We tested EnDED on networks from simulated data of 50 microorganisms. The networks contained on average 50 nodes and 1087 edges, of which 60 were true interactions but 1026 false associations (i.e., environmentally driven or due to chance). Applying each method individually, we detected a moderate to high number of environmentally driven edges-87% sign pattern and overlap, 67% interaction information, and 44% data processing inequality. Combining these methods in an intersection approach resulted in retaining more interactions, both true and false (32% of environmentally driven associations). After validation with the simulated datasets, we applied EnDED on a marine microbial network inferred from 10 years of monthly observations of microbial-plankton abundance. The intersection combination predicted that 8.3% of the associations were environmentally driven, while individual methods predicted 24.8% (data processing inequality), 25.7% (interaction information), and up to 84.6% (sign pattern as well as overlap). The fraction of environmentally driven edges among negative microbial associations in the real network increased rapidly with the number of environmental factors. CONCLUSIONS: To reach accurate hypotheses about ecological interactions, it is important to determine, quantify, and remove environmentally driven associations in marine microbial association networks. For that, EnDED offers up to four individual methods as well as their combination. However, especially for the intersection combination, we suggest using EnDED with other strategies to reduce the number of false associations and consequently the number of potential interaction hypotheses. Video abstract.

The sequenced genomes of nonflowering land plants reveal the innovative evolutionary history of peptide signaling.

An understanding of land plant evolution is a prerequisite for in-depth knowledge of plant biology. Here we extract and explore information hidden in the increasing number of sequenced plant genomes, from bryophytes to angiosperms, to elucidate a specific biological question-how peptide signaling evolved. To conquer land and cope with changing environmental conditions, plants have gone through transformations that must have required innovations in cell-to-cell communication. We discuss peptides mediating endogenous and exogenous changes by interaction with receptors activating intracellular molecular signaling. Signaling peptides were discovered in angiosperms and operate in tissues and organs such as flowers, seeds, vasculature, and 3D meristems that are not universally conserved across land plants. Nevertheless, orthologs of angiosperm peptides and receptors have been identified in nonangiosperms. These discoveries provoke questions regarding coevolution of ligands and their receptors, and whether de novo interactions in peptide signaling pathways may have contributed to generate novel traits in land plants. The answers to such questions will have profound implications for the understanding of the evolution of cell-to-cell communication and the wealth of diversified terrestrial plants. Under this perspective, we have generated, analyzed, and reviewed phylogenetic, genomic, structural, and functional data to elucidate the evolution of peptide signaling.

Niche adaptation promoted the evolutionary diversification of tiny ocean predators.

Unicellular eukaryotic predators play a crucial role in the functioning of the ocean ecosystem by recycling nutrients and energy that are channeled to upper trophic levels. Traditionally, these evolutionarily diverse organisms have been combined into a single functional group (heterotrophic flagellates), overlooking their organismal differences. Here, we investigated four evolutionarily related species belonging to one cosmopolitan group of uncultured marine picoeukaryotic predators: marine stramenopiles (MAST)-4 (species A, B, C, and E). Co-occurrence and distribution analyses in the global surface ocean indicated contrasting patterns in MAST-4A and C, suggesting adaptation to different temperatures. We then investigated whether these spatial distribution patterns were mirrored by MAST-4 genomic content using single-cell genomics. Analyses of 69 single cells recovered 66 to 83% of the MAST-4A/B/C/E genomes, which displayed substantial interspecies divergence. MAST-4 genomes were similar in terms of broad gene functional categories, but they differed in enzymes of ecological relevance, such as glycoside hydrolases (GHs), which are part of the food degradation machinery in MAST-4. Interestingly, MAST-4 species featuring a similar GH composition (A and C) coexcluded each other in the surface global ocean, while species with a different set of GHs (B and C) appeared to be able to coexist, suggesting further niche diversification associated with prey digestion. We propose that differential niche adaptation to temperature and prey type has promoted adaptive evolutionary diversification in MAST-4. We show that minute ocean predators from the same phylogenetic group may have different biogeography and genomic content, which needs to be accounted for to better comprehend marine food webs.

Genome sequencing and de novo assembly of the giant unicellular alga Acetabularia acetabulum using droplet MDA.

The macroscopic single-celled green alga Acetabularia acetabulum has been a model system in cell biology for more than a century. However, no genomic information is available from this species. Since the alga has a long life cycle, is difficult to grow in dense cultures, and has an estimated diploid genome size of almost 2 Gb, obtaining sufficient genomic material for genome sequencing is challenging. Here, we have attempted to overcome these challenges by amplifying genomic DNA using multiple displacement amplification (MDA) combined with microfluidics technology to distribute the amplification reactions across thousands of microscopic droplets. By amplifying and sequencing DNA from five single cells we were able to recover an estimated ~ 7-11% of the total genome, providing the first draft of the A. acetabulum genome. We highlight challenges associated with genome recovery and assembly of MDA data due to biases arising during genome amplification, and hope that our study can serve as a reference for future attempts on sequencing the genome from non-model eukaryotes.

Disentangling the mechanisms shaping the surface ocean microbiota.

BACKGROUND: The ocean microbiota modulates global biogeochemical cycles and changes in its configuration may have large-scale consequences. Yet, the underlying ecological mechanisms structuring it are unclear. Here, we investigate how fundamental ecological mechanisms (selection, dispersal and ecological drift) shape the smallest members of the tropical and subtropical surface-ocean microbiota: prokaryotes and minute eukaryotes (picoeukaryotes). Furthermore, we investigate the agents exerting abiotic selection on this assemblage as well as the spatial patterns emerging from the action of ecological mechanisms. To explore this, we analysed the composition of surface-ocean prokaryotic and picoeukaryotic communities using DNA-sequence data (16S- and 18S-rRNA genes) collected during the circumglobal expeditions Malaspina-2010 and TARA-Oceans. RESULTS: We found that the two main components of the tropical and subtropical surface-ocean microbiota, prokaryotes and picoeukaryotes, appear to be structured by different ecological mechanisms. Picoeukaryotic communities were predominantly structured by dispersal-limitation, while prokaryotic counterparts appeared to be shaped by the combined action of dispersal-limitation, selection and drift. Temperature-driven selection appeared as a major factor, out of a few selected factors, influencing species co-occurrence networks in prokaryotes but not in picoeukaryotes, indicating that association patterns may contribute to understand ocean microbiota structure and response to selection. Other measured abiotic variables seemed to have limited selective effects on community structure in the tropical and subtropical ocean. Picoeukaryotes displayed a higher spatial differentiation between communities and a higher distance decay when compared to prokaryotes, consistent with a scenario of higher dispersal limitation in the former after considering environmental heterogeneity. Lastly, random dynamics or drift seemed to have a more important role in structuring prokaryotic communities than picoeukaryotic counterparts. CONCLUSIONS: The differential action of ecological mechanisms seems to cause contrasting biogeography, in the tropical and subtropical ocean, among the smallest surface plankton, prokaryotes and picoeukaryotes. This suggests that the idiosyncrasy of the main constituents of the ocean microbiota should be considered in order to understand its current and future configuration, which is especially relevant in a context of global change, where the reaction of surface ocean plankton to temperature increase is still unclear. Video Abstract.

Genetic variation and temperature affects hybrid barriers during interspecific hybridization.

Genomic imprinting regulates parent-specific transcript dosage during seed development and is mainly confined to the endosperm. Elucidation of the function of many imprinted genes has been hampered by the lack of corresponding mutant phenotypes, and the role of imprinting is mainly associated with genome dosage regulation or allocation of resources. Disruption of imprinted genes has also been suggested to mediate endosperm-based post-zygotic hybrid barriers depending on genetic variation and gene dosage. Here, we have analyzed the conservation of a clade from the MADS-box type I class transcription factors in the closely related species Arabidopsis arenosa, A. lyrata, and A. thaliana, and show that AGL36-like genes are imprinted and maternally expressed in seeds of Arabidopsis species and in hybrid seeds between outbreeding species. In hybridizations between outbreeding and inbreeding species the paternally silenced allele of the AGL36-like gene is reactivated in the hybrid, demonstrating that also maternally expressed imprinted genes are perturbed during hybridization and that such effects on imprinted genes are specific to the species combination. Furthermore, we also demonstrate a quantitative effect of genetic diversity and temperature on the strength of the post-zygotic hybridization barrier. Markedly, a small decrease in temperature during seed development increases the survival of hybrid F1 seeds, suggesting that abiotic and genetic parameters play important roles in post-zygotic species barriers, pointing at evolutionary scenarios favoring such effects. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA562212. All sequences generated in this study have been deposited in the National Center for Biotechnology Information Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra/) with project number PRJNA562212.

The planktonic protist interactome: where do we stand after a century of research?

Microbial interactions are crucial for Earth ecosystem function, but our knowledge about them is limited and has so far mainly existed as scattered records. Here, we have surveyed the literature involving planktonic protist interactions and gathered the information in a manually curated Protist Interaction DAtabase (PIDA). In total, we have registered ~2500 ecological interactions from ~500 publications, spanning the last 150 years. All major protistan lineages were involved in interactions as hosts, symbionts (mutualists and commensalists), parasites, predators, and/or prey. Predation was the most common interaction (39% of all records), followed by symbiosis (29%), parasitism (18%), and 'unresolved interactions' (14%, where it is uncertain whether the interaction is beneficial or antagonistic). Using bipartite networks, we found that protist predators seem to be 'multivorous' while parasite-host and symbiont-host interactions appear to have moderate degrees of specialization. The SAR supergroup (i.e., Stramenopiles, Alveolata, and Rhizaria) heavily dominated PIDA, and comparisons against a global-ocean molecular survey (TARA Oceans) indicated that several SAR lineages, which are abundant and diverse in the marine realm, were underrepresented among the recorded interactions. Despite historical biases, our work not only unveils large-scale eco-evolutionary trends in the protist interactome, but it also constitutes an expandable resource to investigate protist interactions and to test hypotheses deriving from omics tools.

Seasonal Dynamics of Algae-Infecting Viruses and Their Inferred Interactions with Protists.

Viruses are a highly abundant, dynamic, and diverse component of planktonic communities that have key roles in marine ecosystems. We aimed to reveal the diversity and dynamics of marine large dsDNA viruses infecting algae in the Northern Skagerrak, South Norway through the year by metabarcoding, targeting the major capsid protein (MCP) and its correlation to protist diversity and dynamics. Metabarcoding results demonstrated a high diversity of algal viruses compared to previous metabarcoding surveys in Norwegian coastal waters. We obtained 313 putative algal virus operational taxonomic units (vOTUs), all classified by phylogenetic analyses to either the Phycodnaviridae or Mimiviridae families, most of them in clades without any cultured or environmental reference sequences. The viral community showed a clear temporal variation, with some vOTUs persisting for several months. The results indicate co-occurrences between abundant viruses and potential hosts during long periods. This study gives new insights into the virus-algal host dynamics and provides a baseline for future studies of algal virus diversity and temporal dynamics.

Quantifying long-term recurrence in planktonic microbial eukaryotes.

How much temporal recurrence is present in microbial assemblages is still an unanswered ecological question. Even though marked seasonal changes have been reported for whole microbial communities, less is known on the dynamics and seasonality of individual taxa. Here, we aim at understanding microbial recurrence at three different levels: community, taxonomic group and operational taxonomic units (OTUs). For that, we focused on a model microbial eukaryotic community populating a long-term marine microbial observatory using 18S rRNA gene data from two organismal size fractions: the picoplankton (0.2-3 µm) and the nanoplankton (3-20 µm). We have developed an index to quantify recurrence in particular taxa. We found that community structure oscillated systematically between two main configurations corresponding to winter and summer over the 10 years studied. A few taxonomic groups such as Mamiellophyceae or MALV-III presented clear recurrence (i.e., seasonality), whereas 13%-19% of the OTUs in both size fractions, accounting for ~40% of the relative abundance, featured recurrent dynamics. Altogether, our work links long-term whole community dynamics with that of individual OTUs and taxonomic groups, indicating that recurrent and non-recurrent changes characterize the dynamics of microbial assemblages.

Single Cell Transcriptomics, Mega-Phylogeny, and the Genetic Basis of Morphological Innovations in Rhizaria.

The innovation of the eukaryote cytoskeleton enabled phagocytosis, intracellular transport, and cytokinesis, and is largely responsible for the diversity of morphologies among eukaryotes. Still, the relationship between phenotypic innovations in the cytoskeleton and their underlying genotype is poorly understood. To explore the genetic mechanism of morphological evolution of the eukaryotic cytoskeleton, we provide the first single cell transcriptomes from uncultured, free-living unicellular eukaryotes: the polycystine radiolarian Lithomelissa setosa (Nassellaria) and Sticholonche zanclea (Taxopodida). A phylogenomic approach using 255 genes finds Radiolaria and Foraminifera as separate monophyletic groups (together as Retaria), while Cercozoa is shown to be paraphyletic where Endomyxa is sister to Retaria. Analysis of the genetic components of the cytoskeleton and mapping of the evolution of these on the revised phylogeny of Rhizaria reveal lineage-specific gene duplications and neofunctionalization of α and ß tubulin in Retaria, actin in Retaria and Endomyxa, and Arp2/3 complex genes in Chlorarachniophyta. We show how genetic innovations have shaped cytoskeletal structures in Rhizaria, and how single cell transcriptomics can be applied for resolving deep phylogenies and studying gene evolution in uncultured protist species.

BIR Pipeline for Preparation of Phylogenomic Data.

SUMMARY: We present a pipeline named BIR (Blast, Identify and Realign) developed for phylogenomic analyses. BIR is intended for the identification of gene sequences applicable for phylogenomic inference. The pipeline allows users to apply their own manually curated sequence alignments (seed) in search for homologous genes in sequence databases and available genomes. BIR automatically adds the identified sequences from these databases to the seed alignments and reconstruct a phylogenetic tree from each. The BIR pipeline is an efficient tool for the identification of orthologous gene copies because it expands user-defined sequence alignments and conducts massive parallel phylogenetic reconstruction. The application is also particularly useful for large-scale sequencing projects that require management of a large number of single-gene alignments for gene comparison, functional annotation, and evolutionary analyses. AVAILABILITY: The BIR user manual is available at http://www.bioportal.no/ and can be accessed through Lifeportal at https://lifeportal.uio.no. Access is free but requires a user account registration using the link "Register for BIR access" from the Lifeportal homepage.

Radiolaria associated with large diversity of marine alveolates.

We have isolated cells of unculturable radiolarians from marine coastal waters. Individual cells were subjected to single cell whole genome amplification (SCWGA) and gene-targeted PCR. Using this approach we recover a surprisingly large diversity of sequences related to the enigmatic marine alveolate groups 1 and 2 (MALV I and MALV II) that most likely represent intracellular symbionts or parasites of the radiolarian cells. 18S rDNA phylogeny of the MALV sequences reveals 4 distinct clades of radiolarian associates here named Radiolarian Associated Sequences (RAS) 1-4. One clade of both phaeodarian and radiolarian associates and one clade of only phaeodarian associates are also identified. The MALV sequences cluster according to host type, i.e. sequences from associates identified in radiolarians, fish, copepods, ciliates or dinoflagellates are not intermixed but separated into distinct clades. This implies several independent colonizations of host lineages and links a large diversity of MALV to radiolarian-associated species. This demonstrates that radiolarians may be an important reservoir for MALV, making them a key group for understanding the impact of intracellular symbionts on the marine ecosystem. This study shows that applying SCWGA on unculturable cells is a promising approach to study the vast diversity and interactions of intracellular eukaryote organisms.

Radiolaria divided into Polycystina and Spasmaria in combined 18S and 28S rDNA phylogeny.

Radiolarians are marine planktonic protists that belong to the eukaryote supergroup Rhizaria together with Foraminifera and Cercozoa. Radiolaria has traditionally been divided into four main groups based on morphological characters; i.e. Polycystina, Acantharia, Nassellaria and Phaeodaria. But recent 18S rDNA phylogenies have shown that Phaeodaria belongs within Cerocozoa, and that the previously heliozoan group Taxopodida should be included in Radiolaria. 18S rDNA phylogenies have not yet resolved the sister relationship between the main Radiolaria groups, but nevertheless suggests that Spumellaria, and thereby also Polycystina, are polyphyletic. Very few sequences other than 18S rDNA have so far been generated from radiolarian cells, mostly due to the fact that Radiolaria has been impossible to cultivate and single cell PCR has been hampered by low success rate. Here we have therefore investigated the mutual evolutionary relationship of the main radiolarian groups by using the novel approach of combining single cell whole genome amplification with targeted PCR amplification of the 18S and 28S rDNA genes. Combined 18S and 28S phylogeny of sequences obtained from single cells shows that Radiolaria is divided into two main lineages: Polycystina (Spumellaria+Nassellaria) and Spasmaria (Acantharia+Taxopodida). Further we show with high support that Foraminifera groups within Radiolaria supporting the Retaria hypothesis.