IgA deficiency destabilizes homeostasis toward intestinal microbes and increases systemic immune dysregulation.

The ability of most patients with selective immunoglobulin A (IgA) deficiency (SIgAD) to remain apparently healthy has been a persistent clinical conundrum. Compensatory mechanisms, including IgM, have been proposed, yet it remains unclear how secretory IgA and IgM work together in the mucosal system and, on a larger scale, whether the systemic and mucosal anti-commensal responses are redundant or have unique features. To address this gap in knowledge, we developed an integrated host-commensal approach combining microbial flow cytometry and metagenomic sequencing (mFLOW-Seq) to comprehensively define which microbes induce mucosal and systemic antibodies. We coupled this approach with high-dimensional immune profiling to study a cohort of pediatric patients with SIgAD and household control siblings. We found that mucosal and systemic antibody networks cooperate to maintain homeostasis by targeting a common subset of commensal microbes. In IgA-deficiency, we find increased translocation of specific bacterial taxa associated with elevated levels of systemic IgG targeting fecal microbiota. Associated features of immune system dysregulation in IgA-deficient mice and humans included elevated levels of inflammatory cytokines, enhanced follicular CD4 T helper cell frequency and activation, and an altered CD8 T cell activation state. Although SIgAD is clinically defined by the absence of serum IgA, the symptomatology and immune dysregulation were concentrated in the SIgAD participants who were also fecal IgA deficient. These findings reveal that mucosal IgA deficiency leads to aberrant systemic exposures and immune responses to commensal microbes, which increase the likelihood of humoral and cellular immune dysregulation and symptomatic disease in patients with IgA deficiency.

Arresting microbiome development limits immune system maturation and resistance to infection in mice.

Disruptions to the intestinal microbiome during weaning lead to negative effects on host immune function. However, the critical host-microbe interactions during weaning that are required for immune system development remain poorly understood. We find that restricting microbiome maturation during weaning stunts immune system development and increases susceptibility to enteric infection. We developed a gnotobiotic mouse model of the early-life microbiome Pediatric Community (PedsCom). These mice develop fewer peripheral regulatory T cells and less IgA, hallmarks of microbiota-driven immune system development. Furthermore, adult PedsCom mice retain high susceptibility to Salmonella infection, which is characteristic of young mice and children. Altogether, our work illustrates how the post-weaning transition in microbiome composition contributes to normal immune maturation and protection from infection. Accurate modeling of the pre-weaning microbiome provides a window into the microbial requirements for healthy development and suggests an opportunity to design microbial interventions at weaning to improve immune development in human infants.

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.

Comparative Analysis of Emerging B.1.1.7+E484K SARS-CoV-2 Isolates.

We report the genome of a B.1.1.7+E484K severe acute respiratory syndrome coronavirus 2 from Southeastern Pennsylvania and compare it with all high-coverage B.1.1.7+E484K genomes (n = 235) available. Analyses showed the existence of at least 4 distinct clades of this variant circulating in the United States and the possibility of at least 59 independent acquisitions of the E484K mutation.

Comparative Analysis of Emerging B.1.1.7+E484K SARS-CoV-2 isolates from Pennsylvania.

Rapid whole genome sequencing of SARS-CoV-2 has presented the ability to detect new emerging variants of concern in near real time. Here we report the genome of a virus isolated in Pennsylvania in March 2021 that was identified as lineage B.1.1.7 (VOC-202012/01) that also harbors the E484K spike mutation, which has been shown to promote "escape" from neutralizing antibodies in vitro. We compare this sequence to the only 5 other B.1.1.7+E484K genomes from Pennsylvania, all of which were isolated in mid March. Beginning in February 2021, only a small number (n=60) of isolates with this profile have been detected in the US, and only a total of 253 have been reported globally (first in the UK in December 2020). Comparative genomics of all currently available high coverage B.1.1.7+E484K genomes (n=235) available on GISAID suggested the existence of 7 distinct groups or clonal complexes (CC; as defined by GNUVID) bearing the E484K mutation raising the possibility of 7 independent acquisitions of the E484K spike mutation in each background. Phylogenetic analysis suggested the presence of at least 3 distinct clades of B.1.1.7+E484K circulating in the US, with the Pennsylvanian isolates belonging to two distinct clades. Increased genomic surveillance will be crucial for detection of emerging variants of concern that can escape natural and vaccine induced immunity.

