Next Article in Journal
A Novel Marine Pathogen Isolated from Wild Cunners (Tautogolabrus adspersus): Comparative Genomics and Transcriptome Profiling of Pseudomonas sp. Strain J380
Next Article in Special Issue
Microbial Communities Present in Hydrothermal Sediments from Deception Island, Antarctica
Previous Article in Journal
Stool Serology: Development of a Non-Invasive Immunological Method for the Detection of Enterovirus-Specific Antibodies in Congo Gorilla Faeces
Previous Article in Special Issue
Shift from Carbon Flow through the Microbial Loop to the Viral Shunt in Coastal Antarctic Waters during Austral Summer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 4
10.3390/microorganisms9040811
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bacterial Communities Associated with Poa annua Roots in Central European (Poland) and Antarctic Settings (King George Island)

by
Anna Znój
1,2,†,
Jakub Grzesiak
1,*,†,
Jan Gawor
3,
Robert Gromadka
3 and
Katarzyna J. Chwedorzewska
4
1
Department of Antarctic Biology, Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawińskiego 5A, 02-106 Warsaw, Poland
2
Botanical Garden-Center for Biological Diversity Conservation, Polish Academy of Sciences, Prawdziwka 2, 02-973 Warsaw, Poland
3
Environmental Laboratory of DNA Sequencing and Synthesis, Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawińskiego 5A, 02-106 Warsaw, Poland
4
Department of Agronomy, Warsaw University of Life Sciences-SGGW, Nowoursynowska 166, 02-787 Warsaw, Poland
*
Author to whom correspondence should be addressed.
†
These authors equally contributed to this work.
Microorganisms 2021, 9(4), 811; https://doi.org/10.3390/microorganisms9040811
Submission received: 12 March 2021 / Revised: 8 April 2021 / Accepted: 9 April 2021 / Published: 12 April 2021
(This article belongs to the Special Issue Polar Microbes)
Browse Figures
Versions Notes

Abstract

:
Poa annua (annual bluegrass) is one of the most ubiquitous grass species in the world. In isolated regions of maritime Antarctica, it has become an invasive organism threatening native tundra communities. In this study, we have explored and compared the rhizosphere and root-endosphere dwelling microbial community of P. annua specimens of maritime Antarctic and Central European origin in terms of bacterial phylogenetic diversity and microbial metabolic activity with a geochemical soil background. Our results show that the rhizospheric bacterial community was unique for each sampling site, yet the endosphere communities were similar to each other. However, key plant-associated bacterial taxa such as the Rhizobiaceae family were poorly represented in Antarctic samples, probably due to high salinity and heavy metal concentrations in the soil. Metabolic activity in the Antarctic material was considerably lower than in Central European samples. Antarctic root endosphere showed unusually high numbers of certain opportunistic bacterial groups, which proliferated due to low competition conditions. Thirteen bacterial families were recognized in this study to form a core microbiome of the P. annua root endosphere. The most numerous were the Flavobacteriaceae, suspected to be major contributors to the ecological success of annual bluegrass, especially in harsh, Antarctic conditions.

1. Introduction

Poa annua L. (annual bluegrass) is one of the most ubiquitous grass species in the world. It can be found growing on every continent, having established populations even in the high Arctic and Antarctica [1,2,3]. Its extensive adaptability to a wide range of environmental conditions makes it a stubborn weed and pioneer species. P. annua has developed a series of features that ensure its ecological success. It can deal with biotic and abiotic stressors as well as habitat instability by compact growth habit, a high fraction of biomass allocated into belowground organs, and the presence of both annual and perennial life forms [4,5,6]. It also displays huge phenotypic and genotypic variability, making it difficult in the past to determine its evolutionary origin [1].
These adaptive traits caused P. annua to be a highly competitive, invasive species in many isolated environments of unique biodiversity [7]. Maritime Antarctica, due to its harsh climate, seemed very unlikely to be colonized by annual bluegrass. However, in 1953, specimens of P. annua were recorded near a ruined whaling station on Deception Island, South Shetland Islands, maritime Antarctica [8]. Since then, it was discovered in several other sites within the region [9]. Its expansion was documented in detail in the vicinity of H. Arctowski Polish Antarctic Station (King George Island, South Shetland Islands), where first sightings of the species were made during the austral summer of 1985–1986, next to the stations living quarters [10,11]. Available historical data point toward Poland-exported soils for greenhouse cultivation purposes as the most likely source of viable P. annua diaspores and its subsequent introduction [8]. However, molecular analyses indicate a constant inflow of fresh bluegrass genetic material into the introduced population, mostly due to ongoing human activity [2,12]. Primarily occurring in human-disturbed sites associated with the station’s infrastructure, P. annua has expanded into recently deglaciated glacier forelands, even entering into Antarctic tundra communities [11,13].
There is a consensus within the scientific community about the positive and negative influence of the associated microorganisms on such plant traits as nutrient acquisition, growth, biomass production, and flowering time but also stress tolerance and disease resistance [14,15]. Those multifactorial plant-microbe interactions take place mostly in two major rhizocompartments within the boundaries of the soil-root system mass [16]. The rhizosphere compartment comprises the soil that is in closest proximity to plant roots, where the microbiome is strongly shaped by root exudates and mucilage and contains a metabolically and phylogenetically diverse community [17]. The endosphere, that is, the root interior, hosts a plant-specific microbiome due to selective entry mechanisms and further biochemical exclusion-enrichment processes [18]. In this respect, invasive plant species that thrive in suboptimal conditions are presumed to rely not only on genetic, physiological, and morpho-anatomical adaptations but also on a mutual relationship with the resident microbiome, which enhances their stress tolerance and contributes to the overall plant strength [19].
In this study, we investigated P. annua root-associated microbiomes sampled in Central Europe (Poland) and maritime Antarctica (King George Island, South Shetland Islands) in terms of bacterial phylogenetic diversity (16S rRNA gene-targeted next-generation amplicon sequencing), metabolic traits (community-level physiological profiling by Biolog Ecoplates) and a geochemical background (soil biogenic component concentrations). The main aim of this investigation was to compare microbial communities residing in the rhizosphere and the root endosphere of P. annua specimens that grow in their native European settings with those growing in maritime Antarctica as invasive plant species. Our hypothesis states that the differences between microbial communities of Central European and maritime Antarctic P. annua roots are severe in the rhizospheric compartment yet are minor in the root endosphere, both metabolically and phylogenetically.
This is the first study that sheds light on the microbiome associated with an invasive plant in Antarctica and the first to provide insight into P. annua associated bacterial community structure by a culture-independent approach.

2. Materials and Methods

2.1. Sites and Sampling

Antarctic material sampling was carried out in the vicinity of the H. Arctowski Polish Antarctic Station, located at the western shore of Admiralty Bay, King George Island (South Shetlands Islands, maritime Antarctica). Samples were collected during the austral summer season of 2017–2018 from three sites where P. annua growth was noted. The maximal temperature on the day of sampling was 6.7 °C. Central European material sampling was carried out at the Botanical Garden of the Polish Academy of Sciences (Powsin, city, Poland) in the fall of 2018 (Table 1). The maximal temperature on the day of sampling was 17.3 °C. Several specimens (4–6 per site) were collected with the root-adjacent soil with the use of sterile tools into sterile plastic containers and transported frozen (−20 °C) to the laboratory in the Institute of Biochemistry and Biophysics, Polish Academy of Sciences (IBB PAS). Additionally, bulk soil samples from those sites were gathered in triplicates (approx. 1.5 kg per site) for component analysis and transported in the same conditions.

2.2. Bacterial Extraction

Bacterial cells were extracted from fifteen (three per site) individual P. annua specimen rhizospheric soil and roots. The following method was devised based on the findings of [20] regarding the separation of prokaryotic cells from mineral and organic debris and the guidelines provided by [21] regarding root-associated microbe isolation. To analyze the microbiome of the root-adjacent soil, a sample of the soil was carefully removed from between the roots with a sterile spatula onto a pre-sterilized aluminum foil piece. Approx. 1 g of the soil was weighed and placed in a 50 mL conical tube containing 20 mL of sterile and cool (4 °C) dilution liquid composed of 0.9% (w/v) saline (NaCl) and 10 mM tetrasodium pyrophosphate (Na4P2O7). The suspension was then shaken for 30 min in a TornadoTM Vortexer at 2000 rpm at 4 °C. The tubes were then placed in a VWR Ultrasonic Cleaner USC-TH filled with chilled water and sonicated for 60 s. The tubes were vortexed afterward for 30 s to suspend detached cells. After brief centrifugation (1 min; 1000 rpm; 4 °C), the suspension was submitted to metabolic fingerprinting by the Biolog Ecoplate technique. To detach the rest of the adhering soil, the root system was washed in 60 mL of sterile NaCl/Na4P2O7 solution by shaking for 30 min in the aforementioned shaker (1000 rpm; 4 °C) and then rinsed 3 times in 5 mL of the same sterile and cooled solution by vortexing. Washed roots were sterilized by incubation in a cooled 10% hydrogen peroxide (H2O2) solution for 5 min, then rinsed 3 times with sterile NaCl/Na4P2O7 solution. The so surface-sterilized roots were placed in a pre-cooled sterile mortar. Two and a half ml of sterile NaCl/Na4P2O7 was added with 0.6 g of sterile, sharp garnet sand (lysing matrix A) and gently ground with a pestle, allowing the sharp angular garnet pieces to comminute the roots to an amorphous pulp. The pulp was transferred to a 50 mL conical tube containing 20 mL of sterile and cool (4 °C) NaCl/Na4P2O7 solution and submitted to the above-mentioned procedure (shaking, ultrasonication and vortexing). The resulting supernatant suspension was submitted to metabolic fingerprinting by the Biolog Ecoplate technique and DNA extraction.

