Genetic and Epigenetic Fine Mapping of Complex Trait Associated Loci in the Human Liver

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

A landscape of gene expression regulation for synovium in arthritis

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes.

Proteomic insights into the associations between obesity, lifestyle factors, and coronary artery disease

BACKGROUND

We aimed to investigate the protein pathways linking obesity and lifestyle factors to coronary artery disease (CAD).

METHODS

Summary-level genome-wide association statistics of CAD were obtained from the CARDIoGRAMplusC4D consortium (60,801 cases and 123,504 controls) and the FinnGen study (R8, 39,036 cases and 303,463 controls). Proteome-wide Mendelian randomization (MR) analysis was conducted to identify CAD-associated blood proteins, supplemented by colocalization analysis to minimize potential bias caused by linkage disequilibrium. Two-sample MR analyses were performed to assess the associations of genetically predicted four obesity measures and 13 lifestyle factors with CAD risk and CAD-associated proteins' levels. A two-step network MR analysis was conducted to explore the mediating effects of proteins in the associations between these modifiable factors and CAD.

RESULTS

Genetically predicted levels of 41 circulating proteins were associated with CAD, and 17 of them were supported by medium to high colocalization evidence. PTK7 (protein tyrosine kinase-7), RGMB (repulsive guidance molecule BMP co-receptor B), TAGLN2 (transgelin-2), TIMP3 (tissue inhibitor of metalloproteinases 3), and VIM (vimentin) were identified as promising therapeutic targets. Several proteins were found to mediate the associations between some modifiable factors and CAD, with PCSK9, C1S, AGER (advanced glycosylation end product-specific receptor), and MST1 (mammalian Ste20-like kinase 1) exhibiting highest frequency among the mediating networks.

CONCLUSIONS

This study suggests pathways explaining the associations of obesity and lifestyle factors with CAD from alterations in blood protein levels. These insights may be used to prioritize therapeutic intervention for further study.

Web tools support predicting protein-nucleic acid complexes stability with affinity changes

Numerous biological processes, such as transcription, replication, and translation, rely on protein-nucleic acid interactions (PNIs). Demonstrating the binding stability of protein-nucleic acid complexes is vital to deciphering the code for PNIs. Numerous web-based tools have been developed to attach importance to protein-nucleic acid stability, facilitating the prediction of PNIs characteristics rapidly. However, the data and tools are dispersed and lack comprehensive integration to understand the stability of PNIs better. In this review, we first summarize existing databases for evaluating the stability of protein-nucleic acid binding. Then, we compare and evaluate the pros and cons of web tools for forecasting the interaction energies of protein-nucleic acid complexes. Finally, we discuss the application of combining models and capabilities of PNIs. We may hope these web-based tools will facilitate the discovery of recognition mechanisms for protein-nucleic acid binding stability. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Genome-wide association studies and fine-mapping identify genomic loci for n-3 and n-6 polyunsaturated fatty acids in Hispanic American and African American cohorts

Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acids (PUFAs) play critical roles in human health. Prior genome-wide association studies (GWAS) of n-3 and n-6 PUFAs in European Americans from the CHARGE Consortium have documented strong genetic signals in/near the FADS locus on chromosome 11. We performed a GWAS of four n-3 and four n-6 PUFAs in Hispanic American (n = 1454) and African American (n = 2278) participants from three CHARGE cohorts. Applying a genome-wide significance threshold of P < 5 × 10-8, we confirmed association of the FADS signal and found evidence of two additional signals (in DAGLA and BEST1) within 200 kb of the originally reported FADS signal. Outside of the FADS region, we identified novel signals for arachidonic acid (AA) in Hispanic Americans located in/near genes including TMX2, SLC29A2, ANKRD13D and POLD4, and spanning a > 9 Mb region on chromosome 11 (57.5 Mb ~ 67.1 Mb). Among these novel signals, we found associations unique to Hispanic Americans, including rs28364240, a POLD4 missense variant for AA that is common in CHARGE Hispanic Americans but absent in other race/ancestry groups. Our study sheds light on the genetics of PUFAs and the value of investigating complex trait genetics across diverse ancestry populations.

Characterizing the landscape of gene expression variance in humans

