ABSTRACT Detection of multiple pneumococcal serotype carriage can enhance monitoring of pneumococcal vaccine impact, particularly among high-burden childhood populations. We assessed methods for identifying co-carriage of pneumococcal serotypes from whole-genome sequences. Twenty-four nasopharyngeal samples were collected during community carriage surveillance from healthy children in Blantyre, Malawi, which were then serotyped by microarray. Pneumococcal DNA from culture plate sweeps were sequenced using Illumina MiSeq, and genomic serotyping was carried out using SeroCall and PneumoKITy. Their sensitivity was calculated in reference to the microarray data. Local maxima in the single-nucleotide polymorphism (SNP) density distributions were assessed for their correspondence to the relative abundance of serotypes. Across the 24 individuals, the microarray detected 77 non-unique serotypes, of which 42 occurred at high relative abundance (>10%) (per individual, median, 3; range, 1–6 serotypes). The average sequencing depth was 57X (range: 21X–88X). The sensitivity of SeroCall for identifying high-abundance serotypes was 98% (95% CI, 0.87–1.00), 20% (0.08–0.36) for low abundance (<10%), and 62% (0.50–0.72) overall. PneumoKITy’s sensitivity was 86% (0.72–0.95), 20% (0.06–0.32), and 56% (0.42–0.65), respectively. Local maxima in the SNP frequency distribution were highly correlated with the relative abundance of high-abundance serotypes. Six samples were resequenced, and the pooled runs had an average fourfold increase in sequencing depth. This allowed genomic serotyping of two of the previously undetectable seven low-abundance serotypes. Genomic serotyping is highly sensitive for the detection of high-abundance serotypes in samples with co-carriage. Serotype-associated reads may be identified through SNP frequency, and increased read depth can increase sensitivity for low-abundance serotype detection. IMPORTANCE Pneumococcal carriage is a prerequisite for invasive pneumococcal disease, which is a leading cause of childhood pneumonia. Multiple carriage of unique pneumococcal serotypes at a single time point is prevalent among high-burden childhood populations. This study assessed the sensitivity of different genomic serotyping methods for identifying pneumococcal serotypes during co-carriage. These methods were evaluated against the current gold standard for co-carriage detection. The results showed that genomic serotyping methods have high sensitivity for detecting high-abundance serotypes in samples with co-carriage, and increasing sequencing depth can increase sensitivity for low-abundance serotypes. These results are important for monitoring vaccine impact, which aims to reduce the prevalence of specific pneumococcal serotypes. By accurately detecting and identifying multiple pneumococcal serotypes in carrier populations, we can better evaluate the effectiveness of vaccination programs.

The principal barrier to an HIV cure is the presence of a latent viral reservoir (LVR) made up primarily of latently infected resting CD4+ (rCD4) T-cells. Studies in the United States have shown that the LVR decays slowly (half-life=3.8 years), but this rate in African populations has been understudied. This study examined longitudinal changes in the inducible replication competent LVR (RC-LVR) of ART-suppressed Ugandans living with HIV (n=88) from 2015-2020 using the quantitative viral outgrowth assay, which measures infectious units per million (IUPM) rCD4 T-cells. In addition, outgrowth viruses were examined with site-directed next-generation sequencing to assess for possible ongoing viral evolution. During the study period (2018-19), Uganda instituted a nationwide rollout of first-line ART consisting of Dolutegravir (DTG) with two NRTI, which replaced the previous regimen that consisted of one NNRTI and the same two NRTI. Changes in the RC-LVR were analyzed using two versions of a novel Bayesian model that estimated the decay rate over time on ART as a single, linear rate (model A) or allowing for an inflection at time of DTG initiation (model B). Model A estimated the population-level slope of RC-LVR change as a non-significant positive increase. This positive slope was due to a temporary increase in the RC-LVR that occurred 0-12 months post-DTG initiation (p<0.0001). This was confirmed with model B, which estimated a significant decay pre-DTG initiation with a half-life of 7.7 years, but a significant positive slope post-DTG initiation leading to a transient estimated doubling-time of 8.1 years. There was no evidence of viral failure in the cohort, or consistent evolution in the outgrowth sequences associated with DTG initiation. These data suggest that either the initiation of DTG, or cessation of NNRTI use, is associated with a significant temporary increase in the circulating RC-LVR.

