Using Amplicon Sequencing To Characterize and Monitor Bacterial Diversity in Drinking Water Distribution Systems

Drinking water assessments use a variety of microbial, physical, and chemical indicators to evaluate water treatment efficiency and product water quality. However, these indicators do not allow the complex biological communities, which can adversely impact the performance of drinking water distribution systems (DWDSs), to be characterized. Entire bacterial communities can be studied quickly and inexpensively using targeted metagenomic amplicon sequencing. Here, amplicon sequencing of the 16S rRNA gene region was performed alongside traditional water quality measures to assess the health, quality, and efficiency of two distinct, full-scale DWDSs: (i) a linear DWDS supplied with unfiltered water subjected to basic disinfection before distribution and (ii) a complex, branching DWDS treated by a four-stage water treatment plant (WTP) prior to disinfection and distribution. In both DWDSs bacterial communities differed significantly after disinfection, demonstrating the effectiveness of both treatment regimes. However, bacterial repopulation occurred further along in the DWDSs, and some end-user samples were more similar to the source water than to the postdisinfection water. Three sample locations appeared to be nitrified, displaying elevated nitrate levels and decreased ammonia levels, and nitrifying bacterial species, such as Nitrospira, were detected. Burkholderiales were abundant in samples containing large amounts of monochloramine, indicating resistance to disinfection. Genera known to contain pathogenic and fecal-associated species were also identified in several locations. From this study, we conclude that metagenomic amplicon sequencing is an informative method to support current compliance-based methods and can be used to reveal bacterial community interactions with the chemical and physical properties of DWDSs.     

Changes in Bacterial Communities of the Marine Sponge Mycale laxissima on Transfer into Aquaculture

The changes in bacterial communities associated with the marine sponge Mycale laxissima on transfer to aquaculture were studied using culture-based and molecular techniques. M. laxissima was maintained alive in flowthrough and closed recirculating aquaculture systems for 2 years and 1 year, respectively. The bacterial communities associated with wild and aquacultured sponges, as well as the surrounding water, were assessed using 16S rRNA gene clone library analysis and denaturing gradient gel electrophoresis (DGGE). Bacterial richness and diversity were measured using DOTUR computer software, and clone libraries were compared using S-LIBSHUFF. DGGE analysis revealed that the diversity of the bacterial community of M. laxissima increased when sponges were maintained in aquaculture and that bacterial communities associated with wild and aquacultured M. laxissima were markedly different than those of the corresponding surrounding water. Clone libraries of bacterial 16S rRNA from sponges confirmed that the bacterial communities changed during aquaculture. These communities were significantly different than those of seawater and aquarium water. The diversity of bacterial communities associated with M. laxissima increased significantly in aquaculture. Our work shows that it is important to monitor changes in bacterial communities when examining the feasibility of growing sponges in aquaculture systems because these communities may change. This could have implications for the health of sponges or for the production of bioactive compounds by sponges in cases where these compounds are produced by symbiotic bacteria rather than by the sponges themselves.     

Pyrosequencing analysis of bacterial communities in drinking water biofilters receiving influents of different types

Biological activated carbon (BAC) filters are commonly used in the world for improvement of drinking water quality. The indigenous microbiota in BAC filters can play a crucial role in reduction or biotransformation of contaminants. Molecular analysis can enhance our understanding of ecological functions of the microbial communities in drinking water BAC filters. In this study, three laboratory-scale drinking water BAC filters receiving influents of different types were constructed. Differences of bacterial communities in the three BAC filters were characterized using 454 pyrosequencing analysis. Pyrosequencing analysis illustrated the usefulness in elucidating the bacterial community structure in drinking water biofilter. High bacterial diversity in granular activated carbon (GAC) samples from each BAC biofilter was observed. Proteobacteria was the largest bacterial phylum in each GAC sample, with a marked shift of the proportions of Alpha-, Beta-, and Gammaproteobacteria. The levels of dissolved organic carbon and ammonia nitrogen in the influents could affect the bacterial diversity and community composition in the BAC biofilters. This work might add some new insights into microbial community and its influential factors in drinking water biofilters.     

Microbial diversity in marine biofilms along a water quality gradient on the Great Barrier Reef

