Investigating trait variability of gene co-expression network architecture in brain by controlling for genomic risk of schizophrenia.

The effect of schizophrenia (SCZ) genetic risk on gene expression in brain remains elusive. A popular approach to this problem has been the application of gene co-expression network algorithms (e.g., WGCNA). To improve reliability with this method it is critical to remove unwanted sources of variance while also preserving biological signals of interest. In this WCGNA study of RNA-Seq data from postmortem prefrontal cortex (78 neurotypical donors, EUR ancestry), we tested the effects of SCZ genetic risk on co-expression networks. Specifically, we implemented a novel design in which gene expression was adjusted by linear regression models to preserve or remove variance explained by biological signal of interest (GWAS genomic scores for SCZ risk-(GS-SCZ), and genomic scores- GS of height (GS-Ht) as a negative control), while removing variance explained by covariates of non-interest. We calculated co-expression networks from adjusted expression (GS-SCZ and GS-Ht preserved or removed), and consensus between them (representative of a "background" network free of genomic scores effects). We then tested the overlap between GS-SCZ preserved modules and background networks reasoning that modules with reduced overlap would be most affected by GS-SCZ biology. Additionally, we tested these modules for convergence of SCZ risk (i.e., enrichment in PGC3 SCZ GWAS priority genes, enrichment in SCZ risk heritability and relevant biological ontologies. Our results highlight key aspects of GS-SCZ effects on brain co-expression networks, specifically: 1) preserving/removing SCZ genetic risk alters the co-expression modules; 2) biological pathways enriched in modules affected by GS-SCZ implicate processes of transcription, translation and metabolism that converge to influence synaptic transmission; 3) priority PGC3 SCZ GWAS genes and SCZ risk heritability are enriched in modules associated with GS-SCZ effects. Overall, our results indicate that gene co-expression networks that selectively integrate information about genetic risk can reveal novel combinations of biological pathways involved in schizophrenia.

Gene expression differences associated with alcohol use disorder in human brain.

Excessive alcohol consumption is a leading cause of preventable death worldwide. To improve understanding of neurobiological mechanisms associated with alcohol use disorder (AUD) in humans, we compared gene expression data from deceased individuals with and without AUD across two addiction-relevant brain regions: the nucleus accumbens (NAc) and dorsolateral prefrontal cortex (DLPFC). Bulk RNA-seq data from NAc and DLPFC (N ≥50 with AUD, ≥46 non-AUD) were analyzed for differential gene expression using modified negative binomial regression adjusting for technical and biological covariates. The region-level results were meta-analyzed with those from an independent dataset (NNAc = 28 AUD, 29 non-AUD; NPFC = 66 AUD, 77 non-AUD). We further tested for heritability enrichment of AUD-related phenotypes, gene co-expression networks, gene ontology enrichment, and drug repurposing. We identified 176 differentially expressed genes (DEGs; 12 in both regions, 78 in NAc only, 86 in DLPFC only) for AUD in our new dataset. After meta-analyzing with published data, we identified 476 AUD DEGs (25 in both regions, 29 in NAc only, 422 in PFC only). Of these DEGs, 17 were significant when looked up in GWAS of problematic alcohol use or drinks per week. Gene co-expression analysis showed both concordant and unique gene networks across brain regions. We also identified 29 and 436 drug compounds that target DEGs from our meta-analysis in NAc and PFC, respectively. This study identified robust AUD-associated DEGs, contributing novel neurobiological insights into AUD and highlighting genes targeted by known drug compounds, generating opportunity for drug repurposing to treat AUD.

Characterizing the nuclear and cytoplasmic transcriptomes in developing and mature human cortex uncovers new insight into psychiatric disease gene regulation.

