Regulatory genomic circuitry of human disease loci by integrative epigenomics

Genome-Wide Association Study of Breast Density among Women of African Ancestry

Breast density, the amount of fibroglandular versus fatty tissue in the breast, is a strong breast cancer risk factor. Understanding genetic factors associated with breast density may help in clarifying mechanisms by which breast density increases cancer risk. To date, 50 genetic loci have been associated with breast density, however, these studies were performed among predominantly European ancestry populations. We utilized a cohort of women aged 40-85 years who underwent screening mammography and had genetic information available from the Penn Medicine BioBank to conduct a Genome-Wide Association Study (GWAS) of breast density among 1323 women of African ancestry. For each mammogram, the publicly available "LIBRA" software was used to quantify dense area and area percent density. We identified 34 significant loci associated with dense area and area percent density, with the strongest signals in GACAT3, CTNNA3, HSD17B6, UGDH, TAAR8, ARHGAP10, BOD1L2, and NR3C2. There was significant overlap between previously identified breast cancer SNPs and SNPs identified as associated with breast density. Our results highlight the importance of breast density GWAS among diverse populations, including African ancestry populations. They may provide novel insights into genetic factors associated with breast density and help in elucidating mechanisms by which density increases breast cancer risk.

An atlas of regulatory elements in chicken: A resource for chicken genetics and genomics

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

hipFG: High-throughput harmonization and integration pipeline for functional genomics data

Preparing functional genomic (FG) data with diverse assay types and file formats for integration into analysis workflows that interpret genome-wide association and other studies is a significant and time-consuming challenge. Here we introduce hipFG, an automatically customized pipeline for efficient and scalable normalization of heterogenous FG data collections into standardized, indexed, rapidly searchable analysis-ready datasets while accounting for FG datatypes (e.g., chromatin interactions, genomic intervals, quantitative trait loci).

The ENCODE Imputation Challenge: a critical assessment of methods for cross-cell type imputation of epigenomic profiles

A promising alternative to comprehensively performing genomics experiments is to, instead, perform a subset of experiments and use computational methods to impute the remainder. However, identifying the best imputation methods and what measures meaningfully evaluate performance are open questions. We address these questions by comprehensively analyzing 23 methods from the ENCODE Imputation Challenge. We find that imputation evaluations are challenging and confounded by distributional shifts from differences in data collection and processing over time, the amount of available data, and redundancy among performance measures. Our analyses suggest simple steps for overcoming these issues and promising directions for more robust research.

Characterization of genome-wide STR variation in 6487 human genomes

Short tandem repeats (STRs) are abundant and highly mutagenic in the human genome. Many STR loci have been associated with a range of human genetic disorders. However, most population-scale studies on STR variation in humans have focused on European ancestry cohorts or are limited by sequencing depth. Here, we depicted a comprehensive map of 366,013 polymorphic STRs (pSTRs) constructed from 6487 deeply sequenced genomes, comprising 3983 Chinese samples (~31.5x, NyuWa) and 2504 samples from the 1000 Genomes Project (~33.3x, 1KGP). We found that STR mutations were affected by motif length, chromosome context and epigenetic features. We identified 3273 and 1117 pSTRs whose repeat numbers were associated with gene expression and 3'UTR alternative polyadenylation, respectively. We also implemented population analysis, investigated population differentiated signatures, and genotyped 60 known disease-causing STRs. Overall, this study further extends the scale of STR variation in humans and propels our understanding of the semantics of STRs.

Comparative single cell epigenomic analysis of gene regulatory programs in the rodent and primate neocortex

Sequence divergence of cis- regulatory elements drives species-specific traits, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains to be elucidated. We investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset, and mouse with single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome, and chromosomal conformation profiles from a total of over 180,000 cells. For each modality, we determined species-specific, divergent, and conserved gene expression and epigenetic features at multiple levels. We find that cell type-specific gene expression evolves more rapidly than broadly expressed genes and that epigenetic status at distal candidate cis -regulatory elements (cCREs) evolves faster than promoters. Strikingly, transposable elements (TEs) contribute to nearly 80% of the human-specific cCREs in cortical cells. Through machine learning, we develop sequence-based predictors of cCREs in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Lastly, we show that epigenetic conservation combined with sequence similarity helps uncover functional cis -regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

CRISPR screens identify gene targets at breast cancer risk loci

BACKGROUND

Genome-wide association studies (GWAS) have identified > 200 loci associated with breast cancer risk. The majority of candidate causal variants are in non-coding regions and likely modulate cancer risk by regulating gene expression. However, pinpointing the exact target of the association, and identifying the phenotype it mediates, is a major challenge in the interpretation and translation of GWAS.

RESULTS

Here, we show that pooled CRISPR screens are highly effective at identifying GWAS target genes and defining the cancer phenotypes they mediate. Following CRISPR mediated gene activation or suppression, we measure proliferation in 2D, 3D, and in immune-deficient mice, as well as the effect on DNA repair. We perform 60 CRISPR screens and identify 20 genes predicted with high confidence to be GWAS targets that promote cancer by driving proliferation or modulating the DNA damage response in breast cells. We validate the regulation of a subset of these genes by breast cancer risk variants.

CONCLUSIONS

We demonstrate that phenotypic CRISPR screens can accurately pinpoint the gene target of a risk locus. In addition to defining gene targets of risk loci associated with increased breast cancer risk, we provide a platform for identifying gene targets and phenotypes mediated by risk variants.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Major cell-types in multiomic single-nucleus datasets impact statistical modeling of links between regulatory sequences and target genes

