SARS-CoV-2 evolution during prolonged infection in immunocompromised patients

Prolonged infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in immunocompromised patients provides an opportunity for viral evolution, potentially leading to the generation of new pathogenic variants. To investigate the pathways of viral evolution, we carried out a study on five patients experiencing prolonged SARS-CoV-2 infection (quantitative polymerase chain reaction-positive for 79-203 days) who were immunocompromised due to treatment for lymphoma or solid organ transplantation. For each timepoint analyzed, we generated at least two independent viral genome sequences to assess the heterogeneity and control for sequencing error. Four of the five patients likely had prolonged infection; the fifth apparently experienced a reinfection. The rates of accumulation of substitutions in the viral genome per day were higher in hospitalized patients with prolonged infection than those estimated for the community background. The spike coding region accumulated a significantly greater number of unique mutations than other viral coding regions, and the mutation density was higher. Two patients were treated with monoclonal antibodies (bebtelovimab and sotrovimab); by the next sampled timepoint, each virus population showed substitutions associated with monoclonal antibody resistance as the dominant forms (spike K444N and spike E340D). All patients received remdesivir, but remdesivir-resistant substitutions were not detected. These data thus help elucidate the trends of emergence, evolution, and selection of mutational variants within long-term infected immunocompromised individuals.

IMPORTANCE

SARS-CoV-2 is responsible for a global pandemic, driven in part by the emergence of new viral variants. Where do these new variants come from? One model is that long-term viral persistence in infected individuals allows for viral evolution in response to host pressures, resulting in viruses more likely to replicate efficiently in humans. In this study, we characterize replication in several hospitalized and long-term infected individuals, documenting efficient pathways of viral evolution.

Microbiome Profiling Demonstrates Concordance of Endotracheal Tube Aspirates With Direct Lower Airway Sampling in Intubated Patients

BACKGROUND

Endotracheal aspirates (ETAs) are widely used for microbiologic studies of the respiratory tract in intubated patients. However, they involve sampling through an established endotracheal tube using suction catheters, both of which can acquire biofilms that may confound results.

RESEARCH QUESTION

Does standard clinical ETA in intubated patients accurately reflect the authentic lower airway bacterial microbiome?

STUDY DESIGN AND METHODS

Comprehensive quantitative bacterial profiling using 16S rRNA V1-V2 gene sequencing was applied to compare bacterial populations captured by standard clinical ETA vs contemporaneous gold standard samples acquired directly from the lower airways through a freshly placed sterile tracheostomy tube. The study included 13 patients undergoing percutaneous tracheostomy following prolonged (median, 15 days) intubation. Metrics of bacterial composition, diversity, and relative quantification were applied to samples.

RESULTS

Pre-tracheostomy ETAs closely resembled the gold standard immediate post-tracheostomy airway microbiomes in bacterial composition and community features of diversity and quantification. Endotracheal tube and suction catheter biofilms also resembled cognate ETA and fresh tracheostomy communities.

INTERPRETATION

Unbiased molecular profiling shows that standard clinical ETA sampling has good concordance with the authentic lower airway microbiome in intubated patients.

Tracking gut microbiome and bloodstream infection in critically ill adults

BACKGROUND

The gut microbiome is believed to contribute to bloodstream infection (BSI) via translocation of dominant gut bacteria in vulnerable patient populations. However, conclusively linking gut and blood organisms requires stringent approaches to establish strain-level identity.

METHODS

We enrolled a convenience cohort of critically ill patients and investigated 86 bloodstream infection episodes that occurred in 57 patients. Shotgun metagenomic sequencing was used to define constituents of their gut microbiomes, and whole genome sequencing and assembly was done on 23 unique bloodstream isolates that were available from 21 patients. Whole genome sequences were downloaded from public databases and used to establish sequence-identity distribution and define thresholds for unrelated genomes of BSI species. Gut microbiome reads were then aligned to whole genome sequences of the cognate bloodstream isolate and unrelated database isolates to assess identity.

RESULTS

Gut microbiome constituents matching the bloodstream infection species were present in half of BSI episodes, and represented >30% relative abundance of gut sequences in 10% of episodes. Among the 23 unique bloodstream organisms that were available for whole genome sequencing, 14 were present in gut at the species level. Sequence alignment applying defined thresholds for identity revealed that 6 met criteria for identical strains in blood and gut, but 8 did not. Sequence identity between BSI isolates and gut microbiome reads was more likely when the species was present at higher relative abundance in gut.

