Evolution of H3N2 influenza virus in a guinea pig model

Studies of influenza virus evolution under controlled experimental conditions can provide a better understanding of the consequences of evolutionary processes with and without immunological pressure. Characterization of evolved strains assists in the development of predictive algorithms for both the selection of subtypes represented in the seasonal influenza vaccine and the design of novel immune refocused vaccines. To obtain data on the evolution of influenza in a controlled setting, naïve and immunized Guinea pigs were infected with influenza A/Wyoming/2003 (H3N2). Virus progeny from nasal wash samples were assessed for variation in the dominant and other epitopes by sequencing the hemagglutinin (HA) gene to quantify evolutionary changes. Viral RNA from the nasal washes from infection of naïve and immune animals contained 6% and 24.5% HA variant sequences, respectively. Analysis of mutations relative to antigenic epitopes indicated that adaptive immunity played a key role in virus evolution. HA mutations in immunized animals were associated with loss of glycosylation and changes in charge and hydrophobicity in and near residues within known epitopes. Four regions of HA-1 (75–85, 125–135, 165–170, 225–230) contained residues of highest variability. These sites are adjacent to or within known epitopes and appear to play an important role in antigenic variation. Recognition of the role of these sites during evolution will lead to a better understanding of the nature of evolution which help in the prediction of future strains for selection of seasonal vaccines and the design of novel vaccines intended to stimulated broadened cross-reactive protection to conserved sites outside of dominant epitopes.     

The Wilhelmine E. Key 1999 Invitational lecture. Predicting the evolution of human influenza A

We studied the evolution of the HA1 domain of the H3 hemagglutinin gene from human influenza virus type A. The phylogeny of these genes showed a single dominant lineage persisting over time. We tested the hypothesis that the progenitors of this single evolutionarily successful lineage were viruses carrying mutations at codons at which prior mutations had helped the virus to avoid human immune surveillance. We found evidence that eighteen hemagglutinin codons appeared to have been under positive selection to change the amino acid they encoded in the past. Retrospective tests show that viral lineages undergoing the greatest number of mutations in the positively selected codons were the progenitors of future H3 lineages in nine of eleven recent influenza seasons. Codons under positive selection were associated with antibody combining sites A or B or the sialic acid receptor binding site. However, not all codons in these sites had predictive value. Monitoring new H3 isolates for additional changes in positively selected codons might help identify the most fit extant viral strains that arise during antigenic drift.     

Role of Positive Selection Pressure on the Evolution of H5N1 Hemagglutinin

The surface glycoprotein hemagglutinin (HA) helps the influenza A virus to evade the host immune system by antigenic variation and is a major driving force for viral evolution. In this study, the selection pressure on HA of H5N1 influenza A virus was analyzed using bioinformatics algorithms. Most of the identified positive selection (PS) sites were found to be within or adjacent to epitope sites. Some of the identified PS sites are consistent with previous experimental studies, providing further support to the biological significance of our findings. The highest frequency of PS sites was observed in recent strains isolated during 2005-2007. Phylogenetic analysis was also conducted on HA sequences from various hosts. Viral drift is almost similar in both avian and human species with a progressive trend over the years. Our study reports new mutations in functional regions of HA that might provide markers for vaccine design or can be used to predict isolates of pandemic potential.     

Clonal interference and the periodic selection of new beneficial mutations in Escherichia coli

The conventional model of adaptation in asexual populations implies sequential fixation of new beneficial mutations via rare selective sweeps that purge all variation and preserve the clonal genotype. However, in large populations multiple beneficial mutations may co-occur, causing competition among them, a phenomenon called "clonal interference." Clonal interference is thus expected to lead to longer fixation times and larger fitness effects of mutations that ultimately become fixed, as well as to a genetically more diverse population. Here, we study the significance of clonal interference in populations consisting of mixtures of differently marked wild-type and mutator strains of Escherichia coli that adapt to a minimal-glucose environment for 400 generations. We monitored marker frequencies during evolution and measured the competitive fitness of random clones from each marker state after evolution. The results demonstrate the presence of multiple beneficial mutations in these populations and slower and more erratic invasion of mutants than expected by the conventional model, showing the signature of clonal interference. We found that a consequence of clonal interference is that fitness estimates derived from invasion trajectories were less than half the magnitude of direct estimates from competition experiments, thus revealing fundamental problems with this fitness measure. These results force a reevaluation of the conventional model of periodic selection for asexual microbes.     

