CRISPRi screen of long non-coding RNAs identifies LINC03045 regulating glioblastoma invasion.

INTRODUCTION: Glioblastoma (GBM) invasion studies have focused on coding genes, while few studies evaluate long non-coding RNAs (lncRNAs), transcripts without protein-coding potential, for role in GBM invasion. We leveraged CRISPR-interference (CRISPRi) to evaluate invasive function of GBM-associated lncRNAs in an unbiased functional screen, characterizing and exploring the mechanism of identified candidates. METHODS: We implemented a CRISPRi lncRNA loss-of-function screen evaluating association of lncRNA knockdown (KD) with invasion capacity in Matrigel. Top screen candidates were validated using CRISPRi and oligonucleotide(ASO)-mediated knockdown in three tumor lines. Clinical relevance of candidates was assessed via The Cancer Genome Atlas(TCGA) and Genotype-Tissue Expression(GTEx) survival analysis. Mediators of lncRNA effect were identified via differential expression analysis following lncRNA KD and assessed for tumor invasion using knockdown and rescue experiments. RESULTS: Forty-eight lncRNAs were significantly associated with 33-83% decrease in invasion (p<0.01) upon knockdown. The top candidate, LINC03045, identified from effect size and p-value, demonstrated 82.7% decrease in tumor cell invasion upon knockdown, while LINC03045 expression was significantly associated with patient survival and tumor grade(p<0.0001). RNAseq analysis of LINC03045 knockdown revealed that WASF3, previously implicated in tumor invasion studies, was highly correlated with lncRNA expression, while WASF3 KD was associated with significant decrease in invasion. Finally, WASF3 overexpression demonstrated rescue of invasive function lost with LINC03045 KD. CONCLUSION: CRISPRi screening identified LINC03045, a previously unannotated lncRNA, as critical to GBM invasion. Gene expression is significantly associated with tumor grade and survival. RNA-seq and mechanistic studies suggest that this novel lncRNA may regulate invasion via WASF3.

The complete sequence and comparative analysis of ape sex chromosomes.

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes.

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Human PD-1 agonist treatment alleviates neutrophilic asthma by reprogramming T cells.

BACKGROUND: Neutrophilic asthma is associated with disease severity and corticosteroid insensitivity. Novel therapies are required to manage this life-threatening asthma phenotype. Programmed cell death protein-1 (PD-1) is a key homeostatic modulator of the immune response for T-cell effector functions. OBJECTIVE: We sought to investigate the role of PD-1 in the regulation of acute neutrophilic inflammation in a murine model of airway hyperreactivity (AHR). METHODS: House dust mite was used to induce and compare neutrophilic AHR in wild-type and PD-1 knockout mice. Then, the therapeutic potential of a human PD-1 agonist was tested in a humanized mouse model in which the PD-1 extracellular domain is entirely humanized. Single-cell RNA sequencing and flow cytometry were mainly used to investigate molecular and cellular mechanisms. RESULTS: PD-1 was highly induced on pulmonary T cells in our inflammatory model. PD-1 deficiency was associated with an increased neutrophilic AHR and high recruitment of inflammatory cells to the lungs. Consistently, PD-1 agonist treatment dampened AHR, decreased neutrophil recruitment, and modulated cytokine production in a humanized PD-1 mouse model. Mechanistically, we demonstrated at the transcriptional and protein levels that the inhibitory effect of PD-1 agonist is associated with the reprogramming of pulmonary effector T cells that showed decreased number and activation. CONCLUSIONS: PD-1 agonist treatment is efficient in dampening neutrophilic AHR and lung inflammation in a preclinical humanized mouse model.

Decoding the IGF1 signaling gene regulatory network behind alveologenesis from a mouse model of bronchopulmonary dysplasia.

