Distal Gene Research Articles

Chromatin profiling identifies putative dual roles for H3K27me3 in regulating transposons and cell type-specific genes in choanoflagellates.

Gene expression is tightly controlled during animal development to allow the formation of specialized cell types. Our understanding of how animals evolved this exquisite regulatory control remains elusive, but evidence suggests that changes in chromatin-based mechanisms may have contributed. To investigate this possibility, here we examine chromatin-based gene regulatory features in the closest relatives of animals, choanoflagellates. Using Salpingoeca rosetta as a model system, we examined chromatin accessibility and histone modifications at the genome scale and compared these features to gene expression. We first observed that accessible regions of chromatin are primarily associated with gene promoters and found no evidence of distal gene regulatory elements resembling the enhancers that animals deploy to regulate developmental gene expression. Remarkably, a histone modification deposited by polycomb repressive complex 2, histone H3 lysine 27 trimethylation (H3K27me3), appeared to function similarly in S. rosetta to its role in animals, because this modification decorated genes with cell type-specific expression. Additionally, H3K27me3 marked transposons, retaining what appears to be an ancestral role in regulating these elements. We further uncovered a putative new bivalent chromatin state at cell type-specific genes that consists of H3K27me3 and histone H3 lysine 4 mono-methylation (H3K4me1). Together, our discoveries support the scenario that gene-associated histone modification states that underpin development emerged before the evolution of animal multicellularity.

Multiple regulatory inputs including cell envelope stress orchestrate expression of the Escherichia coli rpoN operon.

The rpoN operon, an important regulatory hub in Enterobacteriaceae, includes rpoN encoding sigma factor σ54, hpf involved in ribosome hibernation, rapZ regulating glucosamine-6-phosphate levels, and two genes encoding proteins of the nitrogen-related phosphotransferase system. Little is known about regulatory mechanisms controlling the abundance of these proteins. This study employs transposon mutagenesis and chemical screens to dissect the complex expression of the rpoN operon. We find that envelope stress conditions trigger read-through transcription into the rpoN operon from a promoter located upstream of the preceding lptA-lptB locus. This promoter is controlled by the envelope stress sigma factor E and response regulator PhoP is required for its full response to a subset of stress signals. σE also stimulates ptsN-rapZ-npr expression using an element downstream of rpoN, presumably by interfering with mRNA processing by RNase E. Additionally, we identify a novel promoter in the 3' end of rpoN that directs transcription of the distal genes in response to ethanol. Finally, we show that translation of hpf and ptsN is individually regulated by the RNA chaperone Hfq, perhaps involving small RNAs. Collectively, our work demonstrates that the rpoN operon is subject to complex regulation, integrating signals related to envelope stress and carbon source quality.

Identification of Arginine-Vasopressin Receptor 1a (Avpr1a/Avpr1a) as a Novel Candidate Gene for Chronic Visceral Pain Sheds Light on the Potential Role of Enteric Neurons in the Development of Visceral Hypersensitivity

