Search

Phenotypic plasticity, the ability to produce multiple phenotypes from a single genotype, represents an excellent model with which to examine the relationship between gene expression and phenotypes. Analyses of the molecular foundations of phenotypic plasticity are challenging, however, especially in the case of complex social phenotypes. Here we apply a machine learning approach to tackle this challenge by analyzing individual-level gene expression profiles of Polistes dominula paper wasps following the loss of a queen. We find that caste-associated gene expression profiles respond strongly to queen loss, and that this change is partly explained by attributes such as age but occurs even in individuals that appear phenotypically unaffected. These results demonstrate that large changes in gene expression may occur in the absence of outwardly detectable phenotypic changes, resulting here in a socially mediated de-differentiation of individuals at the transcriptomic level but not at the levels of ovarian development or behavior. Connecting genotypes to complex social behaviour is challenging. Taylor et al. use machine learning to show a strong response of caste-associated gene expression to queen loss, wherein individual wasp's expression profiles become intermediate between queen and worker states, even in the absence of behavioural changes. [ABSTRACT FROM AUTHOR]

Single-cell RNA sequencing (scRNA-seq) technologies have become essential tools for characterizing cellular landscapes within complex tissues. Large-scale single-cell transcriptomics holds great potential for identifying rare cell types critical to the pathogenesis of diseases and biological processes. Existing methods for identifying rare cell types often rely on one-time clustering using partial or global gene expression. However, these rare cell types may be overlooked during the clustering phase, posing challenges for their accurate identification. In this paper, we propose a Cluster decomposition-based Anomaly Detection method (scCAD), which iteratively decomposes clusters based on the most differential signals in each cluster to effectively separate rare cell types and achieve accurate identification. We benchmark scCAD on 25 real-world scRNA-seq datasets, demonstrating its superior performance compared to 10 state-of-the-art methods. In-depth case studies across diverse datasets, including mouse airway, brain, intestine, human pancreas, immunology data, and clear cell renal cell carcinoma, showcase scCAD's efficiency in identifying rare cell types in complex biological scenarios. Furthermore, scCAD can correct the annotation of rare cell types and identify immune cell subtypes associated with disease, thereby offering valuable insights into disease progression. Identifying rare cells is essential for advancing our understanding of complex biological systems and disease mechanisms. Here, authors propose scCAD, a method that combines cluster decomposition and anomaly detection to effectively identify rare cell types across diverse biological scenarios. [ABSTRACT FROM AUTHOR]

The goal of the National Cancer Institute's (NCI's) Genomic Data Commons (GDC) is to provide the cancer research community with a data repository of uniformly processed genomic and associated clinical data that enables data sharing and collaborative analysis in the support of precision medicine. The initial GDC dataset include genomic, epigenomic, proteomic, clinical and other data from the NCI TCGA and TARGET programs. Data production for the GDC started in June, 2015 using an OpenStack-based private cloud. By June of 2016, the GDC had analyzed more than 50,000 raw sequencing data inputs, as well as multiple other data types. Using the latest human genome reference build GRCh38, the GDC generated a variety of data types from aligned reads to somatic mutations, gene expression, miRNA expression, DNA methylation status, and copy number variation. In this paper, we describe the pipelines and workflows used to process and harmonize the data in the GDC. The generated data, as well as the original input files from TCGA and TARGET, are available for download and exploratory analysis at the GDC Data Portal and Legacy Archive (https://gdc.cancer.gov/). The Genomic Data Commons repository contains genomic, epigenomic, proteomic and clinical data from the TCGA and TARGET datasets. Here, the authors describe the analysis methods for how these divergent datasets were integrated together. [ABSTRACT FROM AUTHOR]

An amendment to this paper has been published and can be accessed via a link at the top of the paper. [ABSTRACT FROM AUTHOR]

An amendment to this paper has been published and can be accessed via a link at the top of the paper. [ABSTRACT FROM AUTHOR]

The neocortex varies in size and complexity among mammals due to the tremendous variability in the number and diversity of neuronal subtypes across species. The increased cellular diversity is paralleled by the expansion of the pool of neocortical progenitors and the emergence of indirect neurogenesis during brain evolution. The molecular pathways that control these biological processes and are disrupted in neurological disorders remain largely unknown. Here we show that the transcription factors BRN1 and BRN2 have an evolutionary conserved function in neocortical progenitors to control their proliferative capacity and the switch from direct to indirect neurogenesis. Functional studies in mice and ferrets show that BRN1/2 act in concert with NOTCH and primary microcephaly genes to regulate progenitor behavior. Analysis of transcriptomics data from genetically modified macaques provides evidence that these molecular pathways are conserved in non-human primates. Our findings thus demonstrate that BRN1/2 are central regulators of gene expression programs in neocortical progenitors critical to determine brain size during evolution. How the diversity and number of neuronal subtypes across species drives neocortical size and complexity remains incompletely understood. Here authors demonstrate BRN1/2 act in concert with NOTCH and primary microcephaly genes to regulate evolutionarily conserved neocortical progenitor behavior across species and determine brain size. [ABSTRACT FROM AUTHOR]

Molecular imaging holds the potential for noninvasive and accurate grading of liver fibrosis. It is limited by the lack of biomarkers that strongly correlate with liver fibrosis grade. Here, we discover the grading potential of fibroblast activation protein alpha (FAPα) for liver fibrosis through transcriptional analysis and biological assays on clinical liver samples. The protein and mRNA expression of FAPα are linearly correlated with fibrosis grade (R2 = 0.89 and 0.91, respectively). A FAPα-responsive MRI molecular nanoprobe is prepared for quantitatively grading liver fibrosis. The nanoprobe is composed of superparamagnetic amorphous iron nanoparticles (AFeNPs) and paramagnetic gadoteric acid (Gd-DOTA) connected by FAPα-responsive peptide chains (ASGPAGPA). As liver fibrosis worsens, the increased FAPα cut off more ASGPAGPA, restoring a higher T1-MRI signal of Gd-DOTA. Otherwise, the signal remains quenched due to the distance-dependent magnetic resonance tuning (MRET) effect between AFeNPs and Gd-DOTA. The nanoprobe identifies F1, F2, F3, and F4 fibrosis, with area under the curve of 99.8%, 66.7%, 70.4%, and 96.3% in patients' samples, respectively. This strategy exhibits potential in utilizing molecular imaging for the early detection and grading of liver fibrosis in the clinic. Molecular imaging holds promise for liver fibrosis grading but limited by a lack of effective biomarker. Here, the authors discover the grading potential of fibroblast activation protein alpha and design a responsive MRI nanoprobe that achieves fibrosis grading in mice models and clinical samples. [ABSTRACT FROM AUTHOR]

The identification effector targets and characterization of their functions are crucial for understanding pathogen infection mechanisms and components of plant immunity. Here, we identify the effector UgsL, a ustilaginoidin synthetase with a key role in regulating virulence of the rice false smut fungus Ustilaginoidea virens. Heterologous expression of UgsL in rice (Oryza sativa) enhances plant susceptibility to multiple pathogens, and host-induced gene silencing of UgsL enhances plant resistance to U. virens, indicating that UgsL inhibits rice immunity. UgsL interacts with STRUBBELIG RECEPTOR KINASE 3 (OsSRF3). Genome editing and overexpression of OsSRF3 demonstrate that OsSRF3 plays a pivotal role in the resistance of rice to multiple pathogens. Remarkably, overexpressing OsSRF3 enhances resistance without adversely affecting plant growth or yield. We show that BRASSINOSTEROID RECEPTOR-ASSOCIATED KINASE 1 (OsBAK1) interacts with and phosphorylates OsSRF3 to activate pathogen-triggered immunity, inducing the mitogen-activated protein kinase cascade, a reactive oxygen species burst, callose deposition, and expression of defense-related genes. UgsL interferes with the phosphorylation of OsSRF3 by OsBAK1. Furthermore, UgsL mediates OsSRF3 degradation by facilitating its association with the ubiquitin-26S proteasome. Our results reveal that OsSRF3 positively regulates immunity in rice and that UgsL mediates its degradation, thereby inhibiting the activation of OsBAK1–OsSRF3-mediated immune pathways. The Ustilaginoidea virens effector UgsL is shown to interact with receptor kinase OsSRF3. OsSRF3 positively regulates immunity in rice and UgsL mediates its degradation, thereby inhibiting the activation of OsBAK1–OsSRF3-mediated immune responses in rice. [ABSTRACT FROM AUTHOR]

Enhancing drought tolerance in crops and understanding the underlying mechanisms have been subject of intense research. The precise function and molecular mechanisms of B-box zinc finger proteins (BBX) remain elusive. Here, we report a natural allele of BBX18 (BBX18TT) that encodes a C-terminal truncated protein. While most wild tomato germplasms contain the BBX18CC allele and show more drought tolerant, modern cultivated tomatoes mostly carry BBX18TT allele and are more drought sensitive. Knockout of BBX18 leads to improved drought tolerance in transgenic plants of cultivated tomato. Ascorbate peroxidase 1 (APX1) is identified as a BBX18-interacting protein that acts as a positive regulator of drought resistance in tomato. Chromatin immunoprecipitation sequencing analyses reveal that BBX18 binds to a unique cis-acting element of the APX1 promoter and represses its gene expression. This study provides insights into the molecular mechanism underlying drought resistance mediated by the BBX18-APX1 module in plants. Elucidation of drought tolerance mechanism would help to develop tolerant crops. Here, the authors reveal the interaction between BBX18, a zinc finger transcription factor, and ascorbate peroxidase (APX1) to mediate drought stress resistance in tomato. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWAS) identified over fifty loci associated with lung cancer risk. However, underlying mechanisms and target genes are largely unknown, as most risk-associated variants might regulate gene expression in a context-specific manner. Here, we generate a barcode-shared transcriptome and chromatin accessibility map of 117,911 human lung cells from age/sex-matched ever- and never-smokers to profile context-specific gene regulation. Identified candidate cis-regulatory elements (cCREs) are largely cell type-specific, with 37% detected in one cell type. Colocalization of lung cancer candidate causal variants (CCVs) with these cCREs combined with transcription factor footprinting prioritize the variants for 68% of the GWAS loci. CCV-colocalization and trait relevance score indicate that epithelial and immune cell categories, including rare cell types, contribute to lung cancer susceptibility the most. A multi-level cCRE-gene linking system identifies candidate susceptibility genes from 57% of the loci, where most loci display cell-category-specific target genes, suggesting context-specific susceptibility gene function. Multiple genetic loci are associated with lung cancer risk, but the underlying genetic mechanisms remain poorly understood. Here, the authors perform single-cell RNA-seq and ATAC-seq analyses of lung cells from ever- and never-smokers; they report candidate cis-regulatory elements that colocalise with candidate causal variants in lung cancer risk loci and potential susceptibility genes. [ABSTRACT FROM AUTHOR]