Background: The ability of SARS-CoV-2 vaccines to protect against infection and onward transmission determines whether immunisation can control global circulation. We estimated the effectiveness of Pfizer-BioNTech mRNA vaccine (BNT162b2) and Oxford AstraZeneca adenovirus vector vaccine (ChAdOx1) vaccines against acquisition and transmission of the Alpha and Delta variants in a prospective household study in England. Methods: Households were recruited based on adult purported index cases testing positive after reverse transcription-quantitative (RT-q)PCR testing of oral-nasal swabs. Purported index cases and their household contacts took oral-nasal swabs on days 1, 3 and 7 after enrolment and a subset of the PCR-positive swabs underwent genomic sequencing conducted on a subset. We used Bayesian logistic regression to infer vaccine effectiveness against acquisition and transmission, adjusted for age, vaccination history and variant. Results: Between 2 February 2021 and 10 September 2021, 213 index cases and 312 contacts were followed up. After excluding households lacking genomic proximity (N=2) or with unlikely serial intervals (N=16), 195 households with 278 contacts remained, of whom 113 (41%) became PCR positive. Delta lineages had 1.53 times the risk (95% Credible Interval: 1.04 – 2.20) of transmission than Alpha; contacts older than 18 years old were 1.48 (1.20 – 1.91) and 1.02 (0.93 – 1.16) times more likely to acquire an Alpha or Delta infection than children. Effectiveness of two doses of BNT162b2 against transmission of Delta was 36% (-1%, 66%) and 49% (18%, 73%) for ChAdOx1, similar to their effectiveness for Alpha. Protection against infection with Alpha was higher than for Delta, 69% (9%, 95%) vs. 18% (-11%, 59%), respectively, for BNT162b2 and 24% (-41%, 72%) vs. 9% (-15%, 42%), respectively, for ChAdOx1. Conclusions: BNT162b2 and ChAdOx1 reduce transmission of the Delta variant from breakthrough infections in the household setting, although their protection against infection within this setting is low.

Abstract The timing of the establishment of the HIV latent viral reservoir (LVR) is of particular interest, as there is evidence that proviruses are preferentially archived at the time of antiretroviral therapy (ART) initiation. Quantitative viral outgrowth assays (QVOAs) were performed using Peripheral Blood Mononuclear Cells (PBMC) collected from Ugandans living with HIV who were virally suppressed on ART for >1 year, had known seroconversion windows, and at least two archived ART-naïve plasma samples. QVOA outgrowth populations and pre-ART plasma samples were deep sequenced for the pol and gp41 genes. The bayroot program was used to estimate the date that each outgrowth virus was incorporated into the reservoir. Bayroot was also applied to previously published data from a South African cohort. In the Ugandan cohort (n = 11), 87.9 per cent pre-ART and 56.3 per cent viral outgrowth sequences were unique. Integration dates were estimated to be relatively evenly distributed throughout viremia in 9/11 participants. In contrast, sequences from the South African cohort (n = 9) were more commonly estimated to have entered the LVR close to ART initiation, as previously reported. Timing of LVR establishment is variable between populations and potentially viral subtypes, which could limit the effectiveness of interventions that target the LVR only at ART initiation.

Background The ability of SARS-CoV-2 vaccines to protect against infection and onward transmission determines whether immunisation can control global circulation. We estimated the effectiveness of Pfizer-BioNTech mRNA vaccine (BNT162b2) and Oxford AstraZeneca adenovirus vector vaccine (ChAdOx1) vaccines against acquisition and transmission of the Alpha and Delta variants in a prospective household study in England. Methods Households were recruited based on adult purported index cases testing positive after reverse transcription-quantitative (RT-q)PCR testing of oral-nasal swabs. Purported index cases and their household contacts took oral-nasal swabs on days 1, 3 and 7 after enrolment and a subset of the PCR-positive swabs underwent genomic sequencing conducted on a subset. We used Bayesian logistic regression to infer vaccine effectiveness against acquisition and transmission, adjusted for age, vaccination history and variant. Results Between 2 February 2021 and 10 September 2021, 213 index cases and 312 contacts were followed up. After excluding households lacking genomic proximity (N=2) or with unlikely serial intervals (N=16), 195 households with 278 contacts remained, of whom 113 (41%) became PCR positive. Delta lineages had 1.53 times the risk (95% Credible Interval: 1.04 – 2.20) of transmission than Alpha; contacts older than 18 years old were 1.48 (1.20 – 1.91) and 1.02 (0.93 – 1.16) times more likely to acquire an Alpha or Delta infection than children. Effectiveness of two doses of BNT162b2 against transmission of Delta was 36% (-1%, 66%) and 49% (18%, 73%) for ChAdOx1, similar to their effectiveness for Alpha. Protection against infection with Alpha was higher than for Delta, 69% (9%, 95%) vs. 18% (-11%, 59%), respectively, for BNT162b2 and 24% (-41%, 72%) vs. 9% (-15%, 42%), respectively, for ChAdOx1. Conclusions BNT162b2 and ChAdOx1 reduce transmission of the Delta variant from breakthrough infections in the household setting, although their protection against infection within this setting is low.

