The organization of spinal neurons: Insights from single cell sequencing

To understand how the spinal cord enacts complex sensorimotor functions, researchers have studied, classified, and functionally probed it's many neuronal populations for over a century. Recent developments in single-cell RNA-sequencing can characterize the gene expression signatures of the entire set of spinal neuron types and can simultaneously provide an unbiased view of their relationships to each other. This approach has revealed that the location of neurons predicts transcriptomic variability, as dorsal spinal neurons become highly distinct over development as ventral spinal neurons become less so. Temporal specification is also a major source of gene expression variation, subdividing many of the canonical embryonic lineage domains. Together, birthdate and cell body location are fundamental organizing features of spinal neuron diversity.

Universal DNA methylation age across mammalian tissues

Aging, often considered a result of random cellular damage, can be accurately estimated using DNA methylation profiles, the foundation of pan-tissue epigenetic clocks. Here, we demonstrate the development of universal pan-mammalian clocks, using 11,754 methylation arrays from our Mammalian Methylation Consortium, which encompass 59 tissue types across 185 mammalian species. These predictive models estimate mammalian tissue age with high accuracy (r > 0.96). Age deviations correlate with human mortality risk, mouse somatotropic axis mutations and caloric restriction. We identified specific cytosines with methylation levels that change with age across numerous species. These sites, highly enriched in polycomb repressive complex 2-binding locations, are near genes implicated in mammalian development, cancer, obesity and longevity. Our findings offer new evidence suggesting that aging is evolutionarily conserved and intertwined with developmental processes across all mammals.

Molecular and spatial profiling of the paraventricular nucleus of the thalamus

The paraventricular nucleus of the thalamus (PVT) is known to regulate various cognitive and behavioral processes. However, while functional diversity among PVT circuits has often been linked to cellular differences, the molecular identity and spatial distribution of PVT cell types remain unclear. To address this gap, here we used single nucleus RNA sequencing (snRNA-seq) and identified five molecularly distinct PVT neuronal subtypes in the mouse brain. Additionally, multiplex fluorescent in situ hybridization of top marker genes revealed that PVT subtypes are organized by a combination of previously unidentified molecular gradients. Lastly, comparing our dataset with a recently published single-cell sequencing atlas of the thalamus yielded novel insight into the PVT's connectivity with the cortex, including unexpected innervation of auditory and visual areas. This comparison also revealed that our data contains a largely non-overlapping transcriptomic map of multiple midline thalamic nuclei. Collectively, our findings uncover previously unknown features of the molecular diversity and anatomical organization of the PVT and provide a valuable resource for future investigations.

A cellular taxonomy of the adult human spinal cord

The mammalian spinal cord functions as a community of cell types for sensory processing, autonomic control, and movement. While animal models have advanced our understanding of spinal cellular diversity, characterizing human biology directly is important to uncover specialized features of basic function and human pathology. Here, we present a cellular taxonomy of the adult human spinal cord using single-nucleus RNA sequencing with spatial transcriptomics and antibody validation. We identified 29 glial clusters and 35 neuronal clusters, organized principally by anatomical location. To demonstrate the relevance of this resource to human disease, we analyzed spinal motoneurons, which degenerate in amyotrophic lateral sclerosis (ALS) and other diseases. We found that compared with other spinal neurons, human motoneurons are defined by genes related to cell size, cytoskeletal structure, and ALS, suggesting a specialized molecular repertoire underlying their selective vulnerability. We include a web resource to facilitate further investigations into human spinal cord biology.

The neurons that restore walking after paralysis

A spinal cord injury interrupts pathways from the brain and brainstem that project to the lumbar spinal cord, leading to paralysis. Here we show that spatiotemporal epidural electrical stimulation (EES) of the lumbar spinal cord1–3 applied during neurorehabilitation4,5 (EESREHAB) restored walking in nine individuals with chronic spinal cord injury. This recovery involved a reduction in neuronal activity in the lumbar spinal cord of humans during walking. We hypothesized that this unexpected reduction reflects activity-dependent selection of specific neuronal subpopulations that become essential for a patient to walk after spinal cord injury. To identify these putative neurons, we modelled the technological and therapeutic features underlying EESREHAB in mice. We applied single-nucleus RNA sequencing6–9 and spatial transcriptomics10,11 to the spinal cords of these mice to chart a spatially resolved molecular atlas of recovery from paralysis. We then employed cell type12,13 and spatial prioritization to identify the neurons involved in the recovery of walking. A single population of excitatory interneurons nested within intermediate laminae emerged. Although these neurons are not required for walking before spinal cord injury, we demonstrate that they are essential for the recovery of walking with EES following spinal cord injury. Augmenting the activity of these neurons phenocopied the recovery of walking enabled by EESREHAB, whereas ablating them prevented the recovery of walking that occurs spontaneously after moderate spinal cord injury. We thus identified a recovery-organizing neuronal subpopulation that is necessary and sufficient to regain walking after paralysis. Moreover, our methodology establishes a framework for using molecular cartography to identify the neurons that produce complex behaviours.

Single cell atlas of spinal cord injury in mice reveals a pro-regenerative signature in spinocerebellar neurons

After spinal cord injury, tissue distal to the lesion contains undamaged cells that could support or augment recovery. Targeting these cells requires a clearer understanding of their injury responses and capacity for repair. Here, we use single nucleus RNA sequencing to profile how each cell type in the lumbar spinal cord changes after a thoracic injury in mice. We present an atlas of these dynamic responses across dozens of cell types in the acute, subacute, and chronically injured spinal cord. Using this resource, we find rare spinal neurons that express a signature of regeneration in response to injury, including a major population that represent spinocerebellar projection neurons. We characterize these cells anatomically and observed axonal sparing, outgrowth, and remodeling in the spinal cord and cerebellum. Together, this work provides a key resource for studying cellular responses to injury and uncovers the spontaneous plasticity of spinocerebellar neurons, uncovering a potential candidate for targeted therapy.

