Search

Background: X-linked acrogigantism (X-LAG; MIM: 300942) is a severe form of pituitary gigantism caused by chromosome Xq26.3 duplications involving GPR101. X-LAG-associated duplications disrupt the integrity of the topologically associating domain (TAD) containing GPR101 and lead to the formation of a neo-TAD that drives pituitary GPR101 misexpression and gigantism. As X-LAG is fully penetrant and heritable, duplications involving GPR101 identified on prenatal screening studies, like amniocentesis, can pose an interpretation challenge for medical geneticists and raise important concerns for patients and families. Therefore, providing robust information on the functional genomic impact of such duplications has important research and clinical value with respect to gene regulation and triplosensitivity traits. Methods: We employed 4C/HiC-seq as a clinical tool to determine the functional impact of incidentally discovered GPR101 duplications on TAD integrity in three families. After defining duplications and breakpoints around GPR101 by clinical-grade and high-density aCGH, we constructed 4C/HiC chromatin contact maps for our study population and compared them with normal and active (X-LAG) controls. Results: We showed that duplications involving GPR101 that preserved the centromeric invariant TAD boundary did not generate a pathogenic neo-TAD and that ectopic enhancers were not adopted. This allowed us to discount presumptive/suspected X-LAG diagnoses and GPR101 misexpression, obviating the need for intensive clinical follow-up. Conclusions: This study highlights the importance of TAD boundaries and chromatin interactions in determining the functional impact of copy number variants and provides proof-of-concept for using 4C/HiC-seq as a clinical tool to acquire crucial information for genetic counseling and to support clinical decision-making in cases of suspected TADopathies. [ABSTRACT FROM AUTHOR]

Histone modifications, known as histone marks, are pivotal in regulating gene expression within cells. The vast array of potential combinations of histone marks presents a considerable challenge in decoding the regulatory mechanisms solely through biological experimental approaches. To overcome this challenge, we have developed a method called CatLearning. It utilizes a modified convolutional neural network architecture with a specialized adaptation Residual Network to quantitatively interpret histone marks and predict gene expression. This architecture integrates long-range histone information up to 500Kb and learns chromatin interaction features without 3D information. By using only one histone mark, CatLearning achieves a high level of accuracy. Furthermore, CatLearning predicts gene expression by simulating changes in histone modifications at enhancers and throughout the genome. These findings help comprehend the architecture of histone marks and develop diagnostic and therapeutic targets for diseases with epigenetic changes. [ABSTRACT FROM AUTHOR]

The myeloma overexpressed gene (MYEOV) has been proposed to be a protooncogene due to high RNA transcript levels found in multiple cancers, including myeloma, breast, lung, pancreas and esophageal cancer. The presence of an open reading frame (ORF) in humans and other primates suggests protein-coding potential. Yet, we still lack evidence of a functional MYEOV protein. It remains undetermined how MYEOV overexpression affects cancerous tissues. In this work, we show that MYEOV has likely originated and may still function as an enhancer, regulating CCND1 and LTO1. Firstly, MYEOV 3' enhancer activity was confirmed in humans using publicly available ATAC-STARR-seq data, performed on B-cell-derived GM12878 cells. We detected enhancer histone marks H3K4me1 and H3K27ac overlapping MYEOV in multiple healthy human tissues, which include B cells, liver and lung tissue. The analysis of 3D genome datasets revealed chromatin interactions between a MYEOV-3'-putative enhancer and the proto-oncogene CCND1. BLAST searches and multi-sequence alignment results showed that DNA sequence from this human enhancer element is conserved from the amphibians/amniotes divergence, with a 273 bp conserved region also found in all mammals, and even in chickens, where it is consistently located near the corresponding CCND1 orthologues. Furthermore, we observed conservation of an active enhancer state in the MYEOV orthologues of four non-human primates, dogs, rats, and mice. When studying this homologous region in mice, where the ORF of MYEOV is absent, we not only observed an enhancer chromatin state but also found interactions between the mouse enhancer homolog and Ccnd1 using 3D-genome interaction data. This is similar to the interaction observed in humans and, interestingly, coincides with CTCF binding sites in both species. Taken together, this suggests that MYEOV is a primatespecific gene with a de novo ORF that originated at an evolutionarily older enhancer region. This deeply conserved putative enhancer element could regulate CCND1 in both humans and mice, opening the possibility of studying MYEOV regulatory functions in cancer using non-primate animal models. [ABSTRACT FROM AUTHOR]

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

MDS1 and EVI1 complex locus (MECOM), a transcription factor encoding several variants, has been implicated in progression of ovarian cancer. The function of regulatory regions in regulating MECOM expression in ovarian cancer is not fully understood. In this study, MECOM expression was evaluated in ovarian cancer cell lines treated with bromodomain and extraterminal (BET) inhibitor JQ-1. Oncogenic phenotypes were assayed using assays of CCK-8, colony formation, wound-healing and transwell. Oncogenic phenotypes were estimated in stable sgRNA-transfected OVCAR3 cell lines. Xenograft mouse model was assayed via subcutaneous injection of enhancer-deleted OVCAR3 cell lines. The results displayed that expression of MECOM is downregulated in cell lines treated with JQ-1. Data from published ChIP-sequencing (H3K27Ac) in 3 ovarian cancer cell lines displayed a potential enhancer around the first exon. mRNA and protein expression were downregulated in OVCAR3 cells after deletion of the MECOM enhancer. Similarly, oncogenic phenotypes both in cells and in the xenograft mouse model were significantly attenuated. This study demonstrates that JQ-1 can inhibit the expression of MECOM and tumorigenesis. Deletion of the enhancer activity of MECOM has an indispensable role in inhibiting ovarian cancer progress, which sheds light on a promising opportunity for ovarian cancer treatment through the application of this non-coding DNA deletion. [ABSTRACT FROM AUTHOR]

Drosophila Hey is a basic helix–loop–helix–orange (bHLH-O) protein with an important role in the establishment of distinct identities of postmitotic cells. We have previously identified Hey as a transcriptional target and effector of Notch signalling during the asymmetric division of neuronal progenitors, generating neurons of two types, and we have shown that Notch-dependent expression of Hey also marks a subpopulation of the newborn enteroendocrine (EE) cells in the midgut primordium of the embryo. Here, we investigate the transcriptional regulation of Hey in neuronal and intestinal tissues. We isolated two genomic regions upstream of the promoter (HeyUP) and in the second intron (HeyIN2) of the Hey gene, based on the presence of binding motifs for Su(H), the transcription factor that mediates Notch activity. We found that both regions can direct the overlapping expression patterns of reporter transgenes recapitulating endogenous Hey expression. Moreover, we showed that while HeyIN2 represents a Notch-dependent enhancer, HeyUP confers both Notch-dependent and independent transcriptional regulation. We induced mutations that removed the Su(H) binding motifs in either region and then studied the enhancer functionality in the respective Hey mutant lines. Our results provide direct evidence that although both enhancers support Notch-dependent regulation of the Hey gene, their role is redundant, as a Hey loss-of-function lethal phenotype is observed only after deletion of all their Su(H) binding motifs by CRISPR/Cas9. [ABSTRACT FROM AUTHOR]

