A taxonomy of transcriptomic cell types across the isocortex and hippocampal formation.

The isocortex and hippocampal formation (HPF) in the mammalian brain play critical roles in perception, cognition, emotion, and learning. We profiled ∼1.3 million cells covering the entire adult mouse isocortex and HPF and derived a transcriptomic cell-type taxonomy revealing a comprehensive repertoire of glutamatergic and GABAergic neuron types. Contrary to the traditional view of HPF as having a simpler cellular organization, we discover a complete set of glutamatergic types in HPF homologous to all major subclasses found in the six-layered isocortex, suggesting that HPF and the isocortex share a common circuit organization. We also identify large-scale continuous and graded variations of cell types along isocortical depth, across the isocortical sheet, and in multiple dimensions in hippocampus and subiculum. Overall, our study establishes a molecular architecture of the mammalian isocortex and hippocampal formation and begins to shed light on its underlying relationship with the development, evolution, connectivity, and function of these two brain structures.

Conserved cell types with divergent features in human versus mouse cortex.

Elucidating the cellular architecture of the human cerebral cortex is central to understanding our cognitive abilities and susceptibility to disease. Here we used single-nucleus RNA-sequencing analysis to perform a comprehensive study of cell types in the middle temporal gyrus of human cortex. We identified a highly diverse set of excitatory and inhibitory neuron types that are mostly sparse, with excitatory types being less layer-restricted than expected. Comparison to similar mouse cortex single-cell RNA-sequencing datasets revealed a surprisingly well-conserved cellular architecture that enables matching of homologous types and predictions of properties of human cell types. Despite this general conservation, we also found extensive differences between homologous human and mouse cell types, including marked alterations in proportions, laminar distributions, gene expression and morphology. These species-specific features emphasize the importance of directly studying human brain.

Shared and distinct transcriptomic cell types across neocortical areas.

The neocortex contains a multitude of cell types that are segregated into layers and functionally distinct areas. To investigate the diversity of cell types across the mouse neocortex, here we analysed 23,822 cells from two areas at distant poles of the mouse neocortex: the primary visual cortex and the anterior lateral motor cortex. We define 133 transcriptomic cell types by deep, single-cell RNA sequencing. Nearly all types of GABA (γ-aminobutyric acid)-containing neurons are shared across both areas, whereas most types of glutamatergic neurons were found in one of the two areas. By combining single-cell RNA sequencing and retrograde labelling, we match transcriptomic types of glutamatergic neurons to their long-range projection specificity. Our study establishes a combined transcriptomic and projectional taxonomy of cortical cell types from functionally distinct areas of the adult mouse cortex.

Pharmacogenomic associations of cyclophosphamide pharmacokinetic candidate genes with event-free survival in intermediate-risk rhabdomyosarcoma: A report from the Children's Oncology Group.

BACKGROUND: In vitro data suggest that the growth of rhabdomyosarcoma (RMS) cells is suppressed in a concentration-dependent manner by 4-hydroxycyclophosphamide (4HCY), the principal precursor to the cytotoxic metabolite of cyclophosphamide (CY). Various retrospective studies on the relationship between genes encoding proteins involved in the formation and elimination of 4HCY (i.e., 4HCY pharmacokinetics) and cyclophosphamide (CY) efficacy and toxicity have been conflicting. PROCEDURES: We evaluated germline pharmacogenetics in 262 patients with newly diagnosed intermediate-risk RMS who participated in one prospective Children's Oncology Group clinical trial, ARST0531. Patients were treated with either vincristine/actinomycin/cyclophosphamide (VAC) or VAC alternating with vincristine/irinotecan (VAC/VI). We analyzed the associations between event-free survival and 394 single-nucleotide polymorphisms (SNP) in 14 drug metabolizing enzymes or transporters involved in 4HCY pharmacokinetics. RESULTS: Eight SNPs were associated (p-value < .05 by univariate analysis) with 3-year event-free survival; no SNPs survived a false discovery rate < 0.05. CONCLUSIONS: Our data suggest that a pharmacogenomic approach to therapy personalization of cyclophosphamide in intermediate-risk rhabdomyosarcoma is not viable. Other methods to personalize therapy should be explored.

Pharmacogenomic associations of cyclophosphamide pharmacokinetic candidate genes with 4hydroxycyclophosphamide formation in children with Cancer.

