High-depth sequencing characterization of viral dynamics across tissues in fatal COVID-19 reveals compartmentalized infection

Fulminant myocarditis induced by SARS-CoV-2 infection without severe lung involvement: insights into COVID-19 pathogenesis

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection often leads to pulmonary complications. Cardiovascular sequelae, including myocarditis and heart failure, have also been reported. Here, the study presents two fulminant myocarditis cases infected by SARS-CoV-2 exhibiting remarkable elevation of cardiac biomarkers without significant pulmonary injury, as determined by imaging examinations. Immunohistochemical staining reveals viral antigen within cardiomyocytes, indicating that SARS-CoV-2 could directly infect myocardium. The full viral genomes from respiratory, anal, and myocardial specimens are obtained via next-generation sequencing. Phylogenetic analyses of the whole genome and spike gene indicate that viruses in the myocardium/pericardial effusion and anal swabs are closely related and cluster together yet diverge from those in the respiratory samples. In addition, unique mutations are found in the anal/myocardial strains compared to the respiratory strains, suggesting tissue-specific virus mutation and adaptation. These findings indicate genetically distinct SARS-CoV-2 variants have infiltrated and disseminated within myocardial tissues, independent of pulmonary injury, and point to different infection routes between the myocardium and respiratory tract, with myocardial infections potentially arising from intestinal infection. These findings highlight the potential for systemic SARS-CoV-2 infection and the importance of a thorough multi-organ assessment in patients for a comprehensive understanding of the pathogenesis of COVID-19.

Persistent SARS-CoV-2 infection: significance and implications

SARS-CoV-2 causes persistent infections in a subset of individuals, which is a major clinical and public health problem that should be prioritised for further investigation for several reasons. First, persistent SARS-CoV-2 infection often goes unrecognised, and therefore might affect a substantial number of people, particularly immunocompromised individuals. Second, the formation of tissue reservoirs (including in non-respiratory tissues) might underlie the pathophysiology of the persistent SARS-CoV-2 infection and require new strategies for diagnosis and treatment. Finally, persistent SARS-CoV-2 replication, particularly in the setting of suboptimal immune responses, is a possible source of new, divergent virus variants that escape pre-existing immunity on the individual and population levels. Defining optimal diagnostic and treatment strategies for patients with persistent virus replication and monitoring viral evolution are therefore urgent medical and public health priorities.

Natural history of Ebola virus disease in rhesus monkeys shows viral variant emergence dynamics and tissue-specific host responses

Ebola virus (EBOV) causes Ebola virus disease (EVD), marked by severe hemorrhagic fever; however, the mechanisms underlying the disease remain unclear. To assess the molecular basis of EVD across time, we performed RNA sequencing on 17 tissues from a natural history study of 21 rhesus monkeys, developing new methods to characterize host-pathogen dynamics. We identified alterations in host gene expression with previously unknown tissue-specific changes, including downregulation of genes related to tissue connectivity. EBOV was widely disseminated throughout the body; using a new, broadly applicable deconvolution method, we found that viral load correlated with increased monocyte presence. Patterns of viral variation between tissues differentiated primary infections from compartmentalized infections, and several variants impacted viral fitness in a EBOV/Kikwit minigenome system, suggesting that functionally significant variants can emerge during early infection. This comprehensive portrait of host-pathogen dynamics in EVD illuminates new features of pathogenesis and establishes resources to study other emerging pathogens.

Differential activation of programmed cell death in patients with severe SARS-CoV-2 infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes severe lower airway disease and death in a subset of patients. Knowledge on the relative contribution of programmed cell death (PCD) to lung pathology is limited to few human autopsy studies with small sample size/scope, in vitro cell culture, and experimental model systems. In this study, we sought to identify, localize, and quantify activation of apoptosis, ferroptosis, pyroptosis, and necroptosis in FFPE lung tissues from patients that died from severe SARS-CoV-2 infection (n = 28) relative to uninfected controls (n = 13). Immunofluorescence (IF) staining, whole-slide imaging, and Image J software was used to localize and quantify expression of SARS-CoV-2 nucleoprotein and the following PCD protein markers: cleaved Caspase-3, pMLKL, cleaved Gasdermin D, and CD71, respectively. IF showed differential activation of each PCD pathway in infected lungs and dichotomous staining for SARS-CoV-2 nucleoprotein enabling distinction between high (n = 9) vs low viral burden (n = 19). No differences were observed in apoptosis and ferroptosis in SARS-CoV-2 infected lungs relative to uninfected controls. However, both pyroptosis and necroptosis were significantly increased in SARS-CoV-2-infected lungs. Increased pyroptosis was observed in SARS-CoV-2 infected lungs, irrespective of viral burden, suggesting an inflammation-driven mechanism. In contrast, necroptosis exhibited a very strong positive correlation with viral burden (R2 = 0.9925), suggesting a direct SARS-CoV-2 mediated effect. These data indicate a possible novel mechanism for viral-mediated necroptosis and a potential role for both lytic programmed cell death pathways, necroptosis and pyroptosis, in mediating infection outcome.

A Simulation Framework for Modeling the Within-Patient Evolutionary Dynamics of SARS-CoV-2

