Characterization of the human RFX transcription factor family by regulatory and target gene analysis

Mapping the transcriptional and epigenetic landscape of organotypic endothelial diversity in the developing and adult mouse

The vascular endothelium features unique molecular and functional properties across different vessel types, such as between arteries, veins and capillaries, as well as between different organs, such as the leaky sinusoidal endothelium of the liver versus the impermeable vessels of the brain. However, the transcriptional networks governing endothelial organ specialization remain unclear. Here we profile the accessible chromatin and transcriptional landscapes of the endothelium from the mouse liver, lung, heart, kidney, brain and retina, across developmental time, to identify potential transcriptional regulators of endothelial heterogeneity. We then determine which of these putative regulators are conserved in human brain endothelial cells, and using single-cell transcriptomic profiling, we define which regulatory networks are active during brain maturation. Finally, we show that the putative transcriptional regulators identified by these three approaches molecularly and functionally reprogram naive endothelial cells. Thus, this resource can be used to identify potential transcriptional regulators controlling the establishment and maintenance of organ-specific endothelial specialization.

RFX2 downregulates RASSF1 expression and YAP phosphorylation through Hippo signaling to promote immune escape in lung adenocarcinoma

OBJECTIVE

Regulatory Factor X (RFX) transcription factors have been implicated in different cancers. Ras association domain family (RASSF) has been shown clinical significance in lung cancer. This paper was to investigate the interaction of RFX2 and RASSF1 in lung adenocarcinoma (LUAD).

METHODS

The transcriptome differences of LUAD patients in GSE32863, GSE43458, and GSE21933 datasets were analyzed. A-549 and NCI-H358 cell lines after overexpression of RFX2 were co-cultured with activated CD8+ T cells, and the release of IFN-γ, GZMB, PRF1 by CD8+ T cells, and PD-L1 in the LUAD cells were detected. Cell viability, invasion, and apoptosis were analyzed by CCK-8, Transwell, and TUNEL assays. Dual-luciferase assay and ChIP were conducted to detect the interaction between RFX2 and RASSF1 promoter. An in vivo tumor model was constructed to monitor tumor growth. YAP protein levels and phosphorylation were measured. A-549 and NCI-H358 cells treated with DMSO or PY-60 after RFX2 overexpression were co-cultured with activated CD8+ T cells.

RESULTS

RFX2 was notably downregulated in LUAD. RFX2 overexpression increased infiltrating CD8+ T cells within transplanted tumors and inhibited immune escape, proliferation, and invasion of LUAD cells. RFX2 was enriched in the RASSF1 promoter, and RFX2 activated RASSF1 transcription by binding to the RASSF1 promoter. RASSF1 knockdown reversed the ability of RFX2 overexpression to inhibit immune escape. RFX2 depletion downregulated RASSF1, which reduced YAP phosphorylation, thus affecting the Hippo pathway to promote the immune escape.

CONCLUSION

RFX2 Loss in LUAD downregulates RASSF1 expression and YAP phosphorylation, thereby promoting immune escape through Hippo signaling.

Multi-omic analysis of the ciliogenic transcription factor RFX3 reveals a role in promoting activity-dependent responses via enhancing CREB binding in human neurons

Heterozygous loss-of-function (LoF) variants in RFX3, a transcription factor known to play key roles in ciliogenesis, result in autism spectrum disorder (ASD) and neurodevelopmental delay. RFX binding motifs are also enriched upstream of genes found to be commonly dysregulated in transcriptomic analyses of brain tissue from individuals with idiopathic ASD. Still, the precise functions of RFX3 in the human brain is unknown. Here, we studied the impact of RFX3 deficiency using human iPSC-derived neurons and forebrain organoids. Biallelic loss of RFX3 disrupted ciliary gene expression and delayed neuronal differentiation, while monoallelic loss of RFX3 did not. Instead, transcriptomic and DNA binding analyses demonstrated that monoallelic RFX3 loss disrupted synaptic target gene expression and diminished neuronal activity-dependent gene expression. RFX3 binding sites co-localized with CREB binding sites near activity-dependent genes, and RFX3 deficiency led to decreased CREB binding and impaired induction of CREB targets in response to neuronal depolarization. This study demonstrates a novel role of the ASD-associated gene RFX3 in shaping neuronal synaptic development and plasticity.

Transcriptional enhancers in human neuronal differentiation provide clues to neuronal disorders

Genome-wide association studies (GWASs) have identified thousands of variants associated with complex phenotypes, including neuropsychiatric disorders. To better understand their pathogenesis, it is necessary to identify the functional roles of these variants, which are largely located in non-coding DNA regions. Here, we employ a human mesencephalic neuronal cell differentiation model, LUHMES, with sensitive and high-resolution methods to discover enhancers (NET-CAGE), perform DNA conformation analysis (Capture Hi-C) to link enhancers to their target genes, and finally validate selected interactions. We expand the number of known enhancers active in differentiating human LUHMES neurons to 47,350, and find overlap with GWAS variants for Parkinson's disease and schizophrenia. Our findings reveal a fine-tuned regulation of human neuronal differentiation, even between adjacent developmental stages; provide a valuable resource for further studies on neuronal development, regulation, and disorders; and emphasize the importance of exploring the vast regulatory potential of non-coding DNA and enhancers.

Suppression of regulatory factor X 7 alleviates airway remodeling and inflammation in childhood asthma

Objective

Childhood asthma is a chronic heterogeneous syndrome composed of distinct disease entities or phenotypes. This study was conducted to characterize regulatory factor X 7 (RFX7) in childhood asthma.

Material and Methods

Two available transcriptome datasets (GSE65204 and GSE27011) were used to analyze regulatory factor X (RFX) family members in childhood asthma. Random forest, logistic regression, and linear support vector machine (SVM) analyses were performed to construct an RFX-based classification model. Airway smooth muscle cells (ASMCs) were induced through platelet-derived growth factor-BB (PDGF-BB) for an asthma in vitro model. RFX7 expression was measured through immunoblotting. RFX7 was knocked out by transfection of RFX7 small-interfering RNAs, and then airway remodeling and inflammation were assayed.

Results

Among RFX family members, RFX3, RFX7, and RFX-associated protein displayed differential expression in childhood asthma versus healthy controls. Thus, SVM, logistic regression, and random forest-based machine learning models were built. The random forest model presented the best diagnostic efficacy (area under the curve [AUC] = 1 and 0.67 in discovery and verification sets). RFX7 was found to be effective in diagnosing childhood asthma (AUC = 0.724 and 0.775 in discovery and verification sets). In addition, RFX7 was overexpressed in PDGF-BB-stimulated ASMCs (✶ ✶ P < 0.01). Silencing RFX7 remarkably attenuated the proliferative and migrative capacities of ASMCs with PDGF-BB stimulation (✶ ✶ P < 0.01). In addition, RFX7 was positively related to neutrophil infiltration in childhood asthma, and its knockdown downregulated the levels of pro-inflammatory cytokines in PDGF-BB-stimulated ASMCs (✶ ✶ P < 0.01).

Conclusion

The findings of this study indicate that RFX7 is a novel molecule that is correlated with airway remodeling and inflammation in childhood asthma, providing insights into the mechanism underlying this disease and its potential clinical importance.

Expanding the clinical phenotype and variant spectrum associated with RFX7

RFX7 encodes a transcription factor that is ubiquitously expressed and important for neural development. Haploinsufficiency of RFX7 is associated with intellectual disability, developmental delay, and diverse malformations of brain structures. Currently, there are only 16 clinically described individuals who have variants in RFX7. A recognizable pattern of malformation associated with mutation in RFX7 has not yet been uncovered. Here we describe the phenotypic presentation of two additional individuals who have novel de novo variants in RFX7. One of the individuals we describe is from an under-represented Afro-Caribbean population.

Single-cell analysis of chromatin and expression reveals age- and sex-associated alterations in the human heart

