Ecology of aerobic anoxygenic phototrophs on a fine-scale taxonomic resolution in Adriatic Sea unravelled by unsupervised neural network.

BACKGROUND: Aerobic anoxygenic phototrophs are metabolically highly active, diverse and widespread polyphyletic members of bacterioplankton whose photoheterotrophic capabilities shifted the paradigm about simplicity of the microbial food chain. Despite their considerable contribution to the transformation of organic matter in marine environments, relatively little is still known about their community structure and ecology at fine-scale taxonomic resolution. Up to date, there is no comprehensive (i.e. qualitative and quantitative) analysis of their community composition in the Adriatic Sea. RESULTS: Analysis was based on pufM gene metabarcoding and quantitative FISH-IR approach with the use of artificial neural network. Significant seasonality was observed with regards to absolute abundances (maximum average abundances in spring 2.136 ± 0.081 × 104 cells mL-1, minimum in summer 0.86 × 104 cells mL-1), FISH-IR groups (Roseobacter clade prevalent in autumn, other Alpha- and Gammaproteobacteria in summer) and pufM sequencing data agglomerated at genus-level. FISH-IR results revealed heterogeneity with the highest average relative contribution of AAPs assigned to Roseobacter clade (37.66%), followed by Gammaproteobacteria (35.25%) and general Alphaproteobacteria (31.15%). Community composition obtained via pufM sequencing was dominated by Gammaproteobacteria clade NOR5/OM60, specifically genus Luminiphilus, with numerous rare genera present in relative abundances below 1%. The use of artificial neural network connected this community to biotic (heterotrophic bacteria, HNA and LNA bacteria, Synechococcus, Prochlorococcus, picoeukaryotes, heterotrophic nanoflagellates, bacterial production) and abiotic environmental factors (temperature, salinity, chlorophyll a and nitrate, nitrite, ammonia, total nitrogen, silicate, and orthophosphate concentration). A type of neural network, neural gas analysis at order-, genus- and ASV-level, resulted in five distinct best matching units (representing particular environments) and revealed that high diversity was generally independent of temperature, salinity, and trophic status of the environment, indicating a potentially dissimilar behaviour of aerobic anoxygenic phototrophs compared to the general bacterioplankton. CONCLUSION: This research represents the first comprehensive analysis of aerobic anoxygenic phototrophs in the Adriatic Sea on a trophic gradient during a year-round period. This study is also one of the first reports of their genus-level ecology linked to biotic and abiotic environmental factors revealed by unsupervised neural network algorithm, paving the way for further research of substantial contribution of this important bacterial functional group to marine ecosystems.

Picoplankton diversity in an oligotrophic and high salinity environment in the central Adriatic Sea.

By combining qualitative 16S metabarcoding and quantitative CARD-FISH methods with neural gas analysis, different patterns of the picoplankton community were revealed at finer taxonomic levels in response to changing environmental conditions in the Adriatic Sea. We present the results of a one-year study carried out in an oligotrophic environment where increased salinity was recently observed. We have shown that the initial state of community structure changes according to environmental conditions and is expressed as qualitative and quantitative changes. A general pattern of increasing diversity under harsh environmental conditions, particularly under the influence of increasing salinity at the expense of community abundance was observed. Considering the trend of changing seawater characteristics due to climate change, this study helps in understanding a possible structural change in the microbial community of the Adriatic Sea that could affect higher levels of the marine food web.

Evaluation of DNA extraction methods and direct PCR in metabarcoding of mock and marine bacterial communities.

Recent advances in new molecular biology methods and next-generation sequencing (NGS) technologies have revolutionized metabarcoding studies investigating complex microbial communities from various environments. The inevitable first step in sample preparation is DNA extraction which introduces its own set of biases and considerations. In this study, we assessed the influence of five DNA extraction methods [B1: phenol/chloroform/isoamyl extraction, B2 and B3: isopropanol and ethanol precipitations, respectively-both modifications of B1, K1: DNeasy PowerWater Kit (QIAGEN), K2: modified DNeasy PowerWater Kit (QIAGEN) and direct PCR approach (P) that completely circumvents this step on community composition and DNA yield of mock and marine sample communities from the Adriatic Sea]. B1-B3 methods generally produced higher DNA yields and more similar microbial communities, but with higher interindividual variability. Each method demonstrated significant differences in a specific community structure, where rare taxa seem to play a crucial role. There was not one superior method closest to the theoretically expected mock community composition, they all demonstrated skewed ratios, but in a similar way which might be attributed to other factors, such as primer bias or 16S rRNA gene count for specific taxa. Direct PCR represents an interesting approach when high throughput in sample processing is required. We emphasize the importance of making a cautious decision about the choice of the extraction method or direct PCR approach, but even more importantly its consistent application throughout the study.

