Next Article in Journal
Rapid and Label-Free Immunosensing of Shiga Toxin Subtypes with Surface Plasmon Resonance Imaging
Next Article in Special Issue
Vaccination with VLPs Presenting a Linear Neutralizing Domain of S. aureus Hla Elicits Protective Immunity
Previous Article in Journal
Accumulation and Elimination of Tetrodotoxin in the Pufferfish Takifugu obscurus by Dietary Administration of the Wild Toxic Gastropod Nassarius semiplicata
Previous Article in Special Issue
Host–Receptor Post-Translational Modifications Refine Staphylococcal Leukocidin Cytotoxicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Staphylococcus aureus Pathogenicity in Cystic Fibrosis Patients—Results from an Observational Prospective Multicenter Study Concerning Virulence Genes, Phylogeny, and Gene Plasticity

1
Institute of Medical Microbiology, University Hospital Münster, 48149 Münster, Germany
2
Department of Medical Microbiology and Infectious Diseases, Erasmus Medical Center Rotterdam, 3015 CN Rotterdam, The Netherlands
3
Institute of Biostatistics and Clinical Research, University Hospital Münster, 48149 Münster, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally as senior authors.
Submission received: 6 April 2020 / Revised: 16 April 2020 / Accepted: 21 April 2020 / Published: 26 April 2020
(This article belongs to the Special Issue Staphylococcus aureus Toxins: Promoter or Handicap during Infection)

Abstract

:
Staphylococcus aureus and cystic fibrosis (CF) are closely interlinked. To date, however, the impact of S. aureus culture in CF airways on lung function and disease progression has only been elucidated to a limited degree. This analysis aims to identify bacterial factors associated to clinical deterioration. Data were collected during an observational prospective multi-center study following 195 patients from 17 centers. The average follow-up time was 80 weeks. S. aureus isolates (n = 3180) were scanned for the presence of 25 virulence genes and agr-types using single and multiplex PCR. The presence of specific virulence genes was not associated to clinical deterioration. For the agr-types 1 and 4, however, a link to the subjects’ clinical status became evident. Furthermore, a significant longitudinal decrease in the virulence gene quantity was observed. Analyses of the plasticity of the virulence genes revealed significantly increased plasticity rates in the presence of environmental stress. The results suggest that the phylogenetic background defines S. aureus pathogenicity rather than specific virulence genes. The longitudinal loss of virulence genes most likely reflects the adaptation process directed towards a persistent and colonizing rather than infecting lifestyle.
Key Contribution: A longitudinal decline in S. aureus virulence gene quantity was observed. The presence of specific virulence genes characteristic for S. aureus was not associated to clinical deterioration.

1. Introduction

Cystic fibrosis (CF) is the most prevalent life-shortening genetic disorder in the Caucasian population, characterized by increased viscosity of exocrine gland secretions, congestion of physiological outflow pathways, and consecutive cystic and fibrotic remodeling of the respective organs [1,2]. Chronic and recurrent bacterial airway infections are the main cause of morbidity and mortality among CF patients [3].
Staphylococcus aureus is one of the first pathogens that CF patients encounter and the bacterium with the highest prevalence in the CF population [4,5]. In most cases, a long-term colonization with individual S. aureus clones is observed [4]. Persistent S. aureus clones tend to form a genetically and particularly phenotypically heterogeneous population with S. aureus acquisition most likely outside the hospital with only minor transmission events [6,7,8]. The genetical diversity is fostered by mutations and horizontal gene transfer with persistence of dominant lineages and a continuous process of intra-host adaptation [6,9]. Both geno- and phenotypic changes occur in a time-dependent manner and accumulate [10,11]. The phenotypic heterogeneity of S. aureus populations in CF patients is intensified by differences in gene expression, mediated by regulatory mechanisms responding to specific milieu conditions [6]. The overall trend of S. aureus adaptation in CF airways is directed towards a persistent and colonizing rather than infecting lifestyle [6,12]. Even in the case of exacerbation, the presence of S. aureus is restricted to the airways and systemic disease is rarely observed [6]. The preponderance of colonizing clones results from the suppression of virulence genes for the sake of an intensified expression of protective factors [6]. Manifold strategies have been described to evade immune response and antibiotic treatment, including biofilm formation [13,14] emergence of small colony variants (SCVs) [15,16,17], intracellular persistence [18], and agr-dysfunction [19], among others [14,20,21,22].
In children, the colonization with S. aureus is associated with increased inflammatory activity in the airways and a worse nutritional status [23,24,25]. For older patients, however, the impact of S. aureus carriage on lung function and disease progression has not yet been elucidated [26]. This deficit in knowledge is reflected by the absence of internationally accepted therapeutic guidelines concerning S. aureus treatment in CF patients [27,28].
In a prospective longitudinal multicenter study, we recently identified the presence of SCVs, high bacterial density in throat cultures, and elevated IL-6 levels as independent risk factors among others for worse lung function in patients who were older than 6 years [26]. Further analyses revealed associations between antibiotic treatment, age, and S. aureus clonal carriage profiles during persistence assessed by spa-sequence typing of almost 4000 S. aureus isolates from 183 CF patients during a 21-month prospective study [29]. Patients receiving intensive antibiotic therapy regimens more likely carried S. aureus isolates during persistence, which belonged to only one spa-type [7], while older patients carried isolates with related spa-types with mutations of the repeat region including deletions, duplications of repeats or point mutations within repeats, or dominant spa-types cultured [7].
The present analysis continues the evaluation of the abovementioned study intending to identify genetic bacterial factors associated to clinical deterioration. We hypothesized that either the presence of specific S. aureus virulence genes or the occurrence of certain phylogenetic backgrounds might be associated with a deterioration of the patients’ clinical status. To test our hypothesis, we scanned every isolate for the presence of 25 virulence genes (cap 5/8, chp, cna, clf A/B, sdrC/D/E, fnbA/B, sasG/H, eap, emp, sea–sej, eta/etb, tst-1, pvl, and hlg) and the agr-specificity groups 1–4, using mono and/or multiplex PCR. Furthermore, the host response to 19 of these virulence genes was assessed by measuring IgG levels against the specific proteins by a bead-based flow cytometry technique (xMAP1; Luminex Corporation, Austin, TX, USA) as previously described [30,31].

2. Results

2.1. Demographics and Statistics

