Next Article in Journal
Grassland Degradation Has Stronger Effects on Soil Fungal Community Than Bacterial Community across the Semi-Arid Region of Northern China
Previous Article in Journal
Association between the Classification of the Genus of Batrachospermaceae (Rhodophyta) and the Environmental Factors Based on Machine Learning
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 24
10.3390/plants11243487
plants-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

Genomic Characterization of Aureimonas altamirensis C2P003—A Specific Member of the Microbiome of Fraxinus excelsior Trees Tolerant to Ash Dieback

by
Regina Becker
1,
Kristina Ulrich
2,
Undine Behrendt
1,
Volker Schneck
2 and
Andreas Ulrich
1,*
1
Microbial Biogeochemistry, Research Area Landscape Functioning, Leibniz Centre for Agricultural Landscape Research (ZALF), 15374 Müncheberg, Germany
2
Institute of Forest Genetics, Johann Heinrich von Thünen Institute, 15377 Waldsieversdorf, Germany
*
Author to whom correspondence should be addressed.
Plants 2022, 11(24), 3487; https://doi.org/10.3390/plants11243487
Submission received: 15 November 2022 / Revised: 7 December 2022 / Accepted: 9 December 2022 / Published: 13 December 2022
(This article belongs to the Section Plant Protection and Biotic Interactions)
Browse Figures
Versions Notes

Abstract

:
Some European ash trees show tolerance towards dieback caused by the invasive pathogen Hymenoscyphus fraxineus. The microbiome of these trees harbours a range of specific bacterial groups. One of these groups belonging to the species Aureimonas altamirensis was studied in detail by genome analysis and a plant inoculation trial. The strain group was shown to be phylogenetically distinct from clinical isolates by 16S rRNA analysis and phylogenomics. Genome analysis of a representative strain C2P003 resulted in a large number of unique gene sequences in comparison to other well-studied strains of the species. A functional analysis of the genome revealed features associated with the synthesis of exopolysaccharides, protein secretion and biofilm production as well as genes for stress adaptation, suggesting the ability of C2P003 to effectively colonize ash leaves. The inoculation of ash seedlings with C2P003 showed a significant positive effect on the plant health of the seedlings that were exposed to H. fraxineus infection. This effect was maintained over a period of three years and was accompanied by a significant shift in the bacterial microbiome composition one year after inoculation. Overall, the results indicate that C2P003 may suppress H. fraxineus in or on ash leaves via colonization resistance or indirectly by affecting the microbiome.

1. Introduction

The plant microbiome is an integral part of its host and comprises a tremendous amount of beneficial, commensal and potentially pathogenic microorganisms. Both plant-microbe interactions and competition among microorganisms affect the plant performance [1,2]. Plants subjected to diverse biotic and abiotic stresses can induce changes in root and leaf exudates and subsequently alter the plant microbiome assembly [3]. A special role in the microbiome is played by plant beneficial bacteria that effectively colonize the plant and contribute to nutrient supply [4,5] and/or defence against pathogens by direct antagonism, competition or plant-supporting metabolic activities [6,7,8]. The abundance of these microorganisms in the microbiome has been linked to plant growth and health in several studies [9,10,11]. The spread of pathogens may have dramatic effects on ecosystem functioning, particularly in forests. This makes it all the more important to successfully control the pathogens of the trees through the application of biocontrol agents [12].
In the last few years, the spread of ash dieback has resulted in a substantial threat to Fraxinus excelsior stands in central and northern Europe [13,14]. Infection with the causal agent Hymenoscyphus fraxineus begins with the germination of wind-borne ascospores on the leaves, after which the fungus penetrates into the intercellular space and then into shoots and twigs, where it causes necrotic bark lesions and blockages of the xylem vessels [15]. The severity of infestations with H. fraxineus might be reduced by resident phyllosphere or endosphere colonizers, which are able to counteract the entry and penetration of the invader [16,17]. Comparative analyses of the bacterial leaf microbiome of F. excelsior in ash dieback-affected forests revealed a significantly increased abundance of individual bacterial groups in tolerant trees belonging to the genera Aureimonas, Luteimonas, Pseudomonas, Bacillus and Paenibacillus [18]. With the exception of Aureimonas, all these genera are known for strains with plant-beneficial traits. Thus, a health-stabilizing effect on ash seedlings was recently demonstrated for a Luteimonas fraxinea strain obtained from the leaves of healthy ashes [19]. Other members of Luteimonas also had stimulating effects based on disease suppression and nitrogen acquisition [20,21,22]. The genus Pseudomonas comprises many strains contributing to plant health by antibiosis, competition and induction of systemic resistance [23,24,25,26]. Finally, strains of the genera Bacillus and Paenibacillus are of exceptional importance as biocontrol agents, especially due to the synthesis of bioactive secondary metabolites [27,28,29].
Unlike these genera, members of Aureimonas have not yet been characterized as plant-beneficial bacteria. At present, the genus contains 18 species originally derived from a wide range of habitats. Aureimonas altamirensis was first isolated from the subterranean environment of the Altamira cave as part of a complex microbial community [30]. In the following years, several publications reported on the pathogenic potential of A. altamirensis strains based on evidence in human clinical samples [31,32,33,34]. However, the clinical relevance of the species remains unclear since infestations are discussed in relation to immune deficiencies and the simultaneous detection of other organisms [32,35,36]. A. altamirensis is therefore considered an opportunistic pathogen [36,37], and further studies are needed to elucidate its epidemiology, pathogenic mechanisms, and clinical significance [38]. Moreover, the species was also repeatedly detected in environmental habitats, e.g., as a cellulolytic strain obtained from vegetable biomass [39].
We discovered A. altamirensis strains as specific colonizers of leaves from dieback tolerant ash trees [18] and hypothesized that the species could play a role in plant resilience. Here we characterize this group of isolates in detail. The phylogenetic analysis of a representative isolate was performed to determine the position within the genus and the species. The genome was further investigated to reveal features that could be related to plant colonization and plant-protective properties. Finally, an inoculation trial on ash seedlings exposed to H. fraxineus was performed to clarify the effects of the studied isolate on plant health and the microbiome.

2. Results and Discussion

Analyses of the leaf microbiome of 80-years-old ash trees under H. fraxineus exposure revealed several bacterial groups with significantly higher abundance in tolerant trees [18]. One of these groups was assigned to A. altamirensis, which was found together with well-known plant colonizers, such as Luteimonas and Pseudomonas. The abundance of the A. altamirensis strain group was increased 34 and 49-fold in tolerant trees as analysed by 16S rRNA gene amplicon sequencing and culturing, respectively. A total of 62 isolates formed this single group based on MALDI-TOF MS (matrix-assisted laser desorption/ionization time-of-flight mass spectrometry), which were isolated almost exclusively from tolerant ash trees.

2.1. Phylogenetic Analysis

To gain insight into the variability of the ash-associated A. altamirensis group, peak profiles derived from MALDI-TOF MS spectra of the 62 isolates were analysed by principal component analysis. The ordination plot indicated clear differences between the isolates (Figure 1a). The isolate C4P030 showed a strong distance to nearly all other group members. In addition, several different clusters could be identified within the A. altamirensis strain group. Seven isolates with clear distances to each other were selected for sequencing of the ribosomal intergenic spacer (IGS) region. For phylogenetic analysis, the IGS sequences of the isolates were compared with those of the type strain DSM 21988 and the clinically derived strains ON-56566 and OT7. The tree shows a separate cluster of all the ash isolates. Within this branch, the isolates obtained from different trees (C2 and C4) could be distinguished. A strong differentiation was found to DSM 21988 and with even higher distances to the clinical strains ON-56566 and OT7 (Figure 1b). The sequencing of the 16S rRNA gene revealed identical sequences for all seven isolates. The isolate C2P003 was subsequently chosen as a representative isolate and subjected to a comprehensive analysis.
A phylogenetic analysis based on the 16S rRNA gene showed a close relationship of isolate C2P003 to the type strain DSM 21988 (isolated from the subterranean environment of a cave), the strain gall3186, which is a member of a gall wasp microflora [40], and SBP73, a cellulolytic strain obtained from vegetable biomass [39] (Figure 2).
The four strains of environmental origin formed a shared cluster that was distinct from the isolates detected in clinical samples [31,32,37,38,41]. Within the genus, A. frigidaquae and “A. fodinaquatilis” were found to be the closest related species. However, the phylogenetic positions of the Aureimonas species within the genus as suggested by the 16S rRNA gene analysis were only slightly supported by the bootstrap values. Consequently, the analysis was complemented by genome phylogeny using the genome sequences of the currently known Aureimonas species. The comparison between C2P003 and the type strain of A. altamirensis DSM 21988 resulted in an average nucleotide identity (ANI) of 96.4% and a digital DNA-DNA hybridization value (dDDH) of 70.0%. Both values were above the respective thresholds (≥95% and ≥70%) and confirmed the affiliation of C2P003 to the species A. altamirensis. Comparable to the 16S rRNA analysis, the phylogenomic tree showed the close relationship between C2P003 and A. altamirensis DSM 21988 with a clear separation to A. altamirensis ON-56566 and Aureimonas sp. OT7, supported by high bootstrap values (Figure S1). Strain OT7, isolated from human skin, and strain ON-56566, isolated from patient blood culture, were reported to belong to the same species clade that is separated from the A. altamirensis type strain [41]. Consistent with the 16S rRNA data, a phylogenomic analysis indicated that A. frigidaquae and “A. fodinaquatilis” were the species most closely related to A. altamirensis.