Gene targeting (GT) allows precise manipulation of genome sequences, such as knock-ins and sequence substitutions, but GT in seed plants remains a challenging task. Engineered sequence-specific nucleases (SSNs) are known to facilitate GT via homology-directed repair (HDR) in organisms. Here, we demonstrate that Cas12a and a temperature-tolerant Cas12a variant (ttCas12a) can efficiently establish precise and heritable GT at two loci in Arabidopsis thaliana (Arabidopsis) through a sequential transformation strategy. As a result, ttCas12a showed higher GT efficiency than unmodified Cas12a. In addition, the efficiency of transcriptional and translational enhancers for GT via sequential transformation strategy was also investigated. These enhancers and their combinations were expected to show an increase in GT efficiency in the sequential transformation strategy, similar to previous reports of all-in-one strategies, but only a maximum twofold increase was observed. These results indicate that the frequency of double strand breaks (DSBs) at the target site is one of the most important factors determining the efficiency of genetic GT in plants. On the other hand, a higher frequency of DSBs does not always lead to higher efficiency of GT, suggesting that some additional factors are required for GT via HDR. Therefore, the increase in DSB can no longer be expected to improve GT efficiency, and a new strategy needs to be established in the future. This research opens up a wide range of applications for precise and heritable GT technology in plants. Key message: The combination of enhancers, Cas12a, and sequential transformation strategies enables efficient and precise heritable gene targeting in Arabidopsis. [ABSTRACT FROM AUTHOR]

Asthma is a widespread airway disorder where GATA3-dependent Type-2 helper T (Th2) cells and group 2 innate lymphoid cells (ILC2s) play vital roles. Asthma-associated single nucleotide polymorphisms (SNPs) are enriched in a region located 926-970 kb downstream from GATA3 in the 10p14 (hG900). However, it is unknown how hG900 affects the pathogenesis of allergic airway inflammation. To investigate the roles of the asthma-associated GATA3 enhancer region in experimental allergic airway inflammation, we first examined the correlation between GATA3 expression and the activation of the hG900 region was analyzed by flow cytometry and ChIP-qPCR. We found that The activation of enhancers in the hG900 region was strongly correlated to the levels of GATA3 in human peripheral T cell subsets. We next generated mice lacking the mG900 region (mG900KO mice) were generated by the CRISPR-Cas9 system, and the development and function of helper T cells and ILCs in mG900KO mice were analyzed in steady-state conditions and allergic airway inflammation induced by papain or house dust mite (HDM). The deletion of the mG900 did not affect the development of lymphocytes in steady-state conditions or allergic airway inflammation induced by papain. However, mG900KO mice exhibited reduced allergic inflammation and Th2 differentiation in the HDM-induced allergic airway inflammation. The analysis of the chromatin conformation around Gata3 by circular chromosome conformation capture coupled to high-throughput sequencing (4C-seq) revealed that the mG900 region interacted with the transcription start site of Gata3 with an influencing chromatin conformation in Th2 cells. These findings indicate that the mG900 region plays a pivotal role in Th2 differentiation and thus enhances allergic airway inflammation. [ABSTRACT FROM AUTHOR]

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation. [ABSTRACT FROM AUTHOR]

The abnormality of the p53 tumor suppressor is crucial in lung cancer development, because p53 regulates target gene promoters to combat cancer. Recent studies have shown extensive p53 binding to enhancer elements. However, whether p53 exerts a tumor suppressor role by shaping the enhancer landscape remains poorly understood. In the current study, we employed several functional genomics approaches to assess the enhancer activity at p53 binding sites throughout the genome based on our established TP53 knockout (KO) human bronchial epithelial cells (BEAS-2B). A total of 943 active regular enhancers and 370 super-enhancers (SEs) disappeared upon the deletion of p53, indicating that p53 modulates the activity of hundreds of enhancer elements. We found that one p53-dependent SE, located on chromosome 9 and designated as KLF4-SE, regulated the expression of the Krüppel-like factor 4 (KLF4) gene. Furthermore, the deletion of p53 significantly decreased the KLF4-SE enhancer activity and the KLF4 expression, but increased colony formation ability in the nitrosamines 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone-induced cell transformation model. Subsequently, in TP53 KO cells, the overexpression of KLF4 partially reversed the increased clonogenic capacity caused by p53 deficiency. Consistently, KLF4 expression also decreased in lung cancer tissues and cell lines. It appeared that overexpression of KLF4 significantly suppressed the proliferation and migration of lung cancer cells. Collectively, our results suggest that the regulation of enhancer formation and activity by p53 is an integral component of the p53 tumor suppressor function. Therefore, our findings offer some novel insights into the regulation mechanism of p53 in lung oncogenesis and introduce a new strategy for screening therapeutic targets. [ABSTRACT FROM AUTHOR]

Abstract Gene targeting (GT) allows precise manipulation of genome sequences, such as knock-ins and sequence substitutions, but GT in seed plants remains a challenging task. Engineered sequence-specific nucleases (SSNs) are known to facilitate GT via homology-directed repair (HDR) in organisms. Here, we demonstrate that Cas12a and a temperature-tolerant Cas12a variant (ttCas12a) can efficiently establish precise and heritable GT at two loci in Arabidopsis thaliana (Arabidopsis) through a sequential transformation strategy. As a result, ttCas12a showed higher GT efficiency than unmodified Cas12a. In addition, the efficiency of transcriptional and translational enhancers for GT via sequential transformation strategy was also investigated. These enhancers and their combinations were expected to show an increase in GT efficiency in the sequential transformation strategy, similar to previous reports of all-in-one strategies, but only a maximum twofold increase was observed. These results indicate that the frequency of double strand breaks (DSBs) at the target site is one of the most important factors determining the efficiency of genetic GT in plants. On the other hand, a higher frequency of DSBs does not always lead to higher efficiency of GT, suggesting that some additional factors are required for GT via HDR. Therefore, the increase in DSB can no longer be expected to improve GT efficiency, and a new strategy needs to be established in the future. This research opens up a wide range of applications for precise and heritable GT technology in plants.

Abstract Background Recognition of enhancer–promoter Interactions (EPIs) is crucial for human development. EPIs in the genome play a key role in regulating transcription. However, experimental approaches for classifying EPIs are too expensive in terms of effort, time, and resources. Therefore, more and more studies are being done on developing computational techniques, particularly using deep learning and other machine learning techniques, to address such problems. Unfortunately, the majority of current computational methods are based on convolutional neural networks, recurrent neural networks, or a combination of them, which don’t take into consideration contextual details and the long-range interactions between the enhancer and promoter sequences. A new transformer-based model called EPI-Trans is presented in this study to overcome the aforementioned limitations. The multi-head attention mechanism in the transformer model automatically learns features that represent the long interrelationships between enhancer and promoter sequences. Furthermore, a generic model is created with transferability that can be utilized as a pre-trained model for various cell lines. Moreover, the parameters of the generic model are fine-tuned using a particular cell line dataset to improve performance. Results Based on the results obtained from six benchmark cell lines, the average AUROC for the specific, generic, and best models is 94.2%, 95%, and 95.7%, while the average AUPR is 80.5%, 66.1%, and 79.6% respectively. Conclusions This study proposed a transformer-based deep learning model for EPI prediction. The comparative results on certain cell lines show that EPI-Trans outperforms other cutting-edge techniques and can provide superior performance on the challenge of recognizing EPI.

Abstract Background Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma. Methods Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3. Results An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4. Conclusions Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment.