Lung development is precisely controlled by underlying gene regulatory networks (GRN). Disruption of genes in the network can interrupt normal development and cause diseases such as bronchopulmonary dysplasia (BPD) - a chronic lung disease in preterm infants with morbid and sometimes lethal consequences characterized by lung immaturity and reduced alveolarization. Here, we generated a transgenic mouse exhibiting a moderate severity BPD phenotype by blocking IGF1 signaling in secondary crest myofibroblasts (SCMF) at the onset of alveologenesis. Using approaches mirroring the construction of the model GRN in sea urchin's development, we constructed the IGF1 signaling network underlying alveologenesis using this mouse model that phenocopies BPD. The constructed GRN, consisting of 43 genes, provides a bird's eye view of how the genes downstream of IGF1 are regulatorily connected. The GRN also reveals a mechanistic interpretation of how the effects of IGF1 signaling are transduced within SCMF from its specification genes to its effector genes and then from SCMF to its neighboring alveolar epithelial cells with WNT5A and FGF10 signaling as the bridge. Consistently, blocking WNT5A signaling in mice phenocopies BPD as inferred by the network. A comparative study on human samples suggests that a GRN of similar components and wiring underlies human BPD. Our network view of alveologenesis is transforming our perspective to understand and treat BPD. This new perspective calls for the construction of the full signaling GRN underlying alveologenesis, upon which targeted therapies for this neonatal chronic lung disease can be viably developed.

Neuron-specific chromosomal megadomain organization is adaptive to recent retrotransposon expansions.

Regulatory mechanisms associated with repeat-rich sequences and chromosomal conformations in mature neurons remain unexplored. Here, we map cell-type specific chromatin domain organization in adult mouse cerebral cortex and report strong enrichment of Endogenous Retrovirus 2 (ERV2) repeat sequences in the neuron-specific heterochromatic B2NeuN+ megabase-scaling subcompartment. Single molecule long-read sequencing and comparative Hi-C chromosomal contact mapping in wild-derived SPRET/EiJ (Mus spretus) and laboratory inbred C57BL/6J (Mus musculus) reveal neuronal reconfigurations tracking recent ERV2 expansions in the murine germline, with significantly higher B2NeuN+ contact frequencies at sites with ongoing insertions in Mus musculus. Neuronal ablation of the retrotransposon silencer Kmt1e/Setdb1 triggers B2NeuN+ disintegration and rewiring with open chromatin domains enriched for cellular stress response genes, along with severe neuroinflammation and proviral assembly with infiltration of dendrites . We conclude that neuronal megabase-scale chromosomal architectures include an evolutionarily adaptive heterochromatic organization which, upon perturbation, results in transcriptional dysregulation and unleashes ERV2 proviruses with strong neuronal tropism.

Preclinical Evaluation of a Novel Dual Targeting PI3Kδ/BRD4 Inhibitor, SF2535, in B-Cell Acute Lymphoblastic Leukemia.

The PI3K/Akt pathway-and in particular PI3Kδ-is known for its role in drug resistant B-cell acute lymphoblastic leukemia (B-ALL) and it is often upregulated in refractory or relapsed B-ALL. Myc proteins are transcription factors responsible for transcribing pro-proliferative genes and c-Myc is often overexpressed in cancers. The chromatin regulator BRD4 is required for expression of c-Myc in hematologic malignancies including B-ALL. Previously, combination of BRD4 and PI3K inhibition with SF2523 was shown to successfully decrease Myc expression. However, the underlying mechanism and effect of dual inhibition of PI3Kδ/BRD4 in B-ALL remains unknown. To study this, we utilized SF2535, a novel small molecule dual inhibitor which can specifically target the PI3Kδ isoform and BRD4. We treated primary B-ALL cells with various concentrations of SF2535 and studied its effect on specific pharmacological on-target mechanisms such as apoptosis, cell cycle, cell proliferation, and adhesion molecules expression usingin vitro and in vivo models. SF2535 significantly downregulates both c-Myc mRNA and protein expression through inhibition of BRD4 at the c-Myc promoter site and decreases p-AKT expression through inhibition of the PI3Kδ/AKT pathway. SF2535 induced apoptosis in B-ALL by downregulation of BCL-2 and increased cleavage of caspase-3, caspase-7, and PARP. Moreover, SF2535 induced cell cycle arrest and decreased cell counts in B-ALL. Interestingly, SF2535 decreased the mean fluorescence intensity (MFI) of integrin α4, α5, α6, and ß1 while increasing MFI of CXCR4, indicating that SF2535 may work through inside-out signaling of integrins. Taken together, our data provide a rationale for the clinical evaluation of targeting PI3Kδ/BRD4 in refractory or relapsed B-ALL using SF2535.