Epigenomic profiling, including ATACseq, is one of the main tools used to define enhancers. Because enhancers are overwhelmingly cell-type specific, inference of their activity is greatly limited in complex tissues. Multiomic assays that probe in the same nucleus both the open chromatin landscape and gene expression levels enable the study of correlations (links) between these two modalities. Current best practices to infer the regulatory effect of candidate cis-regulatory elements (cCREs) in multiomic data involve removing biases associated with GC content by generating null distributions of matched ATACseq peaks drawn from different chromosomes. This strategy has been broadly adopted by popular single-nucleus multiomic workflows such as Signac. Here, we uncovered limitations and confounders of this approach. We found a strong loss of power to detect a regulatory effect for cCREs with high read counts in the dominant cell-type. We showed that this is largely due to cell-type-specific trans-ATACseq peak correlations creating bimodal null distributions. We tested alternative models and concluded that physical distance and/or the raw Pearson correlation coefficients are the best predictors for peak-gene links when compared to predictions from Epimap (e.g. CD14 area under the curve [AUC] = 0.51 with the method implemented in Signac vs. 0.71 with the Pearson correlation coefficients) or validation by CRISPR perturbations (AUC = 0.63 vs. 0.73).

A unifying statistical framework to discover disease genes from GWASs

Genome-wide association studies (GWASs) identify genomic loci associated with complex traits, but it remains a challenge to identify the genes affected by causal genetic variants in these loci. Attempts to solve this challenge are frustrated by a number of compounding problems. Here, we show how to combine solutions to these problems into a unified mathematical framework. From this synthesis, it becomes possible to compute the probability that each gene in the genome is affected by a causal variant, given a particular trait, without making assumptions about the relevant cell types or tissues. We validate each component of the framework individually and in combination. When applied to large GWASs of human disease, the resulting paradigm can rediscover the majority of well-known disease genes. Moreover, it establishes human genetics support for many genes previously implicated only by clinical or preclinical evidence, and it uncovers a plethora of novel disease genes with compelling biological rationale.

Effect of tissue-grouped regulatory variants associated to type 2 diabetes in related secondary outcomes

Genome-wide association studies have identified over five hundred loci that contribute to variation in type 2 diabetes (T2D), an established risk factor for many diseases. However, the mechanisms and extent through which these loci contribute to subsequent outcomes remain elusive. We hypothesized that combinations of T2D-associated variants acting on tissue-specific regulatory elements might account for greater risk for tissue-specific outcomes, leading to diversity in T2D disease progression. We searched for T2D-associated variants acting on regulatory elements and expression quantitative trait loci (eQTLs) in nine tissues. We used T2D tissue-grouped variant sets as genetic instruments to conduct 2-Sample Mendelian Randomization (MR) in ten related outcomes whose risk is increased by T2D using the FinnGen cohort. We performed PheWAS analysis to investigate whether the T2D tissue-grouped variant sets had specific predicted disease signatures. We identified an average of 176 variants acting in nine tissues implicated in T2D, and an average of 30 variants acting on regulatory elements that are unique to the nine tissues of interest. In 2-Sample MR analyses, all subsets of regulatory variants acting in different tissues were associated with increased risk of the ten secondary outcomes studied on similar levels. No tissue-grouped variant set was associated with an outcome significantly more than other tissue-grouped variant sets. We did not identify different disease progression profiles based on tissue-specific regulatory and transcriptome information. Bigger sample sizes and other layers of regulatory information in critical tissues may help identify subsets of T2D variants that are implicated in certain secondary outcomes, uncovering system-specific disease progression.

Cas9-orthologue-mediated cytosine and adenine base editors recognizing NNAAAA PAM sequences

CRISPR/Cas9 system has been applied as an effective genome-targeting technology. By fusing deaminases with Cas9 nickase (nCas9), various cytosine and adenine base editors (CBEs and ABEs) have been successfully developed that can efficiently induce nucleotide conversions and install pathogenic single nucleotide variants (SNVs) in cultured cells and animal models. However, the applications of BEs are frequently limited by the specific protospacer adjacent motif (PAM) sequences and protein sizes. To expand the toolbox for BEs that can recognize novel PAM sequences, we cloned a Cas9 ortholog from Streptococcus sinensis (named as SsiCas9) with a smaller size and constructed it into APOBEC1- or APOBEC3A-composed CBEs and TadA or TadA*-composed ABEs, which yield high editing efficiencies, low off-targeting activities, and low indel rates in human cells. Compared to PAMless SpRY Cas9-composed BE4max, SsiCas9-mediated BE4max displayed higher editing efficiencies for targets with "NNAAAA" PAM sequences. Moreover, SsiCas9-mediated BE4max induced highly efficient C-to-T conversions in the mouse Ar gene (R841C) to introduce a human androgen resistance syndrome-related mutation (AR R820C) in early mouse embryos. Thus, we developed novel BEs mediated by SsiCas9, expanded the toolbox for base conversions, and broadened the range of editable genomes in vitro and in vivo.

Genome-wide association study using whole-genome sequencing identifies risk loci for Parkinson's disease in Chinese population

Genome-wide association studies (GWASs) have identified numerous susceptibility loci for Parkinson's disease (PD), but its genetic architecture remains underexplored in populations of non-European ancestry. To identify genetic variants associated with PD in the Chinese population, we performed a GWAS using whole-genome sequencing (WGS) in 1,972 cases and 2,478 controls, and a replication study in a total of 8209 cases and 9454 controls. We identified one new risk variant rs61204179 (Pcombined = 1.47 × 10-9) with low allele frequency, four previously reported risk variants (NUCKS1/RAB29-rs11557080, SNCA-rs356182, FYN-rs997368, and VPS13C-rs2251086), as well as three risk variants in LRRK2 coding region (A419V, R1628P, and G2385R) with genome-wide significance (P < 5 × 10-8) for PD in Chinese population. Moreover, of the reported genome-wide significant risk variants found mostly in European ancestry populations, the correlation coefficient (rb) of effect size accounting for sampling errors was 0.91 between datasets and 63.6% attained P < 0.05 in Chinese population. Accordingly, we estimated a heritability of 0.14-0.18 for PD, and a moderate genetic correlation between European ancestry and Chinese populations (rg = 0.47, se = 0.21). Polygenic risk score (PRS) analysis revealed that individuals with PRS values in the highest quartile had a 3.9-fold higher risk of developing PD than the lowest quartile. In conclusion, the present GWAS identified PD-associated variants in Chinese population, as well as genetic factors shared among distant populations. Our findings shed light on the genetic homogeneity and heterogeneity of PD in different ethnic groups and suggested WGS might continue to improve our understanding of the genetic architecture of PD.