Components of Selection in the Evolution of the Influenza Virus: Linkage Effects Beat Inherent Selection

The influenza virus is an important human pathogen, with a rapid rate of evolution in the human population. The rate of homologous recombination within genes of influenza is essentially zero. As such, where two alleles within the same gene are in linkage disequilibrium, interference between alleles will occur, whereby selection acting upon one allele has an influence upon the frequency of the other. We here measured the relative importance of selection and interference effects upon the evolution of influenza. We considered time-resolved allele frequency data from the global evolutionary history of the haemagglutinin gene of human influenza A/H3N2, conducting an in-depth analysis of sequences collected since 1996. Using a model that accounts for selection-caused interference between alleles in linkage disequilibrium, we estimated the inherent selective benefit of individual polymorphisms in the viral population. These inherent selection coefficients were in turn used to calculate the total selective effect of interference acting upon each polymorphism, considering the effect of the initial background upon which a mutation arose, and the subsequent effect of interference from other alleles that were under selection. Viewing events in retrospect, we estimated the influence of each of these components in determining whether a mutant allele eventually fixed or died in the global viral population. Our inherent selection coefficients, when combined across different regions of the protein, were consistent with previous measurements of dN/dS for the same system. Alleles going on to fix in the global population tended to be under more positive selection, to arise on more beneficial backgrounds, and to avoid strong negative interference from other alleles under selection. However, on average, the fate of a polymorphism was determined more by the combined influence of interference effects than by its inherent selection coefficient.     

Epistatically Interacting Substitutions Are Enriched during Adaptive Protein Evolution

Most experimental studies of epistasis in evolution have focused on adaptive changes—but adaptation accounts for only a portion of total evolutionary change. Are the patterns of epistasis during adaptation representative of evolution more broadly? We address this question by examining a pair of protein homologs, of which only one is subject to a well-defined pressure for adaptive change. Specifically, we compare the nucleoproteins from human and swine influenza. Human influenza is under continual selection to evade recognition by acquired immune memory, while swine influenza experiences less such selection due to the fact that pigs are less likely to be infected with influenza repeatedly in a lifetime. Mutations in some types of immune epitopes are therefore much more strongly adaptive to human than swine influenza—here we focus on epitopes targeted by human cytotoxic T lymphocytes. The nucleoproteins of human and swine influenza possess nearly identical numbers of such epitopes. However, mutations in these epitopes are fixed significantly more frequently in human than in swine influenza, presumably because these epitope mutations are adaptive only to human influenza. Experimentally, we find that epistatically constrained mutations are fixed only in the adaptively evolving human influenza lineage, where they occur at sites that are enriched in epitopes. Overall, our results demonstrate that epistatically interacting substitutions are enriched during adaptation, suggesting that the prevalence of epistasis is dependent on the underlying evolutionary forces at play.     

Transmission of clonal hepatitis C virus genomes reveals the dominant but transitory role of CD8⁺ T cells in early viral evolution

The RNA genome of the hepatitis C virus (HCV) diversifies rapidly during the acute phase of infection, but the selective forces that drive this process remain poorly defined. Here we examined whether Darwinian selection pressure imposed by CD8(+) T cells is a dominant force driving early amino acid replacement in HCV viral populations. This question was addressed in two chimpanzees followed for 8 to 10 years after infection with a well-defined inoculum composed of a clonal genotype 1a (isolate H77C) HCV genome. Detailed characterization of CD8(+) T cell responses combined with sequencing of recovered virus at frequent intervals revealed that most acute-phase nonsynonymous mutations were clustered in class I epitopes and appeared much earlier than those in the remainder of the HCV genome. Moreover, the ratio of nonsynonymous to synonymous mutations, a measure of positive selection pressure, was increased 50-fold in class I epitopes compared with the rest of the HCV genome. Finally, some mutation of the clonal H77C genome toward a genotype 1a consensus sequence considered most fit for replication was observed during the acute phase of infection, but the majority of these amino acid substitutions occurred slowly over several years of chronic infection. Together these observations indicate that during acute hepatitis C, virus evolution was driven primarily by positive selection pressure exerted by CD8(+) T cells. This influence of immune pressure on viral evolution appears to subside as chronic infection is established and genetic drift becomes the dominant evolutionary force.     