How do open coastal fire episodes' impact sea surface microlayer neuston communities?

Human-induced climate change is expected to increase the frequency and severity of vegetation fires. The Mediterranean region is considered particularly prone to fire episodes in summer. It is well known that pyrogenic particles are an important source of external nutrients for the marine environment, especially in oligotrophic areas. In this study, the plankton components of the sea surface layers were integrated to evaluate, for the first time, their dynamics over six months and their response to fire events in a typical coastal area of the Adriatic Sea. Concentrations of nutrients and organic compounds, together with plankton communities were significantly higher in the sea surface microlayer (SML, < 1 mm thick), than in the underlying water from 1 m depth. The piconeuston community and chlorophyll a responded with extreme abundance and concentration to the most intense fire event that enriched the SML with NH4+. Phytoneuston abundance increased with a delay of 2 weeks, while diversity indices decreased slightly after the fire events. The large abundances of the studied piconeuston parameters could be explained by the high availability of organic compounds and the immediate availability of NH4+, while the phytoneuston community responded to an increased NO3- concentration, triggered by the fire events. We confirmed that fast-acting marine heterotrophs are important members of biogeochemical cycles associated with fire events and that, together with phytoplankton, they are unavoidable parameters to detect environmental changes.

Variabilities of biochemical properties of the sea surface microlayer: Insights to the atmospheric deposition impacts.

Atmospheric deposition (AD) of nutrients and its impact on the sea surface requires consideration of interfacial processes within the sea surface microlayer (SML), the ocean-atmosphere boundary layer of major importance for many global biogeochemical and climate-related processes. This study comprised a comprehensive dataset, including dissolved NO3-, NH4+ and PO43- in ambient aerosol particles, wet deposition and sea surface samples collected from February to July 2019 at a central Adriatic coastal site. The aerosol mean concentration of dissolved nitrogen (DIN = NO3- + NH4+) and PO43- were 48.8 ± 82.8 µmol m-3 and 0.8 ± 0.6 µmol m-3, respectively, while their total fluxes (dry + wet) ranged from 24.2 to 212.3 µmol m-2 d-1 (mean 123.2 ± 53.2 µmol m-2 d-1) and from 1.2 to 2.1 µmol m-2 d-1 (mean 1.5 ± 0.3 µmol m-2 d-1), respectively. Intensive local episodes of open biomass burning (BB) significantly increased aerosol DIN concentrations as well as DIN deposition fluxes, particularly altering the molar DIN/PO43- ratio of atmospheric samples. The DIN temporal patterns showed high variability in the SML (range 0.2-24.6 µmol L-1, mean 5.0 ± 7.1 µmol L-1) in contrast to the underlying water samples (range 0.5-4.2 µmol L-1, mean 1.9 ± 1.2 µmol L-1), with significant increases during BB periods. Variability in abundance of heterotrophic bacteria and autotrophs in the SML along with concentrations of bulk dissolved and particulate organic carbon as well as dissolved and particulate lipids and carbohydrates, gel particles and surfactants followed DIN enhancements with a two-week delay. This study showed that AD can affect the short-term scale enrichments of organic matter in the SML, especially when accompanied by BB emissions typical of the overall Mediterranean coastal environment. This could have strong implications for global air-sea exchange processes, including those of climate relevant gases, mediated by the SML.

Lineage-Specific Growth Curves Document Large Differences in Response of Individual Groups of Marine Bacteria to the Top-Down and Bottom-Up Controls.