Paracrine IL-2 signalling drives the CD8 + T cell expansion and differentiation that allow protection against viral infections, but the underlying molecular events are incompletely understood. Here we show that the transcription factor SRF, a master regulator of cytoskeletal gene expression, is required for effective IL-2 signalling during L. monocytogenes infection. Acting cell-autonomously with its actin-regulated cofactors MRTF-A and MRTF-B, SRF is dispensible for initial TCR-mediated CD8+ T cell proliferation, but is required for sustained IL-2 dependent CD8+ effector T cell expansion, and persistence of memory cells. Following TCR activation, Mrtfab-null CD8+ T cells produce IL-2 normally, but homotypic clustering is impaired both in vitro and in vivo. Expression of cytoskeletal structural and regulatory genes, most notably actins, is defective in Mrtfab-null CD8+ T cells. Activation-induced cell clustering in vitro requires F-actin assembly, and Mrtfab-null cell clusters are small, contain less F-actin, and defective in IL-2 retention. Clustering of Mrtfab-null cells can be partially restored by exogenous actin expression. IL-2 mediated CD8+ T cell proliferation during infection thus depends on the control of cytoskeletal dynamics and actin gene expression by MRTF-SRF signalling. The transcription factor SRF, together with its co-factors MRTF-A and MRTF-B, controls cytoskeletal gene expression. Here authors show that MRTF/SRF inactivation leads to decreased IL-2 mediated CD8 + T cell proliferation during infection, resulting from disrupted homotypic cell clustering and reduced retention of IL-2. [ABSTRACT FROM AUTHOR]

Drug-induced gene expression profiles can identify potential mechanisms of toxicity. We focus on obtaining signatures for cardiotoxicity of FDA-approved tyrosine kinase inhibitors (TKIs) in human induced-pluripotent-stem-cell-derived cardiomyocytes, using bulk transcriptomic profiles. We use singular value decomposition to identify drug-selective patterns across cell lines obtained from multiple healthy human subjects. Cellular pathways affected by cardiotoxic TKIs include energy metabolism, contractile, and extracellular matrix dynamics. Projecting these pathways to published single cell expression profiles indicates that TKI responses can be evoked in both cardiomyocytes and fibroblasts. Integration of transcriptomic outlier analysis with whole genomic sequencing of our six cell lines enables us to correctly reidentify a genomic variant causally linked to anthracycline-induced cardiotoxicity and predict genomic variants potentially associated with TKI-induced cardiotoxicity. We conclude that mRNA expression profiles when integrated with publicly available genomic, pathway, and single cell transcriptomic datasets, provide multiscale signatures for cardiotoxicity that could be used for drug development and patient stratification. Using a new computational pipeline for identification of drug-selective transcriptomic responses and FAERS data, the authors identified potential pathways and genomic variants indicative of cancer drug cardiotoxicity in iPSC-derived cardiomyocytes. [ABSTRACT FROM AUTHOR]

Nitrogen (N) deficiency responses are essential for plant survival and reproduction. Here, via an expression genome-wide association study (eGWAS), we reveal a mechanism that regulates microRNA (miRNA) dynamics necessary for N deficiency responses in Arabidopsis. Differential expression levels of three NAC transcription factor (TF) genes involved in leaf N deficiency responses among Arabidopsis accessions are most significantly associated with polymorphisms in HASTY (HST), which encodes an importin/exportin family protein responsible for the generation of mature miRNAs. HST acts as a negative regulator of N deficiency-induced leaf senescence, and the disruption and overexpression of HST differently modifies miRNA dynamics in response to N deficiency, altering levels of miRNAs targeting transcripts. Interestingly, N deficiency prevents the interaction of HST with HST-interacting proteins, DCL1 and RAN1, and some miRNAs. This suggests that HST-mediated regulation of miRNA dynamics collectively controls regulations mediated by multiple N deficiency response-associated NAC TFs, thereby being central to the N deficiency response network. Here Sakuraba et al. show that HASTY-mediated regulation of miRNA dynamics controls regulations mediated by multiple N deficiency response-associated NAC transcription factors, thereby being central to the N deficiency response network in Arabidopsis. [ABSTRACT FROM AUTHOR]

