A promoter-level mammalian expression atlas

Core promoterome of barley embryo

Precise localization and dissection of gene promoters are key to understanding transcriptional gene regulation and to successful bioengineering applications. The core RNA polymerase II initiation machinery is highly conserved among eukaryotes, leading to a general expectation of equivalent underlying mechanisms. Still, less is known about promoters in the plant kingdom. In this study, we employed cap analysis of gene expression (CAGE) at three embryonic developmental stages in barley to accurately map, annotate, and quantify transcription initiation events. Unsupervised discovery of de novo sequence clusters grouped promoters based on characteristic initiator and position-specific core-promoter motifs. This grouping was complemented by the annotation of transcription factor binding site (TFBS) motifs. Integration with genome-wide epigenomic data sets and gene ontology (GO) enrichment analysis further delineated the chromatin environments and functional roles of genes associated with distinct promoter categories. The TATA-box presence governs all features explored, supporting the general model of two separate genomic regulatory environments. We describe the extent and implications of alternative transcription initiation events, including those that are specific to developmental stages, which can affect the protein sequence or the presence of regions that regulate translation. The generated promoterome dataset provides a valuable genomic resource for enhancing the functional annotation of the barley genome. It also offers insights into the transcriptional regulation of individual genes and presents opportunities for the informed manipulation of promoter architecture, with the aim of enhancing traits of agronomic importance.

A Functional Survey of the Regulatory Landscape of Estrogen Receptor-Positive Breast Cancer Evolution

Only a handful of somatic alterations have been linked to endocrine therapy resistance in hormone-dependent breast cancer, potentially explaining ∼40% of relapses. If other mechanisms underlie the evolution of hormone-dependent breast cancer under adjuvant therapy is currently unknown. In this work, we employ functional genomics to dissect the contribution of cis-regulatory elements (CRE) to cancer evolution by focusing on 12 megabases of noncoding DNA, including clonal enhancers, gene promoters, and boundaries of topologically associating domains. Parallel epigenetic perturbation (CRISPRi) in vitro reveals context-dependent roles for many of these CREs, with a specific impact on dormancy entrance and endocrine therapy resistance. Profiling of CRE somatic alterations in a unique, longitudinal cohort of patients treated with endocrine therapies identifies a limited set of noncoding changes potentially involved in therapy resistance. Overall, our data uncover how endocrine therapies trigger the emergence of transient features which could ultimately be exploited to hinder the adaptive process. Significance: This study shows that cells adapting to endocrine therapies undergo changes in the usage or regulatory regions. Dormant cells are less vulnerable to regulatory perturbation but gain transient dependencies which can be exploited to decrease the formation of dormant persisters.

Extracting regulatory active chromatin footprint from cell-free DNA

Cell-free DNA (cfDNA) has emerged as a pivotal player in precision medicine, revolutionizing the diagnostic and therapeutic landscape. While its clinical applications have significantly increased in recent years, current cfDNA assays have limited ability to identify the active transcriptional programs that govern complex disease phenotypes and capture the heterogeneity of the disease. To address these limitations, we have developed a non-invasive platform to enrich and examine the active chromatin fragments (cfDNAac) in peripheral blood. The deconvolution of the cfDNAac signal from traditional nucleosomal chromatin fragments (cfDNAnuc) yields a catalog of features linking these circulating chromatin signals in blood to specific regulatory elements across the genome, including enhancers, promoters, and highly transcribed genes, mirroring the epigenetic data from the ENCODE project. Notably, these cfDNAac counts correlate strongly with RNA polymerase II activity and exhibit distinct expression patterns for known circadian genes. Additionally, cfDNAac signals across gene bodies and promoters show strong correlations with whole blood gene expression levels defined by GTEx. This study illustrates the utility of cfDNAac analysis for investigating epigenomics and gene expression, underscoring its potential for a wide range of clinical applications in precision medicine.

Angiocrine signaling in sinusoidal homeostasis and liver diseases

The hepatic sinusoids are composed of liver sinusoidal endothelial cells (LSECs), which are surrounded by hepatic stellate cells (HSCs) and contain liver-resident macrophages called Kupffer cells, and other patrolling immune cells. All these cells communicate with each other and with hepatocytes to maintain sinusoidal homeostasis and a spectrum of hepatic functions under healthy conditions. Sinusoidal homeostasis is disrupted by metabolites, toxins, viruses, and other pathological factors, leading to liver injury, chronic liver diseases, and cirrhosis. Alterations in hepatic sinusoids are linked to fibrosis progression and portal hypertension. LSECs are crucial regulators of cellular crosstalk within their microenvironment via angiocrine signaling. This review discusses the mechanisms by which angiocrine signaling orchestrates sinusoidal homeostasis, as well as the development of liver diseases. Here, we summarise the crosstalk between LSECs, HSCs, hepatocytes, cholangiocytes, and immune cells in health and disease and comment on potential novel therapeutic methods for treating liver diseases.

Discovery of NRG1-VII: the myeloid-derived class of NRG1

The growth factor Neuregulin-1 (NRG1) has pleiotropic roles in proliferation and differentiation of the stem cell niche in different tissues. It has been implicated in gut, brain and muscle development and repair. Six isoform classes of NRG1 and over 28 protein isoforms have been previously described. Here we report a new class of NRG1, designated NRG1-VII to denote that these NRG1 isoforms arise from a myeloid-specific transcriptional start site (TSS) previously uncharacterized. Long-read sequencing was used to identify eight high-confidence NRG1-VII transcripts. These transcripts presented major structural differences from one another, through the use of cassette exons and alternative stop codons. Expression of NRG1-VII was confirmed in primary human monocytes and tissue resident macrophages and induced pluripotent stem cell-derived macrophages (iPSC-derived macrophages). Isoform switching via cassette exon usage and alternate polyadenylation was apparent during monocyte maturation and macrophage differentiation. NRG1-VII is the major class expressed by the myeloid lineage, including tissue-resident macrophages. Analysis of public gene expression data indicates that monocytes and macrophages are a primary source of NRG1. The size and structure of class VII isoforms suggests that they may be more diffusible through tissues than other NRG1 classes. However, the specific roles of class VII variants in tissue homeostasis and repair have not yet been determined.