2.2. Genome Analysis

The sequencing and assembly of genomic DNA from the A. altamirensis strain C2P003 resulted in a complete circular genome. However, the annotation of the genome sequence revealed several clusters with plasmid genes (trbIGFLJEDCB, repB) (Table S1). It can therefore not be ruled out that in the current assembly, the chromosome was merged with one or several plasmids. Accordingly, the latest assembly (GCF_021228915.1) should be considered as a draft genome. The genome sequence comprised 4,592,981 bp with a total of 4289 predicted protein-coding genes (Table 1). The genome size, G+C content and the other annotation statistics are comparable to those of the A. altamirensis type strain DSM 21988 and the clinically derived strains ON-56566 and OT7. However, strain C2P003 differs from its closest relatives by a larger genome size and a higher number of protein coding sequences.
The genomes of the four A. altamirensis strains share 3253 orthologous genes (Figure 3), which account for 69.6% of all predicted protein-coding sequences of strain C2P003. A total of 3725 genes of C2P003 (79.7%) had orthologues in A. altamirensis DSM 21988, indicating the close relationship between the two strains. By comparison, lower proportions of genes (73.28 and 74.4%) were shared with ON-56566 and OT7, respectively. In addition, the analysis of orthologous genes again revealed the specificity of C2P003 due to the highest quantity of unique gene sequences (741 vs. 266–398).
Genome annotation did not reveal any genes whose products are clearly associated with pathogenic functions. In addition, a comparison with pathogenic Brucella species that possess a set of virulence factors allowing bacteria to replicate in host cells and to induce persistent infections was performed. A key role in the virulence of Brucella is played by the T4 secretion system (T4SS), which is used to transfer toxic effector proteins into the host cell [42,43]. On the other hand, T4 secretion systems are equally common in nonhuman-pathogenic species, such as symbiotic rhizobia [44]. The genes for T4SS transfer proteins could also be detected in the genomes of the Aureimonas strains C2P003, DSM 21988 and ON-56566 (Table 2 and Table S1). Further virulence factors of Brucella, in particular genes involved in the biosynthesis of the lipopolysaccharide O-side chain (e.g., wboA, wboB, wbkA, wbkE, gmd) [45], could not be identified in the genome of any Aureimonas strains. The phenotypic tests with strain C2P003 did not indicate a pathogenic potential. The strain was not able to grow at human body temperature (37 °C), did not show extracellular DNAse activity and was haemolytic-negative. Overall, there were no indications for pathogenicity.

2.3. Genomic Features

To study the genetic bases for a plant-associated lifestyle and potential plant-enhancing properties of A. altamirensis C2P003, a functional analysis of the annotated genome was performed. The analysis, conducted in comparison to the type strain DSM 21988 and the clinical strain ON-56566, was focused on the detection of gene mediating functions in plant colonization, defence against competitive microorganisms and resistance to abiotic stress.

2.3.1. Host Colonization

Bacteria gain access to plant tissues and cells by motility, chemotaxis, biofilms and the ability to degrade plant polymers [46,47,48]. Numerous genes conferring these traits have been identified in the genomes of plant-colonizing bacteria [49,50,51]. Cells of C2P003 as well as A. altamirensis DSM 21988 [30] were found to be nonmotile and genome analysis of the three A. altamirensis strains did not reveal genes conferring motility and chemotaxis. However, strains C2P003, DSM 21988 and ON-56566 possess several genes related to the synthesis of exopolysaccharides and protein secretion system types I, IV and VII (Table 2, Supplementary Table S1). The secretion of both polysaccharides and proteins promotes the attachment of bacteria to the host cells and thus the formation of biofilms [52,53]. Protein secretion systems, in particular, are involved in numerous biotic interactions, including mutualistic plant colonization but also the infection induction by pathogens [54,55]. In contrast to strains C2P003 and DSM 21988, the genome of ON-56566 lacks genes for secretion systems type I and VII.
Moreover, all analysed A. altamirensis genomes harbour genes coding for plant cell wall-degrading enzymes that are crucial for bacteria to be able to access and proliferate in the intercellular space [46,56]. The detected genes encode for several enzymes involved in the degradation of cellulose, hemicellulose and pectin (Table 3, Supplementary Table S2).

2.3.2. Antimicrobial Substances

Some further genes in the genomes of C2P003, DSM 21988 and ON-56566 encode for enzymes required in the degradation of bacterial and fungal cell walls, such as β-N-acetylglucosaminidase, N-acetylmuramoyl-L-alanine amidase and phage lysozyme R (Table 3, Supplementary Table S2). β-N-acetylglucosaminidase plays an important role in the degradation of chitin, the bulk component of fungal cell walls [57], as well as in the degradation of peptidoglycan, the main component of bacterial cell walls [58,59]. N-acetylmuramoyl-L-alanine amidase is the second most important enzyme required in the lysis of peptidoglucan [58,60], whereas lysozyme contributes to the breakdown of peptidoglucan and partly also of chitin [61,62]. In addition, the three A. altamirensis strains are equipped with antibiotic resistance genes (Table 2). Both the ability to attack microbial competitors and the self-protection from invaders may support the establishment of C2P003 in the host.

2.3.3. Stress Adaptation

Colonizers of the phyllosphere, the interface between the atmosphere and plants, are exposed to strong and unpredictable fluctuations in temperature, moisture and UV levels. Adaptations to osmotic and oxidative stress are therefore basic abilities for survival in this habitat [63].
The most common response to osmotic stress is the accumulation of potassium or osmolytes, such as amino acids, betaines, polyols or sugars [64]. The genome of C2P003 contains genes coding for betaine aldehyde dehydrogenase and choline dehydrogenase, which are involved in the biosynthesis of the osmoprotectant glycine betaine from choline [65,66] and the production of glycine betaine [67], respectively. Furthermore, the presence of the conserved yehZYXW operon encoding the stress-regulated ABC membrane transport system [68] might contribute to the resilience of C2P003 against osmotic stress. Genes related to the osmotic stress response were also found in the other A. altamirensis strains (Table 2, Supplementary Table S1).
To counteract oxidative stress and to restore a redox balance, aerobic organisms possess complex mechanisms regulating the activation or silencing of genes for defensive enzymes, transcription factors and structural proteins [69,70]. Genes coding for antioxidant enzymes, i.e., superoxide dismutase [71] and glutathione S-transferase [72,73], were found in the genomes of C2P003, DSM 21988 and ON-56566 in equal numbers (Table 2, Supplementary Table S1).
In addition to harmful atmospheric effects, heavy metals absorbed by the plants can also damage plant-associated bacteria. In response to this challenge, bacteria have evolved a range of mechanisms for heavy metal tolerance including transport, sequestration and reduction of metal ions [74,75,76]. This enables plant-associated metal-tolerant bacteria to reduce the metal contamination of plants and to enhance plant growth [77,78]. Genes contributing to copper homeostasis, cobalt-zinc-cadmium resistance and resistance to chromium compounds were detected in the genomes of C2P003, DSM 21988 and ON-56566 (Table 2; Supplementary Table S1).
In conclusion, the analysis of the genomic features of A. altamirensis C2P003 revealed numerous genes involved in metabolic functions for survival and competitiveness in the phyllosphere and the plant endosphere.

2.4. Inoculation of Ash Seedlings