Burkholderia cenocepacia is an opportunistic and infective bacterium containing an orphan DNA methyltransferase called M.BceJIV with roles in regulating gene expression and motility of the bacterium. M.BceJIV recognizes a GTWWAC motif (where W can be an adenine or a thymine) and methylates N6 of the adenine at the fifth base position. Here, we present crystal structures of M.BceJIV/DNA/sinefungin ternary complex and allied biochemical, computational, and thermodynamic analyses. Remarkably, the structures show not one, but two DNA substrates bound to the M.BceJIV dimer, with each monomer contributing to the recognition of two recognition sequences. We also show that methylation at the two recognition sequences occurs independently, and that the GTWWAC motifs are enriched in intergenic regions in the genomes of B. cenocepacia strains. We further computationally assess the interactions underlying the affinities of different ligands (SAM, SAH, and sinefungin) for M.BceJIV, as a step towards developing selective inhibitors for limiting B. cenocepacia infection. Crystal structures of DNA methyltransferase M.BceJIV in complex with DNA and sinefungin reveal an unusual mode of DNA binding and methylation, wherein each M.BceJIV monomer contributes to the recognition of two DNA sequences. [ABSTRACT FROM AUTHOR]

Glioblastoma remains one of the deadliest brain malignancies. First-line therapy consists of maximal surgical tumor resection, accompanied by chemotherapy and radiotherapy. Malignant cells escape surgical resection by migrating into the surrounding healthy brain tissue, where they give rise to the recurrent tumor. Based on gene expression, tumor cores can be subtyped into mesenchymal, proneural, and classical tumors, each being associated with differences in genetic alterations and cellular composition. In contrast, the adjacent brain parenchyma where infiltrating malignant cells escape surgical resection is less characterized in patients. Using spatial transcriptomics (n = 11), we show that malignant cells within proneural or mesenchymal tumor cores display spatially organized differences in gene expression, although such differences decrease within the infiltrated brain tissue. Malignant cells residing in infiltrated brain tissue have increased expression of genes related to neurodevelopmental pathways and glial cell differentiation. Our findings provide an updated view of the spatial landscape of glioblastomas and further our understanding of the malignant cells that infiltrate the healthy brain, providing new avenues for the targeted therapy of these cells after surgical resection. In glioblastoma, malignant cells can escape surgical resection by migrating into the surrounding healthy brain tissue. Here, the authors use spatial transcriptomics to characterise the tumour core and infiltrated brain regions in glioblastoma, and show shifts in malignant cell composition and molecular pathways with potential clinical implications. [ABSTRACT FROM AUTHOR]

Progressive supranuclear palsy (PSP), a rare Parkinsonian disorder, is characterized by problems with movement, balance, and cognition. PSP differs from Alzheimer's disease (AD) and other diseases, displaying abnormal microtubule-associated protein tau by both neuronal and glial cell pathologies. Genetic contributors may mediate these differences; however, the genetics of PSP remain underexplored. Here we conduct the largest genome-wide association study (GWAS) of PSP which includes 2779 cases (2595 neuropathologically-confirmed) and 5584 controls and identify six independent PSP susceptibility loci with genome-wide significant (P < 5 × 10−8) associations, including five known (MAPT, MOBP, STX6, RUNX2, SLCO1A2) and one novel locus (C4A). Integration with cell type-specific epigenomic annotations reveal an oligodendrocytic signature that might distinguish PSP from AD and Parkinson's disease in subsequent studies. Candidate PSP risk gene prioritization using expression quantitative trait loci (eQTLs) identifies oligodendrocyte-specific effects on gene expression in half of the genome-wide significant loci, and an association with C4A expression in brain tissue, which may be driven by increased C4A copy number. Finally, histological studies demonstrate tau aggregates in oligodendrocytes that colocalize with C4 (complement) deposition. Integrating GWAS with functional studies, epigenomic and eQTL analyses, we identify potential causal roles for variation in MOBP, STX6, RUNX2, SLCO1A2, and C4A in PSP pathogenesis. The authors present the largest genome-wide association study to date for a rare Parkinsonian disorder, progressive supranuclear palsy (PSP). They include follow-up investigations of the identified susceptibility loci, functional consequences, and cell-specific pathologies, providing insights into genetic and molecular mechanisms underlying PSP. [ABSTRACT FROM AUTHOR]

Naïve pluripotency is sustained by a self-reinforcing gene regulatory network (GRN) comprising core and naïve pluripotency-specific transcription factors (TFs). Upon exiting naïve pluripotency, embryonic stem cells (ESCs) transition through a formative post-implantation-like pluripotent state, where they acquire competence for lineage choice. However, the mechanisms underlying disengagement from the naïve GRN and initiation of the formative GRN are unclear. Here, we demonstrate that phosphorylated AKT acts as a gatekeeper that prevents nuclear localisation of FoxO TFs in naïve ESCs. PTEN-mediated reduction of AKT activity upon exit from naïve pluripotency allows nuclear entry of FoxO TFs, enforcing a cell fate transition by binding and activating formative pluripotency-specific enhancers. Indeed, FoxO TFs are necessary and sufficient for the activation of the formative pluripotency-specific GRN. Our work uncovers a pivotal role for FoxO TFs in establishing formative post-implantation pluripotency, a critical early embryonic cell fate transition. While the core regulators involved in maintenance of naïve and formative pluripotency have been described, less is known about the regulation of the transition between the two states. Here they show that FoxO transcription factors mediate the naïve to formative pluripotency transition and are regulated via AKT phosphorylation. [ABSTRACT FROM AUTHOR]

Genomic imprinting is observed in endosperm, a placenta-like seed tissue, where transposable elements (TEs) and repeat-derived small RNAs (sRNAs) mediate epigenetic changes in plants. In imprinting, uniparental gene expression arises due to parent-specific epigenetic marks on one allele but not on the other. The importance of sRNAs and their regulation in endosperm development or in imprinting is poorly understood in crops. Here we show that a previously uncharacterized CLASSY (CLSY)-family chromatin remodeler named OsCLSY3 is essential for rice endosperm development and imprinting, acting as an upstream player in the sRNA pathway. Comparative transcriptome and genetic analysis indicated its endosperm-preferred expression and its likely paternal imprinted nature. These important features are modulated by RNA-directed DNA methylation (RdDM) of tandemly arranged TEs in its promoter. Upon perturbation of OsCLSY3 in transgenic lines, we observe defects in endosperm development and a loss of around 70% of all sRNAs. Interestingly, well-conserved endosperm-specific sRNAs (siren) that are vital for reproductive fitness in angiosperms are also dependent on OsCLSY3. We observed that many imprinted genes and seed development-associated genes are under the control of OsCLSY3. These results support an essential role of OsCLSY3 in rice endosperm development and imprinting, and propose similar regulatory strategies involving CLSY3 homologs among other cereals. Using a well-resolved spatial transcriptomics coupled with genetic data, this study identified rice CLSY3 chromatin remodeler as an upstream player in endosperm development. OsCLSY3 is imprinted and endosperm-specific, and, curiously, its expression is controlled by RdDM. [ABSTRACT FROM AUTHOR]

Near-infrared (NIR) probes are highly sought after as fluorescent tags for multicolor cellular and in vivo imaging. Here we develop small NIR fluorescent nanobodies, termed NIR-FbLAG16 and NIR-FbLAG30, enabling background-free visualization of various GFP-derived probes and biosensors. We also design a red-shifted variant, NIR-Fb(718), to simultaneously target several antigens within the NIR spectral range. Leveraging the antigen-stabilizing property of the developed NIR-Fbs, we then create two modular systems for precise control of gene expression in GFP-labeled cells. Applying the NIR-Fbs in vivo, we target cells expressing GFP and the calcium biosensor GCaMP6 in the somatosensory cortex of transgenic mice. Simultaneously tracking calcium activity and the reference signal from NIR-FbLAGs bound to GCaMP6 enables ratiometric deep-brain in vivo imaging. Altogether, NIR-FbLAGs present a promising approach for imaging and manipulating various processes in live cells and behaving animals expressing GFP-based probes. Near-infrared (NIR) probes are in high demand as fluorescent tags for multicolour cellular and in vivo imaging. Here authors develop series of NIR nanobodies that bind to GFP-based biosensors, and demonstrate ratiometric deep-brain in vivo imaging of active GCaMP6f-expressing neurons. [ABSTRACT FROM AUTHOR]