STAG2 mutations reshape the cohesin-structured spatial chromatin architecture to drive gene regulation in acute myeloid leukemia

Cohesin shapes the chromatin architecture, including enhancer-promoter interactions. Its components, especially STAG2, but not its paralog STAG1, are frequently mutated in myeloid malignancies. To elucidate the underlying mechanisms of leukemogenesis, we comprehensively characterized genetic, transcriptional, and chromatin conformational changes in acute myeloid leukemia (AML) patient samples. Specific loci displayed altered cohesin occupancy, gene expression, and local chromatin activation, which were not compensated by the remaining STAG1-cohesin. These changes could be linked to disrupted spatial chromatin looping in cohesin-mutated AMLs. Complementary depletion of STAG2 or STAG1 in primary human hematopoietic progenitors (HSPCs) revealed effects resembling STAG2-mutant AML-specific changes following STAG2 knockdown, not invoked by the depletion of STAG1. STAG2-deficient HSPCs displayed impaired differentiation capacity and maintained HSPC-like gene expression. This work establishes STAG2 as a key regulator of chromatin contacts, gene expression, and differentiation in the hematopoietic system and identifies candidate target genes that may be implicated in human leukemogenesis.

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

BLEND: Probabilistic Cellular Deconvolution with Automated Reference Selection

Cellular deconvolution aims to estimate cell type fractions from bulk transcriptomic and other omics data. Most existing deconvolution methods fail to account for the heterogeneity in cell type-specific (CTS) expression across bulk samples, ignore discrepancies between CTS expression in bulk and cell type reference data, and provide no guidance on cell type reference selection or integration. To address these issues, we introduce BLEND, a hierarchical Bayesian method that leverages multiple reference datasets. BLEND learns the most suitable references for each bulk sample by exploring the convex hulls of references and employs a "bag-of-words" representation for bulk count data for deconvolution. To speed up the computation, we provide an efficient EM algorithm for parameter estimation. Notably, BLEND requires no data transformation, normalization, cell type marker gene selection, or reference quality evaluation. Benchmarking studies on both simulated and real human brain data highlight BLEND's superior performance in various scenarios. The analysis of Alzheimer's disease data illustrates BLEND's application in real data and reference resource integration.

Whole-exome and whole-genome sequencing of 1064 individuals with type 1 diabetes reveals novel genes for diabetic kidney disease

AIMS/HYPOTHESIS

Diabetic kidney disease (DKD) is a severe diabetic complication that affects one third of individuals with type 1 diabetes. Although several genes and common variants have been shown to be associated with DKD, much of the predicted inheritance remains unexplained. Here, we performed next-generation sequencing to assess whether low-frequency variants, extending to a minor allele frequency (MAF) ≤10% (single or aggregated) contribute to the missing heritability in DKD.

METHODS

We performed whole-exome sequencing (WES) of 498 individuals and whole-genome sequencing (WGS) of 599 individuals with type 1 diabetes. After quality control, next-generation sequencing data were available for a total of 1064 individuals, of whom 541 had developed either severe albuminuria or end-stage kidney disease, and 523 had retained normal albumin excretion despite a long duration of type 1 diabetes. Single-variant and gene-aggregate tests for protein-altering variants (PAV) and protein-truncating variants (PTV) were performed separately for WES and WGS data and combined in a meta-analysis. We also performed genome-wide aggregate analyses on genomic windows (sliding window), promoters and enhancers using the WGS dataset.

RESULTS

In the single-variant meta-analysis, no variant reached genome-wide significance, but a suggestively associated common THAP7 rs369250 variant (p=1.50 × 10-5, MAF=49%) was replicated in the FinnGen general population genome-wide association study (GWAS) data for chronic kidney disease and DKD phenotypes. The gene-aggregate meta-analysis provided suggestive evidence (p<4.0 × 10-4) at four genes for DKD, of which NAT16 (MAFPAV≤10%) and LTA (also known as TNFβ, MAFPAV≤5%) are replicated in the FinnGen general population GWAS data. The LTA rs2229092 C allele was associated with significantly lower TNFR1, TNFR2 and TNFR3 serum levels in a subset of FinnDiane participants. Of the intergenic regions suggestively associated with DKD, the enhancer on chromosome 18q12.3 (p=3.94 × 10-5, MAFvariants≤5%) showed interaction with the METTL4 gene; the lead variant was replicated, and predicted to alter binding of the MafB transcription factor.

CONCLUSIONS/INTERPRETATION

Our sequencing-based meta-analysis revealed multiple genes, variants and regulatory regions that were suggestively associated with DKD. However, as no variant or gene reached genome-wide significance, further studies are needed to validate the findings.

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

A spatial expression atlas of the adult human proximal small intestine

The mouse small intestine shows profound variability in gene expression along the crypt-villus axis1,2. Whether similar spatial heterogeneity exists in the adult human gut remains unclear. Here we use spatial transcriptomics, spatial proteomics and single-molecule fluorescence in situ hybridization to reconstruct a comprehensive spatial expression atlas of the adult human proximal small intestine. We describe zonated expression and cell type representation for epithelial, mesenchymal and immune cell types. We find that migrating enterocytes switch from lipid droplet assembly and iron uptake at the villus bottom to chylomicron biosynthesis and iron release at the tip. Villus tip cells are pro-immunogenic, recruiting γδ T cells and macrophages to the tip, in contrast to their immunosuppressive roles in mouse. We also show that the human small intestine contains abundant serrated and branched villi that are enriched at the tops of circular folds. Our study presents a detailed resource for understanding the biology of the adult human small intestine.

Potential Transcriptional Enhancers in Coronaviruses: From Infectious Bronchitis Virus to SARS-CoV-2

