Effects of iron supplements and iron-containing micronutrient powders on the gut microbiome in Bangladeshi infants: a randomized controlled trial.

Anemia is highly prevalent globally, especially in young children in low-income countries, where it often overlaps with a high burden of diarrheal disease. Distribution of iron interventions (as supplements or iron-containing multiple micronutrient powders, MNPs) is a key anemia reduction strategy. Small studies in Africa indicate iron may reprofile the gut microbiome towards pathogenic species. We seek to evaluate the safety of iron and MNPs based on their effects on diversity, composition, and function of the gut microbiome in children in rural Bangladesh as part of a large placebo-controlled randomized controlled trial of iron or MNPs given for 3 months (ACTRN12617000660381). In 923 infants, we evaluate the microbiome before, immediately following, and nine months after interventions, using 16S rRNA gene sequencing and shotgun metagenomics in a subset. We identify no increase in diarrhea with either treatment. In our primary analysis, neither iron nor MNPs alter gut microbiome diversity or composition. However, when not adjusting for multiple comparisons, compared to placebo, children receiving iron and MNPs exhibit reductions in commensal species (e.g., Bifidobacterium, Lactobacillus) and increases in potential pathogens, including Clostridium. These increases are most evident in children with baseline iron repletion and are further supported by trend-based statistical analyses.

Castanet: A pipeline for rapid analysis of targeted multi-pathogen genomic data.

MOTIVATION: Target enrichment strategies generate genomic data from multiple pathogens in a single process, greatly improving sensitivity over metagenomic sequencing and enabling cost-effective, high throughput surveillance and clinical applications. However, uptake by research and clinical laboratories is constrained by an absence of computational tools that are specifically designed for the analysis of multi-pathogen enrichment sequence data. Here we present an analysis pipeline, Castanet, for use with multi-pathogen enrichment sequencing data. Castanet is designed to work with short-read data produced by existing targeted enrichment strategies, but can be readily deployed on any BAM file generated by another methodology. Also included are an optional graphical interface and installer script. RESULTS: In addition to genome reconstruction, Castanet reports method-specific metrics that enable quantification of capture efficiency, estimation of pathogen load, differentiation of low-level positives from contamination, and assessment of sequencing quality. Castanet can be used as a traditional end-to-end pipeline for consensus generation, but its strength lies in the ability to process a flexible, pre-defined set of pathogens of interest directly from multi-pathogen enrichment experiments. In our tests, Castanet consensus sequences were accurate reconstructions of reference sequences, including in instances where multiple strains of the same pathogen were present. Castanet performs effectively on standard computers and can process the entire output of a 96-sample enrichment sequencing run (50M reads) using a single batch process command, in $<$2 h. AVAILABILITY AND IMPLEMENTATION: Source code freely available under GPL-3 license at https://github.com/MultipathogenGenomics/castanet, implemented in Python 3.10 and supported in Ubuntu Linux 22.04. SUPPLEMENTARY INFORMATION: Supplementary data available at the journal's web site. Supporting sequencing data available at https://www.ebi.ac.uk/ena/browser/view/PRJEB77004.

Community use of oral antibiotics transiently reprofiles the intestinal microbiome in young Bangladeshi children.

Antibiotics may alter the gut microbiome, and this is one of the mechanisms by which antimicrobial resistance may be promoted. Suboptimal antimicrobial stewardship in Asia has been linked to antimicrobial resistance. We aim to examine the relationship between oral antibiotic use and composition and antimicrobial resistance in the gut microbiome in 1093 Bangladeshi infants. We leverage a trial of 8-month-old infants in rural Bangladesh: 61% of children were cumulatively exposed to antibiotics (most commonly cephalosporins and macrolides) over the 12-month study period, including 47% in the first 3 months of the study, usually for fever or respiratory infection. 16S rRNA amplicon sequencing in 11-month-old infants reveals that alpha diversity of the intestinal microbiome is reduced in children who received antibiotics within the previous 7 days; these samples also exhibit enrichment for Enterococcus and Escherichia/Shigella genera. No effect is seen in children who received antibiotics earlier. Using shotgun metagenomics, overall abundance of antimicrobial resistance genes declines over time. Enrichment for an Enterococcus-related antimicrobial resistance gene is observed in children receiving antibiotics within the previous 7 days, but not earlier. Presence of antimicrobial resistance genes is correlated to microbiome composition. In Bangladeshi children, community use of antibiotics transiently reprofiles the gut microbiome.

Characterising the molecular epidemiology of human parechovirus in young infants in the UK and Canada.

OBJECTIVES: We evaluated the extent of virus heterogeneity in PeV infected infants in the UK, Canada and Australia. METHODS: Samples were collected from PeV infected infants during 2013-16. Next generation sequencing was used to obtain sequencing data and construct phylogenetic trees based on analysis of the VP1 region. Comparison was made with sequencing data available from an outbreak in Australia. RESULTS: We amplified and sequenced 58 samples. All obtained PeV sequences were genotype 3 apart from one UK sample which was PeV-A5. Phylogenetic analysis revealed that all strains clustered together on the same clade and showed no significant genetic variation. We saw no significant evidence of association between sequence and either clinical severity (defined by admission to paediatric intensive care), geographical origin (compared between Canada and U.K) or year of sample collection (samples sequenced during 2013 - 2018). CONCLUSIONS: In this small cohort, sequencing data indicate that PeV circulating in the UK and Canada from 2013 to 18 are derived from a common ancestor. No association between disease severity and genetic sequence was seen in the UK or Canadian cohorts. Larger studies are required to support these findings.

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.