Activated Wnt/β-catenin pathway is a key genetic event in liver cancer development. Solute carrier (SLC) transporters are promising drug targets. Here, we identify SLC13A3 as a drug-targetable effector downstream of β-catenin in liver cancer. SLC13A3 expression is elevated in human liver cancer samples with gain of function (GOF) mutant CTNNB1, the gene encoding β-catenin. Activation of β-catenin up-regulates SLC13A3, leading to intracellular accumulation of endogenous SLC13A3 substrates. SLC13A3 is identified as a low-affinity transporter for glutathione (GSH). Silencing of SLC13A3 downregulates the leucine transporter SLC7A5 via c-MYC signaling, leading to leucine depletion and mTOR inactivation. Furthermore, silencing of SLC13A3 depletes GSH and induces autophagic ferroptosis in β-catenin-activated liver cancer cells. Importantly, both genetic inhibition of SLC13A3 and a small molecule SLC13A3 inhibitor suppress β-catenin-driven hepatocarcinogenesis in mice. Altogether, our study suggests that SLC13A3 could be a promising therapeutic target for treating human liver cancers with GOF CTNNB1 mutations. The therapeutic approaches directly targeting activated β-catenin in liver cancers are restricted by toxicities. Here the authors identify that the solute carrier transporter SLC13A3 is upregulated by activation of β-catenin and silencing of SLC13A3 induces ferroptosis, which could be exploited as a therapeutic opportunity in β-catenin-driven liver cancers. [ABSTRACT FROM AUTHOR]

Histone modifications are associated with distinct transcriptional states, but it is unclear whether they instruct gene expression. To investigate this, we mutate histone H3.3 K9 and K27 residues in mouse embryonic stem cells (mESCs). Here, we find that H3.3K9 is essential for controlling specific distal intergenic regions and for proper H3K27me3 deposition at promoters. The H3.3K9A mutation resulted in decreased H3K9me3 at regions encompassing endogenous retroviruses and induced a gain of H3K27ac and nascent transcription. These changes in the chromatin environment unleash cryptic enhancers, resulting in the activation of distinctive transcriptional programs and culminating in protein expression normally restricted to specialized immune cell types. The H3.3K27A mutant disrupts the deposition and spreading of the repressive H3K27me3 mark, particularly impacting bivalent genes with higher basal levels of H3.3 at promoters. Therefore, H3.3K9 and K27 crucially orchestrate repressive chromatin states at cis-regulatory elements and bivalent promoters, respectively, and instruct proper transcription in mESCs. In this study, the authors mutate histone H3.3 K9 and K27 residues, demonstrating their importance in maintaining repressive chromatin states at endogenous retrovirus-derived cryptic enhancers and bivalent promoters in mouse embryonic stem cells. [ABSTRACT FROM AUTHOR]

In in situ capturing-based spatial transcriptomics, spots of the same size and printed at fixed locations cannot precisely capture the randomly-located single cells, therefore inherently failing to profile transcriptome at the single-cell level. To this end, we present STIE, an Expectation Maximization algorithm that aligns the spatial transcriptome to its matched histology image-based nuclear morphology and recovers missing cells from ~70% gap area, thereby achieving the real single-cell level and whole-slide scale deconvolution, convolution, and clustering for both low- and high-resolution spots. STIE characterizes cell-type-specific gene expression and demonstrates outperforming concordance with true cell-type-specific transcriptomic signatures than the other spot- and subspot-level methods. Furthermore, STIE reveals the single-cell level insights, for instance, lower actual spot resolution than its reported spot size, unbiased evaluation of cell type colocalization, superior power of high-resolution spot in distinguishing nuanced cell types, and spatial cell-cell interactions at the single-cell level other than spot level. In situ capturing-based spatial transcriptomics cannot precisely capture randomly located single cells, regardless of its spot resolution. Here, authors integrate spot-level gene expression with histology images, computationally achieving single-cell level spatial transcriptomics. [ABSTRACT FROM AUTHOR]

Tandem repeats (TRs) are genomic regions that tandemly change in repeat number, which are often multiallelic. Their characteristics and contributions to gene expression and quantitative traits in rice are largely unknown. Here, we survey rice TR variations based on 231 genome assemblies and the rice pan-genome graph. We identify 227,391 multiallelic TR loci, including 54,416 TR variations that are absent from the Nipponbare reference genome. Only 1/3 TR variations show strong linkage with nearby bi-allelic variants (SNPs, Indels and PAVs). Using 193 panicle and 202 leaf transcriptomic data, we reveal 485 and 511 TRs act as QTLs independently of other bi-allelic variations to nearby gene expression, respectively. Using plant height and grain width as examples, we identify and validate TRs contributions to rice agronomic trait variations. These findings would enhance our understanding of the functions of multiallelic variants and facilitate rice molecular breeding. Tandem repeats (TRs) have unique ability to drive a range of phenotype variations. Here, the authors survey rice TR variations based on 231 genome assemblies and the rice pan-genome graph, identify TR variations associated with expressed genes, and reveal expression TRs contributed to rice agronomic trait variations. [ABSTRACT FROM AUTHOR]

The integration of extrinsic signaling with cell-intrinsic transcription factors can direct progenitor cells to differentiate into distinct cell fates. In the developing Drosophila eye, differentiation of photoreceptors R1–R7 requires EGFR signaling mediated by the transcription factor Pointed, and our single-cell RNA-Seq analysis shows that the same photoreceptors require the eye-specific transcription factor Glass. We find that ectopic expression of Glass and activation of EGFR signaling synergistically induce neuronal gene expression in the wing disc in a Pointed-dependent manner. Targeted DamID reveals that Glass and Pointed share many binding sites in the genome of developing photoreceptors. Comparison with transcriptomic data shows that Pointed and Glass induce photoreceptor differentiation through intermediate transcription factors, including the redundant homologs Scratch and Scrape, as well as directly activating neuronal effector genes. Our data reveal synergistic activation of a multi-layered transcriptional network as the mechanism by which EGFR signaling induces neuronal identity in Glass-expressing cells. The authors show that synergistic activation by two transcription factors - Pointed, which is activated by cell–cell signaling, and Glass, which is eye-specific - drives a program of neuronal gene expression in developing Drosophila photoreceptors. [ABSTRACT FROM AUTHOR]

Upon infecting its vertebrate host, the malaria parasite initially invades the liver where it undergoes massive replication, whilst remaining clinically silent. The coordination of host responses across the complex liver tissue during malaria infection remains unexplored. Here, we perform spatial transcriptomics in combination with single-nuclei RNA sequencing over multiple time points to delineate host-pathogen interactions across Plasmodium berghei-infected liver tissues. Our data reveals significant changes in spatial gene expression in the malaria-infected tissues. These include changes related to lipid metabolism in the proximity to sites of Plasmodium infection, distinct inflammation programs between lobular zones, and regions with enrichment of different inflammatory cells, which we term 'inflammatory hotspots'. We also observe significant upregulation of genes involved in inflammation in the control liver tissues of mice injected with mosquito salivary gland components. However, this response is considerably delayed compared to that observed in P. berghei-infected mice. Our study establishes a benchmark for investigating transcriptome changes during host-parasite interactions in tissues, it provides informative insights regarding in vivo study design linked to infection and offers a useful tool for the discovery and validation of de novo intervention strategies aimed at malaria liver stage infection. During malaria transmission, the liver acts as a portal into the vertebrate host and is a major vaccine target. Here, Hildebrandt et al combine spatial and single cell transcriptomics to delineate host-parasite interactions within distinct spatial regions of the tissue. [ABSTRACT FROM AUTHOR]