Coronaviruses constitute a global threat to human and animal health. It is essential to investigate the long-distance RNA-RNA interactions that approximate remote regulatory elements in strategies, including genome circularization, discontinuous transcription, and transcriptional enhancers, aimed at the rapid replication of their large genomes, pathogenicity, and immune evasion. Based on the primary sequences and modeled RNA-RNA interactions of two experimentally defined coronaviral enhancers, we detected via an in silico primary and secondary structural analysis potential enhancers in various coronaviruses, from the phylogenetically ancient avian infectious bronchitis virus (IBV) to the recently emerged SARS-CoV-2. These potential enhancers possess a core duplex-forming region that could transition between closed and open states, as molecular switches directed by viral or host factors. The duplex open state would pair with remote sequences in the viral genome and modulate the expression of downstream crucial genes involved in viral replication and host immune evasion. Consistently, variations in the predicted IBV enhancer region or its distant targets coincide with cases of viral attenuation, possibly driven by decreased open reading frame (ORF)3a immune evasion protein expression. If validated experimentally, the annotated enhancer sequences could inform structural prediction tools and antiviral interventions.

Genetic and selective constraints on the optimization of gene product diversity

RNA and protein expressed from the same gene can have diverse isoforms due to various post-transcriptional and post-translational modifications. For the vast majority of alternative isoforms, It is unknown whether they are adaptive or simply biological noise. As we cannot experimentally probe the function of each isoform, we can ask whether the distribution of isoforms across genes and across species is consistent with expectations from different evolutionary processes. However, there is currently no theoretical framework that can generate such predictions. To address this, we developed a mathematical model where isoform abundances are determined collectively by cis-acting loci, trans-acting factors, gene expression levels, and isoform decay rates to predict isoform abundance distributions across species and genes in the face of mutation, genetic drift, and selection. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. Notably, suboptimal phenotypes are more likely to evolve when the population is small and/or when the number of cis-loci is large. We also explored scenarios where modification processes have both beneficial and detrimental effects, revealing a non-monotonic relationship between effective population size and optimization, demonstrating how opposing selection pressures on cis- and trans-acting loci can constrain the optimization of gene product diversity. As a demonstration of the power of our theory, we compared the expected distribution of A-to-I RNA editing levels in coleoids and found this to be largely consistent with non-adaptive explanations.

Pan-cancer transcriptome analysis reveals widespread regulation through alternative tandem transcription initiation

Abnormal transcription initiation from alternative first exon has been reported to promote tumorigenesis. However, the prevalence and impact of gene expression regulation mediated by alternative tandem transcription initiation were mostly unknown in cancer. Here, we developed a robust computational method to analyze alternative tandem transcription start site (TSS) usage from standard RNA sequencing data. Applying this method to pan-cancer RNA sequencing datasets, we observed widespread dysregulation of tandem TSS usage in tumors, many of which were independent of changes in overall expression level or alternative first exon usage. We showed that the dynamics of tandem TSS usage was associated with epigenomic modulation. We found that significant 5' untranslated region shortening of gene TIMM13 contributed to increased protein production, and up-regulation of TIMM13 by CRISPR-mediated transcriptional activation promoted proliferation and migration of lung cancer cells. Our findings suggest that dysregulated tandem TSS usage represents an addtional layer of cancer-associated transcriptome alterations.

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-located variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-located variants to enhance model building. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these models to large GWAS datasets for breast, prostate, and lung cancers. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene prediction models and identifying disease-associated genes, as evidenced by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

A critical review of the impact of candidate copy number variants on autism spectrum disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder (NDD) influenced by genetic, epigenetic, and environmental factors. Recent advancements in genomic analysis have shed light on numerous genes associated with ASD, highlighting the significant role of both common and rare genetic mutations, as well as copy number variations (CNVs), single nucleotide polymorphisms (SNPs) and unique de novo variants. These genetic variations disrupt neurodevelopmental pathways, contributing to the disorder's complexity. Notably, CNVs are present in 10 %-20 % of individuals with autism, with 3 %-7 % detectable through cytogenetic methods. While the role of submicroscopic CNVs in ASD has been recently studied, their association with genomic loci and genes has not been thoroughly explored. In this review, we focus on 47 CNV regions linked to ASD, encompassing 1632 genes, including protein-coding genes and long non-coding RNAs (lncRNAs), of which 659 show significant brain expression. Using a list of ASD-associated genes from SFARI, we detect 17 regions harboring at least one known ASD-related protein-coding gene. Of the remaining 30 regions, we identify 24 regions containing at least one protein-coding gene with brain-enriched expression and a nervous system phenotype in mouse mutants, and one lncRNA with both brain-enriched expression and upregulation in iPSC to neuron differentiation. This review not only expands our understanding of the genetic diversity associated with ASD but also underscores the potential of lncRNAs in contributing to its etiology. Additionally, the discovered CNVs will be a valuable resource for future diagnostic, therapeutic, and research endeavors aimed at prioritizing genetic variations in ASD.

Proteome-scale tagging and functional screening in mammalian cells by ORFtag

The systematic determination of protein function is a key goal of modern biology, but remains challenging with current approaches. Here we present ORFtag, a versatile, cost-effective and highly efficient method for the massively parallel tagging and functional interrogation of proteins at the proteome scale. ORFtag uses retroviral vectors bearing a promoter, peptide tag and splice donor to generate fusions between the tag and endogenous open reading frames (ORFs). We demonstrate the utility of ORFtag through functional screens for transcriptional activators, repressors and posttranscriptional regulators in mouse embryonic stem cells. Each screen recovers known and identifies new regulators, including long ORFs inaccessible by other methods. Among other hits, we find that Zfp574 is a highly selective transcriptional activator and that oncogenic fusions often function as transactivators.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

Alternative TSS use is widespread in Cryptococcus fungi in response to environmental cues and regulated genome-wide by the transcription factor Tur1

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain.

The annotation of GBA1 has been concealed by its protein-coding pseudogene GBAP1

Mutations in GBA1 cause Gaucher disease and are the most important genetic risk factor for Parkinson's disease. However, analysis of transcription at this locus is complicated by its highly homologous pseudogene, GBAP1. We show that >50% of short RNA-sequencing reads mapping to GBA1 also map to GBAP1. Thus, we used long-read RNA sequencing in the human brain, which allowed us to accurately quantify expression from both GBA1 and GBAP1. We discovered significant differences in expression compared to short-read data and identify currently unannotated transcripts of both GBA1 and GBAP1. These included protein-coding transcripts from both genes that were translated in human brain, but without the known lysosomal function-yet accounting for almost a third of transcription. Analyzing brain-specific cell types using long-read and single-nucleus RNA sequencing revealed region-specific variations in transcript expression. Overall, these findings suggest nonlysosomal roles for GBA1 and GBAP1 with implications for our understanding of the role of GBA1 in health and disease.

