GREGOR: evaluating global enrichment of trait-associated variants in epigenomic features using a systematic, data-driven approach

A compendium of genetic variations associated with promoter usage across 49 human tissues

Promoters play a crucial role in regulating gene transcription. However, our understanding of how genetic variants influence alternative promoter selection is still incomplete. In this study, we implement a framework to identify genetic variants that affect the relative usage of alternative promoters, known as promoter usage quantitative trait loci (puQTLs). By constructing an atlas of human puQTLs across 49 different tissues from 838 individuals, we have identified approximately 76,856 independent loci associated with promoter usage, encompassing 602,009 genetic variants. Our study demonstrates that puQTLs represent a distinct type of molecular quantitative trait loci, effectively uncovering regulatory targets and patterns. Furthermore, puQTLs are regulating in a tissue-specific manner and are enriched with binding sites of epigenetic marks and transcription factors, especially those involved in chromatin architecture formation. Notably, we have also found that puQTLs colocalize with complex traits or diseases and contribute to their heritability. Collectively, our findings underscore the significant role of puQTLs in elucidating the molecular mechanisms underlying tissue development and complex diseases.

Cross-ancestry analysis of brain QTLs enhances interpretation of schizophrenia genome-wide association studies

Research on brain expression quantitative trait loci (eQTLs) has illuminated the genetic underpinnings of schizophrenia (SCZ). Yet most of these studies have been centered on European populations, leading to a constrained understanding of population diversities and disease risks. To address this gap, we examined genotype and RNA-seq data from African Americans (AA, n = 158), Europeans (EUR, n = 408), and East Asians (EAS, n = 217). When comparing eQTLs between EUR and non-EUR populations, we observed concordant patterns of genetic regulatory effect, particularly in terms of the effect sizes of the eQTLs. However, 343,737 cis-eQTLs linked to 1,276 genes and 198,769 SNPs were found to be specific to non-EUR populations. Over 90% of observed population differences in eQTLs could be traced back to differences in allele frequency. Furthermore, 35% of these eQTLs were notably rare in the EUR population. Integrating brain eQTLs with SCZ signals from diverse populations, we observed a higher disease heritability enrichment of brain eQTLs in matched populations compared to mismatched ones. Prioritization analysis identified five risk genes (SFXN2, VPS37B, DENR, FTCDNL1, and NT5DC2) and three potential regulatory variants in known risk genes (CNNM2, MTRFR, and MPHOSPH9) that were missed in the EUR dataset. Our findings underscore that increasing genetic ancestral diversity is more efficient for power improvement than merely increasing the sample size within single-ancestry eQTLs datasets. Such a strategy will not only improve our understanding of the biological underpinnings of population structures but also pave the way for the identification of risk genes in SCZ.

Widespread transposable element dysregulation in human aging brains with Alzheimer's disease

INTRODUCTION

Transposable element (TE) dysregulation is associated with neuroinflammation in Alzheimer's disease (AD) brains. Yet, TE quantitative trait loci (teQTL) have not been well characterized in human aged brains with AD.

METHODS

We leveraged large-scale bulk and single-cell RNA sequencing, whole-genome sequencing (WGS), and xQTL from three human AD brain biobanks to characterize TE expression dysregulation and experimentally validate AD-associated TEs using CRISPR interference (CRISPRi) assays in human induced pluripotent stem cell (iPSC)-derived neurons.

RESULTS

We identified 26,188 genome-wide significant TE expression QTLs (teQTLs) in human aged brains. Subsequent colocalization analysis of teQTLs with AD genetic loci identified AD-associated teQTLs and linked locus TEs. Using CRISPRi assays, we pinpointed a neuron-specific suppressive role of the activated short interspersed nuclear element (SINE; chr11:47608036-47608220) on expression of C1QTNF4 via reducing neuroinflammation in human iPSC-derived neurons.

DISCUSSION

We identified widespread TE dysregulation in human AD brains and teQTLs offer a complementary analytic approach to identify likely AD risk genes.

HIGHLIGHTS

Widespread transposable element (TE) dysregulations are observed in human aging brains with degrees of neuropathology, apolipoprotein E (APOE) genotypes, and neuroinflammation in Alzheimer's disease (AD). A catalog of TE quantitative trait loci (teQTLs) in human aging brains was created using matched RNA sequencing and whole-genome sequencing data. CRISPR interference assays reveal that an upregulated intergenic TE from the MIR family (chr11: 47608036-47608220) suppresses expression of its nearest anti-inflammatory gene C1QTNF4 in human induced pluripotent stem cell-derived neurons.

Genome-wide discovery and integrative genomic characterization of insulin resistance loci using serum triglycerides to HDL-cholesterol ratio as a proxy

Insulin resistance causes multiple epidemic metabolic diseases, including type 2 diabetes, cardiovascular disease, and fatty liver, but is not routinely measured in epidemiological studies. To discover novel insulin resistance genes in the general population, we conducted genome-wide association studies in 382,129 individuals for triglyceride to HDL-cholesterol ratio (TG/HDL), a surrogate marker of insulin resistance calculable from commonly measured serum lipid profiles. We identified 251 independent loci, of which 62 were more strongly associated with TG/HDL compared to TG or HDL alone, suggesting them as insulin resistance loci. Candidate causal genes at these loci were prioritized by fine mapping with directions-of-effect and tissue specificity annotated through analysis of protein coding and expression quantitative trait variation. Directions-of-effect were corroborated in an independent cohort of individuals with directly measured insulin resistance. We highlight two phospholipase encoding genes, PLA2G12A and PLA2G6, which liberate arachidonic acid and improve insulin sensitivity, and VGLL3, a transcriptional co-factor that increases insulin resistance partially through enhanced adiposity. Finally, we implicate the anti-apoptotic gene TNFAIP8 as a sex-dimorphic insulin resistance factor, which acts by increasing visceral adiposity, specifically in females. In summary, our study identifies several candidate modulators of insulin resistance that have the potential to serve as biomarkers and pharmacological targets.

Untangling the genetics of beta cell dysfunction and death in type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D) is a complex multi-system disease which arises from both environmental and genetic factors, resulting in the destruction of insulin-producing pancreatic beta cells. Over the past two decades, human genetic studies have provided new insight into the etiology of T1D, including an appreciation for the role of beta cells in their own demise.

SCOPE OF REVIEW

Here, we outline models supported by human genetic data for the role of beta cell dysfunction and death in T1D. We highlight the importance of strong evidence linking T1D genetic associations to bona fide candidate genes for mechanistic and therapeutic consideration. To guide rigorous interpretation of genetic associations, we describe molecular profiling approaches, genomic resources, and disease models that may be used to construct variant-to-gene links and to investigate candidate genes and their role in T1D.