CONCLUSION

In assessing potential gut source for BSI, stringent sequence-based approaches are essential to determine if organisms responsible for BSI are identical to those in gut: of 14 evaluable patients in which the same species was present in both sites, they were identical in 6/14, but were non-identical in 8/14 and thus inconsistent with gut source. This report demonstrates application of sequencing as a key tool to investigate infection tracking within patients.

Effects of Mask Reuse on the Oropharyngeal, Skin, and Mask Microbiome

BACKGROUND

Face masks have been critical in the coronavirus disease 2019 (COVID-19) pandemic, but supplies were sometimes limited and disposable masks contribute greatly to environmental waste. Studies suggest that filtration capacity is retained with repeated use, and surveys indicate many people reuse surgical masks. However, the impact of mask reuse on the host is understudied.

METHODS

We applied 16S rRNA gene sequencing to investigate the bacterial microbiome of the facial skin and oropharynx of individuals randomized to wearing fresh surgical masks daily versus masks reused for 1 week.

RESULTS

Compared to daily fresh masks, reuse was associated with increased richness (number of taxa) of the skin microbiome and trend towards greater diversity, but no difference in the oropharyngeal microbiome. Used masks had either skin-dominant or oropharynx-dominant bacterial sequences, and reused masks had >100-fold higher bacterial content but no change in composition compared to those used for 1 day.

CONCLUSIONS

One week of mask reuse increased the number of low-abundance taxa on the face but did not impact the upper respiratory microbiome. Thus, face mask reuse has little impact on the host microbiome, although whether minor changes to the skin microbiome might relate to reported skin sequelae of masking (maskne) remains to be determined.

SARS-CoV-2 variants associated with vaccine breakthrough in the Delaware Valley through summer 2021

The severe acute respiratory coronavirus-2 (SARS-CoV-2) is the cause of the global outbreak of COVID-19. Evidence suggests that the virus is evolving to allow efficient spread through the human population, including vaccinated individuals. Here we report a study of viral variants from surveillance of the Delaware Valley, including the city of Philadelphia, and variants infecting vaccinated subjects. We sequenced and analyzed complete viral genomes from 2621 surveillance samples from March 2020 to September 2021 and compared them to genome sequences from 159 vaccine breakthroughs. In the early spring of 2020, all detected variants were of the B.1 and closely related lineages. A mixture of lineages followed, notably including B.1.243 followed by B.1.1.7 (alpha), with other lineages present at lower levels. Later isolations were dominated by B.1.617.2 (delta) and other delta lineages; delta was the exclusive variant present by the last time sampled. To investigate whether any variants appeared preferentially in vaccine breakthroughs, we devised a model based on Bayesian autoregressive moving average logistic multinomial regression to allow rigorous comparison. This revealed that B.1.617.2 (delta) showed three-fold enrichment in vaccine breakthrough cases (odds ratio of 3; 95% credible interval 0.89-11). Viral point substitutions could also be associated with vaccine breakthroughs, notably the N501Y substitution found in the alpha, beta and gamma variants (odds ratio 2.04; 95% credible interval of 1.25-3.18). This study thus provides a detailed picture of viral evolution in the Delaware Valley and a geographically matched analysis of vaccine breakthroughs; it also introduces a rigorous statistical approach to interrogating enrichment of viral variants.

Signatures of COVID-19 Severity and Immune Response in the Respiratory Tract Microbiome

