Cis-regulatory Regions Research Articles

Genome-wide association study of subfoveal choroidal thickness in a longitudinal cohort of older adults

To identify genetic influences on subfoveal choroidal thickness of older adults using a genome-wide association study (GWAS). We recruited 300 participants from the population-based Korean Longitudinal Study on Health and Aging (KLoSHA) and Korean Longitudinal Study on Cognitive Aging and Dementia (KLOSCAD) cohort studies and 500 participants from the Bundang age-related macular degeneration (AMD) cohort study dataset. We conducted a GWAS on older adult populations in the KLoSHA and KLOSCAD cohorts. Single nucleotide polymorphisms (SNPs) associated with choroidal thickness were identified with P values < 1.0 × 10−4 in both the right and left eyes, followed by validation using the Bundang AMD cohort dataset. This association was further confirmed by a functional in vitro study using human umbilical vein endothelial cells (HUVECs). The ages of the cohort participants in the discovery and validation datasets were 73.5 ± 3.3 and 71.3 ± 7.9 years, respectively. In the discovery dataset, three SNPs (rs1916762, rs7587019, and rs13320098) were significantly associated with choroidal thickness in both eyes. This association was confirmed for rs1916762 (genotypes GG, GA, and AA) and rs7587019 (genotypes GG, GA, and AA), but not for rs13320098. The mean choroidal thickness decreased by 56.7 μm (AA, 73.8%) and 31.1 μm (GA, 85.6%) compared with that of the GG genotype of rs1916762, and by 55.4 μm (AA, 74.2%) and 28.2 μm (GA, 86.7%) compared with that of the GG genotype of rs7587019. The SNPs rs1916762 and rs7587019 were located close to the FAM124B gene near its cis-regulatory region. Moreover, FAM124B was highly expressed in vascular endothelial cells. In vitro HUVEC experiments showed that the inhibition of FAM124B was associated with decreased vascular endothelial proliferation, suggesting a potential mechanism of choroidal thinning. FAM124B was identified as a susceptibility gene affecting subfoveal choroidal thickness in older adults. This gene may be involved in mechanisms underlying retinal diseases associated with altered choroidal thickness, such as age-related macular degeneration.

7690 Acute hormonal signaling-induced 3D genome organization controls transcriptional dynamics during adipocyte thermogenesis

Abstract Disclosure: Y. Zhang: None. R. Zheng: None. T. Tsuji: None. Y. Xing: None. C. Wang: None. J. Darcy: None. K. Chen: None. Y. Tseng: None. Brown adipose tissue (BAT) maintains body temperature by dissipating energy as heat, presenting a promising target against obesity and its associated metabolic disorders. When exposed to cold, the β3-adrenergic receptor (β3-AR) in brown adipocytes is stimulated, initiating a sequential signaling cascade and activating a thermogenic gene program essential for heat generation. Despite remarkable progress has been made towards understanding the engagement of transcriptional and epigenetic regulators during thermogenesis in response to β3-AR stimulation, the precise impact of three-dimensional (3D) genome organization, an increasingly important modulator of gene expression, on the activation of thermogenic gene program remains elusive. In this study, we employed high-resolution Micro-C in murine brown adipocytes to interrogate the rapidity of 3D genome reorganization in response to β3-AR hormonal signaling and its impact on transcriptional dynamics during thermogenesis. Our findings unveiled substantial rewiring of chromatin loops within 4 hours of β3-AR stimulation. These dynamic changes in chromatin loops, particularly enhancer-promoter (E-P) interactions, were coupled with gene activation during thermogenesis. Furthermore, we uncovered the mechanisms by which β3-AR hormonal signaling mediates 3D chromatin reorganization to regulate the thermogenic gene program. We found that β3-AR hormonal signaling induced the phosphorylation of p18Hamlet, a subunit of the SRCAP complex, facilitating the deposition of the histone variant H2A.Z into regulatory DNA regions. This nucleosome incorporation of H2A.Z further mediated nucleosome rearrangement, favoring the formation of nucleosome-depleted regions (NDRs) and enhancing the accessibility of regulatory DNA regions. As a result, H2A.Z chromatin incorporation facilitated interactions between cis-regulatory regions and gene promoters, such as E-P interactions, thereby activating thermogenic gene expression in response to β3-AR stimulation. Loss of H2A.Z disrupted loop dynamics, impeding gene activation and BAT thermogenic capacity. The conserved role of H2A.Z in human thermogenesis was validated, and we identified adiposity-related SNPs in human homologous loop anchors that could potentially be targeted for treating obesity. Taken together, our results unveil a previously underappreciated aspect of acute hormonal signaling-induced chromatin dynamics in orchestrating the thermogenic program in brown adipocytes. These findings not only enhance our understanding of molecular regulation of thermogenesis but also shed light on the role of rapid and dynamic chromatin reconfiguration in regulating the cellular response to external stimuli. Our study bridges a critical gap between 3D genome organization and gene activation in the thermogenic program mediated by acute β3-AR hormonal signaling. Presentation: 6/2/2024