Chronic abdominal pain in the absence of ongoing disease is the hallmark of disorders of gut-brain interaction (DGBIs), including irritable bowel syndrome (IBS). While the etiology of DGBIs remains poorly understood, there is evidence that both genetic and environmental factors play a role. In this study, we report the identification and validation of arginine-vasopressin receptor 1A (Avpr1a) as a novel candidate gene for visceral hypersensitivity (VH), a primary peripheral mechanism underlying abdominal pain in DGBI/IBS. Comparing 2 C57BL/6 (BL/6) substrains (C57BL/6NTac and C57BL/6J) revealed differential susceptibility to the development of chronic VH following intrarectal zymosan instillation, a validated preclinical model for postinflammatory IBS. Using whole-genome sequencing, we identified a single-nucleotide polymorphism differentiating the 2 strains in the 5′ intergenic region upstream of Avpr1a, encoding the protein Avpr1a. We used behavioral, histological, and molecular approaches to identify distal colon-specific gene expression and neuronal hyperresponsiveness covarying with Avpr1a genotype and VH susceptibility. While the 2 BL/6 substrains did not differ across other gastrointestinal phenotypes (eg, fecal water retention), VH-susceptible BL/6NTac mice had higher colonic Avpr1a mRNA and protein expression. These results parallel findings that patients' colonic Avpr1a mRNA expression corresponded to higher pain ratings. Moreover, neurons of the enteric nervous system were hyperresponsive to the Avpr1a agonist arginine-vasopressin, suggesting a role for enteric neurons in the pathology underlying VH. Taken together, these findings implicate differential regulation of Avpr1a as a novel mechanism of VH susceptibility as well as a potential therapeutic target specific to VH. PerspectiveThis article presents evidence of Avpr1a as a novel candidate gene for VH in a mouse model of IBS. Avpr1a genotype and/or tissue-specific expression represents a potential biomarker for chronic abdominal pain susceptibility.

Potential global cis and trans regulation of lncRNAs in Saccharomyces cerevisiae subjected to ethanol stress

Long noncoding RNAs (lncRNAs) are regulatory RNAs. Saccharomyces cerevisiae strains transcribe hundreds of lncRNAs. LncRNAs can regulate the expression of adjacent genes (cis-regulation) or distant genes from lncRNAs (trans-regulation). Here, we analyzed the potential global cis and trans-regulation of lncRNAs of yeast subjected to ethanol stress. For potential cis regulation, for BMA641-A and S288C strains, we observed that most lncRNA-neighbor gene pairs increased the expression at a certain point followed by a decrease, and vice versa. Based on the transcriptome profile and triple helix prediction between lncRNAs and promoters of coding genes, we observed nine different ways of potential trans regulation that work in a strain-specific manner. Our data provide an initial landscape of potential cis and trans regulation in yeast, which seems to be strain-specific.

Abstract 151: Single Cell Transcriptional And Epigenetic Map Of Human Coronary Artery Disease In A Multi-ethnic Population

Identifying target genes that are causally linked to coronary artery disease (CAD) remains an elusive problem despite the plethora of CAD genome-wide association study (GWAS) risk loci. Here, we performed single nucleus multiomic (RNA + ATAC) sequencing in 44 human coronary arteries and bulk Hi-C in 16 coronaries from multi-ethnic participants. From sequencing 126,804 nuclei, we identified cell type specific transcriptional signatures, chromatin peaks, and gene regulatory networks through peak to gene linkage analysis. Single-cell quantitative trait loci (QTL) analysis identified cell-specific caQTL’s including cell-type specific caQTLs within CAD GWAS loci. The overlap of CAD GWAS loci with ATAC peaks demonstrated that vascular smooth muscle and endothelial cells harbor the greatest genetic risk. To identify proximal and distal putative target genes for CAD loci in a cell type specific manner, we built a enhancer-gene map by leveraging coronary artery ChiP-seq, cell specific RNA/ATAC-seq, and Hi-C based chromatin contact loops. Using this map, we identified distal target genes linked to CAD GWAS loci. Finally, given the effect of ancestry on CAD GWAS risk interpretation, we leveraged Hi-C data to compare the chromatin architecture in coronary tissue from 8 European Americans vs 8 African Americans. A meta-map, sample level, and local ancestry analysis identified ancestry-specific differences in chromatin looping and topologically associated domains. Collectively, our results provide a cell-type specific transcriptomic and epigenetic map which along with bulk chromatin contacts can be used to interpret and prioritize CAD GWAS candidates and nominate cell specific causal target genes.

A Bidirectional Non-Coding RNA Promoter Mediates Long-Range Gene Expression Regulation.

