Diverse DNA modification in marine prokaryotic and viral communities.

DNA chemical modifications, including methylation, are widespread and play important roles in prokaryotes and viruses. However, current knowledge of these modification systems is severely biased towards a limited number of culturable prokaryotes, despite the fact that a vast majority of microorganisms have not yet been cultured. Here, using single-molecule real-time sequencing, we conducted culture-independent 'metaepigenomic' analyses (an integrated analysis of metagenomics and epigenomics) of marine microbial communities. A total of 233 and 163 metagenomic-assembled genomes (MAGs) were constructed from diverse prokaryotes and viruses, respectively, and 220 modified motifs and 276 DNA methyltransferases (MTases) were identified. Most of the MTase genes were not genetically linked with the endonuclease genes predicted to be involved in defense mechanisms against extracellular DNA. The MTase-motif correspondence found in the MAGs revealed 10 novel pairs, 5 of which showed novel specificities and experimentally confirmed the catalytic specificities of the MTases. We revealed novel alternative specificities in MTases that are highly conserved in Alphaproteobacteria, which may enhance our understanding of the co-evolutionary history of the methylation systems and the genomes. Our findings highlight diverse unexplored DNA modifications that potentially affect the ecology and evolution of prokaryotes and viruses in nature.

Contamination of antibiotics and sul and tet(M) genes in veterinary wastewater, river, and coastal sea in Thailand.

Water systems in Southeast Asia accumulate antibiotics and antibiotic resistance genes (ARGs) from multiple origins, notably including human clinics and animal farms. To ascertain the fate of antibiotics and ARGs in natural water environments, we monitored the concentrations of these items in Thailand. Here, we show high concentrations of tetracyclines (72,156.9 ng/L) and lincomycin (23,968.0 ng/L) in pig farms, followed by nalidixic acid in city canals. The city canals and rivers contained diverse distributions of antibiotics and ARGs. Assessments of targeted ARGs, including sul1, sul2, sul3, and tet(M), showed that freshwater (pig farm wastewater, rivers, and canals) consistently contained these ARGs, but these genes were less abundant in seawater. Although sulfonamides were low concentrations (<170 ng/mL), sul1 and sul2 genes were abundant in freshwater (minimum 4.4 × 10-3-maximum 1.0 × 100 copies/16S), suggesting that sul genes have disseminated over a long period, despite cessation of use of this class of antibiotics. Ubiquitous distribution of sul genes in freshwater appeared to be independent of selection pressure. In contrast, water of the coastal sea in the monitored area was not contaminated by these antibiotics or ARGs. The density of Enterobacteriales was lower in seawater than in freshwater, suggesting that the number of ARG-possessing Enterobacteriales falls after entering seawater. From the pig farms, through rivers/canals, to the coastal sea, the occurrence of tetracyclines and tet(M) exhibited some correlation, although not a strong one. However, no correlations were found between concentrations of total antibiotics and ARGs, nor between sulfonamides and sul genes. This is the first comprehensive study showing Thai features of antibiotics and ARGs contaminations. The pig farm is hot spot of antibiotics and ARGs, and sul genes ubiquitously distribute in freshwater environments, which become less abundant in seawater.

Determining the Distribution of Fluorescent Organic Matter in the Indian Ocean Using in situ Fluorometry.

In order to determine the dynamics of marine fluorescent organic matter (FOM) using high-resolution spatial data, in situ fluorometers have been used in the open ocean. In this study, we measured FOM during the Global Ocean Ship-based Hydrographic Investigations Program (GO-SHIP) expedition from early December 2019 to early February 2020, using an in situ fluorometer at 148 stations along the two meridional transects (at ∼80 and ∼57°E) in the Indian Ocean, covering latitudinal ranges from ∼6°N to ∼20°S and ∼30 to ∼65°S, respectively. The FOM data obtained from the fluorometer were corrected for known temperature dependence and calibrated using FOM data measured onboard by a benchtop fluorometer. Using the relative water mass proportions estimated from water mass analyses, we determined the intrinsic values of FOM and apparent oxygen utilization (AOU) for each of the 12 water masses observed. We then estimated the basin-scale relationship between the intrinsic FOM and the AOU, as well as the turnover time for FOM in the Indian Ocean (410 ± 19 years) in combination with the microbial respiration rate in the dark ocean (>200 m). Consistent to previous estimates in the global tropical and subtropical ocean, the FOM turnover time obtained is of the same order of magnitude as the circulation age of the Indian Ocean, indicating that the FOM is refractory and is a sink for reduced carbon in the dark ocean. A decoupling of FOM and AOU from the basin-scale relationship was also observed in the abyssal waters of the northern Indian Ocean. The local variability may be explained by the effect of sinking organic matter altered by denitrification through the oxygen-deficient zone on enhanced abyssal FOM production relative to oxygen consumption.