Matson et al. performed single nucleus sequencing of the “spared” spinal cord tissue distal to an injury in mice. They found that spinocerebellar neurons expressed a pro-regenerative gene signature and showed axon outgrowth after injury.

Intersectional genetic tools to study skilled reaching in mice

Reaching to grasp is an evolutionarily conserved behavior and a crucial part of the motor repertoire in mammals. As it is studied in the laboratory, reaching has become the prototypical example of dexterous forelimb movements, illuminating key principles of motor control throughout the spinal cord, brain, and peripheral nervous system. Here, we (1) review the motor elements or phases that comprise the reach, grasp, and retract movements of reaching behavior, (2) highlight the role of intersectional genetic tools in linking these movements to their neuronal substrates, (3) describe spinal cord cell types and their roles in skilled reaching, and (4) how descending pathways from the brain and the sensory systems contribute to skilled reaching. We emphasize that genetic perturbation experiments can pin-point the neuronal substrates of specific phases of reaching behavior.

A harmonized atlas of mouse spinal cord cell types and their spatial organization

Single-cell RNA sequencing data can unveil the molecular diversity of cell types. Cell type atlases of the mouse spinal cord have been published in recent years but have not been integrated together. Here, we generate an atlas of spinal cell types based on single-cell transcriptomic data, unifying the available datasets into a common reference framework. We report a hierarchical structure of postnatal cell type relationships, with location providing the highest level of organization, then neurotransmitter status, family, and finally, dozens of refined populations. We validate a combinatorial marker code for each neuronal cell type and map their spatial distributions in the adult spinal cord. We also show complex lineage relationships among postnatal cell types. Additionally, we develop an open-source cell type classifier, SeqSeek, to facilitate the standardization of cell type identification. This work provides an integrated view of spinal cell types, their gene expression signatures, and their molecular organization.

Confronting false discoveries in single-cell differential expression

Differential expression analysis in single-cell transcriptomics enables the dissection of cell-type-specific responses to perturbations such as disease, trauma, or experimental manipulations. While many statistical methods are available to identify differentially expressed genes, the principles that distinguish these methods and their performance remain unclear. Here, we show that the relative performance of these methods is contingent on their ability to account for variation between biological replicates. Methods that ignore this inevitable variation are biased and prone to false discoveries. Indeed, the most widely used methods can discover hundreds of differentially expressed genes in the absence of biological differences. To exemplify these principles, we exposed true and false discoveries of differentially expressed genes in the injured mouse spinal cord.

Cerebellospinal Neurons Regulate Motor Performance and Motor Learning

To understand the neural basis of behavior, it is important to reveal how movements are planned, executed, and refined by networks of neurons distributed throughout the nervous system. Here, we report the neuroanatomical organization and behavioral roles of cerebellospinal (CeS) neurons. Using intersectional genetic techniques, we find that CeS neurons constitute a small minority of excitatory neurons in the fastigial and interpositus deep cerebellar nuclei, target pre-motor circuits in the ventral spinal cord and the brain, and control distinct aspects of movement. CeS neurons that project to the ipsilateral cervical cord are required for skilled forelimb performance, while CeS neurons that project to the contralateral cervical cord are involved in skilled locomotor learning. Together, this work establishes CeS neurons as a critical component of the neural circuitry for skilled movements and provides insights into the organizational logic of motor networks.

Sathyamurthy et al. define the organization, function, and targets of cerebellospinal neurons, revealing a direct link between the deep cerebellar nuclei and motor execution circuits in the spinal cord and demonstrating a role for these neurons in motor control.

Selecting single cell clustering parameter values using subsampling-based robustness metrics

Background

Generating and analysing single-cell data has become a widespread approach to examine tissue heterogeneity, and numerous algorithms exist for clustering these datasets to identify putative cell types with shared transcriptomic signatures. However, many of these clustering workflows rely on user-tuned parameter values, tailored to each dataset, to identify a set of biologically relevant clusters. Whereas users often develop their own intuition as to the optimal range of parameters for clustering on each data set, the lack of systematic approaches to identify this range can be daunting to new users of any given workflow. In addition, an optimal parameter set does not guarantee that all clusters are equally well-resolved, given the heterogeneity in transcriptomic signatures in most biological systems.

Results

Here, we illustrate a subsampling-based approach (chooseR) that simultaneously guides parameter selection and characterizes cluster robustness. Through bootstrapped iterative clustering across a range of parameters, chooseR was used to select parameter values for two distinct clustering workflows (Seurat and scVI). In each case, chooseR identified parameters that produced biologically relevant clusters from both well-characterized (human PBMC) and complex (mouse spinal cord) datasets. Moreover, it provided a simple “robustness score” for each of these clusters, facilitating the assessment of cluster quality.

Conclusion

chooseR is a simple, conceptually understandable tool that can be used flexibly across clustering algorithms, workflows, and datasets to guide clustering parameter selection and characterize cluster robustness.

Massively Parallel Single Nucleus Transcriptional Profiling Defines Spinal Cord Neurons and Their Activity during Behavior

