Spatiotemporal Gene Expression Research Articles

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model.

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

Impeller: a path-based heterogeneous graph learning method for spatial transcriptomic data imputation.

Recent advances in spatial transcriptomics allow spatially resolved gene expression measurements with cellular or even sub-cellular resolution, directly characterizing the complex spatiotemporal gene expression landscape and cell-to-cell interactions in their native microenvironments. Due to technology limitations, most spatial transcriptomic technologies still yield incomplete expression measurements with excessive missing values. Therefore, gene imputation is critical to filling in missing data, enhancing resolution, and improving overall interpretability. However, existing methods either require additional matched single-cell RNA-seq data, which is rarely available, or ignore spatial proximity or expression similarity information. To address these issues, we introduce Impeller, a path-based heterogeneous graph learning method for spatial transcriptomic data imputation. Impeller has two unique characteristics distinct from existing approaches. First, it builds a heterogeneous graph with two types of edges representing spatial proximity and expression similarity. Therefore, Impeller can simultaneously model smooth gene expression changes across spatial dimensions and capture similar gene expression signatures of faraway cells from the same type. Moreover, Impeller incorporates both short- and long-range cell-to-cell interactions (e.g. via paracrine and endocrine) by stacking multiple GNN layers. We use a learnable path operator in Impeller to avoid the over-smoothing issue of the traditional Laplacian matrices. Extensive experiments on diverse datasets from three popular platforms and two species demonstrate the superiority of Impeller over various state-of-the-art imputation methods. The code and preprocessed data used in this study are available at https://github.com/aicb-ZhangLabs/Impeller and https://zenodo.org/records/11212604.

Anterior-posterior Wnt signaling network conservation between indirect developing sea urchin and hemichordate embryos.

How animal body plans evolved and diversified is a major question in evolutionary developmental biology. To address this question, it is important to characterize the exact molecular mechanisms that establish the major embryonic axes which give rise to the adult animal body plan. The anterior-posterior (AP) axis is the first axis to be established in most animal embryos, and in echinoderm sea urchin embryos its formation is governed by an integrated network of three different Wnt signaling pathways: Wnt/β-catenin, Wnt/JNK, and Wnt/PKC pathway. The extent to which this embryonic patterning mechanism is conserved among deuterostomes, or more broadly in metazoans, is an important open question whose answers could lead to a deeper appreciation of the evolution of the AP axis. Because Ambulacrarians (echinoderms and hemichordates) reside in a key phylogenetic position as the sister group to chordates, studies in these animals can help inform on how chordate body plans may have evolved. Here, we assayed the spatiotemporal gene expression of a subset of sea urchin AP Wnt patterning gene orthologs in the hemichordate, Schizocardium californicum. Our results show that positioning of the anterior neuroectoderm (ANE) to a territory around the anterior pole during early AP formation is spatially and temporally similar between indirect developing hemichordates and sea urchins. Furthermore, we show that the expression of wnt8 and frizzled5/8, two known drivers of ANE patterning in sea urchins, is similar in hemichordate embryos. Lastly, our results highlight divergence in embryonic expression of several early expressed Wnt genes (wnt1, wnt2 and wnt4). These results suggest that expression of the sea urchin AP Wnt signaling network is largely conserved in indirect developing hemichordates setting the foundation for future functional studies in S. californicum.

MSL2 variants lead to a neurodevelopmental syndrome with lack of coordination, epilepsy, specific dysmorphisms, and a distinct episignature

Epigenetic dysregulation has emerged as an important etiological mechanism of neurodevelopmental disorders (NDDs). Pathogenic variation in epigenetic regulators can impair deposition of histone post-translational modifications leading to aberrant spatiotemporal gene expression during neurodevelopment. The male-specific lethal (MSL) complex is a prominent multi-subunit epigenetic regulator of gene expression and is responsible for histone 4 lysine 16 acetylation (H4K16ac). Using exome sequencing, here we identify a cohort of 25 individuals with heterozygous de novo variants in MSL complex member MSL2. MSL2 variants were associated with NDD phenotypes including global developmental delay, intellectual disability, hypotonia, and motor issues such as coordination problems, feeding difficulties, and gait disturbance. Dysmorphisms and behavioral and/or psychiatric conditions, including autism spectrum disorder, and to a lesser extent, seizures, connective tissue disease signs, sleep disturbance, vision problems, and other organ anomalies, were observed in affected individuals. As a molecular biomarker, a sensitive and specific DNA methylation episignature has been established. Induced pluripotent stem cells (iPSCs) derived from three members of our cohort exhibited reduced MSL2 levels. Remarkably, while NDD-associated variants in two other members of the MSL complex (MOF and MSL3) result in reduced H4K16ac, global H4K16ac levels are unchanged in iPSCs with MSL2 variants. Regardless, MSL2 variants altered the expression of MSL2 targets in iPSCs and upon their differentiation to early germ layers. Our study defines an MSL2-related disorder as an NDD with distinguishable clinical features, a specific blood DNA episignature, and a distinct, MSL2-specific molecular etiology compared to other MSL complex-related disorders.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Lineage Analysis and Birthdating of Drosophila Central Complex Lineages.