MAJOR CONCLUSIONS

We profile advances in understanding the genetic causes of beta cell dysfunction and death at individual T1D risk loci. We discuss how genetic risk prediction models can be used to address disease heterogeneity. Further, we present areas where investment will be critical for the future use of genetics to address open questions in the development of new treatment and prevention strategies for T1D.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

From GWASs toward Mechanistic Understanding with Case Studies in Dermatogenetics

Many human skin diseases result from the complex interplay of genetic and environmental mechanisms that are largely unknown. GWASs have yielded insight into the genetic aspect of complex disease by highlighting regions of the genome or specific genetic variants associated with disease. Leveraging this information to identify causal genes and cell types will provide insight into fundamental biology, inform diagnostics, and aid drug discovery. However, the etiological mechanisms from genetic variant to disease are still unestablished in most cases. There now exists an unprecedented wealth of data and computational methods for variant interpretation in a functional context. It can be challenging to decide where to start owing to a lack of consensus on the best way to identify causal genetic mechanisms. This article highlights 3 key aspects of genetic variant interpretation: prioritizing causal genes, cell types, and pathways. We provide a practical overview of the main methods and datasets, giving examples from recent atopic dermatitis studies to provide a blueprint for variant interpretation. A collection of resources, including brief description and links to the packages and web tools, is provided for researchers looking to start in silico follow-up genetic analysis of associated genetic variants.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

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

BACKGROUND

As 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.

RESULTS

In 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 .

CONCLUSIONS

GaiaAssociation 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.

Sex-specific genetic architecture of blood pressure

The genetic and genomic basis of sex differences in blood pressure (BP) traits remain unstudied at scale. Here, we conducted sex-stratified and combined-sex genome-wide association studies of BP traits using the UK Biobank resource, identifying 1,346 previously reported and 29 new BP trait-associated loci. Among associated loci, 412 were female-specific (Pfemale ≤ 5 × 10-8; Pmale > 5 × 10-8) and 142 were male-specific (Pmale ≤ 5 × 10-8; Pfemale > 5 × 10-8); these sex-specific loci were enriched for hormone-related transcription factors, in particular, estrogen receptor 1. Analyses of gene-by-sex interactions and sexually dimorphic effects identified four genomic regions, showing female-specific associations with diastolic BP or pulse pressure, including the chromosome 13q34-COL4A1/COL4A2 locus. Notably, female-specific pulse pressure-associated loci exhibited enriched acetylated histone H3 Lys27 modifications in arterial tissues and a female-specific association with fibromuscular dysplasia, a female-biased vascular disease; colocalization signals included Chr13q34: COL4A1/COL4A2, Chr9p21: CDKN2B-AS1 and Chr4q32.1: MAP9 regions. Sex-specific and sex-biased polygenic associations of BP traits were associated with multiple cardiovascular traits. These findings suggest potentially clinically significant and BP sex-specific pleiotropic effects on cardiovascular diseases.

RicePilaf: a post-GWAS/QTL dashboard to integrate pangenomic, coexpression, regulatory, epigenomic, ontology, pathway, and text-mining information to provide functional insights into rice QTLs and GWAS loci

BACKGROUND

As the number of genome-wide association study (GWAS) and quantitative trait locus (QTL) mappings in rice continues to grow, so does the already long list of genomic loci associated with important agronomic traits. Typically, loci implicated by GWAS/QTL analysis contain tens to hundreds to thousands of single-nucleotide polmorphisms (SNPs)/genes, not all of which are causal and many of which are in noncoding regions. Unraveling the biological mechanisms that tie the GWAS regions and QTLs to the trait of interest is challenging, especially since it requires collating functional genomics information about the loci from multiple, disparate data sources.

RESULTS

We present RicePilaf, a web app for post-GWAS/QTL analysis, that performs a slew of novel bioinformatics analyses to cross-reference GWAS results and QTL mappings with a host of publicly available rice databases. In particular, it integrates (i) pangenomic information from high-quality genome builds of multiple rice varieties, (ii) coexpression information from genome-scale coexpression networks, (iii) ontology and pathway information, (iv) regulatory information from rice transcription factor databases, (v) epigenomic information from multiple high-throughput epigenetic experiments, and (vi) text-mining information extracted from scientific abstracts linking genes and traits. We demonstrate the utility of RicePilaf by applying it to analyze GWAS peaks of preharvest sprouting and genes underlying yield-under-drought QTLs.

CONCLUSIONS

RicePilaf enables rice scientists and breeders to shed functional light on their GWAS regions and QTLs, and it provides them with a means to prioritize SNPs/genes for further experiments. The source code, a Docker image, and a demo version of RicePilaf are publicly available at https://github.com/bioinfodlsu/rice-pilaf.

Integrative splicing-quantitative-trait-locus analysis reveals risk loci for non-small-cell lung cancer

Splicing quantitative trait loci (sQTLs) have been demonstrated to contribute to disease etiology by affecting alternative splicing. However, the role of sQTLs in the development of non-small-cell lung cancer (NSCLC) remains unknown. Thus, we performed a genome-wide sQTL study to identify genetic variants that affect alternative splicing in lung tissues from 116 individuals of Chinese ancestry, which resulted in the identification of 1,385 sQTL-harboring genes (sGenes) containing 378,210 significant variant-intron pairs. A comprehensive characterization of these sQTLs showed that they were enriched in actively transcribed regions, genetic regulatory elements, and splicing-factor-binding sites. Moreover, sQTLs were largely distinct from expression quantitative trait loci (eQTLs) and showed significant enrichment in potential risk loci of NSCLC. We also integrated sQTLs into NSCLC GWAS datasets (13,327 affected individuals and 13,328 control individuals) by using splice-transcriptome-wide association study (spTWAS) and identified alternative splicing events in 19 genes that were significantly associated with NSCLC risk. By using functional annotation and experiments, we confirmed an sQTL variant, rs35861926, that reduced the risk of lung adenocarcinoma (rs35861926-T, OR = 0.88, 95% confidence interval [CI]: 0.82-0.93, p = 1.87 × 10-5) by promoting FARP1 exon 20 skipping to downregulate the expression level of the long transcript FARP1-011. Transcript FARP1-011 promoted the migration and proliferation of lung adenocarcinoma cells. Overall, our study provided informative lung sQTL resources and insights into the molecular mechanisms linking sQTL variants to NSCLC risk.

