An immunohistochemical atlas of necroptotic pathway expression.

Necroptosis is a lytic form of regulated cell death reported to contribute to inflammatory diseases of the gut, skin and lung, as well as ischemic-reperfusion injuries of the kidney, heart and brain. However, precise identification of the cells and tissues that undergo necroptotic cell death in vivo has proven challenging in the absence of robust protocols for immunohistochemical detection. Here, we provide automated immunohistochemistry protocols to detect core necroptosis regulators - Caspase-8, RIPK1, RIPK3 and MLKL - in formalin-fixed mouse and human tissues. We observed surprising heterogeneity in protein expression within tissues, whereby short-lived immune barrier cells were replete with necroptotic effectors, whereas long-lived cells lacked RIPK3 or MLKL expression. Local changes in the expression of necroptotic effectors occurred in response to insults such as inflammation, dysbiosis or immune challenge, consistent with necroptosis being dysregulated in disease contexts. These methods will facilitate the precise localisation and evaluation of necroptotic signaling in vivo.

Targeted metagenomics reveals association between severity and pathogen co-detection in infants with respiratory syncytial virus.

Respiratory syncytial virus (RSV) is the leading cause of hospitalisation for respiratory infection in young children. RSV disease severity is known to be age-dependent and highest in young infants, but other correlates of severity, particularly the presence of additional respiratory pathogens, are less well understood. In this study, nasopharyngeal swabs were collected from two cohorts of RSV-positive infants <12 months in Spain, the UK, and the Netherlands during 2017-20. We show, using targeted metagenomic sequencing of >100 pathogens, including all common respiratory viruses and bacteria, from samples collected from 433 infants, that burden of additional viruses is common (111/433, 26%) but only modestly correlates with RSV disease severity. In contrast, there is strong evidence in both cohorts and across age groups that presence of Haemophilus bacteria (194/433, 45%) is associated with higher severity, including much higher rates of hospitalisation (odds ratio 4.25, 95% CI 2.03-9.31). There is no evidence for association between higher severity and other detected bacteria, and no difference in severity between RSV genotypes. Our findings reveal the genomic diversity of additional pathogens during RSV infection in infants, and provide an evidence base for future causal investigations of the impact of co-infection on RSV disease severity.

A risk-reward examination of sample multiplexing reagents for single cell RNA-Seq.

Single-cell RNA sequencing (scRNA-Seq) has emerged as a powerful tool for understanding cellular heterogeneity and function. However the choice of sample multiplexing reagents can impact data quality and experimental outcomes. In this study, we compared various multiplexing reagents, including MULTI-Seq, Hashtag antibody, and CellPlex, across diverse sample types such as human peripheral blood mononuclear cells (PBMCs), mouse embryonic brain and patient-derived xenografts (PDXs). We found that all multiplexing reagents worked well in cell types robust to ex vivo manipulation but suffered from signal-to-noise issues in more delicate sample types. We compared multiple demultiplexing algorithms which differed in performance depending on data quality. We find that minor improvements to laboratory workflows such as titration and rapid processing are critical to optimal performance. We also compared the performance of fixed scRNA-Seq kits and highlight the advantages of the Parse Biosciences kit for fragile samples. Highly multiplexed scRNA-Seq experiments require more sequencing resources, therefore we evaluated CRISPR-based destruction of non-informative genes to enhance sequencing value. Our comprehensive analysis provides insights into the selection of appropriate sample multiplexing reagents and protocols for scRNA-Seq experiments, facilitating more accurate and cost-effective studies.

scPipe: an extended preprocessing pipeline for comprehensive single-cell ATAC-Seq data integration in R/Bioconductor.

scPipe is a flexible R/Bioconductor package originally developed to analyse platform-independent single-cell RNA-Seq data. To expand its preprocessing capability to accommodate new single-cell technologies, we further developed scPipe to handle single-cell ATAC-Seq and multi-modal (RNA-Seq and ATAC-Seq) data. After executing multiple data cleaning steps to remove duplicated reads, low abundance features and cells of poor quality, a SingleCellExperiment object is created that contains a sparse count matrix with features of interest in the rows and cells in the columns. Quality control information (e.g. counts per cell, features per cell, total number of fragments, fraction of fragments per peak) and any relevant feature annotations are stored as metadata. We demonstrate that scPipe can efficiently identify 'true' cells and provides flexibility for the user to fine-tune the quality control thresholds using various feature and cell-based metrics collected during data preprocessing. Researchers can then take advantage of various downstream single-cell tools available in Bioconductor for further analysis of scATAC-Seq data such as dimensionality reduction, clustering, motif enrichment, differential accessibility and cis-regulatory network analysis. The scPipe package enables a complete beginning-to-end pipeline for single-cell ATAC-Seq and RNA-Seq data analysis in R.