Transcriptome compartmentalization by the nuclear membrane provides both stochastic and functional buffering of transcript activity in the cytoplasm, and has recently been implicated in neurodegenerative disease processes. Although many mechanisms regulating transcript compartmentalization are also prevalent in brain development, the extent to which subcellular localization differs as the brain matures has yet to be addressed. To characterize the nuclear and cytoplasmic transcriptomes during brain development, we sequenced both RNA fractions from homogenate prenatal and adult human postmortem cortex using poly(A)+ and Ribo-Zero library preparation methods. We find that while many genes are differentially expressed by fraction and developmental expression changes are similarly detectable in nuclear and cytoplasmic RNA, the compartmented transcriptomes become more distinct as the brain matures, perhaps reflecting increased utilization of nuclear retention as a regulatory strategy in adult brain. We examined potential mechanisms of this developmental divergence including alternative splicing, RNA editing, nuclear pore composition, RNA-binding protein motif enrichment, and RNA secondary structure. Intron retention is associated with greater nuclear abundance in a subset of transcripts, as is enrichment for several splicing factor binding motifs. Finally, we examined disease association with fraction-regulated gene sets and found nuclear-enriched genes were also preferentially enriched in gene sets associated with neurodevelopmental psychiatric disorders. These results suggest that although gene-level expression is globally comparable between fractions, nuclear retention of transcripts may play an underappreciated role in developmental regulation of gene expression in brain, particularly in genes whose dysregulation is related to neuropsychiatric disorders.

Identification and prioritization of gene sets associated with schizophrenia risk by co-expression network analysis in human brain.

Schizophrenia polygenic risk is plausibly manifested by complex transcriptional dysregulation in the brain, involving networks of co-expressed and functionally related genes. The main purpose of this study was to identify and prioritize co-expressed gene sets in a hierarchical manner, based on the strength of the relationships with clinical diagnosis and with polygenic risk for schizophrenia. Weighted Gene Co-expression Network Analysis (WGCNA) was applied to RNA-quality-adjusted DLPFC RNA-Seq data from the LIBD Postmortem Human Brain Repository (90 controls, 74 schizophrenia cases; all Caucasians) to construct co-expression networks and detect "modules" of co-expressed genes. After multiple internal and external validation procedures, modules of selected interest were tested for enrichment in biological ontologies, for association with schizophrenia polygenic risk scores (PRSs) and with diagnosis, and also for enrichment in genes within the significant GWAS loci reported by the Psychiatric Genomic Consortium (PGC2). The association between schizophrenia genetic signals and modules of co-expression converged on one module showing not only a significant association with both diagnosis and PRS but also significant overlap with 36 PGC2 loci genes, deemed the strongest candidates for drug targets. This module contained many genes involved in synaptic signaling and neuroplasticity. Fifty-three PGC2 genes were in modules associated only with diagnosis and 59 in modules unrelated to diagnosis or PRS. Our study highlights complex relationships between gene co-expression networks in the brain and clinical state and polygenic risk for SCZ and provides a strategy for using this information in selecting and prioritizing potentially targetable gene sets for therapeutic drug development.

Developmental effects of maternal smoking during pregnancy on the human frontal cortex transcriptome.

Cigarette smoking during pregnancy is a major public health concern. While there are well-described consequences in early child development, there is very little known about the effects of maternal smoking on human cortical biology during prenatal life. We therefore performed a genome-wide differential gene expression analysis using RNA sequencing (RNA-seq) on prenatal (N = 33; 16 smoking-exposed) as well as adult (N = 207; 57 active smokers) human postmortem prefrontal cortices. Smoking exposure during the prenatal period was directly associated with differential expression of 14 genes; in contrast, during adulthood, despite a much larger sample size, only two genes showed significant differential expression (FDR < 10%). Moreover, 1,315 genes showed significantly different exposure effects between maternal smoking during pregnancy and direct exposure in adulthood (FDR < 10%)-these differences were largely driven by prenatal differences that were enriched for pathways previously implicated in addiction and synaptic function. Furthermore, prenatal and age-dependent differentially expressed genes were enriched for genes implicated in non-syndromic autism spectrum disorder (ASD) and were differentially expressed as a set between patients with ASD and controls in postmortem cortical regions. These results underscore the enhanced sensitivity to the biological effect of smoking exposure in the developing brain and offer insight into how maternal smoking during pregnancy affects gene expression in the prenatal human cortex. They also begin to address the relationship between in utero exposure to smoking and the heightened risks for the subsequent development of neuropsychiatric disorders.