Promises and challenges in pharmacoepigenetics

Pharmacogenetics, the study of how interindividual genetic differences affect drug response, does not explain all observed heritable variance in drug response. Epigenetic mechanisms, such as DNA methylation, and histone acetylation may account for some of the unexplained variances. Epigenetic mechanisms modulate gene expression and can be suitable drug targets and can impact the action of nonepigenetic drugs. Pharmacoepigenetics is the field that studies the relationship between epigenetic variability and drug response. Much of this research focuses on compounds targeting epigenetic mechanisms, called epigenetic drugs, which are used to treat cancers, immune disorders, and other diseases. Several studies also suggest an epigenetic role in classical drug response; however, we know little about this area. The amount of information correlating epigenetic biomarkers to molecular datasets has recently expanded due to technological advances, and novel computational approaches have emerged to better identify and predict epigenetic interactions. We propose that the relationship between epigenetics and classical drug response may be examined using data already available by (1) finding regions of epigenetic variance, (2) pinpointing key epigenetic biomarkers within these regions, and (3) mapping these biomarkers to a drug-response phenotype. This approach expands on existing knowledge to generate putative pharmacoepigenetic relationships, which can be tested experimentally. Epigenetic modifications are involved in disease and drug response. Therefore, understanding how epigenetic drivers impact the response to classical drugs is important for improving drug design and administration to better treat disease.

GENCODE: reference annotation for the human and mouse genomes in 2023

GENCODE produces high quality gene and transcript annotation for the human and mouse genomes. All GENCODE annotation is supported by experimental data and serves as a reference for genome biology and clinical genomics. The GENCODE consortium generates targeted experimental data, develops bioinformatic tools and carries out analyses that, along with externally produced data and methods, support the identification and annotation of transcript structures and the determination of their function. Here, we present an update on the annotation of human and mouse genes, including developments in the tools, data, analyses and major collaborations which underpin this progress. For example, we report the creation of a set of non-canonical ORFs identified in GENCODE transcripts, the LRGASP collaboration to assess the use of long transcriptomic data to build transcript models, the progress in collaborations with RefSeq and UniProt to increase convergence in the annotation of human and mouse protein-coding genes, the propagation of GENCODE across the human pan-genome and the development of new tools to support annotation of regulatory features by GENCODE. Our annotation is accessible via Ensembl, the UCSC Genome Browser and https://www.gencodegenes.org.

Single-cell omics: A new direction for functional genetic research in human diseases and animal models

Over the past decade, with the development of high-throughput single-cell sequencing technology, single-cell omics has been emerged as a powerful tool to understand the molecular basis of cellular mechanisms and refine our knowledge of diverse cell states. They can reveal the heterogeneity at different genetic layers and elucidate their associations by multiple omics analysis, providing a more comprehensive genetic map of biological regulatory networks. In the post-GWAS era, the molecular biological mechanisms influencing human diseases will be further elucidated by single-cell omics. This review mainly summarizes the development and trend of single-cell omics. This involves single-cell omics technologies, single-cell multi-omics technologies, multiple omics data integration methods, applications in various human organs and diseases, classic laboratory cell lines, and animal disease models. The review will reveal some perspectives for elucidating human diseases and constructing animal models.

A framework for group-wise summarization and comparison of chromatin state annotations

MOTIVATION

Genome-wide maps of epigenetic modifications are powerful resources for non-coding genome annotation. Maps of multiple epigenetics marks have been integrated into cell or tissue type-specific chromatin state annotations for many cell or tissue types. With the increasing availability of multiple chromatin state maps for biologically similar samples, there is a need for methods that can effectively summarize the information about chromatin state annotations within groups of samples and identify differences across groups of samples at a high resolution.

RESULTS

We developed CSREP, which takes as input chromatin state annotations for a group of samples. CSREP then probabilistically estimates the state at each genomic position and derives a representative chromatin state map for the group. CSREP uses an ensemble of multi-class logistic regression classifiers that predict the chromatin state assignment of each sample given the state maps from all other samples. The difference in CSREP's probability assignments for the two groups can be used to identify genomic locations with differential chromatin state assignments. Using groups of chromatin state maps of a diverse set of cell and tissue types, we demonstrate the advantages of using CSREP to summarize chromatin state maps and identify biologically relevant differences between groups at a high resolution.

AVAILABILITY AND IMPLEMENTATION

The CSREP source code and generated data are available at http://github.com/ernstlab/csrep.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Characterizing cis-regulatory elements using single-cell epigenomics

Cell type-specific gene expression patterns and dynamics during development or in disease are controlled by cis-regulatory elements (CREs), such as promoters and enhancers. Distinct classes of CREs can be characterized by their epigenomic features, including DNA methylation, chromatin accessibility, combinations of histone modifications and conformation of local chromatin. Tremendous progress has been made in cataloguing CREs in the human genome using bulk transcriptomic and epigenomic methods. However, single-cell epigenomic and multi-omic technologies have the potential to provide deeper insight into cell type-specific gene regulatory programmes as well as into how they change during development, in response to environmental cues and through disease pathogenesis. Here, we highlight recent advances in single-cell epigenomic methods and analytical tools and discuss their readiness for human tissue profiling.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

PEAC-seq adopts Prime Editor to detect CRISPR off-target and DNA translocation