Post-transcriptional reprogramming by thousands of mRNA untranslated regions in trypanosomes

Although genome-wide polycistronic transcription places major emphasis on post-transcriptional controls in trypanosomatids, messenger RNA cis-regulatory untranslated regions (UTRs) have remained largely uncharacterised. Here, we describe a genome-scale massive parallel reporter assay coupled with 3’-UTR-seq profiling in the African trypanosome and identify thousands of regulatory UTRs. Increased translation efficiency was associated with dosage of adenine-rich poly-purine tracts (pPuTs). An independent assessment of native UTRs using machine learning based predictions confirmed the robust correspondence between pPuTs and positive control, as did an assessment of synthetic UTRs. Those 3’-UTRs associated with upregulated expression in bloodstream-stage cells were also enriched in uracil-rich poly-pyrimidine tracts, suggesting a mechanism for developmental activation through pPuT ‘unmasking’. Thus, we describe a cis-regulatory UTR sequence ‘code’ that underpins gene expression control in the context of a constitutively transcribed genome. We conclude that thousands of UTRs post-transcriptionally reprogram gene expression profiles in trypanosomes.

PSVII-1 Widespread allele-specific expression across tissues and cells of the chicken genome

Abstract Chickens are a valuable livestock species, providing an affordable and nutritious food source worldwide through eggs and meat. Additionally, chickens are extensively used in scientific research as model organisms in virology, immunology, epigenetics, development, and conservation biology. Considering this, there are increasing efforts to deepen our understanding of their genome complexity, such as efforts by the Functional Annotation of Animal Genomes (FAANG) Consortium. Toward this goal, we are working to contribute to the annotation of the chicken transcriptome. Allele-specific expression (ASE) is an imbalance of allelic gene expression often driven by genetic and epigenetic changes in cis-regulatory regions. We, therefore, aimed to determine allelic imbalances in gene expression across 20 tissues and cells to shed light on mechanisms regulating gene expression. Using replicate paired-end (PE) RNA-seq samples, we identified around 7 million variants across 20 tissues and cells (e.g., reproductive tissues, intestine tissues, muscle, and immune tissues and cells). We applied quality control measures and filtered out monoallelic and homozygous variants. After that, we calculated read counts per allele to determine allelic expression. We found 365,894 significant ASE SNPs and 11,530 ASE genes in at least one tissue or cell type (FDR ≤ 0.05). We discovered that ASE SNPs are widely distributed throughout the chicken genome, and ASE SNPs and ASE genes showed a varied distribution across tissues and cells. Primary macrophage exhibited the most abundant, while the pectoralis major (breast muscle) had the lowest ASE genes and ASE SNPs, ranging from 7,592 to 32,505 ASE SNPs and 773 to 2,387 ASE genes. Our findings reveal several important pathways affected by allelic imbalance of expression. These pathways included fatty acid biosynthesis and regulation, cell proliferation and differentiation, cell metabolism, adipocytokine signaling, proteolysis and autophagy, and focal adhesion. In summary, our study provides insight into relevant biological processes critical to chickens and can contribute to improving genome annotation, helping to expand the transcriptomic reference source.

Altered interactions between cis-regulatory elements partially resolve BLADE-ON-PETIOLE genetic redundancy in Capsella rubella.