Reconstituting the genus Mycobacterium.

The definition of a genus has wide-ranging implications both in terms of binomial species names and also evolutionary relationships. In recent years, the definition of the genus Mycobacterium has been debated due to the proposed split of this genus into five new genera (Mycolicibacterium, Mycolicibacter, Mycolicibacillus, Mycobacteroides and an emended Mycobacterium). Since this group of species contains many important obligate and opportunistic pathogens, it is important that any renaming of species does not cause confusion in clinical treatment as outlined by the nomen periculosum rule (56a) of the Prokaryotic Code. In this study, we evaluated the proposed and original genus boundaries for the mycobacteria, to determine if the split into five genera was warranted. By combining multiple approaches for defining genus boundaries (16S rRNA gene similarity, amino acid identity index, average nucleotide identity, alignment fraction and percentage of conserved proteins) we show that the original genus Mycobacterium is strongly supported over the proposed five-way split. Thus, we propose that the original genus label be reapplied to all species within this group, with the proposed five genera potentially used as sub-genus complex names.

Creation of a Single Cell RNASeq Meta-Atlas to Define Human Liver Immune Homeostasis.

The liver is unique in both its ability to maintain immune homeostasis and in its potential for immune tolerance following solid organ transplantation. Single-cell RNA sequencing (scRNA seq) is a powerful approach to generate highly dimensional transcriptome data to understand cellular phenotypes. However, when scRNA data is produced by different groups, with different data models, different standards, and samples processed in different ways, it can be challenging to draw meaningful conclusions from the aggregated data. The goal of this study was to establish a method to combine 'human liver' scRNA seq datasets by 1) characterizing the heterogeneity between studies and 2) using the meta-atlas to define the dominant phenotypes across immune cell subpopulations in healthy human liver. Publicly available scRNA seq data generated from liver samples obtained from a combined total of 17 patients and ~32,000 cells were analyzed. Liver-specific immune cells (CD45+) were extracted from each dataset, and immune cell subpopulations (myeloid cells, NK and T cells, plasma cells, and B cells) were examined using dimensionality reduction (UMAP), differential gene expression, and ingenuity pathway analysis. All datasets co-clustered, but cell proportions differed between studies. Gene expression correlation demonstrated similarity across all studies, and canonical pathways that differed between datasets were related to cell stress and oxidative phosphorylation rather than immune-related function. Next, a meta-atlas was generated via data integration and compared against PBMC data to define gene signatures for each hepatic immune subpopulation. This analysis defined key features of hepatic immune homeostasis, with decreased expression across immunologic pathways and enhancement of pathways involved with cell death. This method for meta-analysis of scRNA seq data provides a novel approach to broadly define the features of human liver immune homeostasis. Specific pathways and cellular phenotypes described in this human liver immune meta-atlas provide a critical reference point for further study of immune mediated disease processes within the liver.

Mechanistic basis for chromosomal translocations at the E2A gene and its broader relevance to human B cell malignancies.

Analysis of translocation breakpoints in human B cell malignancies reveals that DNA double-strand breaks at oncogenes most frequently occur at CpG sites located within 20-600 bp fragile zones and depend on activation-induced deaminase (AID). AID requires single-stranded DNA (ssDNA) to act, but it has been unclear why or how this region transiently acquires a ssDNA state. Here, we demonstrate the ssDNA state in the 23 bp E2A fragile zone using several methods, including native bisulfite DNA structural analysis in live human pre-B cells. AID deamination within the E2A fragile zone does not require but is increased upon transcription. High C-string density, nascent RNA tails, and direct DNA sequence repeats prolong the ssDNA state of the E2A fragile zone and increase AID deamination at overlapping AID hotspots that contain the CpG sites at which breaks occur in patients. These features provide key insights into lymphoid fragile zones generally.