Schizophrenia risk variants influence multiple classes of transcripts of sorting nexin 19 (SNX19).

Genome-wide association studies (GWAS) have identified many genomic loci associated with risk for schizophrenia, but unambiguous identification of the relationship between disease-associated variants and specific genes, and in particular their effect on risk conferring transcripts, has proven difficult. To better understand the specific molecular mechanism(s) at the schizophrenia locus in 11q25, we undertook cis expression quantitative trait loci (cis-eQTL) mapping for this 2 megabase genomic region using postmortem human brain samples. To comprehensively assess the effects of genetic risk upon local expression, we evaluated multiple transcript features: genes, exons, and exon-exon junctions in multiple brain regions-dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Genetic risk variants strongly associated with expression of SNX19 transcript features that tag multiple rare classes of SNX19 transcripts, whereas they only weakly affected expression of an exon-exon junction that tags the majority of abundant transcripts. The most prominent class of SNX19 risk-associated transcripts is predicted to be overexpressed, defined by an exon-exon splice junction between exons 8 and 10 (junc8.10) and that is predicted to encode proteins that lack the characteristic nexin C terminal domain. Risk alleles were also associated with either increased or decreased expression of multiple additional classes of transcripts. With RACE, molecular cloning, and long read sequencing, we found a number of novel SNX19 transcripts that further define the set of potential etiological transcripts. We explored epigenetic regulation of SNX19 expression and found that DNA methylation at CpG sites near the primary transcription start site and within exon 2 partially mediate the effects of risk variants on risk-associated expression. ATAC sequencing revealed that some of the most strongly risk-associated SNPs are located within a region of open chromatin, suggesting a nearby regulatory element is involved. These findings indicate a potentially complex molecular etiology, in which risk alleles for schizophrenia generate epigenetic alterations and dysregulation of multiple classes of SNX19 transcripts.

qSVA framework for RNA quality correction in differential expression analysis.

RNA sequencing (RNA-seq) is a powerful approach for measuring gene expression levels in cells and tissues, but it relies on high-quality RNA. We demonstrate here that statistical adjustment using existing quality measures largely fails to remove the effects of RNA degradation when RNA quality associates with the outcome of interest. Using RNA-seq data from molecular degradation experiments of human primary tissues, we introduce a method-quality surrogate variable analysis (qSVA)-as a framework for estimating and removing the confounding effect of RNA quality in differential expression analysis. We show that this approach results in greatly improved replication rates (>3×) across two large independent postmortem human brain studies of schizophrenia and also removes potential RNA quality biases in earlier published work that compared expression levels of different brain regions and other diagnostic groups. Our approach can therefore improve the interpretation of differential expression analysis of transcriptomic data from human tissue.

Comprehensive assessment of multiple biases in small RNA sequencing reveals significant differences in the performance of widely used methods.

BACKGROUND: RNA sequencing offers advantages over other quantification methods for microRNA (miRNA), yet numerous biases make reliable quantification challenging. Previous evaluations of these biases have focused on adapter ligation bias with limited evaluation of reverse transcription bias or amplification bias. Furthermore, evaluations of the quantification of isomiRs (miRNA isoforms) or the influence of starting amount on performance have been very limited. No study had yet evaluated the quantification of isomiRs of altered length or compared the consistency of results derived from multiple moderate starting inputs. We therefore evaluated quantifications of miRNA and isomiRs using four library preparation kits, with various starting amounts, as well as quantifications following removal of duplicate reads using unique molecular identifiers (UMIs) to mitigate reverse transcription and amplification biases. RESULTS: All methods resulted in false isomiR detection; however, the adapter-free method tested was especially prone to false isomiR detection. We demonstrate that using UMIs improves accuracy and we provide a guide for input amounts to improve consistency. CONCLUSIONS: Our data show differences and limitations of current methods, thus raising concerns about the validity of quantification of miRNA and isomiRs across studies. We advocate for the use of UMIs to improve accuracy and reliability of miRNA quantifications.

Integrated DNA methylation and gene expression profiling across multiple brain regions implicate novel genes in Alzheimer's disease.

