Stable antibiotic resistance and rapid human adaptation in livestock-associated MRSA
Abstract
Mobile genetic elements (MGEs) are agents of horizontal gene transfer in bacteria, but can also be vertically inherited by daughter cells. Establishing the dynamics that led to contemporary patterns of MGEs in bacterial genomes is central to predicting the emergence and evolution of novel and resistant pathogens. Methicillin-resistant Staphylococcus aureus (MRSA) clonal-complex (CC) 398 is the dominant MRSA in European livestock and a growing cause of human infections. Previous studies have identified three categories of MGEs whose presence or absence distinguishes livestock-associated CC398 from a closely related and less antibiotic-resistant human-associated population. Here, we fully characterise the evolutionary dynamics of these MGEs using a collection of 1180 CC398 genomes, sampled from livestock and humans, over 27 years. We find that the emergence of livestock-associated CC398 coincided with the acquisition of a Tn916 transposon carrying a tetracycline resistance gene, which has been stably inherited for 57 years. This was followed by the acquisition of a type V SCCmec that carries methicillin, tetracycline, and heavy metal resistance genes, which has been maintained for 35 years, with occasional truncations and replacements with type IV SCCmec. In contrast, a class of prophages that carry a human immune evasion gene cluster and that are largely absent from livestock-associated CC398 have been repeatedly gained and lost in both human- and livestock-associated CC398. These contrasting dynamics mean that when livestock-associated MRSA is transmitted to humans, adaptation to the human host outpaces loss of antibiotic resistance. In addition, the stable inheritance of resistance-associated MGEs suggests that the impact of ongoing reductions in antibiotic and zinc oxide use in European farms on livestock-associated MRSA will be slow to be realised.
Editor's evaluation
The reviewers recognised the importance of understanding where new strains of microbes come from and how they change over time for infection control and prevention. Staphylococcus aureus CC398 is an important strain that 'spills over' from livestock to humans, carrying with it high levels of resistance to antibiotics commonly used in farming. This paper compares more than 1000 genomes of CC398 and concludes that spillover is likely to carry resistance to tetracyclines and other antibiotics into humans that will persist over time.
https://doi.org/10.7554/eLife.74819.sa0eLife digest
Antibiotic-resistant infections are a growing threat to human health. In 2019, these hard-to-treat infections resulted in 4.95 million deaths making them the third leading cause of death that year. Excessive use of antibiotics in humans is likely driving the emergence of drug-resistant bacteria. But there is a concern that use of antibiotics on livestock farms is also contributing. A type of bacteria traced back to livestock is a growing cause of human infections that do not respond to treatment with the antibiotic methicillin in Europe. It is called livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA).
Bacteria can share genes that make them drug resistant or more deadly. These genes are often carried on mobile genetic elements that promote their movement from one bacterial cell to another. The most common type of LA-MRSA in Europe is clonal-complex 398 (CC398). It has two mobile genetic elements carrying antibiotic-resistance genes, but generally lacks a mobile genetic element that helps the bacterium escape the human immune system. Learning more about how LA-MRSA acquired these genetic changes may help scientists develop better strategies to protect the public.
Matuszewska, Murray et al. analyzed the genomes of more than 1,000 samples of CC398 collected from humans, pigs and 13 other animal species in 28 countries over 27 years. They used this data to reconstruct the bacteria’s evolutionary history. Matuszewska, Murray et al. show that two mobile elements containing antibiotic resistance genes in CC398 were gained decades ago. One is more than 50 years old and was likely acquired around the time antibiotic use in livestock became common. While most CC398 in livestock do not have a mobile element that helps LA-MRSA evade the human immune system, they often gain it when they infect humans. This leads to highly drug-resistant human MRSA infections.
The results of this study suggest that LA-MRSA is a serious threat to human health. The resistance of this bacterium has persisted for decades, spreading across different livestock species and different countries. These drug-resistant bacteria in livestock readily infect humans. Current efforts to reduce antibiotic use in farms may take decades to mitigate these risks. Additionally, the ban on zinc-oxide use on livestock in the European Union (coming into force June 2022) may not help reduce LA-MRSA, because the genes conferring resistance to bacteria and zinc treatment are not always linked.
Introduction
Mobile genetic elements (MGEs) play an important role in the evolution of bacterial pathogens. They can move rapidly between bacterial genomes, but can also be vertically inherited through stable integration into a host genome. As MGEs often carry genes associated with virulence and antibiotic resistance (Frost et al., 2005; Rankin et al., 2011), an understanding of the drivers and barriers to their acquisition and maintenance is central to predicting the emergence and evolution of novel and resistant pathogens (Brockhurst et al., 2019).
The emergence and evolution of methicillin-resistant Staphylococcus aureus (MRSA) across different ecological niches and host species are associated with the horizontal transfer of MGEs. Methicillin resistance is carried by the staphylococcal cassette chromosome element SCCmec, and additional MGEs carry resistance to other antibiotics, virulence factors, and host-specific adaptations (Hanssen and Ericson Sollid, 2006; Jamrozy et al., 2017; Haag et al., 2019; Turner et al., 2019; Matuszewska et al., 2020). While most MRSA clonal complexes (CCs) show an association with specific MGEs, their dynamics are not widely understood, leading to a gap in our understanding of the adaptive potential of S. aureus CCs.
Intensification combined with high levels of antibiotic use in farming has led to particular concerns about livestock as reservoirs of antibiotic-resistant human infections (World Health Organization, 2015). CC398 has become the dominant MRSA in European livestock. Its rise has been particularly evident in Danish pig farms where the proportion of MRSA-positive herds has increased from <5% in 2008 to 90% in 2018 (Sieber et al., 2018; DANMAP, 2019), but it has also been observed in other European countries, and other livestock species (Lekkerkerk et al., 2015; Islam et al., 2017; Anjum et al., 2019). Livestock-associated (LA) MRSA CC398 has been associated with increasing numbers of human infections, in both people with and without direct contact with livestock (Larsen et al., 2017; van Alen et al., 2017; Sieber et al., 2019). Understanding the emergence and success of CC398 in European livestock and its capacity to infect the human host is integral to managing the risk that it, and other livestock-associated pathogens, pose to public health.
Previous studies have used genome sequences to reconstruct the evolutionary history of CC398 (Price et al., 2012; Ward et al., 2014; Gonçalves da Silva et al., 2017). They identified a largely methicillin-resistant and livestock-associated clade of CC398, that falls as either sister to Ward et al., 2014 or within (Price et al., 2012; Gonçalves da Silva et al., 2017) a largely methicillin-sensitive and human-associated clade. Through comparing the genomes of isolates from livestock- and human-associated CC398, these studies concluded that the emergence of CC398 in livestock was associated with both the acquisition of antibiotic resistance genes and the loss of genes associated with human immune evasion. While these genes are known to be carried on three categories of MGEs (Tn916 conjugative transposons, SCCmec, and φSa3 prophages), little is known about the dynamics of these MGEs within CC398. Here, we undertake a comprehensive reconstruction of the evolutionary dynamics of these MGEs. We find that while their patterns of presence/absence all show a strong association with the transition to livestock, this is the result of contrasting dynamics. These dynamics can inform predictions about the risk posed by LA-MRSA spillover infections in humans, and the resilience of antibiotic resistance in LA-MRSA to ongoing changes in antibiotic use in farming.
Results
Livestock-associated CC398 emerged between 1957 and 1970
We collected and assembled publicly available whole-genome sequence data from CC398, and sequenced five isolates recently sampled from pig farms in the United Kingdom. Our collection includes high-quality whole genome assemblies of 1180 isolates (including 43 complete reference genomes). This collection spans 15 host species (including humans, pigs, cows, chickens, turkeys, and horses), 28 countries (across Europe, America, Asia, and Australasia), and 27 years (1992–2018) (Figure 1—figure supplement 1 and Supplementary file 1).
We constructed a recombination-stripped maximum likelihood phylogeny of CC398 using reference-mapped assemblies of our collection. We rooted the phylogeny with outgroups from four other S. aureus sequence types (STs) in four separate reconstructions. We also constructed a phylogeny from a recombination-stripped concatenated alignment of core genes extracted from de novo assemblies, with a midpoint rooting. These reconstructions consistently returned the same topology, and one that is described in previous studies: a livestock-associated clade of CC398 (704 isolates) that falls within a more diverse, and other largely human-associated clade (476 isolates) (Figure 1a and Figure 1—figure supplement 2; Price et al., 2012; Gonçalves da Silva et al., 2017).
We used the temporal structure in our collection to date the origin of the livestock-associated clade. Due to the size of our collection, we constructed dated phylogenies from three subsampled datasets, each of which includes 250 isolates. Our estimates of the evolutionary rate were consistent across all reconstructions (1.1–1.6 × 10−6 subs/site/year), and similar to estimates from previous studies of CC398 (Ward et al., 2014) and other S. aureus CCs (Hsu et al., 2015). This led to an estimate of the origin of the livestock-associated clade of approximately 1964 (95% CI: 1957–1970) (Figure 1b, Figure 1—figure supplement 3, and Figure 1—figure supplement 4).
The transition to livestock association is associated with changes in the frequencies of three MGEs with very different dynamics
Comparisons of the genomes of isolates from human- and livestock-associated CC398 have previously indicated that the transition to livestock was associated with the acquisition of genes associated with both tetracycline and methicillin resistance (tetM and mecA), and the loss of genes associated with human immune evasion (the immune evasion gene cluster) (Price et al., 2012). Our analyses of this larger collection are broadly consistent with this. We find that the genes whose presence most strongly distinguishes isolates from the human- and livestock-associated groups are associated with three categories of MGEs: (1) a Tn916 transposon carrying tetM, (2) SCCmec carrying mecA, and (3) φSa3 prophages carrying a human immune evasion gene cluster (Figure 1a, Figure 1—figure supplement 5, and Supplementary file 2). Genes associated with the Tn916 transposon and SCCmec elements are more common in the livestock-associated clade, while the reverse is true of genes associated with φSa3 prophages.
Stable maintenance of a Tn916 transposon carrying tetM
We identified a contiguous assembly of a Tn916 transposon carrying tetM in 699/704 isolates in our collection of livestock-associated CC398 (Figure 2a, Supplementary files 3 and 4; de Vries et al., 2009; Roberts and Mullany, 2009). Several lines of evidence indicate that the presence of the Tn916 transposon in livestock-associated CC398 is the result of a single acquisition event, followed by stable inheritance. First, the location of the transposon in the genome of livestock-associated CC398 is conserved (Tn916 is always found next to the same core gene; WP_000902814 in the published annotation of S0385). Second, an alignment of the coding regions of this element (extracted from de novo assemblies) shows a similar average nucleotide diversity to core genes (Figure 2d). Third, a phylogeny constructed from the genes in this element is entirely congruent with the phylogeny of livestock-associated CC398 (Figure 2b, c).
Our analyses indicate that the Tn916 transposon has been maintained in livestock-associated CC398 since its origin, and therefore for around 57 years (Figure 1b). Nevertheless, its absence from 5/704 livestock-associated CC398 isolates in our collection suggests that it remains capable of excision (this is consistent with experimental studies that have shown that Tn916 in CC398 is a functional conjugative transposon; de Vries et al., 2009). None of the genes associated with this element are present in these five isolates, and in two we were able to identify an intact integration site in the assembled genome (Figure 2—figure supplement 1). These five isolates are broadly distributed across the livestock-associated clade, and not linked to a particular host species or geographic location.
More variable maintenance of a SCCmec carrying mecA, tetK, and czrC
Previous studies suggest that LA-MRSA CC398 emerged from human-associated methicillin-sensitive S. aureus (MSSA) (Price et al., 2012). However, the presence of SCCmec elements in recently sampled human-associated CC398 isolates that fall basal to the livestock-associated clade, including clinical isolates from China (He et al., 2018; Zou et al., 2022), Denmark (Moller et al., 2019), and New Zealand (Gonçalves da Silva et al., 2017), makes the association between methicillin resistance and livestock association less clear (Figure 1a).
SCCmec elements in S. aureus are categorised into several types (Hanssen and Ericson Sollid, 2006). Consistent with previous studies of CC398 (e.g. Price et al., 2012), we observe both type V (76%) and type IV (21%) in CC398 (we were unable to confidently type the remaining 3%; Figure 3c, Figure 3—figure supplement 1, and Supplementary files 3 and 5). Most of the type V SCCmec elements belong to the subtype Vc previously described in livestock-associated CC398 (Li et al., 2011; Price et al., 2012; Vandendriessche et al., 2014). This element includes two additional resistance genes: tetK (tetracycline resistance) and czrC (heavy metal resistance) (Figure 3a and Supplementary file 6). We identified a full-length version of this element in 335 genomes (including some of the most recent isolates in our collection, sampled from UK pig farms in 2018) and shorter type V elements in 204 genomes. Full-length versions are only observed in the livestock-associated clade, while shorter versions are found in both livestock- and human-associated groups of CC398. Shorter versions often lacked tetK (n = 90) and czrC (n = 117) genes (Figure 3—figure supplement 2). Type IV elements are only observed in the livestock-associated clade. They include subtypes IVa and IVc, and all only carry a single resistance gene: mecA.
While type Vc is the most common SCCmec in our collection of livestock-associated CC398, this largely reflects isolates from pigs. Type Vc SCCmec is much more common in pigs (77% of isolates; n = 286) than type IV SCCmec (7% of isolates). In cows (n = 74), the difference is reduced: 45% are type Vc and 32% are type IV. And in isolates from other animal species (n = 94) type IV elements (55%) are more common than type Vc (33%).
Diversity within the type Vc SCCmec element indicates that a full-length type Vc SCCmec was acquired once by livestock-associated CC398, and has been maintained within CC398 largely through vertical transmission. First, low nucleotide diversity within full-length versions of the element is consistent with 329/335 sharing a recent origin common within livestock-associated CC398 (Figure 3d). Second, patterns of diversity are largely congruent with the core genome phylogeny, consistent with vertical inheritance (Figure 3b, c). Third, low nucleotide diversity within genes shared across full- and shorter-length versions of the element are consistent with most shorter-length versions being the result of deletion within livestock-associated CC398 (Figure 3—figure supplement 2). In contrast, diversity in the type IV elements in livestock-associated CC398 supports four independent acquisitions from outside of CC398 (Figure 3—figure supplement 3). Similar to the type Vc element, once acquired, these elements tend to be maintained.
While the SCCmec elements carried by human-associated CC398 are always type V, 70/80 fall within a single clade from a hospital outbreak in Denmark in 2016 (Moller et al., 2019). They show a truncation relative to the full-length type Vc that is also observed in three isolates from livestock-associated CC398 (leading to the absence of czrC, Figure 3—figure supplement 2). Pairwise nucleotide distances between these 70 human-associated CC398 elements and the full-length type Vc in livestock-associated CC398 are consistent with a recent common ancestor within CC398. In contrast, nucleotide diversity within the other 10 type V SCCmec in human-associated CC398 and distances from the livestock-associated CC398 type Vc indicate multiple independent acquisitions from outside of CC398 (Figure 3—figure supplement 2).
Using our dated phylogeny of CC398 and categorisation of SCCmec based on diversity within the elements, we inferred the dynamics of gain and loss within CC398 (Figure 3e and Figure 3—figure supplement 4). These reconstructions consistently estimated that the type Vc SCCmec had been acquired by livestock-associated CC398 by around 1986 (95% CI: 1982–1990), and has therefore been maintained within livestock-associated CC398 for around 35 years. While the diversity within this element indicates a single acquisition by CC398, these reconstructions indicate multiple gains, likely reflecting horizontal transmission within CC398. They also indicate that the acquisition by livestock-associated CC398 predated the acquisition by human-associated CC398, consistent with transmission from livestock-associated CC398 to human-associated CC398. In contrast, type IV elements show evidence of several more recent acquisitions in relatively quick succession between 1997 and 2004 (Figure 3e and Figure 3—figure supplement 4).
Together, our analyses are consistent with LA-MRSA emerging from human-associated MSSA and acquiring the type Vc SCCmec element following its initial diversification. The element has been stably maintained (particularly in pigs), although with several exceptions – including deletions of parts of the element and replacements with smaller type IV SCCmec.
Loss of a φSa3 prophage carrying a human immune evasion gene cluster, but with frequent reacquisition
In contrast to the maintenance of Tn916 and SCCmec type V in the livestock-associated clade, the φSa3 prophage is highly dynamic. φSa3 prophages carry human immune evasion gene clusters that include a variable set of functional genes that encode human-specific virulence factors, including sak, chp, scn, sea, and sep (Gladysheva et al., 2003; Postma et al., 2004; Rooijakkers et al., 2005; van Wamel et al., 2006; Thammavongsa et al., 2015; van Alen et al., 2018). They are temperate prophages that primarily integrate into the hlb gene of S. aureus. While φSa3 prophages are present in 88% of the human-associated CC398 isolates in our collection, we find that this is not a consequence of a stable association between CC398 and one prophage. Nucleotide diversity within genes shared across φSa3 prophages (Figure 4 and Supplementary file 3) suggest at least seven (but likely more) acquisitions of an φSa3 prophage into human-associated CC398 from outside of CC398. The set of functional genes carried by these elements indicate that the elements within human-associated CC398 include types C (n = 285), B (n = 111), E (n = 35), A (n = 4), and D (n = 4), and those carried by livestock-associated CC398 isolates include types B (n = 84), E (n = 8), and A (n = 5) (van Wamel et al., 2006; van Alen et al., 2018).
While φSa3 prophages are rare in livestock-associated CC398 (68/704 isolates), diversity within these elements indicate at least 15 (but likely more) acquisitions of these MGEs by livestock-associated CC398. The majority of these elements (69%) do not share a recent common ancestor with those in human-associated CC398. φSa3 prophages present in livestock-associated CC398 generally show evidence of recent acquisition, with a notable exception being a previously described poultry-associated subclade (n = 51) (Price et al., 2012; Larsen et al., 2016b; Pérez-Moreno et al., 2017; Tang et al., 2017). In our collection, 39% of the isolates in this clade are from poultry (compared to 4% across the rest of the livestock-associated clade). Low nucleotide diversity within the type B φSa3 consistently present in isolates within this clade suggests that it has been maintained since its acquisition approximately 21 years ago (95% CI: 1997–2001; Figure 1b and Figure 4—figure supplement 1).
Distinct route to multi-drug resistance in livestock-associated CC398
Livestock-associated CC398 is frequently multi-drug resistant. We see evidence of this in our dataset where 81% of livestock-associated CC398 isolates carry one or more resistance genes for antibiotic classes other than tetracyclines and β-lactams (Figure 5—figure supplement 1, Figure 5—figure supplement 2, and Supplementary file 7). Sixty-seven percent of livestock-associated CC398 isolates have genes associated with trimethoprim resistance (dfrA, dfrK, or dfrG), 42% have genes associated with macrolide resistance (ermA, ermB, ermC, or ermT), and 26% have genes associated with aminoglycoside resistance (aadA, aphA, or aphD). Not only are resistance genes less common in human-associated CC398 (only 20% of isolates carry tetracycline resistance genes, and 9% carry trimethoprim resistance genes), they also differ in their relative frequencies. In particular, human-associated CC398 isolates more commonly carry genes associated with macrolide resistance (91%).
Spillover of livestock-associated CC398 into humans is associated with the acquisition of φSa3 prophages, but not a loss of resistance genes
φSa3 prophages are more common in human isolates (23%) than in livestock or companion animal isolates (11%) in our collection of livestock-associated CC398. To determine the significance of this association, we identified 70 phylogenetically independent clades that include isolates from both human and livestock or companion animal hosts (Figure 5—figure supplement 3). Comparisons across these groups revealed that isolates from humans are consistently more likely to carry an φSa3 prophage (McNemar’s chi-squared test: p = 1.50 × 10−3; Figure 5).
φSa3 prophages are also less common in isolates from pigs than from other non-human species (only 3% of pig isolates carry one). In 42/70 of our phylogenetically independent groups, the only non-human species was a pig, and none of the pig isolates from these groups carried an φSa3 prophage (while 14% of the human isolates did). In the remaining 28 groups (that included isolates from cows, horses, and poultry), φSa3 prophages were observed in non-human hosts in 17% of groups, but still at a higher frequency in humans (39% of groups) (McNemar’s chi-squared test: p = 0.04).
In contrast, we found no evidence that the spillover of livestock-associated CC398 into humans is associated with the loss (or gain) of individual antibiotic resistance genes (Figure 5). This contrasts with the conclusions of a previous study that suggested that the resistance genes carried by livestock-associated CC398 are likely to be lost in human hosts (Sieber et al., 2019). We find that no resistance gene was significantly more or less common in humans than in livestock species (McNemar’s chi-squared tests: p > 0.1), nor was there a consistent shift in the overall number of resistance genes (there were more resistance genes in isolates from human hosts in 22/70 pairs and more in livestock hosts in 34/70 pairs).
Discussion
We have characterised the evolutionary dynamics of the three classes of MGEs that show the greatest changes in frequency across human- and livestock-associated CC398: the Tn916 transposon and the SCCmec element, which are both common in livestock-associated CC398, and the φSa3 prophage, which is common in human-associated CC398. Despite a consistency in the relative frequencies of these elements across CC398, leading to their strong association with the transition to livestock, these three elements show a broad spectrum of dynamics. The Tn916 transposon carrying tetM shows evidence of stable and consistent vertical transmission in livestock-associated CC398 and absence from human-associated CC398. The type Vc SCCmec, carrying not only mecA, but also tetK and czrC, has also been stably maintained by several lineages of livestock-associated CC398, but by a combination of vertical and horizontal transmission, and with occasional replacement with type IV elements. While type V SCCmec elements are also present in human-associated CC398, there is little evidence of their longer-term maintenance. Finally, while φSa3 prophages carrying a human immune evasion gene cluster are rare in livestock-associated CC398 and common in human-associated CC398, there have been frequent gains and losses in both groups. These contrasting dynamics may reflect variation in the selective benefits and costs of the carriage of these MGEs by CC398, their availability in the environments encountered by CC398, and their mechanisms of transfer.
The three classes of MGEs whose dynamics we describe all employ different mechanisms of horizontal transfer, and this may both influence their intrinsic stability in the CC398 genome and their availability for acquisition from outside of CC398. Experimental studies have found that different types of MGEs vary in their rates of transfer between bacterial cells. In particular, in vitro rates of transfer of phage have been found to be several orders of magnitude higher than of transposons (Humphrey et al., 2021). A 16-day in vivo study of MGE dynamics for two strains of CC398 also observed variation in the mobility of MGEs: the Tn916 transposon, the type V SCCmec and an φSa3 prophage were all stably maintained and not horizontally transferred, while other prophages and plasmids were both gained and lost by daughter cells (McCarthy et al., 2014). While variation in the intrinsic stability of MGEs in their host will influence their long-term dynamics, other factors may dominate. For instance, the long-term stability of Tn916 within CC398 could emerge from a broad spectrum of short-term dynamics. It could result from a very low rate of loss, or frequent loss combined with a high selective cost of loss and a low probability of reacquisition from sources other than very closely related cells. Therefore, the long-term dynamics we describe provide a unique insight into the relationship between these MGEs and CC398.
We find that a Tn916 transposon that carries tetM was acquired once by CC398 and this coincided with its origin in European livestock in 1964 (95% CI: 1957–1970). Our observation of the loss of Tn916 on terminal branches of the phylogeny suggests that while losses do occur, lineages that lose Tn916 are either rapidly outcompeted by those that have maintained it, or they rapidly reacquire it from close relatives. Antibiotics, including tetracyclines, were first licensed for use as growth promoters in livestock in European countries in the 1950s, and were in common use by the end of that decade (Lowbury, 1958). While the use of antibiotics as growth promoters was banned by the European Union in 2006, tetracyclines remain the most commonly used antimicrobial class in livestock farming (World Health Organization, 2015; DANMAP, 2019). Carriage of Tn916 is therefore likely to have been associated with a strong selective benefit for livestock-associated CC398 ever since its emergence. However, as CC398’s exposure to tetracyclines will be intermittent, the long-term stability of Tn916 is also likely to reflect a low selective cost in the absence of treatment or a barrier to reacquisition following loss. Tn916-like elements are found across several bacterial genera (Clewell et al., 1995; Roberts and Mullany, 2009), including other opportunistic pathogens in the respiratory microbiome of pigs (Holden et al., 2009; Hoa et al., 2011) and other S. aureus CCs (de Vries et al., 2009). This suggests that the stability of Tn916 in CC398 is not due to the rarity of the element in the environments encountered by CC398, however it may be due to other barriers to successful transfer.
Previous studies have found evidence that the regulatory system of Tn916 promotes its maintenance in a host through ensuring both that excision only occurs in the presence of tetracycline (or other transcription-limiting cell stress) and that any selective burden in the absence of treatment is minimised (Roberts and Mullany, 2009). Nevertheless, Tn916 has been found to have a much more dynamic association with lineages of other bacterial species (e.g. Streptococcus pneumoniae [D’Aeth et al., 2021]). The fitness costs of MGEs can be host (and even insertion locus) specific and can also be mitigated over time (Starikova et al., 2013; Durão et al., 2018), and therefore the stability of Tn916 in livestock-associated CC398 might reflect a low cost that is specific to this lineage. A combination of high selective benefit, low selective cost and inaccessibility for reacquisition following loss could explain the remarkable stability of Tn916 in livestock-associated CC398, and may in part explain the success of this lineage in livestock.
Our results indicate that the acquisition of an SCCmec occurred later in the expansion of livestock-associated CC398. The most common SCCmec in livestock-associated CC398, type Vc, carries tetK and czrC in addition to mecA. These additional resistance genes might be highly advantageous in livestock, and particularly in pigs. While all livestock-associated CC398 carry tetM, there is evidence that carrying tetK in addition to tetM is associated with increased fitness during exposure to sublethal concentrations of tetracycline (Larsen et al., 2016a). czrC is associated with heavy metal resistance (Cavaco et al., 2010), which is likely to be beneficial in the context of the common supplementation of animal feed with zinc oxide, which in pigs is commonly used to prevent diarrhoea in weaners (Nielsen et al., 2021). Additionally, mecA is likely to be beneficial because of the common use of β-lactams in livestock farming, including third generation cephalosporins (Sjölund et al., 2016; Lekagul et al., 2019). The size of the type Vc element, combined with the replacements and truncations we observe suggests that it may come with a selective cost. In addition, SCCmec type Vc appears to be rare (at least in S. aureus). The element has been found in other staphylococci (S. cohhii in Vervet Monkeys [Hoefer et al., 2021]), but has not been reported in other S. aureus CCs.
Loss of the type Vc SCCmec on internal branches within livestock-associated CC398 is associated with replacement with a type IV SCCmec that only carries mecA. While we find evidence of only a single acquisition of the type Vc SCCmec by livestock-associated CC398, we find evidence of at least four acquisitions of type IV SCCmec. The more recent dates of acquisition, and an apparent association with livestock species other than pigs, might reflect a difference in selective pressures across different livestock species, or a loss of the type Vc element during transmission between livestock populations. The long-term maintenance of the type Vc SCCmec is consistent with a high selective benefit (particularly in pigs where this element is most common). However, its repeated loss and replacement with type IV elements are consistent with a low cost of this replacement, at least in certain contexts (perhaps in other livestock hosts), and might also reflect the rarity of the type V element relative to the type IV.
While the loss of the φSa3 prophage that carries a human immune-evasion gene cluster is associated with the transition to livestock, neither the loss nor the gain of these elements is likely to be a substantial hurdle for the adaptation of CC398 to human or non-human hosts. Their ubiquity in human-associated CC398 and frequent acquisition following transmission of livestock-associated CC398 to humans (consistent with Sieber et al., 2019), is consistent with a strong benefit in the human host (Rohmer and Wolz, 2021). On the other hand, as we find that these elements are frequently lost and generally absent in other host species, these elements may carry a selective burden outside of the human host. The diversity of φSa3 prophages in CC398 suggests they form a large pool of elements within human hosts (van Alen et al., 2018), likely reflecting their ubiquity in human carriage populations of S. aureus (Rohmer and Wolz, 2021) and capacity to transfer between S. aureus lineages.
We observe one clear exception to the pattern of recent acquisitions of φSa3 in livestock-associated CC398 in response to spillover events: the acquisition of an φSa3 prophage at the base of a poultry-associated subclade of livestock-associated CC398 (Price et al., 2012; Larsen et al., 2016b; Pérez-Moreno et al., 2017; Tang et al., 2017). This clade was first described as a hybrid LA-MRSA CC9/CC398 lineage (Price et al., 2012), and was subsequently investigated as a lineage associated with human disease (Larsen et al., 2016b). The maintenance of an φSa3 in this lineage may reflect more frequent transmission via a human host, or an adaptation to poultry (a recent study suggested that φSa3 prophages might aid immune evasion in species other than humans [Jung et al., 2017]). Either way, this might make this lineage a greater immediate threat to public health.
While livestock-associated CC398 is found across a broad range of livestock species, it is most commonly associated with pigs. The dynamics of the SCCmec and φSa3 prophages that we have identified are both consistent with livestock-associated CC398 originating in pig farms and later spreading to other livestock species. The lower frequency of the type Vc SCCmec in species other than pigs could reflect random loss during transmission bottlenecks, or a reduced benefit of either tetK or czrC in these species. Similarly, the higher frequency of φSa3 prophages in other species might reflect an increased benefit or reduced cost, or more recent transmission via human hosts.
Our results reveal that LA-MRSA CC398 is a stably antibiotic-resistant pathogen that is capable of dynamic readaptation to humans. We find that Tn916 and SCCmec are both stably maintained in livestock-associated CC398, across different livestock species and countries, and that neither of these MGEs (or other antibiotic resistance genes) tend to be lost when livestock-associated CC398 is transmitted to humans. This suggests that these MGEs are associated with a low selective cost in both livestock and human hosts, and therefore across variable levels and types of antibiotic exposure. The stability of these two MGEs, combined with the capacity of livestock-associated CC398 to rapidly acquire the φSa3 prophage, underlines the threat posed by LA-MRSA to public health. While SCCmec is less stably maintained than Tn916, our identification of several independent acquisitions of type IV SCCmec suggests both a strong selective benefit, not contingent on the carriage of tetK and czrC, and an availability of type IV elements. These dynamics predict that the impact of gradual reductions in antibiotic consumption on LA-MRSA is likely to be slow to be realised and that the forthcoming EU ban on medical zinc supplementation in pig feed may have a limited impact on LA-MRSA (European Medicines Agency, 2017; European Medicines Agency, 2020 DANMAP, 2019). Further work is, however, required to understand the factors that underlie the acquisition and maintenance of resistance genes within LA-MRSA, and how they differ from human-associated MRSA lineages.
Materials and methods
Data collection
Request a detailed protocolAll available genome sequence data relating to S. aureus CC398 were downloaded from public databases (https://www.ncbi.nlm.nih.gov/sra and https://www.ebi.ac.uk/ena/browser/; accessed 2019), with metadata in some cases obtained by request (Supplementary file 1). We additionally sequenced five isolates sampled from UK pig farms in 2018. All publicly available complete S. aureus genomes assemblies (https://www.ncbi.nlm.nih.gov/; accessed 2019) were MLST typed using Pathogenwatch (https://pathogen.watch), and the 43 genomes identified as CC398 were added to our collection. After exclusion of low-quality assemblies and isolates mis-characterised as CC398, this led to a collection of 1,180 genomes.
Genomic library preparation and sequencing
Request a detailed protocolFor sequencing of the five UK pig farm isolates, genomic DNA was extracted from overnight cultures grown in TSB at 37°C with 200 rpm shaking using the MasterPure Gram Positive DNA Purification Kit (Cambio, UK). Illumina library preparation and Hi-Seq sequencing were carried out as previously described (Harrison et al., 2013).
Genome assembly
Request a detailed protocolWe used sequence data from all isolates to generate de novo assemblies with Spades v.3.12.0 (Bankevich et al., 2012). We removed adapters and low-quality reads with Cutadapt v1.16 (Martin, 2011) and Sickle v1.33 (Joshi and Fass, 2011), and screened for contamination using FastQ Screen (Wingett and Andrews, 2018). Optimal k-mers were identified based on average read lengths for each genome. All assemblies were evaluated using QUAST v.5.0.1 (Gurevich et al., 2013) and we mapped reads back to de novo assemblies to investigate polymorphism (indicative of mixed cultures) using Bowtie2 v1.2.2 (Langmead and Salzberg, 2012). Low-quality genome assemblies were excluded from further analysis (i.e. N50 < 10,000, contigs smaller than 1 kb contributing to >15% of the total assembly length, total assembly length outside of the median sequence length ± one standard deviation, or >1500 polymorphic sites). We identified genomes mischaracterised as CC398 via two approaches and excluded them from further analysis. First, we identified STs with MLST-check (Page et al., 2016) and grouped into CCs using the eBURST algorithm with a single locus variant (Francisco et al., 2009). Second, we constructed a neighbour-joining tree based on a concatenated alignment of MLST genes (arcC, aroE, glpF, gmk, pta, tpi, and yqiL) for our collection and 13 additional reference genomes from other CCs, using the ape package in R and a K80 substitution model (Paradis et al., 2004).
We generated reference-mapped assemblies with Bowtie2 using the reference genome S0385 (GenBank accession no. AM990992). For reference genomes, we generated artificial FASTQ files with ArtificialFastqGenerator (Frampton and Houlston, 2012). Average coverage and number of missing sites in these assemblies were used as an additional quality control measure; genomes with average coverage <×50 or with >10% missing sites were excluded.
We identified recombination in the reference-mapped alignment using both Gubbins v2.3.1 (Croucher et al., 2015) and ClonalFrame (Didelot and Wilson, 2015). We masked all the recombinant sites identified from our alignment. We additionally masked a region of ∼123 kb that was identified as horizontally acquired from an ST9 donor in a previous study (Price et al., 2012).
Genome annotation and identification of homologous genes
Request a detailed protocolWe annotated de novo assemblies with Prokka v2.8.2 (Seemann, 2014) and identified orthologous genes with Roary (Page et al., 2015) using recommended parameter values. We created a core gene alignment with Roary and identified recombinant sites using Gubbins. We identified antibiotic resistance genes using the Pathogenwatch AMR prediction module (Wellcome Sanger Institute), which uses BLASTn (Altschul et al., 1990) with a cut-off of 75% coverage and 80–90% identity threshold (depending on the gene) against a S. aureus AMR database.
Phylogenetic analyses
Request a detailed protocolWe carried out phylogenetic reconstruction for the reference-mapped alignment with RAxML v8.2.4 using the GTR+Γ model and 1000 bootstraps (Stamatakis, 2014). Sites where >0.1% of genomes showed evidence of recombination or had missing data were excluded from the analysis. We constructed dated phylogenies using BEAST v1.10 with a HKY+Γ model, a strict molecular clock, and constant population size coalescent prior, from the coding regions of the reference-mapped alignment (Drummond et al., 2012). We fit separate substitution models and molecular clocks to first/second and third codon positions to reflect differences in selective constraint. We constructed phylogenies for three random subsamples of 250 isolates (200 from livestock-associated CC398 and 50 from human-associated CC398). Subsamples were non-overlapping except for 30 genomes representing the most divergent lineages within the livestock-associated clade, to ensure a consistent description of the origin of this clade. We constructed additional phylogenies from subsamples that included only isolates from the LA clade to establish that consistent rate estimates were returned over different evolutionary depths. We investigated temporal signal in each dataset through a regression of root-to-tip distance against sampling date, and a permutation test of dates over tips (with clustering used to correct for any confounding of temporal and genetic structures) (Murray et al., 2016; Rambaut et al., 2016). Trees are visualised and annotated using ITOL (Letunic and Bork, 2021).
Comparative analyses of MGEs
Request a detailed protocolGenes associated with the transition to livestock were identified through comparing frequencies of homologous genes identified with Roary across human- and livestock-associated CC398. We investigated the association between these genes and MGEs through analysis of (1) physical locations within our de novo assemblies and the reference genomes, (2) correlations in their presence/absence across our collection, and (3) comparison with descriptions in the literature and in online databases of particular MGEs.
We confirmed the identity of the Tn916 transposon by comparison with published descriptions and publicly available annotated sequences. SCCmec elements were initially categorised into types using a BLASTn search of all representative SCCmec types from the SCCmecFinder database in our de novo assemblies (Supplementary file 5; Hanssen and Ericson Sollid, 2006). φSa3 prophages were initially identified and categorised through identification of functional genes associated with the human immune evasion gene cluster. To further characterise diversity within these elements we identified genes associated with each element across our collection. We constructed alignments of genes within MGEs using Clustal Omega v1.2.3 (Sievers and Higgins, 2018) and checked for misalignment by eye.
We analysed variation in both gene content and nucleotide diversity within shared genes for each MGE. We estimated pairwise nucleotide distances in concatenated alignments of shared genes for each type of element using the ape package in R (Paradis et al., 2004), and constructed maximum likelihood trees using RaxML and minimum-spanning trees using GrapeTree (Zhou et al., 2018) to investigate co-phylogeny. We generated confidence intervals for our estimates of mean pairwise nucleotide distances within MGEs by re-estimating the mean distance from 1000 bootstrapped samples of sites. We compared this to sites in the core genome by sampling the same number of sites as were in the MGE alignment from a concatenated alignment of core genes (generated by Roary). For SCCmec we used ancestral state reconstruction in BEAST to infer the evolutionary dynamics of these elements and date their origins within CC398. This involved fitting a discrete traits model to the posterior distributions of trees, with each state representing a version of the element that had been independently acquired by CC398. We used a strict clock model that allowed for asymmetric rates of transitions between states, but we found that the results were robust to use of a symmetric model or a relaxed clock.
Phylogenetically independent groups
Request a detailed protocolTo test the association between spillover into the human host and the presence of φSa3 prophages and antibiotic-resistance genes, we identified 70 phylogenetically independent clades of isolates that were sampled from both human and non-human hosts in the livestock-associated clade. We classed a gene as present in a host if it was observed in any of the isolates from that host in that clade.
Data availability
The Illumina sequences generated in this study have been deposited in the NCBI short read archive under the accession numbers ERR3524650, ERR3524328, ERR3524354, ERR3524446, and ERR3524562. All other sequences used in this study are publicly available and their origins are described in Supplementary file 1.
-
NCBI Sequence Read ArchiveID ERR3524650. Tracking_the_dynamics_of_AMR_genes_within_enteric_bacterial_communities_in_pigs_and_humans.
-
NCBI Sequence Read ArchiveID ERR3524328. Tracking_the_dynamics_of_AMR_genes_within_enteric_bacterial_communities_in_pigs_and_humans.
-
NCBI Sequence Read ArchiveID ERR3524354. Tracking_the_dynamics_of_AMR_genes_within_enteric_bacterial_communities_in_pigs_and_humans.
-
NCBI Sequence Read ArchiveID ERR3524446. Tracking_the_dynamics_of_AMR_genes_within_enteric_bacterial_communities_in_pigs_and_humans.
-
NCBI Sequence Read ArchiveID ERR3524562. Tracking_the_dynamics_of_AMR_genes_within_enteric_bacterial_communities_in_pigs_and_humans.
-
NCBI BioProjectID PRJEB12552. Emergence of CC398 MRSA in New Zealand.
-
NCBI BioProjectID PRJEB18725. Detection and molecular characterisation of Livestock-Associated MRSA in raw meat on retail sale in North West England.
-
NCBI BioProjectID PRJNA347471. Staphylococcus aureus CC398 resequencing data.
-
NCBI BioProjectID PRJEB14187. Studying 3 strains of LA-MRSA from Finland with interesting deletion genotypes.
-
NCBI BioProjectID PRJEB5042. MRSA from rats and humans isolated from the Downtown Eastside of Vancouver, BC, Canada.
-
NCBI BioProjectID PRJEB19362. Prevalence and origin of LA-MRSA CC398 in Danish horses.
-
NCBI BioProjectID PRJNA226567. Staphylococcus aureus.
-
NCBI BioProjectID PRJNA312437. poultry-associated genetic elements in Staphylococcus aureus.
-
NCBI GenBankID PRJNA378150. Staphylococcus aureus strain 08S00974 chromosome, complete genome.
-
NCBI BioProjectID PRJNA508272. Unusual MRSA CC398 hospital outbreak.
-
NCBI BioProjectID PRJEB2655. Diversity_of_MRSA.
-
NCBI BioProjectID PRJEB21015. Tracking_the_dynamics_of_AMR_genes_within_enteric_bacterial_communities_in_pigs_and_humans.
-
NCBI BioProjectID PRJNA274898. Staphylococcus aureus CC398: host adaptation and emergence of methicillin resistance in livestock.
-
NCBI BioProjectID PRJNA430150. Staphylococcus aureus Raw sequence reads.
-
NCBI BioProjectID PRJEB14251. LA-MRSA CC398 from UK animals.
-
NCBI BioProjectID PRJEB25608. Spread of LA-MRSA CC398 in pigs and humans in Denmark.
-
NCBI BioProjectID PRJEB12818. Sequencing of Staphylococcus aureus CC398 from Northern Manhattan.
-
NCBI BioProjectID PRJEB7209. Global transmission and antibiotic resistance dynamics of Staphylococcus aureus CC398 in humans and livestock revealed by whole genome sequence analysis.
-
NCBI BioProjectID PRJEB2755. Staphylococcus_aureus__MSSA__study.
-
NCBI BioProjectID PRJNA660925. The Staphylocossus aureus isolates from central China.
References
-
Basic local alignment search toolJournal of Molecular Biology 215:403–410.https://doi.org/10.1016/S0022-2836(05)80360-2
-
SPAdes: A new genome assembly algorithm and its applications to single-cell sequencingJournal of Computational Biology 19:455–477.https://doi.org/10.1089/cmb.2012.0021
-
The Ecology and Evolution of PangenomesCurrent Biology 29:R1094–R1103.https://doi.org/10.1016/j.cub.2019.08.012
-
Cloning and occurrence of czrC, a gene conferring cadmium and zinc resistance in methicillin-resistant Staphylococcus aureus CC398 isolatesAntimicrobial Agents and Chemotherapy 54:3605–3608.https://doi.org/10.1128/AAC.00058-10
-
Diversity of the tetracycline resistance gene tet(M) and identification of Tn916- and Tn5801-like (Tn6014) transposons in Staphylococcus aureus from humans and animalsThe Journal of Antimicrobial Chemotherapy 64:490–500.https://doi.org/10.1093/jac/dkp214
-
ClonalFrameML: efficient inference of recombination in whole bacterial genomesPLOS Computational Biology 11:e1004041.https://doi.org/10.1371/journal.pcbi.1004041
-
Bayesian phylogenetics with BEAUti and the BEAST 1.7Molecular Biology and Evolution 29:1969–1973.https://doi.org/10.1093/molbev/mss075
-
Evolutionary Mechanisms Shaping the Maintenance of Antibiotic ResistanceTrends in Microbiology 26:677–691.https://doi.org/10.1016/j.tim.2018.01.005
-
Mobile genetic elements: the agents of open source evolutionNature Reviews. Microbiology 3:722–732.https://doi.org/10.1038/nrmicro1235
-
QUAST: quality assessment tool for genome assembliesBioinformatics (Oxford, England) 29:1072–1075.https://doi.org/10.1093/bioinformatics/btt086
-
Staphylococcus aureus in AnimalsMicrobiology Spectrum 7:e0060-2.https://doi.org/10.1128/microbiolspec.GPP3-0060-2019
-
SCCmec in staphylococci: genes on the moveFEMS Immunology and Medical Microbiology 46:8–20.https://doi.org/10.1111/j.1574-695X.2005.00009.x
-
adegenet: a R package for the multivariate analysis of genetic markersBioinformatics (Oxford, England) 24:1403–1405.https://doi.org/10.1093/bioinformatics/btn129
-
SoftwareSickle: A Sliding-Window, Adaptive, Quality-Based Trimming Tool for FastQ FilesSickle.
-
Fast gapped-read alignment with Bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Copresence of tet(K) and tet(M) in Livestock-Associated Methicillin-Resistant Staphylococcus aureus Clonal Complex 398 Is Associated with Increased Fitness during Exposure to Sublethal Concentrations of TetracyclineAntimicrobial Agents and Chemotherapy 60:4401–4403.https://doi.org/10.1128/AAC.00426-16
-
Evidence for Human Adaptation and Foodborne Transmission of Livestock-Associated Methicillin-Resistant Staphylococcus aureusClinical Infectious Diseases 63:1349–1352.https://doi.org/10.1093/cid/ciw532
-
Emergence of Livestock-Associated Methicillin-Resistant Staphylococcus aureus Bloodstream Infections in DenmarkClinical Infectious Diseases 65:1072–1076.https://doi.org/10.1093/cid/cix504
-
Patterns of antibiotic use in global pig production: A systematic reviewVeterinary and Animal Science 7:e100058.https://doi.org/10.1016/j.vas.2019.100058
-
Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotationNucleic Acids Research 49:W293–W296.https://doi.org/10.1093/nar/gkab301
-
Novel types of staphylococcal cassette chromosome mec elements identified in clonal complex 398 methicillin-resistant Staphylococcus aureus strainsAntimicrobial Agents and Chemotherapy 55:3046–3050.https://doi.org/10.1128/AAC.01475-10
-
Symposium on epidemiological risks of antibiotics: hospital infectionProceedings of the Royal Society of Medicine 51:807–810.https://doi.org/10.1177/003591575805101006
-
The Evolutionary Genomics of Host Specificity in Staphylococcus aureusTrends in Microbiology 28:465–477.https://doi.org/10.1016/j.tim.2019.12.007
-
Extensive horizontal gene transfer during Staphylococcus aureus co-colonization in vivoGenome Biology and Evolution 6:2697–2708.https://doi.org/10.1093/gbe/evu214
-
The effect of genetic structure on molecular dating and tests for temporal signalMethods in Ecology and Evolution 7:80–89.https://doi.org/10.1111/2041-210X.12466
-
Antibiotic and medical zinc oxide usage in Danish conventional and welfare-label pig herds in 2016-2018Preventive Veterinary Medicine 189:e105283.https://doi.org/10.1016/j.prevetmed.2021.105283
-
Roary: rapid large-scale prokaryote pan genome analysisBioinformatics (Oxford, England) 31:3691–3693.https://doi.org/10.1093/bioinformatics/btv421
-
APE: Analyses of Phylogenetics and Evolution in R languageBioinformatics (Oxford, England) 20:289–290.https://doi.org/10.1093/bioinformatics/btg412
-
Unusual presence of the immune evasion gene cluster in livestock-associated MRSA of lineage CC398 causing peridural and psoas abscesses in a poultry farmerEnfermedades Infecciosas y Microbiologia Clinica 35:651–654.https://doi.org/10.1016/j.eimc.2016.07.008
-
Chemotaxis inhibitory protein of Staphylococcus aureus binds specifically to the C5a and formylated peptide receptorJournal of Immunology (Baltimore, Md) 172:6994–7001.https://doi.org/10.4049/jimmunol.172.11.6994
-
A modular master on the move: the Tn916 family of mobile genetic elementsTrends in Microbiology 17:251–258.https://doi.org/10.1016/j.tim.2009.03.002
-
The Role of hlb-Converting Bacteriophages in Staphylococcus aureus Host AdaptionMicrobial Physiology 31:109–122.https://doi.org/10.1159/000516645
-
Prokka: rapid prokaryotic genome annotationBioinformatics (Oxford, England) 30:2068–2069.https://doi.org/10.1093/bioinformatics/btu153
-
Clustal Omega for making accurate alignments of many protein sequencesProtein Science 27:135–145.https://doi.org/10.1002/pro.3290
-
RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogeniesBioinformatics (Oxford, England) 30:1312–1313.https://doi.org/10.1093/bioinformatics/btu033
-
Fitness costs of various mobile genetic elements in Enterococcus faecium and Enterococcus faecalisThe Journal of Antimicrobial Chemotherapy 68:2755–2765.https://doi.org/10.1093/jac/dkt270
-
Methicillin-resistant and -susceptible Staphylococcus aureus from retail meat in DenmarkInternational Journal of Food Microbiology 249:72–76.https://doi.org/10.1016/j.ijfoodmicro.2017.03.001
-
Staphylococcal manipulation of host immune responsesNature Reviews. Microbiology 13:529–543.https://doi.org/10.1038/nrmicro3521
-
Methicillin-resistant Staphylococcus aureus: an overview of basic and clinical researchNature Reviews. Microbiology 17:203–218.https://doi.org/10.1038/s41579-018-0147-4
-
High genetic diversity of methicillin-susceptible Staphylococcus aureus (MSSA) from humans and animals on livestock farms and presence of SCCmec remnant DNA in MSSA CC398The Journal of Antimicrobial Chemotherapy 69:355–362.https://doi.org/10.1093/jac/dkt366
-
Time-Scaled Evolutionary Analysis of the Transmission and Antibiotic Resistance Dynamics of Staphylococcus aureus Clonal Complex 398Applied and Environmental Microbiology 80:7275–7282.https://doi.org/10.1128/AEM.01777-14
-
A Survey of Chinese Pig Farms and Human Healthcare Isolates Reveals Separate Human and Animal Methicillin-Resistant Staphylococcus aureus PopulationsAdvanced Science (Weinheim, Baden-Wurttemberg, Germany) 9:e2103388.https://doi.org/10.1002/advs.202103388
Article and author information
Author details
Funding
Wellcome Trust (109385/Z/15/Z)
- Gemma GR Murray
- Lucy A Weinert
Medical Research Council
- Marta Matuszewska
Biotechnology and Biological Sciences Research Council (BB/L018934/1)
- Gemma GR Murray
Newnham College, University of Cambridge
- Gemma GR Murray
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.
Acknowledgements
MM was funded by the Medical Research Council, co-funded by the Raymond and Beverly Sackler Fund. GGRM and LAW were supported by a Sir Henry Dale Fellowship jointly funded by the Wellcome Trust and the Royal Society (109385/Z/15/Z). GGRM was also supported by a ZELS BBSRC award (BB/L018934/1) and a Research Fellowship at Newnham College.
Copyright
© 2022, Matuszewska, Murray et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,847
- views
-
- 827
- downloads
-
- 35
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.