Fishing for the Microbiome of Tropical Tuna.

Although tunas represent a significant part of the global fish economy and a major nutritional resource worldwide, their microbiome still remains poorly documented. Here, we conducted an analysis of the taxonomic composition of the bacterial communities inhabiting the gut, skin, and liver of two most consumed tropical tuna species (skipjack and yellowfin), from individuals caught in the Atlantic and Indian oceans. We hypothesized that each organ harbors a specific microbial assemblage whose composition might vary according to different biotic (sex, species) and/or abiotic (environmental) factors. Our results revealed that the composition of the tuna microbiome was totally independent of fish sex, regardless of the species and ocean considered. Instead, the main determinants of observed diversity were (i) tuna species for the gut and (ii) sampling site for the skin mucus layer and (iii) a combination of both parameters for the liver. Interestingly, 4.5% of all amplicon sequence variants (ASV) were shared by the three organs, highlighting the presence of a core-microbiota whose most abundant representatives belonged to the genera Mycoplasma, Cutibacterium, and Photobacterium. Our study also revealed the presence of a unique and diversified bacterial assemblage within the tuna liver, comprising a substantial proportion of potential histamine-producing bacteria, well known for their pathogenicity and their contribution to fish poisoning cases. These results indicate that this organ is an unexplored microbial niche whose role in the health of both the host and consumers remains to be elucidated.

Resistance of the oyster pathogen Vibrio tasmaniensis LGP32 against grazing by Vannella sp. marine amoeba involves Vsm and CopA virulence factors.

Vibrios are ubiquitous in marine environments and opportunistically colonize a broad range of hosts. Strains of Vibrio tasmaniensis present in oyster farms can thrive in oysters during juvenile mortality events and behave as facultative intracellular pathogen of oyster haemocytes. Herein, we wondered whether V. tasmaniensis LGP32 resistance to phagocytosis is specific to oyster immune cells or contributes to resistance to other phagocytes, like marine amoebae. To address this question, we developed an integrative study, from the first description of amoeba diversity in oyster farms to the characterization of LGP32 interactions with amoebae. An isolate of the Vannella genus, Vannella sp. AP1411, which was collected from oyster farms, is ubiquitous, and belongs to one clade of Vannella that could be found associated with Vibrionaceae. LGP32 was shown to be resistant to grazing by Vannella sp. AP1411 and this phenotype depends on some previously identified virulence factors: secreted metalloprotease Vsm and copper efflux p-ATPase CopA, which act at different steps during amoeba-vibrio interactions, whereas some other virulence factors were not involved. Altogether, our work indicates that some virulence factors can be involved in multi-host interactions of V. tasmaniensis ranging from protozoans to metazoans, potentially favouring their opportunistic behaviour.

Exceptional but vulnerable microbial diversity in coral reef animal surface microbiomes.

Coral reefs host hundreds of thousands of animal species that are increasingly threatened by anthropogenic disturbances. These animals host microbial communities at their surface, playing crucial roles for their fitness. However, the diversity of such microbiomes is mostly described in a few coral species and still poorly defined in other invertebrates and vertebrates. Given the diversity of animal microbiomes, and the diversity of host species inhabiting coral reefs, the contribution of such microbiomes to the total microbial diversity of coral reefs could be important, yet potentially vulnerable to the loss of animal species. Analysis of the surface microbiome from 74 taxa, including teleost fishes, hard and soft corals, crustaceans, echinoderms, bivalves and sponges, revealed that more than 90% of their prokaryotic phylogenetic richness was specific and not recovered in surrounding plankton. Estimate of the total richness associated with coral reef animal surface microbiomes reached up to 2.5% of current estimates of Earth prokaryotic diversity. Therefore, coral reef animal surfaces should be recognized as a hotspot of marine microbial diversity. Loss of the most vulnerable reef animals expected under present-day scenarios of reef degradation would induce an erosion of 28% of the prokaryotic richness, with unknown consequences on coral reef ecosystem functioning.

Geological gas-storage shapes deep life.

