Aspergillus fumigatus Can Display Persistence to the Fungicidal Drug Voriconazole.

Aspergillus fumigatus is a filamentous fungus that can infect the lungs of patients with immunosuppression and/or underlying lung diseases. The mortality associated with chronic and invasive aspergillosis infections remain very high, despite availability of antifungal treatments. In the last decade, there has been a worrisome emergence and spread of resistance to the first-line antifungals, the azoles. The mortality caused by resistant isolates is even higher, and patient management is complicated as the therapeutic options are reduced. Nevertheless, treatment failure is also common in patients infected with azole-susceptible isolates, which can be due to several non-mutually exclusive reasons, such as poor drug absorption. In addition, the phenomena of tolerance or persistence, where susceptible pathogens can survive the action of an antimicrobial for extended periods, have been associated with treatment failure in bacterial infections, and their occurrence in fungal infections already proposed. Here, we demonstrate that some isolates of A. fumigatus display persistence to voriconazole. A subpopulation of the persister isolates can survive for extended periods and even grow at low rates in the presence of supra-MIC of voriconazole and seemingly other azoles. Persistence cannot be eradicated with adjuvant drugs or antifungal combinations and seemed to reduce the efficacy of treatment for certain individuals in a Galleria mellonella model of infection. Furthermore, persistence implies a distinct transcriptional profile, demonstrating that it is an active response. We propose that azole persistence might be a relevant and underestimated factor that could influence the outcome of infection in human aspergillosis. IMPORTANCE The phenomena of antibacterial tolerance and persistence, where pathogenic microbes can survive for extended periods in the presence of cidal drug concentrations, have received significant attention in the last decade. Several mechanisms of action have been elucidated, and their relevance for treatment failure in bacterial infections demonstrated. In contrast, our knowledge of antifungal tolerance and, in particular, persistence is still very limited. In this study, we have characterized the response of the prominent fungal pathogen Aspergillus fumigatus to the first-line therapy antifungal voriconazole. We comprehensively show that some isolates display persistence to this fungicidal antifungal and propose various potential mechanisms of action. In addition, using an alternative model of infection, we provide initial evidence to suggest that persistence may cause treatment failure in some individuals. Therefore, we propose that azole persistence is an important factor to consider and further investigate in A. fumigatus.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations.

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes.

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

High SARS-CoV-2 viral load is associated with a worse clinical outcome of COVID-19 disease.

COVID-19 severity and progression are determined by several host and virological factors that may influence the final outcome of SARS-CoV-2-infected patients. The objective of this work was to determine a possible association between viral load, obtained from nasopharyngeal swabs, and the severity of the infection in a cohort of 448 SARS-CoV-2-infected patients from a hospital in Madrid during the first outbreak of the pandemic in Spain. To perform this, we clinically classified patients as mild, moderate and severe COVID-19 according to a number of clinical parameters such as hospitalization requirement, need of oxygen therapy, admission to intensive care units and/or death. Also, Ct values were determined using SARS-CoV-2-specific oligonucleotides directed to ORF1ab. Here we report a statistically significant association between viral load and disease severity, a high viral load being associated with worse clinical prognosis, independently of several previously identified risk factors such as age, sex, hypertension, cardiovascular disease, diabetes, obesity and lung disease (asthma and chronic obstructive pulmonary disease). The data presented here reinforce viral load as a potential biomarker for predicting disease severity in SARS-CoV-2-infected patients. It is also an important parameter in viral evolution since it relates to the numbers and types of variant genomes present in a viral population, a potential determinant of disease progression.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics.

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus.

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium-revealed by the changing composition of the mutant spectrum-may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.

Amino Acid Substitutions Associated with Treatment Failure for Hepatitis C Virus Infection.

Despite the high virological response rates achieved with current directly acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. The identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultradeep sequencing (UDS) methods, for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide resistance-associated substitutions (RAS), from 220 patients who failed therapy. They were present frequently in basal and posttreatment virus of patients who failed different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. They also have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value for treatment response, are discussed.