Sequence Analysis of 20,453 Severe Acute Respiratory Syndrome Coronavirus 2 Genomes from the Houston Metropolitan Area Identifies the Emergence and Widespread Distribution of Multiple Isolates of All Major Variants of Concern.

Since the beginning of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, there has been international concern about the emergence of virus variants with mutations that increase transmissibility, enhance escape from the human immune response, or otherwise alter biologically important phenotypes. In late 2020, several variants of concern emerged globally, including the UK variant (B.1.1.7), the South Africa variant (B.1.351), Brazil variants (P.1 and P.2), and two related California variants of interest (B.1.429 and B.1.427). These variants are believed to have enhanced transmissibility. For the South Africa and Brazil variants, there is evidence that mutations in spike protein permit it to escape from some vaccines and therapeutic monoclonal antibodies. On the basis of our extensive genome sequencing program involving 20,453 coronavirus disease 2019 patient samples collected from March 2020 to February 2021, we report identification of all six of these SARS-CoV-2 variants among Houston Methodist Hospital (Houston, TX) patients residing in the greater metropolitan area. Although these variants are currently at relatively low frequency (aggregate of 1.1%) in the population, they are geographically widespread. Houston is the first city in the United States in which active circulation of all six current variants of concern has been documented by genome sequencing. As vaccine deployment accelerates, increased genomic surveillance of SARS-CoV-2 is essential to understanding the presence, frequency, and medical impact of consequential variants and their patterns and trajectory of dissemination.

Molecular Architecture of Early Dissemination and Massive Second Wave of the SARS-CoV-2 Virus in a Major Metropolitan Area.

We sequenced the genomes of 5,085 SARS-CoV-2 strains causing two COVID-19 disease waves in metropolitan Houston, Texas, an ethnically diverse region with seven million residents. The genomes were from viruses recovered in the earliest recognized phase of the pandemic in Houston, and an ongoing massive second wave of infections. The virus was originally introduced into Houston many times independently. Virtually all strains in the second wave have a Gly614 amino acid replacement in the spike protein, a polymorphism that has been linked to increased transmission and infectivity. Patients infected with the Gly614 variant strains had significantly higher virus loads in the nasopharynx on initial diagnosis. We found little evidence of a significant relationship between virus genotypes and altered virulence, stressing the linkage between disease severity, underlying medical conditions, and host genetics. Some regions of the spike protein - the primary target of global vaccine efforts - are replete with amino acid replacements, perhaps indicating the action of selection. We exploited the genomic data to generate defined single amino acid replacements in the receptor binding domain of spike protein that, importantly, produced decreased recognition by the neutralizing monoclonal antibody CR30022. Our study is the first analysis of the molecular architecture of SARS-CoV-2 in two infection waves in a major metropolitan region. The findings will help us to understand the origin, composition, and trajectory of future infection waves, and the potential effect of the host immune response and therapeutic maneuvers on SARS-CoV-2 evolution.

Molecular Architecture of Early Dissemination and Massive Second Wave of the SARS-CoV-2 Virus in a Major Metropolitan Area.