The effect of A. altamirensis C2P003 on ash dieback was studied together with two further strains, which were also isolated from ash plants in 2017. Bacillus velezensis A4P130 was characterized as a typical antagonistic strain with a high ability to inhibit the growth of the pathogen H. fraxineus in cocultivation assays [18]. A. altamirensis C2P003 as well as L. fraxinea D4P002 showed no antagonistic activity against H. fraxineus [18]. However, both strains D4P002 and C2P003 were found to have significantly higher abundance (>20-fold) in tolerant ash trees and were conclusively specific for the microbiome of F. excelsior plants tolerant to ash dieback.
The three strains were applied for an inoculation test with ash seedlings. To simulate natural conditions of H. fraxineus infection, the seedlings were exposed to the pathogen via application of infected leaf petioles [79] four weeks after inoculation with the bacterial strains. This resulted in the first dieback symptoms in mid-August, such as leaves with dark patches and withered leaves [15]. The effect of the inoculation was monitored over a period of three years. In September of the first year (141 days after inoculation), all three strains showed a significant positive effect on plant health (Figure 4 and Figure S2). The strongest impact was assessed for C2P003. In the following year, new dieback infection symptoms were rarely observed. The majority of the control plants, which were still weakened by the primary infection, showed growth depressions with nonspecific leaf damage. Despite the lack of infection, there was again a distinct effect of the inoculation; however, in contrast to A4P130 and D4P002, the impact of C2P003 was significant at only one of the two time points of examination. In the third year, a strong fructification of H. fraxineus was again visible on the freshly applied petioles, suggesting a high likelihood of infection. Under these conditions, inoculation with the bacterial strains maintained their advantage over the controls through stronger growth and healthier leaves until the end of the vegetation period (Figure 4). The effects of inoculation with A4130 and D4P002 were significant across all three years of the plant trial [19]. Inoculation with the strain C2P003 also resulted in a distinct effect on the plant health of the ash seedlings, but at two time points this effect was not significant (Figure 4).
In the second year (15 months after inoculation), ash leaves were sampled to compare the microbiomes of the three inoculated and the control variants. Amplicon sequencing revealed a complex bacterial community with 13 phyla dominated by Pseudomonadota (average abundance across all treatments of 68.5%), Actinomycetota (19.0%), Deinococcota (9.5%) and Bacteroidota (2.8%). An overview about the taxonomic composition of the bacterial microbiome is shown in Figure 5. At the phylum level, Pseudomonadota tended to be less abundant in the control plants; however, the difference was not significant. The genera Methylobacterium, Sphingomonas, Deinococcus and Klenkia dominated the bacterial microbiome and had comparable abundances in all variants. Alpha-diversity estimates (Shannon and InvSimpson index) showed significantly lower values for the microbiomes of trees inoculated with C2P003 or D4P002 (Shannon: 4.75 ± 0.34/4.94 ± 0.25; mean ± SD) as compared to the control (5.34 ± 0.26; p < 0.02). The variance of the bacterial microbiome composition of the ash plants was studied by a principle coordinate analysis (Figure 6). As shown in the ordination plot, the first axis clearly differentiated the plants inoculated with C2P003 and D4P002 from the control. Correspondingly, ANOSIM and PERMANOVA analyses revealed significant differences between C2P003 and D4P002 vs. the control (R = 0.47/0.42, p < 0.003; R2 = 0.27/0.24, p = 0.006). In contrast, the microbiome composition of A4P130-inoculated plants did not show significant differences from that of the control plants.
Altogether, the strains specifically found on tolerant ash trees were able to trigger a shift in the bacterial microbiome more than one year after the treatment. This finding is consistent with several recent studies describing that inoculation with individual bacterial strains resulted in changes in the plant microbiome that were associated with improved plant health [81,82,83,84,85]. Berg et al. [86] reviewed the effects of microbial inoculants on the indigenous plant microbiome and suggested microbiome modulations as a novel and efficient mode of action for microbial inoculants that can also be mediated via the plant.
In contrast to C2P003 and D4P002, the antagonistic strain B. velezensis A4P130 was not able to modify the microbiome composition, which fits well with data on the comprehensively studied biocontrol strain B. velezensis FZB42 [87]. Comparable to A4P130, the treatment of lettuce with FZB42 led to improved plant performance, but a lasting shift of the microbiome could not be observed [88]. In this case, the antagonistic capabilities of A4P130 were suggested for the health-protecting effects.

3. Materials and Methods

3.1. Isolation and Classification of Bacteria

A. altamirensis isolates were obtained during a sampling campaign of tolerant and susceptible ash trees conducted in July 2017 in a forest area of Northeast Germany with severe infestation of H. fraxineus [18]. All isolates originated from two tolerant trees without visible symptoms that were located in the district Pennin on plot C (54°15′ N, 13°01′ E) and a susceptible tree located in the district Lendershagen on plot A (54°14′ N, 12°51′ E). The strains were isolated from compound leaves using a method allowing the cultivation of both epi- and endophytic bacteria.
Classification was based on MALDI-TOF MS. Isolates that could not be unambiguously identified by the Bruker database were assigned by the sequencing of almost the complete 16S rRNA gene [18].

3.2. Sequencing of the Ribosomal Intergenic Spacer (IGS) Region

The phylogenetic variability of the selected Aureimonas strains was studied by sequencing the 16S–23S rRNA gene intergenic spacer region. Total DNA extraction followed the procedure of Ulrich et al. [89]. Amplification with the primers 1492f and 115r and subsequent Sanger sequencing were performed as described by Tokajian et al. [90].

3.3. Genome Sequencing

Cells of the A. altamirensis strain C2P003 were cultured in R2 broth for two days at 25 °C and washed two times with 0.3% NaCl. Genomic DNA was extracted by the Genomic-Tip 20 Kit (Qiagen, Hilden, Germany) as described previously [19]. The fragment size, quantity and quality of DNA were assessed on a 1% agarose gel and with a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). DNA was sequenced using the Pacific Biosciences (PacBio) RS II sequencing platform at Eurofins Genomics (Konstanz, Germany). Sequence reads were de novo assembled using the PacBio hierarchical genome assembly process (HGAP4). The assembling resulted in one contig with an average genome coverage of 124×. The genome sequence was circulated with the Circlator v. 1.5.5 [91].

3.4. Phylogenetic and Genome Analyses

The genome of A. altamirensis C2P003 was annotated using RAST server version 2.0 [92] and by the NCBI prokaryotic genome annotation pipeline [93]. The RAST platform was also used for comparisons with genomes of the A. altamirensis strains DSM 21988 and ON-56566 as well as Aureimonas sp. OT7 to avoid bias by different annotation systems. The calculation of orthologous genes was based on the predicted coding sequences with an identity of more than 70% at the amino acid level. A Venn diagram was generated using the R package VennDiagram [94].
For phylogenetic analysis, 16S rRNA gene sequences from closely related strains and species were aligned using the ClustalW algorithm with MEGA X [95], resulting in an alignment of 1386 nt. The phylogenetic tree was constructed using the maximum-likelihood algorithm based on evolutionary distances of the Tamura 3-parameter model (+G+I). In the same way, the isolate group of A. altamirensis was analysed using an alignment of IGS sequences comprising 950 nt. The maximum likelihood tree was calculated using distances of the Hasegawa-Kishino-Yano model (+G). The phylogenomic analysis based on core genome phylogeny [96] was performed as described by Ulrich et al. [50]. A total of 120 bacterial core marker genes were used to form a concatenated amino acid sequence alignment, which was used to calculate a maximum-likelihood tree (LG substitution model with F+G+I) with MEGA. The calculation of the ANI and dDDH values between C2P003 and the type strain of A. altamirensis DSM 21988 was performed as described previously [19].

3.5. Assessment of A. altamirensis C2P003 for Human Pathogenicity

The possible pathogenicity of C2P003 was tested by different approaches allowing a differentiation between pathogenic and nonpathogenic strains. First, C2P003 cells were plated on Columbia blood agar (Oxoid-Thermo Fisher Scientific, Germany) and incubated at 37 °C and 5% CO2. Growth at 37 °C is a simple indication of pathogenicity [97] and has been demonstrated for a clinical A. altamirensis strain [34]. A second approach consisted of checking the activity of extracellular DNase, which is involved in the dissemination of bacteria [98]. Third, haemolytic activity as a possible indication of pathogenicity was tested via incubation on Columbia blood agar as described by the manufacturer.

3.6. Inoculation Test on Ash Seedlings

To study the effect of antagonistic bacteria on ash dieback, ash seedlings with a height of approximately 15 cm were inoculated with the three strains A. altamirensis C2P003, B. velezensis A4P130 and L. fraxinea D4P002 at the beginning of June, 2019. For preparation of the inoculum, the bacteria were grown overnight in R2 broth at 22 °C with gentle shaking at 200 rpm. Cells were collected by centrifugation at 5000× g for 5 min, washed with 1/4 Ringer solution and finally resuspended in 20 mL 1/4 Ringer solution. The cell density of the inocula was adjusted to 1 × 108 cells mL−1. The inoculation was performed by thoroughly dipping the seedlings in the inoculum for approximately 1 min. Control plants were immersed in the same manner in 1/4 Ringer solution. Afterwards, the seedlings were potted into 3 L plastic containers with potting soil and covered with glassine bags for 7 days to maintain high humidity. After 2 weeks, the plants were transferred to the nursery under controlled irrigation and shade. Infection with H. fraxineus was initiated four weeks after inoculation as well as in June of the following two years by application of infected ash leaf petioles at the bottom of the plant containers. The seedlings of each treatment (n = 12) were monitored for their health status based on a 5-point Bonitur scale modified according to Peters et al. [80] over a period of three years. An analysis of significant differences between inoculated and control plants was performed using the Kruskal-Wallis test followed by Dunn’s test of multiple comparisons using rank sums as implemented in R version 4.1.2, library dunn.test [99]. Fifteen months after inoculation, the seedlings were sampled for microbiome analysis. Pieces of all parts of compound leaves were taken from each of the plants and stored at −20 °C until further processing.