2.3. DNA Extraction and Targeted 16S rRNA Gene Amplicon Sequencing

Rhizosphere soil DNA was extracted using the PowerSoil® DNA isolation kit (QIAGEN GmbH, Hilden, Germany) according to manufacturer protocol. An approx. 0.2 g of soil was used in triplicates. DNA solutions were kept at 4 °C for further analysis. The dilution liquid containing endosphere bacteria was passed through a sterile 47 mm Whatman polycarbonate filter (0.22 μm pore size). The DNA from the filter-trapped bacteria was extracted using the PowerWater® DNA isolation kit (QIAGEN, GmbH, Hilden, Germany) according to manufacturer protocol and kept at 4 °C. This resulted in 72 DNA samples. The phylogenetic study was performed by targeted sequencing and analysis of the prokaryotic 16S ribosomal RNA gene. A fragment of the 16S rRNA gene containing the V3 and V4 variable regions was amplified using gene-specific primers: 16S_V3-F and 16S_V4-R positions 341-357F and 785-805R, respectively, according to Escherichia coli 16S rRNA gene reference sequence [22]. Illumina Nextera XT overhang adapter nucleotide sequences were included in addition to the 16S rRNA gene-specific sequences, which allowed sample indexing and pooling. Each PCR amplification was conducted in triplicates using KAPA HiFi PCR kit (Roche, Basel, Switzerland) in a final volume of 20 μL per reaction according to the manufacturer’s instructions. Obtained PCR products were pooled into 10 samples (2 rhizocompartments × 5 sampling sites) in equimolar ratio and indexed using Nextera XT barcodes (Illumina, San Diego, CA, USA). Amplicon libraries were sequenced on Illumina MiSeq instrument (Illumina, San Diego, CA, USA) in the DNA Sequencing and Oligonucleotide Synthesis Laboratory (IBB, PAS). Sequencing was conducted in paired-end mode (2 × 300 bp) with the use of a v.3 (600 cycles) chemistry cartridge, which allowed the generation of long paired reads fully covering 16S V3–V4 amplicons.

2.4. Phenotype Fingerprinting with Biolog EcoPlate™

The EcoPlate Biolog assays assess the ability of a mixed microbial community to use any of 31 carbon compounds as the sole carbon source (plus a single control well with no-carbon). Microbial communities were characterized for their ability to catabolize 10 different carbohydrates, 9 carboxylic and acetic acids, 4 polymers, 6 amino acids, and 2 amines [23]. Root-associated bacterial suspensions were adjusted with sterile 0.9% saline to optical transmittance of 0.9. One hundred microliter aliquots of each suspension were added to each well of EcoPlate microplates (Biolog Inc., Hayward, CA, USA). The plates were incubated in darkness at 10 °C for Antarctic samples and 18 °C for European material. The temperatures were chosen to accommodate the activity range of the respective microbial communities: psychrophiles and psychrotrophes for Antarctic material, psychrotrophes and mesophiles for European material [24]. The color development (absorbance) was read at 590 nm (A590) in a Varioscan plate reader (Thermofisher Scientific, Waltham, MA, USA), and cellular respiration was measured kinetically by determining the colorimetric reduction of tetrazolium dye. Data were collected approximately twice a week over a 65 day period. The prolonged incubation of EcoPlates was based on our previous observations [25,26,27]. Data from the forty-second day (Antarctic samples) and twenty-first day (European samples) of incubation were used as there was no further color development after this date. Final absorbance data were first blanked against the time zero reading and then blanked against the respective control well containing no-carbon source. Readings that had the A590 value of 0.25 or higher were scored as a positive EcoPlate response (PER).

2.5. Measurement of Soil Components

Soil pH (in 1 M KCl) and salinity (in double-distilled water (ddH2O)) were measured with a CPC-411 Elmetron™ multiparameter probe according to [28]. Phosphates and nitrates were determined spectrophotometrically in a Shimadzu UV 1601 spectrophotometer and in an Epoll-Eco 20 spectrophotometer respectively. Other elements were determined by atomic absorption spectroscopy [29].

2.6. Data Analysis

Raw sequencing data were cleaned, aligned, and classified automatically by the EzBioCloud platform using the PKSSU4.0 database [30]. Chimeric, low quality, and non-target (chloroplast, mitochondrial, and archaeal) amplicons were automatically discarded. The operational taxonomic unit was defined as a group of sequences that exhibit greater than 97% similarity to each other. Illumina reads were deposited in the NCBI Sequence Read Archive (SRA) as BioProject PRJNA678861. All results were compiled using Excel (MS Office) 2016 for Windows. A two-sample t-test was applied to compare different data sets. Variance within the sets was assessed using the f-test beforehand. Correlations between biological and geochemical parameters were calculated using Pearson’s correlation coefficient. Principal component analysis was performed using the singular value decomposition method. Data visualization and statistical analysis has been performed using the R software (R v.4.0.2) and the following packages: ggplot2, fmsb, Hmisc, ggpubr, corrplot, and autoplot [31].

3. Results

3.1. Soil Components

The nutrient-poor European soil sample (S1) was characterized by high calcium content (5000 mg/100 g) and relatively high pH (7.8), with subsequently low magnesium (38.4 mg/100 g), manganese (52.8 mg/kg), zinc (5.9 mg/kg), copper (2.7 mg/kg), iron and sodium content (4.9 mg/100 g) as well as relatively low nitrate (17.4 mg/100 g), labile phosphorus (6 mg/100 g), and labile potassium (36.1 mg/100 g). The fertilized European soil sample (S2) had very high nitrate, labile phosphorus, and potassium levels (27.5, 219.2, and 251.8 mg per 100 g of soil, respectively). pH in this sample was neutral (7.0) with the highest manganese (218.1 mg/kg) in all examined samples. Antarctic samples were characterized by overall low nitrate (7.7–16.3 mg/100 g), labile phosphorus (4.9–15.1 mg/100 g), and potassium concentration (31.1–55.8 mg/100 g) values, high salinity (0.4 g NaCl/L), and high magnesium (205–250.7 mg/100 g) and iron contents (4269.8–4670.2 mg/kg). Zinc was low in the glacial foreland sample (9.2 mg/kg) (Table 2).

3.2. Diversity Indices

The highest diversity as assessed by operational taxonomic unit numbers (OTUs) were noted for the fertile Central European soil sample P2S (OTU = 7879), followed by Central European soil sample P1S (OTU = 6813). Antarctic soil samples displayed lower bacterial phylogenetic diversity. Sample P3S showed the highest diversity among them (OTU = 4498), whereas the proglacial soil sample P4S the lowest diversity (OTU = 1322). Bacterial communities residing within the roots displayed lower diversity (av. 2512.4 OTU, sd. 841) than the corresponding rhizosphere soil (av. 5165 OTU, sd. 2083). The most diverse was the bacterial community in soil sample P1R from grass specimens growing in nutrient-poor European soil (OTU = 3694). Samples P2R showed lower values (OTU = 2620). Root-residing bacteriome was diverse in plants growing near the Arctowski Antarctic station: P3R (OTU = 2465) and P5R (OTU = 2461), the one residing in plants growing in postglacial soil (P4R) showed relatively low diversity (OTU = 1322). Other diversity indices are presented in Supplementary File S1. The number of positive EcoPlate responses in rhizospheric soil samples varied between 26 (sample P5S) and 31 (sample P1S), whereas for the root samples, it varied between 17 (sample P5R) and 31 (sample P4R) (Figure 1).

3.3. Bacterial Phylogenetic Diversity