Androgen receptor binding sites enabling genetic prediction of mortality due to prostate cancer in cancer-free subjects

Prostate cancer (PrCa) is the second most common cancer worldwide in males. While strongly warranted, the prediction of mortality risk due to PrCa, especially before its development, is challenging. Here, we address this issue by maximizing the statistical power of genetic data with multi-ancestry meta-analysis and focusing on binding sites of the androgen receptor (AR), which has a critical role in PrCa. Taking advantage of large Japanese samples ever, a multi-ancestry meta-analysis comprising more than 300,000 subjects in total identifies 9 unreported loci including ZFHX3, a tumor suppressor gene, and successfully narrows down the statistically finemapped variants compared to European-only studies, and these variants strongly enrich in AR binding sites. A polygenic risk scores (PRS) analysis restricting to statistically finemapped variants in AR binding sites shows among cancer-free subjects, individuals with a PRS in the top 10% have a strongly higher risk of the future death of PrCa (HR: 5.57, P = 4.2 × 10-10). Our findings demonstrate the potential utility of leveraging large-scale genetic data and advanced analytical methods in predicting the mortality of PrCa.

Leveraging type 1 diabetes human genetic and genomic data in the T1D knowledge portal

To address the challenge of translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight, we have developed the T1D Knowledge Portal (T1DKP), an open-access resource for hypothesis development and target discovery in T1D.

Interpreting non-coding disease-associated human variants using single-cell epigenomics

Genome-wide association studies (GWAS) have linked hundreds of thousands of sequence variants in the human genome to common traits and diseases. However, translating this knowledge into a mechanistic understanding of disease-relevant biology remains challenging, largely because such variants are predominantly in non-protein-coding sequences that still lack functional annotation at cell-type resolution. Recent advances in single-cell epigenomics assays have enabled the generation of cell type-, subtype- and state-resolved maps of the epigenome in heterogeneous human tissues. These maps have facilitated cell type-specific annotation of candidate cis-regulatory elements and their gene targets in the human genome, enhancing our ability to interpret the genetic basis of common traits and diseases.

hGSuite HyperBrowser: A web-based toolkit for hierarchical metadata-informed analysis of genomic tracks

Many high-throughput sequencing datasets can be represented as objects with coordinates along a reference genome. Currently, biological investigations often involve a large number of such datasets, for example representing different cell types or epigenetic factors. Drawing overall conclusions from a large collection of results for individual datasets may be challenging and time-consuming. Meaningful interpretation often requires the results to be aggregated according to metadata that represents biological characteristics of interest. In this light, we here propose the hierarchical Genomic Suite HyperBrowser (hGSuite), an open-source extension to the GSuite HyperBrowser platform, which aims to provide a means for extracting key results from an aggregated collection of high-throughput DNA sequencing data. The hGSuite utilizes a metadata-informed data cube to calculate various statistics across the multiple dimensions of the datasets. With this work, we show that the hGSuite and its associated data cube methodology offers a quick and accessible way for exploratory analysis of large genomic datasets. The web-based toolkit named hGsuite Hyperbrowser is available at https://hyperbrowser.uio.no/hgsuite under a GPLv3 license.

Mapping disease regulatory circuits at cell-type resolution from single-cell multiomics data

Resolving chromatin-remodeling-linked gene expression changes at cell-type resolution is important for understanding disease states. Here we describe MAGICAL (Multiome Accessibility Gene Integration Calling and Looping), a hierarchical Bayesian approach that leverages paired single-cell RNA sequencing and single-cell transposase-accessible chromatin sequencing from different conditions to map disease-associated transcription factors, chromatin sites, and genes as regulatory circuits. By simultaneously modeling signal variation across cells and conditions in both omics data types, MAGICAL achieved high accuracy on circuit inference. We applied MAGICAL to study Staphylococcus aureus sepsis from peripheral blood mononuclear single-cell data that we generated from subjects with bloodstream infection and uninfected controls. MAGICAL identified sepsis-associated regulatory circuits predominantly in CD14 monocytes, known to be activated by bacterial sepsis. We addressed the challenging problem of distinguishing host regulatory circuit responses to methicillin-resistant and methicillin-susceptible S. aureus infections. Although differential expression analysis failed to show predictive value, MAGICAL identified epigenetic circuit biomarkers that distinguished methicillin-resistant from methicillin-susceptible S. aureus infections.

Construction and evaluation of the functional polygenic risk score for gastric cancer in a prospective cohort of the European population

BACKGROUND

A polygenic risk score (PRS) derived from 112 single-nucleotide polymorphisms (SNPs) for gastric cancer has been reported in Chinese populations (PRS-112). However, its performance in other populations is unknown. A functional PRS (fPRS) using functional SNPs (fSNPs) may improve the generalizability of the PRS across populations with distinct ethnicities.

METHODS

We performed functional annotations on SNPs in strong linkage disequilibrium (LD) with the 112 previously reported SNPs to identify fSNPs that affect protein-coding or transcriptional regulation. Subsequently, we constructed an fPRS based on the fSNPs by using the LDpred2-infinitesimal model and then analyzed the performance of the PRS-112 and fPRS in the risk prediction of gastric cancer in 457,521 European participants of the UK Biobank cohort. Finally, the performance of the fPRS in combination with lifestyle factors were evaluated in predicting the risk of gastric cancer.

RESULTS

During 4,582,045 person-years of follow-up with a total of 623 incident gastric cancer cases, we found no significant association between the PRS-112 and gastric cancer risk in the European population (hazard ratio [HR] = 1.00 [95% confidence interval (CI) 0.93-1.09], P = 0.846). We identified 125 fSNPs, including seven deleterious protein-coding SNPs and 118 regulatory non-coding SNPs, and used them to constructed the fPRS-125. Our result showed that the fPRS-125 was significantly associated with gastric cancer risk (HR = 1.11 [95% CI, 1.03-1.20], P = 0.009). Compared to participants with a low fPRS-125 (bottom quintile), those with a high fPRS-125 (top quintile) had a higher risk of incident gastric cancer (HR = 1.43 [95% CI, 1.12-1.84], P = 0.005). Moreover, we observed that participants with both an unfavorable lifestyle and a high genetic risk had the highest risk of incident gastric cancer (HR = 4.99 [95% CI, 1.55-16.10], P = 0.007) compared to those with both a favorable lifestyle and a low genetic risk.

CONCLUSION

These results indicate that the fPRS-125 derived from fSNPs may act as an indicator to measure the genetic risk of gastric cancer in the European population.

Genome-wide association meta-analysis of spontaneous coronary artery dissection identifies risk variants and genes related to artery integrity and tissue-mediated coagulation