Recent evidence suggests that human gene promoters display gene expression regulatory mechanisms beyond the typical single gene local transcription modulation. In mammalian genomes, genes with an associated bidirectional promoter are abundant; bidirectional promoter architecture serves as a regulatory hub for a gene pair expression. However, it has been suggested that its contribution to transcriptional regulation might exceed local transcription initiation modulation. Despite their abundance, the functional consequences of bidirectional promoter architecture remain largely unexplored. This work studies the long-range gene expression regulatory role of a long non-coding RNA gene promoter using chromosome conformation capture methods. We found that this particular bidirectional promoter contributes to distal gene expression regulation in a target-specific manner by establishing promoter-promoter interactions. In particular, we validated that the promoter-promoter interactions of this regulatory element with the promoter of distal gene BBX contribute to modulating the transcription rate of this gene; removing the bidirectional promoter from its genomic context leads to a rearrangement of BBX promoter-enhancer interactions and to increased gene expression. Moreover, long-range regulatory functionality is not directly dependent on its associated non-coding gene pair expression levels.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BackgroundRegulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.ResultsHere we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes.We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.ConclusionBased on these associations of (linearly) distal genes’ transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

A computational approach for deciphering the interactions between proximal and distal gene regulators in GC B-cell response.

Delineating the intricate interplay between promoter-proximal and -distal regulators is crucial for understanding the function of transcriptional mediator complexes implicated in the regulation of gene expression. The present study aimed to develop a computational method for accurately modeling the spatial proximal and distal regulatory interactions. Our method combined regression-based models to identify key regulators through gene expression prediction and a graph-embedding approach to detect coregulated genes. This approach enabled a detailed investigation of the gene regulatory mechanisms for germinal center B cells, accompanied by dramatic rearrangements of the genome structure. We found that while the promoter-proximal regulatory elements were the principal regulators of gene expression, the distal regulators fine-tuned transcription. Moreover, our approach unveiled the presence of modular regulators, such as cofactors and proximal/distal transcription factors, which were co-expressed with their target genes. Some of these modules exhibited abnormal expression patterns in lymphoma. These findings suggest that the dysregulation of interactions between transcriptional and architectural factors is associated with chromatin reorganization failure, which may increase the risk of malignancy. Therefore,our computational approachhelps decipherthe transcriptional cis-regulatory code spatially interacting.

P53 rapidly restructures 3D chromatin organization to trigger a transcriptional response

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

Osteoporosis GWAS-implicated DNM3 locus contextually regulates osteoblastic and chondrogenic fate of mesenchymal stem/progenitor cells through oscillating miR-199a-5p levels.

Genome wide association study (GWAS)-implicated bone mineral density (BMD) signals have been shown to localize in cis-regulatory regions of distal effector genes using 3D genomic methods. Detailed characterization of such genes can reveal novel causal genes for BMD determination. Here, we elected to characterize the "DNM3" locus on chr1q24, where the long non-coding RNA DNM3OS and the embedded microRNA MIR199A2 (miR-199a-5p) are implicated as effector genes contacted by the region harboring variation in linkage disequilibrium with BMD-associated sentinel single nucleotide polymorphism, rs12041600. During osteoblast differentiation of human mesenchymal stem/progenitor cells (hMSC), miR-199a-5p expression was temporally decreased and correlated with the induction of osteoblastic transcription factors RUNX2 and Osterix. Functional relevance of miR-199a-5p downregulation in osteoblastogenesis was investigated by introducing miR-199a-5p mimic into hMSC. Cells overexpressing miR-199a-5p depicted a cobblestone-like morphological change and failed to produce BMP2-dependent extracellular matrix mineralization. Mechanistically, a miR-199a-5p mimic modified hMSC propagated normal SMAD1/5/9 signaling and expressed osteoblastic transcription factors RUNX2 and Osterix but depicted pronounced upregulation of SOX9 and enhanced expression of essential chondrogenic genes ACAN, COMP, and COL10A1. Mineralization defects, morphological changes, and enhanced chondrogenic gene expression associated with miR-199a-5p mimic over-expression were restored with miR-199a-5p inhibitor suggesting specificity of miR-199a-5p in chondrogenic fate specification. The expression of both the DNM3OS and miR-199a-5p temporally increased and correlated with hMSC chondrogenic differentiation. Although miR-199a-5p overexpression failed to further enhance chondrogenesis, blocking miR-199a-5p activity significantly reduced chondrogenic pellet size, extracellular matrix deposition, and chondrogenic gene expression. Taken together, our results indicate that oscillating miR-199a-5p levels dictate hMSC osteoblast or chondrocyte terminal fate. Our study highlights a functional role of miR-199a-5p as a BMD effector gene at the DNM3 BMD GWAS locus, where patients with cis-regulatory genetic variation which increases miR-199a-5p expression could lead to reduced osteoblast activity.