Viral infection of the respiratory tract can be associated with propagating effects on the airway microbiome, and microbiome dysbiosis may influence viral disease. Here, we investigated the respiratory tract microbiome in coronavirus disease 2019 (COVID-19) and its relationship to disease severity, systemic immunologic features, and outcomes. We examined 507 oropharyngeal, nasopharyngeal, and endotracheal samples from 83 hospitalized COVID-19 patients as well as non-COVID patients and healthy controls. Bacterial communities were interrogated using 16S rRNA gene sequencing, and the commensal DNA viruses Anelloviridae and Redondoviridae were quantified by qPCR. We found that COVID-19 patients had upper respiratory microbiome dysbiosis and greater change over time than critically ill patients without COVID-19. Oropharyngeal microbiome diversity at the first time point correlated inversely with disease severity during hospitalization. Microbiome composition was also associated with systemic immune parameters in blood, as measured by lymphocyte/neutrophil ratios and immune profiling of peripheral blood mononuclear cells. Intubated patients showed patient-specific lung microbiome communities that were frequently highly dynamic, with prominence of Staphylococcus. Anelloviridae and Redondoviridae showed more frequent colonization and higher titers in severe disease. Machine learning analysis demonstrated that integrated features of the microbiome at early sampling points had high power to discriminate ultimate level of COVID-19 severity. Thus, the respiratory tract microbiome and commensal viruses are disturbed in COVID-19 and correlate with systemic immune parameters, and early microbiome features discriminate disease severity. Future studies should address clinical consequences of airway dysbiosis in COVID-19, its possible use as biomarkers, and the role of bacterial and viral taxa identified here in COVID-19 pathogenesis. IMPORTANCE COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection of the respiratory tract, results in highly variable outcomes ranging from minimal illness to death, but the reasons for this are not well understood. We investigated the respiratory tract bacterial microbiome and small commensal DNA viruses in hospitalized COVID-19 patients and found that each was markedly abnormal compared to that in healthy people and differed from that in critically ill patients without COVID-19. Early airway samples tracked with the level of COVID-19 illness reached during hospitalization, and the airway microbiome also correlated with immune parameters in blood. These findings raise questions about the mechanisms linking SARS-CoV-2 infection and other microbial inhabitants of the airway, including whether the microbiome might regulate severity of COVID-19 disease and/or whether early microbiome features might serve as biomarkers to discriminate disease severity.

Signatures of COVID-19 severity and immune response in the respiratory tract microbiome

Rationale

Viral infection of the respiratory tract can be associated with propagating effects on the airway microbiome, and microbiome dysbiosis may influence viral disease.

Objective

To define the respiratory tract microbiome in COVID-19 and relationship disease severity, systemic immunologic features, and outcomes.

Methods and Measurements

We examined 507 oropharyngeal, nasopharyngeal and endotracheal samples from 83 hospitalized COVID-19 patients, along with non-COVID patients and healthy controls. Bacterial communities were interrogated using 16S rRNA gene sequencing, commensal DNA viruses Anelloviridae and Redondoviridae were quantified by qPCR, and immune features were characterized by lymphocyte/neutrophil (L/N) ratios and deep immune profiling of peripheral blood mononuclear cells (PBMC).

Main Results

COVID-19 patients had upper respiratory microbiome dysbiosis, and greater change over time than critically ill patients without COVID-19. Diversity at the first time point correlated inversely with disease severity during hospitalization, and microbiome composition was associated with L/N ratios and PBMC profiles in blood. Intubated patients showed patient-specific and dynamic lung microbiome communities, with prominence of Staphylococcus. Anelloviridae and Redondoviridae showed more frequent colonization and higher titers in severe disease. Machine learning analysis demonstrated that integrated features of the microbiome at early sampling points had high power to discriminate ultimate level of COVID-19 severity.

Conclusions

The respiratory tract microbiome and commensal virome are disturbed in COVID-19, correlate with systemic immune parameters, and early microbiome features discriminate disease severity. Future studies should address clinical consequences of airway dysbiosis in COVID-19, possible use as biomarkers, and role of bacterial and viral taxa identified here in COVID-19 pathogenesis.

SARS-CoV-2 Variants Associated with Vaccine Breakthrough in the Delaware Valley through Summer 2021