Late-onset Alzheimer's disease (AD) is a complex age-related neurodegenerative disorder that likely involves epigenetic factors. To better understand the epigenetic state associated with AD, we surveyed 420,852 DNA methylation (DNAm) sites from neurotypical controls (N = 49) and late-onset AD patients (N = 24) across four brain regions (hippocampus, entorhinal cortex, dorsolateral prefrontal cortex and cerebellum). We identified 858 sites with robust differential methylation collectively annotated to 772 possible genes (FDR < 5%, within 10 kb). These sites were overrepresented in AD genetic risk loci (p = 0.00655) and were enriched for changes during normal aging (p < 2.2 × 10-16), and nearby genes were enriched for processes related to cell-adhesion, immunity, and calcium homeostasis (FDR < 5%). To functionally validate these associations, we generated and analyzed corresponding transcriptome data to prioritize 130 genes within 10 kb of the differentially methylated sites. These 130 genes were differentially expressed between AD cases and controls and their expression was associated with nearby DNAm (p < 0.05). This integrated analysis implicates novel genes in Alzheimer's disease, such as ANKRD30B. These results highlight DNAm differences in Alzheimer's disease that have gene expression correlates, further implicating DNAm as an epigenetic mechanism underlying pathological molecular changes associated with AD. Furthermore, our framework illustrates the value of integrating epigenetic and transcriptomic data for understanding complex disease.

Genomic structure and expression of the human serotonin 2A receptor gene (HTR2A) locus: identification of novel HTR2A and antisense (HTR2A-AS1) exons.

BACKGROUND: The serotonin 2A receptor is widely implicated in genetic association studies and remains an important drug target for psychiatric, neurological, and cardiovascular conditions. RNA sequencing redefined the architecture of the serotonin 2A receptor gene (HTR2A), revealing novel mRNA transcript isoforms utilizing unannotated untranslated regions of the gene. Expression of these untranslated regions is modulated by common single nucleotide polymorphisms (SNPs), namely rs6311. Previous studies did not fully capture the complexity of the sense- and antisense-encoded transcripts with respect to novel exons in the HTR2A gene locus. Here, we comprehensively catalogued exons and RNA isoforms for both HTR2A and HTR2A-AS1 using RNA-Seq from human prefrontal cortex and multiple mouse tissues. We subsequently tested associations between expression of newfound gene features and common SNPs in humans. RESULTS: We find that the human HTR2A gene spans ~66 kilobases and consists of 7, rather than 4 exons. Furthermore, the revised human HTR2A-AS1 gene spans ~474 kilobases and consists of 18, rather than 3 exons. Three HTR2A exons directly overlap with HTR2A-AS1 exons, suggesting potential for complementary nucleotide interactions. The repertoire of possible mouse Htr2a splice isoforms is remarkably similar to humans and we also find evidence for overlapping sense-antisense transcripts in the same relative positions as the human transcripts. rs6311 and SNPs in high linkage disequilibrium are associated with HTR2A-AS1 expression, in addition to previously described associations with expression of the extended 5' untranslated region of HTR2A. CONCLUSIONS: Our proposed HTR2A and HTR2A-AS1 gene structures dramatically differ from current annotations, now including overlapping exons on the sense and anti-sense strands. We also find orthologous transcript isoforms expressed in mice, providing opportunities to elucidate the biological roles of the human isoforms using a model system. Associations between rs6311 and expression of HTR2A and HTR2A-AS1 suggest this polymorphism is capable of modulating the expression of the sense or antisense transcripts. Still unclear is whether these SNPs act directly on the expression of the sense or antisense transcripts and whether overlapping exons are capable of interacting through complimentary base-pairing. Additional studies are necessary to determine the extent and nature of interactions between the SNPs and the transcripts prior to interpreting these findings in the context of phenotypes associated with HTR2A.

RNA-Seq-based toxicogenomic assessment of fresh frozen and formalin-fixed tissues yields similar mechanistic insights.