Spontaneous coronary artery dissection (SCAD) is an understudied cause of myocardial infarction primarily affecting women. It is not known to what extent SCAD is genetically distinct from other cardiovascular diseases, including atherosclerotic coronary artery disease (CAD). Here we present a genome-wide association meta-analysis (1,917 cases and 9,292 controls) identifying 16 risk loci for SCAD. Integrative functional annotations prioritized genes that are likely to be regulated in vascular smooth muscle cells and artery fibroblasts and implicated in extracellular matrix biology. One locus containing the tissue factor gene F3, which is involved in blood coagulation cascade initiation, appears to be specific for SCAD risk. Several associated variants have diametrically opposite associations with CAD, suggesting that shared biological processes contribute to both diseases, but through different mechanisms. We also infer a causal role for high blood pressure in SCAD. Our findings provide novel pathophysiological insights involving arterial integrity and tissue-mediated coagulation in SCAD and set the stage for future specific therapeutics and preventions.

FHL5 Controls Vascular Disease-Associated Gene Programs in Smooth Muscle Cells

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with common vascular diseases, such as coronary artery disease, myocardial infarction, and hypertension. However, the lack of mechanistic insights for many GWAS loci limits their translation into the clinic. Among these loci with unknown functions is UFL1-four-and-a-half LIM (LIN-11, Isl-1, MEC-3) domain 5 (FHL5; chr6q16.1), which reached genome-wide significance in a recent coronary artery disease/ myocardial infarction GWAS meta-analysis. UFL1-FHL5 is also associated with several vascular diseases, consistent with the widespread pleiotropy observed for GWAS loci.

METHODS

We apply a multimodal approach leveraging statistical fine-mapping, epigenomic profiling, and ex vivo analysis of human coronary artery tissues to implicate FHL5 as the top candidate causal gene. We unravel the molecular mechanisms of the cross-phenotype genetic associations through in vitro functional analyses and epigenomic profiling experiments in coronary artery smooth muscle cells.

RESULTS

We prioritized FHL5 as the top candidate causal gene at the UFL1-FHL5 locus through expression quantitative trait locus colocalization methods. FHL5 gene expression was enriched in the smooth muscle cells and pericyte population in human artery tissues with coexpression network analyses supporting a functional role in regulating smooth muscle cell contraction. Unexpectedly, under procalcifying conditions, FHL5 overexpression promoted vascular calcification and dysregulated processes related to extracellular matrix organization and calcium handling. Lastly, by mapping FHL5 binding sites and inferring FHL5 target gene function using artery tissue gene regulatory network analyses, we highlight regulatory interactions between FHL5 and downstream coronary artery disease/myocardial infarction loci, such as FOXL1 and FN1 that have roles in vascular remodeling.

CONCLUSIONS

Taken together, these studies provide mechanistic insights into the pleiotropic genetic associations of UFL1-FHL5. We show that FHL5 mediates vascular disease risk through transcriptional regulation of downstream vascular remodeling gene programs. These transacting mechanisms may explain a portion of the heritable risk for complex vascular diseases.

The Type 2 Diabetes Knowledge Portal: An open access genetic resource dedicated to type 2 diabetes and related traits

Associations between human genetic variation and clinical phenotypes have become a foundation of biomedical research. Most repositories of these data seek to be disease-agnostic and therefore lack disease-focused views. The Type 2 Diabetes Knowledge Portal (T2DKP) is a public resource of genetic datasets and genomic annotations dedicated to type 2 diabetes (T2D) and related traits. Here, we seek to make the T2DKP more accessible to prospective users and more useful to existing users. First, we evaluate the T2DKP's comprehensiveness by comparing its datasets with those of other repositories. Second, we describe how researchers unfamiliar with human genetic data can begin using and correctly interpreting them via the T2DKP. Third, we describe how existing users can extend their current workflows to use the full suite of tools offered by the T2DKP. We finally discuss the lessons offered by the T2DKP toward the goal of democratizing access to complex disease genetic results.

Genetic control of RNA editing in neurodegenerative disease