Sex differences and age-related changes in the human heart at the tissue, cell, and molecular level have been well-documented and many may be relevant for cardiovascular disease. However, how molecular programs within individual cell types vary across individuals by age and sex remains poorly characterized. To better understand this variation, we performed single-nucleus combinatorial indexing (sci) ATAC- and RNA-Seq in human heart samples from nine donors. We identify hundreds of differentially expressed genes by age and sex and find epigenetic signatures of variation in ATAC-Seq data in this discovery cohort. We then scale up our single-cell RNA-Seq analysis by combining our data with five recently published single nucleus RNA-Seq datasets of healthy adult hearts. We find variation such as metabolic alterations by sex and immune changes by age in differential expression tests, as well as alterations in abundance of cardiomyocytes by sex and neurons with age. In addition, we compare our adult-derived ATAC-Seq profiles to analogous fetal cell types to identify putative developmental-stage-specific regulatory factors. Finally, we train predictive models of cell-type-specific RNA expression levels utilizing ATAC-Seq profiles to link distal regulatory sequences to promoters, quantifying the predictive value of a simple TF-to-expression regulatory grammar and identifying cell-type-specific TFs. Our analysis represents the largest single-cell analysis of cardiac variation by age and sex to date and provides a resource for further study of healthy cardiac variation and transcriptional regulation at single-cell resolution.

p53 target ANKRA2 cooperates with RFX7 to regulate tumor suppressor genes

The transcription factor regulatory factor X 7 (RFX7) has been identified as a tumor suppressor that is recurrently mutated in lymphoid cancers and appears to be dysregulated in many other cancers. RFX7 is activated by the well-known tumor suppressor p53 and regulates several other known tumor suppressor genes. However, what other factors regulate RFX7 and its target genes remains unclear. Here, reporter gene assays were used to identify that RFX7 regulates the tumor suppressor gene PDCD4 through direct interaction with its X-box promoter motif. We utilized mass spectrometry to identify factors that bind to DNA together with RFX7. In addition to RFX7, we also identified RFX5, RFXAP, RFXANK, and ANKRA2 that bind to the X-box motif in the PDCD4 promoter. We demonstrate that ANKRA2 is a bona fide direct p53 target gene. We used transcriptome analyses in two cell systems to identify genes regulated by ANKRA2, its sibling RFXANK, and RFX7. These results revealed that ANKRA2 functions as a critical cofactor of RFX7, whereas RFXANK regulates largely distinct gene sets.

Pan-cancer investigation of RFX family and associated genes identifies RFX8 as a therapeutic target in leukemia

Background

Several transcription factors and co-factors are encoded by the RFX (Regulatory Factor X) family (RFX1-8) and associated genes (RFXAP and RFXANK). Increasing evidence suggests that the RFX family and associated genes are involved in the development and progression of cancer. However, no prior research has focused on a multi-omic analysis of these genes to evaluate their role in tumor progression.

Methods

Using combined TCGA and GTEx pan-cancer data, we investigated the expression patterns and survival profiles of these ten genes. We then focused on RFX8 to analyze its clinicopathological and therapeutic features. Finally, we conducted experimental validation of RFX8 function in acute myeloid leukemia (AML).

Results

RFX5 and RFXANK showed higher expression levels, while RFX6 showed lower expression levels in most types of cancer, with RFX8 being the most upregulated in LAML. RFX2 and RFXAP demonstrated prognostic significance in eight types of cancer, and RFX8 showed significance in six types of cancer. The expression of these ten genes exhibited specific characteristics in immune subtypes, tumor microenvironment, and stemness. The expression of RFX8 was correlated with various tumor stages, microsatellite instability (MSI), tumor mutation burden (TMB), immune cell infiltration, and immune-checkpoint expression. Additionally, RFX8 was found to regulate tumorigenesis and sensitivity to chelerythrine in AML.

Conclusions

Our work delineated the landscape of the RFX family and associated genes in the pan-cancer context and the specific role of RFX8 in AML. These findings might offer cues for further investigations of these genes in cancer biology.

An Inherited Allele Confers Prostate Cancer Progression and Drug Resistance via RFX6/HOXA10-Orchestrated TGFβ Signaling

Genetic and epigenetic alterations are cancer hallmark characteristics. However, the role of inherited cancer predisposition alleles in co-opting lineage factor epigenetic reprogramming and tumor progression remains elusive. Here the FinnGen cohort phenome-wide analysis, along with multiple genome-wide association studies, has consistently identified the rs339331-RFX6/6q22 locus associated with prostate cancer (PCa) risk across diverse populations. It is uncovered that rs339331 resides in a reprogrammed androgen receptor (AR) binding site in PCa tumors, with the T risk allele enhancing AR chromatin occupancy. RFX6, an AR-regulated gene linked to rs339331, exhibits synergistic prognostic value for PCa recurrence and metastasis. This comprehensive in vitro and in vivo studies demonstrate the oncogenic functions of RFX6 in promoting PCa cell proliferation and metastasis. Mechanistically, RFX6 upregulates HOXA10 that profoundly correlates with adverse PCa outcomes and is pivotal in RFX6-mediated PCa progression, facilitating the epithelial-mesenchymal transition (EMT) and modulating the TGFβ/SMAD signaling axis. Clinically, HOXA10 elevation is associated with increased EMT scores, tumor advancement and PCa recurrence. Remarkably, reducing RFX6 expression restores enzalutamide sensitivity in resistant PCa cells and tumors. This findings reveal a complex interplay of genetic and epigenetic mechanisms in PCa pathogenesis and drug resistance, centered around disrupted prostate lineage AR signaling and abnormal RFX6 expression.

RFX6 haploinsufficiency predisposes to diabetes through impaired beta cell function

AIMS/HYPOTHESIS

Regulatory factor X 6 (RFX6) is crucial for pancreatic endocrine development and differentiation. The RFX6 variant p.His293LeufsTer7 is significantly enriched in the Finnish population, with almost 1:250 individuals as a carrier. Importantly, the FinnGen study indicates a high predisposition for heterozygous carriers to develop type 2 and gestational diabetes. However, the precise mechanism of this predisposition remains unknown.

METHODS

To understand the role of this variant in beta cell development and function, we used CRISPR technology to generate allelic series of pluripotent stem cells. We created two isogenic stem cell models: a human embryonic stem cell model; and a patient-derived stem cell model. Both were differentiated into pancreatic islet lineages (stem-cell-derived islets, SC-islets), followed by implantation in immunocompromised NOD-SCID-Gamma mice.

RESULTS

Stem cell models of the homozygous variant RFX6-/- predictably failed to generate insulin-secreting pancreatic beta cells, mirroring the phenotype observed in Mitchell-Riley syndrome. Notably, at the pancreatic endocrine stage, there was an upregulation of precursor markers NEUROG3 and SOX9, accompanied by increased apoptosis. Intriguingly, heterozygous RFX6+/- SC-islets exhibited RFX6 haploinsufficiency (54.2% reduction in protein expression), associated with reduced beta cell maturation markers, altered calcium signalling and impaired insulin secretion (62% and 54% reduction in basal and high glucose conditions, respectively). However, RFX6 haploinsufficiency did not have an impact on beta cell number or insulin content. The reduced insulin secretion persisted after in vivo implantation in mice, aligning with the increased risk of variant carriers to develop diabetes.

CONCLUSIONS/INTERPRETATION

Our allelic series isogenic SC-islet models represent a powerful tool to elucidate specific aetiologies of diabetes in humans, enabling the sensitive detection of aberrations in both beta cell development and function. We highlight the critical role of RFX6 in augmenting and maintaining the pancreatic progenitor pool, with an endocrine roadblock and increased cell death upon its loss. We demonstrate that RFX6 haploinsufficiency does not affect beta cell number or insulin content but does impair function, predisposing heterozygous carriers of loss-of-function variants to diabetes.

DATA AVAILABILITY