No safe haven: Loss of avirulence in the plant pathogen Leptosphaeria maculans by DNA duplication and repeat‐induced point mutation

AbstractMicrobes can overcome the ability of plant resistance genes to confer protection against disease by mutating their corresponding avirulence genes. The fungus Leptosphaeria maculans causes blackleg disease on canola crops and subverts Brassica napus resistance genes through several DNA mutation mechanisms. One of these is repeat‐induced point (RIP) mutation, which can ‘leak’ into the avirulence genes from the adjacent repetitive sequences that the mutation process is targeting. Here, we identified populations of L. maculans in Australia that have extensive RIP mutations in the avirulence gene AvrLm2 and show that this has been triggered by a duplication of the gene and surrounding DNA that includes the distant (>55 kb in total) AvrLm6 gene. This finding provides another mechanism of mutation by which pathogens can overcome host resistance, and more broadly contributes to understanding the complex balance between gene duplication and genome defence.

The Role of Mechanotransduction in Contact Inhibition of Locomotion and Proliferation.

Contact inhibition (CI) represents a crucial tumor-suppressive mechanism responsible for controlling the unbridled growth of cells, thus preventing the formation of cancerous tissues. CI can be further categorized into two distinct yet interrelated components: CI of locomotion (CIL) and CI of proliferation (CIP). These two components of CI have historically been viewed as separate processes, but emerging research suggests that they may be regulated by both distinct and shared pathways. Specifically, recent studies have indicated that both CIP and CIL utilize mechanotransduction pathways, a process that involves cells sensing and responding to mechanical forces. This review article describes the role of mechanotransduction in CI, shedding light on how mechanical forces regulate CIL and CIP. Emphasis is placed on filamin A (FLNA)-mediated mechanotransduction, elucidating how FLNA senses mechanical forces and translates them into crucial biochemical signals that regulate cell locomotion and proliferation. In addition to FLNA, trans-acting factors (TAFs), which are proteins or regulatory RNAs capable of directly or indirectly binding to specific DNA sequences in distant genes to regulate gene expression, emerge as sensitive players in both the mechanotransduction and signaling pathways of CI. This article presents methods for identifying these TAF proteins and profiling the associated changes in chromatin structure, offering valuable insights into CI and other biological functions mediated by mechanotransduction. Finally, it addresses unanswered research questions in these fields and delineates their possible future directions.

Machine learning-based prediction model for distant metastasis of breast cancer

BackgroundBreast cancer is the most prevalent malignancy in women. Advanced breast cancer can develop distant metastases, posing a severe threat to the life of patients. Because the clinical warning signs of distant metastasis are manifested in the late stage of the disease, there is a need for better methods of predicting metastasis. MethodsFirst, we screened breast cancer distant metastasis target genes by performing difference analysis and weighted gene co-expression network analysis (WGCNA) on the selected datasets, and performed analyses such as GO enrichment analysis on these target genes. Secondly, we screened breast cancer distant metastasis target genes by LASSO regression analysis and performed correlation analysis and other analyses on these biomarkers. Finally, we constructed several breast cancer distant metastasis prediction models based on Logistic Regression (LR) model, Random Forest (RF) model, Support Vector Machine (SVM) model, Gradient Boosting Decision Tree (GBDT) model and eXtreme Gradient Boosting (XGBoost) model, and selected the optimal model from them. ResultsSeveral 21-gene breast cancer distant metastasis prediction models were constructed, with the best performance of the model constructed based on the random forest model. This model accurately predicted the emergence of distant metastases from breast cancer, with an accuracy of 93.6 %, an F1-score of 88.9 % and an AUC value of 91.3 % on the validation set. ConclusionOur findings have the potential to be translated into a point-of-care prognostic analysis to reduce breast cancer mortality.