Around the world, several dozen deep sedimentary aquifers are being used for storage of natural gas. Ad hoc studies of the microbial ecology of some of them have suggested that sulfate reducing and methanogenic microorganisms play a key role in how these aquifers' communities function. Here, we investigate the influence of gas storage on these two metabolic groups by using high-throughput sequencing and show the importance of sulfate-reducing Desulfotomaculum and a new monophyletic methanogenic group. Aquifer microbial diversity was significantly related to the geological level. The distance to the stored natural gas affects the ratio of sulfate-reducing Firmicutes to deltaproteobacteria. In only one aquifer, the methanogenic archaea dominate the sulfate-reducers. This aquifer was used to store town gas (containing at least 50% H2 ) around 50 years ago. The observed decrease of sulfates in this aquifer could be related to stimulation of subsurface sulfate-reducers. These results suggest that the composition of the microbial communities is impacted by decades old transient gas storage activity. The tremendous stability of these gas-impacted deep subsurface microbial ecosystems suggests that in situ biotic methanation projects in geological reservoirs may be sustainable over time.

Deep ocean prokaryotic communities are remarkably malleable when facing long-term starvation.

The bathypelagic ocean is one of the largest ecosystems on Earth and sustains half of the ocean's microbial activity. This microbial activity strongly relies on surface-derived particles, but there is growing evidence that the carbon released through solubilization of these particles may not be sufficient to meet the energy demands of deep ocean prokaryotes. To explore how bathypelagic prokaryotes respond to the absence of external inputs of carbon, we followed the long-term (1 year) dynamics of an enclosed community. Despite the lack of external energy supply, we observed a continuous succession of active prokaryotic phylotypes, which was driven by recruitment of taxa from the seed bank (i.e., initially rare operational taxonomic units [OTUs]). A single OTU belonging to Marine Group I of Thaumarchaeota, which was originally rare, dominated the microbial community for ∼ 4 months and played a fundamental role in this succession likely by introducing new organic carbon through chemolithoautotrophy. This carbon presumably produced a priming effect, because after the decline of Thaumarchaeota, the diversity and metabolic potential of the community increased back to the levels present at the start of the experiment. Our study demonstrates the profound versatility of deep microbial communities when facing organic carbon deprivation.

Metal contaminations impact archaeal community composition, abundance and function in remote alpine lakes.

Using the 16S rRNA and mcrA genes, we investigated the composition, abundance and activity of sediment archaeal communities within 18 high-mountain lakes under contrasted metal levels from different origins (bedrock erosion, past-mining activities and atmospheric depositions). Bathyarchaeota, Euryarchaeota and Woesearchaeota were the major phyla found at the meta-community scale, representing 48%, 18.3% and 15.2% of the archaeal community respectively. Metals were equally important as physicochemical variables in explaining the assemblage of archaeal communities and their abundance. Methanogenesis appeared as a process of central importance in the carbon cycle within sediments of alpine lakes as indicated by the absolute abundance of methanogen 16S rRNA and mcrA gene transcripts (105 to 109 copies g-1 ). We showed that methanogen abundance and activity were significantly reduced with increasing concentrations of Pb and Cd, two indicators of airborne metal contaminations. Considering the ecological importance of methanogenesis in sediment habitats, these metal contaminations may have system wide implications even in remote area such as alpine lakes. Overall, this work was pioneer in integrating the effect of long-range atmospheric depositions on archaeal communities and indicated that metal contamination might significantly compromise the contribution of Archaea to the carbon cycling of the mountain lake sediments.

Analysis of metagenomic data reveals common features of halophilic viral communities across continents.

Microbial communities from hypersaline ponds, dominated by halophilic archaea, are considered specific of such extreme conditions. The associated viral communities have accordingly been shown to display specific features, such as similar morphologies among different sites. However, little is known about the genetic diversity of these halophilic viral communities across the Earth. Here, we studied viral communities in hypersaline ponds sampled on the coast of Senegal (8-36% of salinity) using metagenomics approach, and compared them with hypersaline viromes from Australia and Spain. The specificity of hyperhalophilic viruses could first be demonstrated at a community scale, salinity being a strong discriminating factor between communities. For the major viral group detected in all samples (Caudovirales), only a limited number of halophilic Caudovirales clades were highlighted. These clades gather viruses from different continents and display consistent genetic composition, indicating that they represent related lineages with a worldwide distribution. Non-tailed hyperhalophilic viruses display a greater rate of gene transfer and recombination, with uncharacterized genes conserved across different kind of viruses and plasmids. Thus, hypersaline viral communities around the world appear to form a genetically consistent community that are likely to harbour new genes coding for enzymes specifically adapted to these environments.

Wastewater treatment plant effluents change abundance and composition of ammonia-oxidizing microorganisms in mediterranean urban stream biofilms.