The global impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for rarely acting positive selection are best performed via comparison of empirical data with simulated data wherein commonly acting evolutionary factors, including mutation and recombination, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. Although there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intrahost evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them with existing empirical data. Of these, 592 models (∼5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intrahost SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed toward strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

From viral democratic genomes to viral wild bunch of quasispecies

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

SARS-CoV-2 Bottlenecks and Tissue-Specific Adaptation in the Central Nervous System

Severe COVID-19 and post-acute sequelae of SARS-CoV-2 infection are associated with neurological complications that may be linked to direct infection of the central nervous system (CNS), but the selective pressures ruling neuroinvasion are poorly defined. Here, we assessed SARS-CoV-2 evolution in the lung versus CNS of infected mice. Higher levels of viral diversity were observed in the CNS than the lung after intranasal challenge with a high frequency of mutations in the Spike furin cleavage site (FCS). Deletion of the FCS significantly attenuated virulence after intranasal challenge, with lower viral titers and decreased morbidity compared to the wild-type virus. Intracranial inoculation of the FCS-deleted virus, however, was sufficient to restore virulence. After intracranial inoculation, both viruses established infection in the lung, but this required reversion of the FCS deletion. Cumulatively, these data suggest a critical role for the FCS in determining SARS-CoV-2 tropism and compartmentalization with possible implications for the treatment of neuroinvasive COVID-19.

Intra-host variation in the spike S1/S2 region of a feline coronavirus type-1 in a cat with persistent infection

Feline coronavirus type 1 (FCoV-1) is widely known for causing feline infectious peritonitis (FIP), a systemic infection that is often fatal, with the virus known as the FIPV biotype. However, subclinical disease also occurs, in which cats may not show signs and intermittently shed the virus, including in feces, possibly for long periods of time. This virus is known as the FECV biotype. Progression of FECV to FIPV has been linked to several genomic changes, however a specific region of the viral spike protein at the interface of the spike S1 and S2 domains has been especially implicated. In this study, we followed a cat (#576) for six years from 2017, at which time FCoV-1 was detected in feces and conjunctival swabs, until 2022, when the animal was euthanized based on a diagnosis of alimentary small cell lymphoma. Over this time period, the cat was clinically diagnosed with inflammatory bowel disease and chronic rhinitis, and cardiac problems were also suspected. Using hybridization capture targeting the spike (S) gene of FCoV followed by next-generation sequencing, we screened 27 clinical samples. We detected FCoV-1 in 4 samples taken in 2017 (intestine and nasal tissue, feces, and conjunctiva), and 3 samples taken in 2022 (feces, and intestinal and heart tissue), but not in fecal samples taken in 2019 and 2020. Next, we focused on the S1/S2 region within S, which contains the furin cleavage site (FCS), a key regulator of viral transmission and pathogenesis. We show that the FCoV-1 variants obtained from feces in 2017 and 2022 were identical, while the ones from conjunctiva (2017), heart (2022), and intestine (2017 and 2022) were distinct. Sequence comparison of all the variants obtained showed that most of the non-synonymous changes in the S1/S2 region occur within the FCS. In the heart, we found two variants that differed by a single nucleotide, resulting in distinct FCS motifs that differ in one amino acid. It is predicted that one of these FCS motifs will down-regulate spike cleavability. The variant from the conjunctiva (2017) had a 6-nucleotide in-frame insertion that resulted in a longer and more exposed S1/S2 loop, which is predicted to be more accessible to the furin protease. Our studies indicate that FCoV-1 can independently persist in the gastrointestinal tract and heart of a cat over a long period of time without evidence of typical FIP signs, with intermittent viral shedding from the gastrointestinal and respiratory tracts.

A simulation framework for modeling the within-patient evolutionary dynamics of SARS-CoV-2

The global impact of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for positive selection are best performed via comparison of empirical data to simulated data wherein evolutionary factors, including mutation and recombination rates, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. While there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intra-host evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them to existing empirical data. Of these, 592 models (~5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intra-host SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed towards strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

Differential Activation of Programmed Cell Death in Patients with Severe SARS-CoV-2 Infection

SARS-CoV-2 (SARS-2) causes severe lower airway disease and death in a subset of patients. Knowledge on the relative contribution of programmed cell death (PCD) to lung pathology is limited to few human autopsy studies with small sample size/scope, in vitro cell culture, and experimental model systems. In this study, we sought to identify, localize, and quantify activation of apoptosis, ferroptosis, pyroptosis, and necroptosis in FFPE lung tissues from patients that died from severe SARS-2 infection (n=28) relative to uninfected controls (n=13). Immunofluorescence (IF) staining, whole-slide imaging, and Image J software was used to localize and quantify expression of SARS-2 nucleoprotein and the following PCD protein markers: cleaved Caspase-3, pMLKL, cleaved Gasdermin D, and CD71, respectively. IF showed differential activation of each PCD pathway in SARS-2 infected lungs and dichotomous staining for SARS-2 nucleoprotein enabling distinction between high (n=9) vs low viral burden (n= 19). No differences were observed in apoptosis and ferroptosis in SARS-2 infected lungs relative to uninfected controls. However, both pyroptosis and necroptosis were significantly increased in SARS-2 infected lungs. Increased pyroptosis was observed in SARS-2 infected lungs, irrespective of viral burden, suggesting an inflammation-driven mechanism. In contrast, necroptosis exhibited a very strong positive correlation with viral burden (R2=0.9925), suggesting a direct SARS-2 mediated effect. These data indicate a possible novel mechanism for viral-mediated necroptosis and a potential role for both lytic programmed cell death pathways, necroptosis and pyroptosis, in mediating infection outcome.