Evolution of the Influenza A Virus Genome during Development of Oseltamivir Resistance In Vitro

Influenza A virus (IAV) is a major cause of morbidity and mortality throughout the world. Current antiviral therapies include oseltamivir, a neuraminidase inhibitor that prevents the release of nascent viral particles from infected cells. However, the IAV genome can evolve rapidly, and oseltamivir resistance mutations have been detected in numerous clinical samples. Using an in vitro evolution platform and whole-genome population sequencing, we investigated the population genomics of IAV during the development of oseltamivir resistance. Strain A/Brisbane/59/2007 (H1N1) was grown in Madin-Darby canine kidney cells with or without escalating concentrations of oseltamivir over serial passages. Following drug treatment, the H274Y resistance mutation fixed reproducibly within the population. The presence of the H274Y mutation in the viral population, at either a low or a high frequency, led to measurable changes in the neuraminidase inhibition assay. Surprisingly, fixation of the resistance mutation was not accompanied by alterations of viral population diversity or differentiation, and oseltamivir did not alter the selective environment. While the neighboring K248E mutation was also a target of positive selection prior to H274Y fixation, H274Y was the primary beneficial mutation in the population. In addition, once evolved, the H274Y mutation persisted after the withdrawal of the drug, even when not fixed in viral populations. We conclude that only selection of H274Y is required for oseltamivir resistance and that H274Y is not deleterious in the absence of the drug. These collective results could offer an explanation for the recent reproducible rise in oseltamivir resistance in seasonal H1N1 IAV strains in humans.     

Coadministration of Hedera helix L. Extract Enabled Mice to Overcome Insufficient Protection against Influenza A/PR/8 Virus Infection under Suboptimal Treatment with Oseltamivir

Several anti-influenza drugs that reduce disease manifestation exist, and although these drugs provide clinical benefits in infected patients, their efficacy is limited by the emergence of drug-resistant influenza viruses. In the current study, we assessed the therapeutic strategy of enhancing the antiviral efficacy of an existing neuraminidase inhibitor, oseltamivir, by coadministering with the leaf extract from Hedera helix L, commonly known as ivy. Ivy extract has anti-inflammatory, antibacterial, antifungal, and antihelminthic properties. In the present study, we investigated its potential antiviral properties against influenza A/PR/8 (PR8) virus in a mouse model with suboptimal oseltamivir that mimics a poor clinical response to antiviral drug treatment. Suboptimal oseltamivir resulted in insufficient protection against PR8 infection. Oral administration of ivy extract with suboptimal oseltamivir increased the antiviral activity of oseltamivir. Ivy extract and its compounds, particularly hedrasaponin F, significantly reduced the cytopathic effect in PR8-infected A549 cells in the presence of oseltamivir. Compared with oseltamivir treatment alone, coadministration of the fraction of ivy extract that contained the highest proportion of hedrasaponin F with oseltamivir decreased pulmonary inflammation in PR8-infected mice. Inflammatory cytokines and chemokines, including tumor necrosis factor-alpha and chemokine (C-C motif) ligand 2, were reduced by treatment with oseltamivir and the fraction of ivy extract. Analysis of inflammatory cell infiltration in the bronchial alveolar of PR8-infected mice revealed that CD11b⁺Ly6G⁺ and CD11b⁺Ly6C^int cells were recruited after virus infection; coadministration of the ivy extract fraction with oseltamivir reduced infiltration of these inflammatory cells. In a model of suboptimal oseltamivir treatment, coadministration of ivy extract fraction that includes hedrasaponin F increased protection against PR8 infection that could be explained by its antiviral and anti-inflammatory activities.     

Oseltamivir Resistance in Influenza A(H6N2) Caused by an R292K Substitution in Neuraminidase Is Not Maintained in Mallards without Drug Pressure

Background

Wild waterfowl is the natural reservoir of influenza A virus (IAV); hosted viruses are very variable and provide a source for genetic segments which can reassort with poultry or mammalian adapted IAVs to generate novel species crossing viruses. Additionally, wild waterfowl act as a reservoir for highly pathogenic IAVs. Exposure of wild birds to the antiviral drug oseltamivir may occur in the environment as its active metabolite can be released from sewage treatment plants to river water. Resistance to oseltamivir, or to other neuraminidase inhibitors (NAIs), in IAVs of wild waterfowl has not been extensively studied.

Aim and Methods

In a previous in vivo Mallard experiment, an influenza A(H6N2) virus developed oseltamivir resistance by the R292K substitution in the neuraminidase (NA), when the birds were exposed to oseltamivir. In this study we tested if the resistance could be maintained in Mallards without drug exposure. Three variants of resistant H6N2/R292K virus were each propagated during 17 days in five successive pairs of naïve Mallards, while oseltamivir exposure was decreased and removed. Daily fecal samples were analyzed for viral presence, genotype and phenotype.

Results and Conclusion

Within three days without drug exposure no resistant viruses could be detected by NA sequencing, which was confirmed by functional NAI sensitivity testing. We conclude that this resistant N2 virus could not compete in fitness with wild type subpopulations without oseltamivir drug pressure, and thus has no potential to circulate among wild birds. The results of this study contrast to previous observations of drug induced resistance in an avian H1N1 virus, which was maintained also without drug exposure in Mallards. Experimental observations on persistence of NAI resistance in avian IAVs resemble NAI resistance seen in human IAVs, in which resistant N2 subtypes do not circulate, while N1 subtypes with permissive mutations can circulate without drug pressure. We speculate that the phylogenetic group N1 NAs may easier compensate for NAI resistance than group N2 NAs, though further studies are needed to confirm such conclusions.     

Triple combination antiviral drug (TCAD) composed of amantadine, oseltamivir, and ribavirin impedes the selection of drug-resistant influenza A virus

Widespread resistance among circulating influenza A strains to at least one of the anti-influenza drugs is a major public health concern. A triple combination antiviral drug (TCAD) regimen comprised of amantadine, oseltamivir, and ribavirin has been shown to have synergistic and broad spectrum activity against influenza A strains, including drug resistant strains. Here, we used mathematical modeling along with three different experimental approaches to understand the effects of single agents, double combinations, and the TCAD regimen on resistance in influenza in vitro, including: 1) serial passage at constant drug concentrations, 2) serial passage at escalating drug concentrations, and 3) evaluation of the contribution of each component of the TCAD regimen to the suppression of resistance. Consistent with the modeling which demonstrated that three drugs were required to suppress the emergence of resistance in influenza A, treatment with the TCAD regimen resulted in the sustained suppression of drug resistant viruses, whereas treatment with amantadine alone or the amantadine-oseltamivir double combination led to the rapid selection of resistant variants which comprised ～100% of the population. Furthermore, the TCAD regimen imposed a high genetic barrier to resistance, requiring multiple mutations in order to escape the effects of all the drugs in the regimen. Finally, we demonstrate that each drug in the TCAD regimen made a significant contribution to the suppression of virus breakthrough and resistance at clinically achievable concentrations. Taken together, these data demonstrate that the TCAD regimen was superior to double combinations and single agents at suppressing resistance, and that three drugs at a minimum were required to impede the selection of drug resistant variants in influenza A virus. The use of mathematical modeling with multiple experimental designs and molecular readouts to evaluate and optimize combination drug regimens for the suppression of resistance may be broadly applicable to other infectious diseases.     

Evolutionary novelties and losses in geometric morphometrics: a practical approach through hominin molar morphology

Geometric morphometric techniques may offer a promising methodological approach to analyze evolutionary novelties in a quantitative framework. Nevertheless, and despite continuous improvements to this methodology, the inclusion of novel features in these studies presents some difficulties. In the present study, different methods to explicitly include novel traits in geometric morphometric analyses are compared, including homology-free approaches, landmark-based approaches, and combinations of both techniques. The two-dimensional occlusal morphology of the lower second molar in multiple hominin species was chosen to evaluate these methods, as an example of an anatomical structure including one novelty: a distal fifth cusp is present in earlier hominins, and notably absent in many later Homo species. Results reveal that different approaches provide different results, highlighting that the design of the conformations of landmarks has a high impact on the inferred conclusions. Among diverse methods, a combined approach including landmarks, sliding semilandmarks, and only one landmark related to the studied novelty (an indicator of its absence or presence and of its size, when present), was able to directly discern structures with and without the novel feature, circumventing some of the methodological difficulties associated with these traits. This study demonstrates the ability of geometric morphometric techniques to investigate evolutionary novelties and explores the implications of different methods, providing a reference context for future studies.     

Simultaneous analysis of multiple data types in pharmacogenomic studies using weighted sparse canonical correlation analysis

Variation in drug response results from a combination of factors that include differences in gender, ethnicity, and environment, as well as genetic variation that may result in differences in mRNA and protein expression. This article presents two integrative analytic approaches that make use of both genome-wide SNP and mRNA expression data available on the same set of subjects: a step-wise integrative approach and a comprehensive analysis using sparse canonical correlation analysis (SCCA). In addition to applying standard SCCA, we present a novel modification of SCCA which allows different weighting for the various pair-wise relationships in the SCCA. These integrative approaches are illustrated with both simulated data and data from a pharmacogenomic study of the drug gemcitabine. Results from these analyses found little overlap in terms of genes detected, possibly detecting different biological mechanisms. In addition, we found the proposed weighted SCCA to outperform its unweighted counterpart in detecting associations between the genomic features and phenotype. Further research is needed to develop and assess new integrative methods for pharmacogenomic studies, as these types of analyses may uncover novel insights into the relationship between genomic variation and drug response.     

Cutting Edge: Oseltamivir decreases T cell GM1 expression and inhibits clearance of respiratory syncytial virus: potential role of endogenous sialidase in antiviral immunity

The sialoglycosphingolipid GM1 is important for lipid rafts and immune cell signaling. T cell activation in vitro increases GM1 expression and increases endogenous sialidase activity. GM1 expression has been hypothesized to be regulated by endogenous sialidase. We tested this hypothesis in vivo using a mouse model of respiratory syncytial virus (RSV) infection. RSV infection increased endogenous sialidase activity in lung mononuclear cells. RSV infection increased lung CD8+ T cell surface GM1 expression. Activated CD8+ T cells in the lungs of RSV-infected mice were GM1(high). Treatment of RSV-infected mice with the sialidase/neuraminidase inhibitor oseltamivir decreased T cell surface GM1 levels. Oseltamivir treatment decreased RSV-induced weight loss and inhibited RSV clearance. Our data indicate a novel role for an endogenous sialidase in regulating T cell GM1 expression and antiviral immunity. Also, oseltamivir, an important anti-influenza drug, inhibits the clearance of a respiratory virus that lacks a neuraminidase gene, RSV.     

Impact of Different Oseltamivir Regimens on Treating Influenza A Virus Infection and Resistance Emergence: Insights from a Modelling Study

Several studies have proven oseltamivir to be efficient in reducing influenza viral titer and symptom intensity. However, the usefulness of oseltamivir can be compromised by the emergence and spread of drug-resistant virus. The selective pressure exerted by different oseltamivir therapy regimens have received little attention. Combining models of drug pharmacokinetics, pharmacodynamics, viral kinetics and symptom dynamics, we explored the efficacy of oseltamivir in reducing both symptoms (symptom efficacy) and viral load (virological efficacy). We simulated samples of 1000 subjects using previously estimated between-subject variability in viral and symptom dynamic parameters to describe the observed heterogeneity in a patient population. We simulated random mutations conferring resistance to oseltamivir. We explored the effect of therapy initiation time, dose, intake frequency and therapy duration on influenza infection, illness dynamics, and emergence of viral resistance. Symptom and virological efficacies were strongly associated with therapy initiation time. The proportion of subjects shedding resistant virus was 27-fold higher when prophylaxis was initiated during the incubation period compared with no treatment. It fell to below 1% when treatment was initiated after symptom onset for twice-a-day intakes. Lower doses and prophylaxis regimens led to lower efficacies and increased risk of resistance emergence. We conclude that prophylaxis initiated during the incubation period is the main factor leading to resistance emergence.     

Identification of Potential Inhibitors of H5N1 Influenza A Virus Neuraminidase by Ligand-Based Virtual Screening Approach

The neuraminidase (NA) of the influenza virus is the target of antiviral drug, oseltamivir. Recently, cases were reported that influenza virus becoming resistant to oseltamivir, necessitating the development of new long-acting antiviral compounds. In this report, a novel class of lead molecule with potential NA inhibitory activity was identified using a combination of virtual screening (VS), molecular docking, and molecular dynamic approach. The PubChem database was used to perform the VS analysis by employing oseltamivir as query. Subsequently, the data reduction was carried out by employing molecular docking study. Furthermore, the screened lead molecules were analyzed with respect to the Lipinski rule of five, drug-likeness, toxicity profiles, and other physico-chemical properties of drugs by suitable software program. Final screening was carried out by normal mode analysis and molecular dynamic simulation approach. The result indicates that CID 25145634, deuterium-enriched oseltamivir, become a promising lead compound and be effective in treating oseltamivir sensitive as well as resistant influenza virus strains.     

Computational studies of H5N1 influenza virus resistance to oseltamivir

Influenza A (H5N1) virus is one of the world's greatest pandemic threats. Neuraminidase (NA) inhibitors, oseltamivir and zanamivir, prevent the spread of influenza, but drug‐resistant viruses have reduced their effectiveness. Resistance depends on the binding properties of NA‐drug complexes. Key residue mutations within the active site of NA glycoproteins diminish binding, thereby resulting in drug resistance. We performed molecular simulations and calculations to characterize the mechanisms of H5N1 influenza virus resistance to oseltamivir and predict potential drug‐resistant mutations. We examined two resistant NA mutations, H274Y and N294S, and one non‐drug‐resistant mutation, E119G. Six‐nanosecond unrestrained molecular dynamic simulations with explicit solvent were performed using NA‐oseltamivir complexes containing either NA wild‐type H5N1 virus or a variant. MM_PBSA techniques were then used to rank the binding free energies of these complexes. Detailed analyses indicated that conformational change of E276 in the Pocket 1 region of NA is a key source of drug resistance in the H274Y mutant but not in the N294S mutant.     

Oseltamivir-resistant influenza A viruses are transmitted efficiently among guinea pigs by direct contact but not by aerosol

Influenza viruses resistant to the neuraminidase (NA) inhibitor oseltamivir arise under drug selection pressure both in vitro and in vivo. Several mutations in the active site of the viral NA are known to confer relative resistance to oseltamivir, and influenza viruses with certain oseltamivir resistance mutations have been shown to transmit efficiently among cocaged ferrets. However, it is not known whether NA mutations alter aerosol transmission of drug-resistant influenza virus. Here, we demonstrate that recombinant human influenza A/H3N2 viruses without and with oseltamivir resistance mutations (in which NA carries the mutation E119V or the double mutations E119V I222V) have similar in ovo growth kinetics and infectivity in guinea pigs. These viruses also transmit efficiently by the contact route among cocaged guinea pigs, as in the ferret model. However, in an aerosol transmission model, in which guinea pigs are caged separately, the oseltamivir-resistant viruses transmit poorly or not at all; in contrast, the oseltamivir-sensitive virus transmits efficiently even in the absence of direct contact. The present results suggest that oseltamivir resistance mutations reduce aerosol transmission of influenza virus, which could have implications for public health measures taken in the event of an influenza pandemic.     

Triple Combination of Amantadine,Ribavirin, and Oseltamivir Is Highly Active and Synergistic against Drug Resistant Influenza Virus Strains In Vitro

The rapid emergence and subsequent spread of the novel 2009 Influenza A/H1N1 virus (2009 H1N1) has prompted the World Health Organization to declare the first pandemic of the 21^st century, highlighting the threat of influenza to public health and healthcare systems. Widespread resistance to both classes of influenza antivirals (adamantanes and neuraminidase inhibitors) occurs in both pandemic and seasonal viruses, rendering these drugs to be of marginal utility in the treatment modality. Worldwide, virtually all 2009 H1N1 and seasonal H3N2 strains are resistant to the adamantanes (rimantadine and amantadine), and the majority of seasonal H1N1 strains are resistant to oseltamivir, the most widely prescribed neuraminidase inhibitor (NAI). To address the need for more effective therapy, we evaluated the in vitro activity of a triple combination antiviral drug (TCAD) regimen composed of drugs with different mechanisms of action against drug-resistant seasonal and 2009 H1N1 influenza viruses. Amantadine, ribavirin, and oseltamivir, alone and in combination, were tested against amantadine- and oseltamivir-resistant influenza A viruses using an in vitro infection model in MDCK cells. Our data show that the triple combination was highly synergistic against drug-resistant viruses, and the synergy of the triple combination was significantly greater than the synergy of any double combination tested (P<0.05), including the combination of two NAIs. Surprisingly, amantadine and oseltamivir contributed to the antiviral activity of the TCAD regimen against amantadine- and oseltamivir-resistant viruses, respectively, at concentrations where they had no activity as single agents, and at concentrations that were clinically achievable. Our data demonstrate that the TCAD regimen composed of amantadine, ribavirin, and oseltamivir is highly synergistic against resistant viruses, including 2009 H1N1. The TCAD regimen overcomes baseline drug resistance to both classes of approved influenza antivirals, and thus may represent a highly active antiviral therapy for seasonal and pandemic influenza.     

Nutritional Systems Biology Modeling: From Molecular Mechanisms to Physiology

The use of computational modeling and simulation has increased in many biological fields, but despite their potential these techniques are only marginally applied in nutritional sciences. Nevertheless, recent applications of modeling have been instrumental in answering important nutritional questions from the cellular up to the physiological levels. Capturing the complexity of today''s important nutritional research questions poses a challenge for modeling to become truly integrative in the consideration and interpretation of experimental data at widely differing scales of space and time. In this review, we discuss a selection of available modeling approaches and applications relevant for nutrition. We then put these models into perspective by categorizing them according to their space and time domain. Through this categorization process, we identified a dearth of models that consider processes occurring between the microscopic and macroscopic scale. We propose a “middle-out” strategy to develop the required full-scale, multilevel computational models. Exhaustive and accurate phenotyping, the use of the virtual patient concept, and the development of biomarkers from “-omics” signatures are identified as key elements of a successful systems biology modeling approach in nutrition research—one that integrates physiological mechanisms and data at multiple space and time scales.     

Radiosyntheses of two positron emission tomography probes: [11C]Oseltamivir and its active metabolite [11C]Ro 64-0802

Oseltamivir phosphate (Tamiflu, 1.H(3)PO(4) is an orally active anti-influenza drug, which is hydrolyzed by esterase to its carboxylate metabolite Ro 64-0802 (2) with potent activity inhibiting neuraminidase. In this study, for the first time, we synthesized carbon-11-labeled oseltamivir ([(11)C]1) and Ro 64-0802 ([(11)C]2) as two novel positron emission tomography probes and demonstrated that [(11)C]1 had twofold higher radioactivity concentration in the mouse brains than [(11)C]2.     

Evolutionary potential of morphological traits across different life‐history stages

Despite accumulating examples of selection acting on heritable traits in the wild, predicted evolutionary responses are often different from observed phenotypic trends. Various explanations have been suggested for these mismatches. These include within‐individual changes across lifespan that can create important variation in genetic architecture of traits and selection acting on them, but also potential problems with the methodological approach used to predict evolutionary responses of traits. Here, we used an 8‐year data set on tree swallow (Tachycineta bicolor) to first assess the effects of differences among three nestling life‐history stages on the genetic (co)variances of two morphological traits (body mass and primary feather length) and the selection acting on them over three generations. We then estimated the evolutionary potential of these traits by predicting their evolutionary responses using the breeder's equation and the secondary theorem of selection approaches. Our results showed variation in strength and direction of selection and slight changes in trait variance across ages. Predicted evolutionary responses differed importantly between both approaches for half of the trait–age combinations we studied, suggesting the presence of environmentally induced correlations between focal traits and fitness possibly biasing breeder's equation predictions. Our results emphasize that predictions of evolutionary potential for morphological traits are likely to be highly variable, both in strength and direction, depending on the life stage and method used, thus mitigating our capacity to predict adaptation and persistence of wild populations.     

In silico prediction of drug resistance due to S247R mutation of Influenza H1N1 neuraminidase protein

We present here in silico studies on antiviral drug resistance due to a novel mutation of influenza A/H1N1 neuraminidase (NA) protein. Influenza A/H1N1 virus was responsible for a recent pandemic and is currently circulating among the seasonal influenza strains. M2 and NA are the two major viral proteins related to pathogenesis in humans and have been targeted for drug designing. Among them, NA is preferred because the ligand-binding site of NA is highly conserved between different strains of influenza virus. Different mutations of the NA active site residues leading to drug resistance or susceptibility of the virus were studied earlier. We report here a novel mutation (S247R) in the NA protein that was sequenced earlier from the nasopharyngeal swab from Sri Lanka and Thailand in the year 2009 and 2011, respectively. Another mutation (S247N) was already known to confer resistance to oseltamivir. We did a comparative study of these two mutations vis-a-vis the drug-sensitive wild type NA to understand the mechanism of drug resistance of S247N and to predict the probability of the novel S247R mutation to become resistant to the currently available drugs, oseltamivir and zanamivir. We performed molecular docking- and molecular dynamics-based analysis of both the mutant proteins and showed that mutation of S247R affects drug binding to the protein by positional displacement due to altered active site cavity architecture, which in turn reduces the affinity of the drug molecules to the NA active site. Our analysis shows that S247R may have high probability of being resistant.     

Prevalence of epistasis in the evolution of influenza A surface proteins

The surface proteins of human influenza A viruses experience positive selection to escape both human immunity and, more recently, antiviral drug treatments. In bacteria and viruses, immune-escape and drug-resistant phenotypes often appear through a combination of several mutations that have epistatic effects on pathogen fitness. However, the extent and structure of epistasis in influenza viral proteins have not been systematically investigated. Here, we develop a novel statistical method to detect positive epistasis between pairs of sites in a protein, based on the observed temporal patterns of sequence evolution. The method rests on the simple idea that a substitution at one site should rapidly follow a substitution at another site if the sites are positively epistatic. We apply this method to the surface proteins hemagglutinin and neuraminidase of influenza A virus subtypes H3N2 and H1N1. Compared to a non-epistatic null distribution, we detect substantial amounts of epistasis and determine the identities of putatively epistatic pairs of sites. In particular, using sequence data alone, our method identifies epistatic interactions between specific sites in neuraminidase that have recently been demonstrated, in vitro, to confer resistance to the drug oseltamivir; these epistatic interactions are responsible for widespread drug resistance among H1N1 viruses circulating today. This experimental validation demonstrates the predictive power of our method to identify epistatic sites of importance for viral adaptation and public health. We conclude that epistasis plays a large role in shaping the molecular evolution of influenza viruses. In particular, sites with , which would normally not be identified as positively selected, can facilitate viral adaptation through epistatic interactions with their partner sites. The knowledge of specific interactions among sites in influenza proteins may help us to predict the course of antigenic evolution and, consequently, to select more appropriate vaccines and drugs.     

Characterization of two distinct neuraminidases from avian-origin human-infecting H7N9 influenza viruses

An epidemic of an avian-origin H7N9 influenza virus has recently emerged in China, infecting 134 patients of which 45 have died. This is the first time that an influenza virus harboring an N9 serotype neuraminidase (NA) has been known to infect humans. H7N9 viruses are divergent and at least two distinct NAs and hemagglutinins (HAs) have been found, respectively, from clinical isolates. The prototypes of these viruses are A/Anhui/1/2013 and A/Shanghai/1/2013. NAs from these two viruses are distinct as the A/Shanghai/1/2013 NA has an R294K substitution that can confer NA inhibitor oseltamivir resistance. Oseltamivir is by far the most commonly used anti-influenza drug due to its potency and high bioavailability. In this study, we show that an R294K substitution results in multidrug resistance with extreme oseltamivir resistance (over 100 000-fold) using protein- and virus-based assays. To determine the molecular basis for the inhibitor resistance, we solved high-resolution crystal structures of NAs from A/Anhui/1/2013 N9 (R294-containing) and A/Shanghai/1/2013 N9 (K294-containing). R294K substitution results in an unfavorable E276 conformation for oseltamivir binding, and consequently loss of inhibitor carboxylate interactions, which compromises the binding of all classical NA ligands/inhibitors. Moreover, we found that R294K substitution results in reduced NA catalytic efficiency along with lower viral fitness. This helps to explain why K294 has predominantly been found in clinical cases of H7N9 infection under the selective pressure of oseltamivir treatment and not in the dominant human-infecting viruses. This implies that oseltamivir can still be efficiently used in the treatment of H7N9 infections.     

Artificial selection on larval growth curves in Tribolium: correlated responses and constraints

Body size is often constrained from evolving. Although artificial selection on body size in insects frequently results in a sizable response, these responses usually bear fitness costs. Further, these experiments tend to select only on size at one landmark age, rather than selecting for patterns of growth over the whole larval life stage. To address whether constraints may be caused by larval growth patterns rather than final size, we implemented a function‐valued (FV) trait method of selection, in which entire larval growth curves from Tribolium were artificially selected. The selection gradient function used was previously predicted to give the maximal response and was implemented using a novel selection index in the FV framework. Results indicated a significant response after one generation of selection, but no response in subsequent generations. Correlated responses included increased mortality, increased critical weight, and decreased development time (DT). The lack of response in size and development time after the first generation was likely caused by increased mortality suffered in selected lines; we demonstrated that the selection criterion caused both increased body size and increased mortality. We conclude that artificial selection on continuous traits using FV methods is very efficient and that the constraint of body size evolution is likely caused by a suite of trade‐offs with other traits.