Viral Communities in the Global Deep Ocean Conveyor Belt Assessed by Targeted Viromics.

Viruses are an abundant, diverse and dynamic component of marine and terrestrial ecosystems. In the ocean, viruses play a key role in the biogeochemical cycles and controlling microbial abundance, diversity and evolution. Recent metagenomics studies assessed the structure of the viral community in the upper ocean. However, little is known about the compositional changes in viral communities along the deep ocean conveyor belt. To assess potential changes in the viral community in the global deep-water circulation system, water samples were collected in the core of the North Atlantic Deep Water (NADW) (∼2,500 m) and Pacific Antarctic Bottom Water (∼4,000 m). Microbial and viral abundance were evaluated by flow cytometry. Subsequently, flow cytometry was used to sort virus-like particles and next generation sequencing was applied to build DNA libraries from the sorted virus populations. The viral communities were highly diverse across different oceanic regions with high dissimilarity between samples. Only 18% of the viral protein clusters were shared between the NADW and the Pacific Antarctic Bottom Water. Few viral groups, mainly associated with uncultured environmental and uncultured Mediterranean viruses were ubiquitously distributed along the global deep-water circulation system. Thus, our results point to a few groups of widely distributed abundant viruses in addition to the presence of rare and diverse types of viruses at a local scale.

Highly variable mRNA half-life time within marine bacterial taxa and functional genes.

Messenger RNA can provide valuable insights into the variability of metabolic processes of microorganisms. However, due to uncertainties that include the stability of RNA, its application for activity profiling of environmental samples is questionable. We explored different factors affecting the decay rate of transcripts of three marine bacterial isolates using qPCR and determined mRNA half-life time of specific bacterial taxa and of functional genes by metatranscriptomics of a coastal environmental prokaryotic community. The half-life time of transcripts from 11 genes from bacterial isolates ranged from 1 to 46 min. About 80% of the analysed transcripts exhibited half-live times shorter than 10 min. Significant differences were found in the half-life time between mRNA and rRNA. The half-life time of mRNA obtained from a coastal metatranscriptome ranged from 9 to 400 min. The shortest half-life times of the metatranscriptome corresponded to transcripts from the same clusters of orthologous groups (COGs) in all bacterial classes. The prevalence of short mRNA half-life time in genes related to defence mechanisms and motility indicate a tight connection of RNA decay rate to environmental stressors. The short half-life time of RNA and its high variability needs to be considered when assessing metatranscriptomes especially in environmental samples.

Publisher Correction: Marine microbial metagenomes sampled across space and time.

Due to a typesetting error, 25 rows were omitted from Table 3 in the original version of this Data Descriptor. These missing rows correspond to the following sample names.

Occurrence of sul and tet(M) genes in bacterial community in Japanese marine aquaculture environment throughout the year: Profile comparison with Taiwanese and Finnish aquaculture waters.

The use of antibiotics in aquaculture causes selection pressure for antibiotic-resistant bacteria (ARB). Antibiotic resistance genes (ARGs) may persist in ARB and the environment for long time even after stopping drug administration. Here we show monthly differences in the occurrences of genes conferring resistance to sulfonamides (i.e. sul1, sul2, sul3), and tetracyclines (tet(M)) in Japanese aquaculture seawater accompanied by records of drug administration. sul2 was found to persist throughout the year, whereas the occurrences of sul1, sul3, and tet(M) changed month-to-month. sul3 and tet(M) were detected in natural bacterial assemblages in May and July, but not in colony-forming bacteria, thus suggesting that the sul3 was harbored by the non-culturable fraction of the bacterial community. Comparison of results from Taiwanese, Japanese, and Finnish aquaculture waters reveals that the profile of sul genes and tet(M) in Taiwan resembles that in Japan, but is distinct from that in Finland. To our knowledge, this work represents the first report to use the same method to compare the dynamics of sul genes and tet(M) in aquaculture seawater in different countries.