Twelve major (>1%) phyla were found in the rhizosphere and root interior of the P. annua specimens (Figure 2). Proteobacterial sequences were the most numerous in investigated samples. According to their percentage contribution, Proteobacteria have been enriched in the root-interior compartment in comparison to the respective rhizosphere soil community. The enrichment is approx. 2-fold in every case analyzed, being it European (P1S/P1R = 25.32%/48.14%; P2S/P2R = 30.5%/42.52%) or Antarctic material (P3S/P3R = 26.60%/51.83%; P4S/P4R = 26.80%/55.04%; P5S/P5R = 26.79%/46.87%). Sequences belonging to the Bacteroidetes phylum also were highly abundant in the sampling material. Their enrichment in the root interior was apparent only in some cases (P1S/P1R = 14.53%/24.32%; P2S/P2R = 17.04%/38.68%; P5S/P5R = 13.18%/24.12%), whereas not so in others (P3S/P3R = 23.18%/24.96%; P4S/P4R = 31.90%/33.84%). Acidobacterial presence was reduced in root interior compared to the rhizosphere soil (P1S/P1R = 8.01%/2.77%; P2S/P2R = 8.95%/0.65%; P3S/P3R = 8.72%/2.25%; P4S/P4R = 3.29%/0.19%; P5S/P5R = 7.98%/1.81%). Similarly, Planctomycetes abundance also diminished inside roots (P1S/P1R = 14.38%/2.90%; P2S/P2R = 7.39%/0.61%; P3S/P3R = 5.28%/1.09%; P4S/P4R = 2.92%/0.18%; P5S/P5R = 4.67%/0.65%). Noteworthy are the low Verrucomicrobia numbers in the root interior of specimens growing in postglacial soils (P4R = 0.71%) and high Firmicutes numbers in the corresponding rhizosphere soil (P4S = 13%).
Relative abundance of family-rank groups revealed several discrepancies, but also similarities between European and Antarctic P. annua root-associated bacterial communities (Figure 3, Supplementary File S2). In the rhizosphere community, the most noticeable was the difference in the relative abundances of the family Flavobacteriaceae (Europe—av. 3.6%; Antarctica—av. 13.9%) and also of the Planctomycetaceae (Europe—av. 4.5%; Antarctica—av. 1.1%). Some families displayed a comparable average relative abundance in both sample types, most notably the Chitinophagaceae (Europe—av. 4.5%; Antarctica—av. 4.1%) and the Chthoniobacteraceae (Europe—av. 3.0%; Antarctica—av. 2.9%). High relative abundance in the rhizosphere of a few families was site-specific. Site 4 (proglacial terrain) was especially rich in the following: Clostridiaceae (12.7%), Oxalobacteraceae (8.3%), and Pseudomonadaceae (7.4%). The endosphere community harbored high relative abundances of the family Flavobacteriaceae, both for the European (av. 17.5%) and the Antarctic (av. 17.7%) material. Major differences between European and Antarctic root-endosphere communities could be noticed in the relatively low abundances for the latter of the families: Rhizobiaceae (Europe—av. 6.6%; Antarctica—av. 2.4%) and Comamonadaceae (Europe—av. 7.3%; Antarctica—av. 3.3%). Similar to the rhizosphere, some endospheric families displayed unusually high relative abundances only at specific sites. Site 3 was particularly rich in Pseudomonadaceae sequences (25.8%), Site 4 in Oxalobacteraceae sequences (21.0%), and Site 5 in Sphingomonadaceae (15.3%) and Sphingobacteriaceae sequences (10.2%).

3.4. Community-Level Physiological Profiling by Biolog EcoPlates

Microbial community response intensity in Biolog Ecoplates was considerably higher for the European samples than in those obtained from Antarctic material (Figure 4, Supplementary File S3). In order to be comparable, the data needed to be normalized in terms of intensity value (values within a particular sample type were divided by the highest value in that sample). The highest absorbance value at 590 nm (A590) obtained in the European material was A590 = 4.69, whereas, for the Antarctic samples, A590 = 2.82. D-mannitol was the most intensely metabolized carbon source by all examined microbial communities, especially by root-inhabiting communities from European material (A590 = 4.49–4.69). Apart from that, severe discrepancies in the metabolism of other compounds have been observed. Microbial community residing in nutrient-poor rhizospheric soil (P1S) was most efficient in metabolizing D-cellobiose (A590 = 4.27) and α-D-lactose (A590 = 4.03) as well as D-xylose (A590 = 3.81), L-asparagine (A590 = 3.59), i-erythritol (A590 = 3.60), and 4-hydroxy benzoic acid (A590 = 3.52). The corresponding root community (P1R) efficiently metabolized L-asparagine (A590 = 4.27), γ-hydroxybutyric acid (A590 = 3.57), and also D-cellobiose (A590 = 3.33). The fertile European rhizospheric soil community (P2S) provided the strongest responses to L-arginine (A590 = 2.44), L-serine (A590 = 2.27), L-phenylalanine (A590 = 1.98), α-cyclodextrin (A590 = 1.97), N-acetyl-D-glucosamine (A590 = 2.05), D-cellobiose (A590 = 2.44), and D-malic acid (A590 = 1.94). The corresponding root community metabolized preferably: L-arginine (A590 = 4.21) and L-asparagine (A590 = 4.62) but also Tween40 (A590 = 3.45). Apart from metabolizing D-mannitol, the Antarctic rhizosphere and root communities had an affinity for i-erythritol (A590 = 2.30–2.82), except for the P5S sample (A590 = 1.18). α-cyclodextrin and α-D-lactose, except P3R in both cases, was also readily metabolized by the Antarctic microbial community (A590 = 1.38–2.46 and A590 = 1.57–2.30, respectively). Efficient N-acetyl-D-glucosamine and D-cellobiose metabolism were apparent only in sample P5S (A590 = 1.83 and A590 = 2.44, respectively).

3.5. Correlations between Biological and Geochemical Data

Pearson’s correlation coefficient between the percentile abundance of bacterial families within the rhizosphere and the chemical composition of the soil revealed several significant (p < 0.05) negative correlations (Figure 5A, Supplementary File S4). Salinity displayed the highest number of negative correlations with bacterial family-rank groups: Thermoleophilaceae (r = −0.98, p = 0.002), Rhizobiaceae (r = −0.98, p = 0.003), Planctomycetaceae (r = −0.98, p = 0.003), Caulobacteraceae (r = −0.98, p = 0.004), Verrucomicrobiaceae (r = −0.96, p = 0.01), Micromonosporaceae (r = −0.95, p = 0.01), Bradyrhizobiaceae (−0.95, p = 0.01), Pseudonocardiaceae (r = −0.92, p = 0.03), Iamiaceae (r = −0.92, p = 0.03), Cytophagaceae (r = −0.91, p = 0.03), Blastocatellaceae (r = −0.90, p = 0.04), Gemmataceae (r =−0.90, p = 0.04), and Sinobacteraceae (r = −0.89, p = 0.04). Iron concentration showed negative correlations with the relative abundance of the following families: Rhizobiaceae (r = −0.99, p = 0.001), Thermoleophilaceae (r = −0.94, p = 0.02), Planctomycetaceae (r = −0.93, p = 0.02), Caulobacteraceae (r = −0.92, p = 0.03), and Verrucomicrobiaceae (r = −0.91, p = 0.03). Copper concentrations negatively correlated with the relative abundance of the following bacterial families: Rhizobiaceae (r = −0.96, p = 0.01), Sinobacteraceae (r = −0.93, p = 0.02), Verrucomicrobiaceae (r = −0.90, p = 0.04), and Planctomycetaceae (r = −0.89, p = 0.04). Pearson’s correlation coefficient between the percentile abundance of bacterial families of the P. annua root endosphere highlighted several significant correlations (Figure 5B, Supplementary File S4). Most were obtained for the Flavobacteriaceae, which were negatively correlated with: Acidimicrobiaceae (r = −0.99, p = 0.002), Chthoniobacteraceae (r = −0.98, p = 0.004), Micropepsaceae (r = −0.97, p = 0.01), Polyangiaceae (r = −0.96, p = 0.01), Chitinophagaceae (r = −0.92, p = 0.03), Vicinamibacteraceae (r = −0.92, p = 0.03), Blastocatellaceae (r = −0.89, p = 0.04), and Bradyrhizobiaceae (r = −0.89, p = 0.04). The Pseudomonadaceae abundance displayed negative correlations with the contribution of: Caulobacteraceae (r = −0.97, p = 0.01), Rhizobiaceae (r = −0.94, p = 0.02), Sphingobacteriaceae (r = −0.93, p = 0.02), and Intrasporangiaceae (r = −0.90, p = 0.04). The Oxalobacteraceae displayed a negative correlation with the Iamiaceae (r = −0.93, p = 0.02).

3.6. Principal Component Analysis

Principal component analysis (PCA) based on the abundance of sequences of a family-rank taxon showed a tight clustering of the root-dwelling bacterial communities, both of European and Antarctic origin (Figure 6A). Soil bacterial communities were vastly different, with the proglacial sample P4S bearing the closest resemblance to root communities. PCA based on the Biolog Ecoplate responses showed that European root-interior microbial communities were the most similar despite the differences in their corresponding rhizospheric communities. Antarctic root communities bore more resemblance to their respective rhizosphere communities than to each other or the corresponding European samples (Figure 6B). PCA of soil chemistry data showed a great discrepancy between the nutrient-poor (P1S) and the nutrient-rich European rhizospheric soils (P2S). Antarctic samples clustered closely together, showing great similarity (Figure 6C). PCA clustering made using the combined phylogenetic and functional (Biolog Ecoplate) data revealed four distinct groups: European rhizospheric soil community (P1S, P2S), Antarctic rhizospheric soil community (P3S, P4S, P5S), European root community (P1R, P2R), and Antarctic root community (P3R, P4R, P5R) (Figure 6D).

3.7. Significant Differences between Microbial Parameters

