Distinct transcriptome architectures underlying lupus establishment and exacerbation.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease involving multiple immune cells. To elucidate SLE pathogenesis, it is essential to understand the dysregulated gene expression pattern linked to various clinical statuses with a high cellular resolution. Here, we conducted a large-scale transcriptome study with 6,386 RNA sequencing data covering 27 immune cell types from 136 SLE and 89 healthy donors. We profiled two distinct cell-type-specific transcriptomic signatures: disease-state and disease-activity signatures, reflecting disease establishment and exacerbation, respectively. We then identified candidate biological processes unique to each signature. This study suggested the clinical value of disease-activity signatures, which were associated with organ involvement and therapeutic responses. However, disease-activity signatures were less enriched around SLE risk variants than disease-state signatures, suggesting that current genetic studies may not well capture clinically vital biology. Together, we identified comprehensive gene signatures of SLE, which will provide essential foundations for future genomic and genetic studies.

Whole-body integration of gene expression and single-cell morphology.

Animal bodies are composed of cell types with unique expression programs that implement their distinct locations, shapes, structures, and functions. Based on these properties, cell types assemble into specific tissues and organs. To systematically explore the link between cell-type-specific gene expression and morphology, we registered an expression atlas to a whole-body electron microscopy volume of the nereid Platynereis dumerilii. Automated segmentation of cells and nuclei identifies major cell classes and establishes a link between gene activation, chromatin topography, and nuclear size. Clustering of segmented cells according to gene expression reveals spatially coherent tissues. In the brain, genetically defined groups of neurons match ganglionic nuclei with coherent projections. Besides interneurons, we uncover sensory-neurosecretory cells in the nereid mushroom bodies, which thus qualify as sensory organs. They furthermore resemble the vertebrate telencephalon by molecular anatomy. We provide an integrated browser as a Fiji plugin for remote exploration of all available multimodal datasets.

Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases.

Genetic studies have revealed many variant loci that are associated with immune-mediated diseases. To elucidate the disease pathogenesis, it is essential to understand the function of these variants, especially under disease-associated conditions. Here, we performed a large-scale immune cell gene-expression analysis, together with whole-genome sequence analysis. Our dataset consists of 28 distinct immune cell subsets from 337 patients diagnosed with 10 categories of immune-mediated diseases and 79 healthy volunteers. Our dataset captured distinctive gene-expression profiles across immune cell types and diseases. Expression quantitative trait loci (eQTL) analysis revealed dynamic variations of eQTL effects in the context of immunological conditions, as well as cell types. These cell-type-specific and context-dependent eQTLs showed significant enrichment in immune disease-associated genetic variants, and they implicated the disease-relevant cell types, genes, and environment. This atlas deepens our understanding of the immunogenetic functions of disease-associated variants under in vivo disease conditions.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms.

BACKGROUND: Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots. RESULTS: In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity. CONCLUSIONS: These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

Widespread diversity in the transcriptomes of functionally divergent limb tendons.

KEY POINTS: Tendon is a hypocellular, matrix-rich tissue that has been excluded from comparative transcriptional atlases. These atlases have provided important knowledge about biological heterogeneity between tissues, and our study addresses this important gap. We performed measures on four of the most studied tendons, the Achilles, forepaw flexor, patellar and supraspinatus tendons of both mice and rats. These tendons are functionally distinct and are also among the most commonly injured, and therefore of important translational interest. Approximately one-third of the filtered transcriptome was differentially regulated between Achilles, forepaw flexor, patellar and supraspinatus tendons within either mice or rats. Nearly two-thirds of the transcripts that are expressed in anatomically similar tendons were different between mice and rats. The overall findings from this study identified that although tendons across the body share a common anatomical definition based on their physical location between skeletal muscle and bone, tendon is a surprisingly genetically heterogeneous tissue. ABSTRACT: Tendon is a functionally important connective tissue that transmits force between skeletal muscle and bone. Previous studies have evaluated the architectural designs and mechanical properties of different tendons throughout the body. However, less is known about the underlying transcriptional differences between tendons that may dictate their designs and properties. Therefore, our objective was to develop a comprehensive atlas of the transcriptome of limb tendons in adult mice and rats using systems biology techniques. We selected the Achilles, forepaw digit flexor, patellar, and supraspinatus tendons due to their divergent functions and high rates of injury and tendinopathies in patients. Using RNA sequencing data, we generated the Comparative Tendon Transcriptional Database (CTTDb) that identified substantial diversity in the transcriptomes of tendons both within and across species. Approximately 30% of filtered transcripts were differentially regulated between tendons of a given species, and nearly 60% of the filtered transcripts present in anatomically similar tendons were different between species. Many of the genes that differed between tendons and across species are important in tissue specification and limb morphogenesis, tendon cell biology and tenogenesis, growth factor signalling, and production and maintenance of the extracellular matrix. This study indicates that tendon is a surprisingly heterogenous tissue with substantial genetic variation based on anatomical location and species.

