Spatial transcriptomics reveals novel genes during the remodelling of the embryonic human arterial valves.

Abnormalities of the arterial valves, including bicuspid aortic valve (BAV) are amongst the most common congenital defects and are a significant cause of morbidity as well as predisposition to disease in later life. Despite this, and compounded by their small size and relative inaccessibility, there is still much to understand about how the arterial valves form and remodel during embryogenesis, both at the morphological and genetic level. Here we set out to address this in human embryos, using Spatial Transcriptomics (ST). We show that ST can be used to investigate the transcriptome of the developing arterial valves, circumventing the problems of accurately dissecting out these tiny structures from the developing embryo. We show that the transcriptome of CS16 and CS19 arterial valves overlap considerably, despite being several days apart in terms of human gestation, and that expression data confirm that the great majority of the most differentially expressed genes are valve-specific. Moreover, we show that the transcriptome of the human arterial valves overlaps with that of mouse atrioventricular valves from a range of gestations, validating our dataset but also highlighting novel genes, including four that are not found in the mouse genome and have not previously been linked to valve development. Importantly, our data suggests that valve transcriptomes are under-represented when using commonly used databases to filter for genes important in cardiac development; this means that causative variants in valve-related genes may be excluded during filtering for genomic data analyses for, for example, BAV. Finally, we highlight "novel" pathways that likely play important roles in arterial valve development, showing that mouse knockouts of RBP1 have arterial valve defects. Thus, this study has confirmed the utility of ST for studies of the developing heart valves and broadens our knowledge of the genes and signalling pathways important in human valve development.

Unified rhombic lip origins of group 3 and group 4 medulloblastoma.

Medulloblastoma, a malignant childhood cerebellar tumour, segregates molecularly into biologically distinct subgroups, suggesting that a personalized approach to therapy would be beneficial1. Mouse modelling and cross-species genomics have provided increasing evidence of discrete, subgroup-specific developmental origins2. However, the anatomical and cellular complexity of developing human tissues3-particularly within the rhombic lip germinal zone, which produces all glutamatergic neuronal lineages before internalization into the cerebellar nodulus-makes it difficult to validate previous inferences that were derived from studies in mice. Here we use multi-omics to resolve the origins of medulloblastoma subgroups in the developing human cerebellum. Molecular signatures encoded within a human rhombic-lip-derived lineage trajectory aligned with photoreceptor and unipolar brush cell expression profiles that are maintained in group 3 and group 4 medulloblastoma, suggesting a convergent basis. A systematic diagnostic-imaging review of a prospective institutional cohort localized the putative anatomical origins of group 3 and group 4 tumours to the nodulus. Our results connect the molecular and phenotypic features of clinically challenging medulloblastoma subgroups to their unified beginnings in the rhombic lip in the early stages of human development.

Spatial and cell type transcriptional landscape of human cerebellar development.

The human neonatal cerebellum is one-fourth of its adult size yet contains the blueprint required to integrate environmental cues with developing motor, cognitive and emotional skills into adulthood. Although mature cerebellar neuroanatomy is well studied, understanding of its developmental origins is limited. In this study, we systematically mapped the molecular, cellular and spatial composition of human fetal cerebellum by combining laser capture microscopy and SPLiT-seq single-nucleus transcriptomics. We profiled functionally distinct regions and gene expression dynamics within cell types and across development. The resulting cell atlas demonstrates that the molecular organization of the cerebellar anlage recapitulates cytoarchitecturally distinct regions and developmentally transient cell types that are distinct from the mouse cerebellum. By mapping genes dominant for pediatric and adult neurological disorders onto our dataset, we identify relevant cell types underlying disease mechanisms. These data provide a resource for probing the cellular basis of human cerebellar development and disease.

Spatiotemporal expansion of primary progenitor zones in the developing human cerebellum.

We present histological and molecular analyses of the developing human cerebellum from 30 days after conception to 9 months after birth. Differences in developmental patterns between humans and mice include spatiotemporal expansion of both ventricular and rhombic lip primary progenitor zones to include subventricular zones containing basal progenitors. The human rhombic lip persists longer through cerebellar development than in the mouse and undergoes morphological changes to form a progenitor pool in the posterior lobule, which is not seen in other organisms, not even in the nonhuman primate the macaque. Disruptions in human rhombic lip development are associated with posterior cerebellar vermis hypoplasia and Dandy-Walker malformation. The presence of these species-specific neural progenitor populations refines our insight into human cerebellar developmental disorders.

Recessive gene disruptions in autism spectrum disorder.

Autism spectrum disorder (ASD) affects up to 1 in 59 individuals1. Genome-wide association and large-scale sequencing studies strongly implicate both common variants2-4 and rare de novo variants5-10 in ASD. Recessive mutations have also been implicated11-14 but their contribution remains less well defined. Here we demonstrate an excess of biallelic loss-of-function and damaging missense mutations in a large ASD cohort, corresponding to approximately 5% of total cases, including 10% of females, consistent with a female protective effect. We document biallelic disruption of known or emerging recessive neurodevelopmental genes (CA2, DDHD1, NSUN2, PAH, RARB, ROGDI, SLC1A1, USH2A) as well as other genes not previously implicated in ASD including FEV (FEV transcription factor, ETS family member), which encodes a key regulator of the serotonergic circuitry. Our data refine estimates of the contribution of recessive mutation to ASD and suggest new paths for illuminating previously unknown biological pathways responsible for this condition.

Evolutionarily dynamic alternative splicing of GPR56 regulates regional cerebral cortical patterning.