The t-test calculations showcased significant (p < 0.05) differences between Central European and Antarctic rhizospheric microbial parameters (Figure 7). The relative abundance of the following taxa was higher in the European samples: Planctomycetes, Caulobacteraceae, Planctomycetaceae, Rhizobiaceae, Sinobacteraceae, Thermoleophilaceae, and Verrucomicrobiaceae. The catabolic intensity of several compounds was also higher in the Central European rhizosphere samples, most notably for: β-methyl-D-glucoside, phenylethylamine, glucose-1-phosphate, γ-hydroxybutyric acid, itaconic acid, D-galacturonic acid, and D-xylose. Antarctic material had a significantly higher sequence abundance of the Saccharibacteria phylum and the actinobacterial family of Microbacteriaceae. Significant differences between the root endosphere communities highlighted higher abundances/intensities of several features in the Central European samples. In those samples, a higher relative abundance of the following bacterial families were noted: Comamonadaceae, Cytophagaceae, Micromonosporaceae, Nocardioidaceae, and Rhizobiaceae. As for compound use intensity, it was true for the following substrates: β-methyl-D-glucoside, D-galactonic acid γ-lactone, D-galacturonic acid, L-asparagine, Tween 80, N-acetyl-D-glucosamine, γ-hydroxybutyric acid, D-cellobiose, glucose-1-phosphate, D-malic acid, and putrescine. The only feature significantly higher in the Antarctic material was i-erythritol use intensity. Most of the significant differences were apparent between the rhizosphere and the endosphere of P. annua. The root interior showed significant enrichment in the relative abundance of the following taxa: Proteobacteria, Microbacteriaceae, Sphingobacteriaceae, Hyphomicrobiaceae, Rhizobiaceae, Bacteroidetes, and Comamonadaceae while also displaying significantly lower catabolism intensities of L-phenylalanine, α-cyclodextrin, glycogen, and D-cellobiose.

3.8. Poa annua Core Microbiome

Thirteen bacterial families had an average relative abundance higher than the average abundance of family-rank group in the root endosphere and are therefore considered to be the core microbiome taxa of P. annua in this study (Figure 8). The Flavobacteriaceae had the highest abundance on average (8.0–28.7%, Me = 16.2) followed by Sphingobacteriaceae (3.1–10.2%, Me = 6.8%), Pseudomonadaceae (1.4–22.8%, Me = 6.7%), Rhizobiaceae (0.8–6.9%, Me = 4.9%), Comamonadaceae (1.2–7.3%, Me = 4.6%), Oxalobacteraceae (1.5–21.0%, Me = 4.1%), Microbacteriaceae (2.8–7.9%, Me = 3.5%), Sphingomonadaceae (2.4–15.3%, Me = 3.3%), Caulobacteraceae (1.3–3.3%, Me = 3.0%), Chitinophagaceae (0.7–4.2%, Me = 2.4%), Xanthomonadaceae (0.1–6.0%, Me = 2.0%), Hyphomicrobiaceae (1.0–3.7%, Me = 1.7%), and Cytophagaceae (0.9–2.8%, Me = 1.6%).

4. Discussion

The multiphasic approach applied in this study revealed several phenomena within the P. annua root-associated microbiome.
Rhizospheric soils examined in this study each contained a unique bacterial community in terms of phylogenetic diversity, as shown in the Principal Component Analysis. This diversity was higher in the European rhizospheric soil samples than in the Antarctic ones. Antarctic ecosystems were repeatedly proven to be less diverse and less complex than those in lower latitudes [32,33,34]. Geographical isolation was often made responsible for the low diversity of Antarctic macro-organisms [35,36]. However, this kind of isolation may not directly apply to bacteria, as it was discovered that live bacterial cells can be transported over long geographic distances [37]. Regarding rhizospheric microbial communities, their diversity can be largely controlled by host diversity and abundance [38], which in the case of Antarctic flora is indeed limited due to geographic isolation but also due to harsh climatic conditions [39]. The specific geochemical characteristics of the examined Antarctic soils may be even more responsible for the low bacterial diversity [40]. Those soils shared some features such as high salinity, sodium, magnesium, and heavy metal concentrations. This is mainly due to sea spray influence and the volcanic history of this site [41]. These characteristics, especially salinity, showed significant, negative correlations with members belonging to several soil community bacterial families, including the plant growth-promoting Rhizobiaceae. Salt heavy metal sensitivity, especially that of free rhizobial cells, has been extensively studied, as it impacts the biomass production of several agricultural plants [42,43,44]. Although there were no significant differences in the relative sequence abundance of the majority of bacterial phyla between European and Antarctic rhizospheric soils, they did significantly differ in the abundance of Planctomycetes and Saccharibacteria. European rhizosphere had a higher contribution of Planctomycetes, members of which are connected to plant biomass decomposition, but also to high soil calcium concentrations, which was noted in those soils [45]. Furthermore, European samples displayed significantly higher contributions of key taxa, such as the Rhizobiaceae, a family that contains plant-beneficial nitrogen fixers, and the Sinobacteraceae-harboring rhizospheric ammonia oxidizers [46,47]. The members of other families such as Caulobacteraceae and Verrucomicrobiaceae, known for grass-related rhizospheric competence, have also been found in significantly higher amounts in those samples [48,49]. Antarctic rhizospheric soils had on average a higher contribution of Saccharibacteria sequences. Based on a limited amount of reports about the members of this phylum, they can be suspected to be scavengers or even parasites rather than active biomass degraders [50]. Members of the Microbacteriaceae displayed significantly higher amounts in the P. annua rhizosphere growing in the Antarctic. Rhizospheric enrichment in Microbacteriaceae has rarely been observed but could be influenced by climatic factors rather than soil composition [51]. A curious anomaly was noted in the Antarctic rhizospheric soils of postglacial origin. The postglacial rhizospheric soil sample contained high numbers of sequences belonging to the Firmicutes phylum, the majority of which were identified as Clostridiaceae, a family containing anaerobic, endospore-forming bacteria [52]. This can be connected to the recent deglaciation event [53]. Subglacial habitats have been speculated to be mostly anoxic. The experiment of [54] confirmed that in such conditions, members of the Clostridiaceae could be considerably enriched from Antarctic, supraglacial materials. Furthermore, the spores of the Clostridiaceae can persist in Antarctic settings even after the conditions have changed to oxic [26]. Besides higher phylogenetical diversity, European rhizospheric communities also displayed significantly higher use abilities for several types of compounds, such as carbohydrates, amino and carboxylic acids, lipids, and amines. Together with the high phylogenetic diversity, this suggests a plethora of different bacterial niches and high assemblage complexity, hinting at the presence of a rich reservoir of rhizosphere-competent bacteria [55,56].
The rhizosphere is considered the direct reservoir of root endophytes, albeit a selective entry mechanism exists and is well-described [57,58]. In this respect, the phylogenetic diversity in the endosphere was considerably lower than that in the rhizosphere. Furthermore, rhizosphere diversity had a direct impact on the endosphere diversity. This seems to be a staple when examining plant-associated bacterial communities as it was observed for wild plants [59] as well as a field [60] or greenhouse cultivars [61]. The endospheric communities of P. annua assessed in this study showed much phylogenetic similarity to each other when examined collectively with the rhizosphere community. This was probably due to the depletion in the root of several major soil-dwelling high-ranking taxa such as Acidobacteria and Planctomycetes, along with some other low abundance phyla. Sequences of the Proteobacteria and Bacteroidetes were significantly enriched in the endosphere. Despite this apparent similarity, several differences in the composition of the European and Antarctic endophytic communities could be noticed. Members of the family Rhizobiaceae and also Comamonadaceae were detected to largely contribute to European root endospheric bacteriomes but were significantly less numerous in Antarctic material. Those families hold nitrogen-fixing members, and their presence can largely increase plant biomass production [62,63]. In addition, the endosphere, as well as the rhizosphere of Antarctic P. annua specimens, were low in some of the most pivotal actinobacterial taxa found in well-developed soils in lower latitudes, namely the mycelium forming members of the families Nocardioidaceae and Micromonosporaceae, known for their antimicrobial activities [64,65]. On a family-rank level, there were no significantly enriched groups in the Antarctic endosphere. However, a few bacterial families showed abnormally high contribution in the Antarctic samples, namely: Oxalobacteraceae, Pseudomonadaceae, and Sphingomonadaceae. This is consistent with the findings of [66], where Pseudomonadaceae and Sphingomonadaceae were among high abundance taxa in the endosphere of native Antarctic plants: Colobanthus quitensis and Deschampsia antarctica. However, in those communities, the family Enterobacteriaceae was present in considerable amounts, unlike in the endosphere of examined here P. annua. The Pseudomonadaceae showed negative correlations with several key plant endophytic taxa such as the Rhizobiaceae, Intrasporangiaceae, Caulobacteraceae, and Sphingobacteriaceae, whereas the Oxalobacteraceae were negatively correlated with the Iamiaceae. Microbial interactions are usually very complex, albeit, within the root interior, they are more simplified due to the selection-mediated diversity impoverishment. In this respect, this could be a case of antagonisms or niche overlap. A recent study on the endophytic Oxalobacteraceae member Massilia sp. revealed its copiotrophic and r-strategy based lifestyle while also being competition-sensitive [67]. Such an explanation could be plausible in the Antarctic habitat, especially in the low diversity proglacial site, which houses opportunistic pioneer species [68]. The principal component analysis of the catabolic responses showed an interesting feature of the Antarctic bacterial endophytic community. While the Central European endophytic communities differed from the corresponding rhizospheric communities (and being more metabolically diverse), the Antarctic endophyte communities showed more similarity to the matching rhizospheric communities. This points toward Antarctic soil bacterial communities being catabolically fixed on certain abundant and stable nutrient sources, and the root-incorporated bacteria might not be true endophyte specialists. Furthermore, i-erythritol catabolism intensity was significantly higher in the Antarctic endophytic community. This sugar alcohol, contrary to the popular Mannitol, is not produced by plants but by green algae [69]. Antarctic soils house different species of aeroterrestrial algae, such as the widespread in the Antarctic region Prasiola crispa [10]. This compound’s high catabolism potential in the Antarctic P. annua roots could mean aeroterrestrial algae-associated bacteria incorporation into the endophytic community. Despite the aforementioned dissimilarities between root communities, a core bacterial endophytobiome of P. annua can be extrapolated. The family Flavobacteriaceae has been on average the most abundant in the root interior, despite a wide abundance range. Only in recent years have the members of this family been recognized as important contributors to the root endosphere community [70] and were found to largely contribute to the endosphere of Antarctic native grass D. antarctica [71]. It is suspected that a metabolically diverse pool of Flavobacterium spp. is constantly present in the root interior, shifting in abundance patterns in response to environmental conditions and the physiological state of the host plant [70]. The latter is supported by the findings of [72], where plant-associated factors dictated flavobacterial presence. Investigations on temperate region invasive plants’ root microbiota indicate that they increase the fitness of the host in non-native settings mainly by enhancing the acquisition of nutrients such as phosphorous and nitrogen [73]. However, in polar regions, the alleviation of cold-induced stresses may be of key importance for alien plant establishment [74]. In this respect, studies on the bacterially mediated cold resistance of agricultural crops point toward the action of ACC deaminase [75,76]. This microbially produced enzyme lowers the concentration of plant stress hormone ethylene, a substance that inflicts a significant reduction in plant growth and development [77]. Certain strains of the genus Flavobacterium have displayed ACC deaminase activity [78]. In consequence, Flavobacteriaceae bacteria may play a crucial role in P. annua adaptation to Antarctic conditions by modifying stress responses of the plant during adaptation [77].