Streams affected by wastewater treatment plant (WWTP) effluents are hotspots of nitrification. We analyzed the influence of WWTP inputs on the abundance, distribution, and composition of epilithic ammonia-oxidizing (AO) assemblages in five Mediterranean urban streams by qPCR and amoA gene cloning and sequencing of both archaea (AOA) and bacteria (AOB). The effluents significantly modified stream chemical parameters, and changes in longitudinal profiles of both NH(4)(+) and NO(3)(-) indicated stimulated nitrification activity. WWTP effluents were an allocthonous source of both AOA, essentially from the Nitrosotalea cluster, and mostly of AOB, mainly Nitrosomonas oligotropha, Nitrosomonas communis, and Nitrosospira spp. changing the relative abundance and the natural composition of AO assemblages. Under natural conditions, Nitrososphaera and Nitrosopumilus AOA dominated AO assemblages, and AOB were barely detected. After the WWTP perturbation, epilithic AOB increased by orders of magnitude whereas AOA did not show quantitative changes but a shift in population composition to dominance of Nitrosotalea spp. The foraneous AOB successfully settled in downstream biofilms and probably carried out most of the nitrification activity. Nitrosotalea were only observed downstream and only in biofilms exposed to either darkness or low irradiance. In addition to other potential environmental limitations for AOA distribution, this result suggests in situ photosensitivity as previously reported for Nitrosotalea under laboratory conditions.

Targeting spatiotemporal dynamics of planktonic SAGMGC-1 and segregation of ammonia-oxidizing thaumarchaeota ecotypes by newly designed primers and quantitative polymerase chain reaction.

The annual dynamics of three different ammonia-oxidizing archaea (AOA) ecotypes (amoA gene) and of the SAGMGC-1 (Nitrosotalea-like aquatic Thaumarchaeota) group (16S rRNA gene) were studied by newly designed specific primers and quantitative polymerase chain reaction analysis in a deep oligotrophic high mountain lake (Lake Redon, Limnological Observatory of the Pyrenees, Spain). We observed segregated distributions of the main AOA populations, peaking separately in time and space, and under different ammonia concentrations and irradiance conditions. Strong positive correlation in gene abundances was found along the annual survey between 16S rRNA SAGMAGC-1 and one of the amoA ecotypes suggesting the potential for ammonia oxidation in the freshwater SAGMAGC-1 clade. We also observed dominance of Nitrosotalea-like ecotypes over Nitrosopumilus-like (Marine Group 1.1a) and not the same annual dynamics for the two thaumarchaeotal clades. The fine scale segregation in space and time of the different AOA ecotypes indicated the presence of phylogenetically close but ecologically segregated AOA species specifically adapted to specific environmental conditions. It remains to be elucidated what would be such environmental drivers.

Trace elements in the water column of high-altitude Pyrenean lakes: Impact of local weathering and long-range atmospheric input.

High altitude (alpine) lakes are efficient sentinels of environmental processes, including local pollution and long-range atmospheric transfer, because these lakes are highly vulnerable to ongoing climate changes and increasing anthropogenic pressure. Towards improving the knowledge of trace element geochemistry in the water column of alpine lakes, we assessed 64 physico-chemical parameters, including macro- and micronutrients, major and trace element concentrations in the water column of 18 lakes in the Pyrenees, located along the border between France and Spain. Lake depth, morphology, retention time and watershed rock lithology did not exhibit sizable impact on major and trace element concentrations in the water column. However, acidic (pH = 4.7 ± 0.2) lakes were distinctly different from circumneutral lakes (pH = 6.8 ± 0.5) as they exhibited >10 times higher concentrations of SO42- and trace metals (Fe, Mn, Zn, Cd, Pb, Co, Ni, Be, Al, Ga and REEs). While some of these elements clearly mark the presence of sulphide-rich minerals within the watershed (Fe, Zn, Cd and Pb), the increased mobility of lithogenic elements (Be, Al, Ga and REEs) in acidic lakes may reflect the leaching of these elements from silicate dust derived from atmospheric deposits or surrounding granites. At the same time, compared to circumneutral lakes, acidic lake water displayed lower concentrations of dissolved oxyanions (As, Mo, V, B and W) and elevated SO42- concentrations. The latter could lead to efficient Ba removal from the water column. The exploitation of metal ores within the watershed of three lakes clearly impacted high Zn and Cd concentrations observed in their water column, despite two of these lakes not being acidic. We conclude that local impacts have a greater effect on the water column than long-range atmospheric inputs and that dissolved trace element concentration measurements can be used for revealing sulphide-rich minerals or acid mine drainage within the lakes' watershed.