Vitamin D deficiency exacerbates UV/endorphin and opioid addiction.

The current opioid epidemic warrants a better understanding of genetic and environmental factors that contribute to opioid addiction. Here we report an increased prevalence of vitamin D (VitD) deficiency in patients diagnosed with opioid use disorder and an inverse and dose-dependent association of VitD levels with self-reported opioid use. We used multiple pharmacologic approaches and genetic mouse models and found that deficiencies in VitD signaling amplify exogenous opioid responses that are normalized upon restoration of VitD signaling. Similarly, physiologic endogenous opioid analgesia and reward responses triggered by ultraviolet (UV) radiation are repressed by VitD signaling, suggesting that a feedback loop exists whereby VitD deficiency produces increased UV/endorphin-seeking behavior until VitD levels are restored by cutaneous VitD synthesis. This feedback may carry the evolutionary advantage of maximizing VitD synthesis. However, unlike UV exposure, exogenous opioid use is not followed by VitD synthesis (and its opioid suppressive effects), contributing to maladaptive addictive behavior.

CD200-CD200R immune checkpoint engagement regulates ILC2 effector function and ameliorates lung inflammation in asthma.

The prevalence of asthma and airway hyperreactivity (AHR) is increasing at an alarming rate. Group 2 innate lymphoid cells (ILC2s) are copious producers of type 2 cytokines, which leads to AHR and lung inflammation. Here, we show that mouse ILC2s express CD200 receptor (CD200R) and this expression is inducible. CD200R engagement inhibits activation, proliferation and type 2 cytokine production, indicating an immunoregulatory function for the CD200-CD200R axis on ILC2s. Furthermore, CD200R engagement inhibits both canonical and non-canonical NF-κB signaling pathways in activated ILC2s. Additionally, we demonstrate both preventative and therapeutic approaches utilizing CD200R engagement on ILC2s, which lead to improved airway resistance, dynamic compliance and eosinophilia. These results show CD200R is expressed on human ILC2s, and its engagement ameliorates AHR in humanized mouse models, emphasizing the translational applications for treatment of ILC2-related diseases such as allergic asthma.

In vivo CRISPR screening for novel noncoding RNA functional targets in glioblastoma models.

CRISPR (clustered regularly interspaced short palindromic repeat)-based genetic screens offer unbiased and powerful tools for systematic and specific evaluation of phenotypes associated with specific target genes. CRISPR screens have been utilized heavily in vitro to identify functional coding and noncoding genes in a large number of cell types, including glioblastoma (GB), though no prior study has described the evaluation of CRISPR screening in GB in vivo. Here, we describe a protocol for targeting and transcriptionally repressing GB-specific long noncoding RNAs (lncRNAs) by CRISPR interference (CRISPRi) system in vivo, with tumor growth in the mouse cerebral cortex. Given the target-specific parameters of each individual screen, we list general steps involved in transducing guide RNA libraries into GB tumor lines, maintaining sufficient coverage, as well as cortically injecting and subsequently isolating transduced screen tumor cell populations for analysis. Finally, in order to demonstrate the use of this technique to discern an essential lncRNA, HOTAIR, from a nonessential lncRNA, we injected a 1:1 (HOTAIR:control nonessential lncRNA knockdown) mixture of fluorescently tagged U87 GB cells into the cortex of eight mice, evaluating selective depletion of HOTAIR-tagged cells at 2 weeks of growth. Fluorescently tagged populations were analyzed via flow cytometry for hiBFP (control knockdown) and green fluorescent protein (HOTAIR knockdown), revealing 17% (p = 0.007) decrease in fluorescence associated with HOTAIR knockdown relative to control. The described in vivo CRISPR screening methodology thus appears to be an effective option for identifying noncoding (and coding) genes affecting GB growth within the mouse cortex.

Runx2-Twist1 interaction coordinates cranial neural crest guidance of soft palate myogenesis.