5. Conclusions

In conclusion, bacterial root-associated communities of P. annua differed both phylogenetically and metabolically between those of Central European and those of maritime Antarctic origin. Several key plant-beneficial bacterial groups were less abundant in the Antarctic material, probably due to the geochemical makeup of the soil. Some bacterial families displayed unusually high abundance in the root endosphere of the Antarctic P. annua specimens. They most likely contained competition-sensitive opportunists that proliferated in the absence of antagonistic microbes. Nonetheless, an endophytic core microbiome could be assumed consisting of 13 bacterial families belonging to the Proteobacteria, Bacteroidetes, and Actinobacteria phyla. The Flavobacteriaceae family was the most numerous and most likely to positively influence the adaptation of P. annua to Antarctic conditions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9040811/s1, Table S1: Alpha-diversity indices for the bacterial communities associated with Poa annua roots, Table S2: Relative abundance heatmap of sequences identified on a family-rank taxonomic level, Table S3: Heatmap displaying Poa annua rhizosphere and root community responses on Biolog Ecoplates, Table S4: Correlations of rhizospheric family-rank sequence abundance data and soil chemistry, Table S5: Correlations between root endosphere family-rank sequence abundance data.

Author Contributions

A.Z.: conceptualization, validation, investigation, resources, writing—review and editing, project administration, funding acquisition. J.G. (Jakub Grzesiak): conceptualization, software, formal analysis, data curation, writing—original draft, visualization. J.G. (Jan Gawor): validation, investigation, data curation, writing—review and editing. R.G.: writing—review and editing. K.J.C.: conceptualization, resources, writing—review and editing, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Science Center, Poland (Grant 2016/21/N/NZ9/01536).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Illumina reads were deposited in the NCBI Sequence Read Archive (SRA) as BioProject PRJNA678861.

Acknowledgments

Samples and data were obtained due to the scientific facility of the H. Arctowski Polish Antarctic Station.

Conflicts of Interest

Authors declare no conflict of interest.