Enoblituzumab, an immunotherapeutic agent targeting CD276, shows both safety and efficacy in activating T cells and oligodendrocyte-like cells against various cancers. Preclinical studies and mouse models suggest that therapies targeting CD276 may outperform PD1/PD-L1 blockade. However, data from mouse models indicate a significant non-responsive population to anti-CD276 treatment, with the mechanisms of resistance still unclear. In this study, we evaluate the activity of anti-CD276 antibodies in a chemically-induced murine model of head and neck squamous cell carcinoma. Using models of induced and orthotopic carcinogenesis, we identify ITGB6 as a key gene mediating differential responses to anti-CD276 treatment. Through single-cell RNA sequencing and gene-knockout mouse models, we find that ITGB6 regulates the expression of the tumor-associated chemokine CX3CL1, which recruits and activates PF4+ macrophages that express high levels of CX3CR1. Inhibition of the CX3CL1-CX3CR1 axis suppresses the infiltration and secretion of CXCL16 by PF4+ macrophages, thereby reinvigorating cytotoxic CXCR6+ CD8+ T cells and enhancing sensitivity to anti-CD276 treatment. Further investigations demonstrate that inhibiting ITGB6 restores sensitivity to PD1 antibodies in mice resistant to anti-PD1 treatment. In summary, our research reveals a resistance mechanism associated with immune checkpoint inhibitor therapy and identifies potential targets to overcome resistance in cancer treatment. Response rate to anti-CD276 based immunotherapy remains suboptimal in patients with cancer. Here, in a chemically-induced murine model of head and neck squamous cell carcinoma, the authors show that expression of ITGB6 in tumor cells promotes resistance to anti-CD276 therapy, associated with accumulation of PF4+ macrophages and T cell dysfunction. [ABSTRACT FROM AUTHOR]

Mitochondrial oxidative phosphorylation (OXPHOS) fuels cellular ATP demands. OXPHOS defects lead to severe human disorders with unexplained tissue specific pathologies. Mitochondrial gene expression is essential for OXPHOS biogenesis since core subunits of the complexes are mitochondrial-encoded. COX14 is required for translation of COX1, the central mitochondrial-encoded subunit of complex IV. Here we describe a COX14 mutant mouse corresponding to a patient with complex IV deficiency. COX14M19I mice display broad tissue-specific pathologies. A hallmark phenotype is severe liver inflammation linked to release of mitochondrial RNA into the cytosol sensed by RIG-1 pathway. We find that mitochondrial RNA release is triggered by increased reactive oxygen species production in the deficiency of complex IV. Additionally, we describe a COA3Y72C mouse, affected in an assembly factor that cooperates with COX14 in early COX1 biogenesis, which displays a similar yet milder inflammatory phenotype. Our study provides insight into a link between defective mitochondrial gene expression and tissue-specific inflammation. Dysfunctions of the mitochondrial OXPHOS system lead to severe human disorders. Here, the authors analyzed a mouse mutant that mimics a mitochondrial disorder patient and find that reactive oxygen species trigger RNA release and inflammatory pathway in liver. [ABSTRACT FROM AUTHOR]

HMGA1 is an abundant non-histone chromatin protein that has been implicated in embryonic development, cancer, and cellular senescence, but its specific role remains elusive. Here, we combine functional genomics approaches with graph theory to investigate how HMGA1 genomic deposition controls high-order chromatin networks in an oncogene-induced senescence model. While the direct role of HMGA1 in gene activation has been described previously, we find little evidence to support this. Instead, we show that the heterogeneous linear distribution of HMGA1 drives a specific 3D chromatin organization. HMGA1-dense loci form highly interactive networks, similar to, but independent of, constitutive heterochromatic loci. This, coupled with the exclusion of HMGA1-poor chromatin regions, leads to coordinated gene regulation through the repositioning of genes. In the absence of HMGA1, the whole process is largely reversed, but many regulatory interactions also emerge, amplifying the inflammatory senescence-associated secretory phenotype. Such HMGA1-mediated fine-tuning of gene expression contributes to the heterogeneous nature of senescence at the single-cell level. A similar 'buffer' effect of HMGA1 on inflammatory signalling is also detected in lung cancer cells. Our study reveals a mechanism through which HMGA1 modulates chromatin compartmentalization and gene regulation in senescence and beyond. HMGA1 helps regulate the topology of the chromatin network, supporting compartmentalization. Here, Olan et al. show that, in oncogene-induced senescence, genes are included or excluded from HMGA1 cores, suggesting HMGA1 has a fine-tuning role. [ABSTRACT FROM AUTHOR]

Lipid nanoparticle-assisted mRNA inhalation therapy necessitates addressing challenges such as resistance to shear force damage, mucus penetration, cellular internalization, rapid lysosomal escape, and target protein expression. Here, we introduce the innovative "LOOP" platform with a four-step workflow to develop inhaled lipid nanoparticles specifically for pulmonary mRNA delivery. iLNP-HP08LOOP featuring a high helper lipid ratio, acidic dialysis buffer, and excipient-assisted nebulization buffer, demonstrates exceptional stability and enhanced mRNA expression in the lungs. By incorporating mRNA encoding IL-11 single chain fragment variable (scFv), scFv@iLNP-HP08LOOP effectively delivers and secretes IL-11 scFv to the lungs of male mice, significantly inhibiting fibrosis. This formulation surpasses both inhaled and intravenously injected IL-11 scFv in inhibiting fibroblast activation and extracellular matrix deposition. The HP08LOOP system is also compatible with commercially available ALC0315 LNPs. Thus, the "LOOP" method presents a powerful platform for developing inhaled mRNA nanotherapeutics with potential for treating various respiratory diseases, including idiopathic pulmonary fibrosis. Inhaled lipid nanoparticles face several challenges when delivered to the lungs. Here, the authors show "LOOP"-optimized iLNP-HP08LOOP could withstand shear damage generated during nebulization process and be used to develop mRNA-mediated antibody therapy for treating idiopathic pulmonary fibrosis. [ABSTRACT FROM AUTHOR]

Cell fate specification occurs along invariant species-specific trajectories that define the animal body plan. This process is controlled by gene regulatory networks that regulate the expression of the limited set of transcription factors encoded in animal genomes. Here we globally assess the spatial expression of ~90% of expressed transcription factors during sea urchin development from embryo to larva to determine the activity of gene regulatory networks and their regulatory states during cell fate specification. We show that >200 embryonically expressed transcription factors together define >70 cell fates that recapitulate the morphological and functional organization of this organism. Most cell fate-specific regulatory states consist of ~15–40 transcription factors with similarity particularly among functionally related cell types regardless of developmental origin. Temporally, regulatory states change continuously during development, indicating that progressive changes in regulatory circuit activity determine cell fate specification. We conclude that the combinatorial expression of transcription factors provides molecular definitions that suffice for the unique specification of cell states in time and space during embryogenesis. Gene regulatory networks (GRNs) control developmental processes that define the body plan of multicellular animals. In this study, Valencia and Peter monitored the activity of GRNs during embryogenesis in sea urchins by determining the spatial and temporal changes in the expression of over 200 transcription factors. [ABSTRACT FROM AUTHOR]