Target genes regulated by CLEC16A intronic region associated with common variable immunodeficiency

BackgroundCLEC16A intron 19 has been identified as a candidate locus for common variable immunodeficiency (CVID). ObjectiveThe objective of this study is to elucidate the molecular mechanism by which variants at the CLEC16A intronic locus may contribute to the pathogenesis of CVID. MethodsWe performed fine-mapping of the CLEC16A locus in a CVID cohort, then deleted the candidate functional SNP in T-cell lines by the CRISPR-Cas9 technique and conducted RNA-seq to identify target gene(s). The interactions between the CLEC16A locus and its target genes were identified using Circularized Chromosome Conformation Capture. The transcription factor complexes mediating the chromatin interactions were determined by proteomic approach. The molecular pathways regulated by the CLEC16A locus were examined by RNA-seq and reverse phase protein array. ResultsOur study showed that the CLEC16A locus is an enhancer regulating expression of multiple target genes including a distant gene ATF7IP2 through chromatin interactions. Distinct transcription factor complexes mediate the chromatin interactions in an allele-specific manner. Disruption of the CLEC16A locus affects the AKT signaling pathway, as well as the molecular response of CD4+T cells to immune stimulation. ConclusionThrough multi-omics and targeted experimental approaches, we elucidated the underlying target genes and signaling pathways involved in the genetic association of CLEC16A with CVID, and highlighted plausible molecular targets for developing novel therapeutics.

Cell-specific housekeeping role of lncRNAs in COVID-19-infected and recovered patients.

A plethora of studies have demonstrated the roles of lncRNAs in modulating disease severity and outcomes during infection. However, the spatio-temporal expression of these lncRNAs is poorly understood. In this study, we used single-cell RNA-seq to understand the spatio-temporal expression dynamics of lncRNAs across healthy, SARS-CoV-2-infected, and recovered individuals and their functional role in modulating the disease and recovery. We identified 203 differentially expressed lncRNAs, including cell type-specific ones like MALAT1, NEAT1, ZFAS1, SNHG7, SNHG8, and SNHG25 modulating immune function in classical monocyte, NK T, proliferating NK, plasmablast, naive,and activated B/T cells. Interestingly, we found invariant lncRNAs (no significant change in expression across conditions) regulating essential housekeeping functions (for example, HOTAIR, NRAV, SNHG27, SNHG28, and UCA1) in infected and recovered individuals. Despite similar repeat element abundance, variant lncRNAs displayed higher Alu content, suggesting increased interactions with proximal and distal genes, crucial for immune response modulation. The comparable repeat abundance but distinct expression levels of variant and invariant lncRNAs highlight the significance of investigating the regulatory mechanisms of invariant lncRNAs. Overall, this study offers new insights into the spatio-temporal expression patterns and functional roles of lncRNAs in SARS-CoV-2-infected and recovered individuals while highlighting the importance of invariant lncRNAs in the disease context.

A refactored biosynthetic pathway for the production of glycosylated microbial sunscreens.