We sequenced the genomes of 5,085 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains causing two coronavirus disease 2019 (COVID-19) disease waves in metropolitan Houston, TX, an ethnically diverse region with 7 million residents. The genomes were from viruses recovered in the earliest recognized phase of the pandemic in Houston and from viruses recovered in an ongoing massive second wave of infections. The virus was originally introduced into Houston many times independently. Virtually all strains in the second wave have a Gly614 amino acid replacement in the spike protein, a polymorphism that has been linked to increased transmission and infectivity. Patients infected with the Gly614 variant strains had significantly higher virus loads in the nasopharynx on initial diagnosis. We found little evidence of a significant relationship between virus genotype and altered virulence, stressing the linkage between disease severity, underlying medical conditions, and host genetics. Some regions of the spike protein-the primary target of global vaccine efforts-are replete with amino acid replacements, perhaps indicating the action of selection. We exploited the genomic data to generate defined single amino acid replacements in the receptor binding domain of spike protein that, importantly, produced decreased recognition by the neutralizing monoclonal antibody CR3022. Our report represents the first analysis of the molecular architecture of SARS-CoV-2 in two infection waves in a major metropolitan region. The findings will help us to understand the origin, composition, and trajectory of future infection waves and the potential effect of the host immune response and therapeutic maneuvers on SARS-CoV-2 evolution.IMPORTANCE There is concern about second and subsequent waves of COVID-19 caused by the SARS-CoV-2 coronavirus occurring in communities globally that had an initial disease wave. Metropolitan Houston, TX, with a population of 7 million, is experiencing a massive second disease wave that began in late May 2020. To understand SARS-CoV-2 molecular population genomic architecture and evolution and the relationship between virus genotypes and patient features, we sequenced the genomes of 5,085 SARS-CoV-2 strains from these two waves. Our report provides the first molecular characterization of SARS-CoV-2 strains causing two distinct COVID-19 disease waves.

Carbapenems drive the collateral resistance to ceftaroline in cystic fibrosis patients with MRSA.

Chronic airways infection with methicillin-resistant Staphylococcus aureus (MRSA) is associated with worse respiratory disease cystic fibrosis (CF) patients. Ceftaroline is a cephalosporin that inhibits the penicillin-binding protein (PBP2a) uniquely produced by MRSA. We analyzed 335 S. aureus isolates from CF sputum samples collected at three US centers between 2015-2018. Molecular relationships demonstrated that high-level resistance of preceding isolates to carbapenems were associated with subsequent isolation of ceftaroline resistant CF MRSA. In vitro evolution experiments showed that pre-exposure of CF MRSA to meropenem with further selection with ceftaroline implied mutations in mecA and additional mutations in pbp1 and pbp2, targets of carbapenems; no effects were achieved by other ß-lactams. An in vivo pneumonia mouse model showed the potential therapeutic efficacy of ceftaroline/meropenem combination against ceftaroline-resistant CF MRSA infections. Thus, the present findings highlight risk factors and potential therapeutic strategies offering an opportunity to both prevent and address antibiotic resistance in this patient population.

Treatment of COVID-19 Patients with Convalescent Plasma in Houston, Texas.

BACKGROUND: COVID-19 disease, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spread globally, and no proven treatments are available. Convalescent plasma therapy has been used with varying degrees of success to treat severe microbial infections for more than 100 years. METHODS: Patients (n=25) with severe and/or life-threatening COVID-19 disease were enrolled at the Houston Methodist hospitals from March 28 to April 14, 2020. Patients were transfused with convalescent plasma obtained from donors with confirmed SARS-CoV-2 infection and had been symptom free for 14 days. The primary study outcome was safety, and the secondary outcome was clinical status at day 14 post-transfusion. Clinical improvement was assessed based on a modified World Health Organization 6-point ordinal scale and laboratory parameters. Viral genome sequencing was performed on donor and recipient strains. RESULTS: At baseline, all patients were receiving supportive care, including anti-inflammatory and anti-viral treatments, and all patients were on oxygen support. At day 7 post-transfusion with convalescent plasma, nine patients had at least a 1-point improvement in clinical scale, and seven of those were discharged. By day 14 post-transfusion, 19 (76%) patients had at least a 1-point improvement in clinical status and 11 were discharged. No adverse events as a result of plasma transfusion were observed. The whole genome sequencing data did not identify a strain genotype-disease severity correlation. CONCLUSIONS: The data indicate that administration of convalescent plasma is a safe treatment option for those with severe COVID-19 disease. Randomized, controlled trials are needed to determine its efficacy.

Population Genomic Molecular Epidemiological Study of Macrolide-Resistant Streptococcus pyogenes in Iceland, 1995 to 2016: Identification of a Large Clonal Population with a pbp2x Mutation Conferring Reduced In Vitro ß-Lactam Susceptibility.