Abstract Background Cervical cancer is a human papillomavirus (HPV)-related disease. HPV type 16 (HPV16), which is the predominant cause of cervical cancer, can encode miRNAs (HPV16-miRNAs). However, the role of HPV16-miRNAs in the pathogenesis of cervical cancer remains unclear. Methods Human cervical cancer cell lines SiHa (HPV16-positive) and C33A (HPV-negative), and cervical cancer tissues were collected to investigate the expression levels of two HPV16-miRNAs (HPV16-miR-H1 and HPV16-miR-H6). The overexpression and knockdown of HPV16-miR-H1 and HPV16-miR-H6 were performed using the lentiviral vector system and miRNA inhibitors, respectively. RNA-sequencing (RNA-seq) analysis and H3K27ac chromatin immunoprecipitation and sequencing (CHIP-seq) experiments were utilized to explore the roles of HPV16-miR-H1 and HPV16-miR-H6 facilitated by enhancers. CCK8, EdU, transwell, and wound healing assays were performed to verify the effects of HPV16-miR-H1 and HPV16-miR-H6 on cell proliferation and migration. Results HPV16-miR-H1 and HPV16-miR-H6 were highly expressed in both SiHa cells and tissue samples from HPV16-positive cervical cancer patients. RNA-seq analysis showed that HPV16-miR-H1 and HPV16-miR-H6 induced the upregulation of numerous tumor progression-associated genes. H3K27ac CHIP-seq experiments further revealed that HPV16-miR-H1 and HPV16-miR-H6 modulated the expression of critical genes by regulating their enhancer activity. The functional study demonstrated that HPV16-miR-H1 and HPV16-miR-H6 increased the migratory capacity of SiHa cells. Conclusions Our data shed light on the role of HPV16-encoded miRNAs in cervical cancer, particularly emphasizing their involvement in the miRNA-enhancer-target gene system. This novel regulatory mechanism of HPV16-miRNAs provides new insights and approaches for the development of therapeutic strategies by targeting HPV16-positive cervical cancer.

Transcriptional dysregulation of genes is a hallmark of tumors and can serve as targets for cancer drug development. However, it is extremely challenging to develop small-molecule inhibitors to target abnormally expressed transcription factors (TFs) except for the nuclear receptor family of TFs. Little is known about the interaction between TFs and transcription cofactors in gastroesophageal adenocarcinoma (GEA) or the therapeutic effects of targeting TF and transcription cofactor complexes. In this study, we found that ETS homologous factor (EHF) expression is promoted by a core transcriptional regulatory circuitry (CRC), specifically ELF3-KLF5-GATA6, and interference with its expression suppressed the malignant biological behavior of GEA cells. Importantly, we identified Ajuba LIM protein (AJUBA) as a new coactivator of EHF that cooperatively orchestrates transcriptional network activity in GEA. Furthermore, we identified KRAS signaling as a common pathway downstream of EHF and AJUBA. Applicably, dual targeting of EHF and AJUBA by lipid nanoparticles cooperatively attenuated the malignant biological behaviors of GEA in vitro and in vivo. In conclusion, EHF is upregulated by the CRC and promotes GEA malignancy by interacting with AJUBA through the KRAS pathway. Targeting of both EHF and its coactivator AJUBA through lipid nanoparticles is a novel potential therapeutic strategy.

The Grey allele in horses is causing premature hair greying and susceptibility to melanoma. The causal mutation is a 4.6 kb tandem duplication in intron 6 of the Syntaxin 17 gene. A recent study demonstrated that the most common allele at the Grey locus (G3) involves three tandem copies of this sequence, whilst a more rare allele (G2) has two tandem copies and the wild-type allele (G1) only one copy. The G3 allele is causing fast greying and high incidence of skin melanoma, whereas the G2 allele is causing slow greying and no obvious increase in melanoma incidence. Further somatic copy number expansion has been documented in melanoma tissue from Grey horses. Functional studies showed that this intronic sequence acts as a weak melanocyte-specific enhancer that becomes substantially stronger by the copy number expansion. The Grey mutation is associated with upregulated expression of both Syntaxin 17 and the neighbouring NR4A3 gene in Grey horse melanomas. It is still an open question which of these genes is most important for the phenotypic effects or if causality is due to the combined effect of upregulation of both genes. Interestingly, RNAseq data in the Human Protein Atlas give support for a possible role of NR4A3 because it is particularly upregulated in human skin cancer, and it belongs to a cluster of genes associated with skin cancer and melanin biosynthesis. The Grey mutation and its association with melanoma provide a possibility to study the path to tumour development in numerous Grey horses carrying exactly the same predisposing mutation.

Background: Recognition of enhancer–promoter Interactions (EPIs) is crucial for human development. EPIs in the genome play a key role in regulating transcription. However, experimental approaches for classifying EPIs are too expensive in terms of effort, time, and resources. Therefore, more and more studies are being done on developing computational techniques, particularly using deep learning and other machine learning techniques, to address such problems. Unfortunately, the majority of current computational methods are based on convolutional neural networks, recurrent neural networks, or a combination of them, which don't take into consideration contextual details and the long-range interactions between the enhancer and promoter sequences. A new transformer-based model called EPI-Trans is presented in this study to overcome the aforementioned limitations. The multi-head attention mechanism in the transformer model automatically learns features that represent the long interrelationships between enhancer and promoter sequences. Furthermore, a generic model is created with transferability that can be utilized as a pre-trained model for various cell lines. Moreover, the parameters of the generic model are fine-tuned using a particular cell line dataset to improve performance. Results: Based on the results obtained from six benchmark cell lines, the average AUROC for the specific, generic, and best models is 94.2%, 95%, and 95.7%, while the average AUPR is 80.5%, 66.1%, and 79.6% respectively. Conclusions: This study proposed a transformer-based deep learning model for EPI prediction. The comparative results on certain cell lines show that EPI-Trans outperforms other cutting-edge techniques and can provide superior performance on the challenge of recognizing EPI. [ABSTRACT FROM AUTHOR]

Background: Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma. Methods: Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3. Results: An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4. Conclusions: Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment. [ABSTRACT FROM AUTHOR]

Kaposi’s sarcoma-associated herpesvirus (KSHV) belongs to the gamma-herpesvirus family and is a well-known human oncogenic virus. In infected cells, the viral genome of 165 kbp is circular DNA wrapped in chromatin. The tight control of gene expression is critical for latency, the transition into the lytic phase, and the development of viral-associated malignancies. Distal cis-regulatory elements, such as enhancers and silencers, can regulate gene expression in a position- and orientation-independent manner. Open chromatin is another characteristic feature of enhancers. To systematically search for enhancers, we cloned all the open chromatin regions in the KSHV genome downstream of the luciferase gene and tested their enhancer activity in infected and uninfected cells. A silencer was detected upstream of the latency-associated nuclear antigen promoter. Two constitutive enhancers were identified in the K12p-OriLyt-R and ORF29 Intron regions, where ORF29 Intron is a tissue-specific enhancer. The following promoters: OriLyt-L, PANp, ALTp, and the terminal repeats (TRs) acted as lytically induced enhancers. The expression of the replication and transcription activator (RTA), the master regulator of the lytic cycle, was sufficient to induce the activity of lytic enhancers in uninfected cells. We propose that the TRs that span about 24 kbp region serve as a “viral super-enhancer” that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The silencer and enhancers described here provide an additional layer to the complex gene regulation of herpesviruses. IMPORTANCE In this study, we performed a systematic functional assay to identify cis-regulatory elements within the genome of the oncogenic herpesvirus, Kaposi’s sarcoma-associated herpesvirus (KSHV). Similar to other herpesviruses, KSHV presents both latent and lytic phases. Therefore, our assays were performed in uninfected cells, during latent infection, and under lytic conditions. We identified two constitutive enhancers, one of which seems to be a tissue-specific enhancer. In addition, four lytically induced enhancers, which are all responsive to the replication and transcription activator (RTA), were identified. Furthermore, a silencer was identified between the major latency promoter and the lytic gene locus. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The terminal repeats, spanning a region of about 24 kbp, seem like a “viral super-enhancer” that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA to regulate latency to lytic transition. [ABSTRACT FROM AUTHOR]