Ultra-deep bulk RNA-seq data for pancreatic differentiation stages 3, 5 and 7 of H1 RFX6 genotypes are deposited in the Gene Expression Omnibus database with accession code GSE234289. Original western blot images are deposited at Mendeley ( https://data.mendeley.com/datasets/g75drr3mgw/2 ).

LINC00908 attenuates LUAD tumorigenesis through DEAD-box helicase 54

Lung adenocarcinoma (LUAD) is one of the leading causes of cancer-related death worldwide. We identified a specific long non-coding RNA (LncRNA), LINC00908, which was downregulated in LUAD tissues and associated with good outcome. LINC00908 inhibited glycolysis by regulating the expression of the DEAD-box helicase 54 (DDX54), which was screened by a nine-gene risk signature, where DDX54 showed a positive correlation with several glycolysis-related genes. Experimental verification confirmed that DDX54 regulated nine key glycolytic enzymes, thereby affecting the level of glycolysis in LUAD. Further, the expression of LINC00908 in LUAD tumorigenesis was modulated by a transcription factor, regulatory factor X2 (RFX2). The RFX2/LINC00908/DDX54 axis regulated LUAD tumor growth, migration, invasion, cell apoptosis and glycolysis both in vitro and in vivo. These results demonstrate that this axis may serve as a novel mediator in LUAD progress and offer a novel therapeutic target for more precise diagnosis and treatment of LUAD.

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.

RFX transcription factor in the human-associated yeast Candida albicans regulates adhesion to oral epithelium

Adhesion to mucosal surfaces is a critical step in many bacterial and fungal infections. Here, using a mouse model of oral infection by the human fungal pathobiont Candida albicans, we report the identification of a novel regulator of C. albicans adhesion to the oral mucosa. The regulator is a member of the regulatory factor X (RFX) family of transcription factors, which control cellular processes ranging from genome integrity in model yeasts to tissue differentiation in vertebrates. Mice infected with the C. albicans rfx1 deletion mutant displayed increased fungal burden in tongues compared to animals infected with the reference strain. High-resolution imaging revealed RFX1 transcripts being expressed by C. albicans cells during infection. Concomitant with the increase in fungal burden, the rfx1 mutant elicited an enhanced innate immune response. Transcriptome analyses uncovered HWP1, a gene encoding an adhesion protein that mediates covalent attachment to buccal cells, as a major RFX1-regulated locus. Consistent with this result, we establish that C. albicans adhesion to oral cells is modulated by RFX1 in an HWP1-dependent manner. Our findings expand the repertoire of biological processes controlled by the RFX family and illustrate a mechanism whereby C. albicans can adjust adhesion to the oral epithelium.

Primary cilia promote the differentiation of human neurons through the WNT signaling pathway

BACKGROUND

Primary cilia emanate from most human cell types, including neurons. Cilia are important for communicating with the cell's immediate environment: signal reception and transduction to/from the ciliated cell. Deregulation of ciliary signaling can lead to ciliopathies and certain neurodevelopmental disorders. In the developing brain cilia play well-documented roles for the expansion of the neural progenitor cell pool, while information about the roles of cilia during post-mitotic neuron differentiation and maturation is scarce.

RESULTS

We employed ciliated Lund Human Mesencephalic (LUHMES) cells in time course experiments to assess the impact of ciliary signaling on neuron differentiation. By comparing ciliated and non-ciliated neuronal precursor cells and neurons in wild type and in RFX2 -/- mutant neurons with altered cilia, we discovered an early-differentiation "ciliary time window" during which transient cilia promote axon outgrowth, branching and arborization. Experiments in neurons with IFT88 and IFT172 ciliary gene knockdowns, leading to shorter cilia, confirm these results. Cilia promote neuron differentiation by tipping WNT signaling toward the non-canonical pathway, in turn activating WNT pathway output genes implicated in cyto-architectural changes.

CONCLUSIONS

We provide a mechanistic entry point into when and how ciliary signaling coordinates, promotes and translates into anatomical changes. We hypothesize that ciliary alterations causing neuron differentiation defects may result in "mild" impairments of brain development, possibly underpinning certain aspects of neurodevelopmental disorders.

RFX4 is an intrinsic factor for neuronal differentiation through induction of proneural genes POU3F2 and NEUROD1

Proneural genes play a crucial role in neuronal differentiation. However, our understanding of the regulatory mechanisms governing proneural genes during neuronal differentiation remains limited. RFX4, identified as a candidate regulator of proneural genes, has been reported to be associated with the development of neuropsychiatric disorders. To uncover the regulatory relationship, we utilized a combination of multi-omics data, including ATAC-seq, ChIP-seq, Hi-C, and RNA-seq, to identify RFX4 as an upstream regulator of proneural genes. We further validated the role of RFX4 using an in vitro model of neuronal differentiation with RFX4 knock-in and a CRISPR-Cas9 knock-out system. As a result, we found that RFX4 directly interacts with the promoters of POU3F2 and NEUROD1. Transcriptomic analysis revealed a set of genes associated with neuronal development, which are highly implicated in the development of neuropsychiatric disorders, including schizophrenia. Notably, ectopic expression of RFX4 can drive human embryonic stem cells toward a neuronal fate. Our results strongly indicate that RFX4 serves as a direct upstream regulator of proneural genes, a role that is essential for normal neuronal development. Impairments in RFX4 function could potentially be related to the development of various neuropsychiatric disorders. However, understanding the precise mechanisms by which the RFX4 gene influences the onset of neuropsychiatric disorders requires further investigation through human genetic studies.

Human RFX6 regulates endoderm patterning at the primitive gut tube stage

Transcriptional factor RFX6 is known to be a causal gene of Mitchell-Riley syndrome (MRS), an autosomal recessive neonatal diabetes associated with pancreatic hypoplasia and intestinal atresia/malformation. The morphological defects are limited to posterior foregut and mid-hindgut endodermal lineages and do not occur in the anterior foregut lineage; the mechanism remains to be fully elucidated. In this study, we generated RFX6+/eGFP heterozygous knockin and RFX6eGFP/eGFP homozygous knockin/knockout human-induced pluripotent stem cell (hiPSC) lines and performed in vitro endoderm differentiation to clarify the role of RFX6 in early endoderm development. RFX6 expression was found to surge at the primitive gut tube (PGT) stage in comparison with that in the undifferentiated or definitive endoderm stage. At the PGT stage, the expression of PDX1 and CDX2, posterior foregut and mid-hindgut master regulators, respectively, was decreased by the RFX6 deficit. PDX1+ and CDX2+ cells were mostly green fluorescent protein (GFP)+ in RFX6+/eGFP hiPSCs, but their cell number was markedly decreased in RFX6eGFP/eGFP hiPSCs. The expression of SOX2, an anterior foregut marker, was not affected by the RFX6 deficit. In addition, we found a putative RFX6-binding X-box motif using cap analysis of gene expression-seq and the motif-containing sequences in the enhancer regions of PDX1 and CDX2 bound to RFX6 in vitro. Thus, RFX6 regulates the ParaHox genes PDX1 and CDX2 but does not affect SOX2 in early endodermal differentiation, suggesting that defects in early stage endoderm patterning account for the morphological pathology of MRS.

Vitamin D and K Supplementation Is Associated with Changes in the Methylation Profile of U266-Multiple Myeloma Cells, Influencing the Proliferative Potential and Resistance to Bortezomib

Multiple myeloma (MM) is a plasma cell malignancy that, despite recent advances in therapy, continues to pose a major challenge to hematologists. Currently, different classes of drugs are applied to treat MM, among others, proteasome inhibitors, immunomodulatory drugs, and monoclonal antibodies. Most of them participate in an interplay with the immune system, hijacking its effector functions and redirecting them to anti-MM activity. Therefore, adjuvant therapies boosting the immune system may be potentially beneficial in MM therapy. Vitamin D (VD) and vitamin K (VK) have multiple so called "non-classical" actions. They exhibit various anti-inflammatory and anti-cancer properties. In this paper, we investigated the influence of VD and VK on epigenetic alterations associated with the proliferative potential of MM cells and the development of BTZ resistance. Our results showed that the development of BTZ resistance is associated with a global decrease in DNA methylation. On the contrary, both control MM cells and BTZ-resistant MM cells exposed to VD alone and to the combination of VD and VK exhibit a global increase in methylation. In conclusion, VD and VK in vitro have the potential to induce epigenetic changes that reduce the proliferative potential of plasma cells and may at least partially prevent the development of resistance to BTZ. However, further ex vivo and in vivo studies are needed to confirm the results and introduce new supplementation recommendations as part of adjuvant therapy.

RFX2 promotes tumor cell stemness through epigenetic regulation of PAF1 in spinal ependymoma

PURPOSE

Spinal ependymoma (SE) is a rare tumor that is most commonly low-grade and tends to recur when complete tumor resection is not feasible. We investigated the molecular mechanism induces stem cell features in SE.

METHODS

Immunohistochemical staining was conducted to analyze the expression of RFX2 in tumor tissues of SE patients at different stages. The expression of tumor stemness markers (Netsin and CD133) was analyzed using western blot analysis and IF, and the efficiency of sphere formation in SE cells was analyzed. The biological activities of SE cells were analyzed by EdU proliferation assay, TUNEL, wound healing, and Transwell assays. The regulatory relationship of RFX2 on PAF1 was verified by ChIP-qPCR and the dual-luciferase assay. SE cells were injected into the spinal cord of nude mice for in vivo assays.

RESULTS

RFX2 was higher in the tumor tissues of SE-III patients than in the tumor tissues of SE-I patients. RFX2 knockdown reduced the expression of tumor stemness markers in SE cells and inhibited the sphere formation efficiency. Moreover, RFX2 knockdown ameliorated the malignant progression of SE in nude mice, as manifested by prolonged survival and alleviated SE tumor infiltration. RFX2 bound to the PAF1 promoter to induce its transcription. Overexpression of PAF1 overturned the effects of RFX2 knockdown on stem cell features and biological activities of SE cells, thereby reducing survival in mice.

CONCLUSIONS

RFX2 activates PAF1 transcription, which promotes tumor stemness of SE cells and leads to the malignant progression of SE.

Characterization of human transcription factor function and patterns of gene regulation in HepG2 cells

Transcription factors (TFs) are trans-acting proteins that bind cis-regulatory elements (CREs) in DNA to control gene expression. Here, we analyzed the genomic localization profiles of 529 sequence-specific TFs and 151 cofactors and chromatin regulators in the human cancer cell line HepG2, for a total of 680 broadly termed DNA-associated proteins (DAPs). We used this deep collection to model each TF's impact on gene expression, and identified a cohort of 26 candidate transcriptional repressors. We examine high occupancy target (HOT) sites in the context of three-dimensional genome organization and show biased motif placement in distal-promoter connections involving HOT sites. We also found a substantial number of closed chromatin regions with multiple DAPs bound, and explored their properties, finding that a MAFF/MAFK TF pair correlates with transcriptional repression. Altogether, these analyses provide novel insights into the regulatory logic of the human cell line HepG2 genome and show the usefulness of large genomic analyses for elucidation of individual TF functions.

Characterizing the SREB G protein-coupled receptor family in fish: Brain gene expression and genomic differences in upstream transcription factor binding sites

The SREB (Super-conserved Receptors Expressed in Brain) family of orphan G protein-coupled receptors is highly conserved in vertebrates and consists of three members: SREB1 (orphan designation GPR27), SREB2 (GPR85), and SREB3 (GPR173). SREBs are associated with processes ranging from neuronal plasticity to reproductive control. Relatively little is known about similarities across the entire family, or how mammalian gene expression patterns compare to non-mammalian vertebrates. In fish, this system may be particularly complex, as some species have gained a fourth member (SREB3B) while others have lost genes. To better understand the system, the present study aimed to: 1) use qPCR to characterize sreb and related gene expression patterns in the brains of three fish species with different systems, and 2) identify possible differences in transcriptional regulation among the receptors, using upstream transcription factor binding sites across 70 ray-finned fish genomes. Overall, regional patterns of sreb expression were abundant in forebrain-related areas. However, some species-specific patterns were detected, such as abundant expression of receptors in zebrafish (Danio rerio) hypothalamic-containing sections, and divergence between sreb3a and sreb3b in pufferfish (Dichotomyctere nigroviridis). In addition, a gene possibly related to the system (dkk3a) was spatially correlated with the receptors in all three species. Genomic regions upstream of sreb2 and sreb3b, but largely not sreb1 or sreb3a, contained many highly conserved transcription factor binding sites. These results provide novel information about expression differences and transcriptional regulation across fish that may inform future research to better understand these receptors.