Background: Inference on pneumococcal transmission has mostly relied on longitudinal studies which are costly and resource intensive. Therefore, we conducted a pilot study to test the ability to infer who infected whom from cross-sectional pneumococcal sequences using phylogenetic inference. Methods: Five suspected transmission pairs, for which there was epidemiological evidence of who infected whom, were selected from a household study. For each pair, Streptococcus pneumoniae full genomes were sequenced from nasopharyngeal swabs collected on the same day. The within-host genetic diversity of the pneumococcal population was used to infer the transmission direction and then cross-validated with the direction suggested by the epidemiological records. Results: The pneumococcal genomes clustered into the five households from which the samples were taken. The proportion of concordantly inferred transmission direction generally increased with increasing minimum genome fragment size and single nucleotide polymorphisms. We observed a larger proportion of unique polymorphic sites in the source bacterial population compared to that of the recipient in four of the five pairs, as expected in the case of a transmission bottleneck. The only pair that did not exhibit this effect was also the pair that had consistent discordant transmission direction compared to the epidemiological records suggesting potential misdirection as a result of false-negative sampling. Conclusions: This pilot provided support for further studies to test if the direction of pneumococcal transmission can be reliably inferred from cross-sectional samples if sequenced with sufficient depth and fragment length.

Background: Inference on pneumococcal transmission has mostly relied on longitudinal studies which are costly and resource intensive. Therefore, we conducted a pilot study to test the ability to infer who infected whom from cross-sectional pneumococcal sequences using phylogenetic inference. Methods: Five suspected transmission pairs, for which there was epidemiological evidence of who infected whom, were selected from a household study. For each pair, Streptococcus pneumoniae full genomes were sequenced from nasopharyngeal swabs collected on the same day. The within-host genetic diversity of the pneumococcal population was used to infer the transmission direction and then cross-validated with the direction suggested by the epidemiological records. Results: The pneumococcal genomes clustered into the five households from which the samples were taken. The proportion of concordantly inferred transmission direction generally increased with increasing minimum genome fragment size and single nucleotide polymorphisms. We observed a larger proportion of unique polymorphic sites in the source bacterial population compared to that of the recipient in four of the five pairs, as expected in the case of a transmission bottleneck. The only pair that did not exhibit this effect was also the pair that had consistent discordant transmission direction compared to the epidemiological records suggesting potential misdirection as a result of false-negative sampling. Conclusions: This pilot provided support for further studies to test if the direction of pneumococcal transmission can be reliably inferred from cross-sectional samples if sequenced with sufficient depth and fragment length.

BACKGROUND:The human immunodeficiency virus (HIV)-1 latent reservoir (LR) in resting CD4+ T cells is a barrier to cure. LR measurements are commonly performed on blood samples and therefore may miss latently infected cells residing in tissues, including lymph nodes.METHODS:We determined the frequency of intact HIV-1 proviruses and proviral inducibility in matched peripheral blood (PB) and lymph node (LN) samples from 10 HIV-1-infected patients on antiretroviral therapy (ART) using the intact proviral DNA assay and a novel quantitative viral induction assay. Prominent viral sequences from induced viral RNA were characterized using a next-generation sequencing assay.RESULTS:The frequencies of CD4+ T cells with intact proviruses were not significantly different in PB versus LN (61/106 vs 104/106 CD4+ cells), and they were substantially lower than frequencies of CD4+ T cells with defective proviruses. The frequencies of CD4+ T cells induced to produce high levels of viral RNA were not significantly different in PB versus LN (4.3/106 vs 7.9/106), but they were 14-fold lower than the frequencies of cells with intact proviruses. Sequencing of HIV-1 RNA from induced proviruses revealed comparable sequences in paired PB and LN samples.CONCLUSIONS:These results further support the use of PB as an appropriate proxy for the HIV-1 LR in secondary lymphoid organs.