Analysis of Adaptation Mutants in the Hemagglutinin of the Influenza A(H1N1)pdm09 Virus

Hemagglutinin is the major surface glycoprotein of influenza viruses. It participates in the initial steps of viral infection through receptor binding and membrane fusion events. The influenza pandemic of 2009 provided a unique scenario to study virus evolution. We performed molecular dynamics simulations with four hemagglutinin variants that appeared throughout the 2009 influenza A (H1N1) pandemic. We found that variant 1 (S143G, S185T) likely arose to avoid immune recognition. Variant 2 (A134T), and variant 3 (D222E, P297S) had an increased binding affinity for the receptor. Finally, variant 4 (E374K) altered hemagglutinin stability in the vicinity of the fusion peptide. Variants 1 and 4 have become increasingly predominant, while variants 2 and 3 declined as the pandemic progressed. Our results show some of the different strategies that the influenza virus uses to adapt to the human host and provide an example of how selective pressure drives antigenic drift in viral proteins.     

The genomic rate of molecular adaptation of the human influenza A virus

Quantifying adaptive evolution at the genomic scale is an essential yet challenging aspect of evolutionary biology. Here, we develop a method that extends and generalizes previous approaches to estimate the rate of genomic adaptation in rapidly evolving populations and apply it to a large data set of complete human influenza A virus genome sequences. In accord with previous studies, we observe particularly high rates of adaptive evolution in domain 1 of the viral hemagglutinin (HA1). However, our novel approach also reveals previously unseen adaptation in other viral genes. Notably, we find that the rate of adaptation (per codon per year) is higher in surface residues of the viral neuraminidase than in HA1, indicating strong antibody-mediated selection on the former. We also observed high rates of adaptive evolution in several nonstructural proteins, which may relate to viral evasion of T-cell and innate immune responses. Furthermore, our analysis provides strong quantitative support for the hypothesis that human H1N1 influenza experiences weaker antigenic selection than H3N2. As well as shedding new light on the dynamics and determinants of positive Darwinian selection in influenza viruses, the approach introduced here is applicable to other pathogens for which densely sampled genome sequences are available, and hence is ideally suited to the interpretation of next-generation genome sequencing data.     

Coordinated Evolution of Influenza A Surface Proteins

The surface proteins hemagglutinin (HA) and neuraminidase (NA) of human influenza A virus evolve under selection pressures to escape adaptive immune responses and antiviral drug treatments. In addition to these external selection pressures, some mutations in HA are known to affect the adaptive landscape of NA, and vice versa, because these two proteins are physiologically interlinked. However, the extent to which evolution of one protein affects the evolution of the other one is unknown. Here we develop a novel phylogenetic method for detecting the signatures of such genetic interactions between mutations in different genes – that is, inter-gene epistasis. Using this method, we show that influenza surface proteins evolve in a coordinated way, with mutations in HA affecting subsequent spread of mutations in NA and vice versa, at many sites. Of particular interest is our finding that the oseltamivir-resistance mutations in NA in subtype H1N1 were likely facilitated by prior mutations in HA. Our results illustrate that the adaptive landscape of a viral protein is remarkably sensitive to its genomic context and, more generally, that the evolution of any single protein must be understood within the context of the entire evolving genome.     

Analysis of antigenically important residues in human influenza A virus in terms of B-cell epitopes

In this paper we undertake an analysis of the antigenicity of influenza A virus hemagglutinin. We developed a novel computational approach to the identification of antigenically active regions and showed that the amino acid substitutions between successive predominant seasonal strains form clusters that are consistent, in terms of both their location and their size, with the properties of B-cell epitopes in general and with those epitopes that have been identified experimentally in influenza A virus hemagglutinin to date. Such an interpretation provides a biologically plausible framework for an understanding of the location of antigenically important substitutions that is more specific than the canonical "antigenic site" model and provides an effective basis for deriving models that predict antigenic escape in the H3N2 subtype. Our results support recent indications that antibodies binding to the "stalk" region of hemagglutinin are found in the human population and exert evolutionary pressure on the virus. Our computational approach provides a possible method for identifying antigenic escape through evolution in this region, which in some cases will not be identified by the hemagglutinin inhibition assay.     