3.7. Microbiome Analysis

Total DNA was extracted from 100 mg of the ground plant material using the DNeasy Plant Mini Kit as described by Ulrich et al. [18]. PCR amplification was performed with primers 799F and 1115R which exclude the chloroplast and mitochondrial DNA of the host plant. Library preparation and 300-bp paired-end sequencing on an Illumina MiSeq was performed at LGC Genomics (Berlin, Germany). The sequence reads were processed using the DADA2 pipeline [100]. After quality filtering and the removing of potential chimaeras, residual plastid and mitochondrial sequences from ash trees were excluded from the dataset, which resulted in 2,987,432 high-quality reads. Statistical analyses were performed using the phyloseq, vegan, metagMisc, decipher, phangorn and ggplot2 packages in R. To estimate alpha diversity, the Shannon and InvSimpson indices were calculated. Differences in the bacterial microbiome composition were visualized by applying a principal coordinate analysis (PCoA) based on weighted UniFraq distances of the ASVs. Significant differences between the bacterial communities were tested using Analysis of Similarity (anosim) and PERMANOVA (adonis).

4. Conclusions

The studied group of A. altamirensis strains which were specifically found on ash trees tolerant to H. fraxineus infection represents the first evidence of the species as a plant-associated bacterium. A phylogenetic analysis resulted in a clear distinction from previously described clinical isolates of the species. A functional genome analysis of the representative strain C2P003 revealed the presence of many genes necessary for successful colonization of the phyllosphere and/or the endosphere of plants. Although C2P003 did not possess in vitro antagonistic activity against the fungal pathogen H. fraxineus, an inoculation experiment on ash seedlings demonstrated the plant health-protecting effect of C2P003, which was accompanied by a shift in the leaf microbiome one year after inoculation. It suggests that the strain C2P003 may suppress H. fraxineus via colonization resistance or indirectly by changing the microbiome composition.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11243487/s1, Figure S1: Phylogenomic tree showing the position of A. altamirensis C2P003 among other A. altamirensis strains and reference strains of related Aureimonas species. Figure S2: Ash seedlings inoculated with B. velezensis A4P130 (A), A. altamirensis C2P003 (B) and L. fraxinea D4P002 (C) as well as control plants (D). Plant health was evaluated 141 days after inoculation. Table S1: Protein-encoding genes involved in host colonization and stress adaptation in the genome of A. altamirensis C2P003. Table S2: Protein-encoding genes involved in the degradation of plant and microbe cell walls in the genome of A. altamirensis C2P003.

Author Contributions

Conceptualization, A.U. and V.S.; formal analysis, R.B., U.B., K.U. and A.U.; investigation, R.B. and K.U.; resources, R.B., K.U. and V.S.; data curation, A.U.; writing—original draft preparation, R.B. and K.U.; writing—review and editing, R.B., K.U., U.B., V.S. and A.U.; visualization, R.B., K.U., U.B. and A.U.; project administration, A.U. and V.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding via the Waldklimafonds (WKF) funded by the German Federal Ministry of Food and Agriculture (BMEL) and Federal Ministry for the Environment, Nature Conservation, Nuclear Safety and Consumer Protection (BMUV) administrated by the Agency for Renewable Resources (FNR), grant number 2219WK22I4.

Data Availability Statement

The genome sequence of the strain C2P003 was deposited in the GenBank database under the accession no. CP089311. The paired sequence reads generated for the microbiome analysis were deposited in the Sequence Read Archive under accession number PRJNA809658.

Acknowledgments