Resistance to macrolide antibiotics is a global concern in the treatment of Streptococcus pyogenes (group A Streptococcus [GAS]) infections. In Iceland, since the detection of the first macrolide-resistant isolate in 1998, three epidemic waves of macrolide-resistant GAS infections have occurred, with peaks in 1999, 2004, and 2008. We conducted whole-genome sequencing of all 1,575 available GAS macrolide-resistant clinical isolates of all infection types collected at the national reference laboratory in Reykjavik, Iceland, from 1998 to 2016. Among 1,515 erythromycin-resistant isolates, 90.3% were of only three emm types, emm4 (n = 713), emm6 (n = 324), and emm12 (n = 332), with each being predominant in a distinct epidemic peak. The antibiotic efflux pump genes, mef(A) and msr(D), were present on chimeric mobile genetic elements in 99.3% of the macrolide-resistant isolates of these emm types. Of note, in addition to macrolide resistance, virtually all emm12 isolates had a single amino acid substitution in penicillin-binding protein PBP2X that conferred a 2-fold increased penicillin G and ampicillin MIC among the isolates tested. We conclude that each of the three large epidemic peaks of macrolide-resistant GAS infections occurring in Iceland since 1998 was caused by the emergence and clonal expansion of progenitor strains, with macrolide resistance being conferred predominantly by inducible Mef(A) and Msr(D) drug efflux pumps. The occurrence of emm12 strains with macrolide resistance and decreased beta-lactam susceptibility was unexpected and is of public health concern.

Streptococcus pyogenes genes that promote pharyngitis in primates.

Streptococcus pyogenes (group A streptococcus; GAS) causes 600 million cases of pharyngitis annually worldwide. There is no licensed human GAS vaccine despite a century of research. Although the human oropharynx is the primary site of GAS infection, the pathogenic genes and molecular processes used to colonize, cause disease, and persist in the upper respiratory tract are poorly understood. Using dense transposon mutant libraries made with serotype M1 and M28 GAS strains and transposon-directed insertion sequencing, we performed genome-wide screens in the nonhuman primate (NHP) oropharynx. We identified many potentially novel GAS fitness genes, including a common set of 115 genes that contribute to fitness in both genetically distinct GAS strains during experimental NHP pharyngitis. Targeted deletion of 4 identified fitness genes/operons confirmed that our newly identified targets are critical for GAS virulence during experimental pharyngitis. Our screens discovered many surface-exposed or secreted proteins - substrates for vaccine research - that potentially contribute to GAS pharyngitis, including lipoprotein HitA. Pooled human immune globulin reacted with purified HitA, suggesting that humans produce antibodies against this lipoprotein. Our findings provide new information about GAS fitness in the upper respiratory tract that may assist in translational research, including developing novel vaccines.

Genetic heterogeneity of the Spy1336/R28-Spy1337 virulence axis in Streptococcus pyogenes and effect on gene transcript levels and pathogenesis.

Streptococcus pyogenes is a strict human pathogen responsible for more than 700 million infections annually worldwide. Strains of serotype M28 S. pyogenes are typically among the five more abundant types causing invasive infections and pharyngitis in adults and children. Type M28 strains also have an unusual propensity to cause puerperal sepsis and neonatal disease. We recently discovered that a one-nucleotide indel in an intergenic homopolymeric tract located between genes Spy1336/R28 and Spy1337 altered virulence in a mouse model of infection. In the present study, we analyzed size variation in this homopolymeric tract and determined the extent of heterogeneity in the number of tandemly-repeated 79-amino acid domains in the coding region of Spy1336/R28 in large samples of strains recovered from humans with invasive infections. Both repeat sequence elements are highly polymorphic in natural populations of M28 strains. Variation in the homopolymeric tract results in (i) changes in transcript levels of Spy1336/R28 and Spy1337 in vitro, (ii) differences in virulence in a mouse model of necrotizing myositis, and (iii) global transcriptome changes as shown by RNAseq analysis of isogenic mutant strains. Variation in the number of tandem repeats in the coding sequence of Spy1336/R28 is responsible for size variation of R28 protein in natural populations. Isogenic mutant strains in which genes encoding R28 or transcriptional regulator Spy1337 are inactivated are significantly less virulent in a nonhuman primate model of necrotizing myositis. Our findings provide impetus for additional studies addressing the role of R28 and Spy1337 variation in pathogen-host interactions.