Adaptation by Ancient Horizontal Acquisition of Butyrate Metabolism Genes in Aggregatibacter actinomycetemcomitans.

Like the bacterial residents of the human gut, it is likely that many of the species in the human oral microbiota have evolved to better occupy and persist in their niche. Aggregatibacter actinomycetemcomitans (Aa) is both a common colonizer of the oral cavity and has been implicated in the pathogenesis of periodontal disease. Here, we present a whole-genome phylogenetic analysis of Aa isolates from humans and nonhuman primates that revealed an ancient origin for this species and a long history of association with the Catarrhini, the lineage that includes Old World monkeys (OWM) and humans. Further genomic analysis showed a strong association with the presence of a short-chain fatty acid (SCFA) catabolism locus (atoRDAEB) in many human isolates that was absent in almost all nonhuman OWM isolates. We show that this locus was likely acquired through horizontal gene transfer. When grown under conditions that are similar to those at the subgingival site of periodontitis (anaerobic, SCFA replete), Aa strains with atoRDAEB formed robust biofilms and showed upregulation of genes involved in virulence, colonization, and immune evasion. Both an isogenic deletion mutant and nonhuman primate isolates lacking the ato locus failed to grow in a robust biofilm under these conditions, but grew well under the carbohydrate-rich conditions similar to those found above the gumline. We propose that the acquisition of the ato locus was a key evolutionary step allowing Aa to utilize SCFAs, adapt, and modulate subgingival disease.IMPORTANCE There has been considerable interest in the impact of short-chain fatty acids (SCFAs) on inflammatory effects related to the microbiome. Here, we present evidence that SCFAs may also be important in disease by providing an energy source or disease-associated cue for colonizing pathogens. We propose that SCFAs allow Aggregatibacter actinomycetemcomitans (Aa) to adapt to the subgingival anaerobic environment, which is the site of human periodontitis. Under anaerobic, SCFA-rich conditions, human-derived Aa strains that possess butyrate metabolism genes form strong biofilms and upregulate virulence genes. Our phylogenetic analysis highlights a long history of evolution of Aa with its primate hosts and suggests that the acquisition of butyrate metabolism genes may have been a critical step in allowing Aa to colonize a new niche and cause disease in humans. Overall, this study highlights the important role that horizontal gene transfer may play in microbial adaptation and the evolution of infectious disease.

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.

Diet-induced microbial autofluorescence confounds flow cytometry of ex vivo isolated fecal microbes.

Microbial flow cytometry is a powerful emerging technology with a broad range of applications including the study of complex microbial communities. Immunologists are increasingly using this technology to study antibody responses against pathogenic and commensal microbes. We employed microbial flow cytometry to quantify the proportion of fecal microbes bound by six different Ig isotypes: IgA, IgM, IgG1, IgG2b, IgG2c, and IgG3. In healthy mammals, secretory IgA (sIgA) binds to a subset of commensal microbes in the gut whereas IgG is not typically found in the intestinal tract of healthy mammals. Unexpectedly, fecal microbes isolated from SPF C57BL/6 mice housed in the Hill facility and imported from the vendors The Jackson Laboratory and Taconic Biosciences showed a strong signal in the Brilliant Violet 711 (BV711) channel. Unstained fecal samples from these mice demonstrated that the BV711 signal was due to bacterial autofluorescence. We found that murine diets containing alfalfa induce ex vivo microbial autofluorescence in the far red spectrum, likely due to chlorophyll. Analysis of unstained intestinal microbes is an important step in microbial flow cytometry to identify diet-induced autofluorescence. We recommend fluorophores with emission spectra below 650 nm (e.g. BV421, PE).