CapTrap-seq: a platform-agnostic and quantitative approach for high-fidelity full-length RNA sequencing

Long-read RNA sequencing is essential to produce accurate and exhaustive annotation of eukaryotic genomes. Despite advancements in throughput and accuracy, achieving reliable end-to-end identification of RNA transcripts remains a challenge for long-read sequencing methods. To address this limitation, we develop CapTrap-seq, a cDNA library preparation method, which combines the Cap-trapping strategy with oligo(dT) priming to detect 5' capped, full-length transcripts. In our study, we evaluate the performance of CapTrap-seq alongside other widely used RNA-seq library preparation protocols in human and mouse tissues, employing both ONT and PacBio sequencing technologies. To explore the quantitative capabilities of CapTrap-seq and its accuracy in reconstructing full-length RNA molecules, we implement a capping strategy for synthetic RNA spike-in sequences that mimics the natural 5'cap formation. Our benchmarks, incorporating the Long-read RNA-seq Genome Annotation Assessment Project (LRGASP) data, demonstrate that CapTrap-seq is a competitive, platform-agnostic RNA library preparation method for generating full-length transcript sequences.

Protocol for extracting and isolating porcine bone-marrow-derived macrophages from ribs

Due to anatomical and biological similarities with humans, pigs are increasingly used for inflammation- and immune-related studies in biomedical research, including the field of osteonecrosis and osteoimmunology. Here, we present a protocol for rib extraction, isolation of the bone marrow by centrifugation, and processing to obtain bone-marrow-derived macrophages (BMDMs). Then, we describe the procedures of in vitro experiments to evaluate the cell phenotype. For complete details on the use and execution of this protocol, please refer to Andre et al.1.

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Whole-genome sequencing identifies variants in ANK1, LRRN1, HAS1, and other genes and regulatory regions for stroke in type 1 diabetes

Individuals with type 1 diabetes (T1D) carry a markedly increased risk of stroke, with distinct clinical and neuroimaging characteristics as compared to those without diabetes. Using whole-exome or whole-genome sequencing of 1,051 individuals with T1D, we aimed to find rare and low-frequency genomic variants associated with stroke in T1D. We analysed the genome comprehensively with single-variant analyses, gene aggregate analyses, and aggregate analyses on genomic windows, enhancers and promoters. In addition, we attempted replication in T1D using a genome-wide association study (N = 3,945) and direct genotyping (N = 3,263), and in the general population from the large-scale population-wide FinnGen project and UK Biobank summary statistics. We identified a rare missense variant on SREBF1 exome-wide significantly associated with stroke (rs114001633, p.Pro227Leu, p-value = 7.30 × 10-8), which replicated for hemorrhagic stroke in T1D. Using gene aggregate analysis, we identified exome-wide significant genes: ANK1 and LRRN1 displayed replication evidence in T1D, and LRRN1, HAS1 and UACA in the general population (UK Biobank). Furthermore, we performed sliding-window analyses and identified 14 genome-wide significant windows for stroke on 4q33-34.1, of which two replicated in T1D, and a suggestive genomic window on LINC01500, which replicated in T1D. Finally, we identified a suggestively stroke-associated TRPM2-AS promoter (p-value = 5.78 × 10-6) with borderline significant replication in T1D, which we validated with an in vitro cell-based assay. Due to the rarity of the identified genetic variants, future replication of the genomic regions represented here is required with sequencing of individuals with T1D. Nevertheless, we here report the first genome-wide analysis on stroke in individuals with diabetes.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

RNAirport: a deep neural network-based database characterizing representative gene models in plants

A 5'-leader, known initially as the 5'-untranslated region, contains multiple isoforms due to alternative splicing (aS) and alternative transcription start site (aTSS). Therefore, a representative 5'-leader is demanded to examine the embedded RNA regulatory elements in controlling translation efficiency. Here, we develop a ranking algorithm and a deep-learning model to annotate representative 5'-leaders for five plant species. We rank the intra-sample and inter-sample frequency of aS-mediated transcript isoforms using the Kruskal-Wallis test-based algorithm and identify the representative aS-5'-leader. To further assign a representative 5'-end, we train the deep-learning model 5'leaderP to learn aTSS-mediated 5'-end distribution patterns from cap-analysis gene expression data. The model accurately predicts the 5'-end, confirmed experimentally in Arabidopsis and rice. The representative 5'-leader-contained gene models and 5'leaderP can be accessed at RNAirport (http://www.rnairport.com/leader5P/). The Stage 1 annotation of 5'-leader records 5'-leader diversity and will pave the way to Ribo-Seq open-reading frame annotation, identical to the project recently initiated by human GENCODE.

Transcriptional Hubs Within Cliques in Ensemble Hi-C Chromatin Interaction Networks

Chromatin conformation capture technologies permit the study of chromatin spatial organization on a genome-wide scale at a variety of resolutions. Despite the increasing precision and resolution of high-throughput chromatin conformation capture (Hi-C) methods, it remains challenging to conclusively link transcriptional activity to spatial organizational phenomena. We have developed a clique-based approach for analyzing Hi-C data that helps identify chromosomal hotspots that feature considerable enrichment of chromatin annotations for transcriptional start sites and, building on previously published work, show that these chromosomal hotspots are not only significantly enriched in RNA polymerase II binding sites as identified by the ENCODE project, but also identify a noticeable increase in FANTOM5 and GTEx transcription within our identified cliques across a variety of tissue types. From the obtained data, we surmise that our cliques are a suitable method for identifying transcription factories in Hi-C data, and outline further extensions to the method that may make it useful for locating regions of increased transcriptional activity in datasets where in-depth expression or polymerase data may not be available.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

Alternative isoform expression of key thermogenic genes in human beige adipocytes

Background