Microbial communities are potential indicators for water quality as they respond rapidly to environmental changes. In the Whitsunday Islands, Australia, microbial biofilm communities from two offshore islands were compared to those from two inshore islands subjected to poor water quality. Biofilm community composition was characterized using three culture-independent molecular techniques. The clone libraries indicated high genetic diversity, with somewhat higher scores in the offshore sites (57%) compared to the inshore sites (41%). The majority of microbes in the biofilms were related to Alphaproteobacteria (39.8%), Gammaproteobacteria (14.1%), Bacteroidetes (13.2%), diatoms (8.3%) and Cyanobacteria (3.9%). Redundancy analysis (RDA) for the CARD-FISH data showed distinct microbial assemblages between offshore and inshore communities. Additionally, 5 out of 13 water quality parameters (DIN, Chla, POP, TSS and POC) explained a significant amount of variation in the microbial communities and high values of these were associated with inshore communities. Analysis of variance (ANOVA) indicated that Cyanobacteria (p = 0.01), Bacteroidetes (p = 0.04) and to some extent Alphaproteobacteria (p = 0.07), were significantly more abundant in the offshore biofilm communities. Principal Component Analysis (PCA) of DGGE data showed clear grouping of cyanobacterial communities into inshore and offshore communities. Reasons for community shifts in the bacterial lineages are currently not resolved. One possible causative factor may be that autotrophic primary producers are more dominant in offshore sites due to the higher light availability as well as the limitation by DIN. The trends found in this study are the bases for more detailed research on microbial indicator species for changes in water quality.     

Phyllosphere Bacterial Community of Floating Macrophytes in Paddy Soil Environments as Revealed by Illumina High-Throughput Sequencing

The phyllosphere of floating macrophytes in paddy soil ecosystems, a unique habitat, may support large microbial communities but remains largely unknown. We took Wolffia australiana as a representative floating plant and investigated its phyllosphere bacterial community and the underlying driving forces of community modulation in paddy soil ecosystems using Illumina HiSeq 2000 platform-based 16S rRNA gene sequence analysis. The results showed that the phyllosphere of W. australiana harbored considerably rich communities of bacteria, with Proteobacteria and Bacteroidetes as the predominant phyla. The core microbiome in the phyllosphere contained genera such as Acidovorax, Asticcacaulis, Methylibium, and Methylophilus. Complexity of the phyllosphere bacterial communities in terms of class number and α-diversity was reduced compared to those in corresponding water and soil. Furthermore, the bacterial communities exhibited structures significantly different from those in water and soil. These findings and the following redundancy analysis (RDA) suggest that species sorting played an important role in the recruitment of bacterial species in the phyllosphere. The compositional structures of the phyllosphere bacterial communities were modulated predominantly by water physicochemical properties, while the initial soil bacterial communities had limited impact. Taken together, the findings from this study reveal the diversity and uniqueness of the phyllosphere bacterial communities associated with the floating macrophytes in paddy soil environments.     

Pyrosequencing Analysis of the Bacterial Community in Drinking Water Wells

Wells used for drinking water often have a large biomass and a high bacterial diversity. Current technologies are not always able to reduce the bacterial population, and the threat of pathogen proliferation in drinking water sources is omnipresent. The environmental conditions that shape the microbial communities in drinking water sources have to be elucidated, so that pathogen proliferation can be foreseen. In this work, the bacterial community in nine water wells of a groundwater aquifer in Northern Mexico were characterized and correlated to environmental characteristics that might control them. Although a large variation was observed between the water samples, temperature and iron concentration were the characteristics that affected the bacterial community structure and composition in groundwater wells. Small increases in the concentration of iron in water modified the bacterial communities and promoted the growth of the iron-oxidizing bacteria Acidovorax. The abundance of the genera Flavobacterium and Duganella was correlated positively with temperature and the Acidobacteria Gp4 and Gp1, and the genus Acidovorax with iron concentrations in the well water. Large percentages of Flavobacterium and Pseudomonas bacteria were found, and this is of special concern as bacteria belonging to both genera are often biofilm developers, where pathogens survival increases.     

Dynamic changes in microbial community structure in farming pond water and their effect on the intestinal microbial community profile in juvenile common carp (Cyprinus carpio L.)