The impact of seasonal and year-round transmission regimes on the evolution of influenza A virus

Punctuated antigenic change is believed to be a key element in the evolution of influenza A; clusters of antigenically similar strains predominate worldwide for several years until an antigenically distant mutant emerges and instigates a selective sweep. It is thought that a region of East-Southeast Asia with year-round transmission acts as a source of antigenic diversity for influenza A and seasonal epidemics in temperate regions make little contribution to antigenic evolution. We use a mathematical model to examine how different transmission regimes affect the evolutionary dynamics of influenza over the lifespan of an antigenic cluster. Our model indicates that, in non-seasonal regions, mutants that cause significant outbreaks appear before the peak of the wild-type epidemic. A relatively large proportion of these mutants spread globally. In seasonal regions, mutants that cause significant local outbreaks appear each year before the seasonal peak of the wild-type epidemic, but only a small proportion spread globally. The potential for global spread is strongly influenced by the intensity of non-seasonal circulation and coupling between non-seasonal and seasonal regions. Results are similar if mutations are neutral, or confer a weak to moderate antigenic advantage. However, there is a threshold antigenic advantage, depending on the non-seasonal transmission intensity, beyond which mutants can escape herd immunity in the non-seasonal region and there is a global explosion in diversity. We conclude that non-seasonal transmission regions are fundamental to the generation and maintenance of influenza diversity owing to their epidemiology. More extensive sampling of viral diversity in such regions could facilitate earlier identification of antigenically novel strains and extend the critical window for vaccine development.     

Limited Predictability of Amino Acid Substitutions in Seasonal Influenza Viruses

Seasonal influenza viruses repeatedly infect humans in part because they rapidly change their antigenic properties and evade host immune responses, necessitating frequent updates of the vaccine composition. Accurate predictions of strains circulating in the future could therefore improve the vaccine match. Here, we studied the predictability of frequency dynamics and fixation of amino acid substitutions. Current frequency was the strongest predictor of eventual fixation, as expected in neutral evolution. Other properties, such as occurrence in previously characterized epitopes or high Local Branching Index (LBI) had little predictive power. Parallel evolution was found to be moderately predictive of fixation. Although the LBI had little power to predict frequency dynamics, it was still successful at picking strains representative of future populations. The latter is due to a tendency of the LBI to be high for consensus-like sequences that are closer to the future than the average sequence. Simulations of models of adapting populations, in contrast, show clear signals of predictability. This indicates that the evolution of influenza HA and NA, while driven by strong selection pressure to change, is poorly described by common models of directional selection such as traveling fitness waves.     

Multiplexed screening of natural humoral immunity identifies antibodies at fine specificity for complex and dynamic viral targets

Viral entry targets with therapeutic neutralizing potential are subject to multiple escape mechanisms, including antigenic drift, immune dominance of functionally irrelevant epitopes, and subtle variations in host cell mechanisms. A surprising finding of recent years is that potent neutralizing antibodies to viral epitopes independent of strain exist, but are poorly represented across the diverse human population. Identifying these antibodies and understanding the biology mediating the specific immune response is thus difficult. An effective strategy for meeting this challenge is to incorporate multiplexed antigen screening into a high throughput survey of the memory B cell repertoire from immune individuals. We used this approach to discover suites of cross-clade antibodies directed to conformational epitopes in the stalk region of the influenza A hemagglutinin (HA) protein and to select high-affinity anti-peptide antibodies to the glycoprotein B (gB) of human cytomegalovirus. In each case, our screens revealed a restricted V_H and V_L germline usage, including published and previously unidentified gene families. The in vivo evolution of paratope specificity with optimal neutralizing activity was understandable after correlating biological activities with kinetic binding and epitope recognition. Iterative feedback between antigen probe design based on structure and function information with high throughput multiplexed screening demonstrated a generally applicable strategy for efficient identification of safe, native, finely tuned antibodies with the potential for high genetic barriers to viral escape.     

Multiplexed screening of natural humoral immunity identifies antibodies at fine specificity for complex and dynamic viral targets