Arachis hypogaea gene expression atlas for fastigiata subspecies of cultivated groundnut to accelerate functional and translational genomics applications.

Spatio-temporal and developmental stage-specific transcriptome analysis plays a crucial role in systems biology-based improvement of any species. In this context, we report here the Arachis hypogaea gene expression atlas (AhGEA) for the world's widest cultivated subsp. fastigiata based on RNA-seq data using 20 diverse tissues across five key developmental stages. Approximately 480 million paired-end filtered reads were generated followed by identification of 81 901 transcripts from an early-maturing, high-yielding, drought-tolerant groundnut variety, ICGV 91114. Further, 57 344 genome-wide transcripts were identified with ≥1 FPKM across different tissues and stages. Our in-depth analysis of the global transcriptome sheds light into complex regulatory networks namely gravitropism and photomorphogenesis, seed development, allergens and oil biosynthesis in groundnut. Importantly, interesting insights into molecular basis of seed development and nodulation have immense potential for translational genomics research. We have also identified a set of stable expressing transcripts across the selected tissues, which could be utilized as internal controls in groundnut functional genomics studies. The AhGEA revealed potential transcripts associated with allergens, which upon appropriate validation could be deployed in the coming years to develop consumer-friendly groundnut varieties. Taken together, the AhGEA touches upon various important and key features of cultivated groundnut and provides a reference for further functional, comparative and translational genomics research for various economically important traits.

Whole-organism cellular gene-expression atlas reveals conserved cell types in the ventral nerve cord of Platynereis dumerilii.

The comparative study of cell types is a powerful approach toward deciphering animal evolution. To avoid selection biases, however, comparisons ideally involve all cell types present in a multicellular organism. Here, we use image registration and a newly developed "Profiling by Signal Probability Mapping" algorithm to generate a cellular resolution 3D expression atlas for an entire animal. We investigate three-segmented young worms of the marine annelid Platynereis dumerilii, with a rich diversity of differentiated cells present in relatively low number. Starting from whole-mount expression images for close to 100 neural specification and differentiation genes, our atlas identifies and molecularly characterizes 605 bilateral pairs of neurons at specific locations in the ventral nerve cord. Among these pairs, we identify sets of neurons expressing similar combinations of transcription factors, located at spatially coherent anterior-posterior, dorsal-ventral, and medial-lateral coordinates that we interpret as cell types. Comparison with motor and interneuron types in the vertebrate neural tube indicates conserved combinations, for example, of cell types cospecified by Gata1/2/3 and Tal transcription factors. These include V2b interneurons and the central spinal fluid-contacting Kolmer-Agduhr cells in the vertebrates, and several neuron types in the intermediate ventral ganglionic mass in the annelid. We propose that Kolmer-Agduhr cell-like mechanosensory neurons formed part of the mucociliary sole in protostome-deuterostome ancestors and diversified independently into several neuron types in annelid and vertebrate descendants.

Melonet-DB, a Grand RNA-Seq Gene Expression Atlas in Melon (Cucumis melo L.).