CRISPR technology holds significant promise for biological studies and gene therapies because of its high flexibility and efficiency when applied in mammalian cells. But endonuclease (e.g., Cas9) potentially generates undesired edits; thus, there is an urgent need to comprehensively identify off-target sites so that the genotoxicities can be accurately assessed. To date, it is still challenging to streamline the entire process to specifically label and efficiently enrich the cleavage sites from unknown genomic locations. Here we develop PEAC-seq, in which we adopt the Prime Editor to insert a sequence-optimized tag to the editing sites and enrich the tagged regions with site-specific primers for high throughput sequencing. Moreover, we demonstrate that PEAC-seq could identify DNA translocations, which are more genotoxic but usually overlooked by other off-target detection methods. As PEAC-seq does not rely on exogenous oligodeoxynucleotides to label the editing site, we also conduct in vivo off-target identification as proof of concept. In summary, PEAC-seq provides a comprehensive and streamlined strategy to identify CRISPR off-targeting sites in vitro and in vivo, as well as DNA translocation events. This technique further diversified the toolkit to evaluate the genotoxicity of CRISPR applications in research and clinics.

Integration of 3D genome topology and local chromatin features uncovers enhancers underlying craniofacial-specific cartilage defects

Aberrations in tissue-specific enhancers underlie many developmental defects. Disrupting a noncoding region distal from the human SOX9 gene causes the Pierre Robin sequence (PRS) characterized by the undersized lower jaw. Such a craniofacial-specific defect has been previously linked to enhancers transiently active in cranial neural crest cells (CNCCs). We demonstrate that the PRS region also strongly regulates Sox9 in CNCC-derived Meckel's cartilage (MC), but not in limb cartilages, even after decommissioning of CNCC enhancers. Such an MC-specific regulatory effect correlates with the MC-specific chromatin contacts between the PRS region and Sox9, highlighting the importance of lineage-dependent chromatin topology in instructing enhancer usage. By integrating the enhancer signatures and chromatin topology, we uncovered >10,000 enhancers that function differentially between MC and limb cartilages and demonstrated their association with human diseases. Our findings provide critical insights for understanding the choreography of gene regulation during development and interpreting the genetic basis of craniofacial pathologies.

Associations between NKX2-5 gene polymorphisms and congenital heart disease in the Chinese Tibetan population

BACKGROUND

The pathogenesis of congenital heart disease (CHD) has not been fully elucidated, and this study considers the interaction between inheritance and the environment as the main cause of CHD. Previous studies have found that the incidence of CHD in the Tibetan plateau population is significantly higher than in low-altitude populations. Numerous reports have confirmed that NKX2-5 gene mutations can lead to coronary heart disease, but the relationship between NKX2-5 and Tibetan nationality has not yet been reported.

OBJECTIVE

To explore the relationship between NKX2-5 gene polymorphisms and CHD in Tibetan people.

METHODS

Blood samples were collected retrospectively from Tibetan patients diagnosed with CHD as well as healthy Tibetans, and the exons of NKX2-5 were sequenced. The MassARRAY technique was used to detect and genotype candidate tag single nucleotide polymorphisms (SNPs) in the non-coding regions of NKX2-5.

RESULTS

Exon sequencing revealed no difference in the coding regions of the NKX2-5 gene between the CHD and control groups. In the non-coding regions of NKX2-5, rs6882776 and rs2546741 differed significantly between the two groups. Strong linkage disequilibrium was found between the selected sites of NKX2-5.

CONCLUSIONS

The NKX2-5 exons do not associate with CHD in Tibetans. Rs6882776 and rs2546741 in the non-coding regions of NKX2-5 may protect against CHD in Tibetans. The NKX2-5 haplotype associated with CHD occurrence in the Tibetan population.

Genetics of Primary Biliary Cholangitis

Primary biliary cholangitis (PBC) is a rare disease of the liver characterized by an autoimmune attack on the small bile ducts. PBC is a complex trait, meaning that a large list of genetic factors interacts with environmental agents to determine its onset. Genome-wide association studies have had a huge impact in fostering research in PBC, but many steps need still to be done compared with other autoimmune diseases of similar prevalence. This review presents the state-of-the-art regarding the genetic architecture of PBC and provides some thoughtful reflections about possible future lines of research, which can be helpful to fill the missing heritability gap in PBC.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

From novel discovery tools and biomarkers to precision medicine-basic cardiovascular science highlights of 2021/22

Here, we review the highlights of cardiovascular basic science published in 2021 and early 2022 on behalf of the European Society of Cardiology Council for Basic Cardiovascular Science. We begin with non-coding RNAs which have emerged as central regulators cardiovascular biology, and then discuss how technological developments in single-cell 'omics are providing new insights into cardiovascular development, inflammation, and disease. We also review recent discoveries on the biology of extracellular vesicles in driving either protective or pathogenic responses. The Nobel Prize in Physiology or Medicine 2021 recognized the importance of the molecular basis of mechanosensing and here we review breakthroughs in cardiovascular sensing of mechanical force. We also summarize discoveries in the field of atherosclerosis including the role of clonal haematopoiesis of indeterminate potential, and new mechanisms of crosstalk between hyperglycaemia, lipid mediators, and inflammation. The past 12 months also witnessed major advances in the field of cardiac arrhythmia including new mechanisms of fibrillation. We also focus on inducible pluripotent stem cell technology which has demonstrated disease causality for several genetic polymorphisms in long-QT syndrome and aortic valve disease, paving the way for personalized medicine approaches. Finally, the cardiovascular community has continued to better understand COVID-19 with significant advancement in our knowledge of cardiovascular tropism, molecular markers, the mechanism of vaccine-induced thrombotic complications and new anti-viral therapies that protect the cardiovascular system.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.

Scalable approaches for functional analyses of whole-genome sequencing non-coding variants