Duplicated genes are thought to follow one of three evolutionary trajectories that resolve their redundancy: neofunctionalization, subfunctionalization, or pseudogenization. Differences in expression patterns have been documented for many duplicated gene pairs and interpreted as evidence of subfunctionalization and a loss of redundancy. However, little is known about the functional impact of such differences and about their molecular basis. Here, we investigate the genetic and molecular basis for the partial loss of redundancy between the two BLADE-ON-PETIOLE genes BOP1 and BOP2 in red shepherd's purse (Capsella rubella) compared to Arabidopsis (Arabidopsis thaliana). While both genes remain almost fully redundant in A. thaliana, BOP1 in C. rubella can no longer ensure wild-type floral organ numbers and suppress bract formation, due to an altered expression pattern in the region of the cryptic bract primordium. We use two complementary approaches, transgenic rescue of A. thaliana atbop1 atbop2 double mutants and deletions in the endogenous AtBOP1 promoter, to demonstrate that several BOP1 promoter regions containing conserved noncoding sequences interact in a nonadditive manner to control BOP1 expression in the bract primordium and that changes in these interactions underlie the evolutionary divergence between C. rubella and A. thaliana BOP1 expression and activity. Similarly, altered interactions between cis-regulatory regions underlie the divergence in functional promoter architecture related to the control of floral organ abscission by BOP1. These findings highlight the complexity of promoter architecture in plants and suggest that changes in the interactions between cis-regulatory elements are key drivers for evolutionary divergence in gene expression and the loss of redundancy.

The non-canonically organized members of MIR395 gene family in Brassica juncea are associated with developmentally regulated, sulfate-stress responsive bidirectional promoters that exhibit orientation-dependent differential transcriptional activity

Several MICRORNA genes belonging to same family or different families are often found in homologous or non-homologous clusters. Among the various classes, head-to-head arranged genes form one of the largest categories of non-canonically organized genes. Such head-to-head arranged, non-canonically organized genes possibly share cis-regulatory region with the intergenic sequence having the potential to function as bi-directional promoter (BDP). The transcriptional regulation of head-to-head arranged genes, especially with bidirectional promoters, remains an enigma. In the past, bidirectional promoters have been characterized for a small set of protein-coding gene pairs in plants; however, to the best of our knowledge, no such study has been carried so far for MICRORNA genes. The present study thus functionally characterizes bidirectional promoters associated with members of MIR395 family, which is evolutionary conserved and is most frequently occurring cluster across plant kingdom.In Arabidopsis thaliana, the MIR395 gene family contains six members with two head-to-head arranged gene pairs- MIR395A-B and MIR395E-F. This organization was found to be conserved at seven loci for MIR395A-B, and eleven loci for MIR395E-F in five Brassica sps. Sequence analysis of the putative bidirectional promoters revealed variation in length, GC content and distribution of strict TATA-box. Comparatively higher level of conservation at both the ends of the bidirectional promoters, corresponding to ca. 250 bp upstream of 5’end of the respective MIRNA precursor, was observed. These conserved regions harbour several abiotic stress (nutrient, salt, drought) and hormone (ABA, ethylene) responsive cis-motifs. Functional characterization of putative bidirectional promoters associated with MIR395A-B and MIR395E-F from Arabidopsis and their respective orthologs from Brassica juncea (Bj_A08 MIR395A-B, Bj_B03 MIR395A-B, Bj_A07.1 MIR395E-F and Bj_A07.2 MIR395E-F) was carried out using a dual-reporter vector with β-glucuronidase (GUS) and Green Fluorescent Protein (GFP). Analysis of transcriptional regulation of the two reporter genes - GUS and GFP during developmental stages confirmed their bidirectional nature. Orientation-dependent differential reporter activity indicated asymmetric nature of the promoters. Comparison of the reporter activity amongst orthologs, paralogs and homeologs revealed regulatory diversification, an outcome expected in polyploid genomes. Interestingly, reporter gene activities driven by selected bidirectional promoters were also observed in anther and siliques apart vegetative tissues indicating role of miR395 in anther and fruit development. Finally, we evaluated the activity of reporter genes driven under transcriptional regulation of bidirectional promoters under normal and sulfate-deprived conditions which revealed asymmetric inducibility under sulfate-starvation, in agreement with the known role of miR395 in sulfate homeostasis.

Loss of cold tolerance is conferred by absence of the WRKY34 promoter fragment during tomato evolution