Marine bacterioplankton represent a diverse assembly of species differing largely in their abundance, physiology, metabolic activity, and role in microbial food webs. To analyze their sensitivity to bottom-up and top-down controls, we performed a manipulation experiment where grazers were removed, with or without the addition of phosphate. Using amplicon-reads normalization by internal standard (ARNIS), we reconstructed growth curves for almost 300 individual phylotypes. Grazer removal caused a rapid growth of most bacterial groups, which grew at rates of 0.6 to 3.5 day-1, with the highest rates (>4 day-1) recorded among Rhodobacteraceae, Oceanospirillales, Alteromonadaceae, and Arcobacteraceae. Based on their growth response, the phylotypes were divided into three basic groups. Most of the phylotypes responded positively to both grazer removal as well as phosphate addition. The second group (containing, e.g., Rhodobacterales and Rhizobiales) responded to the grazer removal but not to the phosphate addition. Finally, some clades, such as SAR11 and Flavobacteriaceae, responded only to phosphate amendment but not to grazer removal. Our results show large differences in bacterial responses to experimental manipulations at the phylotype level and document different life strategies of marine bacterioplankton. In addition, growth curves of 130 phylogroups of aerobic anoxygenic phototrophs were reconstructed based on changes of the functional pufM gene. The use of functional genes together with rRNA genes may significantly expand the scientific potential of the ARNIS technique. IMPORTANCE Growth is one of the main manifestations of life. It is assumed generally that bacterial growth is constrained mostly by nutrient availability (bottom-up control) and grazing (top-down control). Since marine bacteria represent a very diverse assembly of species with different metabolic properties, their growth characteristics also largely differ accordingly. Currently, the growth of marine microorganisms is typically evaluated using microscopy in combination with fluorescence in situ hybridization (FISH). However, these laborious techniques are limited in their throughput and taxonomical resolution. Therefore, we combined a classical manipulation experiment with next-generation sequencing to resolve the growth dynamics of almost 300 bacterial phylogroups in the coastal Adriatic Sea. The analysis documented that most of the phylogroups responded positively to both grazer removal and phosphate addition. We observed significant differences in growth kinetics among closely related species, which could not be distinguished by the classical FISH technique.

Artificial neural network analysis of microbial diversity in the central and southern Adriatic Sea.

Bacteria are an active and diverse component of pelagic communities. The identification of main factors governing microbial diversity and spatial distribution requires advanced mathematical analyses. Here, the bacterial community composition was analysed, along with a depth profile, in the open Adriatic Sea using amplicon sequencing of bacterial 16S rRNA and the Neural gas algorithm. The performed analysis classified the sample into four best matching units representing heterogenic patterns of the bacterial community composition. The observed parameters were more differentiated by depth than by area, with temperature and identified salinity as important environmental variables. The highest diversity was observed at the deep chlorophyll maximum, while bacterial abundance and production peaked in the upper layers. The most of the identified genera belonged to Proteobacteria, with uncultured AEGEAN-169 and SAR116 lineages being dominant Alphaproteobacteria, and OM60 (NOR5) and SAR86 being dominant Gammaproteobacteria. Marine Synechococcus and Cyanobium-related species were predominant in the shallow layer, while Prochlorococcus MIT 9313 formed a higher portion below 50 m depth. Bacteroidota were represented mostly by uncultured lineages (NS4, NS5 and NS9 marine lineages). In contrast, Actinobacteriota were dominated by a candidatus genus Ca. Actinomarina. A large contribution of Nitrospinae was evident at the deepest investigated layer. Our results document that neural network analysis of environmental data may provide a novel insight into factors affecting picoplankton in the open sea environment.

Competitive feeding interactions between native Ostrea edulis and non-native Crassostrea gigas with implications of introducing C. gigas into commercial aquaculture in the eastern Adriatic Sea.

In order to detect the possible regulatory effect of non-native C. gigas on the native O. edulis, under aquaculture conditions, feeding interactions between them were investigated in a highly productive environment of Lim Bay (Adriatic Sea). The present study uses a multi-methodological approach, including stomach content, DNA barcoding and stable isotope analysis to elucidate the feeding ecology of two oyster species. The research confirmed a high overlap throughout the year in the feeding traits among native and non-native oyster species. Competition for food was not the only relationship that exists between the investigated species as the presence of O. edulis larvae in C. gigas stomach content was confirmed by DNA analysis. Findings are not in favour of introducing C. gigas to commercial aquaculture in any new areas in the Adriatic Sea and support the need to improve the existing O. edulis aquaculture and conserve its wild stocks.

