Abstract
Experimental models of the central nervous system (CNS) are imperative for developmental and pathophysiological studies of neurological diseases. Among these models, three-dimensional (3D) induced pluripotent stem cell (iPSC)-derived brain organoid models have been successful in mitigating some of the drawbacks of 2D models; however, they are plagued by high organoid-to-organoid variability, making it difficult to compare specific gene regulatory pathways across 3D organoids with those of the native brain. Single-cell RNA sequencing (scRNA-seq) transcriptome datasets have recently emerged as powerful tools to perform integrative analyses and compare variability across organoids. However, transcriptome studies focusing on late-stage neural functionality development have been underexplored. Here, we combine and analyze 8 brain organoid transcriptome databases to study the correlation between differentiation protocols and their resulting cellular functionality across various 3D organoid and exogenous brain models. We utilize dimensionality reduction methods including principal component analysis (PCA) and uniform manifold approximation projection (UMAP) to identify and visualize cellular diversity among 3D models and subsequently use gene set enrichment analysis (GSEA) and developmental trajectory inference to quantify neuronal behaviors such as axon guidance, synapse transmission and action potential. We showed high similarity in cellular composition, cellular differentiation pathways and expression of functional genes in human brain organoids during induction and differentiation phases, i.e., up to 3 months in culture. However, during the maturation phase, i.e., 6-month timepoint, we observed significant developmental deficits and depletion of neuronal and astrocytes functional genes as indicated by our GSEA results. Our results caution against use of organoids to model pathophysiology and drug response at this advanced time point and provide insights to tune in vitro iPSC differentiation protocols to achieve desired neuronal functionality and improve current protocols.