The beneficial effect of thermogenic adipocytes in maintaining body weight and protecting against metabolic disorders has raised interest in understanding the regulatory mechanisms defining white and beige adipocyte identity. Although alternative splicing has been shown to propagate adipose browning signals in mice, this has yet to be thoroughly investigated in human adipocytes.

Methods

We performed parallel white and beige adipogenic differentiation using primary adipose stem cells from 6 unrelated healthy subjects and assessed differential gene and isoform expression in mature adipocytes by RNA sequencing.

Results

We find 777 exon junctions with robust differential usage between white and beige adipocytes in all 6 subjects, mapping to 562 genes. Importantly, only 10% of these differentially spliced genes are also differentially expressed, indicating that alternative splicing constitutes an additional layer of gene expression regulation during beige adipocyte differentiation. Functional classification of alternative isoforms points to a gain of function for key thermogenic transcription factors such as PPARG and CITED1, and enzymes such as PEMT, or LPIN1. We find that a large majority of the splice variants arise from differential TSS usage, with beige-specific TSSs being enriched for PPARγ and MED1 binding compared to white-specific TSSs. Finally, we validate beige specific isoform expression at the protein level for two thermogenic regulators, PPARγ and PEMT.

Discussion

These results suggest that differential isoform expression through alternative TSS usage is an important regulatory mechanism for human adipocyte thermogenic specification.

Histone H3 posttranslational modified enzymes defined neutrophil plasticity and their vulnerability to IL-10 in the course of the inflammation

BACKGROUND

Neutrophils are a heterogeneous population capable of antimicrobial functions associated with pre-activation/activation and tissue regeneration. The specific polarisation of immune cells is mediated by the modification of 'chromatin landscapes', which enables differentiated access and activity of regulatory elements that guarantee their plasticity during inflammation No specific pattern within histone posttranslational modifications (PTMs) controlling this plasticity has been identified.

METHODS

Using the in vitro model of inflammation, reflecting different states of neutrophils from resting, pre-activated cells to activated and reducing tissue regeneration, we have analysed 11 different histone posttranslational modifications (PTMs), PTM enzymes associated with remodelling neutrophil chromatin, and H3K4me3 ChIP-Seq Gene Ontology analysis focusing on the processes related to histone PTMs. These findings were verified by extrapolation to adequate clinical status, using neutrophils derived from the patients with sepsis (systemic septic inflammation with LPS-stimulated neutrophils), neuromyelitis optical spectrum disorders (aseptic inflammation with pre-activated neutrophils) and periodontitis (local self-limiting septic inflammation with IL-10-positive neutrophils).

RESULTS

Physiological activation of neutrophils comprises a pre-activation characterised by histone H3K27ac and H3K4me1, which position enhancers; direct LPS exposure is induced explicitly by H3K4me3 which marked Transcription Start Site (TSS) regions and low-level of H3K9me3, H3K79me2 and H3K27me3 which, in turn, marked repressed genes. Contrary to antimicrobial action, IL-10 positively induced levels of H3S10p and negatively H3K9me3, which characterised processes related to the activation of genes within heterochromatin mediated by CHD1 and H3K9me3 specific demethylase JMJD2A. IL-10 protects changes within histone PTMs induced by TNF or LPS that affected H3K4me3-specific methyltransferase SETD1A and MLL1. Neutrophils previously exposed to inflammatory factors become unvulnerable to IL-10 because previous LPS stimulation interrupts TSS regions marked by H3K4me3 of CHD1 and JMJD2A genes. Therefore, LPS-activated neutrophils are disabled to induce CHD1/JMJD2A enzymes by IL-10, making this process irreversible. Because transcription of JMJD2A and CHD1 also depends on TSS positioning by H3K4me3, neutrophils before LPS stimulation become insensitive to IL-10.

CONCLUSION

Neutrophils, once pre-activated by TNF or directly stimulated by LPS, become insensitive to the anti-inflammatory effects of IL-10, and vice versa; IL-10 protects neutrophils against these proinflammatory stimuli. This phenomenon is responsible for disturbing the natural process of resolving inflammation and tissue regeneration.

Lineage-determining transcription factor-driven promoters regulate cell type-specific macrophage gene expression

Mammalian promoters consist of multifarious elements, which make them unique and support the selection of the proper transcript variants required under diverse conditions in distinct cell types. However, their direct DNA-transcription factor (TF) interactions are mostly unidentified. Murine bone marrow-derived macrophages (BMDMs) are a widely used model for studying gene expression regulation. Thus, this model serves as a rich source of various next-generation sequencing data sets, including a large number of TF cistromes. By processing and integrating the available cistromic, epigenomic and transcriptomic data from BMDMs, we characterized the macrophage-specific direct DNA-TF interactions, with a particular emphasis on those specific for promoters. Whilst active promoters are enriched for certain types of typically methylatable elements, more than half of them contain non-methylatable and prototypically promoter-distal elements. In addition, circa 14% of promoters-including that of Csf1r-are composed exclusively of 'distal' elements that provide cell type-specific gene regulation by specialized TFs. Similar to CG-rich promoters, these also contain methylatable CG sites that are demethylated in a significant portion and show high polymerase activity. We conclude that this unusual class of promoters regulates cell type-specific gene expression in macrophages, and such a mechanism might exist in other cell types too.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

Mouse Genome Informatics: an integrated knowledgebase system for the laboratory mouse

Mouse Genome Informatics (MGI) is a federation of expertly curated information resources designed to support experimental and computational investigations into genetic and genomic aspects of human biology and disease using the laboratory mouse as a model system. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are core MGI databases that share data and system architecture. MGI serves as the central community resource of integrated information about mouse genome features, variation, expression, gene function, phenotype, and human disease models acquired from peer-reviewed publications, author submissions, and major bioinformatics resources. To facilitate integration and standardization of data, biocuration scientists annotate using terms from controlled metadata vocabularies and biological ontologies (e.g. Mammalian Phenotype Ontology, Mouse Developmental Anatomy, Disease Ontology, Gene Ontology, etc.), and by applying international community standards for gene, allele, and mouse strain nomenclature. MGI serves basic scientists, translational researchers, and data scientists by providing access to FAIR-compliant data in both human-readable and compute-ready formats. The MGI resource is accessible at https://informatics.jax.org. Here, we present an overview of the core data types represented in MGI and highlight recent enhancements to the resource with a focus on new data and functionality for MGD and GXD.