From insects to humans, the nervous system generates complex behaviors mediated by distinct neural circuits that are composed of diverse cell types. During development, the spatiotemporal gene expression of the neural progenitors expands the diversity of neuronal and glial subtypes. Various neural stem cell-intrinsic and -extrinsic gene programs have been identified that are thought to play a major role in generating diverse neuronal and glial cell types. Drosophila has served as an excellent model system for discovering the fundamental principles of nervous system development and function. The sophisticated genetic tools allow us to link the origin and birth timing (the time when a particular neuron is born during development) of neuron types to unique neural stem cells (NSCs) and to a developmental time. In Drosophila, a special class of NSCs called Type II NSCs has adopted a more advanced division mode to generate lineages for the higher-order brain center, the central complex, which is an evolutionarily conserved brain region found in all insects. Type II NSCs, similar to the human outer radial glia, generate intermediate neural progenitors (INPs), which divide many times to produce about eight to 10 progeny. Both Type II NSCs and INPs express distinct transcription factors and RNA-binding proteins that have been proposed to regulate the specification of cell types populating the adult central complex. Here, we describe the recently invented lineage filtering system, called cell class-lineage intersection (CLIn), which enables the tracking and birthdating of the Type II NSC lineages. Using CLIn, one can easily generate clones of different Type II NSCs and identify not only the origins of neurons of interest but also their birth time.

272 Integrative Genomics Implicates Primary Impact on Cortical Development in Spina Bifida

INTRODUCTION: Spina bifida is one of the most common malformation of the central nervous system that is associated with poor neurodevelopmental outcomes. Altered neurodevelopment in spina bifida has been attributed to venting of CSF through the neural tube defect. Despite correction of the spinal defect, some patients still develop hydrocephalus with poor neurodevelopmental outcomes, suggesting incomplete understanding of disease pathogenesis. Recent human genetic efforts have identified candidate genes that confer spina bifida risk, however, the mechanisms linking gene mutations to altered neurodevelopment remain poorly understood. METHODS: We curated a list of 439 genes that confer risk to spina bifida from a whole-genome sequencing study of 310 individuals. We then mapped the spatiotemporal and cell-type specific expression of these risk genes onto bulk and single-cell RNA sequencing data of normally developing human brains covering the entire lifespan of human development. RESULTS: We found that spina bifida risk genes are enriched at post-conception weeks 5-6 and 25-38 in the developing human brain. After birth, spina bifida risk genes are enriched at postnatal months 6-10 and then after 60 years of age. At the cell type level, spina bifida risk genes are enriched in proliferative neuroprogenitor cells and microglia in the prenatal human cortex. CONCLUSIONS: Our analyses show that spina bifida risk genes are expressed in the developing human brain, suggesting a primary impact on brain development that is independent from the open caudal end of the neural tube defect. The finding that spina bifida risk genes continue to be expressed postnatally suggests a window for therapeutic intervention to optimize neurodevelopmental outcomes after birth in at least a subset of patients.

Genome-wide identification and functional characterization of trans-IDS genes give new insights into terpene biosynthesis in Syringa oblata

Floral scents play a vital role in plant pollination and ornamental plant industry. Syringa oblata, a typical aromatic and medicinal plant, holds a high terpene concentration as the main fragrant volatiles and plays a crucial role in the breeding in the Syringa genus. The trans-isopentenyl diphosphate synthases (trans-IDSs) were critical enzymes in terpenoid diversity and accumulation. However, genome-wide characterization and systematic analysis of the lilac trans-IDS family are still under-explored. In this study, we identified sixteen trans-IDS genes belonging to five subfamilies on 10 chromosomes of the lilac genome. Each subfamily shared similar exon, intron, and motif distributions. Synteny analysis suggested that the lilac trans-IDS gene family expanded through segmental and tandem duplication events. We also analyzed their expression levels in various organs and flower development stages for terpene biosynthesis functional insights. Furthermore, spatiotemporal gene expression, virus-induced gene silencing (VIGS), transient overexpression, and transgenic Arabidopsis results collectively showed SoGPPS2 enhances terpene volatile synthesis. SoGPPS2 also impacted plant growth and flowering. And potential upstream transcription regulatory networks for SoGPPS2 were identified. In summary, this study enhances our comprehensive understanding of the lilac trans-IDS gene family for the first time, providing crucial insights into terpene synthesis and laying the groundwork for future floral aromatic breeding, medicinal and essential oil industries in lilac.