Water quality parameter dynamics, gut, sediment and water bacteria communities were studied to understand the environmental influence on the gut microbial community of a new strain of Huanghe common carp. A total of 3,384,078 raw tags and 5105 OTUs were obtained for the gut, water and sediment bacteria. The water quality had a stronger influence on the water bacteria community than gut and sediment bacteria communities. The ambient water quality parameters also significantly influenced the water and sediment bacteria communities. Comparing the gut, sediment, and water microbial communities, a relationship was found among them. However, gut bacteria were more closely related to sediment bacterial communities than to water bacteria communities. The results showed that the top three bacterial taxa were identical in gut and sediment samples in the early days of rearing. Interestingly, bacterial communities in the carp gut, water, and sediment had different adaptabilities to variations in environmental factors.     

Host-Specificity and Dynamics in Bacterial Communities Associated with Bloom-Forming Freshwater Phytoplankton

Many freshwater phytoplankton species have the potential to form transient nuisance blooms that affect water quality and other aquatic biota. Heterotrophic bacteria can influence such blooms via nutrient regeneration but also via antagonism and other biotic interactions. We studied the composition of bacterial communities associated with three bloom-forming freshwater phytoplankton species, the diatom Aulacoseira granulata and the cyanobacteria Microcystis aeruginosa and Cylindrospermopsis raciborskii. Experimental cultures incubated with and without lake bacteria were sampled in three different growth phases and bacterial community composition was assessed by 454-Pyrosequencing of 16S rRNA gene amplicons. Betaproteobacteria were dominant in all cultures inoculated with lake bacteria, but decreased during the experiment. In contrast, Alphaproteobacteria, which made up the second most abundant class of bacteria, increased overall during the course of the experiment. Other bacterial classes responded in contrasting ways to the experimental incubations causing significantly different bacterial communities to develop in response to host phytoplankton species, growth phase and between attached and free-living fractions. Differences in bacterial community composition between cyanobacteria and diatom cultures were greater than between the two cyanobacteria. Despite the significance, major differences between phytoplankton cultures were in the proportion of the OTUs rather than in the absence or presence of specific taxa. Different phytoplankton species favoring different bacterial communities may have important consequences for the fate of organic matter in systems where these bloom forming species occur. The dynamics and development of transient blooms may also be affected as bacterial communities seem to influence phytoplankton species growth in contrasting ways.     

Bacterial Communities of Surface Mixed Layer in the Pacific Sector of the Western Arctic Ocean during Sea-Ice Melting

From July to August 2010, the IBRV ARAON journeyed to the Pacific sector of the Arctic Ocean to monitor bacterial variation in Arctic summer surface-waters, and temperature, salinity, fluorescence, and nutrient concentrations were determined during the ice-melting season. Among the measured physicochemical parameters, we observed a strong negative correlation between temperature and salinity, and consequently hypothesized that the melting ice decreased water salinity. The bacterial community compositions of 15 samples, includicng seawater, sea-ice, and melting pond water, were determined using a pyrosequencing approach and were categorized into three habitats: (1) surface seawater, (2) ice core, and (3) melting pond. Analysis of these samples indicated the presence of local bacterial communities; a deduction that was further corroborated by the discovery of seawater- and ice-specific bacterial phylotypes. In all samples, the Alphaproteobacteria, Flavobacteria, and Gammaproteobacteria taxa composed the majority of the bacterial communities. Among these, Alphaproteobacteria was the most abundant and present in all samples, and its variation differed among the habitats studied. Linear regression analysis suggested that changes in salinity could affect the relative proportion of Alphaproteobacteria in the surface water. In addition, the species-sorting model was applied to evaluate the population dynamics and environmental heterogeneity in the bacterial communities of surface mixed layer in the Arctic Ocean during sea-ice melting.     

Bacterial communities in water and sediment shaped by paper mill pollution and indicated bacterial taxa in sediment in Daling River

Effluents from paper mills are highly toxic and are a major source of aquatic pollution. In this study, we collected water and sediment samples to examine the microbial communities using denaturing gradient gel electrophoresis and identified bacterial taxa greatly affected by paper mill pollution using next-generation sequencing data. Our results indicated that bacterial communities in downstream sediments were similar to those in paper mill discharge sites, indicating obvious effects of pollution, while bacterial communities in downstream water samples showed similar profiles to those in upstream sites, both being quite different from the bacterial communities in paper mill discharge sites. This was possibly because of the short contact period. In addition, bacterial communities in the estuary were quite different from those in other water and sediment samples, which was owing to the special habitat type. Considering the storage of paper mill pollutants in sediment and the significant effect on shifts in bacterial communities, we selected Clostridia and Epsilonproteobacteria at the class level and Fusibacter and Desulfobulbus at the genus level as bacterial indicators of paper mill pollution. To monitor the remediation of polluted aquatic environments, we propose Sphingobacteria, Alphaproteobacteria, Actinobacteria, Subdivision3, Planctomycetacia and Verrucomicrobiae at the class level and Bacillus, Steroidobacter, Nocardioides, Terrimonas, Pirellula and Methylibium at the genus level.     