During development, neural stem cells in the cerebral cortex, also known as radial glial cells (RGCs), generate excitatory neurons, followed by production of cortical macroglia and inhibitory neurons that migrate to the olfactory bulb (OB). Understanding the mechanisms for this lineage switch is fundamental for unraveling how proper numbers of diverse neuronal and glial cell types are controlled. We and others recently showed that Sonic Hedgehog (Shh) signaling promotes the cortical RGC lineage switch to generate cortical oligodendrocytes and OB interneurons. During this process, cortical RGCs generate intermediate progenitor cells that express critical gliogenesis genes Ascll, Egfr, and Olig2. The increased Ascll expression and appearance of Egfr+ and Olig2+ cortical progenitors are concurrent with the switch from excitatory neurogenesis to gliogenesis and OB interneuron neurogenesis in the cortex. While Shh signaling promotes Olig2 expression in the developing spinal cord, the exact mechanism for this transcriptional regulation is not known. Furthermore, the transcriptional regulation of Olig2 and Egfr has not been explored. Here, we show that in cortical progenitor cells, multiple regu-latory programs, including Pax6 and Gli3, prevent precocious expression of Olig2, a gene essential for production of cortical oligodendrocytes and astrocytes. We identify multiple enhancers that control Olig2 expression in cortical progenitors and show that the mechanisms for regulating Olig2 expression are conserved between the mouse and human. Our study reveals evolutionarily conserved regulatory logic controlling the lineage switch of cortical neural stem cells. [ABSTRACT FROM AUTHOR]

Background: Cervical cancer is a human papillomavirus (HPV)-related disease. HPV type 16 (HPV16), which is the predominant cause of cervical cancer, can encode miRNAs (HPV16-miRNAs). However, the role of HPV16-miRNAs in the pathogenesis of cervical cancer remains unclear. Methods: Human cervical cancer cell lines SiHa (HPV16-positive) and C33A (HPV-negative), and cervical cancer tissues were collected to investigate the expression levels of two HPV16-miRNAs (HPV16-miR-H1 and HPV16-miR-H6). The overexpression and knockdown of HPV16-miR-H1 and HPV16-miR-H6 were performed using the lentiviral vector system and miRNA inhibitors, respectively. RNA-sequencing (RNA-seq) analysis and H3K27ac chromatin immunoprecipitation and sequencing (CHIP-seq) experiments were utilized to explore the roles of HPV16-miR-H1 and HPV16-miR-H6 facilitated by enhancers. CCK8, EdU, transwell, and wound healing assays were performed to verify the effects of HPV16-miR-H1 and HPV16-miR-H6 on cell proliferation and migration. Results: HPV16-miR-H1 and HPV16-miR-H6 were highly expressed in both SiHa cells and tissue samples from HPV16-positive cervical cancer patients. RNA-seq analysis showed that HPV16-miR-H1 and HPV16-miR-H6 induced the upregulation of numerous tumor progression-associated genes. H3K27ac CHIP-seq experiments further revealed that HPV16-miR-H1 and HPV16-miR-H6 modulated the expression of critical genes by regulating their enhancer activity. The functional study demonstrated that HPV16-miR-H1 and HPV16-miR-H6 increased the migratory capacity of SiHa cells. Conclusions: Our data shed light on the role of HPV16-encoded miRNAs in cervical cancer, particularly emphasizing their involvement in the miRNA-enhancer-target gene system. This novel regulatory mechanism of HPV16-miRNAs provides new insights and approaches for the development of therapeutic strategies by targeting HPV16-positive cervical cancer. [ABSTRACT FROM AUTHOR]

High-resolution melting (HRM) is a cost-efficient tool for targeted DNA methylation analysis. HRM yields the average methylation status across all CpGs in PCR products. Moreover, it provides information on the methylation pattern, e.g., the occurrence of monoallelic methylation. HRM assays have to be calibrated by analyzing DNA methylation standards of known methylation status and mixtures thereof. In general, DNA methylation levels determined by the classical calibration approach, including the whole temperature range in between normalization intervals, are in good agreement with the mean of the DNA methylation status of individual CpGs determined by pyrosequencing (PSQ), the gold standard of targeted DNA methylation analysis. However, the classical calibration approach leads to highly inaccurate results for samples with heterogeneous DNA methylation since they result in more complex melt curves, differing in their shape compared to those of DNA standards and mixtures thereof. Here, we present a novel calibration approach, i.e., temperature-wise calibration. By temperature-wise calibration, methylation profiles over temperature are obtained, which help in finding the optimal calibration range and thus increase the accuracy of HRM data, particularly for heterogeneous DNA methylation. For explaining the principle and demonstrating the potential of the novel calibration approach, we selected the promoter and two enhancers of MGMT, a gene encoding the repair protein MGMT. [ABSTRACT FROM AUTHOR]

Transcriptional dysregulation of genes is a hallmark of tumors and can serve as targets for cancer drug development. However, it is extremely challenging to develop small-molecule inhibitors to target abnormally expressed transcription factors (TFs) except for the nuclear receptor family of TFs. Little is known about the interaction between TFs and transcription cofactors in gastroesophageal adenocarcinoma (GEA) or the therapeutic effects of targeting TF and transcription cofactor complexes. In this study, we found that ETS homologous factor (EHF) expression is promoted by a core transcriptional regulatory circuitry (CRC), specifically ELF3-KLF5-GATA6, and interference with its expression suppressed the malignant biological behavior of GEA cells. Importantly, we identified Ajuba LIM protein (AJUBA) as a new coactivator of EHF that cooperatively orchestrates transcriptional network activity in GEA. Furthermore, we identified KRAS signaling as a common pathway downstream of EHF and AJUBA. Applicably, dual targeting of EHF and AJUBA by lipid nanoparticles cooperatively attenuated the malignant biological behaviors of GEA in vitro and in vivo. In conclusion, EHF is upregulated by the CRC and promotes GEA malignancy by interacting with AJUBA through the KRAS pathway. Targeting of both EHF and its coactivator AJUBA through lipid nanoparticles is a novel potential therapeutic strategy. EHF is upregulated by a core transcriptional regulatory circuitry and promotes malignancy of gastroesophageal adenocarcinoma via AJUBA coactivation through the KRAS pathway, offering a novel dual-targeting therapeutic strategy. [Display omitted] [ABSTRACT FROM AUTHOR]

Background: Cis-regulatory elements (CREs) play crucial roles in regulating gene expression during erythroid cell differentiation. Genome-wide erythroid-specific CREs have not been characterized in chicken erythroid cells, which is an organism model used to study epigenetic regulation during erythropoiesis. Methods: Analysis of public genome-wide accessibility (ATAC-seq) maps, along with transcription factor (TF) motif analysis, CTCF, and RNA Pol II occupancy, as well as transcriptome analysis in fibroblasts and erythroid HD3 cells, were used to characterize erythroid-specific CREs. An a-globin CRE was identified, and its regulatory activity was validated in vitro and in vivo by luciferase activity and genome-editing assays in HD3 cells, respectively. Additionally, circular chromosome conformation capture (UMI-4C) assays were used to distinguish its role in structuring the a-globin domain in erythroid chicken cells. Results: Erythroid-specific CREs displayed occupancy by erythroid TF binding motifs, CTCF, and RNA Pol II, as well as an association with genes involved in hematopoiesis and cell differentiation. An a-globin CRE, referred to as CRE-2, was identified as exhibiting enhancer activity over aD and aA genes in vitro and in vivo. Induction of terminal erythroid differentiation showed that a-globin CRE-2 is required for the induction of aD and aA. Analysis of TF binding motifs at a-globin CRE-2 shows apparent regulation mediated by GATA-1, YY1, and CTCF binding. Conclusion: Our findings demonstrate that cell-specific CREs constitute a key mechanism that contributes to the fine-tuning gene regulation of erythroid cell differentiation and provide insights into the annotation and characterization of CREs in chicken cells. [ABSTRACT FROM AUTHOR]