Bacterial Stress and Mortality may be a Source of Cell-free Enzymatic Activity in the Marine Environment.

Marine microbes play a central role in driving biogeochemical cycles. Microbial extracellular enzymatic activities (EEA) are the 'gatekeeper' of the marine carbon cycle, and these enzymes may be found attached to cells or dissolved (cell-free). Recent studies indicated that the proportion of dissolved enzymatic activity is generally similar to (if not higher than) cell-attached activity. Thus, it is critical to understand the sources and sinks of cell-free EEA in the ocean. We herein empirically tested whether bacterial stress and mortality (induced by mitomycin C) are a source of the cell-free EEA of alkaline phosphatase (APase), beta-glucosidase (BGase), and leucine aminopeptidase (LAPase). We found that bacterial stress and mortality caused relative increases in the proportion of dissolved relative to total EEA of up to 10.5% for APase, 13.5% for BGase, and 7.3% for LAPase. These relative increases in dissolved EEA corresponded to absolute increases in the cell-free pool of 4.8, 7.2, and 3.8% for APase, BGase and LAPase, respectively. Collectively, our results contribute relevant information on the origin of free dissolved extracellular enzymes in marine waters, indicating that bacterial stress and mortality are a source of cell-free enzymatic activity and suggesting a potential link between microbial interactions and the degradation of organic matter via the release of cell-free enzymes.

Teasing apart the different size pools of extracellular enzymatic activity in the ocean.

Extracellular enzymatic activity (EEA) is performed by cell-associated and cell-free (i.e., "dissolved") enzymes. This cell-free fraction is operationally defined as passing through a 0.22 µm filter. The contribution of cell-free to total EEA is comparable to the cell-associated counterpart, so it is critical to understand what controls the relative importance of cell-free versus cell-associated EEA. However, attempts to tease apart the contribution of EEAs in the so-called dissolved fraction (<0.22 µm) in general, and of the nanoparticle size fraction (0.020-0.20 µm) in particular, to the total EEA pool are lacking. Here we performed experiments with Northern and Southern Hemisphere coastal waters to characterize the potential contribution of that nanoparticle fraction to the total EEA fraction of alkaline phosphatase, beta-glucosidase and leucine aminopeptidase. We found a significant contribution (in both hemispheres) of the nanoparticle fraction to the total EEA pool (up to 53%) that differed depending on the enzyme type and location. Collectively, our results indicate that a significant fraction of the so-called "dissolved EEA" is not really dissolved but associated to nanoparticles, colloidal nanogels and/or viruses. Thus, the total marine EEA pool can actually be divided into a cell-associated, undissolved-cell-free (associated to nano-particle of different origins such as viruses and nanogels) and a dissolved-cell-free pools. Our results also imply that the dissolved EEA pool is more complex than thus far anticipated. Future research will be now needed to further characterize the factors controlling the relative importance of these different pools of EEA, which are key in the recycling of organic matter in the ocean.

Differential Response of Cafeteria roenbergensis to Different Bacterial and Archaeal Prey Characteristics.

In the marine environment, the abundance of Bacteria and Archaea is either controlled bottom-up via nutrient availability or top-down via grazing. Heterotrophic nanoflagellates (HNF) are mainly responsible for prokaryotic grazing losses besides viral lysis. However, the grazing specificity of HNF on specific bacterial and archaeal taxa is under debate. Bacteria and Archaea might have different nutritive values and surface properties affecting the growth rates of HNF. In this study, we offered different bacterial and archaeal strains with different morphologic and physiologic characteristics to Cafeteria roenbergensis, one of the most abundant and ubiquitous species of HNF in the ocean. Two Nitrosopumilus maritimus-related strains isolated from the northern Adriatic Sea (Nitrosopumilus adriaticus, Nitrosopumilus piranensis), two Nitrosococcus strains, and two fast growing marine Bacteria (Pseudoalteromonas sp. and Marinobacter sp.) were fed to Cafeteria cultures. Cafeteria roenbergensis exhibited high growth rates when feeding on Pseudoalteromonas sp., Marinobacter sp., and Nitrosopumilus adriaticus, while the addition of the other strains resulted in minimal growth. Taken together, our data suggest that the differences in growth of Cafeteria roenbergensis associated to grazing on different thaumarchaeal and bacterial strains are likely due to the subtle metabolic, cell size, and physiological differences between different bacterial and thaumarchaeal taxa. Moreover, Nitrosopumilus adriaticus experienced a similar grazing pressure by Cafeteria roenbergensis as compared to the other strains, suggesting that other HNF may also prey on Archaea which might have important consequences on the global biogeochemical cycles.