Mycosporine-like amino acids (MAAs) are a family of water-soluble and colorless secondary metabolites, with high extinction coefficients, that function as microbial sunscreens. MAAs share a cyclohexinimine chromophore that is diversified through amino acid substitutions and attachment of sugar moieties. The genetic and enzymatic bases for the chemical diversity of MAAs remain largely unexplored. Here we report a series of structurally distinct MAAs and evidence for an unusual branched biosynthetic pathway from a cyanobacterium isolated from lake sediment. We used a combination of high-resolution liquid chromatography-mass spectrometry (HR-LCMS) analysis and nuclear magnetic resonance (NMR) spectroscopy to identify diglycosylated-palythine-Ser (C22H36N2O15) as the dominant chemical variant in a series of MAAs from Nostoc sp. UHCC 0302 that contained either Ser or Thr. We obtained a complete 9.9 Mb genome sequence to gain insights into the genetic basis for the biosynthesis of these structurally distinct MAAs. We identified MAA biosynthetic genes encoded at two locations on the circular chromosome. Surprisingly, direct pathway cloning and heterologous expression of the complete mysABCJ 1 D 1 G 1 H biosynthetic gene cluster in Escherichia coli (E. coli) led to the production of 450 Da monoglycosylated-palythine-Thr (C18H30N2O11). We reconstructed combinations of the two distant biosynthetic gene clusters in refactored synthetic pathways and expressed them in the heterologous host. These results demonstrated that the MysD1 and MysD2 enzymes displayed a preference for Thr and Ser, respectively. Furthermore, one of the four glycosyltransferases identified, MysG1, was active in E. coli and catalysed the attachment of a hexose moiety to the palythine-Thr intermediate. Together these results provide the first insights into the enzymatic basis for glycosylation of MAAs and demonstrates how paralogous copies of the MysD enzymes allow the simultaneous biosynthesis of specific chemical variants to increase the structural variation in this family of microbial sunscreens.

Genetically corrected RAG2-SCID human hematopoietic stem cells restore V(D)J-recombinase and rescue lymphoid deficiency

AbstractRecombination-activating genes (RAG1 and RAG2) are critical for lymphoid cell development and function by initiating the variable (V), diversity (D), and joining (J) (V(D)J)-recombination process to generate polyclonal lymphocytes with broad antigen specificity. The clinical manifestations of defective RAG1/2 genes range from immune dysregulation to severe combined immunodeficiencies (SCIDs), causing life-threatening infections and death early in life without hematopoietic cell transplantation (HCT). Despite improvements, haploidentical HCT without myeloablative conditioning carries a high risk of graft failure and incomplete immune reconstitution. The RAG complex is only expressed during the G0-G1 phase of the cell cycle in the early stages of T- and B-cell development, underscoring that a direct gene correction might capture the precise temporal expression of the endogenous gene. Here, we report a feasibility study using the CRISPR/Cas9-based “universal gene-correction” approach for the RAG2 locus in human hematopoietic stem/progenitor cells (HSPCs) from healthy donors and RAG2-SCID patient. V(D)J-recombinase activity was restored after gene correction of RAG2-SCID-derived HSPCs, resulting in the development of T-cell receptor (TCR) αβ and γδ CD3+ cells and single-positive CD4+ and CD8+ lymphocytes. TCR repertoire analysis indicated a normal distribution of CDR3 length and preserved usage of the distal TRAV genes. We confirmed the in vivo rescue of B-cell development with normal immunoglobulin M surface expression and a significant decrease in CD56bright natural killer cells. Together, we provide specificity, toxicity, and efficacy data supporting the development of a gene-correction therapy to benefit RAG2-deficient patients.

A cis-regulatory element regulates ERAP2 expression through autoimmune disease risk SNPs