Melon (Cucumis melo L.) is an important Cucurbitaceae crop produced worldwide, exhibiting wide genetic variations and comprising both climacteric and non-climacteric fruit types. The muskmelon cultivar "'Earl's favorite Harukei-3 (Harukei-3)"' known for its sweetness and rich aroma is used for breeding of high-grade muskmelon in Japan. We conducted RNA sequencing (RNA-seq) transcriptome studies in 30 different tissues of the 'Harukei-3' melon. These included root, stems, leaves, flowers, regenerating callus and ovaries, in addition to the flesh and peel sampled at seven stages of fruit development. The expression patterns of 20,752 genes were determined with fragments per kilobase of transcript per million fragments sequenced (FPKM) >1 in at least one tissue. Principal component analysis distinguished 30 melon tissues based on the global gene expression profile and, further, the weighted gene correlation network analysis classified melon genes into 45 distinct coexpression groups. Some coexpression groups exhibited tissue-specific gene expression. Furthermore, we developed and published web application tools designated "'Gene expression map viewer"' and "'Coexpression viewer"' on our website Melonet-DB (http://melonet-db.agbi.tsukuba.ac.jp/) to promote functional genomics research in melon. By using both tools, we analyzed melon homologs of tomato fruit ripening regulators such as E8, RIPENING-INHIBITOR (RIN) and NON-RIPENING (NOR). The "'Coexpression viewer"' clearly distinguished fruit ripening-associated melon RIN/NOR/CNR homologs from those expressed in other tissues. In addition, several other MADS-box, NAM/ATAF/CUC (NAC) and homeobox transcription factor genes were identified as fruit ripening-associated genes. Our tools provide useful information for research not only on melon but also on other fleshy fruit plants.

A gene expression atlas of a bicoid-depleted Drosophila embryo reveals early canalization of cell fate.

In developing embryos, gene regulatory networks drive cells towards discrete terminal fates, a process called canalization. We studied the behavior of the anterior-posterior segmentation network in Drosophila melanogaster embryos by depleting a key maternal input, bicoid (bcd), and measuring gene expression patterns of the network at cellular resolution. This method results in a gene expression atlas containing the levels of mRNA or protein expression of 13 core patterning genes over six time points for every cell of the blastoderm embryo. This is the first cellular resolution dataset of a genetically perturbed Drosophila embryo that captures all cells in 3D. We describe the technical developments required to build this atlas and how the method can be employed and extended by others. We also analyze this novel dataset to characterize the degree and timing of cell fate canalization in the segmentation network. We find that in two layers of this gene regulatory network, following depletion of bcd, individual cells rapidly canalize towards normal cell fates. This result supports the hypothesis that the segmentation network directly canalizes cell fate, rather than an alternative hypothesis whereby cells are initially mis-specified and later eliminated by apoptosis. Our gene expression atlas provides a high resolution picture of a classic perturbation and will enable further computational modeling of canalization and gene regulation in this transcriptional network.

Gene expression atlas of pigeonpea and its application to gain insights into genes associated with pollen fertility implicated in seed formation.

Pigeonpea (Cajanus cajan) is an important grain legume of the semi-arid tropics, mainly used for its protein rich seeds. To link the genome sequence information with agronomic traits resulting from specific developmental processes, a Cajanus cajan gene expression atlas (CcGEA) was developed using the Asha genotype. Thirty tissues/organs representing developmental stages from germination to senescence were used to generate 590.84 million paired-end RNA-Seq data. The CcGEA revealed a compendium of 28 793 genes with differential, specific, spatio-temporal and constitutive expression during various stages of development in different tissues. As an example to demonstrate the application of the CcGEA, a network of 28 flower-related genes analysed for cis-regulatory elements and splicing variants has been identified. In addition, expression analysis of these candidate genes in male sterile and male fertile genotypes suggested their critical role in normal pollen development leading to seed formation. Gene network analysis also identified two regulatory genes, a pollen-specific SF3 and a sucrose-proton symporter, that could have implications for improvement of agronomic traits such as seed production and yield. In conclusion, the CcGEA provides a valuable resource for pigeonpea to identify candidate genes involved in specific developmental processes and to understand the well-orchestrated growth and developmental process in this resilient crop.

Full-length de novo assembly of RNA-seq data in pea (Pisum sativum L.) provides a gene expression atlas and gives insights into root nodulation in this species.

Next-generation sequencing technologies allow an almost exhaustive survey of the transcriptome, even in species with no available genome sequence. To produce a Unigene set representing most of the expressed genes of pea, 20 cDNA libraries produced from various plant tissues harvested at various developmental stages from plants grown under contrasting nitrogen conditions were sequenced. Around one billion reads and 100 Gb of sequence were de novo assembled. Following several steps of redundancy reduction, 46 099 contigs with N50 length of 1667 nt were identified. These constitute the 'Caméor' Unigene set. The high depth of sequencing allowed identification of rare transcripts and detected expression for approximately 80% of contigs in each library. The Unigene set is now available online (http://bios.dijon.inra.fr/FATAL/cgi/pscam.cgi), allowing (i) searches for pea orthologs of candidate genes based on gene sequences from other species, or based on annotation, (ii) determination of transcript expression patterns using various metrics, (iii) identification of uncharacterized genes with interesting patterns of expression, and (iv) comparison of gene ontology pathways between tissues. This resource has allowed identification of the pea orthologs of major nodulation genes characterized in recent years in model species, as a major step towards deciphering unresolved pea nodulation phenotypes. In addition to a remarkable conservation of the early transcriptome nodulation apparatus between pea and Medicago truncatula, some specific features were highlighted. The resource provides a reference for the pea exome, and will facilitate transcriptome and proteome approaches as well as SNP discovery in pea.

A gene expression atlas of early craniofacial development.

We present a gene expression atlas of early mouse craniofacial development. Laser capture microdissection (LCM) was used to isolate cells from the principal critical microregions, whose development, differentiation and signaling interactions are responsible for the construction of the mammalian face. At E8.5, as migrating neural crest cells begin to exit the neural fold/epidermal ectoderm boundary, we examined the cranial mesenchyme, composed of mixed neural crest and paraxial mesoderm cells, as well as cells from adjacent neuroepithelium. At E9.5 cells from the cranial mesenchyme, overlying olfactory placode/epidermal ectoderm, and underlying neuroepithelium, as well as the emerging mandibular and maxillary arches were sampled. At E10.5, as the facial prominences form, cells from the medial and lateral prominences, the olfactory pit, multiple discrete regions of underlying neuroepithelium, the mandibular and maxillary arches, including both their mesenchymal and ectodermal components, as well as Rathke's pouch, were similarly sampled and profiled using both microarray and RNA-seq technologies. Further, we performed single cell studies to better define the gene expression states of the early E8.5 pioneer neural crest cells and paraxial mesoderm. Taken together, and analyzable by a variety of biological network approaches, these data provide a complementing and cross validating resource capable of fueling discovery of novel compartment specific markers and signatures whose combinatorial interactions of transcription factors and growth factors/receptors are responsible for providing the master genetic blueprint for craniofacial development.

An atlas and database of neuropeptide gene expression in the adult zebrafish forebrain.

Zebrafish is a useful model organism in neuroscience; however, its gene expression atlas in the adult brain is not well developed. In the present study, we examined the expression of 38 neuropeptides, comparing with GABAergic and glutamatergic neuron marker genes in the adult zebrafish brain by comprehensive in situ hybridization. The results are summarized as an expression atlas in 19 coronal planes of the forebrain. Furthermore, the scanned data of all brain sections were made publicly available in the Adult Zebrafish Brain Gene Expression Database (https://ssbd.riken.jp/azebex/). Based on these data, we performed detailed comparative neuroanatomical analyses of the hypothalamus and found that several regions previously described as one nucleus in the reference zebrafish brain atlas contain two or more subregions with significantly different neuropeptide/neurotransmitter expression profiles. Subsequently, we compared the expression data in zebrafish telencephalon and hypothalamus obtained in this study with those in mice, by performing a cluster analysis. As a result, several nuclei in zebrafish and mice were clustered in close vicinity. The present expression atlas, database, and anatomical findings will contribute to future neuroscience research using zebrafish.

Functional analysis of the Drosophila eve locus in response to non-canonical combinations of gap gene expression levels.

Transcription factor combinations play a key role in shaping cellular identity. However, the precise relationship between specific combinations and downstream effects remains elusive. Here, we investigate this relationship within the context of the Drosophila eve locus, which is controlled by gap genes. We measure spatiotemporal levels of four gap genes in heterozygous and homozygous gap mutant embryos and correlate them with the striped eve activity pattern. Although changes in gap gene expression extend beyond the manipulated gene, the spatial patterns of Eve expression closely mirror canonical activation levels in wild type. Interestingly, some combinations deviate from the wild-type repertoire but still drive eve activation. Although in homozygous mutants some Eve stripes exhibit partial penetrance, stripes consistently emerge at reproducible positions, even with varying gap gene levels. Our findings suggest a robust molecular canalization of cell fates in gap mutants and provide insights into the regulatory constraints governing multi-enhancer gene loci.

A Multi-Tissue Gene Expression Atlas of Water Buffalo (Bubalus bubalis) Reveals Transcriptome Conservation between Buffalo and Cattle.

We generated 73 transcriptomic data of water buffalo, which were integrated with publicly available data in this species, yielding a large dataset of 355 samples representing 20 major tissue categories. We established a multi-tissue gene expression atlas of water buffalo. Furthermore, by comparing them with 4866 cattle transcriptomic data from the cattle genotype-tissue expression atlas (CattleGTEx), we found that the transcriptomes of the two species exhibited conservation in their overall gene expression patterns, tissue-specific gene expression and house-keeping gene expression. We further identified conserved and divergent expression genes between the two species, with the largest number of differentially expressed genes found in the skin, which may be related to structural and functional differences in the skin of the two species. This work provides a source of functional annotation of the buffalo genome and lays the foundations for future genetic and evolutionary studies in water buffalo.

Screening and identification of multiple abiotic stress responsive candidate genes based on hybrid-sequencing in Vicia sativa.

Common vetch is an important leguminous forage for both livestock fodder and green manure and has a tremendous latent capacity in a sustainable agroecosystem. In the present study, a comprehensive transcriptome analysis of the aboveground leaves and underground roots of common vetch under multiple abiotic stress treatments, including NaCl, drought, cold, and cold drought, was performed using hybrid-sequencing technology, i. e. single-molecule real-time sequencing technology (SMRT) and supplemented by next-generation sequencing (NGS) technology. A total of 485,038 reads of insert (ROIs) with a mean length of 2606 bp and 228,261 full-length nonchimeric (FLNC) reads were generated. After deduplication, 39,709 transcripts were generated. Of these transcripts, we identified 1059 alternative splicing (AS) events, 17,227 simple sequence repeats (SSRs), and 1647 putative transcription factors (TFs). Furthermore, 640 candidates long noncoding RNAs (lncRNAs) and 28,256 complete coding sequences (CDSs) were identified. In gene annotation analyses, a total of 38,826 transcripts (97.78%) were annotated in eight public databases. Finally, seven multiple abiotic stress-responsive candidate genes were obtained through gene expression, annotation information, and protein-protein interaction (PPI) networks. Our research not only enriched the structural information of FL transcripts in common vetch, but also provided useful information for exploring the molecular mechanism of multiple abiotic stress tolerance between aboveground and underground tissues in common vetch and related legumes.

Spatial transcriptomics map of the embryonic mouse brain - a tool to explore neurogenesis.

The developing brain has a well-organized anatomical structure comprising different types of neural and non-neural cells. Stem cells, progenitors and newborn neurons tightly interact with their neighbouring cells and tissue microenvironment, and this intricate interplay ultimately shapes the output of neurogenesis. Given the relevance of spatial cues during brain development, we acknowledge the necessity for a spatial transcriptomics map accessible to the neurodevelopmental community. To fulfil this need, we generated spatially resolved RNA sequencing (RNAseq) data from embryonic day 13.5 mouse brain sections immunostained for mitotic active neural and vascular cells. Unsupervised clustering defined specific cell type populations of diverse lineages and differentiation states. Differential expression analysis revealed unique transcriptional signatures across specific brain areas, uncovering novel features inherent to particular anatomical domains. Finally, we integrated existing single-cell RNAseq datasets into our spatial transcriptomics map, adding tissue context to single-cell RNAseq data. In summary, we provide a valuable tool that enables the exploration and discovery of unforeseen molecular players involved in neurogenesis, particularly in the crosstalk between different cell types.

Achilles Tendons Display Region-Specific Transcriptomic Signatures Associated With Distinct Mechanical Properties.

BACKGROUND: Previous studies have examined the transcriptomes and mechanical properties of whole tendons in different regions of the body. However, less is known about these characteristics within a single tendon. PURPOSE: To develop a regional transcriptomic atlas and evaluate the region-specific mechanical properties of Achilles tendons. STUDY DESIGN: Descriptive laboratory study. METHODS: Achilles tendons from 2-month-old male Sprague Dawley rats were used. Tendons were isolated and divided into proximal, middle, and distal thirds for RNA sequencing (n = 5). For mechanical testing, the Achilles muscle-tendon-calcaneus unit was mounted in a custom-designed materials testing system with the unit clamped over the musculotendinous junction (MTJ) and the calcaneus secured at 90° of dorsiflexion (n = 9). Tendons were stretched to 20 N at a constant speed of 0.0167 mm/s. Cross-sectional area, strain, stress, and Young modulus were determined in each tendon region. RESULTS: An open-access, interactive transcriptional atlas was generated that revealed distinct gene expression signatures in each tendon region. The proximal and distal regions had the largest differences in gene expression, with 2596 genes significantly differentially regulated at least 1.5-fold (q < .01). The proximal tendon displayed increased expression of genes resembling a tendon phenotype and increased expression of nerve cell markers. The distal region displayed increases in genes involved in extracellular matrix synthesis and remodeling, immune cell regulation, and a phenotype similar to cartilage and bone. There was a 3.72-fold increase in Young modulus from the proximal to middle region (P < .01) and an additional 1.34-fold increase from the middle to distal region (P = .027). CONCLUSION: Within a single tendon, there are region-specific transcriptomic signatures and mechanical properties, and there is likely a gradient in the biological and functional phenotype from the proximal origin at the MTJ to the distal insertion at the enthesis. CLINICAL RELEVANCE: These findings improve our understanding of the underlying biological heterogeneity of tendon tissue and will help inform the future targeted use of regenerative medicine and tissue engineering strategies for patients with tendon disorders.

DEEPsc: A Deep Learning-Based Map Connecting Single-Cell Transcriptomics and Spatial Imaging Data.

Single-cell RNA sequencing (scRNA-seq) data provides unprecedented information on cell fate decisions; however, the spatial arrangement of cells is often lost. Several recent computational methods have been developed to impute spatial information onto a scRNA-seq dataset through analyzing known spatial expression patterns of a small subset of genes known as a reference atlas. However, there is a lack of comprehensive analysis of the accuracy, precision, and robustness of the mappings, along with the generalizability of these methods, which are often designed for specific systems. We present a system-adaptive deep learning-based method (DEEPsc) to impute spatial information onto a scRNA-seq dataset from a given spatial reference atlas. By introducing a comprehensive set of metrics that evaluate the spatial mapping methods, we compare DEEPsc with four existing methods on four biological systems. We find that while DEEPsc has comparable accuracy to other methods, an improved balance between precision and robustness is achieved. DEEPsc provides a data-adaptive tool to connect scRNA-seq datasets and spatial imaging datasets to analyze cell fate decisions. Our implementation with a uniform API can serve as a portal with access to all the methods investigated in this work for spatial exploration of cell fate decisions in scRNA-seq data. All methods evaluated in this work are implemented as an open-source software with a uniform interface.

Using the Allen gene expression atlas of the adult mouse brain to gain further insight into the physiological significance of TAG-1/Contactin-2.

The anatomic gene expression atlas (AGEA) of the adult mouse brain of the Allen Institute for Brain Science is a transcriptome-based atlas of the adult C57Bl/6 J mouse brain, based on the extensive in situ hybridization dataset of the Institute. This spatial mapping of the gene expression levels of mice under baseline conditions could assist in the formation of new, reasonable transcriptome-derived hypotheses on brain structure and underlying biochemistry, which could also have functional implications. The aim of this work is to use the data of the AGEA (in combination with Tabula Muris, a compendium of single cell transcriptome data collected from mice, enabling direct and controlled comparison of gene expression among cell types) to provide further insights into the physiology of TAG-1/Contactin-2 and its interactions, by presenting the expression of the corresponding gene across the adult mouse brain under baseline conditions and to investigate any spatial genomic correlations between TAG-1/Contactin-2 and its interacting proteins and markers of mature and immature oligodendrocytes, based on the pre-existing experimental or bibliographical evidence. The across-brain correlation analysis on the gene expression intensities showed a positive spatial correlation of TAG-1/Contactin-2 with the gene expression of Plp1, Myrf, Mbp, Mog, Cldn11, Bace1, Kcna1, Kcna2, App and Nfasc and a negative spatial correlation with the gene expression of Cspg4, Pdgfra, L1cam, Ncam1, Ncam2 and Ptprz1. Spatially correlated genes are mainly expressed by mature oligodendrocytes (like Cntn2), while spatially anticorrelated genes are mainly expressed by oligodendrocyte precursor cells. According to the data presented in this work, we propose that even though Contactin-2 expression during development correlates with high plasticity events, such as neuritogenesis, in adulthood it correlates with pathways characterized by low plasticity.