ASCL1 is a transcription factor that directs neural progenitors towards lineage differentiation. Although many of the molecular mechanisms underlying its action have been described, several of its targets remain unidentified. We identified in the chick genome a putative enhancer (cE1) upstream of the transcription factor Scratch2 (Scrt2) locus with a predicted heterodimerization motif for ASCL1 and POU3F2. In this study, we investigated the role of ASCL1 and this enhancer in regulating the expression of the Scrt2 in the embryonic spinal cord. We confirmed that cE1 region interacted with the Scrt2 promoter. cE1 was sufficient to mediate ASCL1-driven expression in the neural tube through the heterodimerization sites. Moreover, Scrt2 expression was inhibited when we removed cE1 from the genome. These findings strongly indicate that ASCL1 regulates Scrt2 transcription in the neural tube through cE1. [ABSTRACT FROM AUTHOR]

Human endogenous retroviruses (HERVs) are derived from the infection and integration of exogenetic retroviruses. HERVs account for 8% of human genome, and the majority of HERVs are solitary LTRs (solo-LTRs) due to homologous recombination. Multiple findings have showed that solo-LTRs could provide an enormous reservoir of transcriptional regulatory sequences involved in diverse biological processes, especially carcinogenesis and cancer development. The link between solo-LTRs and human diseases still remains poorly understood. This review focuses on the regulatory modules of solo-LTRs, which contribute greatly to the diversification and evolution of human genes. More importantly, although inactivating mutations, insertions and deletions have been identified in solo-LTRs, the inherited regulatory elements of solo-LTRs initiate the expression of chimeric lncRNA transcripts, which have been reported to play crucial roles in human health and disease. These findings provide valuable insights into the evolutionary and functional mechanisms underlying the presence of HERVs in human genome. Taken together, in this review, we will present evidences showing the regulatory and encoding capacity of solo-LTRs as well as the significant impact on various aspects of human biology. [ABSTRACT FROM AUTHOR]

Surgical castration can effectively avoid boar taint and improve pork quality by removing the synthesis of androstenone in the testis, thereby reducing its deposition in adipose tissue. The expression of genes involved in testis-derived hormone metabolism was altered following surgical castration, but the upstream regulatory factors and underlying mechanism remain unclear. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics in liver tissue of castrated and intact full-sibling Yorkshire pigs. First, we identified 897 differentially expressed genes and 6864 differential accessible regions (DARs) using RNA- and ATAC-seq. By integrating the RNA- and ATAC-seq results, 227 genes were identified, and a significant positive correlation was revealed between differential gene expression and the ATAC-seq signal. We constructed a transcription factor regulatory network after motif analysis of DARs and identified a candidate transcription factor (TF) SP1 that targeted the HSD3B1 gene, which was responsible for the metabolism of androstenone. Subsequently, we annotated DARs by incorporating H3K27ac ChIP-seq data, marking 2234 typical enhancers and 245 super enhancers involved in the regulation of all testis-derived hormones. Among these, four typical enhancers associated with HSD3B1 were identified. Furthermore, an in-depth investigation was conducted on the androstenone-related enhancers, and an androstenone-related mutation was identified in a newfound candidatetypical enhancer (andEN) with dual-luciferase assays. These findings provide further insights into how enhancers function as links between phenotypic and non-coding area variations. The discovery of upstream TF and enhancers of HSD3B1 contributes to understanding the regulatory networks of androstenone metabolism and provides an important foundation for improving pork quality. [ABSTRACT FROM AUTHOR]

High metastatic propensity of osteosarcoma leads to its therapeutic failure and poor prognosis. Although nuclear activation miRNAs (NamiRNAs) are reported to activate gene transcription via targeting enhancer and further promote tumor metastasis, it remains uncertain whether NamiRNAs regulate osteosarcoma metastasis and their exact mechanism. Here, we found that extracellular vesicles of the malignant osteosarcoma cells (143B) remarkably increased the migratory abilities of MNNG cells representing the benign osteosarcoma cells by two folds, which attributed to their high miR‐1246 levels. Specially, miR‐1246 located in nucleus could activate the migration gene expression (such as MMP1) to accelerate MNNG cell migration through elevating the enhancer activities via increasing H3K27ac enrichment. Instead, MMP1 expression was dramatically inhibited after Argonaute 2 (AGO2) knockdown. Notably, in vitro assays demonstrated that AGO2 recognized the hybrids of miR‐1246 and its enhancer DNA via PAZ domains to prevent their degradation from RNase H and these protective roles of AGO2 may favor the gene activation by miR‐1246 in vivo. Collectively, our findings suggest that miR‐1246 could facilitate osteosarcoma metastasis through interacting with enhancer to activate gene expression dependent on AGO2, highlighting the nuclear AGO2 as a guardian for NamiRNA‐targeted gene activation and the potential of miR‐1246 for osteosarcoma metastasis therapy. [ABSTRACT FROM AUTHOR]

Human-specific insertions play important roles in human phenotypes and diseases. Here we reported a 446-bp insertion (Insert-446) in intron 11 of the TBC1D8B gene, located on chromosome X, and traced its origin to a portion of intron 6 of the EBF1 gene on chromosome 5. Interestingly, Insert-446 was present in the human Neanderthal and Denisovans genomes, and was fixed in humans after human-chimpanzee divergence. We have demonstrated that Insert-446 acts as an enhancer through binding transcript factors that promotes a higher expression of human TBC1D8B gene as compared with orthologs in macaques. In addition, over-expression TBC1D8B promoted cell proliferation and migration through "a dual finger" catalytic mechanism (Arg538 and Gln573) in the TBC domain in vitro and knockdown of TBC1D8B attenuated tumorigenesis in vivo. Knockout of Insert-446 prevented cell proliferation and migration in cancer and normal cells. Our results reveal that the human-specific Insert-446 promotes cell proliferation and migration by upregulating the expression of TBC1D8B gene. These findings provide a significant insight into the effects of human-specific insertions on evolution. [ABSTRACT FROM AUTHOR]

The β-catenin-dependent canonical Wnt signaling is pivotal in organ development, tissue homeostasis, and cancer. Here, we identified an upstream enhancer of Ctnnb1 – the coding gene for β-catenin, named ieCtnnb1 (intestinal enhancer of Ctnnb1), which is crucial for intestinal homeostasis. ieCtnnb1 is predominantly active in the base of small intestinal crypts and throughout the epithelia of large intestine. Knockout of ieCtnnb1 led to a reduction in Ctnnb1 transcription, compromising the canonical Wnt signaling in intestinal crypts. Single-cell sequencing revealed that ieCtnnb1 knockout altered epithelial compositions and potentially compromised functions of small intestinal crypts. While deletion of ieCtnnb1 hampered epithelial turnovers in physiologic conditions, it prevented occurrence and progression of Wnt/β-catenin-driven colorectal cancers. Human ieCTNNB1 drove reporter gene expression in a pattern highly similar to mouse ieCtnnb1. ieCTNNB1 contains a single-nucleotide polymorphism associated with CTNNB1 expression levels in human gastrointestinal epithelia. The enhancer activity of ieCTNNB1 in colorectal cancer tissues was stronger than that in adjacent normal tissues. HNF4α and phosphorylated CREB1 were identified as key trans-factors binding to ieCTNNB1 and regulating CTNNB1 transcription. Together, these findings unveil an enhancer-dependent mechanism controlling the dosage of Wnt signaling and homeostasis in intestinal epithelia.