Natural evolution has resulted in reduced cold tolerance in cultivated tomato (Solanum lycopersicum). Herein, we perform a combined analysis of ATAC-Seq and RNA-Seq in cold-sensitive cultivated tomato and cold-tolerant wild tomato (S. habrochaites). We identify that WRKY34 has the most significant association with differential chromatin accessibility and expression patterns under cold stress. We find that a 60 bp InDel in the WRKY34 promoter causes differences in its transcription and cold tolerance among 376 tomato accessions. This 60 bp fragment contains a GATA cis-regulatory element that binds to SWIBs and GATA29, which synergistically suppress WRKY34 expression under cold stress. Moreover, WRKY34 interferes with the CBF cold response pathway through regulating transcription and protein levels. Our findings emphasize the importance of polymorphisms in cis-regulatory regions and their effects on chromatin structure and gene expression during crop evolution.

Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish

BackgroundCis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq).ResultsUsing the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions.ConclusionsThus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.

Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs

BackgroundLong-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones. Nevertheless, the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear.ResultsHere, we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical, high-altitude, and frigid regions. Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds. Analysis of selection signals revealed regions under selection for adaptation in tropical (55.5 Mb), high-altitude (43.6 Mb), and frigid (17.72 Mb) regions. The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments, including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments. Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation. These pathways included blood circulation, protein degradation, and inflammation for adaptation to tropical environments; heart and lung development, hypoxia response, and DNA damage repair for high-altitude adaptation; and thermogenesis, cold-induced vasodilation (CIVD), and the cell cycle for adaptation to frigid environments. By examining the chromatin state of the selection signatures, we identified the lung and ileum as two candidate functional tissues for environmental adaptation. Finally, we identified a mutation (chr1: G246,175,129A) in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments.ConclusionsIn this study, we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical, high-altitude, and frigid environments. Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.

Uncoupling histone modification crosstalk by engineering lysine demethylase LSD1.

Biochemical crosstalk between two or more histone modifications is often observed in epigenetic enzyme regulation, but its functional significance in cells has been difficult to discern. Previous enzymatic studies revealed that Lys14 acetylation of histone H3 can inhibit Lys4 demethylation by lysine-specific demethylase 1 (LSD1). In the present study, we engineered a mutant form of LSD1, Y391K, which renders the nucleosome demethylase activity of LSD1 insensitive to Lys14 acetylation. K562 cells with the Y391K LSD1 CRISPR knockin show decreased expression of a set of genes associated with cellular adhesion and myeloid leukocyte activation. Chromatin profiling revealed that the cis-regulatory regions of these silenced genes display a higher level of H3 Lys14 acetylation, and edited K562 cells show diminished H3 mono-methyl Lys4 near these silenced genes, consistent with a role for enhanced LSD1 demethylase activity. These findings illuminate the functional consequences of disconnecting histone modification crosstalk for a key epigenetic enzyme.

From enhancers to genome conformation: complex transcriptional control underlies expression of a single herpesviral gene.

Complex transcriptional control is a conserved feature of both eukaryotes and the viruses that infect them. Here, we illustrate this by combining high-density functional genomics, expression profiling, and viral-specific chromosome conformation capture to define with unprecedented detail the transcriptional regulation of a single gene, ORF68, from Kaposi's sarcoma-associated herpesvirus (KSHV). We first identified seven cis-regulatory regions by densely tiling the ~154 kb KSHV genome with CRISPRi. A parallel Cas9 nuclease screen indicated that three of these regions act as promoters of genes that regulate ORF68. RNA expression profiling demonstrated that three more of these regions act by either repressing or enhancing other distal viral genes involved in ORF68 transcriptional regulation. Finally, we tracked how the 3D structure of the viral genome changes during its lifecycle, revealing that these enhancing regulatory elements are physically closer to their targets when active, and that disrupting some elements caused large-scale changes to the 3D genome. These data enable us to construct a complete model revealing that the mechanistic diversity of this essential regulatory circuit matches that of human genes.

Proteogenomic Reprogramming to a Functional Human Totipotent Stem Cell State via a PARP-DUX4 Regulatory Axis.