Gene expression variance has been linked to organismal function and fitness but remains a commonly neglected aspect of molecular research. As a result, we lack a comprehensive understanding of the patterns of transcriptional variance across genes, and how this variance is linked to context-specific gene regulation and gene function. Here, we use 57 large publicly available RNA-seq data sets to investigate the landscape of gene expression variance. These studies cover a wide range of tissues and allowed us to assess if there are consistently more or less variable genes across tissues and data sets and what mechanisms drive these patterns. We show that gene expression variance is broadly similar across tissues and studies, indicating that the pattern of transcriptional variance is consistent. We use this similarity to create both global and within-tissue rankings of variation, which we use to show that function, sequence variation, and gene regulatory signatures contribute to gene expression variance. Low-variance genes are associated with fundamental cell processes and have lower levels of genetic polymorphisms, have higher gene-gene connectivity, and tend to be associated with chromatin states associated with transcription. In contrast, high-variance genes are enriched for genes involved in immune response, environmentally responsive genes, immediate early genes, and are associated with higher levels of polymorphisms. These results show that the pattern of transcriptional variance is not noise. Instead, it is a consistent gene trait that seems to be functionally constrained in human populations. Furthermore, this commonly neglected aspect of molecular phenotypic variation harbors important information to understand complex traits and disease.

Genome-Wide Association Studies and fine-mapping of genomic loci for n-3 and n-6 Polyunsaturated Fatty Acids in Hispanic American and African American Cohorts

Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acids (PUFAs) play critical roles in human health. Prior genome-wide association studies (GWAS) of n-3 and n-6 PUFAs in European Americans from the CHARGE Consortium have documented strong genetic signals in/near the FADS locus on chromosome 11. We performed a GWAS of four n-3 and four n-6 PUFAs in Hispanic American (n = 1454) and African American (n = 2278) participants from three CHARGE cohorts. Applying a genome-wide significance threshold of P  < 5 x 10 - 8 , we confirmed association of the FADS signal and found evidence of two additional signals (in DAGLA and BEST1 ) within 200 kb of the originally reported FADS signal. Outside of the FADS region, we identified novel signals for arachidonic acid (AA) in Hispanic Americans located in/near genes including TMX2 , SLC29A2 , ANKRD13D and POLD4, and spanning a > 9 Mb region on chromosome 11 (57.5Mb ~ 67.1Mb). Among these novel signals, we found associations unique to Hispanic Americans, including rs28364240, a POLD4 missense variant for AA that is common in CHARGE Hispanic Americans but absent in other race/ancestry groups. Our study sheds light on the genetics of PUFAs and the value of investigating complex trait genetics across diverse ancestry populations.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Systematic analysis of the effects of genetic variants on chromatin accessibility to decipher functional variants in non-coding regions