Viral entry targets with therapeutic neutralizing potential are subject to multiple escape mechanisms, including antigenic drift, immune dominance of functionally irrelevant epitopes, and subtle variations in host cell mechanisms. A surprising finding of recent years is that potent neutralizing antibodies to viral epitopes independent of strain exist, but are poorly represented across the diverse human population. Identifying these antibodies and understanding the biology mediating the specific immune response is thus difficult. An effective strategy for meeting this challenge is to incorporate multiplexed antigen screening into a high throughput survey of the memory B cell repertoire from immune individuals. We used this approach to discover suites of cross-clade antibodies directed to conformational epitopes in the stalk region of the influenza A hemagglutinin (HA) protein and to select high-affinity anti-peptide antibodies to the glycoprotein B (gB) of human cytomegalovirus. In each case, our screens revealed a restricted V_H and V_L germline usage, including published and previously unidentified gene families. The in vivo evolution of paratope specificity with optimal neutralizing activity was understandable after correlating biological activities with kinetic binding and epitope recognition. Iterative feedback between antigen probe design based on structure and function information with high throughput multiplexed screening demonstrated a generally applicable strategy for efficient identification of safe, native, finely tuned antibodies with the potential for high genetic barriers to viral escape.     

Assessing Antigenic Drift of Seasonal Influenza A(H3N2) and A(H1N1)pdm09 Viruses

Under selective pressure from the host immune system, antigenic epitopes of influenza virus hemagglutinin (HA) have continually evolved to escape antibody recognition, termed antigenic drift. We analyzed the genomes of influenza A(H3N2) and A(H1N1)pdm09 virus strains circulating in Thailand between 2010 and 2014 and assessed how well the yearly vaccine strains recommended for the southern hemisphere matched them. We amplified and sequenced the HA gene of 120 A(H3N2) and 81 A(H1N1)pdm09 influenza virus samples obtained from respiratory specimens and calculated the perfect-match vaccine efficacy using the p _epitope model, which quantitated the antigenic drift in the dominant epitope of HA. Phylogenetic analysis of the A(H3N2) HA1 genes classified most strains into genetic clades 1, 3A, 3B, and 3C. The A(H3N2) strains from the 2013 and 2014 seasons showed very low to moderate vaccine efficacy and demonstrated antigenic drift from epitopes C and A to epitope B. Meanwhile, most A(H1N1)pdm09 strains from the 2012–2014 seasons belonged to genetic clades 6A, 6B, and 6C and displayed the dominant epitope mutations at epitopes B and E. Finally, the vaccine efficacy for A(H1N1)pdm09 (79.6–93.4%) was generally higher than that of A(H3N2). These findings further confirmed the accelerating antigenic drift of the circulating influenza A(H3N2) in recent years.     

Positive selection operates continuously on hemagglutinin during evolution of H3N2 human influenza A virus

It has been proposed that antigenic evolution of hemagglutinin 1 (HA1) for H3N2 human influenza A virus was punctuated. In the population genetic analysis, however, it was controversial whether positive selection operated on HA1 in a punctuated manner for the branches of the phylogenetic tree where transitions to new antigenic clusters occurred (C branches), or continuously. In the molecular evolutionary analysis, positive selection was detected for the trunk (T) branches but the relationship between antigenic evolution and positive selection was unclear. Here molecular evolutionary analysis was conducted to examine natural selection operating on HA1 of H3N2 human influenza A virus by dividing HA1 into epitopes A-E and other sites, as well as dividing the phylogenetic tree into the C branches overlapping with the T branches (C-T branches), those not overlapping with the T branches (C-NT branches), the T branches not overlapping with the C branches (NC-T branches), and other branches (NC-NT branches). Positive selection was detected for C, T, and NC-T branches, whereas evolution for the NC-NT branches appeared to be mainly neutral. Positive selection appeared to have operated throughout the trunk, which covered the entire time period of the phylogenetic tree, suggesting that positive selection operated continuously on HA1 during evolution of H3N2 human influenza A virus.     

Structural differences among hemagglutinins of influenza A virus subtypes are reflected in their antigenic architecture: analysis of H9 escape mutants