Development of a small-scale bioreactor method to monitor the molecular diversity and environmental impacts of bacterial biofilm communities from an acid mine drainage impacted creek

Shamokin Creek is a tributary of the Susquehanna River in central Pennsylvania that is heavily impacted by the acid mine drainage (AMD) caused by the oxidation of pyrite from the region's extensive anthracite coal mining industry. Recent studies have begun to characterize the microbial communities present in this and other AMD-impacted waters, but varying environmental conditions have complicated attempts to determine the ecological impacts of individual bacterial species within these communities. This study developed a small-scale biofilm reactor protocol that allowed us to simultaneously monitor the development of bacterial biofilm communities in AMD-impacted creek collected water using terminal restriction fragment length polymorphism (T-RFLP) analysis of 16S rRNA genes, while assessing the impacts that the developing biofilms were having on water quality. Our analysis confirmed that the diversity and composition of these small in situ biofilm communities could be monitored using molecular methods, and indicated the possible presence of many taxa frequently found in AMD environments, including Sulfobacillus, Nitrospira, Desulfovibrio, Geobacter, and Leptothrix species. A significant increase in the total sulfate was observed in the bioreactor, and as most likely due to the accumulation of sulfur-oxidizing bacteria such as Sulfobacillus in the biofilms. This system will allow us to study the microbial ecology of Shamokin Creek through controlled experiments that will ultimately integrate microscopic, molecular, physiological and chemical analyses, and that can be utilized to develop more effective and cost-efficient environmental remediation techniques for AMD-impacted areas.     

Epiphytic bacterial communities of the alga <Emphasis Type="Italic">Fucus vesiculosus</Emphasis> in oil-contaminated water areas of the Barents Sea

Taxonomic compositions of epiphytic bacterial communities in water areas differing in levels of oil pollution were revealed. In total, 82 bacterial genera belonging to 16 classes and 11 phyla were detected. All detected representatives of epiphytic bacterial communities belonged to the phyla Actinobacteria, Bacteroidetes, Planctomycetes, Proteobacteria, Verrucomicrobia, Acidobacteria, Cyanobacteria, Firmicutes, and Fusobacteria and candidate division TM7. The ratio of the phyla in the communities varied depending on the levels of oil pollution. New data on taxonomic composition of uncultivated epiphytic bacterial communities of Fucus vesiculosus were obtained.     

Vertical distribution of bacterial community is associated with the degree of soil organic matter decomposition in the active layer of moist acidic tundra

The increasing temperature in Arctic tundra deepens the active layer, which is the upper layer of permafrost soil that experiences repeated thawing and freezing. The increasing of soil temperature and the deepening of active layer seem to affect soil microbial communities. Therefore, information on soil microbial communities at various soil depths is essential to understand their potential responses to climate change in the active layer soil. We investigated the community structure of soil bacteria in the active layer from moist acidic tundra in Council, Alaska. We also interpreted their relationship with some relevant soil physicochemical characteristics along soil depth with a fine scale (5 cm depth interval). The bacterial community structure was found to change along soil depth. The relative abundances of Acidobacteria, Gammaproteobacteria, Planctomycetes, and candidate phylum WPS-2 rapidly decreased with soil depth, while those of Bacteroidetes, Chloroflexi, Gemmatimonadetes, and candidate AD3 rapidly increased. A structural shift was also found in the soil bacterial communities around 20 cm depth, where two organic (upper Oi and lower Oa) horizons are subdivided. The quality and the decomposition degree of organic matter might have influenced the bacterial community structure. Besides the organic matter quality, the vertical distribution of bacterial communities was also found to be related to soil pH and total phosphorus content. This study showed the vertical change of bacterial community in the active layer with a fine scale resolution and the possible influence of the quality of soil organic matter on shaping bacterial community structure.     