Metabolic pathways of glucose and fructose: II spatiotemporal expression of genes involved in synthesis and transport of lactate in ovine conceptuses.

Lactate, an abundant molecule in fetal fluids and blood of mammalian species is often overlooked as a metabolic waste product generated during pregnancy. Most of the glucose and fructose consumed by ovine conceptuses is converted to lactate, but proteins involved in lactate metabolism and transport have not been investigated. This study characterized total lactate produced by ovine conceptuses throughout gestation, as well as expression of mRNAs and proteins involved in lactate metabolism. Lactate increased in abundance in the uterine lumen during the preimplantation period and was more abundant than pyruvate. The abundance of lactate in allantoic and amniotic fluids increased with advancing days of gestation and most abundant on Day 125 of pregnancy (P < 0.05). Lactate dehydrogenase (LDH) subunits A (converts pyruvate to lactate) and B (converts lactate to pyruvate) were expressed by conceptuses throughout gestation. Lactate is transported via monocarboxylic acid transporters SLC16A1 and SLC16A3, both of which were expressed by the conceptus throughout gestation. Additionally, the interplacentomal chorioallantois from Day 126 expressed SLC16A1 and SLC16A3 and transported lactate across the tissue. Hydrocarboxylic acid receptor 1 (HCAR1), a receptor for lactate, was localized to the uterine luminal and superficial glandular epithelia of pregnant ewes throughout gestation, and conceptus trophectoderm during the peri-implantation period of gestation. These results provide novel insights into the spatiotemporal profiles of enzymes, transporters, and receptor for lactate by ovine conceptuses throughout pregnancy.

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Precise spatio-temporal expression of genes in organogenesis is regulated by the coordinated interplay of DNA elements such as promoter and enhancers present in the regulatory region of a given locus. POU1F1 transcription factor plays a crucial role in the development of somatotrophs, lactotrophs and thyrotrophs in the anterior pituitary gland, and in maintaining high expression of growth hormone, prolactin and TSH. In mouse, expression of POU1F1 is controlled by a region fenced by two CTCF sites, containing 5 upstream enhancer elements, designated E-A (5′ to 3′). Elements C, B and A correspond to elements shown previously to play a role in pituitary development and hormonal expression; functional roles for elements E and D have not been reported. We performed comparative sequence analysis of this regulatory region and discovered that three elements, B, C and E, are present in all vertebrate groups except Agnatha. One very long (>2 kb) element (A) is unique to mammals suggesting a specific change in regulation of the gene in this group. Using DNA accessibility assay (ATAC-seq) we showed that conserved elements in anterior pituitary of four non-mammals are open, suggesting functionality as regulatory elements. We showed that, in many non-mammalian vertebrates, an additional upstream exon closely follows element E, leading to alternatively spliced transcripts. Here, element E functions as an alternative promoter, but in mammals this feature is lost, suggesting conversion of alternative promoter to enhancer. Our work shows that regulation of POU1F1 changed markedly during the course of vertebrate evolution, use of a low number of enhancer elements combined with alternative promoters in non-mammalian vertebrates being replaced by use of a unique combination of regulatory units in mammals. Most importantly, our work suggests that evolutionary conversion of alternate promoter to enhancer could be one of the evolutionary mechanisms of enhancer birth.

Spatiotemporal formation of glands in plants is modulated by MYB-like transcription factors

About one third of vascular plants develop glandular trichomes, which produce defensive compounds that repel herbivores and act as a natural biofactory for important pharmaceuticals such as artemisinin and cannabinoids. However, only a few regulators of glandular structures have been characterized so far. Here we have identified two closely-related MYB-like genes that redundantly inhibit the formation of glandular cells in tomatoes, and they are named as GLAND CELL REPRESSOR (GCR) 1 and 2. The GCR genes highly express in the apical cells of tomato trichomes, with expression gradually diminishing as the cells transition into glands. The spatiotemporal expression of GCR genes is coordinated by a two-step inhibition process mediated by SlTOE1B and GCRs. Furthermore, we demonstrate that the GCR genes act by suppressing Leafless (LFS), a gene that promotes gland formation. Intriguingly, homologous GCR genes from tobacco and petunia also inhibit gland formation, suggesting that the GCR-mediated repression mechanism likely represents a conserved regulatory pathway for glands across different plant species.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