Cranial neural crest (CNC) cells give rise to bone, cartilage, tendons, and ligaments of the vertebrate craniofacial musculoskeletal complex, as well as regulate mesoderm-derived craniofacial muscle development through cell-cell interactions. Using the mouse soft palate as a model, we performed an unbiased single-cell RNA-seq analysis to investigate the heterogeneity and lineage commitment of CNC derivatives during craniofacial muscle development. We show that Runx2, a known osteogenic regulator, is expressed in the CNC-derived perimysial and progenitor populations. Loss of Runx2 in CNC-derivatives results in reduced expression of perimysial markers (Aldh1a2 and Hic1) as well as soft palate muscle defects in Osr2-Cre;Runx2fl/fl mice. We further reveal that Runx2 maintains perimysial marker expression through suppressing Twist1, and that myogenesis is restored in Osr2-Cre;Runx2fl/fl;Twist1fl/+ mice. Collectively, our findings highlight the roles of Runx2, Twist1, and their interaction in regulating the fate of CNC-derived cells as they guide craniofacial muscle development through cell-cell interactions.

Spatiotemporal cellular movement and fate decisions during first pharyngeal arch morphogenesis.

Cranial neural crest (CNC) cells contribute to different cell types during embryonic development. It is unknown whether postmigratory CNC cells undergo dynamic cellular movement and how the process of cell fate decision occurs within the first pharyngeal arch (FPA). Our investigations demonstrate notable heterogeneity within the CNC cells, refine the patterning domains, and identify progenitor cells within the FPA. These progenitor cells undergo fate bifurcation that separates them into common progenitors and mesenchymal cells, which are characterized by Cdk1 and Spry2/Notch2 expression, respectively. The common progenitors undergo further bifurcations to restrict them into osteogenic/odontogenic and chondrogenic/fibroblast lineages. Disruption of a patterning domain leads to specific mandible and tooth defects, validating the binary cell fate restriction process. Different from the compartment model of mandibular morphogenesis, our data redefine heterogeneous cellular domains within the FPA, reveal dynamic cellular movement in time, and describe a sequential series of binary cell fate decision-making process.

Negative feedback control of neuronal activity by microglia.

Microglia, the brain's resident macrophages, help to regulate brain function by removing dying neurons, pruning non-functional synapses, and producing ligands that support neuronal survival1. Here we show that microglia are also critical modulators of neuronal activity and associated behavioural responses in mice. Microglia respond to neuronal activation by suppressing neuronal activity, and ablation of microglia amplifies and synchronizes the activity of neurons, leading to seizures. Suppression of neuronal activation by microglia occurs in a highly region-specific fashion and depends on the ability of microglia to sense and catabolize extracellular ATP, which is released upon neuronal activation by neurons and astrocytes. ATP triggers the recruitment of microglial protrusions and is converted by the microglial ATP/ADP hydrolysing ectoenzyme CD39 into AMP; AMP is then converted into adenosine by CD73, which is expressed on microglia as well as other brain cells. Microglial sensing of ATP, the ensuing microglia-dependent production of adenosine, and the adenosine-mediated suppression of neuronal responses via the adenosine receptor A1R are essential for the regulation of neuronal activity and animal behaviour. Our findings suggest that this microglia-driven negative feedback mechanism operates similarly to inhibitory neurons and is essential for protecting the brain from excessive activation in health and disease.

Single-Cell RNA-Seq Mapping of Human Thymopoiesis Reveals Lineage Specification Trajectories and a Commitment Spectrum in T Cell Development.

The challenges in recapitulating in vivo human T cell development in laboratory models have posed a barrier to understanding human thymopoiesis. Here, we used single-cell RNA sequencing (sRNA-seq) to interrogate the rare CD34+ progenitor and the more differentiated CD34- fractions in the human postnatal thymus. CD34+ thymic progenitors were comprised of a spectrum of specification and commitment states characterized by multilineage priming followed by gradual T cell commitment. The earliest progenitors in the differentiation trajectory were CD7- and expressed a stem-cell-like transcriptional profile, but had also initiated T cell priming. Clustering analysis identified a CD34+ subpopulation primed for the plasmacytoid dendritic lineage, suggesting an intrathymic dendritic specification pathway. CD2 expression defined T cell commitment stages where loss of B cell potential preceded that of myeloid potential. These datasets delineate gene expression profiles spanning key differentiation events in human thymopoiesis and provide a resource for the further study of human T cell development.