Although the role of G-quadruplex (G4) DNA structures has been suggested in chromosomal looping this was not tested directly. Here, to test causal function, an array of G4s, or control sequence that does not form G4s, were inserted within chromatin in cells. In vivo G4 formation of the inserted G4 sequence array, and not the control sequence, was confirmed using G4-selective antibody. Compared to the control insert, we observed a remarkable increase in the number of 3D chromatin looping interactions from the inserted G4 array. This was evident within the immediate topologically associated domain (TAD) and throughout the genome. Locally, recruitment of enhancer histone marks and the transcriptional coactivator p300/Acetylated-p300 increased in the G4-array, but not in the control insertion. Resulting promoter-enhancer interactions and gene activation were clear up to 5 Mb away from the insertion site. Together, these show the causal role of G4s in enhancer function and long-range chromatin interactions. Mechanisms of 3D topology are primarily based on DNA-bound architectural proteins that induce/stabilize long-range interactions. Involvement of the underlying intrinsic DNA sequence/structure in 3D looping shown here therefore throws new light on how long-range chromosomal interactions might be induced or maintained.

The myeloma overexpressed gene (MYEOV) has been proposed to be a proto-oncogene due to high RNA transcript levels found in multiple cancers, including myeloma, breast, lung, pancreas and esophageal cancer. The presence of an open reading frame (ORF) in humans and other primates suggests protein-coding potential. Yet, we still lack evidence of a functional MYEOV protein. It remains undetermined how MYEOV overexpression affects cancerous tissues. In this work, we show that MYEOV has likely originated and may still function as an enhancer, regulating CCND1 and LTO1. Firstly, MYEOV 3′ enhancer activity was confirmed in humans using publicly available ATAC-STARR-seq data, performed on B-cell-derived GM12878 cells. We detected enhancer histone marks H3K4me1 and H3K27ac overlapping MYEOV in multiple healthy human tissues, which include B cells, liver and lung tissue. The analysis of 3D genome datasets revealed chromatin interactions between a MYEOV-3′-putative enhancer and the proto-oncogene CCND1. BLAST searches and multi-sequence alignment results showed that DNA sequence from this human enhancer element is conserved from the amphibians/amniotes divergence, with a 273 bp conserved region also found in all mammals, and even in chickens, where it is consistently located near the corresponding CCND1 orthologues. Furthermore, we observed conservation of an active enhancer state in the MYEOV orthologues of four non-human primates, dogs, rats, and mice. When studying this homologous region in mice, where the ORF of MYEOV is absent, we not only observed an enhancer chromatin state but also found interactions between the mouse enhancer homolog and Ccnd1 using 3D-genome interaction data. This is similar to the interaction observed in humans and, interestingly, coincides with CTCF binding sites in both species. Taken together, this suggests that MYEOV is a primate-specific gene with a de novo ORF that originated at an evolutionarily older enhancer region. This deeply conserved putative enhancer element could regulate CCND1 in both humans and mice, opening the possibility of studying MYEOV regulatory functions in cancer using non-primate animal models.

Background: Children with severe asthma suffer from recurrent symptoms and impaired quality of life despite advanced treatment. Underlying causes of severe asthma are not completely understood, although genetic mechanisms are known to be important. Objective: The aim of this study was to identify gene regulatory enhancers in leukocytes, to describe the role of these enhancers in regulating genes related to severe and mild asthma in children, and to identify known asthma-related SNPs situated in proximity to enhancers. Methods: Gene enhancers were identified and expression of enhancers and genes were measured by Cap Analysis Gene Expression (CAGE) data from peripheral blood leukocytes from children with severe asthma (n = 13), mild asthma (n = 15), and age-matched controls (n = 9). Results: From a comprehensive set of 8,289 identified enhancers, we further defined a robust sub-set of the high-confidence and most highly expressed 4,738 enhancers. Known single nucleotide polymorphisms, SNPs, related to asthma coincided with enhancers in general as well as with specific enhancer-gene interactions. Blocks of enhancer clusters were associated with genes including TGF-beta, PPAR and IL-11 signaling as well as genes related to vitamin A and D metabolism. A signature of 91 enhancers distinguished between children with severe and mild asthma as well as controls. Conclusions: Gene regulatory enhancers were identified in leukocytes with potential roles related to severe and mild asthma in children. Enhancers hosting known SNPs give the opportunity to formulate mechanistic hypotheses about the functions of these SNPs.

Unsaturated fatty acids in avocado oil can help reduce erythema brought on by prolonged UV-B exposure. The effectivity of sunscreen absorption into the skin will be enhanced by the use of isopropyl myristate (IPM) in nanoemulsion. The purpose of this study was to determine the physical characteristics and sunscreen effectiveness of avocado oil nanoemulsion (AVN) modified with IPM. 1% (FI), 3% (FII), and 5% (FIII) IPM variation were used to make AVN with 5% oil. The AVN were tested for physical characteristics such as organoleptic, pH, viscosity, rheology, particle size and polydispersity index (PI). The products were also tested for sunscreen effectivity by in vitro and Minimum Erythemal Dose (MED) method. The data obtained were analyzed statistically. The results showed that the AVN was pale yellow and clear with transmittance percentage were 96%. The rheogram showed that the products were newtonian. The pH values range were from 6.62 to 6.66; viscosity 1.65-1.84 dPa.s; particle size < 17 nm, zeta potential was in range of -30,54±1,72 to -37,85±3,11 and PI < 0.5 for all formula. In vitro SPF values were 16.43 ± 4.50 (FI), 16.27 ± 4.20 (FII) and 17.88 ± 3.20 (FIII) (p >0.05), and categorized as ultra protection. MED value were 12.28 ± 1.34 (FI); 12.51 ± 1.68 (FII); and 13.22 ± 1.84 (FIII) (p< 0.05) and categorized as maximum protection. Isopropyl myristate increased the sunscreen product's MED value without changing its physical characteristics.

Abstract Background Drug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions. Results To identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci. Conclusions We describe “primary T cell crisprQTL” — a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.

Abstract Genome-wide association studies (GWAS) are a powerful tool for detecting variants associated with complex traits and can help risk stratification and prevention strategies against pancreatic ductal adenocarcinoma (PDAC). However, the strict significance threshold commonly used makes it likely that many true risk loci are missed. Functional annotation of GWAS polymorphisms is a proven strategy to identify additional risk loci. We aimed to investigate single-nucleotide polymorphisms (SNP) in regulatory regions [transcription factor binding sites (TFBSs) and enhancers] that could change the expression profile of multiple genes they act upon and thereby modify PDAC risk. We analyzed a total of 12,636 PDAC cases and 43,443 controls from PanScan/PanC4 and the East Asian GWAS (discovery populations), and the PANDoRA consortium (replication population). We identified four associations that reached study-wide statistical significance in the overall meta-analysis: rs2472632(A) (enhancer variant, OR 1.10, 95%CI 1.06,1.13, p = 5.5 × 10−8), rs17358295(G) (enhancer variant, OR 1.16, 95%CI 1.10,1.22, p = 6.1 × 10−7), rs2232079(T) (TFBS variant, OR 0.88, 95%CI 0.83,0.93, p = 6.4 × 10−6) and rs10025845(A) (TFBS variant, OR 1.88, 95%CI 1.50,1.12, p = 1.32 × 10−5). The SNP with the most significant association, rs2472632, is located in an enhancer predicted to target the coiled-coil domain containing 34 oncogene. Our results provide new insights into genetic risk factors for PDAC by a focused analysis of polymorphisms in regulatory regions and demonstrating the usefulness of functional prioritization to identify loci associated with PDAC risk.