The spectrum of pre-mRNA splicing in autism.

Disruptions in spatiotemporal gene expression can result in atypical brain function. Specifically, autism spectrum disorder (ASD) is characterized by abnormalities in pre-mRNA splicing. Abnormal splicing patterns have been identified in the brains of individuals with ASD, and mutations in splicing factors have been found to contribute to neurodevelopmental delays associated with ASD. Here we review studies that shed light on the importance of splicing observed in ASD and that explored the intricate relationship between splicing factors and ASD, revealing how disruptions in pre-mRNA splicing may underlie ASD pathogenesis. We provide an overview of the research regarding all splicing factors associated with ASD and place a special emphasis on five specific splicing factors-HNRNPH2, NOVA2, WBP4, SRRM2, and RBFOX1-known to impact the splicing of ASD-related genes. In the discussion of the molecular mechanisms influenced by these splicing factors, we lay the groundwork for a deeper understanding of ASD's complex etiology. Finally, we discuss the potential benefit of unraveling the connection between splicing and ASD for the development of more precise diagnostic tools and targeted therapeutic interventions. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Evolution and Genomics > Computational Analyses of RNA RNA-Based Catalysis > RNA Catalysis in Splicing and Translation.

Self-organized BMP signaling dynamics underlie the development and evolution of digit segmentation patterns in birds and mammals.

Repeating patterns of synovial joints are a highly conserved feature of articulated digits, with variations in joint number and location resulting in diverse digit morphologies and limb functions across the tetrapod clade. During the development of the amniote limb, joints form iteratively within the growing digit ray, as a population of distal progenitors alternately specifies joint and phalanx cell fates to segment the digit into distinct elements. While numerous molecular pathways have been implicated in this fate choice, it remains unclear how they give rise to a repeating pattern. Here, using single-cell RNA sequencing and spatial gene expression profiling, we investigate the transcriptional dynamics of interphalangeal joint specification in vivo. Combined with mathematical modeling, we predict that interactions within the BMP signaling pathway-between the ligand GDF5, the inhibitor NOGGIN, and the intracellular effector pSMAD-result in a self-organizing Turing system that forms periodic joint patterns. Our model is able to recapitulate the spatiotemporal gene expression dynamics observed in vivo, as well as phenocopy digit malformations caused by BMP pathway perturbations. By contrasting in silico simulations with in vivo morphometrics of two morphologically distinct digits, we show how changes in signaling parameters and growth dynamics can result in variations in the size and number of phalanges. Together, our results reveal a self-organizing mechanism that underpins amniote digit segmentation and its evolvability and, more broadly, illustrate how Turing systems based on a single molecular pathway may generate complex repetitive patterns in a wide variety of organisms.

Three foxg1 paralogues in lampreys and gnathostomes-brothers or cousins?

Foxg1 is a key regulator of the early development of the vertebrate forebrain and sensory organs. In this study, we describe for the first time three foxg1 paralogues in lamprey, representative of one of two basally diverged lineages of vertebrates-the agnathans. We also first describe three foxg1 genes in sterlet-representative of one of the evolutionarily ancient clades of gnathostomes. According to the analysis of local genomic synteny, three foxg1 genes of agnathans and gnathostomes have a common origin as a result of two rounds of genomic duplications in the early evolution of vertebrates. At the same time, it is difficult to reliably establish pairwise orthology between foxg1 genes of agnathans and gnathostomes based on the analysis of phylogeny and local genomic synteny, as well as our studies of the spatiotemporal expression of foxg1 genes in the river lamprey Lampetra fluviatilis and the sterlet Acipenser ruthenus. Thus, the appearance of three foxg1 paralogues in agnathans and gnathostomes could have occurred either as a result of two rounds of duplication of the vertebrate common ancestor genome (2R hypothesis) or as a result of the first common round followed by subsequent independent polyploidizations in two evolutionary lineages (1R hypothesis).

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

Investigating the spatiotemporal expression of CBTS genes lead to the discovery of tobacco root as a cembranoid-producing organ.

