Global Trends in Marine Plankton Diversity across Kingdoms of Life.

The ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation. VIDEO ABSTRACT.

Species detection and delineation in the marine planktonic diatoms Chaetoceros and Bacteriastrum through metabarcoding: making biological sense of haplotype diversity.

High-throughput sequencing (HTS) metabarcoding is commonly applied to assess phytoplankton diversity. Usually, haplotypes are grouped into operational taxonomic units (OTUs) through clustering, whereby the resulting number of OTUs depends on chosen similarity thresholds. We applied, instead, a phylogenetic approach to infer taxa among 18S rDNA V4-metabarcode haplotypes gathered from 48 time-series samples using the marine planktonic diatoms Chaetoceros and Bacteriastrum as test case. The 73 recovered taxa comprised both solitary haplotypes and polytomies, the latter composed each of a highly abundant, dominant haplotype and one to several minor, peripheral haplotypes. The solitary and dominant haplotypes usually matched reference sequences, enabling species assignation of taxa. We hypothesise that the super-abundance of reads in dominant haplotypes results from the homogenization effect of concerted evolution. Reads of populous peripheral haplotypes and dominant haplotypes show comparable distribution patterns over the sample dates, suggesting that they are part of the same population. Many taxa revealed marked seasonality, with closely related ones generally showing distinct periodicity, whereas others occur year-round. Phylogenies inferred from metabarcode haplotypes enable delineation of biologically meaningful taxa, whereas OTUs resulting from clustering algorithms often deviate markedly from such taxa.

Relationship between West African ancestry with lung cancer risk and survival in African Americans.

PURPOSE: African Americans, especially men, have a higher incidence of lung cancer compared with all other racial and ethnic groups in the US. Self-reported race is frequently used in genomic research studies to capture an individual's race or ethnicity. However, it is clear from studies of genetic admixture that human genetic variation does not segregate into the same biologically discrete categories as socially defined categories of race. Previous studies have suggested that the degree of West African ancestry among African Americans can contribute to cancer risk in this population, though few studies have addressed this question in lung cancer. METHODS: Using a genetic ancestry panel of 100 SNPs, we estimated West African, European, and Native American ancestry in 1,407 self-described African Americans and 2,413 European Americans. RESULTS: We found that increasing West African ancestry was associated with increased risk of lung cancer among African American men (ORQ5 vs Q1 = 2.55 (1.45-4.48), p = 0.001), while no association was observed in African American women (ORQ5 vs Q1 = 0.90 (0.51-1.59), p = 0.56). This relationship diminished following adjustment for income and education. CONCLUSIONS: Genetic ancestry is not a major contributor to lung cancer risk or survival disparities.

Green and golden seaweed tides on the rise.

Sudden beaching of huge seaweed masses smother the coastline and form rotting piles on the shore. The number of reports of these events in previously unaffected areas has increased worldwide in recent years. These 'seaweed tides' can harm tourism-based economies, smother aquaculture operations or disrupt traditional artisanal fisheries. Coastal eutrophication is the obvious, ultimate explanation for the increase in seaweed biomass, but the proximate processes that are responsible for individual beaching events are complex and require dedicated study to develop effective mitigation strategies. Harvesting the macroalgae, a valuable raw material, before they beach could well be developed into an effective solution.

Biological Effects of the Azaspiracid-Producing Dinoflagellate Azadinium dexteroporum in Mytilus galloprovincialis from the Mediterranean Sea.