Non-coding genetic variants outside of protein-coding genome regions play an important role in genetic and epigenetic regulation. It has become increasingly important to understand their roles, as non-coding variants often make up the majority of top findings of genome-wide association studies (GWAS). In addition, the growing popularity of disease-specific whole-genome sequencing (WGS) efforts expands the library of and offers unique opportunities for investigating both common and rare non-coding variants, which are typically not detected in more limited GWAS approaches. However, the sheer size and breadth of WGS data introduce additional challenges to predicting functional impacts in terms of data analysis and interpretation. This review focuses on the recent approaches developed for efficient, at-scale annotation and prioritization of non-coding variants uncovered in WGS analyses. In particular, we review the latest scalable annotation tools, databases and functional genomic resources for interpreting the variant findings from WGS based on both experimental data and in silico predictive annotations. We also review machine learning-based predictive models for variant scoring and prioritization. We conclude with a discussion of future research directions which will enhance the data and tools necessary for the effective functional analyses of variants identified by WGS to improve our understanding of disease etiology.

Every gene everywhere all at once: High-precision measurement of 3D chromosome architecture with single-cell Hi-C

The three-dimensional (3D) structure of chromosomes influences essential biological processes such as gene expression, genome replication, and DNA damage repair and has been implicated in many developmental and degenerative diseases. In the past two centuries, two complementary genres of technology-microscopy, such as fluorescence in situ hybridization (FISH), and biochemistry, such as chromosome conformation capture (3C or Hi-C)-have revealed general principles of chromosome folding in the cell nucleus. However, the extraordinary complexity and cell-to-cell variability of the chromosome structure necessitate new tools with genome-wide coverage and single-cell precision. In the past decade, single-cell Hi-C emerges as a new approach that builds upon yet conceptually differs from bulk Hi-C assays. Instead of measuring population-averaged statistical properties of chromosome folding, single-cell Hi-C works as a proximity-based "biochemical microscope" that measures actual 3D structures of individual genomes, revealing features hidden in bulk Hi-C such as radial organization, multi-way interactions, and chromosome intermingling. Single-cell Hi-C has been used to study highly dynamic processes such as the cell cycle, cell-type-specific chromosome architecture ("structure types"), and structure-expression interplay, deepening our understanding of DNA organization and function.

Identifying enhancer properties associated with genetic risk for complex traits using regulome-wide association studies

Genetic risk for complex traits is strongly enriched in non-coding genomic regions involved in gene regulation, especially enhancers. However, we lack adequate tools to connect the characteristics of these disruptions to genetic risk. Here, we propose RWAS (Regulome Wide Association Study), a new application of the MAGMA software package to identify the characteristics of enhancers that contribute to genetic risk for disease. RWAS involves three steps: (i) assign genotyped SNPs to cell type- or tissue-specific regulatory features (e.g., enhancers); (ii) test associations of each regulatory feature with a trait of interest for which genome-wide association study (GWAS) summary statistics are available; (iii) perform enhancer-set enrichment analyses to identify quantitative or categorical features of regulatory elements that are associated with the trait. These steps are implemented as a novel application of MAGMA, a tool originally developed for gene-based GWAS analyses. Applying RWAS to interrogate genetic risk for schizophrenia, we discovered a class of risk-associated AT-rich enhancers that are active in the developing brain and harbor binding sites for multiple transcription factors with neurodevelopmental functions. RWAS utilizes open-source software, and we provide a comprehensive collection of annotations for tissue-specific enhancer locations and features, including their evolutionary conservation, AT content, and co-localization with binding sites for hundreds of TFs. RWAS will enable researchers to characterize properties of regulatory elements associated with any trait of interest for which GWAS summary statistics are available.

WhichTF is functionally important in your open chromatin data?

We present WhichTF, a computational method to identify functionally important transcription factors (TFs) from chromatin accessibility measurements. To rank TFs, WhichTF applies an ontology-guided functional approach to compute novel enrichment by integrating accessibility measurements, high-confidence pre-computed conservation-aware TF binding sites, and putative gene-regulatory models. Comparison with prior sheer abundance-based methods reveals the unique ability of WhichTF to identify context-specific TFs with functional relevance, including NF-κB family members in lymphocytes and GATA factors in cardiac cells. To distinguish the transcriptional regulatory landscape in closely related samples, we apply differential analysis and demonstrate its utility in lymphocyte, mesoderm developmental, and disease cells. We find suggestive, under-characterized TFs, such as RUNX3 in mesoderm development and GLI1 in systemic lupus erythematosus. We also find TFs known for stress response, suggesting routine experimental caveats that warrant careful consideration. WhichTF yields biological insight into known and novel molecular mechanisms of TF-mediated transcriptional regulation in diverse contexts, including human and mouse cell types, cell fate trajectories, and disease-associated cells.

Transcription factors (TFs), a class of DNA binding proteins, regulate tissue- and cell-type-specific expression of genes. Identifying the critical TFs in a given cellular context leads to investigating molecular regulatory mechanisms in development, differentiation, and disease. Because there are more than 1,500 human TFs, experimental measurements of genome-wide occupancy across all TFs have been challenging. While computational approaches play pivotal roles, most existing methods rely on statistical enrichment, focusing either on sequence motif similarity recognized by TFs or the similarity of the genomic region of interest with the previously characterized TF occupancy profile. Here we propose WhichTF as an alternative, incorporating curated biomedical knowledge from ontology and integrating it with the high-confidence prediction of conserved TF binding sites in user-provided genomic regions of interest. We develop a new WhichTF score to rank TFs and demonstrate its applicability across human and mouse cell types, cellular differentiation trajectories, and disease-associated cells.

Shared regulation and functional relevance of local gene co-expression revealed by single cell analysis