Single-nucleotide polymorphisms (SNPs) near the ERAP2 gene are associated with various autoimmune conditions, as well as protection against lethal infections. Due to high linkage disequilibrium, numerous trait-associated SNPs are correlated with ERAP2 expression; however, their functional mechanisms remain unidentified. We show by reciprocal allelic replacement that ERAP2 expression is directly controlled by the splice region variant rs2248374. However, disease-associated variants in the downstream LNPEP gene promoter are independently associated with ERAP2 expression. Allele-specific conformation capture assays revealed long-range chromatin contacts between the gene promoters of LNPEP and ERAP2 and showed that interactions were stronger in patients carrying the alleles that increase susceptibility to autoimmune diseases. Replacing the SNPs in the LNPEP promoter by reference sequences lowered ERAP2 expression. These findings show that multiple SNPs act in concert to regulate ERAP2 expression and that disease-associated variants can convert a gene promoter region into a potent enhancer of a distal gene.

Proximal and Distal Bronchioles Contribute to the Pathogenesis of Non-Cystic Fibrosis Bronchiectasis.

Rationale: Non-cystic fibrosis bronchiectasis (NCFB) may originate in bronchiolar regions of the lung. Accordingly, there is a need to characterize the morphology and molecular characteristics of NCFB bronchioles. Objectives: Test the hypothesis that NCFB exhibits a major component of bronchiolar disease manifest by mucus plugging and ectasia. Methods: Morphologic criteria and region-specific epithelial gene expression, measured histologically and by RNA in situ hybridization and immunohistochemistry, identified proximal and distal bronchioles in excised NCFB lungs. RNA in situ hybridization and immunohistochemistry assessed bronchiolar mucus accumulation and mucin gene expression. CRISPR-Cas9-mediated IL-1R1 knockout in human bronchial epithelial cultures tested IL-1α and IL-1β contributions to mucin production. Spatial transcriptional profiling characterized NCFB distal bronchiolar gene expression. Measurements and Main Results: Bronchiolar perimeters and lumen areas per section area were increased in proximal, but not distal, bronchioles in NCFB versus control lungs, suggesting proximal bronchiolectasis. In NCFB, mucus plugging was observed in ectatic proximal bronchioles and associated nonectatic distal bronchioles in sections with disease. MUC5AC and MUC5B mucins were upregulated in NCFB proximal bronchioles, whereas MUC5B was selectively upregulated in distal bronchioles. Bronchiolar mucus plugs were populated by IL-1β-expressing macrophages. NCFB sterile sputum supernatants induced human bronchial epithelial MUC5B and MUC5AC expression that was >80% blocked by IL-1R1 ablation. Spatial transcriptional profiling identified upregulation of genes associated with secretory cells, hypoxia, interleukin pathways, and IL-1β-producing macrophages in mucus plugs and downregulation of epithelial ciliogenesis genes. Conclusions: NCFB exhibits distinctive proximal and distal bronchiolar disease. Both bronchiolar regions exhibit bronchiolar secretory cell features and mucus plugging but differ in mucin gene regulation and ectasia.

Differential usage of DNA modifications in neurons, astrocytes, and microglia

BackgroundCellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation, DNA modifications in particular, of gene expression between neurons and glia.ResultsAfter validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT-whole genome oxidative bisulfite sequencing (WGoxBS) to assess the neuronal translatome and epigenome in the hippocampus of young mice (4 months old). WGoxBS findings were validated with enzymatic methyl-Seq (EM-Seq) and nanopore sequencing. Comparing neuronal data to microglial and astrocytic data from NuTRAP models, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, rather than proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of gene body mCG and a positive relationship between distal promoter and gene body hmCG with gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions.ConclusionsNeurons, astrocytes, and microglia demonstrate different genome-wide levels of mCG, hmCG, and mCH that are reproducible across analytical methods. However, modification-gene expression relationships are conserved across cell types. Enrichment of differential modifications across cell types in gene bodies and distal regulatory elements, but not proximal promoters, highlights epigenomic patterning in these regions as potentially greater determinants of cell identity. These findings also demonstrate the importance of differentiating between mC and hmC in neuroepigenomic analyses, as up to 30% of what is conventionally interpreted as mCG can be hmCG, which often has a different relationship to gene expression than mCG.