Azaspiracids (AZAs) are marine biotoxins including a variety of analogues. Recently, novel AZAs produced by the Mediterranean dinoflagellate Azadinium dexteroporum were discovered (AZA-54, AZA-55, 3-epi-AZA-7, AZA-56, AZA-57 and AZA-58) and their biological effects have not been investigated yet. This study aimed to identify the biological responses (biomarkers) induced in mussels Mytilus galloprovincialis after the bioaccumulation of AZAs from A. dexteroporum. Organisms were fed with A. dexteroporum for 21 days and subsequently subjected to a recovery period (normal diet) of 21 days. Exposed organisms accumulated AZA-54, 3-epi-AZA-7 and AZA-55, predominantly in the digestive gland. Mussels' haemocytes showed inhibition of phagocytosis activity, modulation of the composition of haemocytic subpopulation and damage to lysosomal membranes; the digestive tissue displayed thinned tubule walls, consumption of storage lipids and accumulation of lipofuscin. Slight genotoxic damage was also observed. No clear occurrence of oxidative stress and alteration of nervous activity was detected in AZA-accumulating mussels. Most of the altered parameters returned to control levels after the recovery phase. The toxic effects detected in M. galloprovincialis demonstrate a clear biological impact of the AZAs produced by A. dexteroporum, and could be used as early indicators of contamination associated with the ingestion of seafood.

Insights into global diatom distribution and diversity in the world's ocean.

Diatoms (Bacillariophyta) constitute one of the most diverse and ecologically important groups of phytoplankton. They are considered to be particularly important in nutrient-rich coastal ecosystems and at high latitudes, but considerably less so in the oligotrophic open ocean. The Tara Oceans circumnavigation collected samples from a wide range of oceanic regions using a standardized sampling procedure. Here, a total of ∼12 million diatom V9-18S ribosomal DNA (rDNA) ribotypes, derived from 293 size-fractionated plankton communities collected at 46 sampling sites across the global ocean euphotic zone, have been analyzed to explore diatom global diversity and community composition. We provide a new estimate of diversity of marine planktonic diatoms at 4,748 operational taxonomic units (OTUs). Based on the total assigned ribotypes, Chaetoceros was the most abundant and diverse genus, followed by Fragilariopsis, Thalassiosira, and Corethron We found only a few cosmopolitan ribotypes displaying an even distribution across stations and high abundance, many of which could not be assigned with confidence to any known genus. Three distinct communities from South Pacific, Mediterranean, and Southern Ocean waters were identified that share a substantial percentage of ribotypes within them. Sudden drops in diversity were observed at Cape Agulhas, which separates the Indian and Atlantic Oceans, and across the Drake Passage between the Atlantic and Southern Oceans, indicating the importance of these ocean circulation choke points in constraining diatom distribution and diversity. We also observed high diatom diversity in the open ocean, suggesting that diatoms may be more relevant in these oceanic systems than generally considered.

Comparison of coastal phytoplankton composition estimated from the V4 and V9 regions of the 18S rRNA gene with a focus on photosynthetic groups and especially Chlorophyta.

We compared the composition of eukaryotic communities using two genetic markers (18S rRNA V4 and V9 regions) at 27 sites sampled during Ocean Sampling Day 2014, with a focus on photosynthetic groups and, more specifically green algae (Chlorophyta). Globally, the V4 and V9 regions of the 18S rRNA gene provided similar images of alpha diversity and ecological patterns. However, V9 provided 20% more OTUs built at 97% identity than V4. 34% of the genera were found with both markers and, of the remnant, 22% were found only with V4 and 44% only with V9. For photosynthetic groups, V4 and V9 performed equally well to describe global communities at different taxonomic levels from the division to the genus and provided similar Chlorophyta distribution patterns. However, at lower taxonomic level, the V9 dataset failed for example to describe the diversity of Dolichomastigales (Chlorophyta, Mamiellophyceae) emphasizing the lack of V9 sequences for this group and the importance of the reference database for metabarcode analysis. We conclude that in order to address questions regarding specific groups (e.g., a given genus), it is necessary to choose the marker based not only on the genetic divergence within this group but also on the existence of reference sequences in databases.

Diatom flagellar genes and their expression during sexual reproduction in Leptocylindrus danicus.