Most human genes are co-expressed with a nearby gene. Previous studies have revealed this local gene co-expression to be widespread across chromosomes and across dozens of tissues. Yet, so far these studies used bulk RNA-seq, averaging gene expression measurements across millions of cells, thus being unclear if this co-expression stems from transcription events in single cells. Here, we leverage single cell datasets in >85 individuals to identify gene co-expression across cells, unbiased by cell-type heterogeneity and benefiting from the co-occurrence of transcription events in single cells. We discover >3800 co-expressed gene pairs in two human cell types, induced pluripotent stem cells (iPSCs) and lymphoblastoid cell lines (LCLs) and (i) compare single cell to bulk RNA-seq in identifying local gene co-expression, (ii) show that many co-expressed genes – but not the majority – are composed of functionally related genes and (iii) using proteomics data, provide evidence that their co-expression is maintained up to the protein level. Finally, using single cell RNA-sequencing (scRNA-seq) and single cell ATAC-sequencing (scATAC-seq) data for the same single cells, we identify gene-enhancer associations and reveal that >95% of co-expressed gene pairs share regulatory elements. These results elucidate the potential reasons for co-expression in single cell gene regulatory networks and warrant a deeper study of shared regulatory elements, in view of explaining disease comorbidity due to affecting several genes. Our in-depth view of local gene co-expression and regulatory element co-activity advances our understanding of the shared regulatory architecture between genes.

Using single-cell data from cell lines, the co-expression of genes and co-activity of regulatory elements is analyzed, providing insight into shared architecture and regulation between genes.

Identification of TMEM129, encoding a ubiquitin-protein ligase, as an effector gene of osteoarthritis genetic risk

BACKGROUND

Osteoarthritis is highly heritable and genome-wide studies have identified single nucleotide polymorphisms (SNPs) associated with the disease. One such locus is marked by SNP rs11732213 (T > C). Genotype at rs11732213 correlates with the methylation levels of nearby CpG dinucleotides (CpGs), forming a methylation quantitative trait locus (mQTL). This study investigated the regulatory activity of the CpGs to identify a target gene of the locus.

METHODS

Nucleic acids were extracted from the articular cartilage of osteoarthritis patients. Samples were genotyped, and DNA methylation was quantified by pyrosequencing at 14 CpGs within a 259-bp interval. CpGs were tested for enhancer effects in immortalised chondrocytes using a reporter gene assay. DNA methylation at the locus was altered using targeted epigenome editing, with the impact on gene expression determined using quantitative polymerase chain reaction.

RESULTS

rs11732213 genotype correlated with DNA methylation at nine CpGs, which formed a differentially methylated region (DMR), with the osteoarthritis risk allele T corresponding to reduced levels of methylation. The DMR acted as an enhancer and demethylation of the CpGs altered expression of TMEM129. Allelic imbalance in TMEM129 expression was identified in cartilage, with under-expression of the risk allele.

CONCLUSIONS

TMEM129 is a target of osteoarthritis genetic risk at this locus. Genotype at rs11732213 impacts DNA methylation at the enhancer, which, in turn, modulates TMEM129 expression. TMEM129 encodes an enzyme involved in protein degradation within the endoplasmic reticulum, a process previously implicated in osteoarthritis. TMEM129 is a compelling osteoarthritis susceptibility target.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Comparative chromatin accessibility upon BDNF stimulation delineates neuronal regulatory elements

Neuronal stimulation induced by the brain-derived neurotrophic factor (BDNF) triggers gene expression, which is crucial for neuronal survival, differentiation, synaptic plasticity, memory formation, and neurocognitive health. However, its role in chromatin regulation is unclear. Here, using temporal profiling of chromatin accessibility and transcription in mouse primary cortical neurons upon either BDNF stimulation or depolarization (KCl), we identify features that define BDNF-specific chromatin-to-gene expression programs. Enhancer activation is an early event in the regulatory control of BDNF-treated neurons, where the bZIP motif-binding Fos protein pioneered chromatin opening and cooperated with co-regulatory transcription factors (Homeobox, EGRs, and CTCF) to induce transcription. Deleting cis-regulatory sequences affect BDNF-mediated Arc expression, a regulator of synaptic plasticity. BDNF-induced accessible regions are linked to preferential exon usage by neurodevelopmental disorder-related genes and the heritability of neuronal complex traits, which were validated in human iPSC-derived neurons. Thus, we provide a comprehensive view of BDNF-mediated genome regulatory features using comparative genomic approaches to dissect mammalian neuronal stimulation.

Analysis of somatic mutations in 131 human brains reveals aging-associated hypermutability

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

Comparison of chromatin accessibility landscapes during early development of prefrontal cortex between rhesus macaque and human

Epigenetic information regulates gene expression and development. However, our understanding of the evolution of epigenetic regulation on brain development in primates is limited. Here, we compared chromatin accessibility landscapes and transcriptomes during fetal prefrontal cortex (PFC) development between rhesus macaques and humans. A total of 304,761 divergent DNase I-hypersensitive sites (DHSs) are identified between rhesus macaques and humans, although many of these sites share conserved DNA sequences. Interestingly, most of the cis-elements linked to orthologous genes with dynamic expression are divergent DHSs. Orthologous genes expressed at earlier stages tend to have conserved cis-elements, whereas orthologous genes specifically expressed at later stages seldom have conserved cis-elements. These genes are enriched in synapse organization, learning and memory. Notably, DHSs in the PFC at early stages are linked to human educational attainment and cognitive performance. Collectively, the comparison of the chromatin epigenetic landscape between rhesus macaques and humans suggests a potential role for regulatory elements in the evolution of differences in cognitive ability between non-human primates and humans.

The evolution of epigenetic regulation of brain development in primates is not well understood. Here, the authors perform a comparative study of epigenetic dynamics of early prefrontal cortex development between human and rhesus macaque, finding divergent regulatory elements that may be related to cognitive capacity.

Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease

More than 800 million people suffer from kidney disease, yet the mechanism of kidney dysfunction is poorly understood. In the present study, we define the genetic association with kidney function in 1.5 million individuals and identify 878 (126 new) loci. We map the genotype effect on the methylome in 443 kidneys, transcriptome in 686 samples and single-cell open chromatin in 57,229 kidney cells. Heritability analysis reveals that methylation variation explains a larger fraction of heritability than gene expression. We present a multi-stage prioritization strategy and prioritize target genes for 87% of kidney function loci. We highlight key roles of proximal tubules and metabolism in kidney function regulation. Furthermore, the causal role of SLC47A1 in kidney disease is defined in mice with genetic loss of Slc47a1 and in human individuals carrying loss-of-function variants. Our findings emphasize the key role of bulk and single-cell epigenomic information in translating genome-wide association studies into identifying causal genes, cellular origins and mechanisms of complex traits.

Allelic imbalance of chromatin accessibility in cancer identifies candidate causal risk variants and their mechanisms

While many germline cancer risk variants have been identified through genome-wide association studies (GWAS), the mechanisms by which these variants operate remain largely unknown. Here we used 406 cancer ATAC-Seq samples across 23 cancer types to identify 7,262 germline allele-specific accessibility QTLs (as-aQTLs). Cancer as-aQTLs had stronger enrichment for cancer risk heritability (up to 145 fold) than any other functional annotation across seven cancer GWAS. Most cancer as-aQTLs directly altered transcription factor (TF) motifs and exhibited differential TF binding and gene expression in functional screens. To connect as-aQTLs to putative risk mechanisms, we introduced the regulome-wide associations study (RWAS). RWAS identified genetically associated accessible peaks at >70% of known breast and prostate loci and discovered new risk loci in all examined cancer types. Integrating as-aQTL discovery, motif analysis and RWAS identified candidate causal regulatory elements and their probable upstream regulators. Our work establishes cancer as-aQTLs and RWAS analysis as powerful tools to study the genetic architecture of cancer risk.

Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity

Disease-associated single-nucleotide polymorphisms (SNPs) generally do not implicate target genes, as most disease SNPs are regulatory. Many SNP-to-gene (S2G) linking strategies have been developed to link regulatory SNPs to the genes that they regulate in cis. Here, we developed a heritability-based framework for evaluating and combining different S2G strategies to optimize their informativeness for common disease risk. Our optimal combined S2G strategy (cS2G) included seven constituent S2G strategies and achieved a precision of 0.75 and a recall of 0.33, more than doubling the recall of any individual strategy. We applied cS2G to fine-mapping results for 49 UK Biobank diseases/traits to predict 5,095 causal SNP-gene-disease triplets (with S2G-derived functional interpretation) with high confidence. We further applied cS2G to provide an empirical assessment of disease omnigenicity; we determined that the top 1% of genes explained roughly half of the SNP heritability linked to all genes and that gene-level architectures vary with variant allele frequency.

Integrating convolution and self-attention improves language model of human genome for interpreting non-coding regions at base-resolution

Interpretation of non-coding genome remains an unsolved challenge in human genetics due to impracticality of exhaustively annotating biochemically active elements in all conditions. Deep learning based computational approaches emerge recently to help interpret non-coding regions. Here, we present LOGO (Language of Genome), a self-attention based contextualized pre-trained language model containing only two self-attention layers with 1 million parameters as a substantially light architecture that applies self-supervision techniques to learn bidirectional representations of the unlabelled human reference genome. LOGO is then fine-tuned for sequence labelling task, and further extended to variant prioritization task via a special input encoding scheme of alternative alleles followed by adding a convolutional module. Experiments show that LOGO achieves 15% absolute improvement for promoter identification and up to 4.5% absolute improvement for enhancer-promoter interaction prediction. LOGO exhibits state-of-the-art multi-task predictive power on thousands of chromatin features with only 3% parameterization benchmarking against the fully supervised model, DeepSEA and 1% parameterization against a recent BERT-based DNA language model. For allelic-effect prediction, locality introduced by one dimensional convolution shows improved sensitivity and specificity for prioritizing non-coding variants associated with human diseases. In addition, we apply LOGO to interpret type 2 diabetes (T2D) GWAS signals and infer underlying regulatory mechanisms. We make a conceptual analogy between natural language and human genome and demonstrate LOGO is an accurate, fast, scalable, and robust framework to interpret non-coding regions for global sequence labeling as well as for variant prioritization at base-resolution.

Sequence-based modeling of three-dimensional genome architecture from kilobase to chromosome scale

To learn how genomic sequence influences multiscale three-dimensional (3D) genome architecture, this manuscript presents a sequence-based deep-learning approach, Orca, that predicts directly from sequence the 3D genome architecture from kilobase to whole-chromosome scale. Orca captures the sequence dependencies of structures including chromatin compartments and topologically associating domains, as well as diverse types of interactions from CTCF-mediated to enhancer-promoter interactions and Polycomb-mediated interactions with cell-type specificity. Orca enables various applications including predicting structural variant effects on multiscale genome organization and it recapitulated effects of experimentally studied variants at varying sizes (300 bp to 90 Mb). Moreover, Orca enables in silico virtual screens to probe the sequence basis of 3D genome organization at different scales. At the submegabase scale, it predicted specific transcription factor motifs underlying cell-type-specific genome interactions. At the compartment scale, virtual screens of sequence activities suggest a model for the sequence basis of chromatin compartments with a prominent role of transcription start sites.

Glioblastoma Embryonic-like Stem Cells Exhibit Immune-Evasive Phenotype

Simple Summary

Most glioblastoma (GBM) patients relapse after an initial response to treatment. These aggressive traits are often associated with the presence of glioma stem cells (GSCs) within the tumor bulk, which are thought to participate in GBM therapy resistance. Given GBM cellular heterogeneity, we hypothesized that GSCs might also display cellular hierarchies associated with different degrees of stemness. Based on single-cell RNAseq data from GBM patients, we identified a subpopulation of GSCs, named core-GSCs (c-GSCs), with a similar profile to embryonic stem cells and downregulation of immune-associated pathways. In addition, we developed an in vitro induced c-GSC (ic-GSC) model resembling their tumor counterpart. The characterization of immune-privileged c-GSCs provides a valuable resource to study immune evasion mechanisms in GBM and to identify potential unexplored targets to improve immunotherapy treatments.