Abstract The breakthrough high-throughput measurement of the cis-regulatory activity of millions of randomly generated promoters provides an unprecedented opportunity to systematically decode the cis-regulatory logic that determines the expression values. We developed an end-to-end transformer encoder architecture named Proformer to predict the expression values from DNA sequences. Proformer used a Macaron-like Transformer encoder architecture, where two half-step feed forward (FFN) layers were placed at the beginning and the end of each encoder block, and a separable 1D convolution layer was inserted after the first FFN layer and in front of the multi-head attention layer. The sliding k-mers from one-hot encoded sequences were mapped onto a continuous embedding, combined with the learned positional embedding and strand embedding (forward strand vs. reverse complemented strand) as the sequence input. Moreover, Proformer introduced multiple expression heads with mask filling to prevent the transformer models from collapsing when training on relatively small amount of data. We empirically determined that this design had significantly better performance than the conventional design such as using the global pooling layer as the output layer for the regression task. These analyses support the notion that Proformer provides a novel method of learning and enhances our understanding of how cis-regulatory sequences determine the expression values.

Abstract Background The vast majority of findings from human genome-wide association studies (GWAS) map to non-coding sequences, complicating their mechanistic interpretations and clinical translations. Non-coding sequences that are evolutionarily conserved and biochemically active could offer clues to the mechanisms underpinning GWAS discoveries. However, genetic effects of such sequences have not been systematically examined across a wide range of human tissues and traits, hampering progress to fully understand regulatory causes of human complex traits. Results Here we develop a simple yet effective strategy to identify functional elements exhibiting high levels of human-mouse sequence conservation and enhancer-like biochemical activity, which scales well to 313 epigenomic datasets across 106 human tissues and cell types. Combined with 468 GWAS of European (EUR) and East Asian (EAS) ancestries, these elements show tissue-specific enrichments of heritability and causal variants for many traits, which are significantly stronger than enrichments based on enhancers without sequence conservation. These elements also help prioritize candidate genes that are functionally relevant to body mass index (BMI) and schizophrenia but were not reported in previous GWAS with large sample sizes. Conclusions Our findings provide a comprehensive assessment of how sequence-conserved enhancer-like elements affect complex traits in diverse tissues and demonstrate a generalizable strategy of integrating evolutionary and biochemical data to elucidate human disease genetics.

Abstract Background Current evidence suggests that cis-regulatory elements controlling gene expression may be the predominant target of natural selection in humans and other species. Detecting selection acting on these elements is critical to understanding evolution but remains challenging because we do not know which mutations will affect gene regulation. Results To address this, we devise an approach to search for lineage-specific selection on three critical steps in transcriptional regulation: chromatin activity, transcription factor binding, and chromosomal looping. Applying this approach to lymphoblastoid cells from 831 individuals of either European or African descent, we find strong signals of differential chromatin activity linked to gene expression differences between ancestries in numerous contexts, but no evidence of functional differences in chromosomal looping. Moreover, we show that enhancers rather than promoters display the strongest signs of selection associated with sites of differential transcription factor binding. Conclusions Overall, our study indicates that some cis-regulatory adaptation may be more easily detected at the level of chromatin than DNA sequence. This work provides a vast resource of genomic interaction data from diverse human populations and establishes a novel selection test that will benefit future study of regulatory evolution in humans and other species.

During the early development of tetrapods, including humans, the embryonic body elongates caudally once the anterior-posterior axis is established. During this process, region-specific vertebral morphogenesis occurs, with the determination of limb positioning along the anterior-posterior axis. We previously reported that Gdf11 functions as an anatomical integration system that determines the positioning of hindlimbs and sacral vertebrae where Gdf11 is expressed. However, the molecular mechanisms underlying induction of Gdf11 expression remain unclear. In this study, we searched for non-coding regions near the Gdf11 locus that were conserved across species to elucidate the regulatory mechanisms of Gdf11 expression. We identified an enhancer of the Gdf11 gene in intron 1 and named it highly conserved region (HCR). In HCR knockout mice, the expression level of endogenous Gdf11 was decreased, and the position of the sacral-hindlimb unit was shifted posteriorly. We also searched for factors upstream of Gdf11 based on the predicted transcription factor binding sites within the HCR. We found that inhibition of FGF signaling increased endogenous Gdf11 expression, suggesting that FGF signaling negatively regulates Gdf11 expression. However, FGF signaling does not regulate HCR activity. Our results suggest that there are species-specific Gdf11 enhancers other than HCR and that FGF signaling regulates Gdf11 expression independent of HCR. [ABSTRACT FROM AUTHOR]

Drosophila matrix metalloproteinase 2 (MMP2) is specifically expressed in posterior follicle cells of stage-14 egg chambers (mature follicles) and is crucial for the breakdown of the follicular wall during ovulation, a process that is highly conserved from flies to mammals. The factors that regulate spatiotemporal expression of MMP2 in follicle cells remain unknown. Here, we demonstrate crucial roles for the ETS-family transcriptional activator Pointed (Pnt) and its endogenous repressor Yan in the regulation of MMP2 expression. We found that Pnt is expressed in posterior follicle cells and overlaps with MMP2 expression in mature follicles. Genetic analysis demonstrated that pnt is both required and sufficient for MMP2 expression in follicle cells. In addition, Yan was temporally upregulated in stage-13 follicle cells to fine-tune Pnt activity and MMP2 expression. Furthermore, we identified a 1.1 kb core enhancer that is responsible for the spatiotemporal expression of MMP2 and contains multiple pnt/yan binding motifs. Mutation of pnt/yan binding sites significantly impaired the Mmp2 enhancer activity. Our data reveal a mechanism of transcriptional regulation of Mmp2 expression in Drosophila ovulation, which could be conserved in other biological systems. [ABSTRACT FROM AUTHOR]

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis -regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis -regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis -regulatory elements in A. thaliana , Drosophila melanogaster , and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis -regulatory elements of A. thaliana , whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens , as well as cis -regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites. [ABSTRACT FROM AUTHOR]

Cis-regulatory elements have an important role in human adaptation to the living environment. However, the lag in population genomic cohort studies and epigenomic studies, hinders the research in the adaptive analysis of cis-regulatory elements in human populations. In this study, we collected 4,013 unrelated individuals and performed a comprehensive analysis of adaptive selection of genome-wide cis-regulatory elements in the Han Chinese. In total, 12.34% of genomic regions are under the influence of adaptive selection, where 1.00% of enhancers and 2.06% of promoters are under positive selection, and 0.06% of enhancers and 0.02% of promoters are under balancing selection. Gene ontology enrichment analysis of these cis-regulatory elements under adaptive selection reveals that many positive selections in the Han Chinese occur in pathways involved in cell–cell adhesion processes, and many balancing selections are related to immune processes. Two classes of adaptive cis-regulatory elements related to cell adhesion were in-depth analyzed, one is the adaptive enhancers derived from neanderthal introgression, leads to lower hyaluronidase level in skin, and brings better performance on UV-radiation resistance to the Han Chinese. Another one is the cis-regulatory elements regulating wound healing, and the results suggest the positive selection inhibits coagulation and promotes angiogenesis and wound healing in the Han Chinese. Finally, we found that many pathogenic alleles, such as risky alleles of type 2 diabetes or schizophrenia, remain in the population due to the hitchhiking effect of positive selections. Our findings will help deepen our understanding of the adaptive evolution of genome regulation in the Han Chinese. [ABSTRACT FROM AUTHOR]