Background The availability of effective, oral direct acting antivirals (DAAs) for hepatitis C virus (HCV) treatment has put elimination of HCV as a public health challenge within reach. However, little is known about the characteristics of transmission networks of people who inject drugs (PWID). Methods Sequencing of a segment of the HCV genome was performed on samples collected from a community-based cohort of PWID between August 2005 and December 2016. Phylogenetic trees were inferred, and clusters were identified (70% bootstrap threshold; 0.04 maximum genetic distance threshold). We describe sex, race, age difference, and HIV infection status of potential transmission partners. Logistic regression was used to assess factors associated with being in an HCV cluster. Results Of 508 HCV genotype 1 viremic PWID, 8% ( n = 41) were grouped into 20 clusters, consisting of 19 pairs and 1 triad. In adjusted analyses, female sex (odds ratio [OR] 2.3 [95% confidence interval (CI) 1.2–4.5]) and HIV infection (OR 5.7 [CI 2.7–11.9]) remained independently associated with being in an HCV infection cluster. Conclusions Molecular epidemiological analysis reveals that, in this cohort of PWID in Baltimore, HIV infection and female sex were associated with HCV clustering. Combination HCV prevention interventions targeting HIV infected PWID and addressing HCV infection prevention needs of women have potential to advance HCV elimination efforts.