Numerous studies have now demonstrated that lncRNAs can influence gene expression programs leading to cell and organismal phenotypes. Typically, lncRNA perturbations and concomitant changes in gene expression are measured on the timescale of many hours to days. Thus, we currently lack a temporally grounded understanding of the primary, secondary, and tertiary relationships of lncRNA-mediated transcriptional and epigenetic regulation—a prerequisite to elucidating lncRNA mechanisms. To begin to address when and where a lncRNA regulates gene expression, we genetically engineer cell lines to temporally induce the lncRNA Firre. Using this approach, we are able to monitor lncRNA transcriptional regulatory events from 15 min to four days. We observe that upon induction, Firre RNA regulates epigenetic and transcriptional states in trans within 30 min. These early regulatory events result in much larger transcriptional changes after 12 h, well before current studies monitor lncRNA regulation. Moreover, Firre-mediated gene expression changes are epigenetically remembered for days. Overall, this study suggests that lncRNAs can rapidly regulate gene expression by establishing persistent epigenetic and transcriptional states. Typically, lncRNA perturbations are measured on the timescale of many hours to days. Here, the authors investigate when and where the lncRNA Firre influences epigenetic and transcriptional states, suggesting that lncRNA-mediated gene regulation occurs within minutes and is retained for days. [ABSTRACT FROM AUTHOR]

Structural variation heavily influences the molecular landscape of cancer, in part by impacting DNA methylation-mediated transcriptional regulation. Here, using multi-omic datasets involving >2400 pediatric brain and central nervous system tumors of diverse histologies from the Children's Brain Tumor Network, we report hundreds of genes and associated CpG islands (CGIs) for which the nearby presence of somatic structural variant (SV) breakpoints is recurrently associated with altered expression or DNA methylation, respectively, including tumor suppressor genes ATRX and CDKN2A. Altered DNA methylation near enhancers associates with nearby somatic SV breakpoints, including MYC and MYCN. A subset of genes with SV-CGI methylation associations also have expression associations with patient survival, including BCOR, TERT, RCOR2, and PDLIM4. DNA methylation changes in recurrent or progressive tumors compared to the initial tumor within the same patient can predict survival in pediatric and adult cancers. Our comprehensive and pan-histology genomic analyses reveal mechanisms of noncoding alterations impacting cancer genes. Somatic structural variants (SVs) in cancer can impact DNA methylation-mediated transcriptional regulation. Here, the authors analyse multi-omics data from over 2400 samples from the Children's Brain Tumor Network and report SVs that are associated with altered gene expression or DNA methylation, including some with prognostic relevance. [ABSTRACT FROM AUTHOR]

Immunochemotherapy has been the mainstay of treatment for newly diagnosed diffuse large B-cell lymphoma (ndDLBCL) yet is inadequate for many patients. In this work, we perform unsupervised clustering on transcriptomic features from a large cohort of ndDLBCL patients and identify seven clusters, one called A7 with poor prognosis, and develop a classifier to identify these clusters in independent ndDLBCL cohorts. This high-risk cluster is enriched for activated B-cell cell-of-origin, low immune infiltration, high MYC expression, and copy number aberrations. We compare and contrast our methodology with recent DLBCL classifiers to contextualize our clusters and show improved prognostic utility. Finally, using pre-clinical models, we demonstrate a mechanistic rationale for IKZF1/3 degraders such as lenalidomide to overcome the low immune infiltration phenotype of A7 by inducing T-cell trafficking into tumors and upregulating MHC I and II on tumor cells, and demonstrate that TCF4 is an important regulator of MYC-related biology in A7. Researchers identified 7 biological group of diffuse large B-cell lymphoma, one with poor prognosis. They developed a gene expression classifier to detect these groups, potentially improving prognosis and future treatments for the high-risk patients. [ABSTRACT FROM AUTHOR]

Multiple reference panels of a given tissue or multiple tissues often exist, and multiple regression methods could be used for training gene expression imputation models for transcriptome-wide association studies (TWAS). To leverage expression imputation models (i.e., base models) trained with multiple reference panels, regression methods, and tissues, we develop a Stacked Regression based TWAS (SR-TWAS) tool which can obtain optimal linear combinations of base models for a given validation transcriptomic dataset. Both simulation and real studies show that SR-TWAS improves power, due to increased training sample sizes and borrowed strength across multiple regression methods and tissues. Leveraging base models across multiple reference panels, tissues, and regression methods, our real studies identify 6 independent significant risk genes for Alzheimer's disease (AD) dementia for supplementary motor area tissue and 9 independent significant risk genes for Parkinson's disease (PD) for substantia nigra tissue. Relevant biological interpretations are found for these significant risk genes. Here the authors develop a Stacked Regression-based TWAS method for obtaining optimal expression imputation models with a given validation transcriptomic dataset, by combining base models trained with multiple reference panels, regression methods, and tissues. [ABSTRACT FROM AUTHOR]

Single-cell analysis across multiple samples and conditions requires quantitative modeling of the interplay between the continuum of cell states and the technical and biological sources of sample-to-sample variability. We introduce GEDI, a generative model that identifies latent space variations in multi-sample, multi-condition single-cell datasets and attributes them to sample-level covariates. GEDI enables cross-sample cell state mapping on par with state-of-the-art integration methods, cluster-free differential gene expression analysis along the continuum of cell states, and machine learning-based prediction of sample characteristics from single-cell data. GEDI can also incorporate gene-level prior knowledge to infer pathway and regulatory network activities in single cells. Finally, GEDI extends all these concepts to previously unexplored modalities that require joint consideration of dual measurements, such as the joint analysis of exon inclusion/exclusion reads to model alternative cassette exon splicing, or spliced/unspliced reads to model the mRNA stability landscapes of single cells. Single-cell analysis of multi-condition cohorts requires modelling the interaction between sample variables and cell states. Here, authors develop GEDI to enable integration, cluster-free differential expression analysis and regulon analysis for both gene expression and alternative splicing modalities. [ABSTRACT FROM AUTHOR]

The ependyma lining the third ventricle (3V) in the mediobasal hypothalamus plays a crucial role in energy balance and glucose homeostasis. It is characterized by a high functional heterogeneity and plasticity, but the underlying molecular mechanisms governing its features are not fully understood. Here, 5481 hypothalamic ependymocytes were cataloged using FACS-assisted scRNAseq from fed, 12h-fasted, and 24h-fasted adult male mice. With standard clustering analysis, typical ependymal cells and β2-tanycytes appear sharply defined, but other subpopulations, β1- and α-tanycytes, display fuzzy boundaries with few or no specific markers. Pseudospatial approaches, based on the 3V neuroanatomical distribution, enable the identification of specific versus shared tanycyte markers and subgroup-specific versus general tanycyte functions. We show that fasting dynamically shifts gene expression patterns along the 3V, leading to a spatial redistribution of cell type-specific responses. Altogether, we show that changes in energy status induce metabolic and functional switches in tanycyte subpopulations, providing insights into molecular and functional diversity and plasticity within the tanycyte population. The ependyma lining the third ventricle in the mediobasal hypothalamus regulates energy balance and glucose homeostasis. Here, using single-cell transcriptomics the authors reveal how hypothalamic ependymocyte subpopulations spatially reorganize functions of lipid processing and neuron communication during energy imbalance. [ABSTRACT FROM AUTHOR]

Epigenetic modifications, such as 5-methylcytosine (5mC), can sometimes be transmitted between generations, provoking speculation that epigenetic changes could play a role in adaptation and evolution. Here, we use experimental evolution to investigate how 5mC levels evolve in populations of biparental insect (Nicrophorus vespilloides) derived from a wild source population and maintained independently under different regimes of parental care in the lab. We show that 5mC levels in the transcribed regions of genes (gene bodies) diverge between populations that have been exposed to different levels of care for 30 generations. These changes in 5mC do not reflect changes in the levels of gene expression. However, the accumulation of 5mC within genes between populations is associated with reduced variability in gene expression within populations. Our results suggest that evolved change in 5mC could contribute to phenotypic evolution by influencing variability in gene expression in invertebrates. It's largely unknown whether heritable, but non-genetic, changes can affect evolution. This study shows that changes in epigenetic modifications that evolve in burying beetles are associated with altered variability in gene expression, potentially influencing adaptation. [ABSTRACT FROM AUTHOR]