Synaptic Activity Causes Minute-scale Changes in BAF Complex Composition and Function

Genes encoding subunits of the SWI/SNF or BAF ATP-dependent chromatin remodeling complex are among the most enriched for deleterious de novo mutations in intellectual disabilities and autism spectrum disorder, but the causative molecular pathways are not fully known 1,2 . Synaptic activity in neurons is critical for learning and memory and proper neural development 3 . Neural activity prompts calcium influx and transcription within minutes, facilitated in the nucleus by various transcription factors (TFs) and chromatin modifiers 4 . While BAF is required for activity-dependent developmental processes such as dendritic outgrowth 5-7 , the immediate molecular consequences of neural activity on BAF complexes and their functions are unknown. Here we mapped minute-scale biochemical consequences of neural activity, modeled by membrane depolarization of embryonic mouse primary cortical neurons, on BAF complexes. We used acute chemical perturbations of BAF ATPase activity and kinase signaling to define the activity-dependent effects on BAF complexes and activity-dependent BAF functions. Our studies found that BAF complexes change in subunit composition and are selectively phosphorylated within 10 minutes of depolarization. Increased levels of the core PBAF subunit Baf200/ Arid2 , uniquely containing an RFX-like DNA-binding domain, are concurrent with ATPase-dependent opening of chromatin at RFX/X-box motifs. Changes in BAF composition and phosphorylation lead to the regulation of chromatin accessibility for critical neurogenesis TFs. These biochemical effects are a convergent phenomenon downstream of multiple growth factor signaling pathways in mouse neurons and fibroblasts suggesting that BAF integrates signaling information from the membrane. In support of such a membrane-to-nucleus signaling cascade, we also identified a BAF-interacting kinase, Dclk2, whose inhibition attenuates BAF phosphorylation selectively. Our findings support a direct role of BAF complexes in responding to synaptic activity to regulate TF binding and transcription.

An RFX transcription factor regulates ciliogenesis in the closest living relatives of animals

Cilia allowed our protistan ancestors to sense and explore their environment, avoid predation, and capture bacterial prey.1,2,3 Regulated ciliogenesis was likely critical for early animal evolution,2,4,5,6 and in modern animals, deploying cilia in the right cells at the right time is crucial for development and physiology. Two transcription factors, RFX and FoxJ1, coordinate ciliogenesis in animals7,8,9 but are absent from the genomes of many other ciliated eukaryotes, raising the question of how the regulation of ciliogenesis in animals evolved.10,11 By comparing the genomes of animals with those of their closest living relatives, the choanoflagellates, we found that the genome of their last common ancestor encoded at least three RFX paralogs and a FoxJ1 homolog. Disruption of the RFX homolog cRFXa in the model choanoflagellate Salpingoeca rosetta resulted in delayed cell proliferation and aberrant ciliogenesis, marked by the collapse and resorption of nascent cilia. In cRFXa mutants, ciliogenesis genes and foxJ1 were significantly downregulated. Moreover, the promoters of S. rosetta ciliary genes are enriched for DNA motifs matching those bound by the cRFXa protein in vitro. These findings suggest that an ancestral cRFXa homolog coordinated ciliogenesis in the progenitors of animals and choanoflagellates and that the selective deployment of the RFX regulatory module may have been necessary to differentiate ciliated from non-ciliated cell types during early animal evolution.

Forkhead transcription factor FKH-8 cooperates with RFX in the direct regulation of sensory cilia in Caenorhabditis elegans

Cilia, either motile or non-motile (a.k.a primary or sensory), are complex evolutionarily conserved eukaryotic structures composed of hundreds of proteins required for their assembly, structure and function that are collectively known as the ciliome. Ciliome gene mutations underlie a group of pleiotropic genetic diseases known as ciliopathies. Proper cilium function requires the tight coregulation of ciliome gene transcription, which is only fragmentarily understood. RFX transcription factors (TF) have an evolutionarily conserved role in the direct activation of ciliome genes both in motile and non-motile cilia cell-types. In vertebrates, FoxJ1 and FoxN4 Forkhead (FKH) TFs work with RFX in the direct activation of ciliome genes, exclusively in motile cilia cell-types. No additional TFs have been described to act together with RFX in primary cilia cell-types in any organism. Here we describe FKH-8, a FKH TF, as a direct regulator of the sensory ciliome genes in Caenorhabditis elegans. FKH-8 is expressed in all ciliated neurons in C. elegans, binds the regulatory regions of ciliome genes, regulates ciliome gene expression, cilium morphology and a wide range of behaviors mediated by sensory ciliated neurons. FKH-8 and DAF-19 (C. elegans RFX) physically interact and synergistically regulate ciliome gene expression. C. elegans FKH-8 function can be replaced by mouse FOXJ1 and FOXN4 but not by other members of other mouse FKH subfamilies. In conclusion, RFX and FKH TF families act jointly as direct regulators of ciliome genes also in sensory ciliated cell types suggesting that this regulatory logic could be an ancient trait predating functional cilia sub-specialization.

Single-cell DNA Methylome and 3D Multi-omic Atlas of the Adult Mouse Brain