Refocusing the microbial rare biosphere concept through a functional lens.

The influential concept of the rare biosphere in microbial ecology has underscored the importance of taxa occurring at low abundances yet potentially playing key roles in communities and ecosystems. Here, we refocus the concept of rare biosphere through a functional trait-based lens and provide a framework to characterize microbial functional rarity, a combination of numerical scarcity across space or time and trait distinctiveness. We demonstrate how this novel interpretation of the rare biosphere, rooted in microbial functions, can enhance our mechanistic understanding of microbial community structure. It also sheds light on functionally distinct microbes, directing conservation efforts towards taxa harboring rare yet ecologically crucial functions.

Metagenomic data of bacterial communities associated with Acropora species from Phu Quoc Islands, Vietnam.

Acropora is one of the most common coral genera found in Phu Quoc Islands, Vietnam. However, the presence of marine snails, such as the coralllivorous gastropod Drupella rugosa, was a potential threat to the survival of many scleractinian species, leading to changes in the health status and bacterial diversity of coral reefs in Phu Quoc Islands. Here, we describe the composition of bacterial communities associated with two species of Acropora (Acropora formosa and Acropora millepora) using the Illumina sequencing technology. This dataset includes 5 coral samples of each status (grazed or healthy), which were collected in Phu Quoc Islands (9°55'20.6″N 104°01'16.4″E) in May 2020. A total of 19 phyla, 34 classes, 98 orders, 216 families and 364 bacterial genera were detected from 10 coral samples. Overall, Proteobacteria and Firmicutes were the two most common bacterial phyla in all samples. Significant differences in the relative abundances of the genera Fusibacter, Halarcobacter, Malaciobacter, and Thalassotalea between grazed and healthy status were observed. However, there was no differences in alpha diversity indices between the two status. Furthermore, the dataset analysis also indicated that Vibrio and Fusibacter were core genera in the grazed samples, whereas Pseudomonas was the core genus in the healthy samples.

Identifying macroplastic pathobiomes and antibiotic resistance in a subtropical fish farm.

Macroplastics are ubiquitous in aquaculture ecosystems. However, to date the potential role of plastics as a support for bacterial biofilm that can include potential human pathogenic bacteria (PHPB) and antibiotic-resistant bacteria (ARB) has been largely overlooked. In this study, we used a combination of metabarcoding and standard antibiotic susceptibility testing to study the pathobiome and resistome of macroplastics, fish guts and the environment in a marine aquaculture farm in Mauritius. Aquaculture macroplastics were found to be higher in PHPB, dominated by the Vibrionaceae family (0.34 % of the total community), compared with environmental samples. Moreover, isolates from aquaculture plastics showed higher significant multiple antibiotic resistance (MAR) compared to non-plastic samples of seawater, sediment and fish guts. These results suggest that plastics act as a reservoir and fomite of PHPB and ARB in aquaculture, potentially threatening the health of farmed fish and human consumers.

Metals and metalloids in high-altitude Pyrenean lakes: sources and distribution in pre-industrial and modern sediments.

High-altitude Pyrenean lakes are ecosystems far from local pollution sources, and thus they are particularly sensitive to the atmospheric deposition of metals and metalloids. This study aims to quantify the effect of human activity in 18 lakes located in both side of the France-Spain frontier. Sediment cores were collected in summer 2013, sampled at a 1cm resolution and the concentration of 24 elements was measured by ICP-MS. Statistic and chemometric analysis of the results highlights the influence of the geographical position and lithogenic features of each lake basin on trapping pollutants. More than the 80% of the lakes showed values of enrichment factor (EF) above 2 for at least one of the elements investigated in at least one core interval, which corroborates the existence of historical anthropogenic inputs of elements in the studied area. The results demonstrate the natural origin of As and Ti in Pyrenees, together with the significant anthropogenic inputs of Cd, Pb, Sb and Sn from ancient times. The data set points mining activities as the main historical source of pollution and illustrate the large impact of the industrial revolution. The regional variability could reflect also differential long-range transport, followed by dry or wet deposition.

Distribution and genomic variation of ammonia-oxidizing archaea in abyssal and hadal surface sediments.