Foam cells in atheroma are engorged with lipid droplets (LDs) that contain esters of regulatory lipids whose metabolism remains poorly understood. LD-associated hydrolase (LDAH) has a lipase structure and high affinity for LDs of foam cells. Using knockout and transgenic mice of both sexes, here we show that LDAH inhibits atherosclerosis development and promotes stable lesion architectures. Broad and targeted lipidomic analyzes of primary macrophages and comparative lipid profiling of atheroma identified a broad impact of LDAH on esterified sterols, including natural liver X receptor (LXR) sterol ligands. Transcriptomic analyzes coupled with rescue experiments show that LDAH modulates the expression of prototypical LXR targets and leads macrophages to a less inflammatory phenotype with a profibrotic gene signature. These studies underscore the role of LDs as reservoirs and metabolic hubs of bioactive lipids, and suggest that LDAH favorably modulates macrophage activation and protects against atherosclerosis via lipolytic mobilization of regulatory sterols. The mechanisms and consequences of bioactive lipids release from lipid droplets remain poorly understood. Here the authors link a lipid droplet enzyme to mobilization of esters of regulatory sterols in foam cells and protection against atherosclerosis. [ABSTRACT FROM AUTHOR]

Gene-gene interactions are crucial to the control of sub-cellular processes but our understanding of their stochastic dynamics is hindered by the lack of simulation methods that can accurately and efficiently predict how the distributions of gene product numbers vary across parameter space. To overcome these difficulties, here we present Holimap (high-order linear-mapping approximation), an approach that approximates the protein or mRNA number distributions of a complex gene regulatory network by the distributions of a much simpler reaction system. We demonstrate Holimap's computational advantages over conventional methods by applying it to predict the stochastic time-dependent dynamics of various gene networks, including transcriptional networks ranging from simple autoregulatory loops to complex randomly connected networks, post-transcriptional networks, and post-translational networks. Holimap is ideally suited to study how the intricate network of gene-gene interactions results in precise coordination and control of gene expression. Understanding gene-gene interactions is key to decoding cellular processes. Here, authors introduce Holimap, a method that accurately and rapidly predicts gene expression dynamics by approximating complex regulatory networks with simpler models. [ABSTRACT FROM AUTHOR]

Hyperuricemia is associated with an increased risk of gout, hypertension, diabetes, and cardiovascular diseases. Most mammals maintain normal serum uric acid (SUA) via urate oxidase (Uox), an enzyme that metabolizes poorly-soluble UA to highly-soluble allantoin. In contrast, Uox became a pseudogene in humans and apes over the long course of evolution. Here we demonstrate an atavistic strategy for treating hyperuricemia based on endogenous expression of Uox in hepatocytes mediated by mRNA (mUox) loaded with an ionizable lipid nanoparticle termed iLAND. mUox@iLAND allows effective transfection and protein expression in vitro. A single dose of mUox@iLAND lowers SUA levels for several weeks in two female murine models, including a novel long-lasting model, which is also confirmed by metabolomics analysis. Together with the excellent safety profiles observed in vivo, the proposed mRNA agent demonstrates substantial potential for hyperuricemia therapy and the prevention of associated conditions. Urate oxidase (Uox) lost its function during evolution. Here the authors propose an atavistic strategy to treat hyperuricemia by using a proprietary lipid nanoparticle to load Uox mRNA, which reduced the serum uric acid levels effectively in two animal models. [ABSTRACT FROM AUTHOR]

Genetic variation in the human leukocyte antigen (HLA) loci is associated with risk of immune-mediated diseases, but the molecular effects of HLA polymorphism are unclear. Here we examined the effects of HLA genetic variation on the expression of 2940 plasma proteins across 45,330 Europeans in the UK Biobank, with replication analyses across multiple ancestry groups. We detected 504 proteins affected by HLA variants (HLA-pQTL), including widespread trans effects by autoimmune disease risk alleles. More than 80% of the HLA-pQTL fine-mapped to amino acid positions in the peptide binding groove. HLA-I and II affected proteins expressed in similar cell types but in different pathways of both adaptive and innate immunity. Finally, we investigated potential HLA-pQTL effects on disease by integrating HLA-pQTL with fine-mapped HLA-disease signals in the UK Biobank. Our data reveal the diverse effects of HLA genetic variation and aid the interpretation of associations between HLA alleles and immune-mediated diseases. Genetic variation in the HLA locus is associated with many traits, including autoimmune diseases. Here, the authors show that HLA genetic variation exerts widespread trans effects on plasma protein expression, aiding interpretation of associations between HLA alleles and immune mediated diseases. [ABSTRACT FROM AUTHOR]

Multiple Sclerosis (MS) is a heterogeneous inflammatory and neurodegenerative disease with an unpredictable course towards progressive disability. Treating progressive MS is challenging due to limited insights into the underlying mechanisms. We examined the molecular changes associated with primary progressive MS (PPMS) using a cross-tissue (blood and post-mortem brain) and multilayered data (genetic, epigenetic, transcriptomic) from independent cohorts. In PPMS, we found hypermethylation of the 1q21.1 locus, controlled by PPMS-specific genetic variations and influencing the expression of proximal genes (CHD1L, PRKAB2) in the brain. Evidence from reporter assay and CRISPR/dCas9 experiments supports a causal link between methylation and expression and correlation network analysis further implicates these genes in PPMS brain processes. Knock-down of CHD1L in human iPSC-derived neurons and knock-out of chd1l in zebrafish led to developmental and functional deficits of neurons. Thus, several lines of evidence suggest a distinct genetic-epigenetic-transcriptional interplay in the 1q21.1 locus potentially contributing to PPMS pathogenesis. This study identifies distinct genetic-epigenetictranscriptional interplay at the 1q21.1 locus in relation to neuronal deficits in primary progressive multiple sclerosis (PPMS). [ABSTRACT FROM AUTHOR]

In plants, the conserved plant-specific photoreceptor UV RESISTANCE LOCUS 8 (UVR8) perceives ultraviolet-B (UV-B) light and mediates UV-B-induced photomorphogenesis and stress acclimation. In this study, we reveal that UV-B light treatment shortens seedlings, increases stem thickness, and enhances UV-B stress tolerance in rice (Oryza sativa) via its two UV-B photoreceptors OsUVR8a and OsUVR8b. Although the rice and Arabidopsis (Arabidopsis thaliana) UVR8 (AtUVR8) photoreceptors all form monomers in response to UV-B light, OsUVR8a, and OsUVR8b function is only partially conserved with respect to AtUVR8 in UV-B-induced photomorphogenesis and stress acclimation. UV-B light and CONSTITUTIVELY PHOTOMORPHOGENIC 1 (COP1) promote the nuclear accumulation of AtUVR8; by contrast, OsUVR8a and OsUVR8b constitutively localize to the nucleus via their own nuclear localization signals, independently of UV-B light and the RING-finger mutation of OsCOP1. We show that OsCOP1 negatively regulates UV-B responses, and shows weak interaction with OsUVR8s, which is ascribed to the N terminus of OsCOP1, which is conserved in several monocots. Furthermore, transcriptome analysis demonstrates that UV-B-responsive gene expression differs globally between Arabidopsis and rice, illuminating the evolutionary divergence of UV-B light signaling pathways between monocot and dicot plants. Differently from Arabidopsis, rice UV-B photoreceptors show a non-classical localization mechanism and protein interaction mode to mediate UV-B responses, suggestive of evolutionary divergence of UV-B light signaling between monocot and dicot plants. [ABSTRACT FROM AUTHOR]