Cytosine DNA methylation is essential in brain development and has been implicated in various neurological disorders. A comprehensive understanding of DNA methylation diversity across the entire brain in the context of the brain's 3D spatial organization is essential for building a complete molecular atlas of brain cell types and understanding their gene regulatory landscapes. To this end, we employed optimized single-nucleus methylome (snmC-seq3) and multi-omic (snm3C-seq1) sequencing technologies to generate 301,626 methylomes and 176,003 chromatin conformation/methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell type taxonomy that contains 4,673 cell groups and 261 cross-modality-annotated subclasses. We identified millions of differentially methylated regions (DMRs) across the genome, representing potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide multiplexed error-robust fluorescence in situ hybridization (MERFISH2) data validated the association of this spatial epigenetic diversity with transcription and allowed the mapping of the DNA methylation and topology information into anatomical structures more precisely than our dissections. Furthermore, multi-scale chromatin conformation diversities occur in important neuronal genes, highly associated with DNA methylation and transcription changes. Brain-wide cell type comparison allowed us to build a regulatory model for each gene, linking transcription factors, DMRs, chromatin contacts, and downstream genes to establish regulatory networks. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a companion whole-brain SMART-seq3 dataset. Our study establishes the first brain-wide, single-cell resolution DNA methylome and 3D multi-omic atlas, providing an unparalleled resource for comprehending the mouse brain's cellular-spatial and regulatory genome diversity.

Loss of RREB1 in pancreatic beta cells reduces cellular insulin content and affects endocrine cell gene expression

AIMS/HYPOTHESIS

Genome-wide studies have uncovered multiple independent signals at the RREB1 locus associated with altered type 2 diabetes risk and related glycaemic traits. However, little is known about the function of the zinc finger transcription factor Ras-responsive element binding protein 1 (RREB1) in glucose homeostasis or how changes in its expression and/or function influence diabetes risk.

METHODS

A zebrafish model lacking rreb1a and rreb1b was used to study the effect of RREB1 loss in vivo. Using transcriptomic and cellular phenotyping of a human beta cell model (EndoC-βH1) and human induced pluripotent stem cell (hiPSC)-derived beta-like cells, we investigated how loss of RREB1 expression and activity affects pancreatic endocrine cell development and function. Ex vivo measurements of human islet function were performed in donor islets from carriers of RREB1 type 2 diabetes risk alleles.

RESULTS

CRISPR/Cas9-mediated loss of rreb1a and rreb1b function in zebrafish supports an in vivo role for the transcription factor in beta cell mass, beta cell insulin expression and glucose levels. Loss of RREB1 also reduced insulin gene expression and cellular insulin content in EndoC-βH1 cells and impaired insulin secretion under prolonged stimulation. Transcriptomic analysis of RREB1 knockdown and knockout EndoC-βH1 cells supports RREB1 as a novel regulator of genes involved in insulin secretion. In vitro differentiation of RREB1KO/KO hiPSCs revealed dysregulation of pro-endocrine cell genes, including RFX family members, suggesting that RREB1 also regulates genes involved in endocrine cell development. Human donor islets from carriers of type 2 diabetes risk alleles in RREB1 have altered glucose-stimulated insulin secretion ex vivo, consistent with a role for RREB1 in regulating islet cell function.

CONCLUSIONS/INTERPRETATION

Together, our results indicate that RREB1 regulates beta cell function by transcriptionally regulating the expression of genes involved in beta cell development and function.

Multi-omics analysis identifies RFX7 targets involved in tumor suppression and neuronal processes

Recurrently mutated in lymphoid neoplasms, the transcription factor RFX7 is emerging as a tumor suppressor. Previous reports suggested that RFX7 may also have a role in neurological and metabolic disorders. We recently reported that RFX7 responds to p53 signaling and cellular stress. Furthermore, we found RFX7 target genes to be dysregulated in numerous cancer types also beyond the hematological system. However, our understanding of RFX7's target gene network and its role in health and disease remains limited. Here, we generated RFX7 knock-out cells and employed a multi-omics approach integrating transcriptome, cistrome, and proteome data to obtain a more comprehensive picture of RFX7 targets. We identify novel target genes linked to RFX7's tumor suppressor function and underscoring its potential role in neurological disorders. Importantly, our data reveal RFX7 as a mechanistic link that enables the activation of these genes in response to p53 signaling.

Human cytomegalovirus pUL97 upregulates SOCS3 expression via transcription factor RFX7 in neural progenitor cells

Congenital human cytomegalovirus (HCMV) infection causes severe damage to the fetal brain, and the underlying mechanisms remain elusive. Cytokine signaling is delicately controlled in the fetal central nervous system to ensure proper development. Here we show that suppressor of cytokine signaling 3 (SOCS3), a negative feedback regulator of the IL-6 cytokine family signaling, was upregulated during HCMV infection in primary neural progenitor cells (NPCs) with a biphasic expression pattern. From viral protein screening, pUL97 emerged as the viral factor responsible for prolonged SOCS3 upregulation. Further, by proteomic analysis of the pUL97-interacting host proteins, regulatory factor X 7 (RFX7) was identified as the transcription factor responsible for the regulation. Depletion of either pUL97 or RFX7 prevented the HCMV-induced SOCS3 upregulation in NPCs. With a promoter-luciferase activity assay, we demonstrated that the pUL97 kinase activity and RFX7 were required for SOCS3 upregulation. Moreover, the RFX7 phosphorylation level was increased by either UL97-expressing or HCMV-infection in NPCs, suggesting that pUL97 induces RFX7 phosphorylation to drive SOCS3 transcription. We further revealed that elevated SOCS3 expression impaired NPC proliferation and migration in vitro and caused NPCs migration defects in vivo. Taken together, these findings uncover a novel regulatory mechanism of sustained SOCS3 expression in HCMV-infected NPCs, which perturbs IL-6 cytokine family signaling, leads to NPCs proliferation and migration defects, and consequently affects fetal brain development.

Do not overwork: cellular communication network factor 3 for life in cartilage

Cellular communication network factor (CCN) 3, which is one of the founding members of the CCN family, displays diverse functions. However, this protein generally represses the proliferation of a variety of cells. Along with skeletal development, CCN3 is produced in cartilaginous anlagen, growth plate cartilage and epiphysial cartilage. Interestingly, CCN3 is drastically induced in the growth plates of mice lacking CCN2, which promotes endochondral ossification. Notably, chondrocytes in these mutant mice with elevated CCN3 production also suffer from impaired glycolysis and energy metabolism, suggesting a critical role of CCN3 in cartilage metabolism. Recently, CCN3 was found to be strongly induced by impaired glycolysis, and in our study, we located an enhancer that mediated CCN3 regulation via starvation. Subsequent investigations specified regulatory factor binding to the X-box 1 (RFX1) as a transcription factor mediating this CCN3 regulation. Impaired glycolysis is a serious problem, resulting in an energy shortage in cartilage without vasculature. CCN3 produced under such starved conditions restricts energy consumption by repressing cell proliferation, leading chondrocytes to quiescence and survival. This CCN3 regulatory system is indicated to play an important role in articular cartilage maintenance, as well as in skeletal development. Furthermore, CCN3 continues to regulate cartilage metabolism even during the aging process, probably utilizing this regulatory system. Altogether, CCN3 seems to prevent "overwork" by chondrocytes to ensure their sustainable life in cartilage by sensing energy metabolism. Similar roles are suspected to exist in relation to systemic metabolism, since CCN3 is found in the bloodstream.

Determinants of renin cell differentiation: a single cell epi-transcriptomics approach

Rationale

Renin cells are essential for survival. They control the morphogenesis of the kidney arterioles, and the composition and volume of our extracellular fluid, arterial blood pressure, tissue perfusion, and oxygen delivery. It is known that renin cells and associated arteriolar cells descend from FoxD1 + progenitor cells, yet renin cells remain challenging to study due in no small part to their rarity within the kidney. As such, the molecular mechanisms underlying the differentiation and maintenance of these cells remain insufficiently understood.

Objective

We sought to comprehensively evaluate the chromatin states and transcription factors (TFs) that drive the differentiation of FoxD1 + progenitor cells into those that compose the kidney vasculature with a focus on renin cells.

Methods and Results

We isolated single nuclei of FoxD1 + progenitor cells and their descendants from FoxD1 ^cre/+ ; R26R-mTmG mice at embryonic day 12 (E12) (n _cells =1234), embryonic day 18 (E18) (n _cells =3696), postnatal day 5 (P5) (n _cells =1986), and postnatal day 30 (P30) (n _cells =1196). Using integrated scRNA-seq and scATAC-seq we established the developmental trajectory that leads to the mosaic of cells that compose the kidney arterioles, and specifically identified the factors that determine the elusive, myo-endocrine adult renin-secreting juxtaglomerular (JG) cell. We confirm the role of Nfix in JG cell development and renin expression, and identified the myocyte enhancer factor-2 (MEF2) family of TFs as putative drivers of JG cell differentiation.