References

  1. Mao, Q.; Huff, D.R. The evolutionary origin of Poa annua L. Crop. Sci. 2012, 52, 1910–1922. [Google Scholar] [CrossRef]
  2. Chwedorzewska, K.J. Poa annua L. in Antarctic: Searching for the source of introduction. Polar Biol. 2008, 31, 263–268. [Google Scholar] [CrossRef]
  3. Heide, O.M. Flowering responses of contrasting ecotypes of Poa annua and their putative ancestors Poa infirma and Poa supina. Ann. Bot. 2001, 87, 795–804. [Google Scholar] [CrossRef] [Green Version]
  4. Galera, H.; Rudak, A.; Czyż, E.A.; Chwedorzewska, K.J.; Znój, A.; Wódkiewicz, M. The role of the soil seed store in the survival of an invasive population of Poa annua at Point Thomas Oasis, King George Island, maritime Antarctica. Glob. Ecol. Conserv. 2019, 19, e00679. [Google Scholar] [CrossRef]
  5. Rudak, A.; Galera, H.; Znój, A.; Chwedorzewska, K.J.; Wódkiewicz, M. Seed germination and invasion success of Poa annua L. in Antarctica. Acta Soc. Bot. Pol. 2018, 87. [Google Scholar] [CrossRef]
  6. Williams, L.K.; Shaw, J.D.; Sindel, B.M.; Wilson, S.C.; Kristiansen, P. Longevity, growth and community ecology of invasive Poa annua across environmental gradients in the subantarctic. Basic Appl. Ecol. 2018, 29, 20–31. [Google Scholar] [CrossRef] [Green Version]
  7. Gardarin, A.; Dürr, C.; Colbach, N. Which model species for weed seedbank and emergence studies? A review. Weed Res. 2009, 49, 117–130. [Google Scholar] [CrossRef]
  8. Skottsberg, C. Antarctic vascular plants. Bot. Tidsskr. 1954, 51, 330–338. [Google Scholar]
  9. Chwedorzewska, K.J.; Giełwanowska, I.; Olech, M.; Molina-Montenegro, M.A.; Wódkiewicz, M.; Galera, H. Poa annua L. in the maritime Antarctic: An overview. Polar Rec. 2015, 51, 637–643. [Google Scholar] [CrossRef]
  10. Olech, M. Human impact on terrestrial ecosystems in west Antarctica. In Proceedings of the NIPR Symposium on Polar Biology; National Institute of Polar Research: Tokyo, Japan, 1996; pp. 299–306. [Google Scholar]
  11. Olech, M.; Chwedorzewska, K.J. Short note: The first appearance and establishment of an alien vascular plant in natural habitats on the forefield of a retreating glacier in Antarctica. Antarct. Sci. 2011, 23, 153. [Google Scholar] [CrossRef]
  12. Chwedorzewska, K.J.; Bednarek, P.T. Genetic and epigenetic variation in a cosmopolitan grass Poa annua from Antarctic and Polish populations. Pol. Polar Res. 2012, 33, 63–80. [Google Scholar] [CrossRef] [Green Version]
  13. Galera, H.; Chwedorzewska, K.J.; Łapiński, S.; Kowalska, M.E.; Pasik, M.; Rajner, M.; Czyż, E.; Bylina PWódkiewicz, M. First step to eradication of Poa annua L. from Arctowski oasis (King George Island, South Shetlands, Antarctica). Polar Biol. 2017, 40, 939–945. [Google Scholar] [CrossRef] [Green Version]
  14. De Zelicourt, A.; Al-Yousif, M.; Hirt, H. Rhizosphere microbes as essential partners for plant stress tolerance. Mol. Plant 2013, 6, 242–245. [Google Scholar] [CrossRef] [Green Version]
  15. Wei, Z.; Yang, T.; Friman, V.P.; Xu, Y.; Shen, Q.; Jousset, A. Trophic network architecture of root-associated bacterial communities determines pathogen invasion and plant health. Nat. Commun. 2015, 6, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Turner, T.R.; Ramakrishnan, K.; Walshaw, J.; Heavens, D.; Alston, M.; Swarbreck, D.; Poole, P.S. Comparative metatranscriptomics reveals kingdom level changes in the rhizosphere microbiome of plants. ISME J. 2013, 7, 2248–2258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Sasse, J.; Martinoia, E.; Northen, T. Feed your friends: Do plant exudates shape the root microbiome? Trends Plant Sci. 2018, 23, 25–41. [Google Scholar] [CrossRef] [Green Version]
  18. Hartmann, A.; Rothballer, M.; Schmid, M. Lorenz Hiltner, a pioneer in rhizosphere microbial ecology and soil bacteriology research. Plant Soil 2008, 312, 7–14. [Google Scholar] [CrossRef]
  19. Trognitz, F.; Hackl, E.; Widhalm, S.; Sessitsch, A. The role of plant–microbiome interactions in weed establishment and control. FEMS Microbiol. Ecol. 2016, 92, fiw138. [Google Scholar] [CrossRef] [Green Version]
  20. Lunau, M.; Lemke, A.; Walther, K.; Martens-Habbena, W.; Simon, M. An improved method for counting bacteria from sediments and turbid environments by epifluorescence microscopy. Environ. Microbiol. 2005, 7, 961–968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Szymańska, S.; Płociniczak, T.; Piotrowska-Seget, Z.; Złoch, M.; Ruppel, S.; Hrynkiewicz, K. Metabolic potential and community structure of endophytic and rhizosphere bacteria associated with the roots of the halophyte Aster tripolium L. Microbiol. Res. 2016, 182, 68–79. [Google Scholar] [CrossRef] [PubMed]
  22. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glöckner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41, e1. [Google Scholar] [CrossRef]
  23. Weber, K.P.; Legge, R.L. One-dimensional metric for tracking bacterial community divergence using sole carbon source utilization patterns. J. Microbiol. Methods 2009, 79, 55–61. [Google Scholar] [CrossRef] [PubMed]
  24. Moyer, C.L.; Morita, R.Y. Psychrophiles and psychrotrophs. eLS 2007. [Google Scholar] [CrossRef]
  25. Grzesiak, J.; Zdanowski, M.K.; Górniak, D.; Świątecki, A.; Aleksandrzak-Piekarczyk, T.; Szatraj, K.; Nieckarz, M. Microbial community changes along the ecology glacier ablation zone (King George Island, Antarctica). Polar Biol. 2015, 38, 2069–2083. [Google Scholar] [CrossRef] [Green Version]
  26. Grzesiak, J.; Kaczyńska, A.; Gawor, J.; Żuchniewicz, K.; Aleksandrzak-Piekarczyk, T.; Gromadka, R.; Zdanowski, M.K. A smelly business: Microbiology of Adélie penguin guano (Point Thomas rookery, Antarctica). Sci. Total Environ. 2020, 714, 136714. [Google Scholar] [CrossRef]
  27. Grzesiak, J.; Woltyńska, A.; Zdanowski, M.K.; Górniak, D.; Świątecki, A.; Olech, M.A.; Aleksandrzak-Piekarczyk, T. Metabolic fingerprinting of the Antarctic cyanolichen Leptogium puberulum–associated bacterial community (Western Shore of Admiralty Bay, King George Island, Maritime Antarctica). Microb. Ecol. 2021, 1–12. [Google Scholar] [CrossRef]
  28. Jones, J.B. Laboratory Guide for Conducting Soil Tests and Plant Analysis (No. BOOK); CRC Press: Boca Raton, FL, USA, 2001. [Google Scholar]
  29. Tatur, A.; Myrcha, A.; Niegodzisz, J. Formation of abandoned penguin rookery ecosystems in the maritime Antarctic. Polar Biol. 1997, 17, 405–417. [Google Scholar] [CrossRef]
  30. Yoon, S.H.; Ha, S.M.; Kwon, S.; Lim, J.; Kim, Y.; Seo, H.; Chun, J. Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 2017, 67, 1613. [Google Scholar] [CrossRef] [PubMed]
  31. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 6 June 2020).
  32. Wynn-Williams, D.D. Antarctic microbial diversity: The basis of polar ecosystem processes. Biodivers. Conserv. 1996, 5, 1271–1293. [Google Scholar] [CrossRef]
  33. Convey, P.; Chown, S.L.; Clarke, A.; Barnes, D.K.; Bokhorst, S.; Cummings, V.; Griffiths, H.J. The spatial structure of Antarctic biodiversity. Ecol. Monogr. 2014, 84, 203–244. [Google Scholar] [CrossRef] [Green Version]
  34. Znój, A.; Chwedorzewska, K.J.; Androsiuk, P.; Cuba-Diaz, M.; Giełwanowska, I.; Koc, J.; Korczak-Abshire, M.; Grzesiak, J.; Zmarz, A. Rapid environmental changes in the Western Antarctic Peninsula region due to climate change and human activity. Appl. Ecol. Environ. Res. 2017, 15, 525–539. [Google Scholar] [CrossRef]
  35. Magalhães, C.; Stevens, M.I.; Cary, S.C.; Ball, B.A.; Storey, B.C.; Wall, D.H.; Ruprecht, U. At limits of life: Multidisciplinary insights reveal environmental constraints on biotic diversity in continental Antarctica. PLoS ONE 2012, 7, e44578. [Google Scholar] [CrossRef] [PubMed]
  36. Barnes, D.K.; Hodgson, D.A.; Convey, P.; Allen, C.S.; Clarke, A. Incursion and excursion of Antarctic biota: Past, present and future. Glob. Ecol. Biogeogr. 2006, 5, 121–142. [Google Scholar] [CrossRef]
  37. Smith, D.J.; Griffin, D.W.; McPeters, R.D.; Ward, P.D.; Schuerger, A.C. Microbial survival in the stratosphere and implications for global dispersal. Aerobiologia 2011, 27, 319–332. [Google Scholar] [CrossRef]
  38. Berendsen, R.L.; Pieterse, C.M.; Bakker, P.A. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012, 17, 478–486. [Google Scholar] [CrossRef]
  39. Alberdi, M.; Bravo, L.A.; Gutiérrez, A.; Gidekel, M.; Corcuera, L.J. Ecophysiology of Antarctic vascular plants. Physiol. Plant 2002, 115, 479–486. [Google Scholar] [CrossRef] [PubMed]
  40. Chong, C.W.; Pearce, D.A.; Convey, P.; Yew, W.C.; Tan, I.K.P. Patterns in the distribution of soil bacterial 16S rRNA gene sequences from different regions of Antarctica. Geoderma 2012, 181, 45–55. [Google Scholar] [CrossRef]
  41. Navas, A.; López-Martínez, J.; Casas, J.; Machín, J.; Durán, J.J.; Serrano, E.; Mink, S. Soil characteristics on varying lithological substrates in the South Shetland Islands, maritime Antarctica. Geoderma 2008, 144, 123–139. [Google Scholar] [CrossRef] [Green Version]
  42. Hua SS, T.; Tsai, V.Y.; Lichens, G.M.; Noma, A.T. Accumulation of amino acids in Rhizobium sp. strain WR1001 in response to sodium chloride salinity. Appl. Environ. Microbiol. 1982, 44, 135–140. [Google Scholar] [PubMed]
  43. Elsheikh, E.A. Effects of salt on rhizobia and bradyrhizobia: A review. Ann. Appl. Biol. 1998, 132, 507–524. [Google Scholar] [CrossRef]
  44. Broos, K.; Beyens, H.; Smolders, E. Survival of rhizobia in soil is sensitive to elevated zinc in the absence of the host plant. Soil Biol. Biochem. 2005, 37, 573–579. [Google Scholar] [CrossRef]
  45. Buckley, D.H.; Huangyutitham, V.; Nelson, T.A.; Rumberger, A.; Thies, J.E. Diversity of planctomycetes in soil in relation to soil history and environmental heterogeneity. Appl. Environ. Microbiol. 2006, 72, 4522–4531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Alves, L.M.C.; de Souza, J.A.M.; de Mello Varani, A.; de Macedo Lemos, E.G. The Family Rhizobiaceae. In The Prokaryotes: Alphaproteobacteria and Betaproteobacteria; Rosenberg, E., DeLong, E.F., Lory, S., Stackebrandt, E., Thompson, F., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 419–437. [Google Scholar]
  47. Liu, Y.; Ludewig, U. Nitrogen dependent bacterial community shifts in root, rhizome and rhizosphere of nutrient efficient Miscanthus x giganteus from long term field trials. GCB Bioenergy 2019, 11, 1334–1347. [Google Scholar] [CrossRef] [Green Version]
  48. Schlatter, D.C.; Yin, C.; Hulbert, S.; Paulitz, T.C. Core rhizosphere microbiomes of dryland wheat are influenced by location and land use history. Appl. Environ. Microbiol. 2020, 86, e02135-19. [Google Scholar] [CrossRef]
  49. Da Rocha, U.N.; van Elsas, J.D.; van Overbeek, L.S. Real-time PCR detection of Holophagae (Acidobacteria) and Verrucomicrobia subdivision 1 groups in bulk and leek (Allium porrum) rhizosphere soils. J. Microbiol. Methods 2010, 83, 141–148. [Google Scholar] [CrossRef] [PubMed]
  50. Figueroa-Gonzalez, P.A.; Bornemann, T.L.; Adam, P.S.; Plewka, J.; Révész, F.; von Hagen, C.A.; Probst, A.J. Saccharibacteria as organic carbon sinks in hydrocarbon-fueled communities. Front. Microbiol. 2020, 11, 3343. [Google Scholar] [CrossRef] [PubMed]
  51. Nuccio, E.E.; Anderson Furgeson, J.; Estera, K.Y.; Pett Ridge, J.; de Valpine, P.; Brodie, E.L.; Firestone, M.K. Climate and edaphic controllers influence rhizosphere community assembly for a wild annual grass. Ecology 2016, 97, 1307–1318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Wiegel, J.; Tanner, R.; Rainey, F.A. An introduction to the family Clostridiaceae. Prokaryotes 2006, 4, 654–678. [Google Scholar]
  53. Pudełko, R.; Angiel, P.J.; Potocki, M.; Jędrejek, A.; Kozak, M. Fluctuation of glacial retreat rates in the eastern part of Warszawa Icefield, King George Island, Antarctica, 1979–2018. Remote Sens. 2018, 10, 892. [Google Scholar] [CrossRef] [Green Version]
  54. Zdanowski, M.K.; Bogdanowicz, A.; Gawor, J.; Gromadka, R.; Wolicka, D.; Grzesiak, J. Enrichment of cryoconite hole anaerobes: Implications for the subglacial microbiome. Microb. Ecol. 2017, 73, 532–538. [Google Scholar] [CrossRef] [PubMed]
  55. Barret, M.; Morrissey, J.P.; O’gara, F. Functional genomics analysis of plant growth-promoting rhizobacterial traits involved in rhizosphere competence. Biol. Fertil. Soils 2011, 47, 729. [Google Scholar] [CrossRef]
  56. Berg, G.; Smalla, K. Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 2009, 68, 1–13. [Google Scholar] [CrossRef] [Green Version]
  57. Compant, S.; Clément, C.; Sessitsch, A. Plant growth-promoting bacteria in the rhizo-and endosphere of plants: Their role, colonization, mechanisms involved and prospects for utilization. Soil Biol. Biochem. 2010, 42, 669–678. [Google Scholar] [CrossRef] [Green Version]
  58. Fournier, J.; Teillet, A.; Chabaud, M.; Ivanov, S.; Genre, A.; Limpens, E.; Barker, D.G. Remodeling of the infection chamber before infection thread formation reveals a two-step mechanism for rhizobial entry into the host legume root hair. Plant Physiol. 2015, 167, 1233–1242. [Google Scholar] [CrossRef]
  59. Tian, X.Y.; Zhang, C.S. Illumina-based analysis of endophytic and rhizosphere bacterial diversity of the coastal halophyte Messerschmidia sibirica. Front. Microbiol. 2017, 8, 2288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Moronta-Barrios, F.; Gionechetti, F.; Pallavicini, A.; Marys, E.; Venturi, V. Bacterial microbiota of rice roots: 16S-based taxonomic profiling of endophytic and rhizospheric diversity, endophytes isolation and simplified endophytic community. Microorganisms 2018, 6, 14. [Google Scholar] [CrossRef] [Green Version]
  61. Dong, C.J.; Wang, L.L.; Li, Q.; Shang, Q.M. Bacterial communities in the rhizosphere, phyllosphere and endosphere of tomato plants. PLoS ONE 2019, 14, e0223847. [Google Scholar] [CrossRef] [PubMed]
  62. Cocking, E.C. Endophytic colonization of plant roots by nitrogen-fixing bacteria. Plant Soil 2003, 252, 169–175. [Google Scholar] [CrossRef]
  63. Bahulikar, R.A.; Chaluvadi, S.R.; Torres-Jerez, I.; Mosali, J.; Bennetzen, J.L.; Udvardi, M. Nitrogen fertilization reduces nitrogen fixation activity of diverse diazotrophs in switchgrass roots. Phytobiomes J. 2020. [Google Scholar] [CrossRef] [Green Version]
  64. Jiang, Z.K.; Tuo, L.; Huang, D.L.; Osterman, I.A.; Tyurin, A.P.; Liu, S.W.; Li, F.N. Diversity, novelty, and antimicrobial activity of endophytic actinobacteria from mangrove plants in Beilun Estuary National Nature Reserve of Guangxi, China. Front. Microbiol. 2018, 9, 868. [Google Scholar] [CrossRef] [PubMed]
  65. Matsumoto, A.; Takahashi, Y. Endophytic actinomycetes: Promising source of novel bioactive compounds. J. Antibiot. 2017, 70, 514–519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Zhang, Q.; Acuña, J.J.; Inostroza, N.G.; Duran, P.; Mora, M.L.; Sadowsky, M.J.; Jorquera, M.A. Niche differentiation in the composition, predicted function, and co-occurrence networks in bacterial communities associated with Antarctic vascular plants. Front. Microbiol. 2020, 11, 1036. [Google Scholar] [CrossRef] [PubMed]
  67. Ofek, M.; Hadar, Y.; Minz, D. Ecology of root colonizing Massilia (Oxalobacteraceae). PLoS ONE 2012, 7, e40117. [Google Scholar] [CrossRef]
  68. Zdanowski, M.K.; Żmuda-Baranowska, M.J.; Borsuk, P.; Świątecki, A.; Górniak, D.; Wolicka, D.; Grzesiak, J. Culturable bacteria community development in postglacial soils of Ecology Glacier, King George Island, Antarctica. Polar Biol. 2013, 36, 511–527. [Google Scholar] [CrossRef] [Green Version]
  69. Gustavs, L.; Görs, M.; Karsten, U. Polyol patterns in biofilm forming aeroterrestrial green algae (Trebouxiophyceae, Chlorophyta) 1. J. Phycol. 2011, 47, 533–537. [Google Scholar] [CrossRef]
  70. Kolton, M.; Erlacher, A.; Berg, G.; Cytryn, E. The Flavobacterium genus in the plant holobiont: Ecological, physiological, and applicative insights. In Microbial Models: From Environmental to Industrial Sustainability; Springer: Singapore, 2016; pp. 189–207. [Google Scholar]
  71. Molina-Montenegro, M.A.; Ballesteros, G.I.; Castro-Nallar, E.; Meneses, C.; Gallardo-Cerda, J.; Torres-Díaz, C. A first insight into the structure and function of rhizosphere microbiota in Antarctic plants using shotgun metagenomic. Polar Biol. 2019, 42, 1825–1835. [Google Scholar] [CrossRef]
  72. Maimaiti, J.; Zhang, Y.; Yang, J.; Cen, Y.P.; Layzell, D.B.; Peoples, M.; Dong, Z. Isolation and characterization of hydrogen oxidizing bacteria induced following exposure of soil to hydrogen gas and their impact on plant growth. Environ. Microbiol. 2007, 9, 435–444. [Google Scholar] [CrossRef] [PubMed]
  73. Coats, V.C.; Rumpho, M.E. The rhizosphere microbiota of plant invaders: An overview of recent advances in the microbiomics of invasive plants. Front. Microbiol. 2014, 5, 368. [Google Scholar] [CrossRef] [PubMed]
  74. Galera, H.; Chwedorzewska, K.J.; Wódkiewicz, M. Response of Poa annua to extreme conditions: Comparison of morphological traits between populations from cold and temperate climate conditions. Polar Biol. 2015, 38, 1657–1666. [Google Scholar] [CrossRef]
  75. Vega-Celedón, P.; Bravo, G.; Velásquez, A.; Cid, F.P.; Valenzuela, M.; Ramírez, I.; Seeger, M. Microbial diversity of psychrotolerant bacteria isolated from wild flora of andes mountains and patagonia of chile towards the selection of plant growth-promoting bacterial consortia to alleviate cold stress in plants. Microorganisms 2021, 9, 538. [Google Scholar] [CrossRef] [PubMed]
  76. Tiryaki, D.; Aydın, İ.; Atıcı, Ö. Psychrotolerant bacteria isolated from the leaf apoplast of cold-adapted wild plants improve the cold resistance of bean (Phaseolus vulgaris L.) under low temperature. Cryobiology 2019, 86, 111–119. [Google Scholar] [CrossRef] [PubMed]
  77. Araya, M.A.; Valenzuela, T.; Inostroza, N.G.; Maruyama, F.; Jorquera, M.A.; Acuña, J.J. Isolation and characterization of cold-tolerant hyper-ACC-degrading bacteria from the rhizosphere, endosphere, and phyllosphere of antarctic vascular plants. Microorganisms 2020, 8, 1788. [Google Scholar] [CrossRef] [PubMed]
  78. Tian, Y.; Gao, L. Bacterial diversity in the rhizosphere of cucumbers grown in soils covering a wide range of cucumber cropping histories and environmental conditions. Microb. Ecol. 2014, 68, 794–806. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 09 00811 g001 550