Validation of an HIV whole genome sequencing method for HIV drug resistance testing in an Australian clinical microbiology laboratory.

Detection of HIV drug resistance (HIVDR) is vital to successful anti-retroviral therapy (ART). HIVDR testing to determine drug-resistance mutations is routinely performed in Australia to guide ART choice in newly diagnosed people living with HIV or in cases of treatment failure. In 2022, our clinical microbiology laboratory sought to validate a next-generation sequencing (NGS)-based HIVDR assay to replace the previous Sanger-sequencing (SS)-based ViroSeq. NGS solutions for HIVDR offer higher throughput, lower costs and higher sensitivity for variant detection. We sought to validate the previously described low-cost probe-based NGS method (veSEQ-HIV) for whole-genome recovery and HIVDR-testing in a diagnostic setting. veSEQ-HIV displayed 100% and 98% accuracy in major and minor mutation detection, respectively, and 100% accuracy of subtyping (provided > 1000 mapped reads were obtained). Pairwise comparison exhibited low inter-and intrarun variability across the whole-genome (Jaccard index [J] = 0.993; J = 0.972) and the Pol gene (J = 0.999; J = 0.999), respectively. veSEQ-HIV met all our pre-set criteria based on WHO recommendations and successfully replaced ViroSeq in our laboratory. Scaling-down veSEQ-HIV to a limited batch size and sequencing on Illumina iSeq. 100, allowed easy implementation of the assay into the workflow of a small sequencing laboratory with minimal staff and equipment and the ability to meet clinically relevant test turn-around times. As HIVDR-testing moves from SS- to NGS-based methods and new ART drugs come to market (particularly those with targets outside the Pol region), whole-genome recovery using veSEQ-HIV provides a robust, cost-effective and "future-proof" NGS method for HIVDR-testing.

Phylogenetic Analysis of Hepatitis C Virus Infections in a Large Belgian Cohort Using Next-Generation Sequencing of Full-Length Genomes.

The hepatitis C virus (HCV) epidemic in Western countries is primarily perpetuated by the sub-populations of men who have sex with men (MSM) and people who inject drugs (PWID). Understanding the dynamics of transmission in these communities is crucial for removing the remaining hurdles towards HCV elimination. We sequenced 269 annotated HCV plasma samples using probe enrichment and next-generation sequencing, obtaining 224 open reading frames of HCV (OR497849-OR498072). Maximum likelihood phylogenies were generated on the four most prevalent subtypes in this study (HCV1a, 1b, 3a, 4d) with a subsequent transmission cluster analysis. The highest rate of clustering was observed for HCV4d samples (13/17 (76.47%)). The second highest rate of clustering was observed in HCV1a samples (42/78 (53.85%)) with significant association with HIV-positive MSM. HCV1b and HCV3a had very low rates of clustering (2/83 (2.41%) and (0/29)). The spread of the prevalent subtype HCV1b appears to have been largely curtailed, and we demonstrate the onwards transmission of HCV1a and HCV4d in the HIV-positive MSM population across municipal borders. More systematic data collection and sequencing is needed to allow a better understanding of the HCV transmission among the community of PWID and overcome the remaining barriers for HCV elimination in Belgium.

Benchmarking long-read RNA-sequencing analysis tools using in silico mixtures.

The lack of benchmark data sets with inbuilt ground-truth makes it challenging to compare the performance of existing long-read isoform detection and differential expression analysis workflows. Here, we present a benchmark experiment using two human lung adenocarcinoma cell lines that were each profiled in triplicate together with synthetic, spliced, spike-in RNAs (sequins). Samples were deeply sequenced on both Illumina short-read and Oxford Nanopore Technologies long-read platforms. Alongside the ground-truth available via the sequins, we created in silico mixture samples to allow performance assessment in the absence of true positives or true negatives. Our results show that StringTie2 and bambu outperformed other tools from the six isoform detection tools tested, DESeq2, edgeR and limma-voom were best among the five differential transcript expression tools tested and there was no clear front-runner for performing differential transcript usage analysis between the five tools compared, which suggests further methods development is needed for this application.