Formalin-fixed, paraffin-embedded (FFPE) pathology specimens represent a potentially vast resource for transcriptomic-based biomarker discovery. We present here a comparison of results from a whole transcriptome RNA-Seq analysis of RNA extracted from fresh frozen and FFPE livers. The samples were derived from rats exposed to aflatoxin B1 (AFB1 ) and a corresponding set of control animals. Principal components analysis indicated that samples were separated in the two groups representing presence or absence of chemical exposure, both in fresh frozen and FFPE sample types. Sixty-five percent of the differentially expressed transcripts (AFB1 vs. controls) in fresh frozen samples were also differentially expressed in FFPE samples (overlap significance: P < 0.0001). Genomic signature and gene set analysis of AFB1 differentially expressed transcript lists indicated highly similar results between fresh frozen and FFPE at the level of chemogenomic signatures (i.e., single chemical/dose/duration elicited transcriptomic signatures), mechanistic and pathology signatures, biological processes, canonical pathways and transcription factor networks. Overall, our results suggest that similar hypotheses about the biological mechanism of toxicity would be formulated from fresh frozen and FFPE samples. These results indicate that phenotypically anchored archival specimens represent a potentially informative resource for signature-based biomarker discovery and mechanistic characterization of toxicity.

Network-wide risk convergence in gene co-expression identifies reproducible genetic hubs of schizophrenia risk.

The omnigenic model posits that genetic risk for traits with complex heritability involves cumulative effects of peripheral genes on mechanistic "core genes," suggesting that in a network of genes, those closer to clusters including core genes should have higher GWAS signals. In gene co-expression networks, we confirmed that GWAS signals accumulate in genes more connected to risk-enriched gene clusters, highlighting across-network risk convergence. This was strongest in adult psychiatric disorders, especially schizophrenia (SCZ), spanning 70% of network genes, suggestive of super-polygenic architecture. In snRNA-seq cell type networks, SCZ risk convergence was strongest in L2/L3 excitatory neurons. We prioritized genes most connected to SCZ-GWAS genes, which showed robust association to a CRISPRa measure of PGC3 regulation and were consistently identified across several brain regions. Several genes, including dopamine-associated ones, were prioritized specifically in the striatum. This strategy thus retrieves current drug targets and can be used to prioritize other potential drug targets.

Mapping recurrent mosaic copy number variation in human neurons.

When somatic cells acquire complex karyotypes, they often are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons are limited by relatively small sample sizes. Here, we develop an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We apply this approach to 2,125 frontal cortical neurons from a neurotypical human brain. SCOVAL identifies 226 CNV neurons, which include a subclass of 65 CNV neurons with highly aberrant karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we find that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contain fewer, but longer, genes.

Sex affects transcriptional associations with schizophrenia across the dorsolateral prefrontal cortex, hippocampus, and caudate nucleus.

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

Analysis of gene expression in the postmortem brain of neurotypical Black Americans reveals contributions of genetic ancestry.

Ancestral differences in genomic variation affect the regulation of gene expression; however, most gene expression studies have been limited to European ancestry samples or adjusted to identify ancestry-independent associations. Here, we instead examined the impact of genetic ancestry on gene expression and DNA methylation in the postmortem brain tissue of admixed Black American neurotypical individuals to identify ancestry-dependent and ancestry-independent contributions. Ancestry-associated differentially expressed genes (DEGs), transcripts and gene networks, while notably not implicating neurons, are enriched for genes related to the immune response and vascular tissue and explain up to 26% of heritability for ischemic stroke, 27% of heritability for Parkinson disease and 30% of heritability for Alzheimer's disease. Ancestry-associated DEGs also show general enrichment for the heritability of diverse immune-related traits but depletion for psychiatric-related traits. We also compared Black and non-Hispanic white Americans, confirming most ancestry-associated DEGs. Our results delineate the extent to which genetic ancestry affects differences in gene expression in the human brain and the implications for brain illness risk.

TrkB-dependent regulation of molecular signaling across septal cell types.

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning, and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

Dopamine signaling enriched striatal gene set predicts striatal dopamine synthesis and physiological activity in vivo.