PARP1 (ARTD1) and Tankyrases (TNKS1/TNKS2; PARP5a/5b) are poly-ADP-ribose polymerases (PARPs) with catalytic and non-catalytic functions that regulate both the genome and proteome during zygotic genome activation (ZGA), totipotent, and pluripotent embryonic stages. Here, we show that primed, conventional human pluripotent stem cells (hPSC) cultured continuously under non-specific TNKS1/TNKS2/PARP1-inhibited chemical naive reversion conditions underwent epigenetic reprogramming to clonal blastomere-like stem cells. TIRN stem cells concurrently expressed hundreds of gene targets of the ZGA-priming pioneer factor DUX4, as well as a panoply of four-cell (4C)-specific (e.g., TPRXL, HOX clusters), eight-cell (8C)-specific (e.g., DUXA, GSC, GATA6), primitive endoderm-specific (e.g., GATA4, SOX17), trophectoderm-specific (e.g., CDX2, TFAP2C), and naive epiblast-specific (e.g., DNMT3L, NANOG, POU5F1(OCT4)) factors; all in a hybrid, combinatorial single-cell manner. Mapping of proteomic and single-cell expressions of TIRN cells against human preimplantation embryo references identified them as relatively homogenous 4C-8C stage populations. Injection of TIRN cells into murine 8C-16C-staged embryos resulted in efficient totipotent-like single cell contributions of human cells to both extra-embryonic (trophectoderm, placenta) and embryonic (neural, fetal liver, hematopoietic) lineages in human-murine blastocyst and fetal chimeras. Pairing of proteome with ubiquitinome analyses of TIRN cells revealed a global shutdown of ADP-ribosylation, and a perturbed TNKS/PARP1 equilibrium which not only impacted the protein levels of hundreds of TNKS/PARP1 substrates via a rewiring of the ubiquitin-proteosome system (UPS), but also de-repressed expression of hundreds of developmental genes associated with PARP1 suppression. ChIP-Seq analysis of core NANOG-SOX2-OCT4 (NSO) pluripotency factors in TIRN cells identified reprogrammed DUX4-accessible distal and cis-regulatory enhancer regions that were co-bound by PARP1 (NSOP). These NSOP enhancer regions possessed co-binding motifs for hundreds of the same ZGA-associated, embryonic, and extraembryonic lineage-specifying pioneer factors (e.g., HOX, FOX, GATA, SOX, TBX, CDX families) that were concurrently co-expressed in TIRN cells; suggesting that PARP1 and DUX4 cooperate with NSO pluripotency core factors to regulate the epigenetic plasticity of a human totipotency program. These findings provide the first demonstration that global, proteome-wide perturbations of post-translational modifications (i.e., ADP-ribosylation, ubiquitination) can regulate epigenetic reprogramming during human embryogenesis. Totipotent TIRN stem cells will provide a valuable cell culture model for studying the proteogenomic regulation of lineage specification from human blastomere stages and may facilitate the efficient generation of human organs in interspecies chimeras.

Comprehensive Identification and Characterization of HML-9 Group in Chimpanzee Genome.

Endogenous retroviruses (ERVs) are related to long terminal repeat (LTR) retrotransposons, comprising gene sequences of exogenous retroviruses integrated into the host genome and inherited according to Mendelian law. They are considered to have contributed greatly to the evolution of host genome structure and function. We previously characterized HERV-K HML-9 in the human genome. However, the biological function of this type of element in the genome of the chimpanzee, which is the closest living relative of humans, largely remains elusive. Therefore, the current study aims to characterize HML-9 in the chimpanzee genome and to compare the results with those in the human genome. Firstly, we report the distribution and genetic structural characterization of the 26 proviral elements and 38 solo LTR elements of HML-9 in the chimpanzee genome. The results showed that the distribution of these elements displayed a non-random integration pattern, and only six elements maintained a relatively complete structure. Then, we analyze their phylogeny and reveal that the identified elements all cluster together with HML-9 references and with those identified in the human genome. The HML-9 integration time was estimated based on the 2-LTR approach, and the results showed that HML-9 elements were integrated into the chimpanzee genome between 14 and 36 million years ago and into the human genome between 18 and 49 mya. In addition, conserved motifs, cis-regulatory regions, and enriched PBS sequence features in the chimpanzee genome were predicted based on bioinformatics. The results show that pathways significantly enriched for ERV LTR-regulated genes found in the chimpanzee genome are closely associated with disease development, including neurological and neurodevelopmental psychiatric disorders. In summary, the identification, characterization, and genomics of HML-9 presented here not only contribute to our understanding of the role of ERVs in primate evolution but also to our understanding of their biofunctional significance.

Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization

BackgroundAs genomic studies continue to implicate non-coding sequences in disease, testing the roles of these variants requires insights into the cell type(s) in which they are likely to be mediating their effects. Prior methods for associating non-coding variants with cell types have involved approaches using linkage disequilibrium or ontological associations, incurring significant processing requirements. GaiaAssociation is a freely available, open-source software that enables thousands of genomic loci implicated in a phenotype to be tested for enrichment at regulatory loci of multiple cell types in minutes, permitting insights into the cell type(s) mediating the studied phenotype.ResultsIn this work, we present Regulatory Landscape Enrichment Analysis (RLEA) by GaiaAssociation and demonstrate its capability to test the enrichment of 12,133 variants across the cis-regulatory regions of 44 cell types. This analysis was completed in 134.0 ± 2.3 s, highlighting the efficient processing provided by GaiaAssociation. The intuitive interface requires only four inputs, offers a collection of customizable functions, and visualizes variant enrichment in cell-type regulatory regions through a heatmap matrix. GaiaAssociation is available on PyPi for download as a command line tool or Python package and the source code can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation.ConclusionsGaiaAssociation is a novel package that provides an intuitive and efficient resource to understand the enrichment of non-coding variants across the cis-regulatory regions of different cells, empowering studies seeking to identify disease-mediating cell types.

Developmental origins shape the paediatric cancer genome.

In the past two decades, technological advances have brought unprecedented insights into the paediatric cancer genome revealing characteristics distinct from those of adult cancer. Originating from developing tissues, paediatric cancers generally have low mutation burden and are driven by variants that disrupt the transcriptional activity, chromatin state, non-coding cis-regulatory regions and other biological functions. Within each tumour, there are multiple populations of cells with varying states, and the lineages of some can be tracked to their fetal origins. Genome-wide genetic screening has identified vulnerabilities associated with both the cell of origin and transcription deregulation in paediatric cancer, which have become a valuable resource for designing new therapeutic approaches including those for small molecules, immunotherapy and targeted protein degradation. In this Review, we present recent findings on these facets of paediatric cancer from a pan-cancer perspective and provide an outlook on future investigations.

Regional-specific calibration enables application of computational evidence for clinical classification of 5′ cis-regulatory variants in Mendelian disease

To date, clinical genetic testing for Mendelian disease variants has focused heavily on exonic coding and intronic gene regions. This multi-step study was undertaken to provide an evidence base for selecting and applying computational approaches for use in clinical classification of 5' cis-regulatory region variants. Curated datasets of clinically reported disease-causing 5' cis-regulatory region variants and variants from matched genomic regions in population controls were used to calibrate six bioinformatic tools as predictors of variant pathogenicity. Likelihood ratio estimates were aligned to code weights following ClinGen recommendations for application of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) classification scheme. Considering code assignment across all reference dataset variants, performance was best for CADD (81.2%) and REMM (81.5%). Optimized thresholds provided moderate evidence toward pathogenicity (CADD, REMM) and moderate (CADD) or supporting (REMM) evidence against pathogenicity. Both sensitivity and specificity of prediction were improved when further categorizing variants based on location in an EPDnew-defined promoter region. Combining predictions (CADD, REMM, and location in a promoter region) increased specificity at the expense of sensitivity. Importantly, the optimal CADD thresholds for assigning ACMG/AMP codes PP3 (≥10) and BP4 (≤8) were vastly different from recommendations for protein-coding variants (PP3≥25.3; BP4≤22.7); CADD <22.7 would incorrectly assign BP4 for >90% of reported disease-causing cis-regulatory region variants. Our results demonstrate the need to consider a tiered approach and tailored score thresholds to optimize bioinformatic impact prediction for clinical classification of 5' cis-regulatory region variants.

Simultaneous profiling of chromatin accessibility and DNA methylation in complete plant genomes using long-read sequencing.