High resolution imaging and five-year tuberculosis contact outcomes.

Background: The evolution of tuberculosis (TB) disease during the clinical latency period remains incompletely understood. Methods: 250 HIV-uninfected, adult household contacts of rifampicin-resistant TB with a negative symptom screen underwent baseline 18F-Fluorodeoxyglucose positron emission and computed tomography (PET/CT), repeated in 112 after 5-15 months. Following South African and WHO guidelines, participants did not receive preventive therapy. All participants had intensive baseline screening with spontaneous, followed by induced, sputum sampling and were then observed for an average of 4.7 years for culture-positive disease. Baseline PET/CT abnormalities were evaluated in relation to culture-positive disease. Results: At baseline, 59 (23.6%) participants had lung PET/CT findings consistent with TB of which 29 (11.6%) were defined as Subclinical TB, and 30 (12%) Subclinical TB-inactive. A further 83 (33.2%) had other lung parenchymal abnormalities and 108 (43.2%) had normal lungs. Over 1107-person years of follow-up 14 cases of culture-positive TB were diagnosed. Six cases were detected by intensive baseline screening, all would have been missed by the South African symptom-based screening strategy and only one detected by a WHO-recommended chest X-Ray screening strategy. Those with baseline Subclinical TB lesions on PET/CT were significantly more likely to be diagnosed with culture-positive TB over the study period, compared to those with normal lung parenchyma (10/29 [34.5%] vs 2/108 [1.9%], Hazard Ratio 22.37 [4.89-102.47, p<0.001]). Conclusions: These findings challenge the latent/active TB paradigm demonstrating that subclinical disease exists up to 4 years prior to microbiological detection and/or symptom onset. There are important implications for screening and management of TB.

Actinobacillus pleuropneumoniae encodes multiple phase-variable DNA methyltransferases that control distinct phasevarions.

Actinobacillus pleuropneumoniae is the cause of porcine pleuropneumonia, a severe respiratory tract infection that is responsible for major economic losses to the swine industry. Many host-adapted bacterial pathogens encode systems known as phasevarions (phase-variable regulons). Phasevarions result from variable expression of cytoplasmic DNA methyltransferases. Variable expression results in genome-wide methylation differences within a bacterial population, leading to altered expression of multiple genes via epigenetic mechanisms. Our examination of a diverse population of A. pleuropneumoniae strains determined that Type I and Type III DNA methyltransferases with the hallmarks of phase variation were present in this species. We demonstrate that phase variation is occurring in these methyltransferases, and show associations between particular Type III methyltransferase alleles and serovar. Using Pacific BioSciences Single-Molecule, Real-Time (SMRT) sequencing and Oxford Nanopore sequencing, we demonstrate the presence of the first ever characterised phase-variable, cytosine-specific Type III DNA methyltransferase. Phase variation of distinct Type III DNA methyltransferase in A. pleuropneumoniae results in the regulation of distinct phasevarions, and in multiple phenotypic differences relevant to pathobiology. Our characterisation of these newly described phasevarions in A. pleuropneumoniae will aid in the selection of stably expressed antigens, and direct and inform development of a rationally designed subunit vaccine against this major veterinary pathogen.

Full-genome next-generation sequencing of hepatitis C virus to assess the accuracy of genotyping by the commercial assay LiPA and the prevalence of resistance-associated substitutions in a Belgian cohort.