The human neocortex has numerous specialized functional areas whose formation is poorly understood. Here, we describe a 15-base pair deletion mutation in a regulatory element of GPR56 that selectively disrupts human cortex surrounding the Sylvian fissure bilaterally including "Broca's area," the primary language area, by disrupting regional GPR56 expression and blocking RFX transcription factor binding. GPR56 encodes a heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptor required for normal cortical development and is expressed in cortical progenitor cells. GPR56 expression levels regulate progenitor proliferation. GPR56 splice forms are highly variable between mice and humans, and the regulatory element of gyrencephalic mammals directs restricted lateral cortical expression. Our data reveal a mechanism by which control of GPR56 expression pattern by multiple alternative promoters can influence stem cell proliferation, gyral patterning, and, potentially, neocortex evolution.

Investigating embryonic expression patterns and evolution of AHI1 and CEP290 genes, implicated in Joubert syndrome.

Joubert syndrome and related diseases (JSRD) are developmental cerebello-oculo-renal syndromes with phenotypes including cerebellar hypoplasia, retinal dystrophy and nephronophthisis (a cystic kidney disease). We have utilised the MRC-Wellcome Trust Human Developmental Biology Resource (HDBR), to perform in-situ hybridisation studies on embryonic tissues, revealing an early onset neuronal, retinal and renal expression pattern for AHI1. An almost identical pattern of expression is seen with CEP290 in human embryonic and fetal tissue. A novel finding is that both AHI1 and CEP290 demonstrate strong expression within the developing choroid plexus, a ciliated structure important for central nervous system development. To test if AHI1 and CEP290 may have co-evolved, we carried out a genomic survey of a large group of organisms across eukaryotic evolution. We found that, in animals, ahi1 and cep290 are almost always found together; however in other organisms either one may be found independent of the other. Finally, we tested in murine epithelial cells if Ahi1 was required for recruitment of Cep290 to the centrosome. We found no obvious differences in Cep290 localisation in the presence or absence of Ahi1, suggesting that, while Ahi1 and Cep290 may function together in the whole organism, they are not interdependent for localisation within a single cell. Taken together these data support a role for AHI1 and CEP290 in multiple organs throughout development and we suggest that this accounts for the wide phenotypic spectrum of AHI1 and CEP290 mutations in man.

Spatiotemporal expression pattern of KIF21A during normal embryonic development and in congenital fibrosis of the extraocular muscles type 1 (CFEOM1).

Congenital fibrosis of the extraocular muscles type 1 (CFEOM1) is a rare inherited strabismus syndrome characterized by non-progressive ophthalmoplegia. We previously identified that CFEOM1 results from heterozygous missense mutations in KIF21A, which encodes a kinesin motor protein. Here we evaluate the expression pattern of KIF21A in human brain and muscles of control and CFEOM1 patients, and during human and mouse embryonic development. KIF21A is expressed in the cell bodies, axons, and dendrites of many neuronal populations including those in the hippocampus, cerebral cortex, cerebellum, striatum, and motor neurons of the oculomotor, trochlear, and abducens nuclei from early development into maturity, and its spatial distribution is not altered in the CFEOM1 tissues available for study. Multiple splice isoforms of KIF21A are identified in human fetal brain, but none of the reported CFEOM1 mutations are located in or near the alternatively spliced exons. KIF21A immunoreactivity is also observed in extraocular and skeletal muscle biopsies of control and CFEOM1 patients, where it co-localizes with triadin, a marker of the excitation-contractile coupling system. The diffuse and widespread expression of KIF21A in the developing human and mouse central and peripheral nervous system as well as in extraocular muscle does not account for the restricted ocular phenotype observed in CFEOM1, nor does it permit the formal exclusion of a myogenic etiology based on expression patterns alone.

Immunophenotyping of diffuse large B-cell lymphoma (DLBCL) defines multiple sub-groups of germinal centre-like tumours displaying different survival characteristics.

Diffuse large B-cell lymphoma (DLBCL) forms a heterogeneous collection of aggressive non-Hodgkin's Lymphoma in which three principle classes of neoplasia have been defined according to gene expression and immunophenotyping studies. The present investigation sought to examine the immunophenotype of proposed subgroups and relate these to patient survival. A series of 155 DLBCL treated uniformly with anthracycline therapy in clinical trials, were stratified upon the basis of common biomarker expression with combination immunophenotype being related to patient overall survival. Stratification of tumours with respect to combined expression profiles of the three biological markers (CD10, Bcl-6 and MUM-1) revealed six groups showing significant differences in survival (p=0.014). The greatest difference resided between distinct populations of germinal centre (GC) cell tumours; the first being CD10-, Bcl-6+, MUM-1- and the second CD10+ Bcl-6+ MUM-1+ (p=0.002). The former group displayed median survival time of 143 months, the latter only 11 months. A third population of GC tumours (CD10+ Bcl-6+ and MUM-1-) also displayed a relative short median survival (32 months). Of the three groups presenting a non-GC or activated B cell (NGC/ABC) phenotype, only one (CD10-, Bcl-6+ and MUM-1+) presented short-term median survival (27 months) comparable with poor prognosis GC sub-populations. Within the remaining ABC tumour groups (CD10- Bcl-6- MUM-1- and CD10- Bcl-6- MUM-1+) patients presented intermediate median survival times of 54 and 58 months, respectively. Thus, the GC phenotype did not act as a universal indicator of good clinical prognosis, but rather multiple groups of GC tumours were associated with distinct overall survival profiles. Ultimately, the data allowed definition of a predictive algorithm defining three groups predicting poor, intermediate and good clinical prognosis. The first of these comprised two patient sub-populations with GC-like tumours together with one sub-population of NGC/ABC, the second two sub-populations of ABC-like tumours, and the final a single group of GC-like tumours associated with optimal long-term survival.