BACKGROUND: Flagella have been lost in the vegetative phase of the diatom life cycle, but they are still present in male gametes of centric species, thereby representing a hallmark of sexual reproduction. This process, besides maintaining and creating new genetic diversity, in diatoms is also fundamental to restore the maximum cell size following its reduction during vegetative division. Nevertheless, sexual reproduction has been demonstrated in a limited number of diatom species, while our understanding of its different phases and of their genetic control is scarce. RESULTS: In the transcriptome of Leptocylindrus danicus, a centric diatom widespread in the world's seas, we identified 22 transcripts related to the flagella development and confirmed synchronous overexpression of 6 flagellum-related genes during the male gamete formation process. These transcripts were mostly absent in the closely related species L. aporus, which does not have sexual reproduction. Among the 22 transcripts, L. danicus showed proteins that belong to the Intra Flagellar Transport (IFT) subcomplex B as well as IFT-A proteins, the latter previously thought to be absent in diatoms. The presence of flagellum-related proteins was also traced in the transcriptomes of several other centric species. Finally, phylogenetic reconstruction of the IFT172 and IFT88 proteins showed that their sequences are conserved across protist species and have evolved similarly to other phylogenetic marker genes. CONCLUSION: Our analysis describes for the first time the diatom flagellar gene set, which appears to be more complete and functional than previously reported based on the genome sequence of the model centric diatom, Thalassiosira pseudonana. This first recognition of the whole set of diatom flagellar genes and of their activation pattern paves the way to a wider recognition of the relevance of sexual reproduction in individual species and in the natural environment.

Mediterranean Azadinium dexteroporum (Dinophyceae) produces six novel azaspiracids and azaspiracid-35: a structural study by a multi-platform mass spectrometry approach.

Azadinium dexteroporum is the first species of the genus described from the Mediterranean Sea and it produces different azaspiracids (AZA). The aims of this work were to characterize the toxin profile of the species and gain structural information on azaspiracids produced by the A. dexteroporum strain SZN-B848 isolated from the Gulf of Naples. Liquid chromatography-mass spectrometry (LC-MS) analyses were carried out on three MS systems having different ion source geometries (ESI, TurboIonSpray®, ESI ION MAX) and different MS analyzers operating either at unit resolution or at high resolution, namely a hybrid triple quadrupole-linear ion trap (Q-Trap MS), a time of flight (TOF MS), and a hybrid linear ion trap Orbitrap XL Fourier transform mass spectrometer (LTQ Orbitrap XL FTMS). As a combined result of these different analyses, A. dexteroporum showed to produce AZA-35, previously reported from Azadinium spinosum, and six compounds that represent new additions to the AZA-group of toxins, including AZA-54 to AZA-58 and 3-epiAZA-7, a stereoisomer of the shellfish metabolite AZA-7. Based on the interpretation of fragmentation patterns, we propose that all these molecules, except AZA-55, have the same A to I ring system as AZA-1, with structural modifications all located in the carboxylic side chain. Considering that none of the azaspiracids produced by the Mediterranean strain of A. dexteroporum is currently regulated by European food safety authorities, monitoring programs of marine biotoxins in the Mediterranean area should take into account the occurrence of the new analogues to avoid an underestimation of the AZA-related risk for seafood consumers. Graphical Abstract A multi-platform MS approach reveals known and new azaspiracids in a Mediterranean strain of Azadinium dexteroporum.

Monoclonal gammopathy-associated pure red cell aplasia.

Pure red cell aplasia (PRCA) is a rare disorder characterized by inhibition of erythroid precursors in the bone marrow and normochromic, normocytic anaemia with reticulocytopenia. Among 51 PRCA patients, we identified 12 (24%) patients having monoclonal gammopathy, monoclonal gammopathy of undetermined significance or smouldering multiple myeloma, with presence of monoclonal protein or abnormal serum free light chains and atypical bone marrow features of clonal plasmacytosis, hypercellularity and fibrosis. Thus far, three patients treated with anti-myeloma based therapeutics have responded with reticulocyte recovery and clinical transfusion independence, suggesting plasma cells play a key role in the pathogenesis of this specific monoclonal gammopathy-associated PRCA.

Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing.