We thank Darline Krebel, Ulrike Klauß and Sonja Wende for providing excellent technical assistance.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Song, C.; Zhu, F.; Carrión, V.J.; Cordovez, V. Beyond plant microbiome composition: Exploiting microbial functions and plant traits via integrated approaches. Front. Bioeng. Biotechnol. 2020, 8, 896. [Google Scholar] [CrossRef] [PubMed]
  2. Babalola, O.O.; Fadiji, A.E.; Enagbonma, B.J.; Alori, E.T.; Ayilara, M.S.; Ayangbenro, A.S. The nexus between plant and plant microbiome: Revelation of the networking strategies. Front. Microbiol. 2020, 11, 548037. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, H.; Brettell, L.E.; Qiu, Z.; Singh, B.K. Microbiome-mediated stress resistance in plants. Trends Plant Sci. 2020, 25, 733–743. [Google Scholar] [CrossRef] [PubMed]
  4. Glick, B.R. Plant growth-promoting bacteria: Mechanisms and applications. Scientifica 2012, 2012, 963401. [Google Scholar] [CrossRef] [Green Version]
  5. Chennappa, G.; Naik, M.K.; Udaykumar, N.; Vidya, M.; Sreenivasa, M.Y.; Amaresh, Y.S.; Mathad, P.F. Chapter 10—Plant growth promoting microbes: A future trend for environmental sustainability. In New and Future Developments in Microbial Biotechnology and Bioengineering; Singh, J.S., Ed.; Elsevier: Amsterdam, The Netherlands, 2019; pp. 163–178. [Google Scholar]
  6. Janakiev, T.; Dimkić, I.; Bojić, S.; Fira, D.; Stanković, S.; Berić, T. Bacterial communities of plum phyllosphere and characterization of indigenous antagonistic Bacillus thuringiensis R3/3 isolate. J. Appl. Microbiol. 2020, 128, 528–543. [Google Scholar] [CrossRef] [PubMed]
  7. Wiraswati, S.; Nawangsih, A.; Rusmana, I.; Wahyudi, A. Rice phyllosphere bacteria producing antifungal compounds as biological control agents of blast disease. Biodivers. J. 2020, 21, 1273–1278. [Google Scholar] [CrossRef]
  8. Muñoz, C.Y.; Zhou, L.; Yi, Y.; Kuipers, O.P. Biocontrol properties from phyllospheric bacteria isolated from Solanum lycopersicum and Lactuca sativa and genome mining of antimicrobial gene clusters. BMC Genom. 2022, 23, 152. [Google Scholar] [CrossRef]
  9. Köberl, M.; Dita, M.; Martinuz, A.; Staver, C.; Berg, G. Members of Gammaproteobacteria as indicator species of healthy banana plants on Fusarium wilt-infested fields in Central America. Sci. Rep. 2017, 7, 45318. [Google Scholar] [CrossRef] [Green Version]
  10. Harman, G.; Khadka, R.; Doni, F.; Uphoff, N. Benefits to plant health and productivity from enhancing plant microbial symbionts. Front. Plant Sci. 2021, 11, 610065. [Google Scholar] [CrossRef]
  11. Zhang, L.; Zhang, M.; Huang, S.; Li, L.; Gao, Q.; Wang, Y.; Zhang, S.; Huang, S.; Yuan, L.; Wen, Y.; et al. A highly conserved core bacterial microbiota with nitrogen-fixation capacity inhabits the xylem sap in maize plants. Nat. Commun. 2022, 13, 3361. [Google Scholar] [CrossRef]
  12. Terhonen, E.; Kovalchuk, A.; Zarsav, A.; Asiegbu, F.O. Biocontrol potential of forest tree endophytes. In Endophytes of Forest Trees: Biology and Applications; Pirttilä, A.M., Frank, A.C., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 283–318. [Google Scholar]
  13. Skovsgaard, J.P.; Wilhelm, G.J.; Thomsen, I.M.; Metzler, B.; Kirisits, T.; Havrdova, L.; Enderle, R.; Dobrowolska, D.; Cleary, M.; Clark, J. Silvicultural strategies for Fraxinus excelsior in response to dieback caused by Hymenoscyphus fraxineus. Forestry 2017, 90, 455–472. [Google Scholar] [CrossRef] [Green Version]
  14. Cleary, M.; Nguyen, D.; Marciulyniene, D.; Berlin, A.; Vasaitis, R.; Stenlid, J. Friend or foe? Biological and ecological traits of the European ash dieback pathogen Hymenoscyphus fraxineus in its native environment. Sci. Rep. 2016, 6, 21895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Gross, A.; Holdenrieder, O.; Pautasso, M.; Queloz, V.; Sieber, T.N. Hymenoscyphus pseudoalbidus, the causal agent of European ash dieback. Mol. Plant Pathol. 2014, 15, 5–21. [Google Scholar] [CrossRef] [PubMed]
  16. Cross, H.; Sønstebø, J.H.; Nagy, N.E.; Timmermann, V.; Solheim, H.; Børja, I.; Kauserud, H.; Carlsen, T.; Rzepka, B.; Wasak, K.; et al. Fungal diversity and seasonal succession in ash leaves infected by the invasive ascomycete Hymenoscyphus fraxineus. New Phytol. 2017, 213, 1405–1417. [Google Scholar] [CrossRef]
  17. Griffiths, S.M.; Galambao, M.; Rowntree, J.; Goodhead, I.; Hall, J.; O’Brien, D.; Atkinson, N.; Antwis, R.E. Complex associations between cross-kingdom microbial endophytes and host genotype in ash dieback disease dynamics. J. Ecol. 2019, 108, 291–309. [Google Scholar] [CrossRef] [Green Version]
  18. Ulrich, K.; Becker, R.; Behrendt, U.; Kube, M.; Ulrich, A. A comparative analysis of ash leaf-colonizing bacterial communities identifies putative antagonists of Hymenoscyphus fraxineus. Front. Microbiol. 2020, 11, 966. [Google Scholar] [CrossRef]
  19. Ulrich, K.; Becker, R.; Behrendt, U.; Kube, M.; Schneck, V.; Ulrich, A. Physiological and genomic characterisation of Luteimonas fraxinea sp. nov., a bacterial species associated with trees tolerant to ash dieback. Syst. Appl. Microbiol. 2022, 45, 126333. [Google Scholar] [CrossRef]
  20. Comby, M.; Gacoin, M.; Robineau, M.; Rabenoelina, F.; Ptas, S.; Dupont, J.; Profizi, C.; Baillieul, F. Screening of wheat endophytes as biological control agents against Fusarium head blight using two different in vitro tests. Microbiol. Res. 2017, 202, 11–20. [Google Scholar] [CrossRef]
  21. Wemheuer, F.; Kaiser, K.; Karlovsky, P.; Daniel, R.; Vidal, S.; Wemheuer, B. Bacterial endophyte communities of three agricultural important grass species differ in their response towards management regimes. Sci. Rep. 2017, 7, 40914. [Google Scholar] [CrossRef] [Green Version]
  22. Kohout, C.E. The Plant Growth-Promoting Potential of Root-Associated Bacteria from Plants Growing in Stressed Environments. Ph.D. Thesis, University of Saskatchewan, Saskatoon, SK, Canada, 2019. [Google Scholar]
  23. Raaijmakers, J.M.; Mazzola, M. Diversity and natural functions of antibiotics produced by beneficial and plant pathogenic bacteria. Annu. Rev. Phytopathol. 2012, 50, 403–424. [Google Scholar] [CrossRef]
  24. Pieterse, C.M.J.; Zamioudis, C.; Berendsen, R.L.; Weller, D.M.; Wees, S.C.M.V.; Bakker, P.A.H.M. Induced systemic resistance by beneficial microbes. Annu. Rev. Phytopathol. 2014, 52, 347–375. [Google Scholar] [CrossRef] [Green Version]
  25. Kamilova, F.; Validov, S.; Azarova, T.; Mulders, I.; Lugtenberg, B. Enrichment for enhanced competitive plant root tip colonizers selects for a new class of biocontrol bacteria. Environ. Microbiol. 2005, 7, 1809–1817. [Google Scholar] [CrossRef] [PubMed]
  26. Biessy, A.; Filion, M. Phenazines in plant-beneficial Pseudomonas spp.: Biosynthesis, regulation, function and genomics. Environ. Microbiol. 2018, 20, 3905–3917. [Google Scholar] [CrossRef] [Green Version]
  27. Fira, D.; Dimkić, I.; Berić, T.; Lozo, J.; Stanković, S. Biological control of plant pathogens by Bacillus species. J. Biotechnol. 2018, 285, 44–55. [Google Scholar] [CrossRef] [PubMed]
  28. Kaspar, F.; Neubauer, P.; Gimpel, M. Bioactive secondary metabolites from Bacillus subtilis: A comprehensive review. J. Nat. Prod. 2019, 82, 2038–2053. [Google Scholar] [CrossRef] [PubMed]
  29. Grady, E.N.; MacDonald, J.; Liu, L.; Richman, A.; Yuan, Z.-C. Current knowledge and perspectives of Paenibacillus: A review. Microb. Cell Factories 2016, 15, 203. [Google Scholar] [CrossRef] [Green Version]
  30. Jurado, V.; Gonzalez, J.M.; Laiz, L.; Saiz-Jimenez, C. Aurantimonas altamirensis sp. nov., a member of the order Rhizobiales isolated from Altamira Cave. Int. J. Syst. Evol. Microbiol. 2006, 56, 2583–2585. [Google Scholar] [CrossRef] [Green Version]
  31. Luong, M.L.; Békal, S.; Vinh, D.C.; Lauzon, D.; Leung, V.; Al-Rawahi, G.N.; Ng, B.; Burdz, T.; Bernard, K. First report of isolation and characterization of Aurantimonas altamirensis from clinical samples. J. Clin. Microbiol. 2008, 46, 2435–2437. [Google Scholar] [CrossRef] [Green Version]
  32. Téllez-Castillo, C.J.; González Granda, D.; Bosch Alepuz, M.; Jurado Lobo, V.; Saiz-Jimenez, C.; Juan, J.L.; Millán Soria, J. Isolation of Aurantimonas altamirensis from pleural effusions. J. Med. Microbiol. 2010, 59, 1126–1129. [Google Scholar] [CrossRef] [Green Version]
  33. García-Lozano, T.; Oroval, E.A.; Bañón, J.L.J. First isolation in Spain of Aurantimonas altamirensis in a blood culture from a port-a-cath in a patient with Bence-Jones type multiple myeloma. Enferm. Infecc. Microbiol. Clin. 2012, 30, 217–218. [Google Scholar] [CrossRef]
  34. Schröttner, P.; Rudolph, W.W.; Taube, F.; Gunzer, F. First report on the isolation of Aureimonas altamirensis from a patient with peritonitis. Int. J. Infect. Dis. 2014, 29, 71–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Mendes, R.E.; Denys, G.A.; Fritsche, T.R.; Jones, R.N. Case report of Aurantimonas altamirensis bloodstream infection. J. Clin. Microbiol. 2009, 47, 514–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Eshaghi, A.; Shahinas, D.; Patel, S.N.; Kus, J.V. First draft genome sequence of Aureimonas altamirensis, isolated from patient blood culture. FEMS Microbiol. Lett. 2015, 362, fnv016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kim, N.; Hwang, J.H.; Cho, Y.G.; Kim, D.S.; Lee, H.S.; Lee, J. First case of Aureimonas altamirensis bacteremia in Korea. Ann. Lab. Med. 2019, 39, 587–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Zhao, S.; Chen, W.; Lan, Y.; Chen, M.; Li, J.; Guan, Y.; Hu, F.; Li, F.; Tang, X.; Li, L. First case of biliary Aureimonas altamirensis infection in a patient with colon cancer in China. Infect. Dis. Res. 2021, 2, 7–10. [Google Scholar] [CrossRef]
  39. Ventorino, V.; Aliberti, A.; Faraco, V.; Robertiello, A.; Giacobbe, S.; Ercolini, D.; Amore, A.; Fagnano, M.; Pepe, O. Exploring the microbiota dynamics related to vegetable biomasses degradation and study of lignocellulose-degrading bacteria for industrial biotechnological application. Sci. Rep. 2015, 5, 8161. [Google Scholar] [CrossRef] [Green Version]
  40. Liu, Y.; Xu, L.; Zhang, Z.; Huang, Z.; Fang, D.; Zheng, X.; Yang, Z.; Lu, M. Isolation, identification, and analysis of potential functions of culturable bacteria associated with an invasive gall wasp, Leptocybe invasa. Microb. Ecol. 2021, 83, 151–166. [Google Scholar] [CrossRef]
  41. Ganzorig, M.; Lee, K. Complete genome sequence of Aureimonas sp. strain OT7, isolated from human skin. Microbiol. Resour. Announc. 2021, 10, e00024-21. [Google Scholar] [CrossRef]
  42. Xiong, X.; Li, B.; Zhou, Z.; Gu, G.; Li, M.; Liu, J.; Jiao, H. The VirB system plays a crucial role in Brucella intracellular infection. Int. J. Mol. Sci. 2021, 22, 13637. [Google Scholar] [CrossRef]
  43. Ke, Y.; Wang, Y.; Li, W.; Chen, Z. Type IV secretion system of Brucella spp. and its effectors. Front. Cell. Infect. Microbiol. 2015, 5, 72. [Google Scholar] [CrossRef]
  44. Ratu, S.T.N.; Teulet, A.; Miwa, H.; Masuda, S.; Nguyen, H.P.; Yasuda, M.; Sato, S.; Kaneko, T.; Hayashi, M.; Giraud, E.; et al. Rhizobia use a pathogenic-like effector to hijack leguminous nodulation signalling. Sci. Rep. 2021, 11, 2034. [Google Scholar] [CrossRef] [PubMed]
  45. Stranahan, L.W.; Arenas-Gamboa, A.M. When the going gets rough: The significance of Brucella lipopolysaccharide phenotype in host–pathogen interactions. Front. Microbiol. 2021, 12, 713157. [Google Scholar] [CrossRef] [PubMed]
  46. 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]
  47. Pinski, A.; Betekhtin, A.; Hupert-Kocurek, K.; Mur, L.A.J.; Hasterok, R. Defining the genetic basis of plant–endophytic bacteria interactions. Int. J. Mol. Sci. 2019, 20, 1947. [Google Scholar] [CrossRef] [Green Version]
  48. Sessitsch, A.; Hardoim, P.R.; Döring, J.; Weilharter, A.; Krause, A.; Woyke, T.; Mitter, B.; Hauberg-Lotte, L.; Friedrich, F.; Rahalkar, M.C.; et al. Functional characteristics of an endophyte community colonizing rice roots as revealed by metagenomic analysis. Mol. Plant Microbe Interact. 2012, 25, 28–36. [Google Scholar] [CrossRef] [Green Version]
  49. del Barrio-Duque, A.; Samad, A.; Nybroe, O.; Antonielli, L.; Sessitsch, A.; Compant, S. Interaction between endophytic Proteobacteria strains and Serendipita indica enhances biocontrol activity against fungal pathogens. Plant Soil 2020, 451, 277–305. [Google Scholar] [CrossRef]
  50. Ulrich, K.; Kube, M.; Becker, R.; Schneck, V.; Ulrich, A. Genomic analysis of the endophytic Stenotrophomonas strain 169 reveals features related to plant-growth promotion and stress tolerance. Front. Microbiol. 2021, 12, 687463. [Google Scholar] [CrossRef]
  51. Chen, X.H.; Koumoutsi, A.; Scholz, R.; Eisenreich, A.; Schneider, K.; Heinemeyer, I.; Morgenstern, B.; Voss, B.; Hess, W.R.; Reva, O.; et al. Comparative analysis of the complete genome sequence of the plant growth–promoting bacterium Bacillus amyloliquefaciens FZB42. Nat. Biotechnol. 2007, 25, 1007–1014. [Google Scholar] [CrossRef] [Green Version]
  52. Muhammad, M.H.; Idris, A.L.; Fan, X.; Guo, Y.; Yu, Y.; Jin, X.; Qiu, J.; Guan, X.; Huang, T. Beyond Risk: Bacterial biofilms and their regulating approaches. Front. Microbiol. 2020, 11, 928. [Google Scholar] [CrossRef]
  53. Poulin, M.B.; Kuperman, L.L. Regulation of biofilm exopolysaccharide production by cyclic di-guanosine monophosphate. Front. Microbiol. 2021, 12, 730980. [Google Scholar] [CrossRef]
  54. Green, E.R.; Mecsas, J. Bacterial secretion systems: An overview. Microbiol. Spectr. 2016, 4, 4.1.13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Tseng, T.-T.; Tyler, B.M.; Setubal, J.C. Protein secretion systems in bacterial-host associations, and their description in the Gene Ontology. BMC Microbiol. 2009, 9, S2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Meneses, C.H.; Rouws, L.F.; Simoes-Araujo, J.L.; Vidal, M.S.; Baldani, J.I. Exopolysaccharide production is required for biofilm formation and plant colonization by the nitrogen-fixing endophyte Gluconacetobacter diazotrophicus. Mol. Plant Microbe Interact. 2011, 24, 1448–1458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Kubota, T.; Miyamoto, K.; Yasuda, M.; Inamori, Y.; Tsujibo, H. Molecular characterization of an intracellular β-n-acetylglucosaminidase involved in the chitin degradation system of Streptomyces thermoviolaceus OPC-520. Biosci. Biotechnol. Biochem. 2004, 68, 1306–1314. [Google Scholar] [CrossRef] [Green Version]
  58. Vermassen, A.; Leroy, S.; Talon, R.; Provot, C.; Popowska, M.; Desvaux, M. Cell wall hydrolases in bacteria: Insight on the diversity of cell wall amidases, glycosidases and peptidases toward peptidoglycan. Front. Microbiol. 2019, 10, 331. [Google Scholar] [CrossRef] [Green Version]
  59. Alcorlo, M.; Martínez-Caballero, S.; Molina, R.; Hermoso, J.A. Carbohydrate recognition and lysis by bacterial peptidoglycan hydrolases. Curr. Opin. Struct. Biol. 2017, 44, 87–100. [Google Scholar] [CrossRef]
  60. Höltje, J.-V. From growth to autolysis: The murein hydrolases in Escherichia coli. Arch. Microbiol. 1995, 164, 243–254. [Google Scholar] [CrossRef]
  61. Herlihey, F.A.; Clarke, A.J. Controlling autolysis during flagella insertion in Gram-negative bacteria. Adv. Exp. Med. Biol. 2017, 925, 41–56. [Google Scholar]
  62. Skujiņś, J.; Puķite, A.; McLaren, A.D. Adsorption and reactions of chitinase and lysozyme on chitin. Mol. Cell. Biochem. 1973, 2, 221–228. [Google Scholar] [CrossRef]
  63. Koskella, B. The phyllosphere. Curr. Biol. 2020, 30, R1143–R1146. [Google Scholar] [CrossRef]
  64. Bremer, E.; Krämer, R. Responses of microorganisms to osmotic stress. Annu. Rev. Microbiol. 2019, 73, 313–334. [Google Scholar] [CrossRef] [PubMed]
  65. Boch, J.; Kempf, B.; Schmid, R.; Bremer, E. Synthesis of the osmoprotectant glycine betaine in Bacillus subtilis: Characterization of the gbsAB genes. J. Bacteriol. 1996, 178, 5121–5129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Boch, J.; Nau-Wagner, G.; Kneip, S.; Bremer, E. Glycine betaine aldehyde dehydrogenase from Bacillus subtilis: Characterization of an enzyme required for the synthesis of the osmoprotectant glycine betaine. Arch. Microbiol. 1997, 168, 282–289. [Google Scholar] [CrossRef]
  67. Lin, J.C.; Gant, N. Chapter 2.3—The Biochemistry of Choline. In Magnetic Resonance Spectroscopy; Stagg, C., Rothman, D., Eds.; Academic Press: San Diego, CA, USA, 2014; pp. 104–110. [Google Scholar]
  68. Herrou, J.; Willett, J.W.; Czyż, D.M.; Babnigg, G.; Kim, Y.; Crosson, S. Conserved ABC transport system regulated by the general stress response pathways of alpha- and gammaproteobacteria. J. Bacteriol. 2017, 199, e00746-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Dalton, T.P.; Shertzer, H.G.; Puga, A. Regulation of gene expression by reactive oxygen. Annu. Rev. Pharmacol. Toxicol. 1999, 39, 67–101. [Google Scholar] [CrossRef]
  70. Imlay, J.A. Cellular defenses against superoxide and hydrogen peroxide. Annu. Rev. Biochem. 2008, 77, 755–776. [Google Scholar] [CrossRef] [Green Version]
  71. Wang, Y.; Branicky, R.; Noë, A.; Hekimi, S. Superoxide dismutases: Dual roles in controlling ROS damage and regulating ROS signaling. J. Cell Biol. 2018, 217, 1915–1928. [Google Scholar] [CrossRef] [Green Version]
  72. Veal, E.A.; Toone, W.M.; Jones, N.; Morgan, B.A. Distinct roles for glutathione S-transferases in the oxidative stress response in Schizosaccharomyces pombe. J. Biol. Chem. 2002, 277, 35523–35531. [Google Scholar] [CrossRef] [Green Version]
  73. Zhuge, X.-L.; Xu, H.; Xiu, Z.-J.; Yang, H.-L. Biochemical functions of glutathione S-transferase family of Salix babylonica. Front. Plant Sci. 2020, 11, 364. [Google Scholar] [CrossRef] [Green Version]
  74. Ianieva, O. Mechanisms of bacteria resistance to heavy metals. Mikrobiol. Z. 2009, 71, 54–65. [Google Scholar]
  75. Nanda, M.; Kumar, V.; Sharma, D.K. Multimetal tolerance mechanisms in bacteria: The resistance strategies acquired by bacteria that can be exploited to ‘clean-up’ heavy metal contaminants from water. Aquat. Toxicol. 2019, 212, 121–136. [Google Scholar] [CrossRef]
  76. Pal, A.; Bhattacharjee, S.; Saha, J.; Sarkar, M.; Mandal, P. Bacterial survival strategies and responses under heavy metal stress: A comprehensive overview. Crit. Rev. Microbiol. 2021, 48, 327–355. [Google Scholar] [CrossRef] [PubMed]
  77. Tiwari, S.; Lata, C. Heavy metal stress, signaling, and tolerance due toplant-associated microbes: An overview. Front. Plant Sci. 2018, 9, 452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Franco-Franklin, V.; Moreno-Riascos, S.; Ghneim-Herrera, T. Are endophytic Bacteria an option for increasing heavy metal tolerance of plants? A meta-analysis of the effect size. Front. Environ. Sci. 2021, 8, 603668. [Google Scholar] [CrossRef]
  79. Queloz, V.; Grünig, C.R.; Berndt, R.; Kowalski, T.; Sieber, T.N.; Holdenrieder, O. Cryptic speciation in Hymenoscyphus albidus. Forest Pathol. 2011, 41, 133–142. [Google Scholar] [CrossRef]
  80. Peters, S.; Langer, G.; Kätzel, R. Eschentriebsterben—Kriterien zur Schadensbonitur an Eschen; Fachagentur Nachwachsende Rohstoffe e.V. (FNR): Gülzow-Prüzen, Germany, 2021. [Google Scholar]
  81. Mitter, B.; Pfaffenbichler, N.; Flavell, R.; Compant, S.; Antonielli, L.; Petric, A.; Berninger, T.; Naveed, M.; Sheibani-Tezerji, R.; von Maltzahn, G.; et al. A new approach to modify plant microbiomes and traits by introducing beneficial bacteria at flowering into progeny seeds. Front. Microbiol. 2017, 8, 11. [Google Scholar] [CrossRef] [Green Version]
  82. Andreote, F.D.; da Rocha, U.N.; Araujo, W.L.; Azevedo, J.L.; van Overbeek, L.S. Effect of bacterial inoculation, plant genotype and developmental stage on root-associated and endophytic bacterial communities in potato (Solanum tuberosum). Anton. Leeuw. Int. J. G. 2010, 97, 389–399. [Google Scholar] [CrossRef]
  83. Elsayed, T.R.; Jacquiod, S.; Nour, E.H.; Sørensen, S.J.; Smalla, K. Biocontrol of bacterial wilt disease through complex interaction between tomato plant, antagonists, the indigenous rhizosphere microbiota, and Ralstonia solanacearum. Front. Microbiol. 2020, 10, 2835. [Google Scholar] [CrossRef]
  84. Ardanov, P.; Sessitsch, A.; Haggman, H.; Kozyrovska, N.; Pirttila, A.M. Methylobacterium-induced endophyte community changes correspond with protection of plants against pathogen attack. PLoS ONE 2012, 7, e46802. [Google Scholar] [CrossRef]
  85. Jiménez, J.A.; Novinscak, A.; Filion, M. Inoculation with the plant-growth-promoting rhizobacterium Pseudomonas fluorescens LBUM677 Impacts the rhizosphere microbiome of three oilseed crops. Front. Microbiol. 2020, 11, 569366. [Google Scholar] [CrossRef]
  86. Berg, G.; Kusstatscher, P.; Abdelfattah, A.; Cernava, T.; Smalla, K. Microbiome modulation—Toward a better understanding of plant microbiome response to microbial inoculants. Front. Microbiol. 2021, 12, 650610. [Google Scholar] [CrossRef] [PubMed]
  87. Chowdhury, S.P.; Hartmann, A.; Gao, X.; Borriss, R. Biocontrol mechanism by root-associated Bacillus amyloliquefaciens FZB42—A review. Front. Microbiol. 2015, 6, 780. [Google Scholar] [CrossRef] [PubMed]
  88. Kröber, M.; Wibberg, D.; Grosch, R.; Eikmeyer, F.; Verwaaijen, B.; Chowdhury, S.P.; Hartmann, A.; Pühler, A.; Schlüter, A. Effect of the strain Bacillus amyloliquefaciens FZB42 on the microbial community in the rhizosphere of lettuce under field conditions analyzed by whole metagenome sequencing. Front. Microbiol. 2014, 5, 252. [Google Scholar] [CrossRef] [PubMed]
  89. Ulrich, K.; Ulrich, A.; Ewald, D. Diversity of endophytic bacterial communities in poplar grown under field conditions. FEMS Microbiol. Ecol. 2008, 63, 169–180. [Google Scholar] [CrossRef] [Green Version]
  90. Tokajian, S.; Issa, N.; Salloum, T.; Ibrahim, J.; Farah, M. 16S-23S rRNA gene intergenic spacer region variability helps resolve closely related sphingomonads. Front. Microbiol. 2016, 7, 149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Hunt, M.; Silva, N.D.; Otto, T.D.; Parkhill, J.; Keane, J.A.; Harris, S.R. Circlator: Automated circularization of genome assemblies using long sequencing reads. Genome Biol. 2015, 16, 294. [Google Scholar] [CrossRef] [Green Version]
  92. Overbeek, R.; Olson, R.; Pusch, G.D.; Olsen, G.J.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Parrello, B.; Shukla, M.; et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014, 42, D206–D214. [Google Scholar] [CrossRef]
  93. Li, W.; O’Neill, K.R.; Haft, D.H.; DiCuccio, M.; Chetvernin, V.; Badretdin, A.; Coulouris, G.; Chitsaz, F.; Derbyshire, M.K.; Durkin, A.S.; et al. RefSeq: Expanding the Prokaryotic Genome Annotation Pipeline reach with protein family model curation. Nucleic Acids Res. 2021, 49, D1020–D1028. [Google Scholar] [CrossRef]
  94. Chen, H. Venndiagram: Generate High-Resolution Venn and Euler Plots. R Package Version 1.6.20. Available online: https://CRAN.R-project.org/package=VennDiagram (accessed on 15 March 2019).
  95. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  96. Parks, D.H.; Chuvochina, M.; Waite, D.W.; Rinke, C.; Skarshewski, A.; Chaumeil, P.-A.; Hugenholtz, P. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat. Biotechnol. 2018, 36, 996–1004. [Google Scholar] [CrossRef]
  97. Berg, G.; Martinez, J.L. Friends or foes: Can we make a distinction between beneficial and harmful strains of the Stenotrophomonas maltophilia complex? Front. Microbiol 2015, 6, 241. [Google Scholar] [CrossRef] [PubMed]
  98. Buchanan, J.T.; Simpson, A.J.; Aziz, R.K.; Liu, G.Y.; Kristian, S.A.; Kotb, M.; Feramisco, J.; Nizet, V. DNase expression allows the pathogen group a Streptococcus to escape killing in neutrophil extracellular traps. Curr. Biol. 2006, 16, 396–400. [Google Scholar] [CrossRef] [PubMed]
  99. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: http://www.R-project.org/ (accessed on 1 December 2021).
  100. Callahan, B.J.; McMurdie, P.J.; Holmes, S.P. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017, 11, 2639–2643. [Google Scholar] [CrossRef] [PubMed]