Distinct patterns of marine bacterial communities in the South and North Pacific Oceans

The study of oceanic microbial communities is crucial for our understanding of the role of microbes in terms of biomass, diversity and ecosystem function. In this study, 16S rRNA gene tag pyrosequencing was used to investigate change in bacterial community structure between summer and winter water masses from Gosung Bay in the South Sea of Korea and Chuuk in Micronesia, located in the North and South Pacific Oceans, respectively. Summer and winter sampling from each water mass revealed highly diverse bacterial communities, containing ~900 Operational Taxonomic Units (OTUs). The microbial distribution and highly heterogeneous composition observed at both sampling sites were different from those of most macroorganisms. The bacterial communities in the seawater at both sites were most abundant in Proteobacteria during the summer in Gosung and in Bacterioidetes during the winter. The proportion of Cyanobacteria was higher in summer than in winter in Chuuk and similar in Gosung. Additionally, the microbial community during summer in Gosung was significantly different from other communities observed based on the unweighted UniFrac distance. These data suggest that in both oceanic areas sampled, the bacterial communities had distinct distribution patterns with spatially- and temporally-heterogeneous distributions.     

Bacterioplankton community structure in the Arctic waters as revealed by pyrosequencing of 16S rRNA genes

Fjords and open oceans are two typical marine ecosystems in the Arctic region, where glacial meltwater and sea ice meltwater have great effects on the bacterioplankton community structure during the summer season. This study aimed to determine the differences in bacterioplankton communities between these two ecosystems in the Arctic region. We conducted a detailed census of microbial communities in Kongsfjorden (Spitsbergen) and the Chukchi Borderland using high-throughput pyrosequencing of the 16S rRNA gene. Gammaproteobacteria and Bacteroidetes were the dominant members of the bacterioplankton community in Kongsfjorden. By contrast, the most abundant bacterial groups in the surface seawater samples from the Chukchi Borderland were Alphaproteobacteria and Actinobacteria. Differences in bacterial communities were found between the surface and subsurface waters in the investigation area of the Chukchi Borderland, and significant differences in bacterial community structure were also observed in the subsurface water between the shelf and deep basin areas. These results suggest the effect of hydrogeographic conditions on bacterial communities. Ubiquitous phylotypes found in all the investigated samples belonged to a few bacterial groups that dominate marine bacterioplankton communities. The sequence data suggested that changes in environmental conditions result in abundant rare phylotypes and reduced amounts of other phylotypes.     

Effect of Soybean Monoculture on the Bacterial Communities Associated with Cysts of Heterodera glycines

The soybean cyst nematode (SCN), Heterodera glycines, can cause significant reductions in soybean yield and quality in many parts of the world. Natural biological control may play an important role in regulating SCN population. In this study the bacterial communities associated with SCN cysts obtained from fields under different lengths of soybean monoculture were explored. Soil samples were collected in 2010 and 2011 from six fields that had been used for soybean monoculture for 2 to 41 yr. SCN population densities were determined and bacterial communities from SCN cysts were investigated by Biolog and PCR-DGGE methods. SCN population densities initially increased in the first 5 yr of soybean monoculture but then declined steeply as years of soybean monoculture increased. Catabolic diversity of bacterial communities associated with cysts tended to decline as number of years of monoculture increased. Some specific PCR-DGGE bands, mainly representing Streptomyces and Rhizobium, were obtained from the cysts collected from the long-term monoculture fields. Principal component analysis of Biolog and PCR-DGGE data revealed that bacterial communities associated with cysts could be divided into two groups: those from cysts obtained from shorter (< 8 yr) vs. longer (> 8 yr) monoculture. This research demonstrates that the composition of the bacterial communities obtained from SCN cysts changes with length of soybean monoculture; the suppressive impact of these bacterial communities to SCN is yet to be determined.     

Bacterial Communities Associated with Production Facilities of Two Newly Drilled Thermogenic Natural Gas Wells in the Barnett Shale (Texas, USA)