To understand the cellular basis of behavior, it is necessary to know the cell types that exist in the nervous system and their contributions to function. Spinal networks are essential for sensory processing and motor behavior and provide a powerful system for identifying the cellular correlates of behavior. Here, we used massively parallel single nucleus RNA sequencing (snRNA-seq) to create an atlas of the adult mouse lumbar spinal cord. We identified and molecularly characterized 43 neuronal populations. Next, we leveraged the snRNA-seq approach to provide unbiased identification of neuronal populations that were active following a sensory and a motor behavior, using a transcriptional signature of neuronal activity. This approach can be used in the future to link single nucleus gene expression data with dynamic biological responses to behavior, injury, and disease.

Decoding Cell Type Diversity Within the Spinal Cord

To understand fundamental mechanisms of mammalian spinal cord function, it is necessary to reveal the diverse array of constituent spinal "cell types" - populations that can be consistently identified because they share a unique and cohesive set of characteristics. Many parameters can contribute to the definition of a spinal cord cell type, including location, morphology, lineage, electrophysiological properties, circuit features, gene expression patterns, and behavioral contribution. While it is not necessary for all of these features to align completely at all times to identify an individual cell type, a correlation of these characteristics paints a rich portrait of cell identity. This review will summarize recent advances in the identification of mammalian spinal cord neuronal cell types and will highlight the power of transcriptional profiling to identify and characterize the cell types of the spinal cord.

Isolation of Adult Spinal Cord Nuclei for Massively Parallel Single-nucleus RNA Sequencing

Probing an individual cell's gene expression enables the identification of cell type and cell state. Single-cell RNA sequencing has emerged as a powerful tool for studying transcriptional profiles of cells, particularly in heterogeneous tissues such as the central nervous system. However, dissociation methods required for single cell sequencing can lead to experimental changes in the gene expression and cell death. Furthermore, these methods are generally restricted to fresh tissue, thus limiting studies on archival and bio-bank material. Single nucleus RNA sequencing (snRNA-Seq) is an appealing alternative for transcriptional studies, given that it accurately identifies cell types, permits the study of tissue that is frozen or difficult to dissociate, and reduces dissociation-induced transcription. Here, we present a high-throughput protocol for rapid isolation of nuclei for downstream snRNA-Seq. This method enables isolation of nuclei from fresh or frozen spinal cord samples and can be combined with two massively parallel droplet encapsulation platforms.

Graded Arrays of Spinal and Supraspinal V2a Interneuron Subtypes Underlie Forelimb and Hindlimb Motor Control

The spinal cord contains neural networks that enable regionally distinct motor outputs along the body axis. Nevertheless, it remains unclear how segment-specific motor computations are processed because the cardinal interneuron classes that control motor neurons appear uniform at each level of the spinal cord. V2a interneurons are essential to both forelimb and hindlimb movements, and here we identify two major types that emerge during development: type I neurons marked by high Chx10 form recurrent networks with neighboring spinal neurons and type II neurons that downregulate Chx10 and project to supraspinal structures. Types I and II V2a interneurons are arrayed in counter-gradients, and this network activates different patterns of motor output at cervical and lumbar levels. Single-cell RNA sequencing (RNA-seq) revealed type I and II V2a neurons are each comprised of multiple subtypes. Our findings uncover a molecular and anatomical organization of V2a interneurons reminiscent of the orderly way motor neurons are divided into columns and pools.

The association between physical activity and cognition in men with and without HIV infection

OBJECTIVES

HIV-associated neurocognitive disorders are highly prevalent, and physical activity (PA) is a modifiable behaviour that may affect neurocognitive function. Our objective was to determine the association between PA and neurocognitive function and the effect of HIV on this association.

METHODS

PA was assessed in the Multicenter AIDS Cohort Study with the International Physical Activity Questionnaire. A neuropsychological test battery assessed global impairment and domain-specific impairment (executive function, speed of processing, working memory, learning, memory, and motor function) every 2 years. Semiannually, the Symbol Digit Modalities Test and Trail Making Test Parts A and B were performed. Adjusted logistic regression models were used to assess the PA-neurocognitive function association. Using longitudinal data, we also assessed the PA category-decline of neurocognitive function association with multivariate simple regression.

RESULTS

Of 601 men, 44% were HIV-infected. Low, moderate, and high PA was reported in 27%, 25%, and 48% of the HIV-infected men vs. 19%, 32% and 49% of the HIV-uninfected men, respectively. High PA was associated with lower odds of impairment of learning, memory, and motor function [odds ratio (OR) ranging from 0.52 to 0.57; P < 0.05 for all]. The high PA-global impairment association OR was 0.63 [95% confidence interval (CI) 0.39, 1.02]. Among HIV-infected men only, across multiple domains, the high PA-impairment association was even more pronounced (OR from 0.27 to 0.49). Baseline high/moderate PA was not associated with decline of any domain score over time. HIV infection was marginally associated with a higher speed of decline in motor function.

CONCLUSIONS

A protective effect of high PA on impairment in neurocognitive domains was observed cross-sectionally. Longitudinal PA measurements are needed to elucidate the PA-neurocognitive function relationship over time.

Satb2 Is Required for the Development of a Spinal Exteroceptive Microcircuit that Modulates Limb Position