Ammonia-oxidizing archaea of the phylum Thaumarchaeota play a central role in the biogeochemical cycling of nitrogen in benthic sediments, at the interface between pelagic and subsurface ecosystems. However, our understanding of their niche separation and of the processes controlling their population structure in hadal and abyssal surface sediments is still limited. Here, we reconstructed 47 AOA metagenome-assembled genomes (MAGs) from surface sediments of the Atacama and Kermadec trench systems. They formed deep-sea-specific groups within the family Nitrosopumilaceae and were assigned to six amoA gene-based clades. MAGs from different clades had distinct distribution patterns along oxygen-ammonium counter gradients in surface sediments. At the species level, MAGs thus seemed to form different ecotypes and follow deterministic niche-based distributions. In contrast, intraspecific population structure, defined by patterns of Single Nucleotide Variants (SNV), seemed to reflect more complex contributions of both deterministic and stochastic processes. Firstly, the bathymetric range had a strong effect on population structure, with distinct populations in abyssal plains and hadal trenches. Then, hadal populations were clearly separated by trench system, suggesting a strong isolation-by-topography effect, whereas abyssal populations were rather controlled by sediment depth or geographic distances, depending on the clade considered. Interestingly, genetic variability between samples was lowest in sediment layers where the mean MAG coverage was highest, highlighting the importance of selective pressure linked with each AOA clade's ecological niche. Overall, our results show that deep-sea AOA genome distributions seem to follow both deterministic and stochastic processes, depending on the genomic variability scale considered.

Shift and homogenization of gut microbiome during invasion in marine fishes.

Biological invasion is one of the main components of global changes in aquatic ecosystems. Unraveling how establishment in novel environments affects key biological features of animals is a key step towards understanding invasion. Gut microbiome of herbivorous animals is important for host health but has been scarcely assessed in invasive species. Here, we characterized the gut microbiome of two invasive marine herbivorous fishes (Siganus rivulatus and Siganus luridus) in their native (Red Sea) and invaded (Mediterranean Sea) ranges. The taxonomic and phylogenetic diversity of the microbiome increased as the fishes move away from the native range and its structure became increasingly different from the native microbiome. These shifts resulted in homogenization of the microbiome in the invaded range, within and between the two species. The shift in microbial diversity was associated with changes in its functions related with the metabolism of short-chain fatty acids. Altogether, our results suggest that the environmental conditions encountered by Siganidae during their expansion in Mediterranean ecosystems strongly modifies the composition of their gut microbiome along with its putative functions. Further studies should pursue to identify the precise determinants of these modifications (e.g. changes in host diet or behavior, genetic differentiation) and whether they participate in the ecological success of these species.

Tracking spoilage bacteria in the tuna microbiome.

Like other seafood products, tuna is highly perishable and sensitive to microbial spoilage. Its consumption, whether fresh or canned, can lead to severe food poisoning due to the activity of specific microorganisms, including histamine-producing bacteria. Yet, many grey areas persist regarding their ecology, conditions of emergence, and proliferation in fish. In this study, we used 16S rRNA barcoding to investigate postmortem changes in the bacteriome of fresh and brine-frozen yellowfin tuna (Thunnus albacares), until late stages of decomposition (i.e. 120 h). The results revealed that despite standard refrigeration storage conditions (i.e. 4°C), a diverse and complex spoilage bacteriome developed in the gut and liver. The relative abundance of spoilage bacterial taxa increased rapidly in both organs, representing 82% of the bacterial communities in fresh yellowfin tuna, and less than 30% in brine-frozen tuna. Photobacterium was identified as one of the dominant bacterial genera, and its temporal dynamics were positively correlated with histamine concentration in both gut and liver samples, which ultimately exceeded the recommended sanitary threshold of 50 ppm in edible parts of tuna. The results from this study show that the sanitary risks associated with the consumption of this widely eaten fish are strongly influenced by postcapture storage conditions.

Seasonal changes of freshwater ammonia-oxidizing archaeal assemblages and nitrogen species in oligotrophic alpine lakes.