Although protists are critical components of marine ecosystems, they are still poorly characterized. Here we analysed the taxonomic diversity of planktonic and benthic protist communities collected in six distant European coastal sites. Environmental deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) from three size fractions (pico-, nano- and micro/mesoplankton), as well as from dissolved DNA and surface sediments were used as templates for tag pyrosequencing of the V4 region of the 18S ribosomal DNA. Beta-diversity analyses split the protist community structure into three main clusters: picoplankton-nanoplankton-dissolved DNA, micro/mesoplankton and sediments. Within each cluster, protist communities from the same site and time clustered together, while communities from the same site but different seasons were unrelated. Both DNA and RNA-based surveys provided similar relative abundances for most class-level taxonomic groups. Yet, particular groups were overrepresented in one of the two templates, such as marine alveolates (MALV)-I and MALV-II that were much more abundant in DNA surveys. Overall, the groups displaying the highest relative contribution were Dinophyceae, Diatomea, Ciliophora and Acantharia. Also, well represented were Mamiellophyceae, Cryptomonadales, marine alveolates and marine stramenopiles in the picoplankton, and Monadofilosa and basal Fungi in sediments. Our extensive and systematic sequencing of geographically separated sites provides the most comprehensive molecular description of coastal marine protist diversity to date.

Altered cytokine and chemokine profiles in multiple myeloma and its precursor disease.

Currently, no reliable biomarkers are available to predict transformation from smoldering myeloma (SMM) to multiple myeloma (MM). Using an ultrasensitive enzyme-linked immunosorbent assay (ELISA) we assessed the levels of a broad range of cytokines and chemokines in the peripheral blood (PB) and bone marrow (BM) supernatant collected from 14 SMM and 38 MM patients and compared to healthy donors. We found significantly increased levels of key cytokines, in particular CXCL8 (IL-8), associated with progressive disease state (controls→SMM→MM). Cytokine profiles were found similar in PB and BM. Five of fourteen SMM patients (36%) progressed to MM. Our findings, although based on a limited number of patients, suggest that serum-based cytokines may have a future role as biomarkers for disease progression and could potentially be assessed as novel targets for treatment.

Oxylipin diversity in the diatom family Leptocylindraceae reveals DHA derivatives in marine diatoms.

Marine planktonic organisms, such as diatoms, are prospective sources of novel bioactive metabolites. Oxygenated derivatives of fatty acids, generally referred to as oxylipins, in diatoms comprise a highly diverse and complex family of secondary metabolites. These molecules have recently been implicated in several biological processes including intra- and inter-cellular signaling as well as in defense against biotic stressors and grazers. Here, we analyze the production and diversity of C20 and C22 non-volatile oxylipins in five species of the family Leptocylindraceae, which constitute a basal clade in the diatom phylogeny. We report the presence of species-specific lipoxygenase activity and oxylipin patterns, providing the first demonstration of enzymatic production of docosahexaenoic acid derivatives in marine diatoms. The differences observed in lipoxygenase pathways among the species investigated broadly reflected the relationships observed with phylogenetic markers, thus providing functional support to the taxonomic diversity of the individual species.

DRD1 suppresses cell proliferation and reduces EGFR activation and PD-L1 expression in NSCLC.

Dopamine (DA) acts in various key neurological and physiological processes as both a neurotransmitter and circulating hormone. Over the past several decades, the DA signaling network has been shown to regulate the progression of several types of solid tumors, and considerable evidence has shown it is a druggable pathway in the cancer cell context. However, the specific activity and effect of these pathway components appears to be tissue-type and cell-context-dependent. In the present study, expression and methylation of dopamine receptor D1 (DRD1) were measured using RNA sequencing (RNAseq) and reverse transcription polymerase chain reaction (RT-PCR) in non-small cell lung cancer (NSCLC) samples, and validated using publicly available datasets, including The Cancer Genome Atlas (TCGA). In vitro and in vivo functional experiments were performed for cell proliferation and tumor growth, respectively. Mechanistic analyses of the transcriptome and kinome in DRD1-modulated cells informed further experiments, which characterized the effects on the epidermal growth factor receptor (EGFR) pathway and programmed cell death 1 ligand 1 (PD-L1) proteins. Through these experiments, we identified the DRD1 gene as a negative regulator of disease progression in NSCLC. We show that DRD1, as well as other DA pathway components, are expressed in normal human lung tissue, and that loss of DRD1 expression through promoter hypermethylation is a common feature in NSCLC patients and is associated with worse survival. At the cellular level, DRD1 affects proliferation by inhibiting the activation of EGFR and mitogen-activated protein kinase 1/2 (ERK1/2). Interestingly, we also found that DRD1 regulates the expression of PD-L1 in lung cancer cells. Taken together, these results suggest that DRD1 methylation may constitute a biomarker of poor prognosis in NSCLC patients while other components of this pathway could be targeted to improve response to EGFR- and PD-L1-targeted therapies.