Figure 1. Operational taxonomic unit (OTU) numbers (upper graphs) and positive EcoPlate response (PER) numbers (lower graphs) for the bacterial communities associated with Poa annua roots. OTU numbers represent a value derived from three pooled samples; positive EcoPlate responses are the mean value of 9 replicates.
Figure 1. Operational taxonomic unit (OTU) numbers (upper graphs) and positive EcoPlate response (PER) numbers (lower graphs) for the bacterial communities associated with Poa annua roots. OTU numbers represent a value derived from three pooled samples; positive EcoPlate responses are the mean value of 9 replicates.
Microorganisms 09 00811 g001
Microorganisms 09 00811 g002 550
Figure 2. Relative abundance by percentile contribution of sequences identified on a phylum-rank taxonomic level. Other—group of bacterial phyla below 1% relative abundance each, S—rhizospheric soil samples, R—root samples, P1–P2—Central European (Poland) samples, P3–P5—Antarctic samples (King George Island).
Figure 2. Relative abundance by percentile contribution of sequences identified on a phylum-rank taxonomic level. Other—group of bacterial phyla below 1% relative abundance each, S—rhizospheric soil samples, R—root samples, P1–P2—Central European (Poland) samples, P3–P5—Antarctic samples (King George Island).
Microorganisms 09 00811 g002
Microorganisms 09 00811 g003 550
Figure 3. Relative abundance heatmap of sequences identified on a family-rank taxonomic level. (A)—Relative abundance according to sequence percentage value; color scale—relative abundance (%), (B)—Scaling within family rows across all examined samples; relative abundance above average—red, relative abundance below average—blue. S—rhizospheric soil samples, R—root samples, P1–P2—Central European (Poland) samples, P3–P5—Antarctic samples (King George Island).
Figure 3. Relative abundance heatmap of sequences identified on a family-rank taxonomic level. (A)—Relative abundance according to sequence percentage value; color scale—relative abundance (%), (B)—Scaling within family rows across all examined samples; relative abundance above average—red, relative abundance below average—blue. S—rhizospheric soil samples, R—root samples, P1–P2—Central European (Poland) samples, P3–P5—Antarctic samples (King George Island).
Microorganisms 09 00811 g003
Microorganisms 09 00811 g004 550
Figure 4. Heatmap displaying Poa annua rhizosphere and root community responses on Biolog Ecoplates. (A)—mean absorbance values from three replicates; color scale—absorbance values at 590 nm, (B)—scaling within a particular carbon source across all examined samples; absorbance above average—red, absorbance below average—blue. S—rhizospheric soil samples, R—root samples, P1–P2—Central European (Poland) samples, P3–P5—Antarctic samples (King George Island).
Figure 4. Heatmap displaying Poa annua rhizosphere and root community responses on Biolog Ecoplates. (A)—mean absorbance values from three replicates; color scale—absorbance values at 590 nm, (B)—scaling within a particular carbon source across all examined samples; absorbance above average—red, absorbance below average—blue. S—rhizospheric soil samples, R—root samples, P1–P2—Central European (Poland) samples, P3–P5—Antarctic samples (King George Island).
Microorganisms 09 00811 g004
Microorganisms 09 00811 g005 550
Figure 5. (A)—Correlogram of rhizospheric family-rank sequence abundance data and soil chemistry. Only significant (p < 0.05) correlations are shown. (B)—Correlogram of root endosphere family-rank sequence abundance data. Only significant (p < 0.05) correlations are shown; the color scale represents correlation coefficient values: dark green—positive correlation, dark red—negative correlation.
Figure 5. (A)—Correlogram of rhizospheric family-rank sequence abundance data and soil chemistry. Only significant (p < 0.05) correlations are shown. (B)—Correlogram of root endosphere family-rank sequence abundance data. Only significant (p < 0.05) correlations are shown; the color scale represents correlation coefficient values: dark green—positive correlation, dark red—negative correlation.
Microorganisms 09 00811 g005
Microorganisms 09 00811 g006 550
Figure 6. Principal component analysis (PCA) of biological and chemical data. (A)—PCA based on percentage contribution of bacterial sequences identified on a family-rank level. (B)—PCA based on normalized responses obtained for bacterial communities by the Biolog Ecoplate method. (C)—PCA based on soil chemical data. (D)—PCA based on a combination of family-rank bacterial sequence percentile contribution and normalized community responses on Biolog Ecoplates.
Figure 6. Principal component analysis (PCA) of biological and chemical data. (A)—PCA based on percentage contribution of bacterial sequences identified on a family-rank level. (B)—PCA based on normalized responses obtained for bacterial communities by the Biolog Ecoplate method. (C)—PCA based on soil chemical data. (D)—PCA based on a combination of family-rank bacterial sequence percentile contribution and normalized community responses on Biolog Ecoplates.
Microorganisms 09 00811 g006
Microorganisms 09 00811 g007 550
Figure 7. Statistically significant differences (p < 0.05) between different groups (see legends) of Poa annua root-associated communities: (A)—sequence contribution identified on a phylum taxonomic level; (B)—sequence contribution identified on a family-rank taxonomic level; (C)—community responses on Biolog Ecoplates based on absorbance values at 590 nm (A_590).
Figure 7. Statistically significant differences (p < 0.05) between different groups (see legends) of Poa annua root-associated communities: (A)—sequence contribution identified on a phylum taxonomic level; (B)—sequence contribution identified on a family-rank taxonomic level; (C)—community responses on Biolog Ecoplates based on absorbance values at 590 nm (A_590).
Microorganisms 09 00811 g007
Microorganisms 09 00811 g008 550
Figure 8. Boxplots displaying the range of relative sequence abundance of major Poa annua root-dwelling bacterial families. * significantly enriched in the endosphere compared to the rhizosphere.
Figure 8. Boxplots displaying the range of relative sequence abundance of major Poa annua root-dwelling bacterial families. * significantly enriched in the endosphere compared to the rhizosphere.
Microorganisms 09 00811 g008
Table
Table 1. Sampling site characteristics; m. a. s. l.—meters above sea level.
Table 1. Sampling site characteristics; m. a. s. l.—meters above sea level.
Sampling SiteGeographical CoordinatesDistance to the SeaAltitudeStructure of VegetationLandform and Habitat
P152°06′22″ N
21°05′56″ E
Powsin, Poland		285 km90 m.a.s.l.Poa annua (90%) and Poa pratensis (10%).Rubble path, sand, and sandstone.
Dry site, sheltered away from sources of nutrients.
P252°06′29″ N
21°05′33″ E
Powsin, Poland		285 km95 m.a.s.l.Poa annua (30%), Malus domestica (10%), Trifolium arvense (30%), Sonchus oleraceus L. (10%), Sonchus arvensis (20%).Soil, mechanically altered by human activities.
Wet site, fertile, extra fertilization once a year with an animal post-food mass.
P362°09′35″ S
58°28′26″ W
Arctowski Station,
Antarctica		120 m0.5 m.a.s.l.,Colobanthus quitensis (15%), Deschampsia antarctica (25%), Poa annua (10%), mosses (25%) and fruticose, foliose and crustose lichens (25%).Base station, soil mechanically altered by human activities; Skeletic Eutric Fluvisol (Turbic)
The site is strongly influenced by marine aerosols, moist, with large human influence.
P462°10′05″ S
58°27′46″ W
Ecology Glacier foreland, Antarctica		20 m0.5 m.a.s.l.,Mosses (40%), fruticose and foliose lichens (40%), Colobanthus quitensis (10%) and Deschampsia antarctica (10%).Fluted moraine; Eutric Skeletic Protic Regosol (Turbic).
Dry site, sheltered away from sources of nutrients with a little influence of marine aerosols.
P562°09′33″ S
58°28′25″ W
Arctowski Station, Antarctica		100 m0.5 m.a.s.l.,Colobanthus quitensis (15%), Deschampsia antarctica (25%), Poa annua (10%), mosses (25%) and fruticose, foliose and crustose lichens (25%).Base station, soil, mechanically altered by human activities; Skeletic Eutric Fluvisol (Turbic)
The site is medium influenced by marine aerosols, moist, with large human influence, especially by using big vehicles.
Table
Table 2. Soil component concentrations.
Table 2. Soil component concentrations.
Soil ComponentsS1S2S3S4S5
N (mgNO3/100 g)17.427.58.816.37.7
P (mgP2O5/100 g)6.0219.25.615.14.9
K (mg K2O/100 g)36.1252.831.155.831.9
Mg (mg Mg/100 g)38.4185.4244.7205250.7
Ca (mg Ca/100 g)>50002024.6498.01020.1533.3
Na (mg Na/100 g)4.94.0>150>150>150
Salinity (g NaCl/L)0.20.30.40.40.4
pH7.87.07.47.57.4
Mn (mg Mn/kg)52.8218.1147.5132.1155.5
Zn (mg Zn/kg)5.9>50>509.2>50
Cu (mg Cu/kg)2.75.721.8>3029.8
Fe (mg Fe/kg)389.41502.54412.64670.24269.8
N—nitrate nitrogen, P—labile phosphorus, K—labile potassium, S1–S2—European samples, S3–S5—Antarctic samples.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Znój, A.; Grzesiak, J.; Gawor, J.; Gromadka, R.; Chwedorzewska, K.J. Bacterial Communities Associated with Poa annua Roots in Central European (Poland) and Antarctic Settings (King George Island). Microorganisms 2021, 9, 811. https://doi.org/10.3390/microorganisms9040811

AMA Style

Znój A, Grzesiak J, Gawor J, Gromadka R, Chwedorzewska KJ. Bacterial Communities Associated with Poa annua Roots in Central European (Poland) and Antarctic Settings (King George Island). Microorganisms. 2021; 9(4):811. https://doi.org/10.3390/microorganisms9040811

Chicago/Turabian Style

Znój, Anna, Jakub Grzesiak, Jan Gawor, Robert Gromadka, and Katarzyna J. Chwedorzewska. 2021. "Bacterial Communities Associated with Poa annua Roots in Central European (Poland) and Antarctic Settings (King George Island)" Microorganisms 9, no. 4: 811. https://doi.org/10.3390/microorganisms9040811

APA Style

Znój, A., Grzesiak, J., Gawor, J., Gromadka, R., & Chwedorzewska, K. J. (2021). Bacterial Communities Associated with Poa annua Roots in Central European (Poland) and Antarctic Settings (King George Island). Microorganisms, 9(4), 811. https://doi.org/10.3390/microorganisms9040811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop