Does Plant Breeding for Antioxidant-Rich Foods Have an Impact on Human Health?

Given the general beneficial effects of antioxidants-rich foods on human health and disease prevention, there is a continuous interest in plant secondary metabolites conferring attractive colors to fruits and grains and responsible, together with others, for nutraceutical properties. Cereals and Solanaceae are important components of the human diet, thus, they are the main targets for functional food development by exploitation of genetic resources and metabolic engineering. In this review, we focus on the impact of antioxidants-rich cereal and Solanaceae derived foods on human health by analyzing natural biodiversity and biotechnological strategies aiming at increasing the antioxidant level of grains and fruits, the impact of agronomic practices and food processing on antioxidant properties combined with a focus on the current state of pre-clinical and clinical studies. Despite the strong evidence in in vitro and animal studies supporting the beneficial effects of antioxidants-rich diets in preventing diseases, clinical studies are still not sufficient to prove the impact of antioxidant rich cereal and Solanaceae derived foods on human.

Unveiling gene expression regulation of the Bacillus thuringiensis Cry3Aa toxin receptor ADAM10 by the potato dietary miR171c in Colorado potato beetle.

BACKGROUND: The Colorado potato beetle (CPB) is a worldwide devastating pest of potato plants and other Solanaceae characterized by its remarkable ability to evolve resistance to insecticides. Bacillus thuringiensis (Bt) Cry3Aa toxin represents an environmentally safe alternative for CPB control but larvae susceptibility to this toxin has been reported to vary depending on the host plant on which larvae feed. To gain more insight into how nutrition mediates Bt tolerance through effects on gene expression, here we explored the post-transcriptional regulation by microRNAs (miRNAs) of the CPB-ADAM10 gene encoding the Cry3Aa toxin functional receptor ADAM10. RESULTS: The lower CPB-ADAM10 gene expression in CPB larvae fed on potato plants cv. Vivaldi than those fed on potato cv. Monalisa or tomato plants was inversely related to Cry3Aa toxicity. By high-throughput sequencing we identified seven CPB miRNAs and one potato miRNA predicted to base pair with the CPB-ADAM10 messenger RNA. No differential expression of the endogenous lde-miR1175-5p was found in larvae feeding on any of the two potato plant varieties. However, statistically significant increased amounts of potato stu-miR171c-5p were detected in CPB larvae fed on potato cv. Vivaldi compared to larvae fed on potato cv. Monalisa. CONCLUSION: Our results support a role for dietary miRNAs in Bt toxicity by regulating the CPB-ADAM10 gene encoding the Cry3Aa toxin receptor ADAM10 in CPB larvae and opening up the possibility of exploiting plant natural variation in miRNAs to provide more sustainable potato crop protection against CPB. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Dual Antimicrobial and Antiproliferative Activity of TcPaSK Peptide Derived from a Tribolium castaneum Insect Defensin.

Antimicrobial peptides (AMPs) found in the innate immune system of a wide range of organisms might prove useful to fight infections, due to the reported slower development of resistance to AMPs. Increasing the cationicity and keeping moderate hydrophobicity of the AMPs have been described to improve antimicrobial activity. We previously found a peptide derived from the Tribolium castaneum insect defensin 3, exhibiting antrimicrobial activity against several human pathogens. Here, we analyzed the effect against Staphyloccocus aureus of an extended peptide (TcPaSK) containing two additional amino acids, lysine and asparagine, flanking the former peptide fragment in the original insect defensin 3 protein. TcPaSK peptide displayed higher antimicrobial activity against S. aureus, and additionally showed antiproliferative activity against the MDA-MB-231 triple negative breast cancer cell line. A SWATH proteomic analysis revealed the downregulation of proteins involved in cell growth and tumor progression upon TcPaSK cell treatment. The dual role of TcPaSK peptide as antimicrobial and antiproliferative agent makes it a versatile molecule that warrants exploration for its use in novel therapeutic developments as an alternative approach to overcome bacterial antibiotic resistance and to increase the efficacy of conventional cancer treatments.

The Histone Marks Signature in Exonic and Intronic Regions Is Relevant in Early Response of Tomato Genes to Botrytis cinerea and in miRNA Regulation.

Research into the relationship between epigenetic regulation and resistance to biotic stresses provides alternatives for plant protection and crop improvement. To unravel the mechanisms underlying tomato responses to Botrytis cinerea, we performed a chromatin immunoprecipitation (ChIP) analysis showing the increase in H3K9ac mark along the early induced genes SlyDES, SlyDOX1, and SlyLoxD encoding oxylipin-pathway enzymes, and SlyWRKY75 coding for a transcriptional regulator of hormonal signaling. This histone mark showed a more distinct distribution than the previously studied H3K4me3. The RNAPol-ChIP analysis reflected the actual gene transcription associated with increased histone modifications. A different pattern of marks in the oxylipin-related genes against P. syringae supported a pathogen-specific profile, while no significant differences occurred in SlyWRKY75. The epigenetic regulation of SlyWRKY75 by the intron-binding miR1127-3p was supported by the presence of SlyWRKY75 pre-mRNA in control plants. Interestingly, mRNA was found to be accumulated in response to B. cinerea and P. syringae, while reduction in miRNA only occurred against B. cinerea. The intronic region presented a similar pattern of marks than the rest of the gene in both pathosystems, except for H3K4me3 in the miRNA binding site upon B. cinerea. We located the gene encoding Sly-miR1127-3p, which presented reduced H3K4me3 on its promoter against B. cinerea.

Proteomic insights into the immune response of the Colorado potato beetle larvae challenged with Bacillus thuringiensis.

Bacillus thuringiensis (Bt) toxins constitute effective, environmentally safe biopesticides. Nevertheless, insects' tolerance to Bt is influenced by environmental factors affecting immunity. To understand larval immune response in the devastating coleopteran insect pest Colorado potato beetle (CPB), we undertook a proteomic analysis of hemolymph of non-treated control larvae and larvae consuming non-lethal doses of spore-crystal mixtures containing the coleopteran-active Cry3Aa toxin. Results revealed lower amount of proteins involved in insect growth and higher amount of immune response-related proteins in challenged insects, sustaining the larval weight loss observed. Additionally, we found a potential regulatory role of the evolutionary conserved miR-8 in the insect's immune response relying on antimicrobial peptides (AMPs) production. Upon toxin challenge, different patterns of hemolymph AMPs expression and phenoloxidase activity were observed in CPB larvae reared on different Solanaceae plants. This suggests that diet and diet-associated insect midgut microbiota might modulate this insects' tolerance to non-lethal doses of Bt.

Expression of miR159 Is Altered in Tomato Plants Undergoing Drought Stress.

In a scenario of global climate change, water scarcity is a major threat for agriculture, severely limiting crop yields. Therefore, alternatives are urgently needed for improving plant adaptation to drought stress. Among them, gene expression reprogramming by microRNAs (miRNAs) might offer a biotechnologically sound strategy. Drought-responsive miRNAs have been reported in many plant species, and some of them are known to participate in complex regulatory networks via their regulation of transcription factors involved in water stress signaling. We explored the role of miR159 in the response of Solanum lycopersicum Mill. plants to drought stress by analyzing the expression of sly-miR159 and its target SlMYB transcription factor genes in tomato plants of cv. Ailsa Craig grown in deprived water conditions or in response to mechanical damage caused by the Colorado potato beetle, a devastating insect pest of Solanaceae plants. Results showed that sly-miR159 regulatory function in the tomato plants response to distinct stresses might be mediated by differential stress-specific MYB transcription factor targeting. sly-miR159 targeting of SlMYB33 transcription factor transcript correlated with accumulation of the osmoprotective compounds proline and putrescine, which promote drought tolerance. This highlights the potential role of sly-miR159 in tomato plants' adaptation to water deficit conditions.

Colorado potato beetle chymotrypsin genes are differentially regulated in larval midgut in response to the plant defense inducer hexanoic acid or the Bacillus thuringiensis Cry3Aa toxin.

When Colorado potato beetle larvae ingested potato plants treated with the plant defense inducer compound hexanoic acid, midgut chymotrypsin enzyme activity increased, and the corresponding chymotrypsin genes were differentially expressed, evidence of the larval digestive proteolytic system's plasticity. We previously reported increased susceptibility to Cry3Aa toxin in larvae fed hexanoic acid treated plants. Here we show that the most expressed chymotrypsin gene in larvae fed hexanoic acid treated plants, CTR6, was dramatically downregulated in Cry3Aa intoxicated larvae. lde-miR-965-5p and lde-miR-9a-5p microRNAs, predicted to target CTR6, might be involved in regulating the response to hexanoic acid but not to Cry3Aa toxin.

Identification of Stress Associated microRNAs in Solanum lycopersicum by High-Throughput Sequencing.

Tomato (Solanum lycopersicum) is one of the most important crops around the world and also a model plant to study response to stress. High-throughput sequencing was used to analyse the microRNA (miRNA) profile of tomato plants undergoing five biotic and abiotic stress conditions (drought, heat, P. syringae infection, B. cinerea infection, and herbivore insect attack with Leptinotarsa decemlineata larvae) and one chemical treatment with a plant defence inducer, hexanoic acid. We identified 104 conserved miRNAs belonging to 37 families and we predicted 61 novel tomato miRNAs. Among those 165 miRNAs, 41 were stress-responsive. Reverse transcription quantitative PCR (RT-qPCR) was used to validate high-throughput expression analysis data, confirming the expression profiles of 10 out of 11 randomly selected miRNAs. Most of the differentially expressed miRNAs were stress-specific, except for sly-miR167c-3p upregulated in B. cinerea and P. syringae infection, sly-newmiR26-3p upregulated in drought and Hx treatment samples, and sly-newmiR33-3p, sly-newmiR6-3p and sly-newmiR8-3p differentially expressed both in biotic and abiotic stresses. From mature miRNAs sequences of the 41 stress-responsive miRNAs 279 targets were predicted. An inverse correlation between the expression profiles of 4 selected miRNAs (sly-miR171a, sly-miR172c, sly-newmiR22-3p and sly-miR167c-3p) and their target genes (Kinesin, PPR, GRAS40, ABC transporter, GDP and RLP1) was confirmed by RT-qPCR. Altogether, our analysis of miRNAs in different biotic and abiotic stress conditions highlight the interest to understand the functional role of miRNAs in tomato stress response as well as their putative targets which could help to elucidate plants molecular and physiological adaptation to stress.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae).

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

Epigenetic regulation of the expression of WRKY75 transcription factor in response to biotic and abiotic stresses in Solanaceae plants.

KEY MESSAGE: SlyWRKY75: gene expression was induced in response to biotic stresses, especially in Botrytis cinerea-infected tomato plants, in which Sly-miR1127-3p is a putative SlyWRKY75 regulator and epigenetic marks were detected. WRKY75 transcription factor involved in Pi homeostasis was recently found also induced in defense against necrotrophic pathogens. In this study, we analyzed by RT-qPCR the expression of SlyWRKY75 gene in tomato plants in response to abiotic stresses (drought or heat) and biotic stresses (Colorado potato beetle larvae infestation, Pseudomonas syringae or Botrytis cinerea infection) being only differentially expressed following biotic stresses, especially upon B. cinerea infection (55-fold induction). JA and JA-Ile levels were significantly increased in tomato plants under biotic stresses compared with control plants, indicating that SlyWRKY75 might be a transcriptional regulator of the JA pathway. The contribution of miRNAs and epigenetic molecular mechanisms to the regulation of this gene in B. cinerea-infected tomato plants was explored. We identified a putative Sly-miR1127-3p miRNA predicted to bind the intronic region of the SlyWRKY75 genomic sequence. Sly-miR1127-3p miRNA was repressed in infected plants (0.4-fold) supporting that it might act as an epigenetic regulation factor of SlyWRKY75 gene expression rather than via the post-transcriptional mechanisms of canonical miRNAs. It has been proposed that certain miRNAs can mediate DNA methylation in the plant nucleus broadening miRNA functions with transcriptional gene silencing by targeting intron-containing pre-mRNAs. Histone modifications analysis by chromatin immunoprecipitation (ChIP) demonstrated the presence of the activator histone modification H3K4me3 on SlyWRKY75 transcription start site and gene body. The induction of this gene in response to B. cinerea correlates with the presence of an activator mark. Thus, miRNAs and chromatin modifications might cooperate as epigenetic factors to modulate SlyWRKY75 gene expression.

Oxylipin mediated stress response of a miraculin-like protease inhibitor in Hexanoic acid primed eggplant plants infested by Colorado potato beetle.

Insect-plant interactions are governed by a complex equilibrium between the mechanisms through which plant recognize insect attack and orchestrate downstream signaling events that trigger plant defense responses, and the mechanisms by which insects overcome plant defenses. Due to this tight and dynamic interplay, insight into the nature of the plant defense response can be gained by analyzing changes in the insect herbivores digestive system upon plant feeding. In this work we have identified a Solanum melongena miraculin-like protease inhibitor in the midgut juice of Colorado potato larvae feeding on eggplant plants treated with the natural inducer of plant defenses hexanoic acid. We analyzed the corresponding gene expression by qRT-PCR and our results showed that this eggplant miraculin-like gene enhanced induction contributes to the hexanoic acid priming effect in this Solanaceae species. Moreover, our data evidencing that OPDA might be involved in this gene regulation highlights its potential as biomarker in eggplant plant responses to stress mediated this oxylipin signaling pathway.

The Hammerhead Ribozyme: A Long History for a Short RNA.

Small nucleolytic ribozymes are a family of naturally occurring RNA motifs that catalyse a self-transesterification reaction in a highly sequence-specific manner. The hammerhead ribozyme was the first reported and the most extensively studied member of this family. However, and despite intense biochemical and structural research for three decades since its discovery, the history of this model ribozyme seems to be far from finished. The hammerhead ribozyme has been regarded as a biological oddity typical of small circular RNA pathogens of plants. More recently, numerous and new variations of this ribozyme have been found to inhabit the genomes of organisms from all life kingdoms, although their precise biological functions are not yet well understood.

Combining hexanoic acid plant priming with Bacillus thuringiensis insecticidal activity against Colorado potato beetle.

Interaction between insect herbivores and host plants can be modulated by endogenous and exogenous compounds present in the source of food and might be successfully exploited in Colorado potato beetle (CPB) pest management. Feeding tests with CPB larvae reared on three solanaceous plants (potato, eggplant and tomato) resulted in variable larval growth rates and differential susceptibility to Bacillus thuringiensis Cry3Aa toxin as a function of the host plant. An inverse correlation with toxicity was observed in Cry3Aa proteolytic patterns generated by CPB midgut brush-border membrane vesicles (BBMV) from Solanaceae-fed larvae, being the toxin most extensively proteolyzed on potato, followed by eggplant and tomato. We found that CPB cysteine proteases intestains may interact with Cry3Aa toxin and, in CPB BBMV from larvae fed all three Solanaceae, the toxin was able to compete for the hydrolysis of a papain substrate. In response to treatment with the JA-dependent plant inducer Hexanoic acid (Hx), we showed that eggplant reduced OPDA basal levels and both, potato and eggplant induced JA-Ile. CPB larvae feeding on Hx-induced plants exhibited enhanced Cry3Aa toxicity, which correlated with altered papain activity. Results indicated host-mediated effects on B. thuringiensis efficacy against CPB that can be enhanced in combination with Hx plant induction.

Intronic hammerhead ribozymes in mRNA biogenesis.

Small self-cleaving ribozymes are a group of natural RNAs that are capable of catalyzing their own and sequence-specific endonucleolytic cleavage. One of the most studied members is the hammerhead ribozyme (HHR), a catalytic RNA originally discovered in subviral plant pathogens but recently shown to reside in a myriad of genomes along the tree of life. In eukaryotes, most of the genomic HHRs seem to be related to short interspersed retroelements, with the main exception of a group of strikingly conserved ribozymes found in the genomes of all amniotes (reptiles, birds and mammals). These amniota HHRs occur in the introns of a few specific genes, and clearly point to a preserved biological role during pre-mRNA biosynthesis. More specifically, bioinformatic analysis suggests that these intronic ribozymes could offer a new form of splicing regulation of the mRNA of higher vertebrates. We review here the latest advances in the discovery and biological characterization of intronic HHRs of vertebrates, including new conserved examples in the genomes of the primitive turtle and coelacanth fish.

Functional significance of membrane associated proteolysis in the toxicity of Bacillus thuringiensis Cry3Aa toxin against Colorado potato beetle.

Bacillus thuringiensis Cry toxins are widely used as biocontrol agents in bioinsecticides and transgenic plants. In the three domain-Cry toxins, domain II has been identified as an important determinant of their highly specific activity against insects. In this work, we assessed the role in membrane associated proteolysis and toxicity in Colorado potato beetle (CPB) of a previously reported ADAM recognition motif present in Cry3Aa toxin domain II. We used site-directed mutagenesis to modify the Bacillus thuringiensis cry3A gene in amino acid residues 344, 346, 347, 351 and 353 of the ADAM recognition motif in Cry3Aa toxin. Cry3Aa toxin mutants displayed decreased toxicity when compared to the wild type toxin and impaired ability to compete CPB brush border membrane associated cleavage of an ADAM fluorogenic substrate. Although the proteolytic profile of Cry3Aa toxin mutants generated by brush border membrane associated proteases was similar to that of Cry3Aa toxin, the metalloprotease inhibitor 1,10-phenanthroline was less efficient on the proteolysis of mutants than on that of the wild type toxin. The relevance of the Cry3Aa-ADAM interaction through the predicted recognition sequence was further confirmed by analyzing the effect of membrane integrity disturbance on Cry3Aa toxin membrane associated proteolysis and CPB larvae toxicity. Data support that Cry3Aa proteolysis, as a result of the interaction with ADAM through the Cry3Aa recognition motif, is essential for Cry3Aa toxic action in CPB. Detailed knowledge of Cry3Aa interaction with CPB midgut membrane should facilitate the development of more effective Bt based products against this devastating pest and other Coleoptera.

Identification of a calmodulin-binding site within the domain I of Bacillus thuringiensis Cry3Aa toxin.

Bacillus thuringiensis Cry3Aa toxin is a coleopteran specific toxin highly active against Colorado Potato Beetle (CPB).We have recently shown that Cry3Aa toxin is proteolytically cleaved by CPB midgut membrane associated metalloproteases and that this cleavage is inhibited by ADAM metalloprotease inhibitors. In the present study, we investigated whether the Cry3Aa toxin is a calmodulin (CaM) binding protein, as it is the case of several different ADAM shedding substrates. In pull-down assays using agarose beads conjugated with CaM, we demonstrated that Cry3Aa toxin specifically binds to CaM in a calcium-independent manner. Furthermore, we used gel shift assays and (1)H NMR spectra to demonstrate that CaM binds to a 16-amino acid synthetic peptide corresponding to residues N256-V271 within the domain I of Cry3Aa toxin. Finally, to investigate whether CaM has any effect on Cry3Aa toxin CPB midgut membrane associated proteolysis, cleavage assays were performed in the presence of the CaM-specific inhibitor trifluoperazine. We showed that trifluoperazine significantly increased Cry3Aa toxin proteolysis and also decreased Cry3Aa larval toxicity.

Ubiquitous presence of the hammerhead ribozyme motif along the tree of life.

Examples of small self-cleaving RNAs embedded in noncoding regions already have been found to be involved in the control of gene expression, although their origin remains uncertain. In this work, we show the widespread occurrence of the hammerhead ribozyme (HHR) motif among genomes from the Bacteria, Chromalveolata, Plantae, and Metazoa kingdoms. Intergenic HHRs were detected in three different bacterial genomes, whereas metagenomic data from Galapagos Islands showed the occurrence of similar ribozymes that could be regarded as direct relics from the RNA world. Among eukaryotes, HHRs were detected in the genomes of three water molds as well as 20 plant species, ranging from unicellular algae to vascular plants. These HHRs were very similar to those previously described in small RNA plant pathogens and, in some cases, appeared as close tandem repetitions. A parallel situation of tandemly repeated HHR motifs was also detected in the genomes of lower metazoans from cnidarians to invertebrates, with special emphasis among hematophagous and parasitic organisms. Altogether, these findings unveil the HHR as a widespread motif in DNA genomes, which would be involved in new forms of retrotransposable elements.

Intronic hammerhead ribozymes are ultraconserved in the human genome.

Small ribozymes have been regarded as living fossils of a prebiotic RNA world that would have remained in the genomes of modern organisms. In this study, we report the ultraconserved occurrence of hammerhead ribozymes in Amniota genomes (reptiles, birds and mammals, including humans), similar to those described previously in amphibians and platyhelminth parasites. The ribozymes mapped to intronic regions of different genes, such as the tumour suppressor RECK in birds and mammals, a mammalian tumour antigen and the dystrobrevin beta in lizards and birds. In vitro characterization confirmed a high self-cleavage activity, whereas analysis of RECK-expressed sequence tags revealed fusion events between the in vivo self-cleaved intron and U5 or U6 small nuclear RNA fragments. Together, these results suggest a conserved role for these ribozymes in messenger RNA biogenesis.

The role of positive selection in hepatitis C virus.

Hepatitis C virus (HCV) is a major health problem worldwide, infecting an estimated 170 million people. In this study, we have employed a large data set of sequences (14,654 sequences from between 25 and 100 clone sequences per analyzed region and per patient) from 67 patients infected with HCV genotype 1 (23 subtype 1a and 44 subtype 1b). For all patients, a sample prior to combined therapy with alpha interferon plus ribavirin was available, whereas for some patients additional samples after 6 or 12 months of treatment were also available. Twenty-seven patients responded to treatment (12 subtype 1a and 15 subtype 1b) and forty patients did not respond to treatment (11 subtype 1a vs. 29 subtype 1b). Two regions of the HCV genome were analyzed, one compressing the hypervariable regions (HVR1, HVR2 and HVR3) of the envelope 2 glycoprotein and another one including the interferon sensitive determining region (ISDR) and the V3 domain of the NS5A protein. Previously (Cuevas, J.M., Torres-Puente, M., Jiménez-Hernández, N., Bracho, M.A., García-Robles, I., Wrobel, B., Carnicer, F., del Olmo, J., Ortega, E., Moya, A., González-Candelas, F., 2008b. Genetic variability of hepatitis C virus before and after combined therapy of interferon plus ribavirin. Plos One 3 (8), e3058), several amino acid positions in both regions analyzed were detected to be under positive selection. Here, we have compared the amino acid composition of each positively selected position between responder and non-responder patients for both subtypes. If we exclude some non-conclusive cases, no clear differences were detected in any case. In conclusion, identifying specific positions as completely discriminatory of treatment response seems to be a difficult task. Our results, in concordance with previous studies, suggest that HCV evasion strategies are more likely based on a global increased variability, which would yield combinations of mutations with an increased resistance, than on the fixation of specific amino acids conferring resistance to antiviral treatment or immune response. In this sense, the particular systemic response from each patient could play an essential role in determining the outcome of the antiviral treatment.

Combined therapy of interferon plus ribavirin promotes multiple adaptive solutions in hepatitis C virus.

Hepatitis C virus (HCV) presents several regions involved potentially in evading antiviral treatment and host immune system. Two regions, known as PKR-BD and V3 domains, have been proposed to be involved in resistance to interferon. Additionally, hypervariable regions in the envelope E2 glycoprotein are also good candidates to participate in evasion from the immune system. In this study, we have used a cohort of 22 non-responder patients to combined therapy (interferon alpha-2a plus ribavirin) for which samples obtained just before initiation of therapy and after 6 or/and 12 months of treatment were available. A range of 25-100 clones per patient, genome region and time sample were obtained. The predominant amino acid sequences for each time sample and patient were determined. Next, the sequences of the PKR-BD and V3 domains and the hypervariable regions from different time samples were compared for each patient. The highest levels of variability were detected at the three hypervariable regions of the E2 protein and, to a lower extent, at the V3 domain of the NS5A protein. However, no clear patterns of adaptation to the host immune system or to antiviral treatment were detected. In summary, although high levels of variability are correlated to viral adaptive response, antiviral treatment does not seem to promote convergent adaptive changes. Consequently, other regions must be involved in evasion strategies likely based on a combination of multiple mechanisms, in which pools of changes along the HCV genome could confer viruses the ability to overcome strong selective pressures.