Abstract

Background: Glioma stem cells (GSCs) have self-renewal and tumor-initiating capacities involved in drug resistance and immune evasion mechanisms in glioblastoma (GBM). Methods: Core-GSCs (c-GSCs) were identified by selecting cells co-expressing high levels of embryonic stem cell (ESC) markers from a single-cell RNA-seq patient-derived GBM dataset (n = 28). Induced c-GSCs (ic-GSCs) were generated by reprogramming GBM-derived cells (GBM-DCs) using induced pluripotent stem cell (iPSC) technology. The characterization of ic-GSCs and GBM-DCs was conducted by immunostaining, transcriptomic, and DNA methylation (DNAm) analysis. Results: We identified a GSC population (4.22% ± 0.59) exhibiting concurrent high expression of ESC markers and downregulation of immune-associated pathways, named c-GSCs. In vitro ic-GSCs presented high expression of ESC markers and downregulation of antigen presentation HLA proteins. Transcriptomic analysis revealed a strong agreement of enriched biological pathways between tumor c-GSCs and in vitro ic-GSCs (κ = 0.71). Integration of our epigenomic profiling with 833 functional ENCODE epigenetic maps identifies increased DNA methylation on HLA genes’ regulatory regions associated with polycomb repressive marks in a stem-like phenotype. Conclusions: This study unravels glioblastoma immune-evasive mechanisms involving a c-GSC population. In addition, it provides a cellular model with paired gene expression, and DNA methylation maps to explore potential therapeutic complements for GBM immunotherapy.

Single-Cell Epigenomics Reveals Mechanisms of Cancer Progression

Cancer initiation is driven by the cooperation between genetic and epigenetic aberrations that disrupt gene regulatory programs critical to maintaining specialized cellular functions. After initiation, cells acquire additional genetic and epigenetic alterations influenced by tumor-intrinsic and -extrinsic mechanisms, which increase intratumoral heterogeneity, reshape the cell's underlying gene regulatory networks and promote cancer evolution. Furthermore, environmental or therapeutic insults drive the selection of heterogeneous cell states, with implications for cancer initiation, maintenance, and drug resistance. The advancement of single-cell genomics has begun to uncover the full repertoire of chromatin and gene expression states (cell states) that exist within individual tumors. These single-cell analyses suggest that cells diversify in their regulatory states upon transformation by co-opting damage-induced and nonlineage regulatory programs that can lead to epigenomic plasticity. Here, we review these recent studies related to regulatory state changes in cancer progression and highlight the growing single-cell epigenomics toolkit poised to address unresolved questions in the field.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

MNX1-HNF1B Axis Is Indispensable for Intraductal Papillary Mucinous Neoplasm Lineages

BACKGROUND & AIMS

Chromatin architecture governs cell lineages by regulating the specific gene expression; however, its role in the diversity of cancer development remains unknown. Among pancreatic cancers, pancreatic ductal adenocarcinoma (PDAC) and intraductal papillary mucinous neoplasms (IPMN) with an associated invasive carcinoma (IPMNinv) arise from 2 distinct precursors, and their fundamental differences remain obscure. Here, we aimed to assess the difference of chromatin architecture regulating the transcriptional signatures or biological features in pancreatic cancers.

METHODS

We established 28 human organoids from distinct subtypes of pancreatic tumors, including IPMN, IPMNinv, and PDAC. We performed exome sequencing (seq), RNA-seq, assay for transposase-accessible chromatin-seq, chromatin immunoprecipitation-seq, high-throughput chromosome conformation capture, and phenotypic analyses with short hairpin RNA or clustered regularly interspaced short palindromic repeats interference.

RESULTS

Established organoids successfully reproduced the histology of primary tumors. IPMN and IPMNinv organoids harbored GNAS, RNF43, or KLF4 mutations and showed the distinct expression profiles compared with PDAC. Chromatin accessibility profiles revealed the gain of stomach-specific open regions in IPMN and the pattern of diverse gastrointestinal tissues in IPMNinv. In contrast, PDAC presented an impressive loss of accessible regions compared with normal pancreatic ducts. Transcription factor footprint analysis and functional assays identified that MNX1 and HNF1B were biologically indispensable for IPMN lineages. The upregulation of MNX1 was specifically marked in the human IPMN lineage tissues. The MNX1-HNF1B axis governed a set of genes, including MYC, SOX9, and OLFM4, which are known to be essential for gastrointestinal stem cells. High-throughput chromosome conformation capture analysis suggested the HNF1B target genes to be 3-dimensionally connected in the genome of IPMNinv.

CONCLUSIONS

Our organoid analyses identified the MNX1-HNF1B axis to be biologically significant in IPMN lineages.

Synthesis and Profiling of Highly Selective Inhibitors of Methyltransferase DOT1L Based on Carbocyclic C-Nucleosides

Histone methyltransferase DOT1L is an attractive therapeutic target for the treatment of hematological malignancies. Here, we report the design, synthesis, and profiling of new DOT1L inhibitors based on nonroutine carbocyclic C-nucleoside scaffolds. The experimentally observed SAR was found to be nontrivial as seemingly minor changes of individual substituents resulted in significant changes in the affinity to DOT1L. Molecular modeling suggested that these trends could be related to significant conformational changes of the protein upon interaction with the inhibitors. The compounds 22 and (-)-53 (MU1656), carbocyclic C-nucleoside analogues of the natural nucleoside derivative EPZ004777, and the clinical candidate EPZ5676 (pinometostat) potently and selectively inhibit DOT1L in vitro as well as in the cell. The most potent compound MU1656 was found to be more metabolically stable and significantly less toxic in vivo than pinometostat itself.