Conclusions

We provide the first developmental trajectory of renin cell differentiation as they become JG cells in a single-cell atlas of kidney vascular open chromatin and highlighted novel factors important for their stage-specific differentiation. This improved understanding of the regulatory landscape of renin expressing JG cells is necessary to better learn the control and function of this rare cell population as overactivation or aberrant activity of the RAS is a key factor in cardiovascular and kidney pathologies.

Phenotype expansion and neurological manifestations of neurobehavioural disease caused by a variant in RFX7

The RFX7 gene is one of eight genes within the regulatory factor X family. RFX7 is highly expressed in the brain and plays an important role in cell maturation and differentiation. It has only recently been implicated in disease in humans. Reports from 15 individuals have described RFX-associated phenotype as a neurobehavioural disease, manifesting primarily with global developmental delay and intellectual disability. Autism spectrum disorder and attention deficit hyperactivity disorder have also been described in some children. Here we report a case of a 19-month-old with a de novo missense variant in RFX7 resulting in severe global developmental delay including significant speech delay, microcephaly, dyskinetic movements, and failure to thrive. This is the first association between variants in RFX7 and failure to thrive, expanding the phenotype of this newly described gene. In this report we will also show RFX7 associated progressive central nervous system involvement through serial brain imaging.

Co-segregation analysis and functional trial in vivo of candidate genes for monogenic diabetes

INTRODUCTION

The aim of this study was to perform familial co-segregation analysis and functional trial in vivo during mixed meal tolerance test (MMTT) of novel variants in diabetes candidate genes.

RESEARCH DESIGN AND METHODS

It is a continuation of the project "Genetic diabetes in Lithuania" with the cohort of 1209 patients with diabetes. Prior screening for autoimmune markers confirmed type 1 diabetes (T1D) diagnosis in 88.1% (n=1065) of patients, and targeted next-generation sequencing identified 3.5% (n=42) pathogenic variants in MODY genes. Subsequently, 102 patients were classified as having diabetes of unknown etiology. 12/102 were found to have novel variants in potential diabetes genes (RFX2, RREB1, SLC5A1 (3 patients with variants in this gene), GCKR, MC4R, CASP10, TMPRSS6, HGFAC, DACH1, ZBED3). Co-segregation analysis and MMTT were carried out in order to study beta-cell function in subjects with specific variants.

RESULTS

MMTT analysis showed that probands with variants in MC4R, CASP10, TMPRSS6, HGFAC, and SLC5A1 (c.1415T>C) had sufficient residual beta-cell function with stimulated C-peptide (CP) >200 pmol/L. Seven individuals with variants in RFX2, RREB1, GCKR, DACH1, ZBED3 and SLC5A1 (c.1415T>C, and c.932A>T) presented with complete beta-cell failure. No statistical differences were found between patients with sufficient CP production and those with complete beta-cell failure when comparing age at the onset and duration of diabetes. Nineteen family members were included in co-segregation analysis; no diabetes cases were reported among them. Only in patient with the variant c.1894G>A in RFX2 gene, none of the family members were affected by proband's variant.

CONCLUSIONS

Functional beta-cell study in vivo allowed to select five most probable genes for monogenic diabetes. Familial co-segregation analysis showed that novel variant in RFX2 gene could be a possible cause of diabetes. Future functional analysis in vitro is necessary to support or rule out the genetic background as a cause of diabetes.

Deregulated transcription factors and poor clinical outcomes in cancer patients

Transcription factors are a group of proteins, which possess DNA-binding domains, bind to DNA strands of promoters or enhancers, and initiate transcription of genes with cooperation of RNA polymerase and other co-factors. They play crucial roles in regulating transcription during embryogenesis and development. Their physiological status in different cell types is also important to maintain cellular homeostasis. Therefore, any deregulation of transcription factors will lead to the development of cancer cells and tumor progression. Based on their functions in cancer cells, transcription factors could be either oncogenic or tumor suppressive. Furthermore, transcription factors have been shown to modulate cancer stem cells, epithelial-mesenchymal transition (EMT) and drug response; therefore, measuring deregulated transcription factors is hypothesized to predict treatment outcomes of patients with cancers and targeting deregulated transcription factors could be an encouraging strategy for cancer therapy. Here, we summarize the current knowledge of major deregulated transcription factors and their effects on causing poor clinical outcome of patients with cancer. The information presented here will help to predict the prognosis and drug response and to design novel drugs and therapeutic strategies for the treatment of cancers by targeting deregulated transcription factors.

Endotype reversal as a novel strategy for screening drugs targeting familial Alzheimer's disease

While amyloid-β (Aβ) plaques are considered a hallmark of Alzheimer's disease, clinical trials focused on targeting gamma secretase, an enzyme involved in aberrant Aβ peptide production, have not led to amelioration of AD symptoms or synaptic dysregulation. Screening strategies based on mechanistic, multi-omics approaches that go beyond pathological readouts can aid in the evaluation of therapeutics. Using early-onset Alzheimer's (EOFAD) disease patient lineage PSEN1^A246E iPSC-derived neurons, we performed RNA-seq to characterize AD-associated endotypes, which are in turn used as a screening evaluation metric for two gamma secretase drugs, the inhibitor Semagacestat and the modulator BPN-15606. We demonstrate that drug treatment partially restores the neuronal state while concomitantly inhibiting cell cycle re-entry and dedifferentiation endotypes to different degrees depending on the mechanism of gamma secretase engagement. Our endotype-centric screening approach offers a new paradigm by which candidate AD therapeutics can be evaluated for their overall ability to reverse disease endotypes.

The C. elegans regulatory factor X (RFX) DAF-19M module: A shift from general ciliogenesis to cell-specific ciliary and behavioral specialization

Cilia are important for the interaction with environments and the proper function of tissues. While the basic structure of cilia is well conserved, ciliated cells have various functions. To understand the distinctive identities of ciliated cells, the identification of cell-specific proteins and its regulation is essential. Here, we report the mechanism that confers a specific identity on IL2 neurons in Caenorhabditis elegans, neurons important for the dauer larva-specific nictation behavior. We show that DAF-19M, an isoform of the sole C. elegans RFX transcription factor DAF-19, heads a regulatory subroutine, regulating target genes through an X-box motif variant under the control of terminal selector proteins UNC-86 and CFI-1 in IL2 neurons. Considering the conservation of DAF-19M module in IL2 neurons for nictation and in male-specific neurons for mating behavior, we propose the existence of an evolutionarily adaptable, hard-wired genetic module for distinct behaviors that share the feature "recognizing the environment."

The First Report of a Missense Variant in RFX2 Causing Non-Syndromic Tooth Agenesis in a Consanguineous Pakistani Family

Background: The syndromic and non-syndromic congenital missing teeth phenotype is termed tooth agenesis. Since tooth agenesis is a heterogeneous disorder hence, the patients show diverse absent teeth phenotypes. Thus identifying novel genes involved in the morphogenesis of ectodermal appendages, including teeth, paves the way for establishing signaling pathways.

Methods and Results: We have recruited an autosomal recessive non-syndromic tooth agenesis family with two affected members. The exome sequencing technology identified a novel missense sequence variant c.1421T &gt; C; p.(Ile474Thr) in a regulatory factor X (RFX) family member (RFX2, OMIM: 142,765). During the data analysis eight rare variants on various chromosomal locations were identified, but the co-segregation analysis using Sanger sequencing confirmed the segregation of only two variants RFX2: c.1421T &gt; C; p.(Ile474Thr), DOHH: c.109C &gt; G; p.(Pro37Ala) lying in a common 7.1 MB region of homozygosity on chromosome 19p13.3. Furthermore, the online protein prediction algorithms and protein modeling analysis verified the RFX2 variant as a damaging genetic alteration and ACMG pathogenicity criteria classified it as likely pathogenic. On the other hand, the DOHH variant showed benign outcomes.