We monitored the bacterial communities in the gas–water separator and water storage tank of two newly drilled natural gas wells in the Barnett Shale in north central Texas, using a 16S rRNA gene pyrosequencing approach over a period of 6?months. Overall, the communities were composed mainly of moderately halophilic and halotolerant members of the phyla Firmicutes and Proteobacteria (classes Βeta-, Gamma-, and Epsilonproteobacteria) in both wells at all sampling times and locations. Many of the observed lineages were encountered in prior investigations of microbial communities from various fossil fluid formations and production facilities. In all of the samples, multiple H₂S-producing lineages were encountered; belonging to the sulfate- and sulfur-reducing class Deltaproteobacteria, order Clostridiales, and phylum Synergistetes, as well as the thiosulfate-reducing order Halanaerobiales. The bacterial communities from the separator and tank samples bore little resemblance to the bacterial communities in the drilling mud and hydraulic-fracture waters that were used to drill these wells, suggesting the in situ development of the unique bacterial communities in such well components was in response to the prevalent geochemical conditions present. Conversely, comparison of the bacterial communities on temporal and spatial scales suggested the establishment of a core microbial community in each sampled location. The results provide the first overview of bacterial dynamics and colonization patterns in newly drilled, thermogenic natural gas wells and highlights patterns of spatial and temporal variability observed in bacterial communities in natural gas production facilities.     

Bacterial Community Composition in Central European Running Waters Examined by Temperature Gradient Gel Electrophoresis and Sequence Analysis of 16S rRNA Genes

The bacterial community composition in small streams and a river in central Germany was examined by temperature gradient gel electrophoresis (TGGE) with PCR products of 16S rRNA gene fragments and sequence analysis. Complex TGGE band patterns suggested high levels of diversity of bacterial species in all habitats of these environments. Cluster analyses demonstrated distinct differences among the communities in stream and spring water, sandy sediments, biofilms on stones, degrading leaves, and soil. The differences between stream water and sediment were more significant than those between sites within the same habitat along the stretch from the stream source to the mouth. TGGE data from an entire stream course suggest that, in the upper reach of the stream, a special suspended bacterial community is already established and changes only slightly downstream. The bacterial communities in water and sediment in an acidic headwater with a pH below 5 were highly similar to each other but deviated distinctly from the communities at the other sites. As ascertained by nucleotide sequence analysis, stream water communities were dominated by Betaproteobacteria (one-third of the total bacteria), whereas sediment communities were composed mainly of Betaproteobacteria and members of the Fibrobacteres/Acidobacteria group (each accounting for about 25% of bacteria). Sequences obtained from bacteria from water samples indicated the presence of typical cosmopolitan freshwater organisms. TGGE bands shared between stream and soil samples, as well as sequences found in bacteria from stream samples that were related to those of soil bacteria, demonstrated the occurrence of some species in both stream and soil habitats. Changes in bacterial community composition were correlated with geographic distance along a stream, but in comparisons of different streams and rivers, community composition was correlated only with environmental conditions.     

Bacterial Communities in Semen from Men of Infertile Couples: Metagenomic Sequencing Reveals Relationships of Seminal Microbiota to Semen Quality

Some previous studies have identified bacteria in semen as being a potential factor in male infertility. However, only few types of bacteria were taken into consideration while using PCR-based or culturing methods. Here we present an analysis approach using next-generation sequencing technology and bioinformatics analysis to investigate the associations between bacterial communities and semen quality. Ninety-six semen samples collected were examined for bacterial communities, measuring seven clinical criteria for semen quality (semen volume, sperm concentration, motility, Kruger''s strict morphology, antisperm antibody (IgA), Atypical, and leukocytes). Computer-assisted semen analysis (CASA) was also performed. Results showed that the most abundant genera among all samples were Lactobacillus (19.9%), Pseudomonas (9.85%), Prevotella (8.51%) and Gardnerella (4.21%). The proportion of Lactobacillus and Gardnerella was significantly higher in the normal samples, while that of Prevotella was significantly higher in the low quality samples. Unsupervised clustering analysis demonstrated that the seminal bacterial communities were clustered into three main groups: Lactobacillus, Pseudomonas, and Prevotella predominant group. Remarkably, most normal samples (80.6%) were clustered in Lactobacillus predominant group. The analysis results showed seminal bacteria community types were highly associated with semen health. Lactobacillus might not only be a potential probiotic for semen quality maintenance, but also might be helpful in countering the negative influence of Prevotella and Pseudomonas. In this study, we investigated whole seminal bacterial communities and provided the most comprehensive analysis of the association between bacterial community and semen quality. The study significantly contributes to the current understanding of the etiology of male fertility.