Molecular mechanisms of Guadecitabine induced FGFR4 down regulation in alveolar rhabdomyosarcomas.

Fibroblast growth factor receptor 4 (FGFR4) aberrant expression and activity have been linked to the pathogenesis of a variety of cancers including rhabdomyosarcomas (RMS). We found that treatment of alveolar rhabdomyosarcoma (aRMS) cells with Guadecitabine (SGI-110), a next-generation DNA methyltransferase inhibitor (DNMTi), resulted in a significant reduction of FGFR4 protein levels, 5 days post treatment. Chromatin immunoprecipitation-sequencing (ChIP-seq) in aRMS cells revealed attenuation of the H3K4 mono-methylation across the FGFR4 super enhancer without changes in tri-methylation of either H3K4 or H3K27. These changes were associated with a significant reduction in FGFR4 transcript levels in treated cells. These decreases in H3K4me1 in the FGFR4 super enhancer were also associated with a 240-fold increase in KDM5B (JARID1B) mRNA levels. Immunoblot and immunofluorescent studies also revealed a significant increase in the KDM5B protein levels after treatment in these cells. KDM5B is the only member of KDM5 (JARID1) family of histone lysine demethylases that catalyzes demethylation of H3K4me1. These data together suggest a pleiotropic effect of DNMTi therapy in aRMS cells, converging to significantly lower FGFR4 protein levels in these cells.

Sex-Specific Role for the Long Non-coding RNA LINC00473 in Depression.

Depression is a common disorder that affects women at twice the rate of men. Here, we report that long non-coding RNAs (lncRNAs), a recently discovered class of regulatory transcripts, represent about one-third of the differentially expressed genes in the brains of depressed humans and display complex region- and sex-specific patterns of regulation. We identified the primate-specific, neuronal-enriched gene LINC00473 as downregulated in prefrontal cortex (PFC) of depressed females but not males. Using viral-mediated gene transfer to express LINC00473 in adult mouse PFC neurons, we mirrored the human sex-specific phenotype by inducing stress resilience solely in female mice. This sex-specific phenotype was accompanied by changes in synaptic function and gene expression selectively in female mice and, along with studies of human neuron-like cells in culture, implicates LINC00473 as a CREB effector. Together, our studies identify LINC00473 as a female-specific driver of stress resilience that is aberrant in female depression.

Shared Transcriptional Signatures in Major Depressive Disorder and Mouse Chronic Stress Models.

BACKGROUND: Most of our knowledge of the biological basis of major depressive disorder (MDD) is derived from studies of chronic stress models in rodents. While these models capture certain aspects of the behavioral and neuroendocrine features of MDD, the extent to which they reproduce the molecular pathology of the human syndrome remains unknown. METHODS: We systematically compared transcriptional signatures in two brain regions implicated in depression-medial prefrontal cortex and nucleus accumbens-of humans with MDD and of 3 chronic stress models in mice: chronic variable stress, adult social isolation, and chronic social defeat stress. We used differential expression analysis combined with weighted gene coexpression network analysis to create interspecies gene networks and assess the capacity of each stress paradigm to recapitulate the transcriptional organization of gene networks in human MDD. RESULTS: Our results show significant overlap between transcriptional alterations in medial prefrontal cortex and nucleus accumbens in human MDD and the 3 mouse chronic stress models, with each of the chronic stress paradigms capturing distinct aspects of MDD abnormalities. Chronic variable stress and adult social isolation better reproduce differentially expressed genes, while chronic social defeat stress and adult social isolation better reproduce gene networks characteristic of human MDD. We also identified several gene networks and their constituent genes that are most significantly associated with human MDD and mouse stress models. CONCLUSIONS: This study demonstrates the ability of 3 chronic stress models in mice to recapitulate distinct aspects of the broad molecular pathology of human MDD, with no one mouse model apparently better than another.