Conclusion:RFX2 regulates the Hedgehog and fibroblast growth factor signaling pathways, which are involved in the epithelial and mesenchymal interactions during tooth development. Prior animal model studies have confirmed the expression of rfx2 at a developmental stage governing mouth formation. Moreover, its regulatory role and close association with ciliary and non-ciliary genes causing various dental malformations makes it a potential candidate gene for tooth agenesis phenotype. Further studies will contribute to exploring the direct role of RFX2 in human tooth development.

RFX1-mediated CCN3 induction that may support chondrocyte survival under starved conditions

Cellular communication network factor (CCN) family members are multifunctional matricellular proteins that manipulate and integrate extracellular signals. In our previous studies investigating the role of CCN family members in cellular metabolism, we found three members that might be under the regulation of energy metabolism. In this study, we confirmed that CCN2 and CCN3 are the only members that are tightly regulated by glycolysis in human chondrocytic cells. Interestingly, CCN3 was induced under a variety of impaired glycolytic conditions. This CCN3 induction was also observed in two breast cancer cell lines with a distinct phenotype, suggesting a basic role of CCN3 in cellular metabolism. Reporter gene assays indicated a transcriptional regulation mediated by an enhancer in the proximal promoter region. As a result of analyses in silico, we specified regulatory factor binding to the X-box 1 (RFX1) as a candidate that mediated the transcriptional activation by impaired glycolysis. Indeed, the inhibition of glycolysis induced the expression of RFX1, and RFX1 silencing nullified the CCN3 induction by impaired glycolysis. Subsequent experiments with an anti-CCN3 antibody indicated that CCN3 supported the survival of chondrocytes under impaired glycolysis. Consistent with these findings in vitro, abundant CCN3 production by chondrocytes in the deep zones of developing epiphysial cartilage, which are located far away from the synovial fluid, was confirmed in vivo. Our present study uncovered that RFX1 is the mediator that enables CCN3 induction upon cellular starvation, which may eventually assist chondrocytes in retaining their viability, even when there is an energy supply shortage.

Single-nuclei chromatin profiling of ventral midbrain reveals cell identity transcription factors and cell-type-specific gene regulatory variation

BACKGROUND

Cell types in ventral midbrain are involved in diseases with variable genetic susceptibility, such as Parkinson's disease and schizophrenia. Many genetic variants affect regulatory regions and alter gene expression in a cell-type-specific manner depending on the chromatin structure and accessibility.

RESULTS

We report 20,658 single-nuclei chromatin accessibility profiles of ventral midbrain from two genetically and phenotypically distinct mouse strains. We distinguish ten cell types based on chromatin profiles and analysis of accessible regions controlling cell identity genes highlights cell-type-specific key transcription factors. Regulatory variation segregating the mouse strains manifests more on transcriptome than chromatin level. However, cell-type-level data reveals changes not captured at tissue level. To discover the scope and cell-type specificity of cis-acting variation in midbrain gene expression, we identify putative regulatory variants and show them to be enriched at differentially expressed loci. Finally, we find TCF7L2 to mediate trans-acting variation selectively in midbrain neurons.

CONCLUSIONS

Our data set provides an extensive resource to study gene regulation in mesencephalon and provides insights into control of cell identity in the midbrain and identifies cell-type-specific regulatory variation possibly underlying phenotypic and behavioural differences between mouse strains.

Transcription factor RFX7 governs a tumor suppressor network in response to p53 and stress

Despite its prominence, the mechanisms through which the tumor suppressor p53 regulates most genes remain unclear. Recently, the regulatory factor X 7 (RFX7) emerged as a suppressor of lymphoid neoplasms, but its regulation and target genes mediating tumor suppression remain unknown. Here, we identify a novel p53-RFX7 signaling axis. Integrative analysis of the RFX7 DNA binding landscape and the RFX7-regulated transcriptome in three distinct cell systems reveals that RFX7 directly controls multiple established tumor suppressors, including PDCD4, PIK3IP1, MXD4, and PNRC1, across cell types and is the missing link for their activation in response to p53 and stress. RFX7 target gene expression correlates with cell differentiation and better prognosis in numerous cancer types. Interestingly, we find that RFX7 sensitizes cells to Doxorubicin by promoting apoptosis. Together, our work establishes RFX7’s role as a ubiquitous regulator of cell growth and fate determination and a key node in the p53 transcriptional program.

Genome-Wide Association Studies Reveal Neurological Genes for Dog Herding, Predation, Temperament, and Trainability Traits

Genome-wide association study (GWAS) using dog breed standard values as phenotypic measurements is an efficient way to identify genes associated with morphological and behavioral traits. As a result of strong human purposeful selections, several specialized behavioral traits such as herding and hunting have been formed in different modern dog breeds. However, genetic analyses on this topic are rather limited due to the accurate phenotyping difficulty for these complex behavioral traits. Here, 268 dog whole-genome sequences from 130 modern breeds were used to investigate candidate genes underlying dog herding, predation, temperament, and trainability by GWAS. Behavioral phenotypes were obtained from the American Kennel Club based on dog breed standard descriptions or groups (conventional categorization of dog historical roles). The GWAS results of herding behavior (without body size as a covariate) revealed 44 significantly associated sites within five chromosomes. Significantly associated sites on CFA7, 9, 10, and 20 were located either in or near neuropathological or neuronal genes including THOC1, ASIC2, MSRB3, LLPH, RFX8, and CHL1. MSRB3 and CHL1 genes were reported to be associated with dog fear. Since herding is a restricted hunting behavior by removing killing instinct, 36 hounds and 55 herding dogs were used to analyze predation behavior. Three neuronal-related genes (JAK2, MEIS1, and LRRTM4) were revealed as candidates for predation behavior. The significantly associated variant of temperament GWAS was located within ACSS3 gene. The highest associated variant in trainability GWAS is located on CFA22, with no variants detected above the Bonferroni threshold. Since dog behaviors are correlated with body size, we next incorporate body mass as covariates into GWAS; and significant signals around THOC1, MSRB3, LLPH, RFX8, CHL1, LRRTM4, and ACSS3 genes were still detected for dog herding, predation, and temperament behaviors. In humans, these candidate genes are either involved in nervous system development or associated with mental disorders. In conclusion, our results imply that these neuronal or psychiatric genes might be involved in biological processes underlying dog herding, predation, and temperament behavioral traits.

A Chinese hamster transcription start site atlas that enables targeted editing of CHO cells

Chinese hamster ovary (CHO) cells are widely used for producing biopharmaceuticals, and engineering gene expression in CHO is key to improving drug quality and affordability. However, engineering gene expression or activating silent genes requires accurate annotation of the underlying regulatory elements and transcription start sites (TSSs). Unfortunately, most TSSs in the published Chinese hamster genome sequence were computationally predicted and are frequently inaccurate. Here, we use nascent transcription start site sequencing methods to revise TSS annotations for 15 308 Chinese hamster genes and 3034 non-coding RNAs based on experimental data from CHO-K1 cells and 10 hamster tissues. We further capture tens of thousands of putative transcribed enhancer regions with this method. Our revised TSSs improves upon the RefSeq annotation by revealing core sequence features of gene regulation such as the TATA box and the Initiator and, as exemplified by targeting the glycosyltransferase gene Mgat3, facilitate activating silent genes by CRISPRa. Together, we envision our revised annotation and data will provide a rich resource for the CHO community, improve genome engineering efforts and aid comparative and evolutionary studies.