Hox genes encode transcription factors whose roles in patterning animal body plans during embryonic development are well-documented. Multiple studies demonstrate that Hox genes continue to act in adult cells, in normal differentiation, in regenerative processes, and, with abnormal expression, in diverse types of cancers. However, surprisingly little is known about the regulatory mechanisms that govern Hox gene expression in specific cell types, as they differentiate during late embryonic development, and in the adult organism. The murine Hoxc8 gene determines the identity of multiple skeletal elements in the lower thoracic and lumbar region and continues to play a role in the proliferation and differentiation of cells in cartilage as the skeleton matures. This study was undertaken to identify regulatory elements in the Hoxc8 gene that control transcriptional activity, specifically in cartilage-producing chondrocytes. We report that an enhancer comprising two 416 and 224 bps long interacting DNA elements produces reporter gene activity when assayed on a heterologous transcriptional promoter in transgenic mice. This enhancer is distinct in spatial, temporal, and molecular regulation from previously identified regulatory sequences in the Hoxc8 gene that control its expression in early development. The identification of a tissue-specific Hox gene regulatory element now allows mechanistic investigations into Hox transcription factor expression and function in differentiating cell types and adult tissues and to specifically target these cells during repair processes and regeneration. [ABSTRACT FROM AUTHOR]

A major goal of modern neuroscience is to understand the functions of the varied neuronal types that comprise the mammalian brain. Toward this end, some types of neurons can be targeted and manipulated with enhancer-bearing AAV vectors. These vectors hold great promise to advance basic and translational neuroscience, but to realize this potential, their selectivitymust be characterized. In this study, we investigated the selectivity of AAV vectors carrying an enhancer of the murine Dlx5 and Dlx6 genes. Vectors were injected into the visual cortex of two macaque monkeys, the frontal cortex of two others, and the somatosensory/motor cortex of three rats. Post-mortem immunostaining revealed that parvalbumin-expressing neurons were transduced efficiently in all cases but calretinin-expressing neurons were not. We speculate that this specificity is a consequence of differential activity of this DLX5/6 enhancer in adult neurons of different developmental lineages. [ABSTRACT FROM AUTHOR]

The breakthrough high-throughput measurement of the cis-regulatory activity of millions of randomly generated promoters provides an unprecedented opportunity to systematically decode the cis-regulatory logic that determines the expression values. We developed an end-to-end transformer encoder architecture named Proformer to predict the expression values from DNA sequences. Proformer used a Macaron-like Transformer encoder architecture, where two half-step feed forward (FFN) layers were placed at the beginning and the end of each encoder block, and a separable 1D convolution layer was inserted after the first FFN layer and in front of the multi-head attention layer. The sliding k-mers from one-hot encoded sequences were mapped onto a continuous embedding, combined with the learned positional embedding and strand embedding (forward strand vs. reverse complemented strand) as the sequence input. Moreover, Proformer introduced multiple expression heads with mask filling to prevent the transformer models from collapsing when training on relatively small amount of data. We empirically determined that this design had significantly better performance than the conventional design such as using the global pooling layer as the output layer for the regression task. These analyses support the notion that Proformer provides a novel method of learning and enhances our understanding of how cis-regulatory sequences determine the expression values. [ABSTRACT FROM AUTHOR]

Critical genes involved in embryonic development are often transcription factors, regulating many downstream genes. LMX1B is a homeobox gene that is involved in formation of the limbs, eyes and kidneys, heterozygous loss‐of‐function sequence variants and deletions cause Nail‐Patella syndrome. Most of the reported variants are localised within the gene's coding sequence, however, approximately 5%–10% of affected individuals do not have a pathogenic variant identified within this region. In this study, we present a family with four affected individuals across two generations with a deletion spanning a conserved upstream LMX1B‐binding sequence. This deletion is de novo in the mother of three affected children. Furthermore, in this family, the manifestations appear limited to the nails and limbs, and therefore may reflect an attenuated phenotype of the classic Nail‐Patella phenotype that includes ophthalmological and renal manifestations. [ABSTRACT FROM AUTHOR]

Chromatin remodeling is central to the dynamic changes in gene expression that drive cell fate determination. During development, the sets of enhancers that are accessible for use change globally as cells transition between stages. While transcription factors and nucleosome remodelers are known to work together to control enhancer accessibility, it is unclear how the short stretches of DNA that they individually unmask yield the kilobase-sized accessible regions characteristic of active enhancers. Here, we performed a genetic screen to investigate the role of nucleosome remodelers in control of dynamic enhancer activity. We find that the Drosophila Switch/Sucrose Non-Fermenting complex, BAP, is required for repression of a temporally dynamic enhancer, brdisc. Contrary to expectations, we find that the BAP-specific subunit Osa is dispensable for mediating changes in chromatin accessibility between the early and late stages of wing development. Instead, we find that Osa is required to constrain the levels of brdisc activity when the enhancer is normally active. Genome-wide profiling reveals that Osa directly binds brdisc as well as thousands of other developmentally dynamic regulatory sites, including multiple genes encoding components and targets of the Notch signaling pathway. Transgenic reporter analyses demonstrate that Osa is required for activation and for constraint of different sets of target enhancers in the same cells. Moreover, Osa loss results in hyperactivation of the Notch ligand Delta and development of ectopic sensory structures patterned by Notch signaling early in development. Together, these findings indicate that proper constraint of enhancer activity is necessary for regulation of dose-dependent developmental events. [ABSTRACT FROM AUTHOR]

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. Key messages: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801. [ABSTRACT FROM AUTHOR]

Tauopathies are a group of neurodegenerative diseases defined by abnormal aggregates of tau, a microtubule-associated protein encoded by MAPT. MAPT expression is near absent in neural progenitor cells (NPCs) and increases during differentiation. This temporally dynamic expression pattern suggests that MAPT expression could be controlled by transcription factors and cis -regulatory elements specific to differentiated cell types. Given the relevance of MAPT expression to neurodegeneration pathogenesis, identification of such elements is relevant to understanding disease risk and pathogenesis. Here, we performed chromatin conformation assays (HiC & Capture-C), single-nucleus multiomics (RNA-seq+ATAC-seq), bulk ATAC-seq, and ChIP-seq for H3K27ac and CTCF in NPCs and differentiated neurons to nominate candidate cis -regulatory elements (cCREs). We assayed these cCREs using luciferase assays and CRISPR interference (CRISPRi) experiments to measure their effects on MAPT expression. Finally, we integrated cCRE annotations into an analysis of genetic variation in neurodegeneration-affected individuals and control subjects. We identified both proximal and distal regulatory elements for MAPT and confirmed the regulatory function for several regions, including three regions centromeric to MAPT beyond the H1/H2 haplotype inversion breakpoint. We also found that rare and predicted damaging genetic variation in nominated CREs was nominally depleted in dementia-affected individuals relative to control subjects, consistent with the hypothesis that variants that disrupt MAPT enhancer activity, and thereby reduced MAPT expression, may be protective against neurodegenerative disease. Overall, this study provides compelling evidence for pursuing detailed knowledge of CREs for genes of interest to permit better understanding of disease risk. This study uses functional genomics approaches to assess candidate cis -regulatory elements (cCREs) for MAPT in neurons. The study suggests that impactful rare non-coding variants in MAPT cCREs (hypothesized to result in lower tau expression) may be protective against neurodegeneration, emphasizing the potential importance of rare non-coding variants in disease risk. [ABSTRACT FROM AUTHOR]