Annotation of nuclear lncRNAs based on chromatin interactions

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Molecular mechanisms reconstruction from single-cell multi-omics data with HuMMuS

MOTIVATION

The molecular identity of a cell results from a complex interplay between heterogeneous molecular layers. Recent advances in single-cell sequencing technologies have opened the possibility to measure such molecular layers of regulation.

RESULTS

Here, we present HuMMuS, a new method for inferring regulatory mechanisms from single-cell multi-omics data. Differently from the state-of-the-art, HuMMuS captures cooperation between biological macromolecules and can easily include additional layers of molecular regulation. We benchmarked HuMMuS with respect to the state-of-the-art on both paired and unpaired multi-omics datasets. Our results proved the improvements provided by HuMMuS in terms of transcription factor (TF) targets, TF binding motifs and regulatory regions prediction. Finally, once applied to snmC-seq, scATAC-seq and scRNA-seq data from mouse brain cortex, HuMMuS enabled to accurately cluster scRNA profiles and to identify potential driver TFs.

AVAILABILITY AND IMPLEMENTATION

HuMMuS is available at https://github.com/cantinilab/HuMMuS.

Multiplex, single-cell CRISPRa screening for cell type specific regulatory elements

CRISPR-based gene activation (CRISPRa) is a promising therapeutic approach for gene therapy, upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis- regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to a library of 493 gRNAs targeting candidate cis- regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically upregulating intended target genes and no other neighboring genes within 1 Mb, including gRNAs yielding upregulation of six autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes in neurons. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans- acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate therapeutically relevant genes in a cell type-specific manner.

Differences in 5'untranslated regions highlight the importance of translational regulation of dosage sensitive genes

BACKGROUND

Untranslated regions (UTRs) are important mediators of post-transcriptional regulation. The length of UTRs and the composition of regulatory elements within them are known to vary substantially across genes, but little is known about the reasons for this variation in humans. Here, we set out to determine whether this variation, specifically in 5'UTRs, correlates with gene dosage sensitivity.

RESULTS

We investigate 5'UTR length, the number of alternative transcription start sites, the potential for alternative splicing, the number and type of upstream open reading frames (uORFs) and the propensity of 5'UTRs to form secondary structures. We explore how these elements vary by gene tolerance to loss-of-function (LoF; using the LOEUF metric), and in genes where changes in dosage are known to cause disease. We show that LOEUF correlates with 5'UTR length and complexity. Genes that are most intolerant to LoF have longer 5'UTRs, greater TSS diversity, and more upstream regulatory elements than their LoF tolerant counterparts. We show that these differences are evident in disease gene-sets, but not in recessive developmental disorder genes where LoF of a single allele is tolerated.

CONCLUSIONS

Our results confirm the importance of post-transcriptional regulation through 5'UTRs in tight regulation of mRNA and protein levels, particularly for genes where changes in dosage are deleterious and lead to disease. Finally, to support gene-based investigation we release a web-based browser tool, VuTR, that supports exploration of the composition of individual 5'UTRs and the impact of genetic variation within them.

Sequence basis of transcription initiation in the human genome

Transcription initiation is a process that is essential to ensuring the proper function of any gene, yet we still lack a unified understanding of sequence patterns and rules that explain most transcription start sites in the human genome. By predicting transcription initiation at base-pair resolution from sequences with a deep learning-inspired explainable model called Puffin, we show that a small set of simple rules can explain transcription initiation at most human promoters. We identify key sequence patterns that contribute to human promoter activity, each activating transcription with distinct position-specific effects. Furthermore, we explain the sequence basis of bidirectional transcription at promoters, identify the links between promoter sequence and gene expression variation across cell types, and explore the conservation of sequence determinants of transcription initiation across mammalian species.

Massively parallel reporter assays and mouse transgenic assays provide complementary information about neuronal enhancer activity

Genetic studies find hundreds of thousands of noncoding variants associated with psychiatric disorders. Massively parallel reporter assays (MPRAs) and in vivo transgenic mouse assays can be used to assay the impact of these variants. However, the relevance of MPRAs to in vivo function is unknown and transgenic assays suffer from low throughput. Here, we studied the utility of combining the two assays to study the impact of non-coding variants. We carried out an MPRA on over 50,000 sequences derived from enhancers validated in transgenic mouse assays and from multiple fetal neuronal ATAC-seq datasets. We also tested over 20,000 variants, including synthetic mutations in highly active neuronal enhancers and 177 common variants associated with psychiatric disorders. Variants with a high impact on MPRA activity were further tested in mice. We found a strong and specific correlation between MPRA and mouse neuronal enhancer activity including changes in neuronal enhancer activity in mouse embryos for variants with strong MPRA effects. Mouse assays also revealed pleiotropic variant effects that could not be observed in MPRA. Our work provides a large catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

The dynamic landscape of enhancer-derived RNA during mouse early embryo development

Enhancer-derived RNAs (eRNAs) play critical roles in diverse biological processes by facilitating their target gene expression. However, the abundance and function of eRNAs in early embryos are not clear. Here, we present a comprehensive eRNA atlas by systematically integrating publicly available datasets of mouse early embryos. We characterize the transcriptional and regulatory network of eRNAs and show that different embryo developmental stages have distinct eRNA expression and regulatory profiles. Paternal eRNAs are activated asymmetrically during zygotic genome activation (ZGA). Moreover, we identify an eRNA, MZGAe1, which plays an important function in regulating mouse ZGA and early embryo development. MZGAe1 knockdown leads to a developmental block from 2-cell embryo to blastocyst. We create an online data portal, M2ED2, to query and visualize eRNA expression and regulation. Our study thus provides a systematic landscape of eRNA and reveals the important role of eRNAs in regulating mouse early embryo development.

Transcriptome analysis of long non-coding RNAs in Mycobacterium avium complex-infected macrophages