The severe acute respiratory coronavirus-2 (SARS-CoV-2) is the cause of the global outbreak of COVID-19. Evidence suggests that the virus is evolving to allow efficient spread through the human population, including vaccinated individuals. Here, we report a study of viral variants from surveillance of the Delaware Valley, including the city of Philadelphia, and variants infecting vaccinated subjects. We sequenced and analyzed complete viral genomes from 2621 surveillance samples from March 2020 to September 2021 and compared them to genome sequences from 159 vaccine breakthroughs. In the early spring of 2020, all detected variants were of the B.1 and closely related lineages. A mixture of lineages followed, notably including B.1.243 followed by B.1.1.7 (alpha), with other lineages present at lower levels. Later isolations were dominated by B.1.617.2 (delta) and other delta lineages; delta was the exclusive variant present by the last time sampled. To investigate whether any variants appeared preferentially in vaccine breakthroughs, we devised a model based on Bayesian autoregressive moving average logistic multinomial regression to allow rigorous comparison. This revealed that B.1.617.2 (delta) showed 3-fold enrichment in vaccine breakthrough cases (odds ratio of 3; 95% credible interval 0.89-11). Viral point substitutions could also be associated with vaccine breakthroughs, notably the N501Y substitution found in the alpha, beta and gamma variants (odds ratio 2.04; 95% credible interval of1.25-3.18). This study thus overviews viral evolution and vaccine breakthroughs in the Delaware Valley and introduces a rigorous statistical approach to interrogating enrichment of breakthrough variants against a changing background. IMPORTANCE SARS-CoV-2 vaccination is highly effective at reducing viral infection, hospitalization and death. However, vaccine breakthrough infections have been widely observed, raising the question of whether particular viral variants or viral mutations are associated with breakthrough. Here, we report analysis of 2621 surveillance isolates from people diagnosed with COVID-19 in the Delaware Valley in southeastern Pennsylvania, allowing rigorous comparison to 159 vaccine breakthrough case specimens. Our best estimate is a 3-fold enrichment for some lineages of delta among breakthroughs, and enrichment of a notable spike substitution, N501Y. We introduce statistical methods that should be widely useful for evaluating vaccine breakthroughs and other viral phenotypes.

SARS-CoV-2 Genomic Variation in Space and Time in Hospitalized Patients in Philadelphia

The severe acute respiratory coronavirus 2 (SARS-CoV-2) is the cause of the global outbreak of COVID-19. The epidemic accelerated in Philadelphia, PA, in the spring of 2020, with the city experiencing a first peak of infections on 15 April, followed by a decline through midsummer. Here, we investigate spread of the epidemic in the first wave in Philadelphia using full-genome sequencing of 52 SARS-CoV-2 samples obtained from 27 hospitalized patients collected between 30 March and 17 July 2020. Sequences most commonly resembled lineages circulating at earlier times in New York, suggesting transmission primarily from this location, though a minority of Philadelphia genomes matched sequences from other sites, suggesting additional introductions. Multiple genomes showed even closer matches to other Philadelphia isolates, suggestive of ongoing transmission within Philadelphia. We found that all of our isolates contained the D614G substitution in the viral spike and belong to lineages variously designated B.1, Nextstrain clade 20A or 20C, and GISAID clade G or GH. There were no viral sequence polymorphisms detectably associated with disease outcome. For some patients, genome sequences were determined longitudinally or concurrently from multiple body sites. In both cases, some comparisons showed reproducible polymorphisms, suggesting initial seeding with multiple variants and/or accumulation of polymorphisms after infection. These results thus provide data on the sources of SARS-CoV-2 infection in Philadelphia and begin to explore the dynamics within hospitalized patients.IMPORTANCE Understanding how SARS-CoV-2 spreads globally and within infected individuals is critical to the development of mitigation strategies. We found that most lineages in Philadelphia had resembled sequences from New York, suggesting infection primarily but not exclusively from this location. Many genomes had even nearer neighbors within Philadelphia, indicating local spread. Multiple genome sequences were available for some subjects and in a subset of cases could be shown to differ between time points and body sites within an individual, indicating heterogeneous viral populations within individuals and raising questions on the mechanisms responsible. There was no evidence that different lineages were associated with different outcomes in patients, emphasizing the importance of individual-specific vulnerability.

Microbial Lineages in Sarcoidosis. A Metagenomic Analysis Tailored for Low-Microbial Content Samples

RATIONALE

The etiology of sarcoidosis is unknown, but microbial agents are suspected as triggers.

OBJECTIVES

We sought to identify bacterial, fungal, or viral lineages in specimens from patients with sarcoidosis enriched relative to control subjects using metagenomic DNA sequencing. Because DNA from environmental contamination contributes disproportionately to samples with low authentic microbial content, we developed improved methods for filtering environmental contamination.

METHODS

We analyzed specimens from subjects with sarcoidosis (n = 93), control subjects without sarcoidosis (n = 72), and various environmental controls (n = 150). Sarcoidosis specimens consisted of two independent sets of formalin-fixed, paraffin-embedded lymph node biopsies, BAL, Kveim reagent, and fresh granulomatous spleen from a patient with sarcoidosis. All specimens were analyzed by bacterial 16S and fungal internal transcribed spacer ribosomal RNA gene sequencing. In addition, BAL was analyzed by shotgun sequencing of fractions enriched for viral particles, and Kveim and spleen were subjected to whole-genome shotgun sequencing.