Environmental determinants of the distribution of planktonic diplonemids and kinetoplastids in the oceans.

We analysed a widely used barcode, the V9 region of the 18S rRNA gene, to study the effect of environmental conditions on the distribution of two related heterotrophic protistan lineages in marine plankton, kinetoplastids and diplonemids. We relied on a major published dataset (Tara Oceans) where samples from the mesopelagic zone were available from just 32 of 123 locations, and both groups are most abundant in this zone. To close sampling gaps and obtain more information from the deeper ocean, we collected 57 new samples targeting especially the mesopelagic zone. We sampled in three geographic regions: the Arctic, two depth transects in the Adriatic Sea, and the anoxic Cariaco Basin. In agreement with previous studies, both protist groups are most abundant and diverse in the mesopelagic zone. In addition to that, we found that their abundance, richness, and community structure also depend on geography, oxygen concentration, salinity, temperature, and other environmental variables reflecting the abundance of algae and nutrients. Both groups studied here demonstrated similar patterns, although some differences were also observed. Kinetoplastids and diplonemids prefer tropical regions and nutrient-rich conditions and avoid high oxygen concentration, high salinity, and high density of algae.

Light and Primary Production Shape Bacterial Activity and Community Composition of Aerobic Anoxygenic Phototrophic Bacteria in a Microcosm Experiment.

Phytoplankton is a key component of aquatic microbial communities, and metabolic coupling between phytoplankton and bacteria determines the fate of dissolved organic carbon (DOC). Yet, the impact of primary production on bacterial activity and community composition remains largely unknown, as, for example, in the case of aerobic anoxygenic phototrophic (AAP) bacteria that utilize both phytoplankton-derived DOC and light as energy sources. Here, we studied how reduction of primary production in a natural freshwater community affects the bacterial community composition and its activity, focusing primarily on AAP bacteria. The bacterial respiration rate was the lowest when photosynthesis was reduced by direct inhibition of photosystem II and the highest in ambient light condition with no photosynthesis inhibition, suggesting that it was limited by carbon availability. However, bacterial assimilation rates of leucine and glucose were unaffected, indicating that increased bacterial growth efficiency (e.g., due to photoheterotrophy) can help to maintain overall bacterial production when low primary production limits DOC availability. Bacterial community composition was tightly linked to light intensity, mainly due to the increased relative abundance of light-dependent AAP bacteria. This notion shows that changes in bacterial community composition are not necessarily reflected by changes in bacterial production or growth and vice versa. Moreover, we demonstrated for the first time that light can directly affect bacterial community composition, a topic which has been neglected in studies of phytoplankton-bacteria interactions.IMPORTANCE Metabolic coupling between phytoplankton and bacteria determines the fate of dissolved organic carbon in aquatic environments, and yet how changes in the rate of primary production affect the bacterial activity and community composition remains understudied. Here, we experimentally limited the rate of primary production either by lowering light intensity or by adding a photosynthesis inhibitor. The induced decrease had a greater influence on bacterial respiration than on bacterial production and growth rate, especially at an optimal light intensity. This suggests that changes in primary production drive bacterial activity, but the effect on carbon flow may be mitigated by increased bacterial growth efficiencies, especially of light-dependent AAP bacteria. Bacterial activities were independent of changes in bacterial community composition, which were driven by light availability and AAP bacteria. This direct effect of light on composition of bacterial communities has not been documented previously.

Changes in the Trophic Pathways within the Microbial Food Web in the Global Warming Scenario: An Experimental Study in the Adriatic Sea.

A recent analysis of the Mediterranean Sea surface temperature showed significant annual warming. Since small picoplankton microorganisms play an important role in all major biogeochemical cycles, fluxes and processes occurring in marine systems (the changes at the base of the food web) as a response to human-induced temperature increase, could be amplified through the trophic chains and could also significantly affect different aspects of the structure and functioning of marine ecosystems. In this study, manipulative laboratory growth/grazing experiments were performed under in situ simulated conditions to study the structural and functional changes within the microbial food web after a 3 °C increase in temperature. The results show that a rise in temperature affects the changes in: (1) the growth and grazing rates of picoplankton, (2) their growth efficiency, (3) carrying capacities, (4) sensitivity of their production and grazing mortality to temperature, (5) satisfying protistan grazer carbon demands, (6) their preference in the selection of prey, (7) predator niche breadth and their overlap, (8) apparent uptake rates of nutrients, and (9) carbon biomass flow through the microbial food web. Furthermore, temperature affects the autotrophic and heterotrophic components of picoplankton in different ways.