Basin scale variability of Ostreopsis spp. blooms provides evidence of effectiveness of an integrated sampling approach.

Ostreopsis spp. blooms have been occurring in the last two decades in the Mediterranean Sea in association with a variety of biotic and abiotic substrata (macroalgae, seagrasses, benthic invertebrates, sand, pebbles and rocks). Cells proliferate attached to the surfaces through mucilaginous trichocysts, which lump together microalgal cells, and can also be found in the plankton and on floating aggregates: such tychoplanktonic behavior makes the quantitative assessment of blooms more difficult than planktonic or benthic ones. Different techniques have been so far applied for quantifying cell abundances of benthic microalgae for research, monitoring and risk assessment purposes. In this context, the Benthic Dinoflagellates Integrator (BEDI), a non-destructive quantification method for benthic dinoflagellate abundances, was developed and tested within the EU ENPI-CBCMED project M3-HABs. This device allows mechanical detachment of cells without collecting the benthic substrate, providing an integrated assessment of both epiphytic and planktonic cells, i.e. of the number of cells potentially made available in the water volume from "resuspension" which could have harmful effects on other organisms (including humans). The present study confirms the effectiveness of the BEDI sampling device across different environments across the Mediterranean Sea and constitutes the first large-scale study of Ostreopsis spp. blooms magnitude in function of different macro- and meso­habitat features across the basin.

Marine phycotoxin levels in shellfish-14 years of data gathered along the Italian coast.

Along the Italian coasts, toxins of algal origin in wild and cultivated shellfish have been reported since the 1970s. In this study, we used data gathered by the Veterinary Public Health Institutes (IZS) and the Italian Environmental Health Protection Agencies (ARPA) from 2006 to 2019 to investigate toxicity events along the Italian coasts and relate them to the distribution of potentially toxic species. Among the detected toxins (OA and analogs, YTXs, PTXs, STXs, DAs, AZAs), OA and YTX were those most frequently reported. Levels exceeding regulatory limits in the case of OA (≤2,448 µg equivalent kg-1) were associated with high abundances of Dinophysis spp., and in the case of YTXs (≤22 mg equivalent kg-1) with blooms of Gonyaulax spinifera, Lingulodinium polyedra, and Protoceratium reticulatum. Seasonal blooms of Pseudo-nitzschia spp. occur all along the Italian coast, but DA has only occasionally been detected in shellfish at concentrations always below the regulatory limit (≤18 mg kg-1). Alexandrium spp. were recorded in several areas, although STXs (≤13,782 µg equivalent kg-1) rarely and only in few sites exceeded the regulatory limit in shellfish. Azadinium spp. have been sporadically recorded, and AZAs have been sometimes detected but always in low concentrations (≤7 µg equivalent kg-1). Among the emerging toxins, PLTX-like toxins (≤971 µg kg-1 OVTX-a) have often been detected mainly in wild mussels and sea urchins from rocky shores due to the presence of Ostreopsis cf. ovata. Overall, Italian coastal waters harbour a high number of potentially toxic species, with a few HAB hotspots mainly related to DSP toxins. Nevertheless, rare cases of intoxications have occurred so far, reflecting the whole Mediterranean Sea conditions.

A mechanistic rationale for MEK inhibitor therapy in myeloma based on blockade of MAF oncogene expression.

