Evolution of gene expression across brain regions in behaviourally divergent deer mice.

The evolution of innate behaviours is ultimately due to genetic variation likely acting in the nervous system. Gene regulation may be particularly important because it can evolve in a modular brain-region specific fashion through the concerted action of cis- and trans-regulatory changes. Here, to investigate transcriptional variation and its regulatory basis across the brain, we perform RNA sequencing (RNA-Seq) on ten brain subregions in two sister species of deer mice (Peromyscus maniculatus and P. polionotus)-which differ in a range of innate behaviours, including their social system-and their F1 hybrids. We find that most of the variation in gene expression distinguishes subregions, followed by species. Interspecific differential expression (DE) is pervasive (52-59% of expressed genes), whereas the number of DE genes between sexes is modest overall (~3%). Interestingly, the identity of DE genes varies considerably across brain regions. Much of this modularity is due to cis-regulatory divergence, and while 43% of genes were consistently assigned to the same gene regulatory class across subregions (e.g. conserved, cis- or trans-regulatory divergence), a similar number were assigned to two or more different gene regulatory classes. Together, these results highlight the modularity of gene expression differences and divergence in the brain, which may be key to explain how the evolution of brain gene expression can contribute to the astonishing diversity of animal behaviours.

The genetic basis of parental care evolution in monogamous mice.

Parental care is essential for the survival of mammals, yet the mechanisms underlying its evolution remain largely unknown. Here we show that two sister species of mice, Peromyscus polionotus and Peromyscus maniculatus, have large and heritable differences in parental behaviour. Using quantitative genetics, we identify 12 genomic regions that affect parental care, 8 of which have sex-specific effects, suggesting that parental care can evolve independently in males and females. Furthermore, some regions affect parental care broadly, whereas others affect specific behaviours, such as nest building. Of the genes linked to differences in nest-building behaviour, vasopressin is differentially expressed in the hypothalamus of the two species, with increased levels associated with less nest building. Using pharmacology in Peromyscus and chemogenetics in Mus, we show that vasopressin inhibits nest building but not other parental behaviours. Together, our results indicate that variation in an ancient neuropeptide contributes to interspecific differences in parental care.

A phenotypic screening platform for chronic pain therapeutics using all-optical electrophysiology.

ABSTRACT: Chronic pain associated with osteoarthritis (OA) remains an intractable problem with few effective treatment options. New approaches are needed to model the disease biology and to drive discovery of therapeutics. We present an in vitro model of OA pain, where dorsal root ganglion (DRG) sensory neurons were sensitized by a defined mixture of disease-relevant inflammatory mediators, here called Sensitizing PAin Reagent Composition or SPARC. Osteoarthritis-SPARC components showed synergistic or additive effects when applied in combination and induced pain phenotypes in vivo. To measure the effect of OA-SPARC on neural firing in a scalable format, we used a custom system for high throughput all-optical electrophysiology. This system enabled light-based membrane voltage recordings from hundreds of neurons in parallel with single cell and single action potential resolution and a throughput of up to 500,000 neurons per day. A computational framework was developed to construct a multiparameter OA-SPARC neuronal phenotype and to quantitatively assess phenotype reversal by candidate pharmacology. We screened ∼3000 approved drugs and mechanistically focused compounds, yielding data from over 1.2 million individual neurons with detailed assessment of functional OA-SPARC phenotype rescue and orthogonal "off-target" effects. Analysis of confirmed hits revealed diverse potential analgesic mechanisms including ion channel modulators and other mechanisms including MEK inhibitors and tyrosine kinase modulators. Our results suggest that the Raf-MEK-ERK axis in DRG neurons may integrate the inputs from multiple upstream inflammatory mediators found in osteoarthritis patient joints, and MAPK pathway activation in DRG neurons may contribute to chronic pain in patients with osteoarthritis.

Developing antisense oligonucleotides for a TECPR2 mutation-induced, ultra-rare neurological disorder using patient-derived cellular models.

Mutations in the TECPR2 gene are the cause of an ultra-rare neurological disorder characterized by intellectual disability, impaired speech, motor delay, and hypotonia evolving to spasticity, central sleep apnea, and premature death (SPG49 or HSAN9; OMIM: 615031). Little is known about the biological function of TECPR2, and there are currently no available disease-modifying therapies for this disease. Here we describe implementation of an antisense oligonucleotide (ASO) exon-skipping strategy targeting TECPR2 c.1319delT (p.Leu440Argfs∗19), a pathogenic variant that results in a premature stop codon within TECPR2 exon 8. We used patient-derived fibroblasts and induced pluripotent stem cell (iPSC)-derived neurons homozygous for the p.Leu440Argfs∗19 mutation to model the disease in vitro. Both patient-derived fibroblasts and neurons showed lack of TECPR2 protein expression. We designed and screened ASOs targeting sequences across the TECPR2 exon 8 region to identify molecules that induce exon 8 skipping and thereby remove the premature stop signal. TECPR2 exon 8 skipping restored in-frame expression of a TECPR2 protein variant (TECPR2ΔEx8) containing 1,300 of 1,411 amino acids. Optimization of ASO sequences generated a lead candidate (ASO-005-02) with ∼27 nM potency in patient-derived fibroblasts. To examine potential functional rescue induced by ASO-005-02, we used iPSC-derived neurons to analyze the neuronal localization of TECPR2ΔEx8 and showed that this form of TECPR2 retains the distinct, punctate neuronal expression pattern of full-length TECPR2. Finally, ASO-005-02 had an acceptable tolerability profile in vivo following a single 20-mg intrathecal dose in cynomolgus monkeys, showing some transient non-adverse behavioral effects with no correlating histopathology. Broad distribution of ASO-005-02 and induction of TECPR2 exon 8 skipping was detected in multiple central nervous system (CNS) tissues, supporting the potential utility of this therapeutic strategy for a subset of patients suffering from this rare disease.

Pol II docking and pausing at growth and stress genes in C. elegans.

Fluctuations in nutrient availability profoundly impact gene expression. Previous work revealed postrecruitment regulation of RNA polymerase II (Pol II) during starvation and recovery in Caenorhabditis elegans, suggesting that promoter-proximal pausing promotes rapid response to feeding. To test this hypothesis, we measured Pol II elongation genome wide by two complementary approaches and analyzed elongation in conjunction with Pol II binding and expression. We confirmed bona fide pausing during starvation and also discovered Pol II docking. Pausing occurs at active stress-response genes that become downregulated in response to feeding. In contrast, "docked" Pol II accumulates without initiating upstream of inactive growth genes that become rapidly upregulated upon feeding. Beyond differences in function and expression, these two sets of genes have different core promoter motifs, suggesting alternative transcriptional machinery. Our work suggests that growth and stress genes are both regulated postrecruitment during starvation but at initiation and elongation, respectively, coordinating gene expression with nutrient availability.