Temperature and phosphorus interacts in controlling the picoplankton carbon flux in the Adriatic Sea: an experimental versus field study.

Temperature and phosphorus positively interacted in controlling picoplankton biomass production and its transfer towards higher trophic levels. Two complementary approaches (experimental and field study) indicated several coherent patterns: (1) the impact of temperature on heterotrophic bacteria was high at temperatures lower than 16°C and levelled off at higher temperatures, whereas this impact on autotrophic picoplankton was linear along the entire range of the investigated temperatures; (2) the addition of phosphorus increased the values of picoplankton production and grazing, but did not change the nature of their relationships with temperature substantially; (3) the picoplankton carbon flux towards higher trophic levels was larger during the warmer months (grazing by HNF dominated during the warmer period and by ciliates during the colder period) and also strengthened in conditions without phosphorus limitation; (4) the hypothesis that the available phosphorus can be better utilized at higher temperatures was confirmed for both autotrophic and heterotrophic picoplankton; (5) the hypothesis that the rise in temperature stimulates growth only in conditions of sufficient phosphorus was confirmed only for heterotrophic bacteria. Therefore, in the global warming scenario, an increase of the picoplankton carbon flux towards higher trophic levels can be expected in the Adriatic Sea, particularly under unlimited phosphorus conditions.

Relations between mercury fractions and microbial community components in seawater under the presence and absence of probable phosphorus limitation conditions.

Microbial transformations of toxic monomethylmercury (MMHg) and dissolved gaseous mercury (DGM) at the lower levels of the marine food web are not well understood, especially in oligotrophic and phosphorus-limited seas. To examine the effects of probable phosphorus limitation (PP-limitation) on relations between mercury (Hg) fractions and microorganisms, we determined the total mercury (THg), total methylated mercury (MeHg), DGM, and microbiological and chemical parameters in the Central Adriatic Sea. Using statistical analysis, we assessed the potential microbial effects on Hg transformations and bioaccumulation. Only in the absence of PP-limitation conditions (NO-PP-limitation) is MeHg significantly related to most chemical and microbial parameters, indicating metabolism-dependent Hg transformations. The heterotrophic activity of low nucleic acid bacteria (abundant in oligotrophic regions) seems responsible for most of Hg methylation under NO-PP-limitation. Under these conditions, DGM is strongly related to microbial fractions and chlorophyll a, indicating biological DGM production, which is probably not metabolically induced, as most of these relations are also observed under PP-limitation. MMHg biomagnification was observed through an increased bioaccumulation factor from microseston to mesozooplankton. Our results indicate that Hg transformations and uptake might be enhanced under NO-PP-limitation conditions, emphasizing their impact on the transfer of Hg to higher trophic levels.

The effect of temperature increase on microbial carbon fluxes in the Adriatic Sea: an experimental approach.

An assessment of the temperature increase effect on processes within the microbial food web provides a better insight into the carbon transfer and energy flow processes in marine environments in the global warming perspective. Modified laboratory dilution experiments that allow simultaneous estimates of protozoan grazing and viral lysis on picoplankton groups (bacteria, Prochlorococcus, Synechococcus and pico-eukaryotic algae) under in situ and 3°C above in situ temperatures were performed at seasonal scale. Picoplankton mortality due to grazing was generally higher than that caused by viral lysis, especially in the cold months. The largest part of HNF carbon demand was satisfied by grazing on bacteria throughout the year. Although ciliates satisfied their carbon demand predominantly through grazing on HNF and bacteria, the role of autotrophic picoplankton (APP) as their prey increased significantly in the cold months. Bacteria constituted the most important host for viruses throughout the year. However, during the warm months, APP groups were also significant hosts for viral infection. Under the warming condition the amount of picoplankton biomass transferred to protozoan grazers exceeded the lysed biomass, suggesting that global warming could further increase picoplankton carbon flow toward higher trophic levels in the Adriatic Sea.