Motor behaviors such as walking or withdrawing the limb from a painful stimulus rely upon integrative multimodal sensory circuitry to generate appropriate muscle activation patterns. Both the cellular components and the molecular mechanisms that instruct the assembly of the spinal sensorimotor system are poorly understood. Here we characterize the connectivity pattern of a sub-population of lamina V inhibitory sensory relay neurons marked during development by the nuclear matrix and DNA binding factor Satb2 (ISR(Satb2)). ISR(Satb2) neurons receive inputs from multiple streams of sensory information and relay their outputs to motor command layers of the spinal cord. Deletion of the Satb2 transcription factor from ISR(Satb2) neurons perturbs their cellular position, molecular profile, and pre- and post-synaptic connectivity. These alterations are accompanied by abnormal limb hyperflexion responses to mechanical and thermal stimuli and during walking. Thus, Satb2 is a genetic determinant that mediates proper circuit development in a core sensory-to-motor spinal network.

Involvement of tumor suppressors PTEN and p53 in the formation of multiple subtypes of liposarcoma

Liposarcoma (LPS) is a type of soft tissue sarcoma that mostly occurs in adults, and in humans is characterized by amplifications of MDM2 and CDK4. The molecular pathogenesis of this malignancy is still poorly understood and, therefore, we developed a mouse model with conditional inactivation of PTEN and p53 to investigate these pathways in the progression of the disease. We show that deletion of these two tumor suppressors cooperate in the formation of multiple subtypes of LPS (from well-differentiated LPS to pleomorphic LPS). In addition, progression of the tumors is further characterized by the expression of D cyclins and CDK4/6, which allow for continued cell division. Microarray analysis also revealed novel genes that are differentially expressed between different subtypes of LPS, which could aid in understanding the disease and to unravel potential new therapeutic targets.

Identification of a cellular node for motor control pathways

The rich behavioral repertoire of animals is encoded in the CNS as a set of motorneuron activation patterns, also called 'motor synergies'. However, the neurons that orchestrate these motor programs as well as their cellular properties and connectivity are poorly understood. Here we identify a population of molecularly defined motor synergy encoder (MSE) neurons in the mouse spinal cord that may represent a central node in neural pathways for voluntary and reflexive movement. This population receives direct inputs from the motor cortex and sensory pathways and, in turn, has monosynaptic outputs to spinal motorneurons. Optical stimulation of MSE neurons drove reliable patterns of activity in multiple motor groups, and we found that the evoked motor patterns varied on the basis of the rostrocaudal location of the stimulated MSE. We speculate that these neurons comprise a cellular network for encoding coordinated motor output programs.

MicroRNAs, miR-154, miR-299-5p, miR-376a, miR-376c, miR-377, miR-381, miR-487b, miR-485-3p, miR-495 and miR-654-3p, mapped to the 14q32.31 locus, regulate proliferation, apoptosis, migration and invasion in metastatic prostate cancer cells

miRNAs act as oncogenes or tumor suppressors in a wide variety of human cancers, including prostate cancer (PCa). We found a severe and consistent downregulation of miRNAs, miR-154, miR-299-5p, miR-376a, miR-376c, miR-377, miR-381, miR-487b, miR-485-3p, miR-495 and miR-654-3p, mapped to the 14q32.31 region in metastatic cell lines as compared with normal prostatic epithelial cells (PrEC). In specimens of human prostate (28 normals, 99 primary tumors and 13 metastases), lower miRNA levels correlated significantly with a higher incidence of metastatic events and higher prostate specific antigen (PSA) levels, with similar trends observed for lymph node invasion and the Gleason score. We transiently transfected 10 members of the 14q32.31 cluster in normal prostatic epithelial cell lines and characterized their affect on malignant cell behaviors, including proliferation, apoptosis, migration and invasion. Finally, we identified FZD4, a gene important for epithelial-to-mesenchymal transition in (PCa), as a target of miR-377.

DNA methylation silences miR-132 in prostate cancer

Silencing of microRNAs (miRNAs) by promoter CpG island methylation may be an important mechanism in prostate carcinogenesis. To screen for epigenetically silenced miRNAs in prostate cancer (PCa), we treated prostate normal epithelial and carcinoma cells with 5-aza-2'-deoxycytidine (AZA) and subsequently examined expression changes of 650 miRNAs by megaplex stemloop reverse transcription-quantitative PCR. After applying a selection strategy, we analyzed the methylation status of CpG islands upstream to a subset of miRNAs by methylation-specific PCR. The CpG islands of miR-18b, miR-132, miR-34b/c, miR-148a, miR-450a and miR-542-3p showed methylation patterns congruent with their expression modulations in response to AZA. Methylation analysis of these CpG islands in a panel of 50 human prostate carcinoma specimens and 24 normal controls revealed miR-132 to be methylated in 42% of human cancer cases in a manner positively correlated to total Gleason score and tumor stage. Expression analysis of miR-132 in our tissue panel confirmed its downregulation in methylated tumors. Re-expression of miR-132 in PC3 cells induced cell detachment followed by cell death (anoikis). Two pro-survival proteins-heparin-binding epidermal growth factor and TALIN2-were confirmed as direct targets of miR-132. The results of this study point to miR-132 as a methylation-silenced miRNA with an antimetastatic role in PCa controlling cellular adhesion.

Proliferation and tissue remodeling in cancer: the hallmarks revisited

Although cancers are highly heterogeneous at the genomic level, they can manifest common patterns of gene expression. Here, we use gene expression signatures to interrogate two major processes in cancer, proliferation and tissue remodeling. We demonstrate that proliferation and remodeling signatures are partially independent and result in four distinctive cancer subtypes. Cancers with the proliferation signature are characterized by signatures of p53 and PTEN inactivation and concomitant Myc activation. In contrast, remodeling correlates with RAS, HIF-1α and NFκB activation. From the metabolic point of view, proliferation is associated with upregulation of glycolysis and serine/glycine metabolism, whereas remodeling is characterized by a downregulation of oxidative phosphorylation. Notably, the proliferation signature correlates with poor outcome in lung, prostate, breast and brain cancer, whereas remodeling increases mortality rates in colorectal and ovarian cancer.