PURPOSE: 4-hydroxycyclophosphamide (4HCY) is the principal precursor to the cytotoxic metabolite of cyclophosphamide (CY), which is often used as first-line treatment of children with cancer. There is conflicting data regarding the relationship between CY efficacy, toxicity, and pharmacokinetics with the genes encoding proteins involved in 4HCY pharmacokinetics, specifically its formation and elimination. METHODS: We evaluated germline pharmacogenetics in children with various malignancies receiving their first CY dose. Using linear regression, we analyzed the associations between two pharmacokinetic outcomes - how fast a child cleared CY (i.e., CY clearance) and the ratio of the 4HCY/CY exposure, specifically area under the plasma concentration-time curve (AUC), and 372 single nucleotide polymorphisms (SNP) in 14 drug-metabolizing transporters or enzymes involved in 4HCY formation or elimination. RESULTS: Age was associated with the ratio of 4HCY/CY AUC (P = 0.004); Chemotherapy regimen was associated with CY clearance (P = 0.003). No SNPs were associated with CY clearance or the ratio of 4HCY/CY AUC after controlling for a false discovery rate. CONCLUSION: Age and chemotherapy regimen, but not germline pharmacogenomics, were associated with CY clearance or the ratio of 4HCY/CY AUC. Other methods, such as metabolomics or lipidomics, should be explored.

COX-1 (PTGS1) and COX-2 (PTGS2) polymorphisms, NSAID interactions, and risk of colon and rectal cancers in two independent populations.