Epigenetic regulations, including chromatin accessibility, nucleosome positioningand DNA methylation intricately shape genome function. However, current chromatin profiling techniques relying on short-read sequencing technologies fail to characterise highly repetitive genomic regions and cannot detect multiple chromatin features simultaneously. Here, we performed Simultaneous Accessibility and DNA Methylation Sequencing (SAM-seq) of purified plant nuclei. Thanks to the use of long-read nanopore sequencing, SAM-seq enables high-resolution profiling of m6A-tagged chromatin accessibility together with endogenous cytosine methylation in plants. Analysis of naked genomic DNA revealed significant sequence preference biases of m6A-MTases, controllable through a normalisation step. By applying SAM-seq to Arabidopsis and maize nuclei we obtained fine-grained accessibility and DNA methylation landscapes genome-wide. We uncovered crosstalk between chromatin accessibility and DNA methylation within nucleosomes of genes, TEs, and centromeric repeats. SAM-seq also detects DNA footprints over cis-regulatory regions. Furthermore, using the single-molecule information provided by SAM-seq we identified extensive cellular heterogeneity at chromatin domains with antagonistic chromatin marks, suggesting that bivalency reflects cell-specific regulations. SAM-seq is a powerful approach to simultaneously study multiple epigenetic features over unique and repetitive sequences, opening new opportunities for the investigation of epigenetic mechanisms.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

Polymorphic locus rs8023580 NR2F2-AS1 is associated with breast cancer risk in residents of Central Russia

Breast cancer (BC) is the most common hormone-dependent genetically determined cancer among women. The level of “active” sex hormones in the body, the connection of which with breast cancer is beyond doubt, is determined by the content of the protein transporting sex hormones (SHBG).Aim: To study associations of polymorphic loci linked with the level of SHBG at the full-genomic level of significance with the risk of developing breast cancer.Material and Methods. The work was carried out on a sample of 1,498 women, 358 of whom were BC patients, the control was 1140 individuals. Genotyping of four single nucleotide polymorphic loci (SNP) associated with the level of SHBG was performed according to previously performed genome-wide studies (GWAS): rs7910927 JMJD1C, rs4149056 SLCO1B1, rs8023580 NR2F2-AS1, rs12150660 SHBG. The method of logistic regression was used to search for associations.Results and Discussion. SNP rs8023580 of the NR2F2-AS1 gene is associated with the risk of BC developing. The presence of the CC genotype rs8023580 NR2F2-AS1in a woman has a protective value in the formation of the disease (CCvsTC+TT [recessive model]; OR = 0.58; 95%CI = 0.35–0.96; p = 0.033; pperm = 0.042). SNP rs8023580 NR2F2-AS1 is functionally significant in the liver: it is localized in the enhancer region, affects the level of methylation of the cg01739960 (hg38) genome region, affects the expression of the RP11-327J17.2 gene, determines the interaction of DNA with six transcription factors (Fo xd1,Foxl1,Foxq1,Mef2,PLZF,STAT), which are located in the cis-regulatory region of RNA polymerase II, specific for the DNA binding sequence, and determine the activity of DNA-binding transcription factors specific for RNA polymerase II, are involved in the processes of cellular differentiation and tissue development.Conclusion. The SHBG-enhancing CC genotype rs8023580 of the NR2F2-AS1 gene is a protective factor in the development of BC.

Reprogramming of the LXRα Transcriptome Sustains Macrophage Secondary Inflammatory Responses.

Macrophages regulate essential aspects of innate immunity against pathogens. In response to microbial components, macrophages activate primary and secondary inflammatory gene programs crucial for host defense. The liver X receptors (LXRα, LXRβ) are ligand-dependent nuclear receptors that direct gene expression important for cholesterol metabolism and inflammation, but little is known about the individual roles of LXRα and LXRβ in antimicrobial responses. Here, the results demonstrate that induction of LXRα transcription by prolonged exposure to lipopolysaccharide (LPS) supports inflammatory gene expression in macrophages. LXRα transcription is induced by NF-κB and type-I interferon downstream of TLR4 activation. Moreover, LPS triggers a reprogramming of the LXRα cistrome that promotes cytokine and chemokine gene expression through direct LXRα binding to DNA consensus sequences within cis-regulatory regions including enhancers. LXRα-deficient macrophages present fewer binding of p65 NF-κB and reduced histone H3K27 acetylation at enhancers of secondary inflammatory response genes. Mice lacking LXRα in the hematopoietic compartment show impaired responses to bacterial endotoxin in peritonitis models, exhibiting reduced neutrophil infiltration and decreased expansion and inflammatory activation of recruited F4/80lo-MHC-IIhi peritoneal macrophages. Together, these results uncover a previously unrecognized function for LXRα-dependent transcriptional cis-activation of secondary inflammatory gene expression in macrophages and the host response to microbial ligands.