BACKGROUND: Although most currently used regimens for Hepatitis C virus (HCV) infections can be initiated without prior knowledge of genotype and subtype, genotyping is still useful to identify patients who might benefit from a personalized treatment due to resistance to direct-acting antivirals (DAA). OBJECTIVES: To assess the utility of full-genome next-generation sequencing (FG-NGS) for HCV genotyping. STUDY DESIGN: 138 HCV plasma samples previously genotyped by VERSANT HCV Genotype Assay (LiPA) were subjected to FG-NGS and phylogenetically genotyped Genome Detective. Consensuses were analysed by HCV-GLUE for resistance-associated substitutions (RASs) and their impact on treatment response was investigated. RESULTS: 102/138 (73.9%) samples were sequenced to a genome coverage and depth of >90% of the HCV open reading frame covered by >100 reads/site. Concordant genotype and subtype results were assigned in 97.1% and 79.4% of samples, respectively. FG-NGS resolved the subtype of 13.7% samples that had ambiguous calls by LiPA and identified one dual infection and one recombinant strain. At least one RAS was found for the HCV genes NS3, NS5A, and NS5B in 2.91%, 36.98% and 27.3% samples, respectively. Irrespective of the observed RAS, all patients responded well to DAA treatment, except for HCV1b-infected patients treated with Zepatier (33.3% failure rate (5/15)). CONCLUSION: While LiPA and FG-NGS showed overall good concordance, FG-NGS improved specificity for subtypes, recombinant and mixed infections. FG-NGS enabled the detection of RAS, but its predictive value for treatment outcome in DAA-naïve patients remains uncertain. With additional refinements, FG-NGS may be the way forward for HCV genotyping.

The interplay between community and hospital Enterococcus faecium clones within health-care settings: a genomic analysis.

BACKGROUND: The genomic relationships among Enterococcus faecium isolates are the subject of ongoing research that seeks to clarify the origins of observed lineages and the extent of horizontal gene transfer between them, and to robustly identify links between genotypes and phenotypes. E faecium is considered to form distinct groups-A and B-corresponding to isolates derived from patients who were hospitalised (A) and isolates from humans in the community (B). The additional separation of A into the so-called clades A1 and A2 remains an area of uncertainty. We aimed to investigate the relationships between A1 and non-A1 groups and explore the potential role of non-A1 isolates in shaping the population structure of hospital E faecium. METHODS: We collected short-read sequence data from invited groups that had previously published E faecium genome data. This hospital-based isolate collection could be separated into three groups (or clades, A1, A2, and B) by augmenting the study genomes with published sequences derived from human samples representing the previously defined genomic clusters. We performed phylogenetic analyses, by constructing maximum-likelihood phylogenetic trees, and identified historical recombination events. We assessed the pan-genome, did resistome analysis, and examined the genomic data to identify mobile genetic elements. Each genome underwent chromosome painting by use of ChromoPainter within FineSTRUCTURE software to assess ancestry and identify hybrid groups. We further assessed highly admixed regions to infer recombination directionality. FINDINGS: We assembled a collection of 1095 hospital E faecium sequences from 34 countries, further augmented by 33 published sequences. 997 (88%) of 1128 genomes clustered as A1, 92 (8%) as A2, and 39 (4%) as B. We showed that A1 probably emerged as a clone from within A2 and that, because of ongoing gene flow, hospital isolates currently identified as A2 represent a genetic continuum between A1 and community E faecium. This interchange of genetic material between isolates from different groups results in the emergence of hybrid genomes between clusters. Of the 1128 genomes, 49 (4%) hybrid genomes were identified: 33 previously labelled as A2 and 16 previously labelled as A1. These interactions were fuelled by a directional pattern of recombination mediated by mobile genetic elements. By contrast, the contribution of B group genetic material to A1 was limited to a few small regions of the genome and appeared to be driven by genomic sweep events. INTERPRETATION: A2 and B isolates coming into the hospital form an important reservoir for ongoing A1 adaptation, suggesting that effective long-term control of the effect of E faecium could benefit from strategies to reduce these genomic interactions, such as a focus on reducing the acquisition of hospital A1 strains by patients entering the hospital. FUNDING: Wellcome Trust.

Dual RNA sequencing reveals dendritic cell reprogramming in response to typhoidal Salmonella invasion.

Salmonella enterica represent a major disease burden worldwide. S. enterica serovar Typhi (S. Typhi) is responsible for potentially life-threatening Typhoid fever affecting 10.9 million people annually. While non-typhoidal Salmonella (NTS) serovars usually trigger self-limiting diarrhoea, invasive NTS bacteraemia is a growing public health challenge. Dendritic cells (DCs) are key professional antigen presenting cells of the human immune system. The ability of pathogenic bacteria to subvert DC functions and prevent T cell recognition contributes to their survival and dissemination within the host. Here, we adapted dual RNA-sequencing to define how different Salmonella pathovariants remodel their gene expression in tandem with that of infected DCs. We find DCs harness iron handling pathways to defend against invading Salmonellas, which S. Typhi is able to circumvent by mounting a robust response to nitrosative stress. In parallel, we uncover the alternative strategies invasive NTS employ to impair DC functions.