Environmental pH and peptide signaling control virulence of Streptococcus pyogenes via a quorum-sensing pathway.

Bacteria control gene expression in concert with their population density by a process called quorum sensing, which is modulated by bacterial chemical signals and environmental factors. In the human pathogen Streptococcus pyogenes, production of secreted virulence factor SpeB is controlled by a quorum-sensing pathway and environmental pH. The quorum-sensing pathway consists of a secreted leaderless peptide signal (SIP), and its cognate receptor RopB. Here, we report that the SIP quorum-sensing pathway has a pH-sensing mechanism operative through a pH-sensitive histidine switch located at the base of the SIP-binding pocket of RopB. Environmental acidification induces protonation of His144 and reorganization of hydrogen bonding networks in RopB, which facilitates SIP recognition. The convergence of two disparate signals in the SIP signaling pathway results in induction of SpeB production and increased bacterial virulence. Our findings provide a model for investigating analogous crosstalk in other microorganisms.

Integrated analysis of population genomics, transcriptomics and virulence provides novel insights into Streptococcus pyogenes pathogenesis.

Streptococcus pyogenes causes 700 million human infections annually worldwide, yet, despite a century of intensive effort, there is no licensed vaccine against this bacterium. Although a number of large-scale genomic studies of bacterial pathogens have been published, the relationships among the genome, transcriptome, and virulence in large bacterial populations remain poorly understood. We sequenced the genomes of 2,101 emm28 S. pyogenes invasive strains, from which we selected 492 phylogenetically diverse strains for transcriptome analysis and 50 strains for virulence assessment. Data integration provided a novel understanding of the virulence mechanisms of this model organism. Genome-wide association study, expression quantitative trait loci analysis, machine learning, and isogenic mutant strains identified and confirmed a one-nucleotide indel in an intergenic region that significantly alters global transcript profiles and ultimately virulence. The integrative strategy that we used is generally applicable to any microbe and may lead to new therapeutics for many human pathogens.

Gene fitness landscape of group A streptococcus during necrotizing myositis.

Necrotizing fasciitis and myositis are devastating infections characterized by high mortality. Group A streptococcus (GAS) is a common cause of these infections, but the molecular pathogenesis is poorly understood. We report a genome-wide analysis using serotype M1 and M28 strains that identified GAS genes contributing to necrotizing myositis in nonhuman primates (NHP), a clinically relevant model. Using transposon-directed insertion-site sequencing (TraDIS), we identified 126 and 116 GAS genes required for infection by serotype M1 and M28 organisms, respectively. For both M1 and M28 strains, more than 25% of the GAS genes required for necrotizing myositis encode known or putative transporters. Thirteen GAS transporters contributed to both M1 and M28 strain fitness in NHP myositis, including putative importers for amino acids, carbohydrates, and vitamins and exporters for toxins, quorum-sensing peptides, and uncharacterized molecules. Targeted deletion of genes encoding 5 transporters confirmed that each isogenic mutant strain was significantly (P < 0.05) impaired in causing necrotizing myositis in NHPs. Quantitative reverse-transcriptase PCR (qRT-PCR) analysis showed that these 5 genes are expressed in infected NHP and human skeletal muscle. Certain substrate-binding lipoproteins of these transporters, such as Spy0271 and Spy1728, were previously documented to be surface exposed, suggesting that our findings have translational research implications.