We used a panel of monoclonal antibodies to H9 hemagglutinin to select 18 escape mutants of mouse-adapted influenza A/Swine/Hong Kong/9/98 (H9N2) virus. Cross-reactions of the mutants with the antibodies and the sequencing of hemagglutinin genes revealed two minimally overlapping epitopes. We mapped the amino acid changes to two areas of the recently reported three-dimensional structure of A/Swine/Hong Kong/9/98 hemagglutinin. The grouping of the antigenically relevant amino acid positions in H9 hemagglutinin differs from the pattern observed in H3 and H5 hemagglutinins. Several positions in site B of H3 hemagglutinin are distributed in two sites of H9 hemagglutinin. Unlike any subtype analyzed so far, H9 hemagglutinin does not contain an antigenic site corresponding to site A in H3 hemagglutinin. Positions 145 and 193 (H3 numbering), which in H3 hemagglutinin belong to sites A and B, respectively, are within one site in H9 hemagglutinin. This finding is consistent with the peculiarity of the three-dimensional structure of the H9 molecule, that is, the absence from H9 hemagglutinin of the lateral loop that forms site A in H3 and the equivalent site in H5 hemagglutinins. The escape mutants analyzed displayed phenotypic variations, including decreased virulence for mice and changes in affinity for sialyl substrates. Our results demonstrate a correlation between intersubtype differences in three-dimensional structure and variations among subtypes in the distribution of antigenic areas. Our findings also suggest that covariation and pleiotropic effects of antibody-selected mutations may be important in the evolution of H9 influenza virus, a possible causative agent of a future pandemic.     

The First Steps of Adaptation of Escherichia coli to the Gut Are Dominated by Soft Sweeps

The accumulation of adaptive mutations is essential for survival in novel environments. However, in clonal populations with a high mutational supply, the power of natural selection is expected to be limited. This is due to clonal interference - the competition of clones carrying different beneficial mutations - which leads to the loss of many small effect mutations and fixation of large effect ones. If interference is abundant, then mechanisms for horizontal transfer of genes, which allow the immediate combination of beneficial alleles in a single background, are expected to evolve. However, the relevance of interference in natural complex environments, such as the gut, is poorly known. To address this issue, we have developed an experimental system which allows to uncover the nature of the adaptive process as Escherichia coli adapts to the mouse gut. This system shows the invasion of beneficial mutations in the bacterial populations and demonstrates the pervasiveness of clonal interference. The observed dynamics of change in frequency of beneficial mutations are consistent with soft sweeps, where different adaptive mutations with similar phenotypes, arise repeatedly on different haplotypes without reaching fixation. Despite the complexity of this ecosystem, the genetic basis of the adaptive mutations revealed a striking parallelism in independently evolving populations. This was mainly characterized by the insertion of transposable elements in both coding and regulatory regions of a few genes. Interestingly, in most populations we observed a complete phenotypic sweep without loss of genetic variation. The intense clonal interference during adaptation to the gut environment, here demonstrated, may be important for our understanding of the levels of strain diversity of E. coli inhabiting the human gut microbiota and of its recombination rate.     

Escape from a dominant HLA-B*15-restricted CD8+ T cell response against hepatitis C virus requires compensatory mutations outside the epitope

Antiviral CD8(+) T cells are a key component of the adaptive immune system against hepatitis C virus (HCV). For the development of immune therapies, it is essential to understand how CD8(+) T cells contribute to clearance of infection and why they fail so often. A mechanism for secondary failure is mutational escape of the virus. However, some substitutions in viral epitopes are associated with fitness costs and often require compensatory mutations. We hypothesized that compensatory mutations may point toward epitopes under particularly strong selection pressure that may be beneficial for vaccine design because of a higher genetic barrier to escape. We previously identified two HLA-B*15-restricted CD8(+) epitopes in NS5B (LLRHHNMVY(2450-2458) and SQRQKKVTF(2466-2474)), based on sequence analysis of a large HCV genotype 1b outbreak. Both epitopes are targeted in about 70% of HLA-B*15-positive individuals exposed to HCV. Reproducible selection of escape mutations was confirmed in an independent multicenter cohort in the present study. Interestingly, mutations were also selected in the epitope flanking region, suggesting that compensatory evolution may play a role. Covariation analysis of sequences from the database confirmed a significant association between escape mutations inside one of the epitopes (H2454R and M2456L) and substitutions in the epitope flanking region (S2439T and K2440Q). Functional analysis with the subgenomic replicon Con1 confirmed that the primary escape mutations impaired viral replication, while fitness was restored by the additional substitutions in the epitope flanking region. We concluded that selection of escape mutations inside an HLA-B*15 epitope requires secondary substitutions in the epitope flanking region that compensate for fitness costs.