Viral genome wide association study identifies novel hepatitis C virus polymorphisms associated with sofosbuvir treatment failure.

Persistent hepatitis C virus (HCV) infection is a major cause of chronic liver disease, worldwide. With the development of direct-acting antivirals, treatment of chronically infected patients has become highly effective, although a subset of patients responds less well to therapy. Sofosbuvir is a common component of current de novo or salvage combination therapies, that targets the HCV NS5B polymerase. We use pre-treatment whole-genome sequences of HCV from 507 patients infected with HCV subtype 3a and treated with sofosbuvir containing regimens to detect viral polymorphisms associated with response to treatment. We find three common polymorphisms in non-targeted HCV NS2 and NS3 proteins are associated with reduced treatment response. These polymorphisms are enriched in post-treatment HCV sequences of patients unresponsive to treatment. They are also associated with lower reductions in viral load in the first week of therapy. Using in vitro short-term dose-response assays, these polymorphisms do not cause any reduction in sofosbuvir potency, suggesting an indirect mechanism of action in decreasing sofosbuvir efficacy. The identification of polymorphisms in NS2 and NS3 proteins associated with poor treatment outcomes emphasises the value of systematic genome-wide analyses of viruses in uncovering clinically relevant polymorphisms that impact treatment.

Distinct patterns of within-host virus populations between two subgroups of human respiratory syncytial virus.

Human respiratory syncytial virus (RSV) is a major cause of lower respiratory tract infection in young children globally, but little is known about within-host RSV diversity. Here, we characterised within-host RSV populations using deep-sequencing data from 319 nasopharyngeal swabs collected during 2017-2020. RSV-B had lower consensus diversity than RSV-A at the population level, while exhibiting greater within-host diversity. Two RSV-B consensus sequences had an amino acid alteration (K68N) in the fusion (F) protein, which has been associated with reduced susceptibility to nirsevimab (MEDI8897), a novel RSV monoclonal antibody under development. In addition, several minor variants were identified in the antigenic sites of the F protein, one of which may confer resistance to palivizumab, the only licensed RSV monoclonal antibody. The differences in within-host virus populations emphasise the importance of monitoring for vaccine efficacy and may help to explain the different prevalences of monoclonal antibody-escape mutants between the two subgroups.

High-resolution transcriptional landscape of xeno-free human induced pluripotent stem cell-derived cerebellar organoids.

Current protocols for producing cerebellar neurons from human pluripotent stem cells (hPSCs) often rely on animal co-culture and mostly exist as monolayers, limiting their capability to recapitulate the complex processes in the developing cerebellum. Here, we employed a robust method, without the need for mouse co-culture to generate three-dimensional cerebellar organoids from hPSCs that display hallmarks of in vivo cerebellar development. Single-cell profiling followed by comparison to human and mouse cerebellar atlases revealed the presence and maturity of transcriptionally distinct populations encompassing major cerebellar cell types. Encapsulation with Matrigel aimed to provide more physiologically-relevant conditions through recapitulation of basement-membrane signalling, influenced both growth dynamics and cellular composition of the organoids, altering developmentally relevant gene expression programmes. We identified enrichment of cerebellar disease genes in distinct cell populations in the hPSC-derived cerebellar organoids. These findings ascertain xeno-free human cerebellar organoids as a unique model to gain insight into cerebellar development and its associated disorders.

Targeted single-cell RNA sequencing of transcription factors enhances the identification of cell types and trajectories.

Single-cell RNA sequencing (scRNA-seq) is a widely used method for identifying cell types and trajectories in biologically heterogeneous samples, but it is limited in its detection and quantification of lowly expressed genes. This results in missing important biological signals, such as the expression of key transcription factors (TFs) driving cellular differentiation. We show that targeted sequencing of ∼1000 TFs (scCapture-seq) in iPSC-derived neuronal cultures greatly improves the biological information garnered from scRNA-seq. Increased TF resolution enhanced cell type identification, developmental trajectories, and gene regulatory networks. This allowed us to resolve differences among neuronal populations, which were generated in two different laboratories using the same differentiation protocol. ScCapture-seq improved TF-gene regulatory network inference and thus identified divergent patterns of neurogenesis into either excitatory cortical neurons or inhibitory interneurons. Furthermore, scCapture-seq revealed a role for of retinoic acid signaling in the developmental divergence between these different neuronal populations. Our results show that TF targeting improves the characterization of human cellular models and allows identification of the essential differences between cellular populations, which would otherwise be missed in traditional scRNA-seq. scCapture-seq TF targeting represents a cost-effective enhancement of scRNA-seq, which could be broadly applied to improve scRNA-seq resolution.