The annual changes in the composition and abundance of ammonia-oxidizing archaea (AOA) were analyzed monthly in surface waters of three high mountain lakes within the Limnological Observatory of the Pyrenees (LOOP; northeast Spain) using both 16S rRNA and functional (ammonia monooxygenase gene, amoA) gene sequencing as well as quantitative PCR amplification. The set of biological data was related to changes in nitrogen species and to other relevant environmental variables. The whole archaeal assemblage was dominated by phylotypes closely related to the crenarchaeal 1.1a group (58% ± 18% of total 16S rRNA gene sequences), and consistent structural changes were detected during the study. Water temperature was the environmental variable that better explained spring, summer, and winter (ice-covered lakes) archaeal assemblage structure. The amoA gene was detected year round, and seasonal changes in amoA gene composition were well correlated with changes in the archaeal 16S rRNA gene pool. In addition, copy numbers of both the specific 1.1a group 16 rRNA and archaeal amoA genes were well correlated, suggesting that most freshwater 1.1a Crenarchaeota had the potential to carry out ammonia oxidation. Seasonal changes in the diversity and abundance of AOA (i.e., amoA) were better explained by temporal changes in ammonium, the substrate for nitrification, and mostly nitrite, the product of ammonia oxidation. Lacustrine amoA gene sequences grouped in coherent freshwater phylogenetic clusters, suggesting that freshwater habitats harbor typical amoA-containing ecotypes, which is different from soils and seas. We observed within the freshwater amoA gene sequence pool a high genetic divergence (translating to up to 32% amino acid divergence) between the spring and the remaining AOA assemblages. This suggests that different AOA ecotypes are adapted to different temporal ecological niches in these lakes.

Phylogenetic ecology of widespread uncultured clades of the Kingdom Euryarchaeota.

Despite its widespread distribution and high levels of phylogenetic diversity, microbes are poorly understood creatures. We applied a phylogenetic ecology approach in the Kingdom Euryarchaeota (Archaea) to gain insight into the environmental distribution and evolutionary history of one of the most ubiquitous and largely unknown microbial groups. We compiled 16S rRNA gene sequences from our own sequence libraries and public genetic databases for two of the most widespread mesophilic Euryarchaeota clades, Lake Dagow Sediment (LDS) and Rice Cluster-V (RC-V). The inferred population history indicated that both groups have undergone specific nonrandom evolution within environments, with several noteworthy habitat transition events. Remarkably, the LDS and RC-V groups had enormous levels of genetic diversity when compared with other microbial groups, and proliferation of sequences within each single clade was accompanied by significant ecological differentiation. Additionally, the freshwater Euryarchaeota counterparts unexpectedly showed high phylogenetic diversity, possibly promoted by their environmental adaptability and the heterogeneous nature of freshwater ecosystems. The temporal phylogenetic diversification pattern of these freshwater Euryarchaeota was concentrated both in early times and recently, similarly to other much less diverse but deeply sampled archaeal groups, further stressing that their genetic diversity is a function of environment plasticity. For the vast majority of living beings on Earth (i.e. the uncultured microorganisms), how they differ in the genetic or physiological traits used to exploit the environmental resources is largely unknown. Inferring population history from 16S rRNA gene-based molecular phylogenies under an ecological perspective may shed light on the intriguing relationships between lineage, environment, evolution and diversity in the microbial world.

Seasonality of archaeal proteorhodopsin and associated Marine Group IIb ecotypes (Ca. Poseidoniales) in the North Western Mediterranean Sea.

The Archaea Marine Group II (MGII) is widespread in the world's ocean where it plays an important role in the carbon cycle. Despite recent discoveries on the group's metabolisms, the ecology of this newly proposed order (Candidatus Poseidoniales) remains poorly understood. Here we used a combination of time-series metagenome-assembled genomes (MAGs) and high-frequency 16S rRNA data from the NW Mediterranean Sea to test if the taxonomic diversity within the MGIIb family (Candidatus Thalassarchaeaceae) reflects the presence of different ecotypes. The MAGs' seasonality revealed a MGIIb family composed of different subclades that have distinct lifestyles and physiologies. The vitamin metabolisms were notably different between ecotypes with, in some, a possible link to sunlight's energy. Diverse archaeal proteorhodopsin variants, with unusual signature in key amino acid residues, had distinct seasonal patterns corresponding to changing day length. In addition, we show that in summer, archaea, as opposed to bacteria, disappeared completely from surface waters. Our results shed light on the diversity and the distribution of the euryarchaeotal proteorhodopsin, and highlight that MGIIb is a diverse ecological group. The work shows that time-series based studies of the taxonomy, seasonality, and metabolisms of marine prokaryotes is critical to uncover their diverse role in the ocean.