PCA and clustering reveal alternate mtDNA phylogeny of N and M clades

Phylogenetic trees based on mtDNA polymorphisms are often used to infer the history of recent human migrations. However, there is no consensus on which method to use. Most methods make strong assumptions which may bias the choice of polymorphisms and result in computational complexity which limits the analysis to a few samples/polymorphisms. For example, parsimony minimizes the number of mutations, which biases the results to minimizing homoplasy events. Such biases may miss the global structure of the polymorphisms altogether, with the risk of identifying a "common" polymorphism as ancient without an internal check on whether it either is homoplasic or is identified as ancient because of sampling bias (from oversampling the population with the polymorphism). A signature of this problem is that different methods applied to the same data or the same method applied to different datasets results in different tree topologies. When the results of such analyses are combined, the consensus trees have a low internal branch consensus. We determine human mtDNA phylogeny from 1737 complete sequences using a new, direct method based on principal component analysis (PCA) and unsupervised consensus ensemble clustering. PCA identifies polymorphisms representing robust variations in the data and consensus ensemble clustering creates stable haplogroup clusters. The tree is obtained from the bifurcating network obtained when the data are split into k = 2,3,4,...,kmax clusters, with equal sampling from each haplogroup. Our method assumes only that the data can be clustered into groups based on mutations, is fast, is stable to sample perturbation, uses all significant polymorphisms in the data, works for arbitrary sample sizes, and avoids sample choice and haplogroup size bias. The internal branches of our tree have a 90% consensus accuracy. In conclusion, our tree recreates the standard phylogeny of the N, M, L0/L1, L2, and L3 clades, confirming the African origin of modern humans and showing that the M and N clades arose in almost coincident migrations. However, the N clade haplogroups split along an East-West geographic divide, with a "European R clade" containing the haplogroups H, V, H/V, J, T, and U and a "Eurasian N subclade" including haplogroups B, R5, F, A, N9, I, W, and X. The haplogroup pairs (N9a, N9b) and (M7a, M7b) within N and M are placed in nonnearest locations in agreement with their expected large TMRCA from studies of their migrations into Japan. For comparison, we also construct consensus maximum likelihood, parsimony, neighbor joining, and UPGMA-based trees using the same polymorphisms and show that these methods give consistent results only for the clade tree. For recent branches, the consensus accuracy for these methods is in the range of 1-20%. From a comparison of our haplogroups to two chimp and one bonobo sequences, and assuming a chimp-human coalescent time of 5 million years before present, we find a human mtDNA TMRCA of 206,000 +/- 14,000 years before present.

Nonequilibrium mechanics and dynamics of motor-activated gels

The mechanics of cells is strongly affected by molecular motors that generate forces in the cellular cytoskeleton. We develop a model for cytoskeletal networks driven out of equilibrium by molecular motors exerting transient contractile stresses. Using this model we show how motor activity can dramatically increase the network's bulk elastic moduli. We also show how motor binding kinetics naturally leads to enhanced low-frequency stress fluctuations that result in nonequilibrium diffusive motion within an elastic network, as seen in recent in vitro and in vivo experiments.

Discrepancies between self-report and objective measures for stimulant drug use in HIV: cognitive, medication adherence and psychological correlates

While it has long been recognized that self-reported drug use may be at variance with objectively obtained evidence such as urine toxicology assays, few studies have explored the behavioral correlates of such discrepancies. Here we compared self-reported and objective measures of stimulant drug use for 162 HIV infected individuals and identified a sub-group with discrepancies between data obtained via the two methods. Results showed poorer neurocognitive performance (attention, learning/memory) and lower medication adherence rates for the discrepant group as compared to those who either acknowledged their drug use or accurately denied recent stimulant use. Using the Millon Clinical Multiaxial Inventory-III, it was also found that those in the discrepant group were more hesitant to reveal psychopathology. Comparisons of self-reported and objectively measured medication adherence data are also discussed.

The P53 pathway: what questions remain to be explored?

The p53 pathway is composed of hundreds of genes and their products that respond to a wide variety of stress signals. These responses to stress include apoptosis, cellular senescence or cell cycle arrest. In addition the p53-regulated genes produce proteins that communicate these stress signals to adjacent cells, prevent and repair damaged DNA and create feedback loops that enhance or attenuate p53 activity and communicate with other signal transduction pathways. Many questions remain to be explored in our understanding of how this network of genes plays a role in protection from cancers, therapy and integrating the homeostatic mechanisms of stress management and fidelity in a cell and organism. The goal of this chapter is to elucidate some of those questions and suggest new directions for this area of research.

Proposal of a model of mammalian neural induction

How does the vertebrate embryo make a nervous system? This complex question has been at the center of developmental biology for many years. The earliest step in this process - the induction of neural tissue - is intimately linked to patterning of the entire early embryo, and the molecular and embryological of basis these processes are beginning to emerge. Here, we analyze classic and cutting-edge findings on neural induction in the mouse. We find that data from genetics, tissue explants, tissue grafting, and molecular marker expression support a coherent framework for mammalian neural induction. In this model, the gastrula organizer of the mouse embryo inhibits BMP signaling to allow neural tissue to form as a default fate-in the absence of instructive signals. The first neural tissue induced is anterior and subsequent neural tissue is posteriorized to form the midbrain, hindbrain, and spinal cord. The anterior visceral endoderm protects the pre-specified anterior neural fate from similar posteriorization, allowing formation of forebrain. This model is very similar to the default model of neural induction in the frog, thus bridging the evolutionary gap between amphibians and mammals.