Single cell genomes of Prochlorococcus, Synechococcus, and sympatric microbes from diverse marine environments.

Prochlorococcus and Synechococcus are the dominant primary producers in marine ecosystems and perform a significant fraction of ocean carbon fixation. These cyanobacteria interact with a diverse microbial community that coexists with them. Comparative genomics of cultivated isolates has helped address questions regarding patterns of evolution and diversity among microbes, but the fraction that can be cultivated is miniscule compared to the diversity in the wild. To further probe the diversity of these groups and extend the utility of reference sequence databases, we report a data set of single cell genomes for 489 Prochlorococcus, 50 Synechococcus, 9 extracellular virus particles, and 190 additional microorganisms from a diverse range of bacterial, archaeal, and viral groups. Many of these uncultivated single cell genomes are derived from samples obtained on GEOTRACES cruises and at well-studied oceanographic stations, each with extensive suites of physical, chemical, and biological measurements. The genomic data reported here greatly increases the number of available Prochlorococcus genomes and will facilitate studies on evolutionary biology, microbial ecology, and biological oceanography.

Marine microbial metagenomes sampled across space and time.

Recent advances in understanding the ecology of marine systems have been greatly facilitated by the growing availability of metagenomic data, which provide information on the identity, diversity and functional potential of the microbial community in a particular place and time. Here we present a dataset comprising over 5 terabases of metagenomic data from 610 samples spanning diverse regions of the Atlantic and Pacific Oceans. One set of metagenomes, collected on GEOTRACES cruises, captures large geographic transects at multiple depths per station. The second set represents two years of time-series data, collected at roughly monthly intervals from 3 depths at two long-term ocean sampling sites, Station ALOHA and BATS. These metagenomes contain genomic information from a diverse range of bacteria, archaea, eukaryotes and viruses. The data's utility is strengthened by the availability of extensive physical, chemical, and biological measurements associated with each sample. We expect that these metagenomes will facilitate a wide range of comparative studies that seek to illuminate new aspects of marine microbial ecosystems.

Metagenomic insights into zooplankton-associated bacterial communities.

Zooplankton and microbes play a key role in the ocean's biological cycles by releasing and consuming copious amounts of particulate and dissolved organic matter. Additionally, zooplankton provide a complex microhabitat rich in organic and inorganic nutrients in which bacteria thrive. In this study, we assessed the phylogenetic composition and metabolic potential of microbial communities associated with crustacean zooplankton species collected in the North Atlantic. Using Illumina sequencing of the 16S rRNA gene, we found significant differences between the microbial communities associated with zooplankton and those inhabiting the surrounding seawater. Metagenomic analysis of the zooplankton-associated microbial community revealed a highly specialized bacterial community able to exploit zooplankton as microhabitat and thus, mediating biogeochemical processes generally underrepresented in the open ocean. The zooplankton-associated bacterial community is able to colonize the zooplankton's internal and external surfaces using a large set of adhesion mechanisms and to metabolize complex organic compounds released or exuded by the zooplankton such as chitin, taurine and other complex molecules. Moreover, the high number of genes involved in iron and phosphorus metabolisms in the zooplankton-associated microbiome suggests that this zooplankton-associated bacterial community mediates specific biogeochemical processes (through the proliferation of specific taxa) that are generally underrepresented in the ambient waters.

Abundance of Common Aerobic Anoxygenic Phototrophic Bacteria in a Coastal Aquaculture Area.