Mycobacterium avium complex (MAC) is a non-tuberculous mycobacterium widely distributed in the environment. Even though MAC infection is increasing in older women and immunocompromised patients, to our knowledge there has been no comprehensive analysis of the MAC-infected host-cell transcriptome-and particularly of long non-coding RNAs (lncRNAs). By using in vitro-cultured primary mouse bone-marrow-derived macrophages (BMDMs) and Cap analysis of gene expression, we analyzed the transcriptional and kinetic landscape of macrophage genes, with a focus on lncRNAs, during MAC infection. MAC infection of macrophages induced the expression of immune/inflammatory response genes and other genes similar to those involved in M1 macrophage activation, consistent with previous reports, although Nos2 (M1 activation) and Arg1 (M2 activation) had distinct expression profiles. We identified 31 upregulated and 30 downregulated lncRNA promoters corresponding respectively to 18 and 26 lncRNAs. Upregulated lncRNAs were clustered into two groups-early and late upregulated-predicted to be associated with immune activation and the immune response to infection, respectively. Furthermore, an Ingenuity Pathway Analysis revealed canonical pathways and upstream transcription regulators associated with differentially expressed lncRNAs. Several differentially expressed lncRNAs reported elsewhere underwent expressional changes upon M1 or M2 preactivation and subsequent MAC infection. Finally, we showed that expressional change of lncRNAs in MAC-infected BMDMs was mediated by toll-like receptor 2, although there may be other mechanisms that sense MAC infection. We identified differentially expressed lncRNAs in MAC-infected BMDMs, revealing diverse features that imply the distinct roles of these lncRNAs in MAC infection and macrophage polarization.

Conditional splicing system for tight control of viral overlapping genes

Viral genomes frequently harbor overlapping genes, complicating the development of virus-vectored vaccines and gene therapies. This study introduces a novel conditional splicing system to precisely control the expression of such overlapping genes through recombinase-mediated conditional splicing. We refined site-specific recombinase (SSR) conditional splicing systems and explored their mechanisms. The systems demonstrated exceptional inducibility (116,700-fold increase) with negligible background expression, facilitating the conditional expression of overlapping genes in adenovirus-associated virus (AAV) and human immunodeficiency virus type 1. Notably, this approach enabled the establishment of stable AAV producer cell lines, encapsulating all necessary packaging genes. Our findings underscore the potential of the SSR-conditional splicing system to significantly advance vector engineering, enhancing the efficacy and scalability of viral-vector-based therapies and vaccines.

IMPORTANCE

Regulating overlapping genes is vital for gene therapy and vaccine development using viral vectors. The regulation of overlapping genes presents challenges, including cytotoxicity and impacts on vector capacity and genome stability, which restrict stable packaging cell line development and broad application. To address these challenges, we present a "loxp-splice-loxp"-based conditional splicing system, offering a novel solution for conditional expression of overlapping genes and stable cell line establishment. This system may also regulate other cytotoxic genes, representing a significant advancement in cell engineering and gene therapy as well as biomass production.

An integrated toolkit for human microglia functional genomics

BACKGROUND

Microglia, the brain's resident immune cells, play vital roles in brain development, and disorders like Alzheimer's disease (AD). Human iPSC-derived microglia (iMG) provide a promising model to study these processes. However, existing iMG generation protocols face challenges, such as prolonged differentiation time, lack of detailed characterization, and limited gene function investigation via CRISPR-Cas9.

METHODS

Our integrated toolkit for in-vitro microglia functional genomics optimizes iPSC differentiation into iMG through a streamlined two-step, 20-day process, producing iMG with a normal karyotype. We confirmed the iMG's authenticity and quality through single-cell RNA sequencing, chromatin accessibility profiles (ATAC-Seq), proteomics and functional tests. The toolkit also incorporates a drug-dependent CRISPR-ON/OFF system for temporally controlled gene expression. Further, we facilitate the use of multi-omic data by providing online searchable platform that compares new iMG profiles to human primary microglia: https://sherlab.shinyapps.io/IPSC-derived-Microglia/ .

RESULTS

Our method generates iMG that closely align with human primary microglia in terms of transcriptomic, proteomic, and chromatin accessibility profiles. Functionally, these iMG exhibit Ca2 + transients, cytokine driven migration, immune responses to inflammatory signals, and active phagocytosis of CNS related substrates including synaptosomes, amyloid beta and myelin. Significantly, the toolkit facilitates repeated iMG harvesting, essential for large-scale experiments like CRISPR-Cas9 screens. The standalone ATAC-Seq profiles of our iMG closely resemble primary microglia, positioning them as ideal tools to study AD-associated single nucleotide variants (SNV) especially in the genome regulatory regions.

CONCLUSIONS

Our advanced two-step protocol rapidly and efficiently produces authentic iMG. With features like the CRISPR-ON/OFF system and a comprehensive multi-omic data platform, our toolkit equips researchers for robust microglial functional genomic studies. By facilitating detailed SNV investigation and offering a sustainable cell harvest mechanism, the toolkit heralds significant progress in neurodegenerative disease drug research and therapeutic advancement.

Bayesian estimation of gene constraint from an evolutionary model with gene features

Measures of selective constraint on genes have been used for many applications including clinical interpretation of rare coding variants, disease gene discovery, and studies of genome evolution. However, widely-used metrics are severely underpowered at detecting constraint for the shortest ∼25% of genes, potentially causing important pathogenic mutations to be overlooked. We developed a framework combining a population genetics model with machine learning on gene features to enable accurate inference of an interpretable constraint metric, shet. Our estimates outperform existing metrics for prioritizing genes important for cell essentiality, human disease, and other phenotypes, especially for short genes. Our new estimates of selective constraint should have wide utility for characterizing genes relevant to human disease. Finally, our inference framework, GeneBayes, provides a flexible platform that can improve estimation of many gene-level properties, such as rare variant burden or gene expression differences.

OVCH1 Antisense RNA 1 is differentially expressed between non-frail and frail old adults