A single nucleotide polymorphism in the p53 pathway interacts with gender, environmental stresses and tumor genetics to influence cancer in humans

Cancer biology finds itself in a post-genomic era and the hopes of using inherited genetic variants to improve prevention and treatment strategies are widespread. One of the largest types of inherited genetic variation is the single nucleotide polymorphism (SNP), of which there are at least 4.5 million. The challenge now becomes how to discover which polymorphisms alter cancer in humans and how to begin to understand their mechanism of action. In this report, a series of recent publications will be reviewed that have studied a polymorphism in the p53 tumor suppressor pathway, MDM2 SNP309. These reports have lent insights into how germline genetic variants of the p53 pathway could interact with gender, environmental stresses and tumor genetics to affect cancer in humans. Importantly, these observations have also exposed potential nodes of intervention, which could prove valuable in both the prevention and treatment of this disease in humans.

Imaging the Sublimation Dynamics of Colloidal Crystallites

We studied the kinetics of sublimating crystals with single-particle resolution by experiments with colloidal spheres and by computer simulations. A short-range attraction between spheres led to crystallites one to three layers thick. The spheres were tracked with optical microscopy while the attraction was reduced and the crystals sublimated. Large crystallites sublimated by escape of particles from the perimeter. The rate of shrinkage was greatly enhanced, however, when the size decreased to less than 20 to 50 particles, depending on the location in the phase diagram. At this size, the crystallites transformed into a dense amorphous structure, which rapidly vaporized. The enhancement of kinetics by metastable or unstable phases may play a major role in the melting, freezing, and annealing of crystals.

p53-Mdm2 loop controlled by a balance of its feedback strength and effective dampening using ATM and delayed feedback

When the genomic integrity of a cell is challenged, its fate is determined in part by signals conveyed by the p53 tumour suppressor protein. It was observed recently that such signals are not simple gradations of p53 concentration, but rather a counter-intuitive limit-cycle behaviour. Based on a careful mathematical interpretation of the experimental body of knowledge, we propose a model for the p53 signalling network and characterise the p53 stability and oscillatory dynamics. In our model, ATM, a protein that senses DNA damage, activates p53 by phosphorylation. In its active state, p53 has a decreased degradation rate and an enhanced transactivation of Mdm2, a gene whose protein product Mdm2 tags p53 for degradation. Thus the p53-Mdm2 system forms a negative feedback loop. However, the feedback in this loop is delayed, as the pool of Mdm2 molecules being induced by p53 at a given time will mark for degradation the pool of p53 molecules at some later time, after the Mdm2 molecules have been transcribed, exported out of the nucleus, translated and transported back into the nucleus. The analysis of our model demonstrates how this time lag combines with the ATM-controlled feedback strength and effective dampening of the negative feedback loop to produce limit-cycle oscillations. The picture that emerges is that ATM, once activated by DNA damage, makes the p53-Mdm2 oscillator undergo a supercritical Hopf bifurcation. This approach yields an improved understanding of the global dynamics and bifurcation structure of our time-delayed, negative feedback model and allows for predictions of the behaviour of the p53 system under different perturbations.

Correlated fluctuations of microparticles in viscoelastic solutions: quantitative measurement of material properties by microrheology in the presence of optical traps

The Brownian motions of microscopic particles in viscous or viscoelastic fluids can be used to measure rheological properties. This is the basis of recently developed one- and two-particle microrheology techniques. For increased temporal and spatial resolution, some microrheology techniques employ optical traps, which introduce additional forces on the particles. We have systematically studied the effect that confinement of particles by optical traps has on their auto- and cross-correlated fluctuations. We show that trapping causes anticorrelations in the motion of two particles at low frequencies. We demonstrate how these anticorrelations depend on trap strength and the shear modulus of viscoelastic media. We present a method to account for the effects of optical traps, which permits the quantitative measurement of viscoelastic properties in one- and two-particle microrheology over an extended frequency range in a variety of viscous and viscoelastic media.

GDF3 at the crossroads of TGF-beta signaling

Embryonic stem (ES) cells present an excellent system for addressing the relevance of our current knowledge about how cell fate is determined and how cells integrate multiple signals into a single outcome as a function of time. Many of the factors that mediate these phenomena have been discovered through classical embryological experiments and are organized into several major signal transduction pathways including TGF-beta/BMP, Jak-STAT, Hedgehog, Wnt, Notch and FGF/MAPK. This review will summarize the current understanding of TGF-beta signaling in ES and focus on early embryological roles of the TGF-beta member, GDF-3. GDF-3 is associated with the undifferentiated state of ES cells and two recent and contradictory reports examined the function and mechanism of GDF-3 in the context of both stem cells and early embryonic differentiations. While Levine and Brivanlou found that GDF-3 inhibits its own subfamily members (the BMPs), Chen and colleagues found that GDF-3 acts as a nodal-like TGF-beta ligand. These combined findings raise the intriguing possibility that GDF-3 acts as a bi-functional protein, to regulate the balance between the two modes of TGF-beta signaling.

GDF3, a BMP inhibitor, regulates cell fate in stem cells and early embryos

The TGFbeta superfamily of ligands plays key functions in development and disease. In both human and mouse embryonic stem cells, a member of this family, GDF3, is specifically expressed in the pluripotent state. We show that GDF3 is an inhibitor of its own subfamily, blocks classic BMP signaling in multiple contexts, interacts with BMP proteins and is expressed specifically in the node during gastrulation in a pattern consistent with BMP inhibition. Furthermore, we use gain- and reduction-of-function to show that in a species-specific manner, GDF3 regulates both of the two major characteristics of embryonic stem cells: the ability to maintain the undifferentiated state and the ability to differentiate into the full spectrum of cell types.