Aerobic anoxygenic phototrophic bacteria (AAnPB) rely on not only heterotrophic but also phototrophic energy gain. AAnPB are known to have high abundance in oligotrophic waters and are the major portion of the bacterial carbon stock in the environment. In a yearlong study in an aquaculture area in the Uwa Sea, Japan, AAnPB, accounted for 4.7 to 24% of the total bacteria by count. Since the cell volume of AAnPB is 2.23 ± 0.674 times larger than the mean for total bacteria, AAnPB biomass is estimated to account for 10-53% of the total bacterial assemblage. By examining pufM gene sequence, a common phylogenetic AAnPB species was found in all sampling sites through the year. The common species and other season-specific species were phylogenetically close to unculturable clones recorded in the Sargasso Sea and Pacific Ocean. The present study suggests that the common species may be a cosmopolitan species with worldwide distribution that is abundant not only in the oligotrophic open ocean but also in eutrophic aquaculture areas.

Distribution and Niche Separation of Planktonic Microbial Communities in the Water Columns from the Surface to the Hadal Waters of the Japan Trench under the Eutrophic Ocean.

The Japan Trench is located under the eutrophic Northwestern Pacific while the Mariana Trench that harbors the unique hadal planktonic biosphere is located under the oligotrophic Pacific. Water samples from the sea surface to just above the seafloor at a total of 11 stations including a trench axis station, were investigated several months after the Tohoku Earthquake in March 2011. High turbidity zones in deep waters were observed at most of the sampling stations. The small subunit (SSU) rRNA gene community structures in the hadal waters (water depths below 6000 m) at the trench axis station were distinct from those in the overlying meso-, bathy and abyssopelagic waters (water depths between 200 and 1000 m, 1000 and 4000 m, and 4000 and 6000 m, respectively), although the SSU rRNA gene sequences suggested that potential heterotrophic bacteria dominated in all of the waters. Potential niche separation of nitrifiers, including ammonia-oxidizing archaea (AOA), was revealed by quantitative PCR analyses. It seems likely that Nitrosopumilus-like AOAs respond to a high flux of electron donors and dominate in several zones of water columns including shallow and very deep waters. This study highlights the effects of suspended organic matter, as induced by seafloor deformation, on microbial communities in deep waters and confirm the occurrence of the distinctive hadal biosphere in global trench environments hypothesized in the previous study.

Large-scale distribution of microbial and viral populations in the South Atlantic Ocean.

Viruses are abundant, diverse and dynamic components of the marine environments and play a significant role in the ocean biogeochemical cycles. To assess potential variations in the relation between viruses and microbes in different geographic regions and depths, viral and microbial abundance and production were determined throughout the water column along a latitudinal transect in the South Atlantic Ocean. Path analysis was used to examine the relationships between several abiotic and biotic parameters and the different microbial and viral populations distinguished by flow cytometry. The depth-integrated contribution of microbial and viral abundance to the total microbial and viral biomass differed significantly among the different provinces. Additionally, the virus-to-microbe ratio increased with depth and decreased laterally towards the more productive regions. Our data revealed that the abundance of phytoplankton and microbes is the main controlling factor of the viral populations in the euphotic and mesopelagic layers, whereas in the bathypelagic realm, viral abundance was only weakly related to the biotic and abiotic variables. The relative contribution of the three viral populations distinguished by flow cytometry showed a clear geographical pattern throughout the water column, suggesting that these populations are composed of distinct taxa able to infect specific hosts. Overall, our data indicate the presence of distinct microbial patterns along the latitudinal transect. This variability is not limited to the euphotic layer but also detectable in the meso- and bathypelagic layers.

Extracellular proteases are released by ciliates in defined seawater microcosms.

The biodegradation of proteins in seawater requires various proteases which are commonly thought to be mainly derived from heterotrophic bacteria. We, however, found that protists showed a high protease activity and continuously produced trypsin-type enzymes. The free-living marine heterotrophic ciliate Paranophrys marina together with an associated bacterium was isolated and used for microcosm incubation with different concentrations of killed bacteria as food for 10 days. The results showed that the co-existence of the ciliate with its associated bacterium produced a significant protease activity in both cell-associated and cell-free fractions while that in the associated bacterium only microcosm was negligible. The protease profiles are different between cell-associated and cell-free fractions, and a trypsin-type enzyme hydrolyzing Boc-Val-Leu-Lys-MCA was detected throughout the period in the presence of ciliates. This suggests that ciliates release proteases into the surrounding environment which could play a role in protein digestion outside cells. It has been previously suggested that bacteria are the major transformers in seawater. We here present additional data which indicates that protists, or at least ciliates with their specific enzymes, are a potential player in organic matter degradation in water columns.