Genome-wide association study (GWAS) has identified thousands of single nucleotide polymorphisms (SNPs) associated with complex diseases and traits. However, deciphering the functions of these SNPs still faces challenges. Recent studies have shown that SNPs could alter chromatin accessibility and result in differences in tumor susceptibility between individuals. Therefore, systematically analyzing the effects of SNPs on chromatin accessibility could help decipher the functions of SNPs, especially those in non-coding regions. Using data from The Cancer Genome Atlas (TCGA), chromatin accessibility quantitative trait locus (caQTL) analysis was conducted to estimate the associations between genetic variants and chromatin accessibility. We analyzed caQTLs in 23 human cancer types and identified 9,478 caQTLs in breast carcinoma (BRCA). In BRCA, these caQTLs tend to alter the binding affinity of transcription factors, and open chromatin regions regulated by these caQTLs are enriched in regulatory elements. By integrating with eQTL data, we identified 141 caQTLs showing a strong signal for colocalization with eQTLs. We also identified 173 caQTLs in genome-wide association studies (GWAS) loci and inferred several possible target genes of these caQTLs. By performing survival analysis, we found that ~10% caQTLs potentially influence the prognosis of patients. To facilitate access to relevant data, we developed a user-friendly data portal, BCaQTL (http://gong_lab.hzau.edu.cn/caqtl_database), for data searching and downloading. Our work may facilitate fine-map regulatory mechanisms underlying risk loci of cancer and discover the biomarkers or therapeutic targets for cancer prognosis. The BCaQTL database will be an important resource for genetic and epigenetic studies.

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.

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.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Cis-regulatory architecture of human ESC-derived hypothalamic neuron differentiation aids in variant-to-gene mapping of relevant complex traits

The hypothalamus regulates metabolic homeostasis by influencing behavior and endocrine systems. Given its role governing key traits, such as body weight and reproductive timing, understanding the genetic regulation of hypothalamic development and function could yield insights into disease pathogenesis. However, given its inaccessibility, studying human hypothalamic gene regulation has proven challenging. To address this gap, we generate a high-resolution chromatin architecture atlas of an established embryonic stem cell derived hypothalamic-like neuron model across three stages of in vitro differentiation. We profile accessible chromatin and identify physical contacts between gene promoters and putative cis-regulatory elements to characterize global regulatory landscape changes during hypothalamic differentiation. Next, we integrate these data with GWAS loci for various complex traits, identifying multiple candidate effector genes. Our results reveal common target genes for these traits, potentially affecting core developmental pathways. Our atlas will enable future efforts to determine hypothalamic mechanisms influencing disease susceptibility.

Identification of 22 susceptibility loci associated with testicular germ cell tumors

Testicular germ cell tumors (TGCT) are the most common tumor in young white men and have a high heritability. In this study, the international Testicular Cancer Consortium assemble 10,156 and 179,683 men with and without TGCT, respectively, for a genome-wide association study. This meta-analysis identifies 22 TGCT susceptibility loci, bringing the total to 78, which account for 44% of disease heritability. Men with a polygenic risk score (PRS) in the 95th percentile have a 6.8-fold increased risk of TGCT compared to men with median scores. Among men with independent TGCT risk factors such as cryptorchidism, the PRS may guide screening decisions with the goal of reducing treatment-related complications causing long-term morbidity in survivors. These findings emphasize the interconnected nature of two known pathways that promote TGCT susceptibility: male germ cell development within its somatic niche and regulation of chromosomal division and structure, and implicate an additional biological pathway, mRNA translation.

Genetic effects on liver chromatin accessibility identify disease regulatory variants

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

Transcriptional Silencers: Driving Gene Expression with the Brakes On

Silencers are regulatory DNA elements that reduce transcription from their target promoters; they are the repressive counterparts of enhancers. Although discovered decades ago, and despite evidence of their importance in development and disease, silencers have been much less studied than enhancers. Recently, however, a series of papers have reported systematic studies of silencers in various model systems. Silencers are often bifunctional regulatory elements that can also act as enhancers, depending on cellular context, and are enriched for expression quantitative trait loci (eQTLs) and disease-associated variants. There is not yet evidence of a 'silencer chromatin signature', in the distribution of histone modifications or associated proteins, that is common to all silencers; instead, silencers may fall into various subclasses, acting by distinct (and possibly overlapping) mechanisms.

Colocalization analysis of polycystic ovary syndrome to identify potential disease-mediating genes and proteins

Polycystic ovary syndrome (PCOS) is a common complex disease in women with a strong genetic component and downstream consequences for reproductive, metabolic and psychological health. There are currently 19 known PCOS risk loci, primarily identified in women of Han Chinese or European ancestry, and 14 of these risk loci were identified or replicated in a genome-wide association study of PCOS performed in up to 10,074 cases and 103,164 controls of European descent. However, for most of these loci the gene responsible for the association is unknown. We therefore use a Bayesian colocalization approach (Coloc) to highlight genes in PCOS-associated regions that may have a role in mediating the disease risk. We evaluated the posterior probabilities of evidence consistent with shared causal variants between 14 PCOS genetic risk loci and intermediate cellular phenotypes in one protein (N = 3301) and two expression quantitative trait locus datasets (N = 31,684 and N = 80–491). Through these analyses, we identified seven proteins or genes with evidence of a possibly shared causal variant for almost 30% of known PCOS signals, including follicle stimulating hormone and ERBB3, IKZF4, RPS26, SUOX, ZFP36L2, and C8orf49. Several of these potential effector proteins and genes have been implicated in the hypothalamic–pituitary–gonadal signalling pathway and provide an avenue for functional follow-up in order to demonstrate a causal role in PCOS pathophysiology.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

Exploiting the GTEx resources to decipher the mechanisms at GWAS loci

The resources generated by the GTEx consortium offer unprecedented opportunities to advance our understanding of the biology of human diseases. Here, we present an in-depth examination of the phenotypic consequences of transcriptome regulation and a blueprint for the functional interpretation of genome-wide association study-discovered loci. Across a broad set of complex traits and diseases, we demonstrate widespread dose-dependent effects of RNA expression and splicing. We develop a data-driven framework to benchmark methods that prioritize causal genes and find no single approach outperforms the combination of multiple approaches. Using colocalization and association approaches that take into account the observed allelic heterogeneity of gene expression, we propose potential target genes for 47% (2519 out of 5385) of the GWAS loci examined.

Biological constraints on GWAS SNPs at suggestive significance thresholds reveal additional BMI loci

To uncover novel significant association signals (p<5×10-8), genome-wide association studies (GWAS) requires increasingly larger sample sizes to overcome statistical correction for multiple testing. As an alternative, we aimed to identify associations among suggestive signals (5 × 10-8≤p<5×10-4) in increasingly powered GWAS efforts using chromatin accessibility and direct contact with gene promoters as biological constraints. We conducted retrospective analyses of three GIANT BMI GWAS efforts using ATAC-seq and promoter-focused Capture C data from human adipocytes and embryonic stem cell (ESC)-derived hypothalamic-like neurons. This approach, with its extremely low false-positive rate, identified 15 loci at p<5×10-5 in the 2010 GWAS, of which 13 achieved genome-wide significance by 2018, including at NAV1, MTIF3, and ADCY3. Eighty percent of constrained 2015 loci achieved genome-wide significance in 2018. We observed similar results in waist-to-hip ratio analyses. In conclusion, biological constraints on sub-significant GWAS signals can reveal potentially true-positive loci for further investigation in existing data sets without increasing sample size.

Epigenetics of Hepatic Insulin Resistance

Insulin resistance (IR) is largely recognized as a unifying feature that underlies metabolic dysfunction. Both lifestyle and genetic factors contribute to IR. Work from recent years has demonstrated that the epigenome may constitute an interface where different signals may converge to promote IR gene expression programs. Here, we review the current knowledge of the role of epigenetics in hepatic IR, focusing on the roles of DNA methylation and histone post-translational modifications. We discuss the broad epigenetic changes observed in the insulin resistant liver and its associated pathophysiological states and leverage on the wealth of 'omics' studies performed to discuss efforts in pinpointing specific loci that are disrupted by these changes. We envision that future studies, with increased genomic resolution and larger cohorts, will further the identification of biomarkers of early onset hepatic IR and assist the development of targeted interventions. Furthermore, there is growing evidence to suggest that persistent epigenetic marks may be acquired over prolonged exposure to disease or deleterious exposures, highlighting the need for preventative medicine and long-term lifestyle adjustments to avoid irreversible or long-term alterations in gene expression.

Novel approaches to liver disease diagnosis and modeling

Lack of a prompt and accurate diagnosis remains on top of the list of challenges faced by patients with rare liver diseases. Although rare liver diseases affect a significant percentage of the population as a group, when taken singularly they represent unique diseases and the approaches used for diagnosis of common liver diseases are insufficient. However, the development of new methods for the acquisition of molecular and clinical data (i.e., genomic, proteomics, metabolomics) and computational tools for their analysis and integration, together with advances in modeling diseases using stem cell-based technology [i.e., induced pluripotent stem cells (iPSCs) and tissue organoids] represent a promising and powerful tool to improve the clinical management of these patients. This is the goal of precision medicine, a novel approach of modern medicine that aims at delivering a specific treatment based on disease-specific biological insights and individual profile. This review will discuss the application and advances of these technologies and how they represent a new opportunity in hepatology.

From genome sequencing to the discovery of potential biomarkers in liver disease

Chronic liver disease progresses through several stages, fatty liver, steatohepatitis, cirrhosis, and eventually, it leads to hepatocellular carcinoma (HCC) over a long period of time. Since a large proportion of patients with HCC are accompanied by cirrhosis, it is considered to be an important factor in the diagnosis of liver cancer. This is because cirrhosis leads to an irreversible harmful effect, but the early stages of chronic liver disease could be reversed to a healthy state. Therefore, the discovery of biomarkers that could identify the early stages of chronic liver disease is important to prevent serious liver damage. Biomarker discovery at liver cancer and cirrhosis has enhanced the development of sequencing technology. Next generation sequencing (NGS) is one of the representative technical innovations in the biological field in the recent decades and it is the most important thing to design for research on what type of sequencing methods are suitable and how to handle the analysis steps for data integration. In this review, we comprehensively summarized NGS techniques for identifying genome, transcriptome, DNA methylome and 3D/4D chromatin structure, and introduced framework of processing data set and integrating multi-omics data for uncovering biomarkers.

Systematically Sifting Big Data to Identify Novel Causal Genes for Human Traits

Widespread technological advances have propelled human genetics into a "big data" era, in which genome-wide data from extremely large cohorts can be integrated with other "-omics" datasets from humans and model systems. Li et al. (2020) demonstrate the power of applying multiple computational analyses to publicly available data to prioritize the study of genes with novel trait associations.

Identification of PIEZO1 polymorphisms for human bone mineral density

Bone mineral density (BMD) is a key indicator for diagnosis and treatment for osteoporosis; the reduction of BMD could increase the risk of osteoporotic fracture. It was very recently found that Piezo1 mediated mechanically evoked responses in bone and further participated in bone formation in mice. Here, we performed cross phenotype meta-analysis for human BMD at lumbar spine (LS), femoral neck (FN), distal radius/forearm (FA) and heel and screened out 14 top SNPs for PIEZO1, these SNPs were overlapped with putative enhancers, DNase-I hypersensitive sites and active promoter flanking regions. We found that the signal of the best SNP rs62048221 was mainly from heel ultrasound estimated BMD (-0.02 SD per T allele, P = 8.50E-09), where calcaneus supported most of the mechanical force of body when standing, walking and doing physical exercises. Each copy of the effect allele T of SNP rs62048221 was associated with a decrease of 0.0035 g/cm² BMD (P = 4.6E-27, SE = 0.0003) in UK Biobank data within 477,760 samples. SNP rs62048221 was located at the enhancer region (HEDD enhancer ID 2331049) of gene PIEZO1, site-directed ChIP assays in human mesenchymal stem cells (hMSCs) showed significant enrichment of H3K4me1 and H3K27ac in this region, luciferase assays showed that rs62048221 could significantly affect the activity of the enhancer where it resides. Our results first suggested that SNP rs62048221 might mediate the PIEZO1 expression level via modulating the activity of cis-regulatory elements and then further affect the BMD.

Liver gene regulatory networks: Contributing factors to nonalcoholic fatty liver disease

Metabolic diseases such as nonalcoholic fatty liver disease (NAFLD) result from complex interactions between intrinsic and extrinsic factors, including genetics and exposure to obesogenic environments. These risk factors converge in aberrant gene expression patterns in the liver, which are underlined by altered cis-regulatory networks. In homeostasis and in disease states, liver cis-regulatory networks are established by coordinated action of liver-enriched transcription factors (TFs), which define enhancer landscapes, activating broad gene programs with spatiotemporal resolution. Recent advances in DNA sequencing have dramatically expanded our ability to map active transcripts, enhancers and TF cistromes, and to define the 3D chromatin topology that contains these elements. Deployment of these technologies has allowed investigation of the molecular processes that regulate liver development and metabolic homeostasis. Moreover, genomic studies of NAFLD patients and NAFLD models have demonstrated that the liver undergoes pervasive regulatory rewiring in NAFLD, which is reflected by aberrant gene expression profiles. We have therefore achieved an unprecedented level of detail in the understanding of liver cis-regulatory networks, particularly in physiological conditions. Future studies should aim to map active regulatory elements with added levels of resolution, addressing how the chromatin landscapes of different cell lineages contribute to and are altered in NAFLD and NAFLD-associated metabolic states. Such efforts would provide additional clues into the molecular factors that trigger this disease. This article is categorized under: Biological Mechanisms > Metabolism Biological Mechanisms > Regulatory Biology Laboratory Methods and Technologies > Genetic/Genomic Methods.