Short and long-read genome sequencing methodologies for somatic variant detection; genomic analysis of a patient with diffuse large B-cell lymphoma.

Recent advances in throughput and accuracy mean that the Oxford Nanopore Technologies PromethION platform is a now a viable solution for genome sequencing. Much of the validation of bioinformatic tools for this long-read data has focussed on calling germline variants (including structural variants). Somatic variants are outnumbered many-fold by germline variants and their detection is further complicated by the effects of tumour purity/subclonality. Here, we evaluate the extent to which Nanopore sequencing enables detection and analysis of somatic variation. We do this through sequencing tumour and germline genomes for a patient with diffuse B-cell lymphoma and comparing results with 150 bp short-read sequencing of the same samples. Calling germline single nucleotide variants (SNVs) from specific chromosomes of the long-read data achieved good specificity and sensitivity. However, results of somatic SNV calling highlight the need for the development of specialised joint calling algorithms. We find the comparative genome-wide performance of different tools varies significantly between structural variant types, and suggest long reads are especially advantageous for calling large somatic deletions and duplications. Finally, we highlight the utility of long reads for phasing clinically relevant variants, confirming that a somatic 1.6 Mb deletion and a p.(Arg249Met) mutation involving TP53 are oriented in trans.

The global dissemination of hospital clones of Enterococcus faecium.

BACKGROUND: The hospital-adapted A1 group of Enterococcus faecium remains an organism of significant concern in the context of drug-resistant hospital-associated infections. How this pathogen evolves and disseminates remains poorly understood. METHODS: A large, globally representative collection of short-read genomic data from the hospital-associated A1 group of Enterococcus faecium was assembled (n = 973). We analysed, using a novel analysis approach, global diversity in terms of both the dynamics of the accessory genome and homologous recombination among conserved genes. RESULTS: Two main modes of genomic evolution continue to shape E. faecium: the acquisition and loss of genes, including antimicrobial resistance genes, through mobile genetic elements including plasmids, and homologous recombination of the core genome. These events lead to new clones emerging at the local level, followed by the erosion of signals of clonality through recombination, and in some identifiable cases producing new clonal clusters. These patterns lead to new, emerging lineages which are able to spread globally over relatively short timeframes. CONCLUSIONS: The ability of A1 E. faecium to continually present new combinations of genes for potential selection suggests that controlling this pathogen will remain challenging but establishing a framework for understanding genomic evolution is likely to aid in tracking the threats posed by newly emerging lineages.

Targeted capture and sequencing of Orientia tsutsugamushi genomes from chiggers and humans.

Scrub typhus is a febrile disease caused by Orientia tsutsugamushi, transmitted by larval stage Trombiculid mites (chiggers), whose primary hosts are small mammals. The phylogenomics of O. tsutsugamushi in chiggers, small mammals and humans remains poorly understood. To combat the limitations imposed by the low relative quantities of pathogen DNA in typical O. tsutsugamushi clinical and ecological samples, along with the technical, safety and cost limitations of cell culture, a novel probe-based target enrichment sequencing protocol was developed. The method was designed to capture variation among conserved genes and facilitate phylogenomic analysis at the scale of population samples. A whole-genome amplification step was incorporated to enhance the efficiency of sequencing by reducing duplication rates. This resulted in on-target capture rates of up to 93% for a diverse set of human, chigger, and rodent samples, with the greatest success rate in samples with real-time PCR Ct values below 35. Analysis of the best-performing samples revealed phylogeographic clustering at local, provincial and international scales. Applying the methodology to a comprehensive set of samples could yield a more complete understanding of the ecology, genomic evolution and population structure of O. tsutsugamushi and other similarly challenging organisms, with potential benefits in the development of diagnostic tests and vaccines.