Modulating aberrant transcription of oncogenes is a relatively unexplored opportunity in cancer therapeutics. In approximately 10% of multiple myelomas, the initiating oncogenic event is translocation of musculoaponeurotic fibrosarcoma oncogene homolog (MAF), a transcriptional activator of key target genes, including cyclinD2. Our prior work showed that MAF is up-regulated in an additional 30% of multiple myeloma cases. The present study describes a common mechanism inducing MAF transcription in both instances. The second mode of MAF transcription occurred in myelomas with multiple myeloma SET domain (MMSET) translocation. MMSET knockdown decreased MAF transcription and cell viability. A small-molecule screen found an inhibitor of mitogen-activated protein kinase kinase (MEK), which activates extracellular signal-regulated kinase (ERK)-MAP kinases, reduced MAF mRNA in cells representing MMSET or MAF subgroups. ERK activates transcription of FOS, part of the AP-1 transcription factor. By chromatin immunoprecipitation, FOS bound the MAF promoter, and MEK inhibition decreased this interaction. MEK inhibition selectively induced apoptosis in MAF-expressing myelomas, and FOS inactivation was similarly toxic. Reexpression of MAF rescued cells from death induced by MMSET depletion, MEK inhibition, or FOS inactivation. The data presented herein demonstrate that the MEK-ERK pathway regulates MAF transcription, providing molecular rationale for clinical evaluation of MEK inhibitors in MAF-expressing myeloma.

Immunoparesis and monoclonal gammopathy of undetermined significance are disassociated in advanced age.

Immunoparesis and a skewed serum free light chain (FLC) ratio are indicators of immune dysfunction predictive of progression from monoclonal gammopathy of undetermined significance (MGUS) to multiple myeloma (MM). Previous studies have reported increased prevalence of MGUS by age, but no study has examined the relationship between immunoparesis and abnormal FLC ratios in the elderly. We screened 453 older adults (median age, 80 years; range, 65-96) to characterize the patterns of immunoparesis and abnormal FLC ratio in relation to MGUS. We defined MGUS in 4.4% of the subjects; the prevalence was 12.5% among individuals of >90 years. In MGUS (vs. non-MGUS) cases, immunoparesis and abnormal FLC ratios were detected in 70.0% (vs. 49.0%; P = 0.07) and 50.0% (vs. 12.9%; P = 0.0001), respectively. Based on small numbers, MGUS patients with abnormal FLC ratio were borderline (P = 0.07) more likely to have immunoparesis. Overall, the prevalence of immunoparesis varied in a nonlinear fashion, with lowest frequencies in the youngest and oldest groups. Our observed disassociation between MGUS prevalence and impaired immunoglobulin production suggests that separate mechanisms are involved in the development of MGUS and immunoparesis in advanced age. These findings emphasize the need for molecularly defined methods to characterize myeloma precursor states and better predict progression to MM.

Diverse eukaryotic phytoplankton from around the Marquesas Islands documented by combined microscopy and molecular techniques.

Oceanic phytoplankton serve as a base for the food webs within the largest planetary ecosystem. Despite this, surprisingly little is known about species composition, function and ecology of phytoplankton communities, especially for vast areas of the open ocean. In this study we focus on the marine phytoplankton microflora from the vicinity of the Marquesas Islands in the Southern Pacific Ocean collected during the Tara Oceans expedition. Multiple samples from four sites and two depths were studied in detail using light microscopy, scanning electron microscopy, and automated confocal laser scanning microscopy. In total 289 taxa were identified, with Dinophyceae and Bacillariophyceae contributing 60% and 32% of taxa, respectively, to phytoplankton community composition. Notwithstanding, a large number of cells could not be assigned to any known species. Coccolithophores and other flagellates together contributed less than 8% to the species list. Observed cell densities were generally low, but at sites of high autotrophic biomass, diatoms reached the highest cell densities (1.26 × 104 cells L-1). Overall, 18S rRNA metabarcode-based community compositions matched microscopy-based estimates, particularly for the main diatom taxa, indicating consistency and complementarity between different methods, while the wide range of microscopy-based methods permitted several unknown and poorly studied taxa to be revealed and identified.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes.

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.