While some old adults stay healthy and non-frail up to late in life, others experience multimorbidity and frailty often accompanied by a pro-inflammatory state. The underlying molecular mechanisms for those differences are still obscure. Here, we used gene expression analysis to understand the molecular underpinning between non-frail and frail individuals in old age. Twenty-four adults (50% non-frail and 50% frail) from InCHIANTI study were included. Total RNA extracted from whole blood was analyzed by Cap Analysis of Gene Expression (CAGE). CAGE identified transcription start site (TSS) and active enhancer regions. We identified a set of differentially expressed (DE) TSS and enhancer between non-frail and frail and male and female participants. Several DE TSSs were annotated as lncRNA (XIST and TTTY14) and antisense RNAs (ZFX-AS1 and OVCH1 Antisense RNA 1). The promoter region chr6:366,786,54-366,787,97;+ was DE and overlapping the longevity CDKN1A gene. GWAS-LD enrichment analysis identifies overlapping LD-blocks with the DE regions with reported traits in GWAS catalog (isovolumetric relaxation time and urinary tract infection frequency). Furthermore, we used weighted gene co-expression network analysis (WGCNA) to identify changes of gene expression associated with clinical traits and identify key gene modules. We performed functional enrichment analysis of the gene modules with significant trait/module correlation. One gene module is showing a very distinct pattern in hub genes. Glycogen Phosphorylase L (PYGL) was the top ranked hub gene between non-frail and frail. We predicted transcription factor binding sites (TFBS) and motif activity. TF involved in age-related pathways (e.g., FOXO3 and MYC) shows different expression patterns between non-frail and frail participants. Expanding the study of OVCH1 Antisense RNA 1 and PYGL may help understand the mechanisms leading to loss of homeostasis that ultimately causes frailty.

Progress and limitations in engineering cellular adhesion for research and therapeutics

Intercellular interactions form the cornerstone of multicellular biology. Despite advances in protein engineering, researchers artificially directing physical cell interactions still rely on endogenous cell adhesion molecules (CAMs) alongside off-target interactions and unintended signaling. Recently, methods for directing cellular interactions have been developed utilizing programmable domains such as coiled coils (CCs), nanobody-antigen, and single-stranded DNA (ssDNA). We first discuss desirable molecular- and systems-level properties in engineered CAMs, using the helixCAM platform as a benchmark. Next, we propose applications for engineered CAMs in immunology, developmental biology, tissue engineering, and neuroscience. Biologists in various fields can readily adapt current engineered CAMs to establish control over cell interactions, and their utilization in basic and translational research will incentivize further expansion in engineered CAM capabilities.

Unique Role of Vimentin in the Intermediate Filament Proteins Family

Intermediate filaments (IFs), being traditionally the least studied component of the cytoskeleton, have begun to receive more attention in recent years. IFs are found in different cell types and are specific to them. Accumulated data have shifted the paradigm about the role of IFs as structures that merely provide mechanical strength to the cell. In addition to this role, IFs have been shown to participate in maintaining cell shape and strengthening cell adhesion. The data have also been obtained that point out to the role of IFs in a number of other biological processes, including organization of microtubules and microfilaments, regulation of nuclear structure and activity, cell cycle control, and regulation of signal transduction pathways. They are also actively involved in the regulation of several aspects of intracellular transport. Among the intermediate filament proteins, vimentin is of particular interest for researchers. Vimentin has been shown to be associated with a range of diseases, including cancer, cataracts, Crohn's disease, rheumatoid arthritis, and HIV. In this review, we focus almost exclusively on vimentin and the currently known functions of vimentin intermediate filaments (VIFs). This is due to the structural features of vimentin, biological functions of its domains, and its involvement in the regulation of a wide range of basic cellular functions, and its role in the development of human diseases. Particular attention in the review will be paid to comparing the role of VIFs with the role of intermediate filaments consisting of other proteins in cell physiology.

AIRE relies on Z-DNA to flag gene targets for thymic T cell tolerization

AIRE is an unconventional transcription factor that enhances the expression of thousands of genes in medullary thymic epithelial cells and promotes clonal deletion or phenotypic diversion of self-reactive T cells1-4. The biological logic of AIRE's target specificity remains largely unclear as, in contrast to many transcription factors, it does not bind to a particular DNA sequence motif. Here we implemented two orthogonal approaches to investigate AIRE's cis-regulatory mechanisms: construction of a convolutional neural network and leveraging natural genetic variation through analysis of F1 hybrid mice5. Both approaches nominated Z-DNA and NFE2-MAF as putative positive influences on AIRE's target choices. Genome-wide mapping studies revealed that Z-DNA-forming and NFE2L2-binding motifs were positively associated with the inherent ability of a gene's promoter to generate DNA double-stranded breaks, and promoters showing strong double-stranded break generation were more likely to enter a poised state with accessible chromatin and already-assembled transcriptional machinery. Consequently, AIRE preferentially targets genes with poised promoters. We propose a model in which Z-DNA anchors the AIRE-mediated transcriptional program by enhancing double-stranded break generation and promoter poising. Beyond resolving a long-standing mechanistic conundrum, these findings suggest routes for manipulating T cell tolerance.

C19ORF84 connects piRNA and DNA methylation machineries to defend the mammalian germ line

In the male mouse germ line, PIWI-interacting RNAs (piRNAs), bound by the PIWI protein MIWI2 (PIWIL4), guide DNA methylation of young active transposons through SPOCD1. However, the underlying mechanisms of SPOCD1-mediated piRNA-directed transposon methylation and whether this pathway functions to protect the human germ line remain unknown. We identified loss-of-function variants in human SPOCD1 that cause defective transposon silencing and male infertility. Through the analysis of these pathogenic alleles, we discovered that the uncharacterized protein C19ORF84 interacts with SPOCD1. DNMT3C, the DNA methyltransferase responsible for transposon methylation, associates with SPOCD1 and C19ORF84 in fetal gonocytes. Furthermore, C19ORF84 is essential for piRNA-directed DNA methylation and male mouse fertility. Finally, C19ORF84 mediates the in vivo association of SPOCD1 with the de novo methylation machinery. In summary, we have discovered a conserved role for the human piRNA pathway in transposon silencing and C19ORF84, an uncharacterized protein essential for orchestrating piRNA-directed DNA methylation.