The CF patients participating in the study on average were younger than the German CF population and showed better lung function [26]. Furthermore, there was an overrepresentation of nonF508del genotypes and male sex [26]. Fifty-eight patients (30.5%) suffered an exacerbation during the study period, whereas the clinical status of the remaining 132 (69.5%) subjects did not deteriorate considerably. Subjects suffering one or more exacerbations showed a poorer lung function (mean 75% FEV1% predicted) than healthier patients (mean 84% FEV1% predicted) (p = 0.034 *, Mann–Whitney U test). In total, 3180 S. aureus isolates were analyzed (1312 from throat cultures, 1065 from nose cultures, and 803 from sputum cultures). On average, 3.22 (range: 1–12) different S. aureus clones as determined by spa-typing (Supplementary Materials Table S1 were retrieved from the subjects’ airways. The majority of patients (108; 57.8%) carried at least one dominant clone that was detected on every outpatient clinic visit. Six spa-types were most prevalent within this S. aureus population belonging to the spa-types t084, t012, t091, t015, and t002 (Supplementary Materials Table S2).

2.2. Quantity of Virulence Genes and Clinical Parameters

The starting point for this analysis was the hypothesis that a broad virulence gene profile might be associated to a deterioration of the patients’ clinical parameters. However, there was no statistical difference in mean virulence gene quantity between patients suffering from an exacerbation during the study period (mean of 14.09 virulence genes) and subjects not witnessing exacerbations (mean of 14.29 virulence genes).

2.3. Time-Dependent Quantity of Virulence Genes

The number of virulence genes present in the S. aureus clones populating the patients’ airways was documented at every outpatient visit, enabling a longitudinal analysis concerning time-dependent changes in the quantity of virulence genes exhibited by the patients’ S. aureus clones. At the beginning of the study period, the patients’ S. aureus clones carried on average 14.51 of the 25 virulence genes, whereas at the end, a mean of 14.17 virulence genes was present. A generalized estimating equation model was applied using Poisson distribution with log link function. Appropriate model terms for repeating observations were included using an independent working correlation matrix. Within the model, we could show a significant decrease in virulence gene quantity (p = 0.029 *, Figure 1) over time. On average, patients lost 5.6 × 10 4 genes per day. Extrapolated (beyond the maximal follow-up of this study), this would translate to a loss of one gene every five years (exactly, 1.022/5 years).

2.4. Specific Virulence Genes and Clinical Parameters

The next step consisted in the analysis of a possible association between single virulence genes and clinical parameters. Therefore, the mean prevalence counts of the individual virulence genes in different patient subgroups were compared. Patients suffering an exacerbation during the study period (n = 58) showed similar prevalence counts of the different virulence genes compared to healthier subjects (n = 132) (Figure 2a–c). The only exception was patients with isolates belonging to agr-type 4, showing a significant overrepresentation among patients with an exacerbation (2.34% vs. 7.72%) (p = 0.008 **, Mann–Whitney U test, Figure 2c). Agr-type 1 showed a considerable underrepresentation among patients witnessing an exacerbation (p = 0.051, Mann–Whitney U test). To narrow the spectrum of considered subjects to patients with the presumably most severe disease, patients with intensive antibiotic therapy and exacerbation (n = 35) were compared to the other subjects. In this subgroup, isolates belonging to agr-type 1 were significantly underrepresented (p = 0.011 *, Mann–Whitney U test, Figure 2d).

2.5. Complementary Data Concerning a Humoral Response Towards S. aureus Antigens

S. aureus toxin directed antibody levels in CF patients have been associated to the severity of exacerbations and decline in FEV1% [32]. Therefore, we measured IgG levels against 19 virulence factors in serum samples of 182 subjects and compared them to antibody levels of 53 healthy adults with persistent nasal S. aureus colonization (Supplementary Materials Table S3). While IgG levels against six antigens were significantly elevated in patients compared to controls (Table 1), IgG levels against seven virulence factors were significantly associated to a decline in FEV1% predicted (Table 1).

2.6. Phylogeny and Virulence Genes

In the preceding analyses, it became obvious that neither the number of virulence genes present nor the presence of specific virulence genes was associated to changes in the patients´ clinical parameters. Nevertheless, the significant differences of the mean prevalence counts of the agr-types 1 and 4 among patients with an exacerbation suggest an influence of the phylogenetic background on varying degrees of S. aureus pathogenicity in CF patients. To further pursue these hints, we tested for a link between agr-type and phylogeny.
Agr-type and virulence gene profile: Clones, which belonged to varying agr-types, differed significantly concerning the counts of virulence genes (p = 0.001 ***, ANOVA, Supplementary Materials Figure S1). Agr-type 4 isolates carried less virulence genes compared to isolates belonging to agr-type 1. Next, the prevalence counts of the specific virulence genes in patients whose S. aureus clones predominantly belonged to one of the four agr-types were compared. In this analysis, clear and significant differences in almost all considered virulence genes were found (Kruskal–Wallis test, Figure 3a). Benjamini and Hochberg (False discovery rate) correction was applied to adjust for multiple testing. In a final analysis, it was tested for associations between the most prevalent spa- and agr-types. This revealed that each of the most prevalent spa-types could be clearly assigned to one agr-type (p < 0.001, Chi-squared test, Figure 3b).
Agr-type and longitudinal virulence gene quantity: In our preceding analysis, we showed that the counts of virulence genes present in the patients’ S. aureus clones demonstrated a time-dependent decrease (Figure 1). Reconducting this analysis and taking into account the agr-type of the patients’ isolates revealed clear and significant differences in the process of time-dependent alterations of the virulence gene quantity between isolates belonging to the different agr-types (p = 0.0083 **, Wald test). Isolates belonging to agr-type 1 and 3 showed a longitudinal decline in virulence gene quantity, whereas the virulence gene number in isolates belonging to agr-type 2 slightly increased and agr-type 4 isolates remained stable in virulence gene numbers (Figure 4).

2.7. Cluster Analyses

Since the isolates belonging to the different agr-types significantly varied in the prevalence of many virulence genes (Figure 2), it was considered that certain combinations of virulence genes rather than single virulence genes might be associated to a deterioration in clinical parameters. To test this hypothesis, a distance matrix was created deploying Euclidian distance and then hierarchical clustering via ward linkage was applied. The virulence gene profiles were grouped into four association clusters, implying that the virulence gene profiles of patients within a cluster are more similar to each other than to virulence gene profiles of other clusters. In a second step, we analyzed for associations between these clusters and the different agr-types, rendering significant results (p = 0.001 ***, Chi-squared test, Figure 5). Next, it was investigated whether exacerbation rates differed significantly between the four virulence gene clusters, which was not the case.

2.8. Virulence Gene Plasticity and Environmental Stress

Environmental stress consisting of antibiotic therapy and/or immune response has been shown to be a powerful inducer of hgt, in particular of bacteriophage mediated transduction [33,34,35]. Since this process can entail clinical consequences and can be induced iatrogenically, the data were checked for hints pointing towards stress-induced virulence gene plasticity in the study population. The patients were divided using the “environmental stress high/low” parameter. Patients assigned to the “environmental stress high” group showed at least one exacerbation during the study period and received intensive antibiotic treatment; 633 S. aureus isolates of 35 patients were exposed to high amounts of environmental stress. There were significant differences in plasticity for three out of 25 tested virulence genes (chp, eta, etb), all of them showing increased plasticity in case of high environmental stress (Mann–Whitney U test, Figure 6).

3. Discussion

A prospective longitudinal multi-centre study consisting of 195 patients receiving treatment in 17 specialized CF centres was conducted to identify genetic bacterial factors associated to clinical deterioration in the patient cohort. It was tested whether the presence of specific S. aureus virulence genes or the occurrence of certain phylogenetic backgrounds were linked to differing clinical status of the subjects. Furthermore, the data were analysed for longitudinal changes in the S. aureus virulence gene composition and stress dependent virulence gene plasticity.
The analyses revealed several interesting findings that will be discussed in the following.

3.1. Specific Virulence Genes and Clinical Parameters

None of the virulence genes characteristic of S. aureus (tst-1 and lukS/F, among others) was associated to exacerbations. This finding supports the hypothesis that, given no susceptible host, the presence of single virulence genes is not associated with a significant increase in the overall severity of S. aureus infections and, consequently, that a screening for such “super weapons” is not required [36,37,38,39]. These results, however, only apply to S. aureus disease in CF patients and neither could nor should be generalized. Moreover, there are significant regional differences in the prevalence counts of the various virulence genes. For instance, community-acquired methicillin-resistant S. aureus (CA-MRSA) strains, which mostly carry the Panton-Valentine leucocidin (PVL, lukS/F), are rare in Germany while widespread in the US [40,41,42]. The significant differences in the prevalence counts of the agr-types 1 and 4 between patients with and without exacerbations most likely did not arise from differing intrinsic gene regulatory activities of the varying agr-types, since the agr-system has been described to be inactive in CF patients most of the times [19,43]. Rather, the observation of widespread agr-inactivity or dysfunction in CF airways suggests a link between the agr-system and S. aureus phylogenetic background as an explanation for differing gene prevalence counts in exacerbating patients [39].

3.2. Phylogeny and Pathogenicity

The identification of the different agr-types as the only parameter that was associated to the subjects’ clinical status supports the hypothesis that the basic unit of S. aureus pathogenicity is not the presence of specific single virulence genes but most likely certain combinations of genes linked to a clone’s phylogenetic background [39,44]. The phylogenetic indicator function of the agr-system is highlighted by further analyses, revealing considerably varying virulence gene profiles between patients whose S. aureus clones differed concerning their agr-types (Figure 3a) and the observation that each prevalent spa-type could be clearly assigned to one agr-type (Figure 3b). The results of the cluster analyses added further coherence to these findings: the different agr-types showed a significantly unequal distribution between the four virulence gene clusters (Figure 5).

3.3. Longitudinal Decline in Virulence Gene Quantity

The analyses revealed a significant time-dependent decline in the number of virulence genes exhibited by the subjects’ S. aureus clones (Figure 1). The decrease may seem negligible at first glance; however, it has to be taken into account that this analysis is highly limited in two dimensions. On the one hand, it only covers a short cutout of the longstanding interaction between S. aureus and its human host. On the other hand, only a small fraction of the extensive S. aureus virulence genome is included. Thus, a hypothetical extrapolation concerning follow-up time and considered virulence genes would put the results into a different perspective. Hence, this finding is considered as a clear hint towards a longitudinal decline of S. aureus virulence gene quantity in CF patients and thus indicates that S. aureus shows not only high phenotypic diversity and increased degrees of gene plasticity in CF airways but also long-term changes in its (virulence) gene composition. The finding is also consistent with the persistent and colonizing rather than infectious lifestyle, which S. aureus displays in CF patients [6]. Similar processes have also been proven for chronic P. aeruginosa infection in CF patients [45,46]. The gradual loss of virulence genes could be beneficial for the bacterium concerning the evasion of host defenses, since many virulence factors are targets for the immune system as obvious by mounting of specific IgG against such virulence factors (Table 1) [45]. Future research has to show whether there is evidence for this trend on the level of a more comprehensive virulence gene profile analysis using whole genome sequencing (WGS). While Azarian et al. used WGS to investigate intra-host diversity of sequential MRSA isolates in four CF patients, the authors focused in their analysis, for which they used sequencing data of a mean of 21 isolates per patient from a mean period of 1403 days on the evidence of switched intergenic regions, allele frequency, and single nucleotide polymorphisms but not the presence or absence of special virulence genes [9]. Therefore, it is not possible to compare their data with our analysis.

3.4. Virulence Gene Plasticity and Environmental Stress

In this analysis, the virulence gene plasticity was used as a surrogate parameter to estimate the potential for hgt. It is considered a reasonable parameter for this purpose since genome plasticity forms a necessary precondition for hgt. Furthermore, surroundings that are associated to increased rates of virulence gene plasticity in the data have also been identified as inductors of hgt [47,48]. The analyses revealed significantly increased plasticity rates for several genes in the presence of high environmental stress (Figure 6). These results are in line with research having identified environmental stress as an inducer of hgt and suggest an iatrogenic and host-related impact on the S. aureus virulence genome [34,35]. According to the “insurance hypothesis”, genomic plasticity is a necessary prerequisite for adaptation to changing and challenging surroundings [49]. That these surroundings themselves increase the potential for hgt is an important fact that treating physicians should be aware of. Particularly sub-inhibitory concentrations of antibiotics, which are often observed to occur in CF patients, have been shown to be powerful inducers of hgt [50].

3.5. Strengths and Limitations of the Study and Analyses

The study population differs from the German CF population concerning age, sex, lung function, and genotype [26]. The “No P. aeruginosa” inclusion criterion is considered to be the reason for the subjects’ young age, since with increasing life span patients become more prone to colonization and infection with P. aeruginosa [5]. The young age in turn presumably accounts for the subjects’ better lung function compared to the German CF population, since it is an age-correlated parameter characterized by a longitudinal decline [51,52]. The “No P. aeruginosa” criterion presumably also accounts for the increased prevalence of non-F508del genotypes since those genotypes have been shown to exhibit a delayed colonization by P. aeruginosa compared to F508del homozygous patients [53]. There are no clear hints that these demographic incongruities go along with a limited transferability of the findings to the basic CF population. Both groups show a similar annual decline in lung function and thus do not seem to differ substantially in the course of the disease [26,51,52]. A shortcoming of the study is the rather short follow-up time of app. 80 weeks, which only covers a narrow cutout of the long-lasting interaction between S. aureus and its human host. Moreover, the analyses are primarily based on genetic data. Complementary proteomic data would be desirable, since differences in gene expression contribute to the high level of phenotypic diversity among S. aureus clones [54]. Especially the diversity of the S. aureus exoproteome is considered to be relevantly shaped by gene regulation mechanisms [55,56].

4. Conclusions

The results of our prospective multicenter study suggest that the phylogenetic background defines S. aureus pathogenicity rather than specific virulence genes. The significantly increased gene plasticity rates in the presence of environmental stress suggest an iatrogenic and host influence on S. aureus genome plasticity. The longitudinal loss of virulence genes might reflect the adaptation process directed towards a persistent and colonizing rather than infecting behavior of S. aureus in CF patients.

5. Material and Methods

5.1. Study Design

The data underlying this analysis were generated during an prospective multicenter study (ClinicalTrials.gov, NCT00669760). An ethical statement was obtained at the main study center with the central laboratory in Münster, Germany (2007-496-f-S) on the 1 June 2009. Written informed consent was obtained from all patients and from parents if patients were younger than 18 years.
The study population consisted of 195 CF patients receiving treatment in 16 German and one Austrian CF center. CF patients older than six years whose airways were persistently colonized by S. aureus during the study period (more than 50% of specimen positive for S. aureus) were included into the trial. The follow-up time was on average 80 weeks, with a mean of seven outpatient clinic visits per patient, usually performed at an interval of three months. Inclusion criteria were the presence of a minimum of two airway cultures positive for S. aureus during an interval of six months in the year before study enrollment, and if more than four cultures were available, at least 50% of them had to harbor S. aureus. To prevent possible confounding caused by the presence of bacteria other than S. aureus, potential subjects showing specimens positive for Pseudomonas aeruginosa and/or Burkholderia cepacia complex (BCC) were to be ruled out a priori. At each outpatient visit, nasal and throat and/or sputum samples were retrieved and sent to the central study laboratory in Münster, where the material was cultivated following standardized procedures for CF airway specimens as described elsewhere [57]. Furthermore, antibiotic regimens and clinical parameters were monitored and documented longitudinally. For a more detailed description of the study, see Junge et al. [26]. The obtained S. aureus isolates (n = 3893) were distinguished deploying the spa-typing method [58]. The “Based upon repeat pattern” method (BURP, Ridom StaphType software) was used to analyze the clonal relatedness of different spa-types differentiating 266 spa types (Supplementary Materials Tables S1 and S2) [29,59]. In this analysis, every isolate was scanned for the presence of 25 virulence genes and the 4 agr-specificity groups, using single and multiplex PCR.
In the present study, 3180 S. aureus isolates obtained from 190 patients were analysed; five subjects were excluded because no persistent S. aureus has been detected in their specimens; and 713 isolates of the former 3893 isolates [33] were disregarded due to identity of isolates in regard of the same visit, sample (nose or throat swab or sputum), spa-type, and virulence gene profile.

5.2. Virulence Gene Profiles

The obtained S. aureus isolates (n = 3180) were scanned for the presence of 25 virulence genes and the four agr-specificity groups (cap 5/8, chip, cna, clf A/B, sdrC/D/E, fnbA/B, sasG/H, eap, emp, seasej, eta/etb, tst-1, pvl, hlg, and agr), using single and multiplex PCRs as described elsewhere (Supplementary Materials Table S4 [60,61,62,63,64,65,66,67]).

5.3. Parameters

The study provided an extensive data set, consisting on the one hand of the patients’ clinical data and on the other hand of the results of the microbiological analyses. Since the aim of this analysis was to screen for correlations between these two different data domains, reasoned parameters were to be established to render the data accessible for statistical analysis, while at the same time, oversimplification had to be avoided to generate data that is of use for clinicians. The parameters were chosen and defined as follows:

5.3.1. Exacerbation

The treating physicians monitored the patient’s clinical status according to the Fuchs criteria (change in sputum/sinus discharge, new or increased dyspnea/cough/hemoptysis/malaise/weight loss/sinus tenderness, temperature over 38 °C, radiographic changes indicative of lung infection, and decrease in pulmonary function by 10%) [68]. If at one visit four or more Fuchs criteria were met, the respective patient was assigned to the “exacerbation” group.

5.3.2. Lung Function

Spirometric lung function testing was conducted on every outpatient clinic visit. Lung function was expressed as forced expiratory volume (1 s) (FEV1) percent of normal. Comparison with a healthy cohort was achieved by using the Quanjer prediction model [69]. During the study’s follow-up, lung function changes over time were negligible. Thus, the mean lung function (FEV1%) throughout all outpatient clinic visits was calculated for every patient.

5.3.3. Virulence Gene Prevalence Counts

To analyze for associations between the presence of specific virulence genes and a patient’s clinical outcome, the calculation of individual virulence gene prevalence counts was necessary, which display the percentage of isolates in which the respective virulence gene was present throughout the study period.

5.3.4. Quantity of Virulence Genes

In addition to a possible clinical impact of single virulence genes, it was also evaluated whether a broad or narrow spectrum of virulence genes in general correlated with changes in clinical outcomes. Thus, it was reviewed for every visit how many of the 25 virulence genes were present in the patients’ S. aureus clones. This count was conducted in a complementary way, meaning that, at every visit, each proven virulence gene was only counted once per clone. This parameter also enabled a longitudinal analysis of a possible time-dependent change in the number of virulence genes incorporated by the S. aureus clones populating the patients’ airways.

5.3.5. Agr-Type

The accessory gene regulator (agr) not only represents one of S. aureus’s most important virulence regulators but also is tightly linked to the bacteria’s phylogenetic background [39,44,70,71]. To elicit possible relations between agr-type, virulence gene prevalence counts, and spa-types, the isolates of patients were analyzed in terms of dominance of a specific agr-type. To belong to the group with dominant agr-types, at least a two-third majority of patients’ isolates had to belong to one of the four agr-types, which was the case for 154 patients. Patients without a dominance of one agr-type (36 subjects) were excluded from further analyses regarding this parameter.

5.3.6. Environmental Stress

Previous research has shown an impact of antibiotic therapy and host immune response on genomic variability and horizontal gene transfer (hgt) [33,35,50]. Thus, the parameter “environmental stress” was established, which includes all patients who suffered an exacerbation during the study period and received intensive antibiotic therapy. The latter was defined as ≥7 cycles of antibiotic therapy during the study period. The environmental stress parameter was used to compare gene plasticity and hence the potential for hgt in S. aureus isolates that were exposed to high amounts of environmental stress to isolates subsisting in a less challenging environment.

5.3.7. Virulence Gene Plasticity and Horizontal Gene Transfer

To quantify the plasticity of the virulence genes, the virulence profiles of the patients’ S. aureus clones were compared from visit to visit. Events of gain and loss of virulence genes were counted, resulting in a value for every virulence gene of every clone of each subject, representing its variability in this specific clone of this specific patient. To receive a variability value for every virulence gene for every patient independent of the different clones, the cumulative sum of genetic switch events of each clone of the patient was calculated, resulting in a variability value for every virulence gene for each patient, which later was set into relation to other parameters. The virulence gene plasticity was used as a surrogate parameter to estimate the potential for hgt.

5.3.8. Antistaphylococcal Antibodies

IgG antibody levels were quantified applying a bead-based flow cytometry technique (xMAP ®; Luminex Corporation) as previously described [30,31]. For further information, see Junge et al. [26].

5.4. Statistical Analyses

Study data were analyzed using descriptive and inductive statistical methods. For data description, absolute and relative frequencies were used for categorical variables while continuous parameters were described by mean and standard deviation or by median and interquartile range. To check for normal distribution, histogram plots and the Kolmogorov–Smirnov test were applied. For group comparisons, t-tests and ANOVA were applied. In the case of non-normal distributions nonparametric methods (e.g., Mann–Whitney U tests) were applied. Associations between categorical variables were tested with Chi-squared tests. Regression coefficients were tested using the Wald-test statistic. Further details of the applied statistical methods are given in line with the results of the presented analyses.
The statistical analyses were conducted using SPSS (IBM Corp. Released 2017. IBM SPSS Statistics for Windows, Version 25.0. Armonk, NY, USA: IBM Corp.). The significance level for all performed statistical tests is 0.05. All statistical tests have to be considered exploratory. Results, thus, might need to be validated in prospective trials. Significant p-values, thus, indicate noticeable effects and not significant results in a confirmatory sense. For illustrations, significance levels are indicated by *, **, and ***.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/5/279/s1, Figure S1: Agr-type vs. virulence gene count, Table S1: spa-types of all S. aureus isolates; Table S2: Prevalent spa-types; Table S3: Prevalent spa-types; Table S4: Primers for single and multiplex PCRs.

Author Contributions

Conceptualization, B.C.K.; Methodology, B.C.K., J.L., W.v.W., M.d.R. and D.G.; Software, J.L. and D.G.; Validation, J.L., M.d.R., W.v.W., D.G., and B.C.K.; Formal analysis, J.L., K.H. (Kathrin Heidenreich), K.H. (Katharina Higelin), K.D., V.M., C.R., M.d.R., W.v.W., D.G., and B.C.K.; Investigation, J.L., K.H. (Kathrin Heidenreich), K.H. (Katharina Higelin), K.D., V.M., C.R., D.G., W.v.W., M.d.R. and B.C.K.; resources, B.C.K.; Data curation, J.L., K.H. (Kathrin Heidenreich), K.H. (Katharina Higelin), K.D., V.M., C.R., D.G., and B.C.K.; Writing—Original draft preparation, J.L., D.G., and B.C.K.; Writing—Review and editing, J.L., M.d.R., W.v.W., D.G., and B.C.K.; Visualization, J.L. and D.G.; Supervision, W.v.W., D.G. and B.C.K.; Project administration, B.C.K.; Funding acquisition, B.C.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a grant from the Mukoviszidose e.V. (S05/07), Germany, and partly by grants of the Interdisciplinary Center for Clinical Research (IZKF Münster; Kah2/016/16) and the Transregional Collaborative Research Center 34 (C7) to B.C.K. We acknowledge support from the Open Access Publication Fund of the University of Muenster.

Acknowledgments

We thank the colleagues of the participating centers for recruiting patients and the CF patients and their parents for participating in our observational study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elborn, J.S. Cystic fibrosis. Lancet 2016, 388, 2519–2531. [Google Scholar] [CrossRef]
  2. Farrell, P.M. The prevalence of cystic fibrosis in the European Union. J. Cyst. Fibros. 2008, 7, 450–453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Rosenfeld, M.; Gibson, R.L.; McNamara, S.; Emerson, J.; Burns, J.L.; Castile, R.; Hiatt, P.; McCoy, K.; Wilson, C.B.; Inglis, A.; et al. Early pulmonary infection, inflammation, and clinical outcomes in infants with cystic fibrosis. Pediatr. Pulmonol. 2001, 32, 356–366. [Google Scholar] [CrossRef] [PubMed]
  4. Kahl, B.C. Impact of Staphylococcus aureus on the pathogenesis of chronic cystic fibrosis lung disease. Int. J. Med. Microbiol. 2010, 300, 514–519. [Google Scholar] [CrossRef]
  5. Barley, M.; Mcnally, J.; Marshall, B.; Faro, A.; Elbert, A.; Fink, A.; Sewall, A.; Loeffler, D.; Petren, K.; O’Neil, T.; et al. Annual Data Report 2016 Cystic Fibrosis Foundation Patient Registry. Cyst. Fibros. Found. Patient Regist. 2016, 2016, 1–92. [Google Scholar]
  6. Goerke, C.; Wolz, C. Adaptation of Staphylococcus aureus to the cystic fibrosis lung. Int. J. Med. Microbiol. 2010, 300, 520–525. [Google Scholar] [CrossRef]
  7. Westphal, C.; Görlich, D.; Kampmeier, S.; Herzog, S.; Braun, N.; Hitschke, C.; Mellmann, A.; Peters, G.; Kahl, B.C.; Junge, S.; et al. Antibiotic Treatment and Age Are Associated with Staphylococcus aureus Carriage Profiles During Persistence in the Airways of Cystic Fibrosis Patients. Front. Microbiol. 2020, 11, 230. [Google Scholar] [CrossRef]
  8. Ankrum, A.; Hall, B.G. Population dynamics of Staphylococcus aureus in cystic fibrosis patients to determine transmission events by use of whole-genome sequencing. J. Clin. Microbiol. 2017, 55, 2143–2152. [Google Scholar] [CrossRef] [Green Version]
  9. Azarian, T.; Ridgway, J.P.; Yin, Z.; David, M.Z. Long-Term Intrahost Evolution of Methicillin Resistant Staphylococcus aureus Among Cystic Fibrosis Patients With Respiratory Carriage. Front. Genet. 2019, 10, 1–14. [Google Scholar] [CrossRef] [Green Version]
  10. Hirschhausen, N.; Block, D.; Bianconi, I.; Bragonzi, A.; Birtel, J.; Lee, J.C.; Dübbers, A.; Küster, P.; Kahl, J.; Peters, G.; et al. Extended Staphylococcus aureus persistence in cystic fibrosis is associated with bacterial adaptation. Int. J. Med. Microbiol. 2013, 303, 685–692. [Google Scholar] [CrossRef]
  11. Langhanki, L.; Berger, P.; Treffon, J.; Catania, F.; Kahl, B.C.; Mellmann, A. In vivo competition and horizontal gene transfer among distinct Staphylococcus aureus lineages as major drivers for adaptational changes during long-term persistence in humans. BMC Microbiol. 2018, 18, 152. [Google Scholar] [CrossRef] [PubMed]
  12. Goerke, C.; Gressinger, M.; Endler, K.; Breitkopf, C.; Wardecki, K.; Stern, M.; Wolz, C.; Kahl, B.C. High phenotypic diversity in infecting but not in colonizing Staphylococcus aureus populations. Environ. Microbiol. 2007, 9, 3134–3142. [Google Scholar] [CrossRef] [PubMed]
  13. Donlan, R.M.; Costerton, J.W. Biofilms: Survival Mechanisms of Clinically Relevant Microorganisms. Clin. Microbiol. Rev. 2002, 15, 167–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Schwartbeck, B.; Birtel, J.; Treffon, J.; Langhanki, L.; Mellmann, A.; Kale, D.; Kahl, J.; Hirschhausen, N.; Neumann, C.; Lee, J.C.; et al. Dynamic in vivo mutations within the ica operon during persistence of Staphylococcus aureus in the airways of cystic fibrosis patients. PLoS Pathog. 2016, 12, e1006024. [Google Scholar] [CrossRef] [PubMed]
  15. Kahl, B.; Herrmann, M.; Everding, A.S.; Koch, H.G.; Becker, K.; Harms, E.; Proctor, R.A.; Peters, G. Persistent Infection with Small Colony Variant Strains of Staphylococcus aureus in Patients with Cystic Fibrosis. J. Infect. Dis. 1998, 177, 1023–1029. [Google Scholar] [CrossRef]
  16. Wolter, D.J.; Emerson, J.C.; McNamara, S.; Buccat, A.M.; Qin, X.; Cochrane, E.; Houston, L.S.; Rogers, G.B.; Marsh, P.; Prehar, K.; et al. Staphylococcus aureus Small-Colony Variants Are Independently Associated with Worse Lung Disease in Children with Cystic Fibrosis. Clin. Infect. Dis. 2013, 57, 384–391. [Google Scholar] [CrossRef]
  17. Besier, S.; Ludwig, A.; Ohlsen, K.; Brade, V.; Wichelhaus, T.A. Molecular analysis of the thymidine-auxotrophic small colony variant phenotype of Staphylococcus aureus. Int. J. Med. Microbiol. 2007, 297, 217–225. [Google Scholar] [CrossRef]
  18. Garzoni, C.; Kelley, W.L. Staphylococcus aureus: New evidence for intracellular persistence. Trends Microbiol. 2009, 17, 59–65. [Google Scholar] [CrossRef]
  19. Goerke, C.; Campana, S.; Bayer, M.G.; Döring, G.; Botzenhart, K.; Wolz, C.; Do, G.; Botzenhart, K.; Wolz, C. Direct quantitative transcript analysis of the agr regulon of Staphylococcus aureus during human infection in comparison to the expression profile in vitro. Infect. Immun. 2000, 68, 1304–1311. [Google Scholar] [CrossRef] [Green Version]
  20. Herzog, S.; Dach, F.; de Buhr, N.; Niemann, S.; Schlagowski, J.; Chaves-Moreno, D.; Neumann, C.; Goretzko, J.; Schwierzeck, V.; Mellmann, A.; et al. High nuclease activity of long persisting Staphylococcus aureus isolates within the airways of cystic fibrosis patients protects against NET-mediated killing. Front. Immunol. 2019, 10, 2552. [Google Scholar] [CrossRef] [Green Version]
  21. Treffon, J.; Block, D.; Moche, M.; Reiss, S.; Fuchs, S.; Engelmann, S.; Becher, D.; Langhanki, L.; Mellmann, A.; Peters, G.; et al. Adaptation of Staphylococcus aureus to airway environments in patients with cystic fibrosis by upregulation of superoxide dismutase M and iron-scavenging proteins. J. Infect. Dis. 2018, 217, 1453–1461. [Google Scholar] [CrossRef] [PubMed]
  22. Treffon, J.; Chaves-Moreno, D.; Niemann, S.; Pieper, D.H.; Vogl, T.; Roth, J.; Kahl, B.C. Importance of superoxide dismutases A and M for protection of Staphylococcus aureus in the oxidative stressful environment of cystic fibrosis airways. Cell. Microbiol. 2020, 22, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Gangell, C.; Gard, S.; Douglas, T.; Park, J.; De Klerk, N.; Keil, T.; Brennan, S.; Ranganathan, S.; Robins-Browne, R.; Sly, P.D. Inflammatory responses to individual microorganisms in the lungs of children with cystic fibrosis. Clin. Infect. Dis. 2011, 53, 425–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Armstrong, D.S.; Hook, S.M.; Jamsen, K.M.; Nixon, G.M.; Carzino, R.; Carlin, J.B.; Robertson, C.F.; Grimwood, K. Lower airway inflammation in infants with cystic fibrosis detected by newborn screening. Pediatr. Pulmonol. 2005, 40, 500–510. [Google Scholar] [CrossRef]
  25. Sagel, S.D.; Gibson, R.L.; Emerson, J.; McNamara, S.; Burns, J.L.; Wagener, J.S.; Ramsey, B.W.; Konstan, M.; Chatfield, B.; Retsch-Bogart, G.; et al. Impact of Pseudomonas and Staphylococcus Infection on Inflammation and Clinical Status in Young Children with Cystic Fibrosis. J. Pediatr. 2009, 154, 183–188. [Google Scholar] [CrossRef] [Green Version]
  26. Junge, S.; Görlich, D.; Den Reijer, M.; Wiedemann, B.; Tümmler, B.; Ellemunter, H.; Dübbers, A.; Küster, P.; Ballmann, M.; Koerner-Rettberg, C.; et al. Factors associated with worse lung function in cystic fibrosis patients with persistent Staphylococcus aureus. PLoS ONE 2016, 11, e0166220. [Google Scholar] [CrossRef]
  27. Smyth, A. Prophylactic antibiotics in cystic fibrosis: A conviction without evidence? Pediatr. Pulmonol. 2005, 40, 471–476. [Google Scholar] [CrossRef]
  28. Hurley, M.N.; Smyth, A.R. Staphylococcus aureus in cystic fibrosis: Pivotal role or bit part actor? Curr. Opin. Pulm. Med. 2018, 24, 586–591. [Google Scholar] [CrossRef]
  29. Westphal, C.; Görlich, D.; Herzog, S.; Braun, N.; Vogel, C.; Peters, G.; Kahl, B. Antibiotic treatment and age have an impact on Staphylococcus aureus clonality during persistence in the airways of cystic fibrosis patients. Front Microbiol. 2020, 11, 230. [Google Scholar] [CrossRef]
  30. den Reijer, P.M.; Lemmens-den Toom, N.; Kant, S.; Snijders, S.V.; Boelens, H.; Tavakol, M.; Verkaik, N.J.; van Belkum, A.; Verbrugh, H.A.; van Wamel, W.J.B. Characterization of the Humoral Immune Response during Staphylococcus aureus Bacteremia and Global Gene Expression by Staphylococcus aureus in Human Blood. PLoS ONE 2013, 8, e53391. [Google Scholar] [CrossRef] [Green Version]
  31. Verkaik, N.J.; de Vogel, C.P.; Boelens, H.A.; Grumann, D.; Hoogenboezem, T.; Vink, C.; Hooijkaas, H.; Foster, T.J.; Verbrugh, H.A.; van Belkum, A.; et al. Anti-Staphylococcal Humoral Immune Response in Persistent Nasal Carriers and Noncarriers of Staphylococcus aureus. J. Infect. Dis. 2009, 199, 625–632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Chadha, A.D.; Thomsen, I.P.; Jimenez-Truque, N.; Soper, N.R.; Jones, L.S.; Sokolow, A.G.; Torres, V.J.; Creech, C.B. Host response to Staphylococcus aureus cytotoxins in children with cystic fibrosis. J. Cyst. Fibros. 2016, 15, 597–604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Goerke, C.; Koller, J.; Wolz, C. Ciprofloxacin and trimethoprim cause phage induction and virulence modulation in Staphylococcus aureus. Antimicrob. Agents Chemother. 2006, 50, 171–177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Goerke, C.; Matias y Papenberg, S.; Dasbach, S.; Dietz, K.; Ziebach, R.; Kahl, B.C.; Wolz, C. Increased frequency of genomic alterations in Staphylococcus aureus during chronic infection is in part due to phage mobilization. J. Infect. Dis. 2004, 189, 724–734. [Google Scholar] [CrossRef] [Green Version]
  35. Úbeda, C.; Maiques, E.; Knecht, E.; Lasa, Í.; Novick, R.P.; Penadés, J.R. Antibiotic-induced SOS response promotes horizontal dissemination of pathogenicity island-encoded virulence factors in staphylococci. Mol. Microbiol. 2005, 56, 836–844. [Google Scholar] [CrossRef]
  36. Lindsay, J.A.; Moore, C.E.; Day, N.P.; Peacock, S.J.; Witney, A.A.; Stabler, R.A.; Husain, S.E.; Butcher, P.D.; Hinds, J. Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes. J. Bacteriol. 2006, 188, 669–676. [Google Scholar] [CrossRef] [Green Version]
  37. Shallcross, L.J.; Fragaszy, E.; Johnson, A.M.; Hayward, A.C. The role of the Panton-Valentine leucocidin toxin in staphylococcal disease: A systematic review and meta-analysis. Lancet Infect. Dis. 2013, 13, 43–54. [Google Scholar] [CrossRef] [Green Version]
  38. Laabei, M.; Uhlemann, A.C.; Lowy, F.D.; Austin, E.D.; Yokoyama, M.; Ouadi, K.; Feil, E.; Thorpe, H.A.; Williams, B.; Perkins, M.; et al. Evolutionary Trade-Offs Underlie the Multi-faceted Virulence of Staphylococcus aureus. PLoS Biol. 2015, 13, e1002229. [Google Scholar] [CrossRef] [Green Version]
  39. Jarraud, S.; Mougel, C.; Thioulouse, J.; Lina, G.; Meugnier, H.; Forey, F.; Nesme, X.; Etienne, J.; Vandenesch, F. Relationships between Staphylococcus aureus genetic background, virulence factors, agr groups (alleles), and human disease. Infect. Immun. 2002, 70, 631–641. [Google Scholar] [CrossRef] [Green Version]
  40. Cafiso, V.; Bertuccio, T.; Spina, D.; Campanile, F.; Bongiorno, D.; Santagati, M.; Sciacca, A.; Sciuto, C.; Stefani, S. Methicillin resistance and vancomycin heteroresistance in Staphylococcus aureus in cystic fibrosis patients. Eur. J. Clin. Microbiol. Infect. Dis. 2010, 29, 1277–1285. [Google Scholar] [CrossRef]
  41. Shilo, N.; Quach, C. Pulmonary Infections and Community Associated Methicillin Resistant Staphylococcus aureus: A Dangerous Mix? Paediatr. Respir. Rev. 2011, 12, 182–189. [Google Scholar] [CrossRef] [PubMed]
  42. Glikman, D.; Siegel, J.D.; David, M.Z.; Okoro, N.M.; Boyle-Vavra, S.; Dowell, M.L.; Daum, R.S. Complex Molecular Epidemiology of Methicillin-Resistant Staphylococcus aureus Isolates from Children with Cystic Fibrosis in the Era of Epidemic Community-Associated Methicillin-Resistant S aureus. Chest 2008, 133, 1381–1387. [Google Scholar] [CrossRef] [PubMed]
  43. Goerke, C.; Kümmel, M.; Dietz, K.; Wolz, C. Evaluation of intraspecies interference due to agr polymorphism in Staphylococcus aureus during infection and colonization. J. Infect. Dis. 2003, 188, 250–256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wright, J.S.; Traber, K.E.; Corrigan, R.; Benson, S.A.; Musser, J.M.; Novick, R.P. The agr Radiation: An Early Event in the Evolution of Staphylococci. J. Bacteriol. 2005, 187, 5585–5594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Nguyen, D.; Singh, P.K. Evolving stealth: Genetic adaptation of Pseudomonas aeruginosa during cystic fibrosis infections. Proc. Natl. Acad. Sci. USA 2006, 103, 8305–8306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Smith, E.E.; Buckley, D.G.; Wu, Z.; Saenphimmachak, C.; Hoffman, L.R.; D’Argenio, D.A.; Miller, S.I.; Ramsey, B.W.; Speert, D.P.; Moskowitz, S.M.; et al. Genetic adaptation by Pseudomonas aeruginosa to the airways of cystic fibrosis patients. Proc. Natl. Acad. Sci. USA 2006, 103, 8487–8492. [Google Scholar] [CrossRef] [Green Version]
  47. Knöppel, A.; Lind, P.A.; Lustig, U.; Näsvall, J.; Andersson, D.I. Minor fitness costs in an experimental model of horizontal gene transfer in bacteria. Mol. Biol. Evol. 2014, 31, 1220–1227. [Google Scholar] [CrossRef] [Green Version]
  48. Sorek, R.; Zhu, Y.; Creevey, C.J.; Francino, M.P.; Bork, P.; Rubin, E.M. Genome-Wide Experimental Determination of Barriers to Horizontal Gene Transfer. Science 2007, 318, 1449–1452. [Google Scholar] [CrossRef] [Green Version]
  49. Yachi, S.; Loreau, M. Biodiversity and ecosystem productivity in a fluctuating environment: The insurance hypothesis. Proc. Natl. Acad. Sci. USA 1999, 96, 1463–1468. [Google Scholar] [CrossRef] [Green Version]
  50. Hodille, E.; Warre, R.; Diep, B.A.; Goutelle, S.; Lina, G.; Dumitrescu, O. The role of antibiotics in modulating virulence in Staphylococcus aureus. Clin. Microbiol. Rev. 2017, 30, 887–917. [Google Scholar] [CrossRef] [Green Version]
  51. Liou, T.G.; Elkin, E.P.; Pasta, D.J.; Jacobs, J.R.; Konstan, M.W.; Morgan, W.J.; Wagener, J.S. Year-to-year changes in lung function in individuals with cystic fibrosis. J. Cyst. Fibros. 2010, 9, 250–256. [Google Scholar] [CrossRef] [Green Version]
  52. Kerem, E.; Viviani, L.; Zolin, A.; MacNeill, S.; Hatziagorou, E.; Ellemunter, H.; Drevinek, P.; Gulmans, V.; Krivec, U.; Olesen, H. Factors associated with FEV1 decline in cystic fibrosis: Analysis of the ECFS patient registry. Eur. Respir. J. 2014, 43, 125–133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. McKone, E.F.; Emerson, S.S.; Edwards, K.L.; Aitken, M.L. Effect of genotype on phenotype and mortality in cystic fibrosis: A retrospective cohort study. Lancet 2003, 361, 1671–1676. [Google Scholar] [CrossRef]
  54. Wolf, C.; Kusch, H.; Monecke, S.; Albrecht, D.; Holtfreter, S.; von Eiff, C.; Petzl, W.; Rainard, P.; Bröker, B.M.; Engelmann, S. Genomic and proteomic characterization of Staphylococcus aureus mastitis isolates of bovine origin. Proteomics 2011, 11, 2491–2502. [Google Scholar] [CrossRef] [PubMed]
  55. Ziebandt, A.K.; Kusch, H.; Degner, M.; Jaglitz, S.; Sibbald, M.J.J.B.; Arends, J.P.; Chlebowicz, M.A.; Albrecht, D.; Pantuček, R.; Doškar, J.; et al. Proteomics uncovers extreme heterogeneity in the Staphylococcus aureus exoproteome due to genomic plasticity and variant gene regulation. Proteomics 2010, 10, 1634–1644. [Google Scholar] [CrossRef]
  56. Sibbald, M.J.J.B.; Ziebandt, A.K.; Engelmann, S.; Hecker, M.; de Jong, A.; Harmsen, H.J.M.; Raangs, G.C.; Stokroos, I.; Arends, J.P.; Dubois, J.Y.F.; et al. Mapping the Pathways to Staphylococcal Pathogenesis by Comparative Secretomics. Microbiol. Mol. Biol. Rev. 2006, 70, 755–788. [Google Scholar] [CrossRef] [Green Version]
  57. Hogardt, M.; Haeußler, S.; Balke, B.; Kahl, B.C.; Schmoldt, S.; Leitritz, L.; Jäger, G.; Kappler, M.; Suerbaum, S.; Heesemann, J. MIQ 24: Atemwegsinfektionen bei Mukoviszidose Qualitätsstandards in der Mikrobiologisch-Infektiologischen Diagnostik; Elsevier, Urban Fischer: München, Germany, 2006. [Google Scholar]
  58. Harmsen, D.; Claus, H.; Witte, W.; Rothganger, J.; Claus, H.; Turnwald, D.; Vogel, U. Typing of Methicillin-Resistant Staphylococcus aureus in a University Hospital Setting by Using Novel Software for spa Repeat Determination and Database Management. J. Clin. Microbiol. 2003, 41, 5442–5448. [Google Scholar] [CrossRef] [Green Version]
  59. Mellmann, A.; Weniger, T.; Berssenbrügge, C.; Rothgänger, J.; Sammeth, M.; Stoye, J.; Harmsen, D. Based Upon Repeat Pattern (BURP): An algorithm to characterize the long-term evolution of Staphylococcus aureus populations based on spa polymorphisms. BMC Microbiol. 2007, 7. [Google Scholar] [CrossRef] [Green Version]
  60. van Wamel, W.J.B.; Rooijakkers, S.H.M.; Ruyken, M.; van Kessel, K.P.M.; van Strijp, J.A.G. The Innate Immune Modulators Staphylococcal Complement Inhibitor and Chemotaxis Inhibitory Protein of Staphylococcus aureus Are Located on -Hemolysin-Converting Bacteriophages. J. Bacteriol. 2006, 188, 1310–1315. [Google Scholar] [CrossRef] [Green Version]
  61. Gomes, A.R.; Vinga, S.; Zavolan, M.; Lencastre, H. De Analysis of the Genetic Variability of Virulence-Related Loci in Epidemic Clones of Methicillin-Resistant Staphylococcus aureus Analysis of the Genetic Variability of Virulence-Related Loci in Epidemic Clones of Methicillin-Resistant Staphylococcus aureus. Antimicrob. Agents Chemother. 2005, 49, 366–379. [Google Scholar]
  62. Goerke, C.; Esser, S.; Kümmel, M.; Wolz, C. Staphylococcus aureus strain designation by agr and cap polymorphism typing and delineation of agr diversification by sequence analysis. Int. J. Med. Microbiol. 2005, 295, 67–75. [Google Scholar] [CrossRef] [PubMed]
  63. Roche, F.M.; Massey, R.; Peacock, S.J.; Day, N.P.J.; Visai, L.; Speziale, P.; Lam, A.; Pallen, M.; Foster, T.J. Characterization of novel LPXTG-containing proteins of Staphylococcus aureus identified from genome sequences. Microbiology 2003, 149, 643–654. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Becker, K.; Friedrich, A.W.; Lubritz, G.; Weilert, M.; Peters, G.; Eiff, C. Von Prevalence of Genes Encoding Pyrogenic Toxin Superantigens and Exfoliative Toxins among Strains of Staphylococcus aureus isolated from blood and nasal specimens. J. Clin. Microbiol. 2003, 41, 1434–1439. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Lina, G.; Boutite, F.; Tristan, A.; Bes, M.; Etienne, J.; Vandenesch, F. Bacterial competition for human nasal cavity colonization: Role of Staphylococcal agr alleles. Appl. Environ. Microbiol. 2003, 69, 18–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Hussain, M.; Von Eiff, C.; Sinha, B.; Joost, I.; Herrmann, M.; Peters, G.; Becker, K. eap gene as novel target for specific identification of Staphylococcus aureus. J. Clin. Microbiol. 2008, 46, 470–476. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Hussain, M.; Becker, K.; Von Eiff, C.; Schrenzel, J.; Peters, G.; Herrmann, M. Identification and characterization of a novel 38.5-Kilodalton cell surface protein of Staphylococcus aureus with extended-spectrum binding activity for extracellular matrix and plasma proteins. J. Bacteriol. 2001, 183, 6778–6786. [Google Scholar] [CrossRef] [Green Version]
  68. Fuchs, H.J.; Borowitz, D.S.; Christiansen, D.H.; Morris, E.M.; Nash, M.L.; Ramsey, B.W.; Rosenstein, B.J.; Smith, A.L.; Wohl, M.E. Effect of Aerosolized Recombinant Human DNase on Exacerbations of Respiratory Symptoms and on Pulmonary Function in Patients with Cystic Fibrosis. N. Engl. J. Med. 1994, 331, 637–642. [Google Scholar] [CrossRef]
  69. Quanjer, P.H.; Stanojevic, S.; Cole, T.J.; Baur, X.; Hall, G.L.; Culver, B.H.; Enright, P.L.; Hankinson, J.L.; Ip, M.S.M.; Zheng, J.; et al. Multi-ethnic reference values for spirometry for the 3–95-yr age range: The global lung function 2012 equations. Eur. Respir. J. 2012, 40, 1324–1343. [Google Scholar] [CrossRef]
  70. Recsei, P.; Kreiswirth, B.; O’Reilly, M.; Schlievert, P.; Gruss, A.; Novick, R.P. Regulation of exoprotein gene expression in Staphylococcus aureus by agr. MGG Mol. Gen. Genet. 1986, 202, 58–61. [Google Scholar] [CrossRef]
  71. Yarwood, J.M.; Schlievert, P.M. Quorum sensing in Staphylococcus infections. J. Clin. Invest. 2003, 112, 1620–1625. [Google Scholar] [CrossRef]
Figure 1. Time-dependent quantity of virulence-genes. The red line indicates the averaged loss of virulence genes over time for the study population. The grey dots indicate the virulence gene quantity of the individual S. aureus clones at different dates during the study. Significance levels are indicated by *, **, and ***.
Figure 1. Time-dependent quantity of virulence-genes. The red line indicates the averaged loss of virulence genes over time for the study population. The grey dots indicate the virulence gene quantity of the individual S. aureus clones at different dates during the study. Significance levels are indicated by *, **, and ***.
Toxins 12 00279 g001
Figure 2. Prevalence of virulence genes and agr-types in the investigated CF S. aureus isolates. (a) Prevalence of cap, chip, and adhesin-genes in exacerbation vs. no-exacerbation patients; patients without exacerbation are indicated by blue bars, and patients with exacerbation are indicated by red bars. (b) Prevalence of toxin-genes in exacerbation vs. no-exacerbation patients; (c) prevalence of agr-types in exacerbation vs. no-exacerbation patients (Mann–Whitney U test); and (d) prevalence of the agr-type 1 among subjects without exacerbation compared to exacerbation patients who also received intensive antibiotic treatment (Mann–Whitney U test). Significance levels are indicated by *, **, and ***.
Figure 2. Prevalence of virulence genes and agr-types in the investigated CF S. aureus isolates. (a) Prevalence of cap, chip, and adhesin-genes in exacerbation vs. no-exacerbation patients; patients without exacerbation are indicated by blue bars, and patients with exacerbation are indicated by red bars. (b) Prevalence of toxin-genes in exacerbation vs. no-exacerbation patients; (c) prevalence of agr-types in exacerbation vs. no-exacerbation patients (Mann–Whitney U test); and (d) prevalence of the agr-type 1 among subjects without exacerbation compared to exacerbation patients who also received intensive antibiotic treatment (Mann–Whitney U test). Significance levels are indicated by *, **, and ***.
Toxins 12 00279 g002
Figure 3. Association of agr-types and virulence genes or spa-types. (a) Association between agr-types and virulence gene prevalence counts (Kruskal–Wallis test); the 4 agr-types are indicated as follows: 1 = red, 2 = green, 3 = blue, and 4 = purple. (b) Association between agr-types and prevalent spa-types (Chi-squared test). Significance levels are indicated by *, **, and ***.
Figure 3. Association of agr-types and virulence genes or spa-types. (a) Association between agr-types and virulence gene prevalence counts (Kruskal–Wallis test); the 4 agr-types are indicated as follows: 1 = red, 2 = green, 3 = blue, and 4 = purple. (b) Association between agr-types and prevalent spa-types (Chi-squared test). Significance levels are indicated by *, **, and ***.
Toxins 12 00279 g003
Figure 4. Longitudinal development of virulence gene counts depending on agr-type (Wald Test). Colored lines indicate the averaged loss of virulence genes over time for patient groups showing a predominance of an agr-type (blue = agr-type 1, red = agr-type 2, green = agr type 3, and yellow = agr-type 4). The colored dots indicate the virulence gene quantity of the individual S. aureus clones at different dates during the study. Significance levels are indicated by *, **, and ***.
Figure 4. Longitudinal development of virulence gene counts depending on agr-type (Wald Test). Colored lines indicate the averaged loss of virulence genes over time for patient groups showing a predominance of an agr-type (blue = agr-type 1, red = agr-type 2, green = agr type 3, and yellow = agr-type 4). The colored dots indicate the virulence gene quantity of the individual S. aureus clones at different dates during the study. Significance levels are indicated by *, **, and ***.
Toxins 12 00279 g004
Figure 5. Association between the different agr-types and the four virulence gene clusters (Chi-squared test; red = cluster 1, green = cluster 2, blue = cluster 3, and purple = cluster 4). Significance levels are indicated by *, **, and ***.
Figure 5. Association between the different agr-types and the four virulence gene clusters (Chi-squared test; red = cluster 1, green = cluster 2, blue = cluster 3, and purple = cluster 4). Significance levels are indicated by *, **, and ***.
Toxins 12 00279 g005
Figure 6. Plasticity of virulence genes depending on the presence of environmental stress (Mann–Whitney U test). Events of loss and gain of virulence genes in isolates exposed to environmental stress are indicated in red, and genetical switch events occurring in isolates residing in less-challenging environments are indicated in blue. Significance levels are indicated by *, **, and ***.
Figure 6. Plasticity of virulence genes depending on the presence of environmental stress (Mann–Whitney U test). Events of loss and gain of virulence genes in isolates exposed to environmental stress are indicated in red, and genetical switch events occurring in isolates residing in less-challenging environments are indicated in blue. Significance levels are indicated by *, **, and ***.
Toxins 12 00279 g006
Table 1. Significant IgG levels against S. aureus antigens in patients vs. healthy controls and estimated effect on forced expiratory volume (1 s) FEV1% a.
Table 1. Significant IgG levels against S. aureus antigens in patients vs. healthy controls and estimated effect on forced expiratory volume (1 s) FEV1% a.
AntigenMean IgG Level All Patients (±SE) bMean IgG Level Controls (±SE) bp Value c,dEstimated Effect on FEV1% ep Value d
CHIPS11492 (±220) 11019 (±349)0.104−0.000730.0185
ClfB4587 (±222)4092 (±338)0.38090.002079<0.0001
Eta3999 (±447)2043 (±424)0.1180−0.000610.0033
Etb613 (±113)320 (±94)0.0032−0.00700<0.0001
HlgB13621 (±174)9878 (±407)<0.0001−0.001090.0088
LukF4079 (±133)2718 (±243)<0.0001−0.00305<0.0001
LukS14097 (±133) 7134 (±524)<0.0001−0.00169<0.0001
SdrE3139 (±201)1992 (±221)0.0253−0.00179<0.0001
SED1106 (±118)1292 (±268)0.04980.0014900.0377
SEH2290 (±291)2174 (±359)0.0113−0.000470.0986
a Data are partially retrieved from our earlier study [26]; b significant difference of results between 182 patients and 53 healthy controls; c p values of differences between patient and controls groups (Mann–Whitney U test); d adjusted p-values (Bonferroni correction); and e IgG levels are modelled as continuous factors. Estimated effects are therefore interpreted as mean change in FEV1% predicted.

Share and Cite

MDPI and ACS Style

Lange, J.; Heidenreich, K.; Higelin, K.; Dyck, K.; Marx, V.; Reichel, C.; Wamel, W.v.; Reijer, M.d.; Görlich, D.; Kahl, B.C. Staphylococcus aureus Pathogenicity in Cystic Fibrosis Patients—Results from an Observational Prospective Multicenter Study Concerning Virulence Genes, Phylogeny, and Gene Plasticity. Toxins 2020, 12, 279. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12050279

AMA Style

Lange J, Heidenreich K, Higelin K, Dyck K, Marx V, Reichel C, Wamel Wv, Reijer Md, Görlich D, Kahl BC. Staphylococcus aureus Pathogenicity in Cystic Fibrosis Patients—Results from an Observational Prospective Multicenter Study Concerning Virulence Genes, Phylogeny, and Gene Plasticity. Toxins. 2020; 12(5):279. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12050279

Chicago/Turabian Style

Lange, Jonas, Kathrin Heidenreich, Katharina Higelin, Kristina Dyck, Vanessa Marx, Christian Reichel, Willem van Wamel, Martijn den Reijer, Dennis Görlich, and Barbara C. Kahl. 2020. "Staphylococcus aureus Pathogenicity in Cystic Fibrosis Patients—Results from an Observational Prospective Multicenter Study Concerning Virulence Genes, Phylogeny, and Gene Plasticity" Toxins 12, no. 5: 279. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12050279

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