A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ∼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (∼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.

Y disruption, autosomal hypomethylation and poor male lung cancer survival

Lung cancer is the most frequent cause of cancer death worldwide. It affects more men than women, and men generally have worse survival outcomes. We compared gene co-expression networks in affected and unaffected lung tissue from 126 consecutive patients with Stage IA-IV lung cancer undergoing surgery with curative intent. We observed marked degradation of a sex-associated transcription network in tumour tissue. This disturbance, detected in 27.7% of male tumours in the discovery dataset and 27.3% of male tumours in a further 123-sample replication dataset, was coincident with partial losses of the Y chromosome and extensive autosomal DNA hypomethylation. Central to this network was the epigenetic modifier and regulator of sexually dimorphic gene expression, KDM5D. After accounting for prognostic and epidemiological covariates including stage and histology, male patients with tumour KDM5D deficiency showed a significantly increased risk of death (Hazard Ratio [HR] 3.80, 95% CI 1.40-10.3, P = 0.009). KDM5D deficiency was confirmed as a negative prognostic indicator in a further 1100 male lung tumours (HR 1.67, 95% CI 1.4-2.0, P = 1.2 × 10^-10). Our findings identify tumour deficiency of KDM5D as a prognostic marker and credible mechanism underlying sex disparity in lung cancer.

RFX1: a promising therapeutic arsenal against cancer

Regulatory factor X1 (RFX1) is an evolutionary conserved transcriptional factor that influences a wide range of cellular processes such as cell cycle, cell proliferation, differentiation, and apoptosis, by regulating a number of target genes that are involved in such processes. On a closer look, these target genes also play a key role in tumorigenesis and associated events. Such observations paved the way for further studies evaluating the role of RFX1 in cancer. These studies were indispensable due to the failure of conventional chemotherapeutic drugs to target key cellular hallmarks such as cancer stemness, cellular plasticity, enhanced drug efflux, de-regulated DNA repair machinery, and altered pathways evading apoptosis. In this review, we compile significant evidence for the tumor-suppressive activities of RFX1 while also analyzing its oncogenic potential in some cancers. RFX1 induction decreased cellular proliferation, modulated the immune system, induced apoptosis, reduced chemoresistance, and sensitized cancer stem cells for chemotherapy. Thus, our review discusses the pleiotropic function of RFX1 in multitudinous gene regulations, decisive protein–protein interactions, and also its role in regulating key cell signaling events in cancer. Elucidation of these regulatory mechanisms can be further utilized for RFX1 targeted therapy.

Disruption of RFX family transcription factors causes autism, attention-deficit/hyperactivity disorder, intellectual disability, and dysregulated behavior

PURPOSE

We describe a novel neurobehavioral phenotype of autism spectrum disorder (ASD), intellectual disability, and/or attention-deficit/hyperactivity disorder (ADHD) associated with de novo or inherited deleterious variants in members of the RFX family of genes. RFX genes are evolutionarily conserved transcription factors that act as master regulators of central nervous system development and ciliogenesis.

METHODS

We assembled a cohort of 38 individuals (from 33 unrelated families) with de novo variants in RFX3, RFX4, and RFX7. We describe their common clinical phenotypes and present bioinformatic analyses of expression patterns and downstream targets of these genes as they relate to other neurodevelopmental risk genes.

RESULTS

These individuals share neurobehavioral features including ASD, intellectual disability, and/or ADHD; other frequent features include hypersensitivity to sensory stimuli and sleep problems. RFX3, RFX4, and RFX7 are strongly expressed in developing and adult human brain, and X-box binding motifs as well as RFX ChIP-seq peaks are enriched in the cis-regulatory regions of known ASD risk genes.

CONCLUSION

These results establish a likely role of deleterious variation in RFX3, RFX4, and RFX7 in cases of monogenic intellectual disability, ADHD and ASD, and position these genes as potentially critical transcriptional regulators of neurobiological pathways associated with neurodevelopmental disease pathogenesis.

Long non-coding RNA PCA3 inhibits lipid accumulation and atherosclerosis through the miR-140-5p/RFX7/ABCA1 axis

OBJECTIVE

The purpose of this study was to explore the role of long noncoding RNA (lncRNA) prostate cancer antigen 3 (PCA3) in atherosclerosis and the underlying mechanism.

METHODS

The Gene Expression Omnibus (GEO) datasets were used to divide differentially expressed lncRNAs, microRNAs (miRNAs), and mRNAs. The expression of PCA3, miR-140-5p, RFX7 and ABCA1 were determined by qPCR or Western blot in ox-LDL-treated macrophages. Macrophage lipid accumulation s was evaluated using the Oil Red O staining and high-performance liquid chromatography. Target relationships among PCA3, miR-140-5p, RFX7, and ABCA1 promoter area were validated via dual-luciferase reporter gene assay or chromatin immunoprecipitation assay. The apoE^-/- mouse model in vivo was designed to evaluate the effect of PCA3 on the reverse cholesterol transport (RCT) and atherosclerosis.

RESULTS

PCA3 was down-regulated in foam cells, whereas miR-140-5p was highly expressed. Overexpression of PCA3 promoted ABCA1-mediated cholesterol efflux and reduced lipid accumulation in macrophages. Besides, RFX7 bound to the ABCA1 promoter and increased ABCA1 expression. Targeted relationships and interactions on the expression between miR-140-5p and PCA3 or RFX7 were elucidated. PCA3 up-regulated ABCA1 expression by binding to miR-140-5p to up-regulate RFX7 and ABCA1 expression in macrophages. PCA3 promoted RCT and impeded the progression of atherosclerosis by sponging miR-140-5p in apoE^-/- mice. Meanwhile, miR-140-5p also inhibit ABCA1 expression via downregulation of RFX7 to impede RCT and aggravate atherosclerosis.

CONCLUSIONS

lncRNA PCA3 promotes ABCA1-mediated cholesterol efflux to inhibit atherosclerosis through sponging miR-140-5p and up-regulating RFX7.

Differentiation of ciliated human midbrain-derived LUHMES neurons

Many human cell types are ciliated, including neural progenitors and differentiated neurons. Ciliopathies are characterized by defective cilia and comprise various disease states, including brain phenotypes, where the underlying biological pathways are largely unknown. Our understanding of neuronal cilia is rudimentary, and an easy-to-maintain, ciliated human neuronal cell model is absent. The Lund human mesencephalic (LUHMES) cell line is a ciliated neuronal cell line derived from human fetal mesencephalon. LUHMES cells can easily be maintained and differentiated into mature, functional neurons within one week. They have a single primary cilium as proliferating progenitor cells and as postmitotic, differentiating neurons. These developmental stages are completely separable within one day of culture condition change. The sonic hedgehog (SHH) signaling pathway is active in differentiating LUHMES neurons. RNA-sequencing timecourse analyses reveal molecular pathways and gene-regulatory networks critical for ciliogenesis and axon outgrowth at the interface between progenitor cell proliferation, polarization and neuronal differentiation. Gene expression dynamics of cultured LUHMES neurons faithfully mimic the corresponding in vivo dynamics of human fetal midbrain. In LUHMES cells, neuronal cilia biology can be investigated from proliferation through differentiation to mature neurons.

Interplay of RFX transcription factors 1, 2 and 3 in motile ciliogenesis

Cilia assembly is under strict transcriptional control during animal development. In vertebrates, a hierarchy of transcription factors (TFs) are involved in controlling the specification, differentiation and function of multiciliated epithelia. RFX TFs play key functions in the control of ciliogenesis in animals. Whereas only one RFX factor regulates ciliogenesis in C. elegans, several distinct RFX factors have been implicated in this process in vertebrates. However, a clear understanding of the specific and redundant functions of different RFX factors in ciliated cells remains lacking. Using RNA-seq and ChIP-seq approaches we identified genes regulated directly and indirectly by RFX1, RFX2 and RFX3 in mouse ependymal cells. We show that these three TFs have both redundant and specific functions in ependymal cells. Whereas RFX1, RFX2 and RFX3 occupy many shared genomic loci, only RFX2 and RFX3 play a prominent and redundant function in the control of motile ciliogenesis in mice. Our results provide a valuable list of candidate ciliary genes. They also reveal stunning differences between compensatory processes operating in vivo and ex vivo.

Complex Genetics in Pancreatitis: Insights Gained From a New Candidate Locus Panel

OBJECTIVES

Chronic pancreatitis is the end stage of a pathologic inflammatory syndrome with multiple etiological factors, including genetic. We hypothesized that some pancreatitis etiology originates in pancreatic acinar or duct cells and requires both injury and compensatory mechanism failure.

METHODS

One hundred pancreatitis patients were assessed using a DNA sequencing panel for pancreatitis. Cooccurrence of variants within and between genes was measured. Gene coexpression was confirmed via published single-cell RNA sequencing.

RESULTS

One hundred and twenty-one variants were identified in 2 or more patients, 15 of which were enriched compared with reference populations. Single cell RNA-sequencing data verified coexpression of GGT1, CFTR, and PRSS1 in duct cells, PRSS1, CPA1, CEL, CTRC, and SPINK1 in acinar cells, and UBR1 in both. Multiple-risk variants with injury/stress effects (CEL, CFTR, CPA1, PRSS1) and impaired cell protection (CTRC, GGT1, SPINK1, UBR1) cooccur within duct cells, acinar cells, or both.

CONCLUSIONS

Pancreatitis is a complex disorder with genetic interactions across genes and cell types. These findings suggest a new, non-Mendelian genetic risk/etiology paradigm where a combination of nonpathogenic genetic risk variants in groups of susceptibility genes and injury/dysfunction response genes contribute to acquired pancreatic disease.