Mechanical response of semiflexible networks to localized perturbations

Previous research on semiflexible polymers including cytoskeletal networks in cells has suggested the existence of distinct regimes of elastic response, in which the strain field is either uniform (affine) or nonuniform (nonaffine) under external stress. Associated with these regimes, it has been further suggested that a mesoscopic length scale emerges, which characterizes the scale for the crossover from nonaffine to affine deformations. Here, we extend these studies by probing the response to localized forces and force dipoles. We show that the previously identified nonaffinity length [D. A. Head, Phys. Rev. E 68, 061907 (2003)] controls the mesoscopic response to point forces and the crossover to continuum elastic behavior at large distances.

p53 tumor suppressor protein regulates the levels of huntingtin gene expression

The p53 protein is a transcription factor that integrates various cellular stress signals. The accumulation of the mutant huntingtin protein with an expanded polyglutamine tract plays a central role in the pathology of human Huntington's disease. We found that the huntingtin gene contains multiple putative p53-responsive elements and p53 binds to these elements both in vivo and in vitro. p53 activation in cultured human cells, either by a temperature-sensitive mutant p53 protein or by gamma-irradiation (gamma-irradiation), increases huntingtin mRNA and protein expression. Similarly, murine huntingtin also contains multiple putative p53-responsive elements and its expression is induced by p53 activation in cultured cells. Moreover, gamma-irradiation, which activates p53, increases huntingtin gene expression in the striatum and cortex of mouse brain, the major pathological sites for Huntington's disease, in p53+/+ but not the isogenic p53-/- mice. These results demonstrate that p53 protein can regulate huntingtin expression at transcriptional level, and suggest that a p53 stress response could be a modulator of the process of Huntington's disease.

TGFbeta/activin/nodal signaling is necessary for the maintenance of pluripotency in human embryonic stem cells

Human embryonic stem cells (hESCs) self-renew indefinitely and give rise to derivatives of all three primary germ layers, yet little is known about the signaling cascades that govern their pluripotent character. Because it plays a prominent role in the early cell fate decisions of embryonic development, we have examined the role of TGFbeta superfamily signaling in hESCs. We found that, in undifferentiated cells, the TGFbeta/activin/nodal branch is activated (through the signal transducer SMAD2/3) while the BMP/GDF branch (SMAD1/5) is only active in isolated mitotic cells. Upon early differentiation, SMAD2/3 signaling is decreased while SMAD1/5 signaling is activated. We next tested the functional role of TGFbeta/activin/nodal signaling in hESCs and found that it is required for the maintenance of markers of the undifferentiated state. We extend these findings to show that SMAD2/3 activation is required downstream of WNT signaling, which we have previously shown to be sufficient to maintain the undifferentiated state of hESCs. Strikingly, we show that in ex vivo mouse blastocyst cultures, SMAD2/3 signaling is also required to maintain the inner cell mass (from which stem cells are derived). These data reveal a crucial role for TGFbeta signaling in the earliest stages of cell fate determination and demonstrate an interconnection between TGFbeta and WNT signaling in these contexts.

Single-nucleotide polymorphisms in the p53 pathway

A cell culture assay has been developed that detects and validates single-nucleotide polymorphisms (SNPs) in genes that populate the p53 pathway. One hundred thirteen EBV-transformed human B-lymphocyte cell lines obtained from a diverse population were employed to measure the apoptotic response to gamma radiation. Each cell line undergoes a reproducible, characteristic frequency of apoptosis, and the response of the population forms a normal distribution around a median of 35.5% apoptosis with a range from 12% to 58% apoptosis. Polymorphisms in the AKT1 and Perp genes significantly affect the frequency of apoptosis. The assay can detect both racial and sexual dimorphisms in these genes and has the ability to demonstrate epistatic relationships within the p53 pathway. The cell lines used in this assay provide biological materials to explore the molecular basis of the polymorphisms.

Distinct regimes of elastic response and deformation modes of cross-linked cytoskeletal and semiflexible polymer networks

Semiflexible polymers such as filamentous actin (F-actin) play a vital role in the mechanical behavior of cells, yet the basic properties of cross-linked F-actin networks remain poorly understood. To address this issue, we have performed numerical studies of the linear response of homogeneous and isotropic two-dimensional networks subject to an applied strain at zero temperature. The elastic moduli are found to vanish for network densities at a rigidity percolation threshold. For higher densities, two regimes are observed: one in which the deformation is predominately affine and the filaments stretch and compress; and a second in which bending modes dominate. We identify a dimensionless scalar quantity, being a combination of the material length scales, that specifies to which regime a given network belongs. A scaling argument is presented that approximately agrees with this crossover variable. By a direct geometric measure, we also confirm that the degree of affinity under strain correlates with the distinct elastic regimes. We discuss the implications of our findings and suggest possible directions for future investigations.

Derivation and validation of a clinical scoring system to predict the need for an intra-aortic balloon pump in patients undergoing adult cardiac surgery