A-to-I RNA editing diversifies human transcriptome to confer its functional effects on the downstream genes or regulations, potentially involving in neurodegenerative pathogenesis. Its variabilities are attributed to multiple regulators, including the key factor of genetic variants. To comprehensively investigate the potentials of neurodegenerative disease-susceptibility variants from the view of A-to-I RNA editing, we analyzed matched genetic and transcriptomic data of 1596 samples across nine brain tissues and whole blood from two large consortiums, Accelerating Medicines Partnership-Alzheimer's Disease and Parkinson's Progression Markers Initiative. The large-scale and genome-wide identification of 95 198 RNA editing quantitative trait loci revealed the preferred genetic effects on adjacent editing events. Furthermore, to explore the underlying mechanisms of the genetic controls of A-to-I RNA editing, several top RNA-binding proteins were pointed out, such as EIF4A3, U2AF2, NOP58, FBL, NOP56 and DHX9, since their regulations on multiple RNA-editing events were probably interfered by these genetic variants. Moreover, these variants may also contribute to the variability of other molecular phenotypes associated with RNA editing, including the functions of 3 proteins, expressions of 277 genes and splicing of 449 events. All the analyses results shown in NeuroEdQTL (https://relab.xidian.edu.cn/NeuroEdQTL/) constituted a unique resource for the understanding of neurodegenerative pathogenesis from genotypes to phenotypes related to A-to-I RNA editing.

Translating non-coding genetic associations into a better understanding of immune-mediated disease

Genome-wide association studies have identified hundreds of genetic loci that are associated with immune-mediated diseases. Most disease-associated variants are non-coding, and a large proportion of these variants lie within enhancers. As a result, there is a pressing need to understand how common genetic variation might affect enhancer function and thereby contribute to immune-mediated (and other) diseases. In this Review, we first describe statistical and experimental methods to identify causal genetic variants that modulate gene expression, including statistical fine-mapping and massively parallel reporter assays. We then discuss approaches to characterise the mechanisms by which these variants modulate immune function, such as clustered regularly interspaced short palindromic repeats (CRISPR)-based screens. We highlight examples of studies that, by elucidating the effects of disease variants within enhancers, have provided important insights into immune function and uncovered key pathways of disease.

Genetic adaptation to pathogens and increased risk of inflammatory disorders in post-Neolithic Europe

Ancient genomics can directly detect human genetic adaptation to environmental cues. However, it remains unclear how pathogens have exerted selective pressures on human genome diversity across different epochs and affected present-day inflammatory disease risk. Here, we use an ancestry-aware approximate Bayesian computation framework to estimate the nature, strength, and time of onset of selection acting on 2,879 ancient and modern European genomes from the last 10,000 years. We found that the bulk of genetic adaptation occurred after the start of the Bronze Age, <4,500 years ago, and was enriched in genes relating to host-pathogen interactions. Furthermore, we detected directional selection acting on specific leukocytic lineages and experimentally demonstrated that the strongest negatively selected candidate variant in immunity genes, lipopolysaccharide-binding protein (LBP) D283G, is hypomorphic. Finally, our analyses suggest that the risk of inflammatory disorders has increased in post-Neolithic Europeans, possibly because of antagonistic pleiotropy following genetic adaptation to pathogens.

Leveraging type 1 diabetes human genetic and genomic data in the T1D Knowledge Portal

Translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight can reveal novel biology and therapeutic targets but remains a major challenge. We developed the T1D Knowledge Portal (T1DKP), a disease-specific resource of genetic and functional annotation data that enables users to develop hypotheses for T1D-based research and target discovery. The T1DKP can be used to query genes and genomic regions for genetic associations, identify epigenomic features, access results of bioinformatic analyses, and obtain expert-curated resources. The T1DKP is available at http://t1d.hugeamp.org .

Loci for insulin processing and secretion provide insight into type 2 diabetes risk

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

3D chromatin structure in chondrocytes identifies putative osteoarthritis risk genes

Genome-wide association studies have identified over 100 loci associated with osteoarthritis risk, but the majority of osteoarthritis risk variants are noncoding, making it difficult to identify the impacted genes for further study and therapeutic development. To address this need, we used a multiomic approach and genome editing to identify and functionally characterize potential osteoarthritis risk genes. Computational analysis of genome-wide association studies and ChIP-seq data revealed that chondrocyte regulatory loci are enriched for osteoarthritis risk variants. We constructed a chondrocyte-specific regulatory network by mapping 3D chromatin structure and active enhancers in human chondrocytes. We then intersected these data with our previously collected RNA-seq dataset of chondrocytes responding to fibronectin fragment, a known osteoarthritis trigger. Integration of the 3 genomic datasets with recently reported osteoarthritis genome-wide association study variants revealed a refined set of putative causal osteoarthritis variants and their potential target genes. One of the putative target genes identified was SOCS2, which was connected to a putative causal variant by a 170-kb loop and is differentially regulated in response to fibronectin fragment. CRISPR-Cas9-mediated deletion of SOCS2 in primary human chondrocytes from 3 independent donors led to heightened expression of inflammatory markers after fibronectin fragment treatment. These data suggest that SOCS2 plays a role in resolving inflammation in response to cartilage matrix damage and provides a possible mechanistic explanation for its influence on osteoarthritis risk. In total, we identified 56 unique putative osteoarthritis risk genes for further research and potential therapeutic development.

Integration of multimodal data in the developing tooth reveals candidate regulatory loci driving human odontogenic phenotypes

Human odontogenic aberrations such as abnormal tooth number and delayed tooth eruption can occur as a symptom of rare syndromes or, more commonly, as nonsyndromic phenotypes. These phenotypes can require extensive and expensive dental treatment, posing a significant burden. While many dental phenotypes are heritable, most nonsyndromic cases have not been linked to causal genes. We demonstrate the novel finding that common sequence variants associated with human odontogenic phenotypes are enriched in developmental craniofacial enhancers conserved between human and mouse. However, the bulk nature of these samples obscures if this finding is due to the tooth itself or the surrounding tissues. We therefore sought to identify enhancers specifically active in the tooth anlagen and quantify their contribution to the observed genetic enrichments. We systematically identified 22,001 conserved enhancers active in E13.5 mouse incisors using ChIP-seq and machine learning pipelines and demonstrated biologically relevant enrichments in putative target genes, transcription factor binding motifs, and in vivo activity. Multi-tissue comparisons of human and mouse enhancers revealed that these putative tooth enhancers had the strongest enrichment of odontogenic phenotype-associated variants, suggesting a role for dysregulation of tooth developmental enhancers in human dental phenotypes. The large number of these regions genome-wide necessitated prioritization of enhancer loci for future investigations. As enhancers modulate gene expression, we prioritized regions based on enhancers' putative target genes. We predicted these target genes and prioritized loci by integrating chromatin state, bulk gene expression and coexpression, GWAS variants, and cell type resolved gene expression to generate a prioritized list of putative odontogenic phenotype-driving loci active in the developing tooth. These genomic regions are of particular interest for downstream experiments determining the role of specific dental enhancer:gene pairs in odontogenesis.

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.

Population-level variation in enhancer expression identifies disease mechanisms in the human brain

Identification of risk variants for neuropsychiatric diseases within enhancers underscores the importance of understanding population-level variation in enhancer function in the human brain. Besides regulating tissue-specific and cell-type-specific transcription of target genes, enhancers themselves can be transcribed. By jointly analyzing large-scale cell-type-specific transcriptome and regulome data, we cataloged 30,795 neuronal and 23,265 non-neuronal candidate transcribed enhancers. Examination of the transcriptome in 1,382 brain samples identified robust expression of transcribed enhancers. We explored gene-enhancer coordination and found that enhancer-linked genes are strongly implicated in neuropsychiatric disease. We identified expression quantitative trait loci (eQTLs) for both genes and enhancers and found that enhancer eQTLs mediate a substantial fraction of neuropsychiatric trait heritability. Inclusion of enhancer eQTLs in transcriptome-wide association studies enhanced functional interpretation of disease loci. Overall, our study characterizes the gene-enhancer regulome and genetic mechanisms in the human cortex in both healthy and diseased states.

Genetic regulation of RNA splicing in human pancreatic islets

BACKGROUND

Non-coding genetic variants that influence gene transcription in pancreatic islets play a major role in the susceptibility to type 2 diabetes (T2D), and likely also contribute to type 1 diabetes (T1D) risk. For many loci, however, the mechanisms through which non-coding variants influence diabetes susceptibility are unknown.

RESULTS

We examine splicing QTLs (sQTLs) in pancreatic islets from 399 human donors and observe that common genetic variation has a widespread influence on the splicing of genes with established roles in islet biology and diabetes. In parallel, we profile expression QTLs (eQTLs) and use transcriptome-wide association as well as genetic co-localization studies to assign islet sQTLs or eQTLs to T2D and T1D susceptibility signals, many of which lack candidate effector genes. This analysis reveals biologically plausible mechanisms, including the association of T2D with an sQTL that creates a nonsense isoform in ERO1B, a regulator of ER-stress and proinsulin biosynthesis. The expanded list of T2D risk effector genes reveals overrepresented pathways, including regulators of G-protein-mediated cAMP production. The analysis of sQTLs also reveals candidate effector genes for T1D susceptibility such as DCLRE1B, a senescence regulator, and lncRNA MEG3.

CONCLUSIONS

These data expose widespread effects of common genetic variants on RNA splicing in pancreatic islets. The results support a role for splicing variation in diabetes susceptibility, and offer a new set of genetic targets with potential therapeutic benefit.

Temporal analysis of enhancers during mouse cerebellar development reveals dynamic and novel regulatory functions

We have identified active enhancers in the mouse cerebellum at embryonic and postnatal stages which provides a view of novel enhancers active during cerebellar development. The majority of cerebellar enhancers have dynamic activity between embryonic and postnatal development. Cerebellar enhancers were enriched for neural transcription factor binding sites with temporally specific expression. Putative gene targets displayed spatially restricted expression patterns, indicating cell-type specific expression regulation. Functional analysis of target genes indicated that enhancers regulate processes spanning several developmental epochs such as specification, differentiation and maturation. We use these analyses to discover one novel regulator and one novel marker of cerebellar development: Bhlhe22 and Pax3, respectively. We identified an enrichment of de novo mutations and variants associated with autism spectrum disorder in cerebellar enhancers. Furthermore, by comparing our data with relevant brain development ENCODE histone profiles and cerebellar single-cell datasets we have been able to generalize and expand on the presented analyses, respectively. We have made the results of our analyses available online in the Developing Mouse Cerebellum Enhancer Atlas (https://goldowitzlab.shinyapps.io/developing_mouse_cerebellum_enhancer_atlas/), where our dataset can be efficiently queried, curated and exported by the scientific community to facilitate future research efforts. Our study provides a valuable resource for studying the dynamics of gene expression regulation by enhancers in the developing cerebellum and delivers a rich dataset of novel gene-enhancer associations providing a basis for future in-depth studies in the cerebellum.

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.

Whole genome sequence association analysis of fasting glucose and fasting insulin levels in diverse cohorts from the NHLBI TOPMed program

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits.

Genome-wide association study identifies Sjögren's risk loci with functional implications in immune and glandular cells

Sjögren's disease is a complex autoimmune disease with twelve established susceptibility loci. This genome-wide association study (GWAS) identifies ten novel genome-wide significant (GWS) regions in Sjögren's cases of European ancestry: CD247, NAB1, PTTG1-MIR146A, PRDM1-ATG5, TNFAIP3, XKR6, MAPT-CRHR1, RPTOR-CHMP6-BAIAP6, TYK2, SYNGR1. Polygenic risk scores yield predictability (AUROC = 0.71) and relative risk of 12.08. Interrogation of bioinformatics databases refine the associations, define local regulatory networks of GWS SNPs from the 95% credible set, and expand the implicated gene list to >40. Many GWS SNPs are eQTLs for genes within topologically associated domains in immune cells and/or eQTLs in the main target tissue, salivary glands.

A multiancestry genome-wide association study of unexplained chronic ALT elevation as a proxy for nonalcoholic fatty liver disease with histological and radiological validation

Nonalcoholic fatty liver disease (NAFLD) is a growing cause of chronic liver disease. Using a proxy NAFLD definition of chronic elevation of alanine aminotransferase (cALT) levels without other liver diseases, we performed a multiancestry genome-wide association study (GWAS) in the Million Veteran Program (MVP) including 90,408 cALT cases and 128,187 controls. Seventy-seven loci exceeded genome-wide significance, including 25 without prior NAFLD or alanine aminotransferase associations, with one additional locus identified in European American-only and two in African American-only analyses (P < 5 × 10-8). External replication in histology-defined NAFLD cohorts (7,397 cases and 56,785 controls) or radiologic imaging cohorts (n = 44,289) replicated 17 single-nucleotide polymorphisms (SNPs) (P < 6.5 × 10-4), of which 9 were new (TRIB1, PPARG, MTTP, SERPINA1, FTO, IL1RN, COBLL1, APOH and IFI30). Pleiotropy analysis showed that 61 of 77 multiancestry and all 17 replicated SNPs were jointly associated with metabolic and/or inflammatory traits, revealing a complex model of genetic architecture. Our approach integrating cALT, histology and imaging reveals new insights into genetic liability to NAFLD.

Genetic variants underlying differences in facial morphology in East Asian and European populations

Facial morphology-a conspicuous feature of human appearance-is highly heritable. Previous studies on the genetic basis of facial morphology were performed mainly in European-ancestry cohorts (EUR). Applying a data-driven phenotyping and multivariate genome-wide scanning protocol to a large collection of three-dimensional facial images of individuals with East Asian ancestry (EAS), we identified 244 variants in 166 loci (62 new) associated with typical-range facial variation. A newly proposed polygenic shape analysis indicates that the effects of the variants on facial shape in EAS can be generalized to EUR. Based on this, we further identified 13 variants related to differences between facial shape in EUR and EAS populations. Evolutionary analyses suggest that the difference in nose shape between EUR and EAS populations is caused by a directional selection, due mainly to a local adaptation in Europeans. Our results illustrate the underlying genetic basis for facial differences across populations.

Epigenomic priming of immune genes implicates oligodendroglia in multiple sclerosis susceptibility

Multiple sclerosis (MS) is characterized by a targeted attack on oligodendroglia (OLG) and myelin by immune cells, which are thought to be the main drivers of MS susceptibility. We found that immune genes exhibit a primed chromatin state in single mouse and human OLG in a non-disease context, compatible with transitions to immune-competent states in MS. We identified BACH1 and STAT1 as transcription factors involved in immune gene regulation in oligodendrocyte precursor cells (OPCs). A subset of immune genes presents bivalency of H3K4me3/H3K27me3 in OPCs, with Polycomb inhibition leading to their increased activation upon interferon gamma (IFN-γ) treatment. Some MS susceptibility single-nucleotide polymorphisms (SNPs) overlap with these regulatory regions in mouse and human OLG. Treatment of mouse OPCs with IFN-γ leads to chromatin architecture remodeling at these loci and altered expression of interacting genes. Thus, the susceptibility for MS may involve OLG, which therefore constitutes novel targets for immunological-based therapies for MS.

CoRE-ATAC: A deep learning model for the functional classification of regulatory elements from single cell and bulk ATAC-seq data

Cis-Regulatory elements (cis-REs) include promoters, enhancers, and insulators that regulate gene expression programs via binding of transcription factors. ATAC-seq technology effectively identifies active cis-REs in a given cell type (including from single cells) by mapping accessible chromatin at base-pair resolution. However, these maps are not immediately useful for inferring specific functions of cis-REs. For this purpose, we developed a deep learning framework (CoRE-ATAC) with novel data encoders that integrate DNA sequence (reference or personal genotypes) with ATAC-seq cut sites and read pileups. CoRE-ATAC was trained on 4 cell types (n = 6 samples/replicates) and accurately predicted known cis-RE functions from 7 cell types (n = 40 samples) that were not used in model training (mean average precision = 0.80, mean F1 score = 0.70). CoRE-ATAC enhancer predictions from 19 human islet samples coincided with genetically modulated gain/loss of enhancer activity, which was confirmed by massively parallel reporter assays (MPRAs). Finally, CoRE-ATAC effectively inferred cis-RE function from aggregate single nucleus ATAC-seq (snATAC) data from human blood-derived immune cells that overlapped with known functional annotations in sorted immune cells, which established the efficacy of these models to study cis-RE functions of rare cells without the need for cell sorting. ATAC-seq maps from primary human cells reveal individual- and cell-specific variation in cis-RE activity. CoRE-ATAC increases the functional resolution of these maps, a critical step for studying regulatory disruptions behind diseases.

Perspective of the GEMSTONE Consortium on Current and Future Approaches to Functional Validation for Skeletal Genetic Disease Using Cellular, Molecular and Animal-Modeling Techniques

The availability of large human datasets for genome-wide association studies (GWAS) and the advancement of sequencing technologies have boosted the identification of genetic variants in complex and rare diseases in the skeletal field. Yet, interpreting results from human association studies remains a challenge. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary. Multiple unknowns exist for putative causal genes, including cellular localization of the molecular function. Intermediate traits ("endophenotypes"), e.g. molecular quantitative trait loci (molQTLs), are needed to identify mechanisms of underlying associations. Furthermore, index variants often reside in non-coding regions of the genome, therefore challenging for interpretation. Knowledge of non-coding variance (e.g. ncRNAs), repetitive sequences, and regulatory interactions between enhancers and their target genes is central for understanding causal genes in skeletal conditions. Animal models with deep skeletal phenotyping and cell culture models have already facilitated fine mapping of some association signals, elucidated gene mechanisms, and revealed disease-relevant biology. However, to accelerate research towards bridging the current gap between association and causality in skeletal diseases, alternative in vivo platforms need to be used and developed in parallel with the current -omics and traditional in vivo resources. Therefore, we argue that as a field we need to establish resource-sharing standards to collectively address complex research questions. These standards will promote data integration from various -omics technologies and functional dissection of human complex traits. In this mission statement, we review the current available resources and as a group propose a consensus to facilitate resource sharing using existing and future resources. Such coordination efforts will maximize the acquisition of knowledge from different approaches and thus reduce redundancy and duplication of resources. These measures will help to understand the pathogenesis of osteoporosis and other skeletal diseases towards defining new and more efficient therapeutic targets.

Genetic investigation of fibromuscular dysplasia identifies risk loci and shared genetics with common cardiovascular diseases

Fibromuscular dysplasia (FMD) is an arteriopathy associated with hypertension, stroke and myocardial infarction, affecting mostly women. We report results from the first genome-wide association meta-analysis of six studies including 1556 FMD cases and 7100 controls. We find an estimate of SNP-based heritability compatible with FMD having a polygenic basis, and report four robustly associated loci (PHACTR1, LRP1, ATP2B1, and LIMA1). Transcriptome-wide association analysis in arteries identifies one additional locus (SLC24A3). We characterize open chromatin in arterial primary cells and find that FMD associated variants are located in arterial-specific regulatory elements. Target genes are broadly involved in mechanisms related to actin cytoskeleton and intracellular calcium homeostasis, central to vascular contraction. We find significant genetic overlap between FMD and more common cardiovascular diseases and traits including blood pressure, migraine, intracranial aneurysm, and coronary artery disease.

TIGER: The gene expression regulatory variation landscape of human pancreatic islets

Genome-wide association studies (GWASs) identified hundreds of signals associated with type 2 diabetes (T2D). To gain insight into their underlying molecular mechanisms, we have created the translational human pancreatic islet genotype tissue-expression resource (TIGER), aggregating >500 human islet genomic datasets from five cohorts in the Horizon 2020 consortium T2DSystems. We impute genotypes using four reference panels and meta-analyze cohorts to improve the coverage of expression quantitative trait loci (eQTL) and develop a method to combine allele-specific expression across samples (cASE). We identify >1 million islet eQTLs, 53 of which colocalize with T2D signals. Among them, a low-frequency allele that reduces T2D risk by half increases CCND2 expression. We identify eight cASE colocalizations, among which we found a T2D-associated SLC30A8 variant. We make all data available through the TIGER portal (http://tiger.bsc.es), which represents a comprehensive human islet genomic data resource to elucidate how genetic variation affects islet function and translates into therapeutic insight and precision medicine for T2D.

Genetic variations of DNA bindings of FOXA1 and co-factors in breast cancer susceptibility

Identifying transcription factors (TFs) whose DNA bindings are altered by genetic variants that regulate susceptibility genes is imperative to understand transcriptional dysregulation in disease etiology. Here, we develop a statistical framework to analyze extensive ChIP-seq and GWAS data and identify 22 breast cancer risk-associated TFs. We find that, by analyzing genetic variations of TF-DNA bindings, the interaction of FOXA1 with co-factors such as ESR1 and E2F1, and the interaction of TFs with chromatin features (i.e., enhancers) play a key role in breast cancer susceptibility. Using genetic variants occupied by the 22 TFs, transcriptome-wide association analyses identify 52 previously unreported breast cancer susceptibility genes, including seven with evidence of essentiality from functional screens in breast relevant cell lines. We show that FOXA1 and co-factors form a core TF-transcriptional network regulating the susceptibility genes. Our findings provide additional insights into genetic variations of TF-DNA bindings (particularly for FOXA1) underlying breast cancer susceptibility.

The identification of transcription factors (TFs) whose binding sites are affected by risk genetic variants remains crucial. Here, the authors develop a statistical framework to analyse ChIP-seq and GWAS data, identify 22 breast cancer risk-associated TFs and a core TF-transcriptional network for FOXA1 and co-factors.

Mapping the genetic architecture of human traits to cell types in the kidney identifies mechanisms of disease and potential treatments

The functional interpretation of genome-wide association studies (GWAS) is challenging due to the cell-type-dependent influences of genetic variants. Here, we generated comprehensive maps of expression quantitative trait loci (eQTLs) for 659 microdissected human kidney samples and identified cell-type-eQTLs by mapping interactions between cell type abundances and genotypes. By partitioning heritability using stratified linkage disequilibrium score regression to integrate GWAS with single-cell RNA sequencing and single-nucleus assay for transposase-accessible chromatin with high-throughput sequencing data, we prioritized proximal tubules for kidney function and endothelial cells and distal tubule segments for blood pressure pathogenesis. Bayesian colocalization analysis nominated more than 200 genes for kidney function and hypertension. Our study clarifies the mechanism of commonly used antihypertensive and renal-protective drugs and identifies drug repurposing opportunities for kidney disease.

Disease-specific eQTL screening reveals an anti-fibrotic effect of AGXT2 in non-alcoholic fatty liver disease

BACKGROUND & AIMS

Non-alcoholic fatty liver disease (NAFLD) poses an increasing clinical burden. Genome-wide association studies have revealed a limited contribution of genomic variants to the disease, requiring alternative but robust approaches to identify disease-associated variants and genes. We carried out a disease-specific expression quantitative trait loci (eQTL) screen to identify novel genetic factors that specifically act on NAFLD progression on the basis of genotype.

METHODS

We recruited 125 Korean patients (83 with biopsy-proven NAFLD and 42 without NAFLD) and performed eQTL analyses using 21,272 transcripts and 3,234,941 genotyped and imputed single nucleotide polymorphisms. We then selected eQTLs that were detected only in the NAFLD group, but not in the control group (i.e., NAFLD-eQTLs). An additional cohort of 162 Korean individuals with NAFLD was used for replication. The function of the selected eQTL toward NAFLD development was validated using HepG2, primary hepatocytes and NAFLD mouse models.

RESULTS

The NAFLD-specific eQTL screening yielded 242 loci. Among them, AGXT2, encoding alanine-glyoxylate aminotransferase 2, displayed decreased expression in patients with NAFLD homozygous for the non-reference allele of rs2291702, compared to no-NAFLD individuals with the same genotype (p = 4.79 × 10^-6). This change was replicated in an additional 162 individuals, yielding a combined p value of 8.05 × 10^-8 from a total of 245 patients with NAFLD and 42 controls. Knockdown of AGXT2 induced palmitate-overloaded hepatocyte death by increasing endoplasmic reticulum stress, and exacerbated NAFLD diet-induced liver fibrosis in mice, while overexpression of AGXT2 attenuated liver fibrosis and steatosis.

CONCLUSIONS

We identified a new molecular role for AGXT2 in NAFLD. Our overall approach will serve as an efficient tool for uncovering novel genetic factors that contribute to liver steatosis and fibrosis in patients with NAFLD.

LAY SUMMARY

Elucidating causal genes for non-alcoholic fatty liver disease (NAFLD) has been challenging due to limited tissue availability and the polygenic nature of the disease. Using liver and blood samples from 125 Korean individuals (83 with NAFLD and 42 without NAFLD), we devised a new analytic method to identify causal genes. Among the candidates, we found that AGXT2-rs2291702 protects against liver fibrosis in a genotype-dependent manner with the potential for therapeutic interventions. Our approach enables the discovery of causal genes that act on the basis of genotype.

Motif-Raptor: a cell type-specific and transcription factor centric approach for post-GWAS prioritization of causal regulators

MOTIVATION

Genome-wide association studies (GWASs) have identified thousands of common trait-associated genetic variants but interpretation of their function remains challenging. These genetic variants can overlap the binding sites of transcription factors (TFs) and therefore could alter gene expression. However, we currently lack a systematic understanding on how this mechanism contributes to phenotype.

RESULTS

We present Motif-Raptor, a TF-centric computational tool that integrates sequence-based predictive models, chromatin accessibility, gene expression datasets and GWAS summary statistics to systematically investigate how TF function is affected by genetic variants. Given trait-associated non-coding variants, Motif-Raptor can recover relevant cell types and critical TFs to drive hypotheses regarding their mechanism of action. We tested Motif-Raptor on complex traits such as rheumatoid arthritis and red blood cell count and demonstrated its ability to prioritize relevant cell types, potential regulatory TFs and non-coding SNPs which have been previously characterized and validated.

AVAILABILITY AND IMPLEMENTATION

Motif-Raptor is freely available as a Python package at: https://github.com/pinellolab/MotifRaptor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

vSampler: fast and annotation-based matched variant sampling tool

SUMMARY

Sampling of control variants having matched properties with input variants is widely used in enrichment analysis of genome-wide association studies/quantitative trait loci and negative data construction for pathogenic/regulatory variant prediction methods. Spurious enrichment results because of confounding factors, such as minor allele frequency and linkage disequilibrium pattern, can be avoided by calibration of statistical significance based on matched controls. Here, we presented vSampler which can generate sets of randomly drawn variants with comprehensive choices of matching properties, such as tissue/cell type-specific epigenomic features. Importantly, the development of a novel data structure and sampling algorithms for vSampler makes it significantly fast than existing tools.

AVAILABILITY AND IMPLEMENTATION

vSampler web server and local program are available at http://mulinlab.org/vsampler.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dynamic CTCF binding directly mediates interactions among cis-regulatory elements essential for hematopoiesis

While constitutive CCCTC-binding factor (CTCF)-binding sites are needed to maintain relatively invariant chromatin structures, such as topologically associating domains, the precise roles of CTCF to control cell-type-specific transcriptional regulation remain poorly explored. We examined CTCF occupancy in different types of primary blood cells derived from the same donor to elucidate a new role for CTCF in gene regulation during blood cell development. We identified dynamic, cell-type-specific binding sites for CTCF that colocalize with lineage-specific transcription factors. These dynamic sites are enriched for single-nucleotide polymorphisms that are associated with blood cell traits in different linages, and they coincide with the key regulatory elements governing hematopoiesis. CRISPR-Cas9-based perturbation experiments demonstrated that these dynamic CTCF-binding sites play a critical role in red blood cell development. Furthermore, precise deletion of CTCF-binding motifs in dynamic sites abolished interactions of erythroid genes, such as RBM38, with their associated enhancers and led to abnormal erythropoiesis. These results suggest a novel, cell-type-specific function for CTCF in which it may serve to facilitate interaction of distal regulatory emblements with target promoters. Our study of the dynamic, cell-type-specific binding and function of CTCF provides new insights into transcriptional regulation during hematopoiesis.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.