Plants 11 03487 g001 550
Figure 1. Variability of the A. altamirensis isolates obtained from ash trees: (a) Ordination plot of MALDI-TOF MS spectra, calculated by a principal component analysis. (b) Maximum-likelihood tree of the IGS gene sequences of selected ash-associated A. altamirensis isolates and two clinically derived strains (ON-56566, OT7) as well as the type strain DSM 21988. Numbers at branch nodes refer to bootstrap values >50%.
Figure 1. Variability of the A. altamirensis isolates obtained from ash trees: (a) Ordination plot of MALDI-TOF MS spectra, calculated by a principal component analysis. (b) Maximum-likelihood tree of the IGS gene sequences of selected ash-associated A. altamirensis isolates and two clinically derived strains (ON-56566, OT7) as well as the type strain DSM 21988. Numbers at branch nodes refer to bootstrap values >50%.
Plants 11 03487 g001
Plants 11 03487 g002 550
Figure 2. Maximum-likelihood tree of the 16S rRNA gene sequences showing the position of A. altamirensis C2P003 among other A. altamirensis strains and reference strains of related Aureimonas species. Phylogenetic analysis was performed using the HKY+G+I model. Rhizobium leguminosarum was used as the outgroup. Numbers at branch nodes refer to bootstrap values >50%. Bar: substitutions per nucleotide site. Accession numbers or locus tags (NCBI or IMG database) are indicated in brackets.
Figure 2. Maximum-likelihood tree of the 16S rRNA gene sequences showing the position of A. altamirensis C2P003 among other A. altamirensis strains and reference strains of related Aureimonas species. Phylogenetic analysis was performed using the HKY+G+I model. Rhizobium leguminosarum was used as the outgroup. Numbers at branch nodes refer to bootstrap values >50%. Bar: substitutions per nucleotide site. Accession numbers or locus tags (NCBI or IMG database) are indicated in brackets.
Plants 11 03487 g002
Plants 11 03487 g003 550
Figure 3. Venn diagram showing orthologous genes between the genomes of A. altamirensis C2P003, DSM 21988, ON-56566 and OT7.
Figure 3. Venn diagram showing orthologous genes between the genomes of A. altamirensis C2P003, DSM 21988, ON-56566 and OT7.
Plants 11 03487 g003
Plants 11 03487 g004 550
Figure 4. Impacts of inoculation with A. altamirensis C2P003, B. velezensis A4P130 and L. fraxinea D4P002 on ash seedlings. The effects on ash dieback were monitored by the health status of H. fraxineus-infected seedlings. Health status and ash dieback symptoms were estimated as described by Peters et al. [80]. Bonitur scale: 1: healthy, 2: 10–25% leaf damage, 3: 26–60% leaf damage, 4: 61-99 leaf damage, 5: dead plants. Stars indicate significant differences compared to the control (p < 0.05).
Figure 4. Impacts of inoculation with A. altamirensis C2P003, B. velezensis A4P130 and L. fraxinea D4P002 on ash seedlings. The effects on ash dieback were monitored by the health status of H. fraxineus-infected seedlings. Health status and ash dieback symptoms were estimated as described by Peters et al. [80]. Bonitur scale: 1: healthy, 2: 10–25% leaf damage, 3: 26–60% leaf damage, 4: 61-99 leaf damage, 5: dead plants. Stars indicate significant differences compared to the control (p < 0.05).
Plants 11 03487 g004
Plants 11 03487 g005 550
Figure 5. Composition of the bacterial microbiome of ash seedlings inoculated with A. altamirensis C2P003, B. velezensis A4P130 and L. fraxinea D4P002 as well as control plants. Two different taxonomic levels are demonstrated.
Figure 5. Composition of the bacterial microbiome of ash seedlings inoculated with A. altamirensis C2P003, B. velezensis A4P130 and L. fraxinea D4P002 as well as control plants. Two different taxonomic levels are demonstrated.
Plants 11 03487 g005
Plants 11 03487 g006 550
Figure 6. Ordination plot visualizing the bacterial microbiome composition of ash seedlings inoculated with A. altamirensis C2P003, B. velezensis A4P130 and L. fraxinea D4P002 as well as control plants. A principle coordinate analysis (PCoA) was applied based on the weighted UniFraq distance matrix.
Figure 6. Ordination plot visualizing the bacterial microbiome composition of ash seedlings inoculated with A. altamirensis C2P003, B. velezensis A4P130 and L. fraxinea D4P002 as well as control plants. A principle coordinate analysis (PCoA) was applied based on the weighted UniFraq distance matrix.
Plants 11 03487 g006
Table
Table 1. Genome annotation statistics of the Aureimonas altamirensis strains C2P003, DSM 21988, ON-56566 and OT7.
Table 1. Genome annotation statistics of the Aureimonas altamirensis strains C2P003, DSM 21988, ON-56566 and OT7.
C2P003DSM 21988ON-56566OT7
Contigs111131
Genome size (bp)4,592,9814,190,9654,200,0474,181,223
G+C content (%)64.364.865.265.0
Number of protein-coding genes 14289396739763871
rRNAs 16S (total)3 (9)2 (10)1 (5)3 (9)
tRNAs50464851
1 The data were derived from the NCBI and IMG databases, and the accession numbers are CP089311, 2585427611, 2627854086 and CP062167, respectively.
Table
Table 2. Summary of protein-coding genes involved in host colonization and stress adaptation in the genome of A. altamirensis C2P003 and corresponding genes in the genomes of DSM 21988 and ON-56566.
Table 2. Summary of protein-coding genes involved in host colonization and stress adaptation in the genome of A. altamirensis C2P003 and corresponding genes in the genomes of DSM 21988 and ON-56566.
Category/SubsystemNumber of Genes
C2P003DSM 21988ON-56566
Capsular and extracellular polysaccharides887
Membrane Transport
Protein secretion system, Type I330
Protein secretion system, Type IV45935
Protein secretion system, Type VII330
Stress response: Osmotic stress
Osmoprotectant ABC transporter YehZYXW of Enterobacteriales444
Synthesis of osmoregulated periplasmic glucans333
Choline and betaine uptake and betaine biosynthesis999
Stress response: Oxidative stress
Oxidative stress787
Glutathione: Biosynthesis and gamma-glutamyl cycle635
Glutathione: Nonredox reactions999
Glutathione: Redox cycle333
Resistance to toxic compounds
Copper homeostasis/tolerance81011
Cobalt-zinc-cadmium resistance888
Resistance to chromium compounds111
Resistance to antibiotics
Resistance to fluoroquinolones222
Beta-lactamase111
Table
Table 3. Predicted plant and microbial cell wall-degrading enzymes of strains A. altamirensis C2P003, DSM 21988 and ON-56566.
Table 3. Predicted plant and microbial cell wall-degrading enzymes of strains A. altamirensis C2P003, DSM 21988 and ON-56566.
CAZy FamilySubstrateAnnotationEnzyme CodeNumber of Genes
C2P003DSM 21988ON-56566
CE4PeptidoglycansPeptidoglycan N-acetylglucosamine deacetylaseEC 3.5.1.-100
CE4PolysaccharidesPolysaccharide deacetylase 200
PeptidoglycansD-alanyl-D-alanine carboxypeptidaseEC 3.4.16.4444
GH73PeptidoglycansN-acetylmuramoyl-L-alanine amidaseEC 3.5.1.28111
CE9Polysaccharides/ChitooligosaccharidesN-acetylglucosamine-6-phosphate deacetylaseEC 3.5.1.25111
PolysaccharidesGlucosamine-6-phosphate deaminaseEC 3.5.99.6111
GH13Polysaccharidesα-amylaseEC 3.2.1.1111
GH15PolysaccharidesGlucoamylaseEC 3.2.1.3221
GT35PolysaccharidesGlycogen phosphorylaseEC 2.4.1.1111
GH43Hemicelluloseα-L-arabinofuranosidaseEC 3.2.1.55101
GH9Celluloseβ-1,4-glucanaseEC 3.2.1.4111
GH28Polysaccharides/PectinPectin degradation protein 110
GH24Peptidoglycans/ChitinPhage lysozyme REC 3.2.1.17211
GH73Peptidoglycans/Chitinβ-N-acetylglucosaminidaseEC 3.2.1.52111
GH15PolysaccharidesGlucan 1,4-α-glucosidaseEC 3.2.1.3110
GH16Cellulose, HemicelluloseEndo-β-1,3-1,4 glucanase (licheninase)EC 3.2.1.73111
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Becker, R.; Ulrich, K.; Behrendt, U.; Schneck, V.; Ulrich, A. Genomic Characterization of Aureimonas altamirensis C2P003—A Specific Member of the Microbiome of Fraxinus excelsior Trees Tolerant to Ash Dieback. Plants 2022, 11, 3487. https://doi.org/10.3390/plants11243487

AMA Style

Becker R, Ulrich K, Behrendt U, Schneck V, Ulrich A. Genomic Characterization of Aureimonas altamirensis C2P003—A Specific Member of the Microbiome of Fraxinus excelsior Trees Tolerant to Ash Dieback. Plants. 2022; 11(24):3487. https://doi.org/10.3390/plants11243487

Chicago/Turabian Style

Becker, Regina, Kristina Ulrich, Undine Behrendt, Volker Schneck, and Andreas Ulrich. 2022. "Genomic Characterization of Aureimonas altamirensis C2P003—A Specific Member of the Microbiome of Fraxinus excelsior Trees Tolerant to Ash Dieback" Plants 11, no. 24: 3487. https://doi.org/10.3390/plants11243487

APA Style

Becker, R., Ulrich, K., Behrendt, U., Schneck, V., & Ulrich, A. (2022). Genomic Characterization of Aureimonas altamirensis C2P003—A Specific Member of the Microbiome of Fraxinus excelsior Trees Tolerant to Ash Dieback. Plants, 11(24), 3487. https://doi.org/10.3390/plants11243487

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