Abstract The availability of effective, oral direct acting antivirals for hepatitis C virus (HCV) treatment has fueled optimism for HCV elimination through treatment as prevention (TasP) among people who inject drugs (PWID). Identifying characteristics of individuals in transmission networks would provide critical information for the development and implementation of effective, targeted HCV TasP strategies. The AIDS linked to the IntraVenous Experience (ALIVE) cohort has followed PWID in Baltimore since 1988. Sequencing of the HCV core/E1 region (342 nucleotides) was performed on HCV viremic samples from the most recent study visit attended by ALIVE participants between August, 2005 and December, 2016. Outgroup sequences were retrieved from GenBank through a BLAST search for HCV sequences similar to study sequences to support identification of ‘local clusters’ and were aligned to study sequences using Clustal O. Phylogenetic trees were inferred for each of HCV subtype 1a and 1b separately through maximum likelihood analysis implemented in the MEGA X software using the Tamura-Nei model with gamma distribution and invariant sites. Nucleotide substitution model selection was based on the corrected Akaike information criterion scores of various models in MEGA. Robustness of the resulting tree was assessed by bootstrapping with 1,000 replicates. Clusters were identified using ClusterPicker software (70% bootstrap threshold and 0.05 maximum genetic distance threshold). Sensitivity analyses were performed by varying the genetic distance threshold between 0.025 and 0.05 to determine the effect on identification of factors associated with clustering. HCV infection clustering was defined as > 2 participants with HCV genome sequences satisfying 70 per cent bootstrap and 0.05 genetic threshold distance requirement for sequence similarity. Logistic regression was used to assess sociodemographic factors associated with being in an HCV cluster. Among 512 HCV genotype 1 viremic PWID, HCV subtype prevalence was 83 per cent genotype 1a and 17 per cent genotype 1b. The median age of participants was 54 years, 68 per cent male, 87 per cent Black, and 38 per cent HIV infected. Overall, 9 per cent (n = 44) were grouped into 21 clusters, consisting of 20 pairs and 1 triad. Of the 425 genotype 1a and 87 genotype 1b samples evaluated, 8 per cent (n = 33) and 13 per cent (n = 11) respectively, were in clusters. In unadjusted analyses, membership in a cluster, was associated with younger age (odds ratio (OR) 1.5 [95% confidence interval (CI) 1.1–2.1] per 10 year age decrease); female sex (OR 2.8 [95% CI 1.5–5.3]), HIV infection (OR 4.9 [95% CI 2.5–9.9]), and living in East Baltimore (versus outside East Baltimore, OR 2.0 [95% CI 1.0–3.9]). In adjusted analyses, female sex (OR 2.0 [95% CI 1.0–3.9] and HIV infection (OR 5.4 [95% CI 2.6–11.1] remained independently associated with being in an HCV infection cluster. HIV-infected PWID and their networks should be prioritized for HCV treatment and prevention interventions given an increased likelihood of transmission in these groups.

Abstract Hepatitis C virus (HCV) infection occurs in 30–90 per cent of people who inject drugs (PWID). Although cure rates can exceed 95 per cent, treatment access is limited and approximately 400,000 people die each year due to complications of chronic infection. A temporal analysis of cluster networks among PWID can be used to inform strategies to interdict transmission. In Baltimore, PWID have been recruited for The AIDS Linked to the IntraVenous Experience (ALIVE) cohort. A demographic questionnaire was administered and recorded for baseline and recent participants. Viral RNA underwent PCR with primers targeting the core and envelope-1 protein (CE1) and sequenced via Sanger sequencing. Sequences with > 400 bp reads and Q-scores >370 were used for downstream analysis resulting in 322 ALIVE baseline participants (1988–9) and 548 recently diagnosed subjects enrolled approximately two decades later (2005–16). Cluster networks were rendered with a threshold of 4 per cent in MicrobeTRACE, and statistical analyses were performed in R Studio. Of the 1988–9 subjects, the majority (259/317, 81.7%) were a part of cluster. There were nine clusters and fifty-eight singletons, with two large clusters containing most sequences of genotype 1a (73.5%). Two decades later, a minority of recently diagnosed individuals (235/512, 44.1%) were part of a cluster. There were seventeen clusters with 286 singletons with two large clusters containing 1a genotype individuals (21.5%). Additional clustering was done by parsing the two datasets by subtype 1a (n = 714) and 1b (n = 151). The genotype 1a network demonstrates a majority, 65.8 per cent, of participants in clusters. Moreover, two large clusters can be observed with baseline participants towards the center and recent participants on the outskirts indicative of high linkage at baseline. The genotype 1b network produced a single large cluster but subclusters were observed. The sequences between the two time points co-mingled but subclusters were also observed. Interestingly, the two large clusters from 1988 to 1989 were still evident in the 2005–16 viral sequences. We observed greater cluster diversity in more recently diagnosed individuals, indicative of a less connected network of individuals sharing transmission risk, though major viral strains did persist over time in this cohort.

This study examines correlates of hepatitis C virus (HCV) genetic clustering among community-recruited people who inject drugs enrolled in the AIDS Linked to the IntraVenous Experience cohort in Baltimore between 1988 and 1989. HCV RNA was extracted and the core/envelope-1 region was sequenced. Clusters were identified from maximum likelihood trees with 1000 bootstrap replicates using a 70% aLRT and a 4% genetic-distance threshold in Cluster Picker. Overall, 46% of participants were in a cluster, including 122 genotype-1a and 36 genotype-1b clusters with an average of 2-3 genetically linked HCV infections. The largest cluster consists of 9 participants. In univariable analysis, black race (PR= 1.66 [95% CI: 1.12-2.45]), age < 35 years (PR= 1.18 [95% CI: 1.02-1.37]), and injection drug use of cocaine alone (PR= 1.30 [95% CI: 1.02-1.65]) were significantly associated with being in a cluster. Conversely, a history of medication-associated treatment (MAT) was negatively associated with being in a cluster (PR= 0.82 [95% CI: 0.71-0.95]). In multivariable analysis, black race (APR = 1.62 [95% CI: 1.11-2.38]) remained independently associated being in a cluster while MAT (APR = 0.85 [95% CI: 0.74-0.99]) remained negatively associated with clustering. Our findings suggest strong locally-propagated transmission networks during the early epidemic that was driven by younger PWID. In light of the current opioid epidemic in the US, these findings suggest an urgent need for preventive interventions to mitigate the growth of large HCV transmission networks.

Using human fecal samples and including sequencing for both bacterial and fungal taxa, this study compared the conventional antibiotics used to treat C. difficile infection (CDI) from the perspective of the microbiome, which is particularly relevant, given the relationship between dysbiotic states and the development of CDI. Sequencing and imputed functional analyses suggest that C. difficile -directed antibiotics are associated with distinct forms of dysbiosis that may be influential in the course of CDI. Further, a role for fungal organisms in the perpetuation of the causal dysbiosis of CDI is discussed, suggesting a previously unappreciated, clinically relevant transkingdom interaction that warrants further study.

This study sought to characterize the bacterial and fungal microbiota changes associated with Clostridium difficile infection (CDI) among inpatients with diarrhea, in order to further explain the pathogenesis of this infection as well as to potentially guide new CDI therapies. Twenty-four inpatients with diarrhea were enrolled, 12 of whom had CDI. Each patient underwent stool testing for CDI prior to being treated with difficile-directed antibiotics, when appropriate. Clinical data was obtained from the medical record, while each stool sample underwent 16S rRNA and ITS sequencing for bacterial and fungal elements. An analysis of microbial community structures distinct to the CDI population was also performed. The results demonstrated no difference between the CDI and non-CDI cohorts with respect to any previously reported CDI risk factors. Butyrogenic bacteria were enriched in both CDI and non-CDI patients. A previously unreported finding of increased numbers of Akkermansia muciniphila in CDI patients was observed, an organism which degrades mucin and which therefore may provide a selective advantage toward CDI. Fungal elements of the genus Penicillium were predominant in CDI; these organisms produce antibacterial chemicals which may resist recovery of healthy microbiota. The most frequent CDI microbial community networks involved Peptostreptococcaceae and Enterococcus, with decreased population density of Bacteroides. These results suggest that the development of CDI is associated with microbiota changes which are consistently associated with CDI in human subjects. These gut taxa contribute to the intestinal dysbiosis associated with C. difficile infection.