The spectrum of patients receiving cardiac surgery are increasing in age and severity of illness. With the reduction of complications caused by the placement of an intra-aortic balloon pump (IABP) there is increasing interest in the placement of an IABP prophylactically. We sought to derive a scoring system to guide the placement of IABPs. A total of 3927 patients from the Blackpool Victoria Open Heart Registry were used to derive a range of clinical decision scores using a range of established and novel statistical techniques. This database included 127 patients who received an IABP. The derived scores and rules were then validated on the North Staffordshire Open Heart Registry, containing 3070 patients, and 161 patients who received an IABP. We derived and validated a clinical score that has a sensitivity of 50% and a specificity of 96.5% in the prediction of those patients requiring an IABP. This was robust in the validation dataset and outperformed the Parsonnet score in this context. Our validated clinical scoring system will be useful both to guide individual clinical decision making and to compare variation of IABP usage among institutions.

Nonuniversality of elastic exponents in random bond-bending networks

We numerically investigate the rigidity percolation transition in two-dimensional flexible, random rod networks with freely rotating cross links. Near the transition, networks are dominated by bending modes and the elastic modulii vanish with an exponent f=3.0+/-0.2, in contrast with central force percolation which shares the same geometric exponents. This indicates that universality for geometric quantities does not imply universality for elastic ones. The implications of this result for actin-fiber networks is discussed.

Rheological microscopy: local mechanical properties from microrheology

We demonstrate how tracer microrheology methods can be extended to study submicron scale variations in the viscoelastic response of soft materials; in particular, a semidilute solution of lambda-DNA. The polymer concentration is depleted near the surfaces of the tracer particles, within a distance comparable to the polymer correlation length. The rheology of this microscopic layer alters the tracers' motion and can be precisely quantified using one- and two-point microrheology. Interestingly, we found this mechanically distinct layer to be twice as thick as the layer of depleted concentration, likely due to solvent drainage through the locally perturbed polymer structure.

Fluorescent labeling of endothelial cells allows in vivo, continuous characterization of the vascular development of Xenopus laevis

Appropriate blood supply and vascular development are necessary in development and in cancer, heart disease, and diabetes. Here, we report the use of DiI-labeled acetylated low-density lipoprotein (DiI-Ac-LDL) to label endothelial cells and characterize the vasculature of live Xenopus embryos. The atlas we have created provides a detailed map of normal vascular development against which perturbations of normal patterning can be compared. By following the development of the intersomitic vessels in real-time, we show that, while rostrocaudal gradient of maturing intersomitic vessels occurs, it is not absolute. In addition, the comparative study of the ontogeny of nerve bundles from the spinal cord of transgenic Xenopus embryos expressing green fluorescent protein in the nervous system and blood vessels demonstrates a strong anatomical correlation in neurovascular development. These studies provide the basis for understanding how the vascular system forms and assumes its complicated stereotypical pattern in normal development and in disease.

The p53MH algorithm and its application in detecting p53-responsive genes

A computer algorithm, p53MH, was developed, which identifies putative p53 transcription factor DNA-binding sites on a genomewide scale with high power and versatility. With the sequences from the human and mouse genomes, putative p53 DNA-binding elements were identified in a scan of 2,583 human genes and 1,713 mouse orthologs based on the experimental data of el-Deiry et al. [el-Deiry, W. S., Kern, S. E., Pietenpol, J. A., Kinzler, K. W. & Vogelstein, B. (1992) Nat. Genet. 1, 45-49] and Funk et al. [Funk, W. D., Pak, D. T., Karas, R. H., Wright, W. E. & Shay, J. W. (1992) Mol. Cell. Biol. 12, 2866-2871] (http://linkage.rockefeller.edu/p53). The p53 DNA-binding motif consists of a 10-bp palindrome and most commonly a second related palindrome linked by a spacer region. By scanning from the 5' to 3' end of each gene with an additional 10-kb nucleotide sequence appended at each end (most regulatory DNA elements characterized in the literature are in these regions), p53MH computes the binding likelihood for each site under a discrete discriminant model and then outputs ordered scores, corresponding site positions, sequences, and related information. About 300 genes receiving scores greater than a theoretical cut-off value were identified as potential p53 targets. Semiquantitative reverse transcription-PCR experiments were performed in 2 cell lines on 16 genes that were previously unknown regarding their functional relationship to p53 and were found to have high scores in either proximal promoter or possible distal enhancer regions. Ten (approximately 63%) of these genes responded to the presence of p53.

Serum 25-hydroxyvitamin D, dietary calcium intake, and distal colorectal adenoma risk

Vitamin D has recently emerged as a potentially protective agent against colorectal neoplasia. We assessed the associations between dietary vitamin D, plasma 25-hydroxyvitamin D [25(OH)D], dietary calcium, and colorectal adenomas in a large screening sigmoidoscopy-based case-control study in Southern California. Because conversion of serum 25(OH)D to serum 1,25-vitamin D is highly regulated by serum calcium, we also assessed modification of the 25(OH)D-adenoma association by calcium intake. Cases were 473 subjects with a primary adenoma, and controls were 507 subjects who had no adenomas at sigmoidoscopy and no history of adenomas. Compared with those in the lowest quartile of intake, those in the highest quartile of dietary vitamin D had an adjusted odds ratio (OR) of 0.83 [95% confidence interval (CI) = 0.49-1.41] and those in the highest quartile of dietary calcium had an OR of 0.82 (95% CI = 0.49-1.25). There was a suggestion that plasma 25(OH)D may be protective in this population (OR for highest vs. lowest quartile = 0.74, 95% CI = 0.51-1.09). A significant protective effect of 25(OH)D was clearly evident only in those with calcium intakes below (OR = 0.40 for highest vs. lowest quartile, 95% CI = 0.22-0.71, p for trend = 0.005) and above (OR = 1.17, 95% CI = 0.69-1.99, p for trend = 0.94) the median calcium intake.