The increasing utilization of mouse models in human neuroscience research places higher demands on computational methods to translate findings from the mouse brain to the human one. In this study, we develop BrainAlign, a self-supervised learning approach, for the whole brain alignment of spatial transcriptomics (ST) between humans and mice. BrainAlign encodes spots and genes simultaneously in two separated shared embedding spaces by a heterogeneous graph neural network. We demonstrate that BrainAlign could integrate cross-species spots into the embedding space and reveal the conserved brain regions supported by ST information, which facilitates the detection of homologous regions between humans and mice. Genomic analysis further presents gene expression connections between humans and mice and reveals similar expression patterns for marker genes. Moreover, BrainAlign can accurately map spatially similar homologous regions or clusters onto a unified spatial structural domain while preserving their relative positions. Comparative transcriptomics of whole brains across species is vital in neuroscience. Here, authors develop a deep learning method, BrainAlign, to align spatial transcriptomics across human and mouse brains. BrainAlign identifies conserved brain regions and uncovers similar patterns for marker genes. [ABSTRACT FROM AUTHOR]

Gene expression is regulated by changes in chromatin architecture intrinsic to cellular differentiation and as an active response to environmental stimuli. Chromatin dynamics are a major driver of phenotypic diversity, regulation of development, and manifestation of disease. Remarkably, we know little about the evolutionary dynamics of chromatin reorganisation through time, data essential to characterise the impact of environmental stress during the ongoing biodiversity extinction crisis (20th–21st century). Linking the disparate fields of chromatin biology and museum science through their common use of the preservative formaldehyde (a constituent of formalin), we have generated historical chromatin profiles in museum specimens up to 117 years old. Historical chromatin profiles are reproducible, tissue-specific, sex-specific, and environmental condition-dependent in vertebrate specimens. Additionally, we show that over-fixation modulates differential chromatin accessibility to enable semi-quantitative estimates of relative gene expression in vertebrates and a yeast model. Our approach transforms formalin-fixed biological collections into an accurate, comprehensive, and global record of environmental impact on gene expression and phenotype. Formaldehyde-preserved museum specimens have produced genetic data. Here, the authors generate chromatin profiles from museum specimens 117 years old and experimentally demonstrate chromatin profile presence in formalin-fixed mouse and yeast models. [ABSTRACT FROM AUTHOR]

Engineering metabolism to efficiently produce chemicals from multi-step pathways requires optimizing multi-gene expression programs to achieve enzyme balance. CRISPR-Cas transcriptional control systems are emerging as important tools for programming multi-gene expression, but poor predictability of guide RNA folding can disrupt expression control. Here, we correlate efficacy of modified guide RNAs (scRNAs) for CRISPR activation (CRISPRa) in E. coli with a computational kinetic parameter describing scRNA folding rate into the active structure (rS = 0.8). This parameter also enables forward design of scRNAs, allowing us to design a system of three synthetic CRISPRa promoters that can orthogonally activate (>35-fold) expression of chosen outputs. Through combinatorial activation tuning, we profile a three-dimensional design space expressing two different biosynthetic pathways, demonstrating variable production of pteridine and human milk oligosaccharide products. This RNA design approach aids combinatorial optimization of metabolic pathways and may accelerate routine design of effective multi-gene regulation programs in bacterial hosts. Guide RNA folding affects functionality of CRISPR-Cas transcriptional control systems. Here, the authors report computational gRNA design together with creation of synthetic CRISPRa promoters for orthogonal expression control and demonstrate the application in pteridine and human milk oligosaccharide production in E. coli. [ABSTRACT FROM AUTHOR]

Engineered DNA will slow the growth of a host cell if it redirects limiting resources or otherwise interferes with homeostasis. Escape mutants that alleviate this burden can rapidly evolve and take over cell populations, making genetic engineering less reliable and predictable. Synthetic biologists often use genetic parts encoded on plasmids, but their burden is rarely characterized. We measured how 301 BioBrick plasmids affected Escherichia coli growth and found that 59 (19.6%) were burdensome, primarily because they depleted the limited gene expression resources of host cells. Overall, no BioBricks reduced the growth rate of E. coli by >45%, which agreed with a population genetic model that predicts such plasmids should be unclonable. We made this model available online for education (https://barricklab.org/burden-model) and added our burden measurements to the iGEM Registry. Our results establish a fundamental limit on what DNA constructs and genetic modifications can be successfully engineered into cells. Engineered DNA will slow the growth of a host cell if it redirects limiting resources or otherwise interferes with homeostasis. Here the authors measure how 301 BioBrick plasmids affected Escherichia coli growth and found that 19.6% were burdensome, primarily because they depleted the limited gene expression resources of host cells. [ABSTRACT FROM AUTHOR]

Mutations in the Cockayne Syndrome group B (CSB) gene cause cancer in mice, but premature aging and severe neurodevelopmental defects in humans. CSB, a member of the SWI/SNF family of chromatin remodelers, plays diverse roles in regulating gene expression and transcription-coupled nucleotide excision repair (TC-NER); however, these functions do not explain the distinct phenotypic differences observed between CSB-deficient mice and humans. During investigating Cockayne Syndrome-associated genome instability, we uncover an intrinsic mechanism that involves elongating RNA polymerase II (RNAPII) undergoing transient pauses at internal T-runs where CSB is required to propel RNAPII forward. Consequently, CSB deficiency retards RNAPII elongation in these regions, and when coupled with G-rich sequences upstream, exacerbates genome instability by promoting R-loop formation. These R-loop prone motifs are notably abundant in relatively long genes related to neuronal functions in the human genome, but less prevalent in the mouse genome. These findings provide mechanistic insights into differential impacts of CSB deficiency on mice versus humans and suggest that the manifestation of the Cockayne Syndrome phenotype in humans results from the progressive evolution of mammalian genomes. CSB deficiency is associated with induction of R-loops and genome instability, impacting long genes linked to neuronal functions in humans. Here the authors provide mechanistic understanding on how mutations in the CSB gene affects RNAPII elongation and genome instability. [ABSTRACT FROM AUTHOR]

Duchenne muscular dystrophy (DMD) affecting 1 in 3500–5000 live male newborns is the frequently fatal genetic disease resulted from various mutations in DMD gene encoding dystrophin protein. About 70% of DMD-causing mutations are exon deletion leading to frameshift of open reading frame and dystrophin deficiency. To facilitate translating human DMD-targeting CRISPR therapeutics into patients, we herein establish a genetically humanized mouse model of DMD by replacing exon 50 and 51 of mouse Dmd gene with human exon 50 sequence. This humanized mouse model recapitulats patient's DMD phenotypes of dystrophin deficiency and muscle dysfunction. Furthermore, we target splicing sites in human exon 50 with adenine base editor to induce exon skipping and robustly restored dystrophin expression in heart, tibialis anterior and diaphragm muscles. Importantly, systemic delivery of base editor via adeno-associated virus in the humanized male mouse model improves the muscle function of DMD mice to the similar level of wildtype ones, indicating the therapeutic efficacy of base editing strategy in treating most of DMD types with exon deletion or point mutations via exon-skipping induction. Duchenne muscular dystrophy (DMD) results from mutations in the DMD. Here the authors report using adenine base editing of exon splice sites to restore dystrophin expression and muscle function in a humanized mouse model of DMD with exon 50 and 51 of mouse DMD gene replaced with human exon 50. [ABSTRACT FROM AUTHOR]