Distribution of aerobic anoxygenic phototrophs in the Eastern Adriatic Sea.

The spatial patterns of aerobic anoxygenic phototrophs abundances were investigated, for the first time, in the Adriatic Sea. Also, the spatial patterns of the whole picoplankton community as well as the environmental factors that potentially influence these patterns were highlighted. AAP abundances was in average 66.9 ± 66.8 × 103 cell mL-1, and their proportion in total bacteria was 7.3 ± 4.3%. These values are in the upper range of AAP abundances observed in marine environments. Multivariate analyses proved that environmental factors influenced the picoplankton community interdependently. Chl a was the main driving factor for the picoplankton community, accounting for 33.3% of picoplankton community variance, followed by NO2- (17.9% of variance explained) and temperature (14.2% of variance explained). Chl a showed stronger correlation with AAPs, non-pigmented bacteria and Picoeucaryotes than with cyanobacteria. Abundance of cyanobacteria was stronger correlated to salinity and the N:P ratio than to nutrient concentrations.

The mitochondrial T16189C polymorphism is associated with coronary artery disease in Middle European populations.

BACKGROUND: The pivotal role of mitochondria in energy production and free radical generation suggests that the mitochondrial genome could have an important influence on the expression of multifactorial age related diseases. Substitution of T to C at nucleotide position 16189 in the hypervariable D-loop of the control region (CR) of mitochondrial DNA (mtDNA) has attracted research interest because of its suspected association with various multifactorial diseases. The aim of the present study was to compare the frequency of this polymorphism in the CR of mtDNA in patients with coronary artery disease (CAD, n = 482) and type 2 diabetes mellitus (T2DM, n = 505) from two study centers, with healthy individuals (n = 1481) of Middle European descent in Austria. METHODOLOGY AND PRINCIPAL FINDINGS: CR polymorphisms and the nine major European haplogroups were identified by DNA sequencing and primer extension analysis, respectively. Frequencies and Odds Ratios for the association between cases and controls were calculated. Compared to healthy controls, the prevalence of T16189C was significantly higher in patients with CAD (11.8% vs 21.6%), as well as in patients with T2DM (11.8% vs 19.4%). The association of CAD, but not the one of T2DM, with T16189C remained highly significant after correction for age, sex and body mass index (BMI) and was independent of the two study centers. CONCLUSIONS AND SIGNIFICANCE: Our results show for the first time a significant association of T16189C with CAD in a Middle European population. As reported in other studies, in patients with T2DM an association with T16189C in individuals of European decent remains questionable.

Mitochondrial DNA haplogroup T is associated with coronary artery disease and diabetic retinopathy: a case control study.

BACKGROUND: There is strong and consistent evidence that oxidative stress is crucially involved in the development of atherosclerotic vascular disease. Overproduction of reactive oxygen species (ROS) in mitochondria is an unifying mechanism that underlies micro- and macrovascular atherosclerotic disease. Given the central role of mitochondria in energy and ROS production, mitochondrial DNA (mtDNA) is an obvious candidate for genetic susceptibility studies on atherosclerotic processes. We therefore examined the association between mtDNA haplogroups and coronary artery disease (CAD) as well as diabetic retinopathy. METHODS: This study of Middle European Caucasians included patients with angiographically documented CAD (n = 487), subjects with type 2 diabetes mellitus with (n = 149) or without (n = 78) diabetic retinopathy and control subjects without clinical manifestations of atherosclerotic disease (n = 1527). MtDNA haplotyping was performed using multiplex PCR and subsequent multiplex primer extension analysis for determination of the major European haplogroups. Haplogroup frequencies of patients were compared to those of control subjects without clinical manifestations of atherosclerotic disease. RESULTS: Haplogroup T was significantly more prevalent among patients with CAD than among control subjects (14.8% vs 8.3%; p = 0.002). In patients with type 2 diabetes, the presence of diabetic retinopathy was also significantly associated with a higher prevalence of haplogroup T (12.1% vs 5.1%; p = 0.046). CONCLUSION: Our data indicate that the mtDNA haplogroup T is associated with CAD and diabetic retinopathy in Middle European Caucasian populations.