PURPOSE: Nonsteroidal anti-inflammatory drugs (NSAIDs) target the prostaglandin H synthase enzymes, cyclooxygenase (COX)-1 and COX-2, and reduce colorectal cancer risk. Genetic variation in the genes encoding these enzymes may be associated with changes in colon and rectal cancer risk and in NSAID efficacy. METHODS: We genotyped candidate polymorphisms and tag SNPs in PTGS1 (COX-1) and PTGS2 (COX-2) in a population-based case­control study (Diet, Activity and Lifestyle Study, DALS) of colon cancer (n = 1,470 cases/1,837 controls) and rectal cancer (n = 583/775), and independently among cases and controls from the Colon Cancer Family Registry (CCFR; colon n = 959/1,535, rectal n = 505/839). RESULTS: In PTGS2, a functional polymorphism (-765G[C; rs20417) was associated with a twofold increased rectal cancer risk (p = 0.05) in the DALS. This association replicated with a significant nearly fivefold increased risk of rectal cancer in the CCFR study (ORCC vs. GG = 4.88; 95 % CI 1.54­15.45; ORGC vs. GG = 1.36; 95 %CI 0.95­1.94). Genotype­NSAID interactions were observed in the DALS for PTGS1 and rectal cancer risk and for PTGS2 and colon cancer risk, but were no longer significant after correcting for multiple comparisons and did not replicate in the CCFR. No significant associations between PTGS1 polymorphisms and colon or rectal cancer risk were observed.

Comparative transcriptomics reveals human-specific cortical features.

The cognitive abilities of humans are distinctive among primates, but their molecular and cellular substrates are poorly understood. We used comparative single-nucleus transcriptomics to analyze samples of the middle temporal gyrus (MTG) from adult humans, chimpanzees, gorillas, rhesus macaques, and common marmosets to understand human-specific features of the neocortex. Human, chimpanzee, and gorilla MTG showed highly similar cell-type composition and laminar organization as well as a large shift in proportions of deep-layer intratelencephalic-projecting neurons compared with macaque and marmoset MTG. Microglia, astrocytes, and oligodendrocytes had more-divergent expression across species compared with neurons or oligodendrocyte precursor cells, and neuronal expression diverged more rapidly on the human lineage. Only a few hundred genes showed human-specific patterning, suggesting that relatively few cellular and molecular changes distinctively define adult human cortical structure.

Transcriptomic cytoarchitecture reveals principles of human neocortex organization.

Variation in cytoarchitecture is the basis for the histological definition of cortical areas. We used single cell transcriptomics and performed cellular characterization of the human cortex to better understand cortical areal specialization. Single-nucleus RNA-sequencing of 8 areas spanning cortical structural variation showed a highly consistent cellular makeup for 24 cell subclasses. However, proportions of excitatory neuron subclasses varied substantially, likely reflecting differences in connectivity across primary sensorimotor and association cortices. Laminar organization of astrocytes and oligodendrocytes also differed across areas. Primary visual cortex showed characteristic organization with major changes in the excitatory to inhibitory neuron ratio, expansion of layer 4 excitatory neurons, and specialized inhibitory neurons. These results lay the groundwork for a refined cellular and molecular characterization of human cortical cytoarchitecture and areal specialization.

A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain.

The mammalian brain is composed of millions to billions of cells that are organized into numerous cell types with specific spatial distribution patterns and structural and functional properties. An essential step towards understanding brain function is to obtain a parts list, i.e., a catalog of cell types, of the brain. Here, we report a comprehensive and high-resolution transcriptomic and spatial cell type atlas for the whole adult mouse brain. The cell type atlas was created based on the combination of two single-cell-level, whole-brain-scale datasets: a single-cell RNA-sequencing (scRNA-seq) dataset of ~7 million cells profiled, and a spatially resolved transcriptomic dataset of ~4.3 million cells using MERFISH. The atlas is hierarchically organized into five nested levels of classification: 7 divisions, 32 classes, 306 subclasses, 1,045 supertypes and 5,200 clusters. We systematically analyzed the neuronal, non-neuronal, and immature neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell type organization in different brain regions, in particular, a dichotomy between the dorsal and ventral parts of the brain: the dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. We also systematically characterized cell-type specific expression of neurotransmitters, neuropeptides, and transcription factors. The study uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types across the brain, suggesting they mediate a myriad of modes of intercellular communications. Finally, we found that transcription factors are major determinants of cell type classification in the adult mouse brain and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole-mouse-brain transcriptomic and spatial cell type atlas establishes a benchmark reference atlas and a foundational resource for deep and integrative investigations of cell type and circuit function, development, and evolution of the mammalian brain.

Enhancer-AAVs allow genetic access to oligodendrocytes and diverse populations of astrocytes across species.

Proper brain function requires the assembly and function of diverse populations of neurons and glia. Single cell gene expression studies have mostly focused on characterization of neuronal cell diversity; however, recent studies have revealed substantial diversity of glial cells, particularly astrocytes. To better understand glial cell types and their roles in neurobiology, we built a new suite of adeno-associated viral (AAV)-based genetic tools to enable genetic access to astrocytes and oligodendrocytes. These oligodendrocyte and astrocyte enhancer-AAVs are highly specific (usually > 95% cell type specificity) with variable expression levels, and our astrocyte enhancer-AAVs show multiple distinct expression patterns reflecting the spatial distribution of astrocyte cell types. To provide the best glial-specific functional tools, several enhancer-AAVs were: optimized for higher expression levels, shown to be functional and specific in rat and macaque, shown to maintain specific activity in epilepsy where traditional promoters changed activity, and used to drive functional transgenes in astrocytes including Cre recombinase and acetylcholine-responsive sensor iAChSnFR. The astrocyte-specific iAChSnFR revealed a clear reward-dependent acetylcholine response in astrocytes of the nucleus accumbens during reinforcement learning. Together, this collection of glial enhancer-AAVs will enable characterization of astrocyte and oligodendrocyte populations and their roles across species, disease states, and behavioral epochs.

Single-nucleus and single-cell transcriptomes compared in matched cortical cell types.

Transcriptomic profiling of complex tissues by single-nucleus RNA-sequencing (snRNA-seq) affords some advantages over single-cell RNA-sequencing (scRNA-seq). snRNA-seq provides less biased cellular coverage, does not appear to suffer cell isolation-based transcriptional artifacts, and can be applied to archived frozen specimens. We used well-matched snRNA-seq and scRNA-seq datasets from mouse visual cortex to compare cell type detection. Although more transcripts are detected in individual whole cells (~11,000 genes) than nuclei (~7,000 genes), we demonstrate that closely related neuronal cell types can be similarly discriminated with both methods if intronic sequences are included in snRNA-seq analysis. We estimate that the nuclear proportion of total cellular mRNA varies from 20% to over 50% for large and small pyramidal neurons, respectively. Together, these results illustrate the high information content of nuclear RNA for characterization of cellular diversity in brain tissues.

UGT1A6 polymorphism and salicylic acid glucuronidation following aspirin.

OBJECTIVES: In vivo, aspirin (acetylsalicylic acid) is rapidly deacetylated to form salicylic acid, which then undergoes primary or secondary glucuronidation catalyzed by UDP-glucuronosyltransferases (UGTs). The variant UGT1A6*2 (T181A, R184S) is associated with altered enzyme function. Our objective was to compare salicylic acid glucuronidation in individuals with different UGT1A6 genotypes. METHODS: Following orally dosing with 650 mg aspirin, saliva and urine samples were collected over a period of 24 h from healthy individuals with homozygous wild-type UGT1A6 *1/*1 (n=19) and homozygous variant UGT1A6 *2/*2 (T181A, R184S) (n=9) genotypes. RESULTS: No statistically significant differences were observed in salivary pharmacokinetic parameters. Urinary excretion of the sum of aspirin and its metabolites (salicyluric acid, salicyluric acid phenolic glucuronide, salicyl phenolic glucuronide, salicyl acyl glucuronide, salicylic acid) during the early period of 2-4 h of collection was significantly lower in UGT1A6 *1/*1 than in UGT1A6 *2/*2 individuals. Further, UGT1A6 *1/*1 individuals excreted a lower percentage of aspirin and its metabolites in the first 12 h and a greater percentage after 12 h than UGT1A6 *2/*2 individuals. CONCLUSIONS: The variant UGT1A6*2 or polymorphisms in other UGTs that are in linkage disequilibrium with UGT1A6*2 may confer more rapid glucuronidation of salicylic acid than the wild-type UGT1A6 *1/*1.