The polygenic architecture of schizophrenia implicates several molecular pathways involved in synaptic function. However, it is unclear how polygenic risk funnels through these pathways to translate into syndromic illness. Using tensor decomposition, we analyze gene co-expression in the caudate nucleus, hippocampus, and dorsolateral prefrontal cortex of post-mortem brain samples from 358 individuals. We identify a set of genes predominantly expressed in the caudate nucleus and associated with both clinical state and genetic risk for schizophrenia that shows dopaminergic selectivity. A higher polygenic risk score for schizophrenia parsed by this set of genes predicts greater dopamine synthesis in the striatum and greater striatal activation during reward anticipation. These results translate dopamine-linked genetic risk variation into in vivo neurochemical and hemodynamic phenotypes in the striatum that have long been implicated in the pathophysiology of schizophrenia.

Individual variation in the emergence of anterior-to-posterior neural fates from human pluripotent stem cells.

Variability between human pluripotent stem cell (hPSC) lines remains a challenge and opportunity in biomedicine. In this study, hPSC lines from multiple donors were differentiated toward neuroectoderm and mesendoderm lineages. We revealed dynamic transcriptomic patterns that delineate the emergence of these lineages, which were conserved across lines, along with individual line-specific transcriptional signatures that were invariant throughout differentiation. These transcriptomic signatures predicted an antagonism between SOX21-driven forebrain fates and retinoic acid-induced hindbrain fates. Replicate lines and paired adult tissue demonstrated the stability of these line-specific transcriptomic traits. We show that this transcriptomic variation in lineage bias had both genetic and epigenetic origins, aligned with the anterior-to-posterior structure of early mammalian development, and was present across a large collection of hPSC lines. These findings contribute to developing systematic analyses of PSCs to define the origin and consequences of variation in the early events orchestrating individual human development.

Central neurocytoma exhibits radial glial cell signatures with FGFR3 hypomethylation and overexpression.

We explored the genomic events underlying central neurocytoma (CN), a rare neoplasm of the central nervous system, via multiomics approaches, including whole-exome sequencing, bulk and single-nuclei RNA sequencing, and methylation sequencing. We identified FGFR3 hypomethylation leading to FGFR3 overexpression as a major event in the ontogeny of CN that affects crucial downstream events, such as aberrant PI3K-AKT activity and neuronal development pathways. Furthermore, we found similarities between CN and radial glial cells based on analyses of gene markers and CN tumor cells and postulate that CN tumorigenesis is due to dysregulation of radial glial cell differentiation into neurons. Our data demonstrate the potential role of FGFR3 as one of the leading drivers of tumorigenesis in CN.

Human archetypal pluripotent stem cells differentiate into trophoblast stem cells via endogenous BMP5/7 induction without transitioning through naive state.

Primary human trophoblast stem cells (TSCs) and TSCs derived from human pluripotent stem cells (hPSCs) can potentially model placental processes in vitro. Yet, the pluripotent states and factors involved in the differentiation of hPSCs to TSCs remain poorly understood. In this study, we demonstrate that the primed pluripotent state can generate TSCs by activating pathways such as Epidermal Growth Factor (EGF) and Wingless-related integration site (WNT), and by suppressing tumor growth factor beta (TGFß), histone deacetylases (HDAC), and Rho-associated protein kinase (ROCK) signaling pathways, all without the addition of exogenous Bone morphogenetic protein 4 (BMP4)-a condition we refer to as the TS condition. We characterized this process using temporal single-cell RNA sequencing to compare TS conditions with differentiation protocols involving BMP4 activation alone or BMP4 activation in conjunction with WNT inhibition. The TS condition consistently produced a stable, proliferative cell type that closely mimics first-trimester placental cytotrophoblasts, marked by the activation of endogenous retroviral genes and the absence of amnion expression. This was observed across multiple cell lines, including various primed induced pluripotent stem cell (iPSC) and embryonic stem cell (ESC) lines. Primed-derived TSCs can proliferate for over 30 passages and further specify into multinucleated syncytiotrophoblasts and extravillous trophoblast cells. Our research establishes that the differentiation of primed hPSCs to TSC under TS conditions triggers the induction of TMSB4X, BMP5/7, GATA3, and TFAP2A without progressing through a naive state. These findings propose that the primed hPSC state is part of a continuum of potency with the capacity to differentiate into TSCs through multiple routes.