Tobacco cembranoids, known for their anti-insect and antifungal properties, were shown to be mainly present on the surface of leaves and flowers, being biosynthesized by their trichomes. It remains unclear whether they could be biosynthesized in other organs without trichomes. Cembratrien-ol synthases (CBTSs) catalyze the conversion of GGPP to CBT-ols and thus play an important role in cembranoid biosynthesis. This study identified the CBTS family genes in tobacco and examined their spatiotemporal expression patterns. The CBTS genes showed diverse expression patterns in tobacco organs, with the majority highly expressed in leaves and a few highly expressed in flowers. The expression of CBTS genes were also correlated with the development of tobacco plants, and most of them showed the highest expression level at the budding stage. Furthermore, their expression is mediated by the JA (jasmonate) signaling in all tobacco organs. Several CBTS genes were found to be highly expressed in tobacco roots that have no trichomes, which prompted us to determine the cembranoid production in roots and other organs. GC-MS and UPLC assays revealed that cembranoids were produced in all tobacco organs, which was supported by the bioactivity assay results that almost all these CBTS enzymes could catalyze CBT-ol biosyntheis in yeast, and that the content ratio of CBT-ols and CBT-diols in tobacco roots was different to that in leaves. This work sheds insights into the expression profiles of tobacco CBTS genes and provides a feasibility to engineer tobacco roots for industrial production of cembranoids.

3D organization of enhancers in MuSCs.

Skeletal muscle stem cells (MuSCs), also known as satellite cells, are essential for muscle growth and injury induced regeneration. In healthy adult muscle, MuSCs remain in a quiescent state located in a specialized niche beneath the basal lamina. Upon injury, these dormant MuSCs can quickly activate to re-enter the cell cycle and differentiate into new myofibers, while a subset undergoes self-renewal and returns to quiescence to restore the stem cell pool. The myogenic lineage progression is intricately controlled by complex intrinsic and extrinsic cues and coupled with dynamic transcriptional programs. In transcriptional regulation, enhancers are key regulatory elements controlling spatiotemporal gene expression through physical contacting promoters of target genes. The three-dimensional (3D) chromatin architecture is known to orchestrate the establishment of proper enhancer-promoterinteractions throughout development and aging. However, studies dissecting the 3D organization of enhancers in MuSCs are just emerging. Here, we provide an overview of the general properties of enhancers and newly developed methods for assessing their activity. In particular, we summarize recent discoveries regarding the 3D rewiring of enhancers during MuSC specification, lineage progression as well as aging.

Spatiotemporal expression and coexpression patterns of SRPK1 in the human brain: A neurodevelopmental perspective.

SRPK1 is a splicing-related protein that plays an important role in the development and function of the human brain. This article presents evidence that SRPK1 has distinct spatiotemporal expression patterns enriched in processes related to neurodevelopmental disorders across development. We used the BrainSpan growing mammalian brain transcriptome to evaluate the distribution of SRPK1 throughout the entire brain. RNA-sequencing data were gathered from 524 tissue samples recovered from 41 postmortem brains of physiologically normal individuals spanning early developing fetus (8 postconception weeks, PCW) to later life (40 years of age). Using the Allen Human Brain Atlas (AHBA) dataset, we analyzed the spatial gene expression of 15 adult human brains. Using Toppgene, we identified genes that exhibit significant coexpression with SRPK1. We found evidence that analyzing the spatiotemporal gene expression profile and identifying coexpressed genes reveals that SRPK1 expression is involved in various neurodevelopmental and somatic events throughout the lifetime. Our findings highlight the importance of detailed maps of gene expression in the human brain for improved human-to-human translation and illustrate differences in SRPK1 expression across anatomical areas and developmental stages in healthy human brain tissue.

Aberrant epithelial cell interaction promotes esophageal squamous-cell carcinoma development and progression

Epithelial-mesenchymal transition (EMT) and proliferation play important roles in epithelial cancer formation and progression, but what molecules and how they trigger EMT is largely unknown. Here we performed spatial transcriptomic and functional analyses on samples of multistage esophageal squamous-cell carcinoma (ESCC) from mice and humans to decipher these critical issues. By investigating spatiotemporal gene expression patterns and cell–cell interactions, we demonstrated that the aberrant epithelial cell interaction via EFNB1-EPHB4 triggers EMT and cell cycle mediated by downstream SRC/ERK/AKT signaling. The aberrant epithelial cell interaction occurs within the basal layer at early precancerous lesions, which expands to the whole epithelial layer and strengthens along the cancer development and progression. Functional analysis revealed that the aberrant EFNB1-EPHB4 interaction is caused by overexpressed ΔNP63 due to TP53 mutation, the culprit in human ESCC tumorigenesis. Our results shed new light on the role of TP53-TP63/ΔNP63-EFNB1-EPHB4 axis in EMT and cell proliferation in epithelial cancer formation.