Biodiversity and multifunctionality in a microbial community: a novel theoretical approach to quantify functional redundancy.

Ecosystems have a limited buffering capacity of multiple ecosystem functions against biodiversity loss (i.e. low multifunctional redundancy). We developed a novel theoretical approach to evaluate multifunctional redundancy in a microbial community using the microbial genome database (MBGD) for comparative analysis. In order to fully implement functional information, we defined orthologue richness in a community, each of which is a functionally conservative evolutionary unit in genomes, as an index of community multifunctionality (MF). We constructed a graph of expected orthologue richness in a community (MF) as a function of species richness (SR), fit the power function to SR (i.e. MF = cSR(a)), and interpreted the higher exponent a as the lower multifunctional redundancy. Through a microcosm experiment, we confirmed that MF defined by orthologue richness could predict the actual multiple functions. We simulated random and non-random community assemblages using full genomic data of 478 prokaryotic species in the MBGD, and determined that the exponent in microbial communities ranged from 0.55 to 0.75. This exponent range provided a quantitative estimate that a 6.6-8.9% loss limit in SR occurred in a microbial community for an MF reduction no greater than 5%, suggesting a non-negligible initial loss effect of microbial diversity on MF.

Comparison between MICRO­CARD­FISH and 16S rRNA gene clone libraries to assess the active versus total bacterial community in the coastal Arctic.

We collected surface- and deep-water samples (maximum depth 300 m) during the spring­summer transition in the coastal Arctic along a transect in the Kongsfjorden (Ny-Ålesund, Spitsbergen, Norway) to determine the structure of the active versus total marine bacterioplankton community using different approaches. Catalysed reporter deposition­ fluorescence in situ hybridization combined with microautoradiography (MICRO­CARD­FISH) was used to determine the abundance and activity of different bacterial groups. The bacterial communities were dominated by members of Alphaproteobacteria followed by Bacteroidetes, whereas Gammaproteobacteria were present at low abundance but exhibited a high percentage of active cells taking up leucine. The clone libraries of 16S rRNA genes (16S rDNA) and 16S rRNA from two different depths were used to decipher the bacterial community structure. Independently of the type of clone libraries analysed (16S rDNA- or 16S rRNA-based), four major and four minor taxonomic groups were detected. The bacterioplankton community was mainly dominated at both the DNA and the RNA levels by Alphaproteobacteria followed by Gammaproteobacteria. The Rhodobacteriaceae were the most abundant members of the Alphaproteobacteria in both DNA and RNA clone libraries, followed by the SAR11 clade, which was only detectable at the 16S

Spatial patterns of bacterial and archaeal communities along the Romanche Fracture Zone (tropical Atlantic).

The composition of prokaryotic communities was determined in the meso- and bathypelagic waters funneled through the Romanche Fracture Zone (RFZ, 2°7'S, 31°79'W to 0°6'N, 14°33'W) in the tropical Atlantic. Distinct water masses were identified based on their physical and chemical characteristics. The bacterial and archaeal communities were depth-stratified with a total of 116 and 25 operational taxonomic units (OTUs), respectively, distributed among the distinct water masses as revealed by terminal restriction fragment length polymorphism, and cloning and sequencing. The relative abundance of Thaumarchaeota, determined by catalyzed reporter deposition-fluorescence in situ hybridization, was significantly higher in deeper layers (Antarctic Bottom Water, AABW, > 4000 m depth), contributing up to 31% to the total prokaryotic community, than in the mesopelagic and lower euphotic layer. Although the contribution of SAR11 to bacterial abundance did not increase with depth, SAR202, SAR324, SAR406 and Alteromonas did increase with depth. Terminal restriction fragment length polymorphism analysis revealed successional changes in the bacterial and archaeal community composition of the North Atlantic Deep Water (NADW) with a passage time through the RFZ of c. 4 months but not in the under- and overlying water masses. Our results indicate that specific water masses harbor distinct bacterial and archaeal communities and that the prokaryotic community of the NADW undergoes successional changes in this conduit between the western and eastern Atlantic basin. Apparently, in the absence of major input of organic matter to specific deep-water masses, the indigenous prokaryotic community adapts to subtle physical and biogeochemical changes in the water mass within a time frame of weeks, similar to the reported seasonal changes in surface water prokaryotic communities.