MEASUREMENTS AND MAIN RESULTS

In one tissue set, fungi in the Cladosporiaceae family were enriched in sarcoidosis compared with nonsarcoidosis tissues; in the other tissue set, we detected enrichment of several bacterial lineages in sarcoidosis but not Cladosporiaceae. BAL showed limited enrichment of Aspergillus fungi. Several microbial lineages were detected in Kveim and spleen, including Cladosporium. No microbial lineage was enriched in more than one sample type after correction for multiple comparisons.

CONCLUSIONS

Metagenomic sequencing revealed enrichment of microbes in single types of sarcoidosis samples but limited concordance across sample types. Statistical analysis accounting for environmental contamination was essential to avoiding false positives.

Lung-enriched organisms and aberrant bacterial and fungal respiratory microbiota after lung transplant

RATIONALE

Long-term survival after lung transplantation is limited by infectious complications and by bronchiolitis obliterans syndrome (BOS), a form of chronic rejection linked in part to microbial triggers.

OBJECTIVES

To define microbial populations in the respiratory tract of transplant patients comprehensively using unbiased high-density sequencing.

METHODS

Lung was sampled by bronchoalveolar lavage (BAL) and upper respiratory tract by oropharyngeal wash (OW). Bacterial 16S rDNA and fungal internal transcribed spacer sequencing was used to profile organisms present. Outlier analysis plots defining taxa enriched in lung relative to OW were used to identify bacteria enriched in lung against a background of oropharyngeal carryover.

MEASUREMENTS AND MAIN RESULTS

Lung transplant recipients had higher bacterial burden in BAL than control subjects, frequent appearance of dominant organisms, greater distance between communities in BAL and OW indicating more distinct populations, and decreased respiratory tract microbial richness and diversity. Fungal populations were typically dominated by Candida in both sites or by Aspergillus in BAL but not OW. 16S outlier analysis identified lung-enriched taxa indicating bacteria replicating in the lower respiratory tract. In some cases this confirmed respiratory cultures but in others revealed enrichment by anaerobic organisms or mixed outgrowth of upper respiratory flora and provided quantitative data on relative abundances of bacteria found by culture.

CONCLUSIONS

Respiratory tract microbial communities in lung transplant recipients differ in structure and composition from healthy subjects. Outlier analysis can identify specific bacteria replicating in lung. These findings provide novel approaches to address the relationship between microbial communities and transplant outcome and aid in assessing lung infections.

Topographical continuity of bacterial populations in the healthy human respiratory tract

RATIONALE

Defining the biogeography of bacterial populations in human body habitats is a high priority for understanding microbial-host relationships in health and disease. The healthy lung was traditionally considered sterile, but this notion has been challenged by emerging molecular approaches that enable comprehensive examination of microbial communities. However, studies of the lung are challenging due to difficulties in working with low biomass samples.

OBJECTIVES

Our goal was to use molecular methods to define the bacterial microbiota present in the lungs of healthy individuals and assess its relationship to upper airway populations.

METHODS

We sampled respiratory flora intensively at multiple sites in six healthy individuals. The upper tract was sampled by oral wash and oro-/nasopharyngeal swabs. Two bronchoscopes were used to collect samples up to the glottis, followed by serial bronchoalveolar lavage and lower airway protected brush. Bacterial abundance and composition were analyzed by 16S rDNA Q-PCR and deep sequencing.

MEASUREMENTS AND MAIN RESULTS

Bacterial communities from the lung displayed composition indistinguishable from the upper airways, but were 2 to 4 logs lower in biomass. Lung-specific sequences were rare and not shared among individuals. There was no unique lung microbiome.

CONCLUSIONS

In contrast to other organ systems, the respiratory tract harbors a homogenous microbiota that decreases in biomass from upper to lower tract. The healthy lung does not contain a consistent distinct microbiome, but instead contains low levels of bacterial sequences largely indistinguishable from upper respiratory flora. These findings establish baseline data for healthy subjects and sampling approaches for sequence-based analysis of diseases.