A mega-analysis of expression quantitative trait loci (eQTL) provides insight into the regulatory architecture of gene expression variation in liver

Levels of the HtrA1 Protein in Serum and Vitreous Humor Are Independent of Genetic Risk for Age-Related Macular Degeneration at the 10q26 Locus

Purpose

The purpose of this study was to determine if levels of the HtrA1 protein in serum or vitreous humor are influenced by genetic risk for age-related macular degeneration (AMD) at the 10q26 locus, age, sex, AMD status, and/or AMD disease severity, and, therefore, to determine the contribution of systemic and ocular HtrA1 to the AMD disease process.

Methods

A custom-made sandwich ELISA assay (SCTM ELISA) for detection of the HtrA1 protein was designed and compared with three commercial assays (R&D Systems, MyBiosource 1 and MyBiosource 2) using 65 serum samples. Concentrations of HtrA1 were thereafter determined in serum and vitreous samples collected from 248 individuals and 145 human donor eyes, respectively.

Results

The SCTM ELISA demonstrated high specificity, good recovery, and parallelism within its linear detection range and performed comparably to the R&D Systems assay. In contrast, we were unable to demonstrate the specificity of the two assays from MyBioSource using either recombinant or native HtrA1. Analyses of concentrations obtained using the validated SCTM assay revealed that genetic risk at the 10q26 locus, age, sex, or AMD status are not significantly associated with altered levels of the HtrA1 protein in serum or in vitreous humor (P > 0.05).

Conclusions

HtrA1 levels in serum and vitreous do not reflect the risk for AMD associated with the 10q26 locus or disease status. Localized alteration in HTRA1 expression in the retinal pigment epithelium, rather than systemic changes in HtrA1, is the most likely driver of elevated risk for developing AMD among individuals with risk variants at the 10q26 locus.

Transcriptome-wide association study of the plasma proteome reveals cis and trans regulatory mechanisms underlying complex traits

Regulation of transcription and translation are mechanisms through which genetic variants affect complex traits. Expression quantitative trait locus (eQTL) studies have been more successful at identifying cis-eQTL (within 1 Mb of the transcription start site) than trans-eQTL. Here, we tested the cis component of gene expression for association with observed plasma protein levels to identify cis- and trans-acting genes that regulate protein levels. We used transcriptome prediction models from 49 Genotype-Tissue Expression (GTEx) Project tissues to predict the cis component of gene expression and tested the predicted expression of every gene in every tissue for association with the observed abundance of 3,622 plasma proteins measured in 3,301 individuals from the INTERVAL study. We tested significant results for replication in 971 individuals from the Trans-omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA). We found 1,168 and 1,210 cis- and trans-acting associations that replicated in TOPMed (FDR < 0.05) with a median expected true positive rate (π1) across tissues of 0.806 and 0.390, respectively. The target proteins of trans-acting genes were enriched for transcription factor binding sites and autoimmune diseases in the GWAS catalog. Furthermore, we found a higher correlation between predicted expression and protein levels of the same underlying gene (R = 0.17) than observed expression (R = 0.10, p = 7.50 × 10-11). This indicates the cis-acting genetically regulated (heritable) component of gene expression is more consistent across tissues than total observed expression (genetics + environment) and is useful in uncovering the function of SNPs associated with complex traits.

Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single hematopoietic stem cell lineage. Using whole-genome sequencing of 45,510 Icelandic and 130,709 UK Biobank participants combined with a mutational barcode method, we identified 16,306 people with CH. Prevalence approaches 50% in elderly participants. Smoking demonstrates a dosage-dependent impact on risk of CH. CH associates with several smoking-related diseases. Contrary to published claims, we find no evidence that CH is associated with cardiovascular disease. We provide evidence that CH is driven by genes that are commonly mutated in myeloid neoplasia and implicate several new driver genes. The presence and nature of a driver mutation alters the risk profile for hematological disorders. Nevertheless, most CH cases have no known driver mutations. A CH genome-wide association study identified 25 loci, including 19 not implicated previously in CH. Splicing, protein and expression quantitative trait loci were identified for CD164 and TCL1A.

Insights into the liver-eyes connections, from epidemiological, mechanical studies to clinical translation

Maintenance of internal homeostasis is a sophisticated process, during which almost all organs get involved. Liver plays a central role in metabolism and involves in endocrine, immunity, detoxification and storage, and therefore it communicates with distant organs through such mechanisms to regulate pathophysiological processes. Dysfunctional liver is often accompanied by pathological phenotypes of distant organs, including the eyes. Many reviews have focused on crosstalk between the liver and gut, the liver and brain, the liver and heart, the liver and kidney, but with no attention paid to the liver and eyes. In this review, we summarized intimate connections between the liver and the eyes from three aspects. Epidemiologically, we suggest liver-related, potential, protective and risk factors for typical eye disease as well as eye indicators connected with liver status. For molecular mechanism aspect, we elaborate their inter-organ crosstalk from metabolism (glucose, lipid, proteins, vitamin, and mineral), detoxification (ammonia and bilirubin), and immunity (complement and inflammation regulation) aspect. In clinical application part, we emphasize the latest advances in utilizing the liver-eye axis in disease diagnosis and therapy, involving artificial intelligence-deep learning-based novel diagnostic tools for detecting liver disease and adeno-associated viral vector-based gene therapy method for curing blinding eye disease. We aim to focus on and provide novel insights into liver and eyes communications and help resolve existed clinically significant issues.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing

BACKGROUND

A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation.

RESULTS

In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10^-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10^-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10^-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework.

CONCLUSIONS

Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

Genetic variability shapes the alternative pathway complement activity and predisposition to complement-related diseases

The implementation of next-generation sequencing technologies has provided a sharp picture of the genetic variability in the components and regulators of the alternative pathway (AP) of the complement system and has revealed the association of many AP variants with different rare and common diseases. An important finding that has emerged from these analyses is that each of these complement-related diseases associate with genetic variants altering specific aspects of the activation and regulation of the AP. These genotype-phenotype correlations have provided valuable insights into their pathogenic mechanisms with important diagnostic and therapeutic implications. While genetic variants in coding regions and structural variants are reasonably well characterized and occasionally have been instrumental to uncover unknown features of the complement proteins, data about complement expressed quantitative trait loci are still very limited. A crucial task for future studies will be to identify these quantitative variations and to determine their impact in the overall activity of the AP. This is fundamental as it is now clear that the consequences of genetic variants in the AP are additive and that susceptibility or resistance to disease is the result of specific combinations of genetic variants in different complement components and regulators ("complotypes").

Human variation in population-wide gene expression data predicts gene perturbation phenotype

Population-scale datasets of healthy individuals capture genetic and environmental factors influencing gene expression. The expression variance of a gene of interest (GOI) can be exploited to set up a quasi loss- or gain-of-function "in population" experiment. We describe here an approach, huva (human variation), taking advantage of population-scale multi-layered data to infer gene function and relationships between phenotypes and expression. Within a reference dataset, huva derives two experimental groups with LOW or HIGH expression of the GOI, enabling the subsequent comparison of their transcriptional profile and functional parameters. We demonstrate that this approach robustly identifies the phenotypic relevance of a GOI allowing the stratification of genes according to biological functions, and we generalize this concept to almost 16,000 genes in the human transcriptome. Additionally, we describe how huva predicts monocytes to be the major cell type in the pathophysiology of STAT1 mutations, evidence validated in a clinical cohort.

The Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel

Background

Unraveling the relationship between genetic variation and phenotypic traits remains a fundamental challenge in biology. Mapping variants underlying complex traits while controlling for confounding environmental factors is often problematic. To address this, we establish a vertebrate genetic resource specifically to allow for robust genotype-to-phenotype investigations. The teleost medaka (Oryzias latipes) is an established genetic model system with a long history of genetic research and a high tolerance to inbreeding from the wild.

Results

Here we present the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel: the first near-isogenic panel of 80 inbred lines in a vertebrate model derived from a wild founder population. Inbred lines provide fixed genomes that are a prerequisite for the replication of studies, studies which vary both the genetics and environment in a controlled manner, and functional testing. The MIKK panel will therefore enable phenotype-to-genotype association studies of complex genetic traits while allowing for careful control of interacting factors, with numerous applications in genetic research, human health, drug development, and fundamental biology.

Conclusions

Here we present a detailed characterization of the genetic variation across the MIKK panel, which provides a rich and unique genetic resource to the community by enabling large-scale experiments for mapping complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02623-z.

Interpreting Functional Impact of Genetic Variations by Network QTL for Genotype–Phenotype Association Study

An enormous challenge in the post-genome era is to annotate and resolve the consequences of genetic variation on diverse phenotypes. The genome-wide association study (GWAS) is a well-known method to identify potential genetic loci for complex traits from huge genetic variations, following which it is crucial to identify expression quantitative trait loci (eQTL). However, the conventional eQTL methods usually disregard the systematical role of single-nucleotide polymorphisms (SNPs) or genes, thereby overlooking many network-associated phenotypic determinates. Such a problem motivates us to recognize the network-based quantitative trait loci (QTL), i.e., network QTL (nQTL), which is to detect the cascade association as genotype → network → phenotype rather than conventional genotype → expression → phenotype in eQTL. Specifically, we develop the nQTL framework on the theory and approach of single-sample networks, which can identify not only network traits (e.g., the gene subnetwork associated with genotype) for analyzing complex biological processes but also network signatures (e.g., the interactive gene biomarker candidates screened from network traits) for characterizing targeted phenotype and corresponding subtypes. Our results show that the nQTL framework can efficiently capture associations between SNPs and network traits (i.e., edge traits) in various simulated data scenarios, compared with traditional eQTL methods. Furthermore, we have carried out nQTL analysis on diverse biological and biomedical datasets. Our analysis is effective in detecting network traits for various biological problems and can discover many network signatures for discriminating phenotypes, which can help interpret the influence of nQTL on disease subtyping, disease prognosis, drug response, and pathogen factor association. Particularly, in contrast to the conventional approaches, the nQTL framework could also identify many network traits from human bulk expression data, validated by matched single-cell RNA-seq data in an independent or unsupervised manner. All these results strongly support that nQTL and its detection framework can simultaneously explore the global genotype–network–phenotype associations and the underlying network traits or network signatures with functional impact and importance.

Chromatin accessibility and gene expression during adipocyte differentiation identify context-dependent effects at cardiometabolic GWAS loci

Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and mechanisms at genome-wide association study (GWAS) loci. To identify regulatory elements that display differential activity across adipocyte differentiation, we performed ATAC-seq and RNA-seq in a human cell model of preadipocytes and adipocytes at days 4 and 14 of differentiation. For comparison, we created a consensus map of ATAC-seq peaks in 11 human subcutaneous adipose tissue samples. We identified 58,387 context-dependent chromatin accessibility peaks and 3,090 context-dependent genes between all timepoint comparisons (log2 fold change>1, FDR<5%) with 15,919 adipocyte- and 18,244 preadipocyte-dependent peaks. Adipocyte-dependent peaks showed increased overlap (60.1%) with Roadmap Epigenomics adipocyte nuclei enhancers compared to preadipocyte-dependent peaks (11.5%). We linked context-dependent peaks to genes based on adipocyte promoter capture Hi-C data, overlap with adipose eQTL variants, and context-dependent gene expression. Of 16,167 context-dependent peaks linked to a gene, 5,145 were linked by two or more strategies to 1,670 genes. Among GWAS loci for cardiometabolic traits, adipocyte-dependent peaks, but not preadipocyte-dependent peaks, showed significant enrichment (LD score regression P<0.005) for waist-to-hip ratio and modest enrichment (P < 0.05) for HDL-cholesterol. We identified 659 peaks linked to 503 genes by two or more approaches and overlapping a GWAS signal, suggesting a regulatory mechanism at these loci. To identify variants that may alter chromatin accessibility between timepoints, we identified 582 variants in 454 context-dependent peaks that demonstrated allelic imbalance in accessibility (FDR<5%), of which 55 peaks also overlapped GWAS variants. At one GWAS locus for palmitoleic acid, rs603424 was located in an adipocyte-dependent peak linked to SCD and exhibited allelic differences in transcriptional activity in adipocytes (P = 0.003) but not preadipocytes (P = 0.09). These results demonstrate that context-dependent peaks and genes can guide discovery of regulatory variants at GWAS loci and aid identification of regulatory mechanisms.

Cardiovascular and metabolic diseases are widespread, and an increased understanding of genetic mechanisms behind these diseases could improve treatment. Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and genetic mechanisms for disease traits. A relevant context for cardiovascular and metabolic disease traits is adipocyte differentiation. To identify regulatory elements and genes that display differences in activity during adipocyte differentiation, we profiled chromatin accessibility and gene expression in a human cell model of preadipocytes and adipocytes. We identified chromatin regions that change accessibility during differentiation and predicted genes they may affect. We also linked these chromatin regions to genetic variants associated with risk of disease. At one genomic region linked to fatty acids, a chromatin region more accessible in adipocytes linked to a fatty acid synthesis gene and exhibited allelic differences in transcriptional activity in adipocytes but not preadipocytes. These results demonstrate that chromatin regions and genes that change during cell context can guide discovery of regulatory variants and aid identification of disease mechanisms.

Assigning Co-Regulated Human Genes and Regulatory Gene Clusters

Elucidating the role of genetic variation in the regulation of gene expression is key to understanding the pathobiology of complex diseases which, in consequence, is crucial in devising targeted treatment options. Expression quantitative trait locus (eQTL) analysis correlates a genetic variant with the strength of gene expression, thus defining thousands of regulated genes in a multitude of human cell types and tissues. Some eQTL may not act independently of each other but instead may be regulated in a coordinated fashion by seemingly independent genetic variants. To address this issue, we combined the approaches of eQTL analysis and colocalization studies. Gene expression was determined in datasets comprising 49 tissues from the Genotype-Tissue Expression (GTEx) project. From about 33,000 regulated genes, over 14,000 were found to be co-regulated in pairs and were assembled across all tissues to almost 15,000 unique clusters containing up to nine regulated genes affected by the same eQTL signal. The distance of co-regulated eGenes was, on average, 112 kilobase pairs. Of 713 genes known to express clinical symptoms upon haploinsufficiency, 231 (32.4%) are part of at least one of the identified clusters. This calls for caution should treatment approaches aim at an upregulation of a haploinsufficient gene. In conclusion, we present an unbiased approach to identifying co-regulated genes in and across multiple tissues. Knowledge of such common effects is crucial to appreciate implications on biological pathways involved, specifically when a treatment option targets a co-regulated disease gene.

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.

Associations between pancreatic expression quantitative traits and risk of pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is among the most lethal cancers. Its poor prognosis is predominantly due to the fact that most patients remain asymptomatic until the disease reaches an advanced stage, alongside the lack of early markers and screening strategies. A better understanding of PDAC risk factors is essential for the identification of groups at high risk in the population. Genome-wide association studies (GWAS) have been a powerful tool for detecting genetic variants associated with complex traits, including pancreatic cancer. By exploiting functional and GWAS data, we investigated the associations between polymorphisms affecting gene function in the pancreas (expression quantitative trait loci, eQTLs) and PDAC risk. In a two-phase approach, we analysed 13 713 PDAC cases and 43 784 controls and identified a genome-wide significant association between the A allele of the rs2035875 polymorphism and increased PDAC risk (P = 7.14 × 10-10). This allele is known to be associated with increased expression in the pancreas of the keratin genes KRT8 and KRT18, whose increased levels have been reported to correlate with various tumour cell characteristics. Additionally, the A allele of the rs789744 variant was associated with decreased risk of developing PDAC (P = 3.56 × 10-6). This single nucleotide polymorphism is situated in the SRGAP1 gene and the A allele is associated with higher expression of the gene, which in turn inactivates the cyclin-dependent protein 42 (CDC42) gene expression, thus decreasing the risk of PDAC. In conclusion, we present here a functional-based novel PDAC risk locus and an additional strong candidate supported by significant associations and plausible biological mechanisms.

Common haplotypes at the CFH locus and low-frequency variants in CFHR2 and CFHR5 associate with systemic FHR concentrations and age-related macular degeneration

Age-related macular degeneration (AMD) is the principal cause of blindness in the elderly population. A strong effect on AMD risk has been reported for genetic variants at the CFH locus, encompassing complement factor H (CFH) and the complement-factor-H-related (CFHR) genes, but the underlying mechanisms are not fully understood. We aimed to dissect the role of factor H (FH) and FH-related (FHR) proteins in AMD in a cohort of 202 controls and 216 individuals with AMD. We detected elevated systemic levels of FHR-1 (p = 1.84 × 10-6), FHR-2 (p = 1.47 × 10-4), FHR-3 (p = 1.05 × 10-5) and FHR-4A (p = 1.22 × 10-2) in AMD, whereas FH concentrations remained unchanged. Common AMD genetic variants and haplotypes at the CFH locus strongly associated with FHR protein concentrations (e.g., FH p.Tyr402His and FHR-2 concentrations, p = 3.68 × 10-17), whereas the association with FH concentrations was limited. Furthermore, in an International AMD Genomics Consortium cohort of 17,596 controls and 15,894 individuals with AMD, we found that low-frequency and rare protein-altering CFHR2 and CFHR5 variants associated with AMD independently of all previously reported genome-wide association study (GWAS) signals (p = 5.03 × 10-3 and p = 2.81 × 10-6, respectively). Low-frequency variants in CFHR2 and CFHR5 led to reduced or absent FHR-2 and FHR-5 concentrations (e.g., p.Cys72Tyr in CFHR2 and FHR-2, p = 2.46 × 10-16). Finally, we showed localization of FHR-2 and FHR-5 in the choriocapillaris and in drusen. Our study identifies FHR proteins as key proteins in the AMD disease mechanism. Consequently, therapies that modulate FHR proteins might be effective for treating or preventing progression of AMD. Such therapies could target specific individuals with AMD on the basis of their genotypes at the CFH locus.

GWAS of serum ALT and AST reveals an association of SLC30A10 Thr95Ile with hypermanganesemia symptoms

Understanding mechanisms of hepatocellular damage may lead to new treatments for liver disease, and genome-wide association studies (GWAS) of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) serum activities have proven useful for investigating liver biology. Here we report 100 loci associating with both enzymes, using GWAS across 411,048 subjects in the UK Biobank. The rare missense variant SLC30A10 Thr95Ile (rs188273166) associates with the largest elevation of both enzymes, and this association replicates in the DiscovEHR study. SLC30A10 excretes manganese from the liver to the bile duct, and rare homozygous loss of function causes the syndrome hypermanganesemia with dystonia-1 (HMNDYT1) which involves cirrhosis. Consistent with hematological symptoms of hypermanganesemia, SLC30A10 Thr95Ile carriers have increased hematocrit and risk of iron deficiency anemia. Carriers also have increased risk of extrahepatic bile duct cancer. These results suggest that genetic variation in SLC30A10 adversely affects more individuals than patients with diagnosed HMNDYT1.

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.

Epistatic interactions of genetic loci associated with age-related macular degeneration

The currently largest genome-wide association study (GWAS) for age-related macular degeneration (AMD) defines disease association with genome-wide significance for 52 independent common and rare genetic variants across 34 chromosomal loci. Overall, these loci contain over 7200 variants and are enriched for genes with functions indicating several shared cellular processes. Still, the precise mechanisms leading to AMD pathology are largely unknown. Here, we exploit the phenomenon of epistatic interaction to identify seemingly independent AMD-associated variants that reveal joint effects on gene expression. We focus on genetic variants associated with lipid metabolism, organization of extracellular structures, and innate immunity, specifically the complement cascade. Multiple combinations of independent variants were used to generate genetic risk scores allowing gene expression in liver to be compared between low and high-risk AMD. We identified genetic variant combinations correlating significantly with expression of 26 genes, of which 19 have not been associated with AMD before. This study defines novel targets and allows prioritizing further functional work into AMD pathobiology.

TreeMap: a structured approach to fine mapping of eQTL variants

Motivation

Expression quantitative trait loci (eQTL) harbor genetic variants modulating gene transcription. Fine mapping of regulatory variants at these loci is a daunting task due to the juxtaposition of causal and linked variants at a locus as well as the likelihood of interactions among multiple variants. This problem is exacerbated in genes with multiple cis-acting eQTL, where superimposed effects of adjacent loci further distort the association signals.

Results

We developed a novel algorithm, TreeMap, that identifies putative causal variants in cis-eQTL accounting for multisite effects and genetic linkage at a locus. Guided by the hierarchical structure of linkage disequilibrium, TreeMap performs an organized search for individual and multiple causal variants. Via extensive simulations, we show that TreeMap detects co-regulating variants more accurately than current methods. Furthermore, its high computational efficiency enables genome-wide analysis of long-range eQTL. We applied TreeMap to GTEx data of brain hippocampus samples and transverse colon samples to search for eQTL in gene bodies and in 4 Mbps gene-flanking regions, discovering numerous distal eQTL. Furthermore, we found concordant distal eQTL that were present in both brain and colon samples, implying long-range regulation of gene expression.

Availability and implementation

TreeMap is available as an R package enabled for parallel processing at https://github.com/liliulab/treemap.

Supplementary information

Supplementary data are available at Bioinformatics online.

A systematic comparison of normalization methods for eQTL analysis

Expression quantitative trait loci (eQTL) analysis has been widely used in interpreting disease-associated loci through correlating genetic variant loci with the expression of specific genes. RNA-sequencing (RNA-Seq), which can quantify gene expression at the genome-wide level, is often used in eQTL identification. Since different normalization methods of gene expression have substantial impacts on RNA-seq downstream analysis, it is of great necessity to systematically compare the effects of these methods on eQTL identification. Here, by using RNA-seq and genotype data of four different cancers in The Cancer Genome Atlas (TCGA) database, we comprehensively evaluated the effect of eight commonly used normalization methods on eQTL identification. Our results showed that the application of different methods could cause 20-30% differences in the final results of eQTL identification. Among these methods, COUNT, Median of Ratio (MED) and Trimmed Mean of M-values (TMM) generated similar results for identifying eQTLs, while Fragments Per Kilobase Million (FPKM) or RANK produced more differential results compared with other methods. Based on the accuracy and receiver operating characteristic (ROC) curve, the TMM method was found to be the optimal method for normalizing gene expression data in eQTLs analysis. In addition, we also evaluated the performance of different pairwise combinations of these methods. As a result, compared with single normalization methods, the combination of methods can not only identify more cis-eQTLs, but also improve the performance of the ROC curve. Overall, this study provides a comprehensive comparison of normalization methods for identifying eQTLs from RNA-seq data, and proposes some practical recommendations for diverse scenarios.

MRLocus: Identifying causal genes mediating a trait through Bayesian estimation of allelic heterogeneity

Expression quantitative trait loci (eQTL) studies are used to understand the regulatory function of non-coding genome-wide association study (GWAS) risk loci, but colocalization alone does not demonstrate a causal relationship of gene expression affecting a trait. Evidence for mediation, that perturbation of gene expression in a given tissue or developmental context will induce a change in the downstream GWAS trait, can be provided by two-sample Mendelian Randomization (MR). Here, we introduce a new statistical method, MRLocus, for Bayesian estimation of the gene-to-trait effect from eQTL and GWAS summary data for loci with evidence of allelic heterogeneity, that is, containing multiple causal variants. MRLocus makes use of a colocalization step applied to each nearly-LD-independent eQTL, followed by an MR analysis step across eQTLs. Additionally, our method involves estimation of the extent of allelic heterogeneity through a dispersion parameter, indicating variable mediation effects from each individual eQTL on the downstream trait. Our method is evaluated against other state-of-the-art methods for estimation of the gene-to-trait mediation effect, using an existing simulation framework. In simulation, MRLocus often has the highest accuracy among competing methods, and in each case provides more accurate estimation of uncertainty as assessed through interval coverage. MRLocus is then applied to five candidate causal genes for mediation of particular GWAS traits, where gene-to-trait effects are concordant with those previously reported. We find that MRLocus's estimation of the causal effect across eQTLs within a locus provides useful information for determining how perturbation of gene expression or individual regulatory elements will affect downstream traits. The MRLocus method is implemented as an R package available at https://mikelove.github.io/mrlocus.

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.

Implications of genetic variation in the complement system in age-related macular degeneration

Age-related macular degeneration (AMD) is the main cause of vision loss among the elderly in the Western world. While AMD is a multifactorial disease, the complement system was identified as one of the main pathways contributing to disease risk. The strong link between the complement system and AMD was demonstrated by genetic associations, and by elevated complement activation in local eye tissue and in the systemic circulation of AMD patients. Several complement inhibitors have been and are being explored in clinical trials, but thus far with limited success, leaving the majority of AMD patients without treatment options to date. This indicates that there is still a gap of knowledge regarding the functional implications of the complement system in AMD pathogenesis and how to bring these towards clinical translation. Many different experimental set-ups and disease models have been used to study complement activation in vivo and in vitro, and recently emerging patient-derived induced pluripotent stem cells and genome-editing techniques open new opportunities to study AMD disease mechanisms and test new therapeutic strategies in the future. In this review we provide an extensive overview of methods employed to understand the molecular processes of complement activation in AMD pathogenesis. We discuss the findings, advantages and challenges of each approach and conclude with an outlook on how recent, exciting developments can fill in current knowledge gaps and can aid in the development of effective complement-targeting therapeutic strategies in AMD.

geno5mC: A Database to Explore the Association between Genetic Variation (SNPs) and CpG Methylation in the Human Genome

Genetic variation, gene expression and DNA methylation influence each other in a complex way. To study the impact of sequence variation and DNA methylation on gene expression, we generated geno⁵mC, a database that contains statistically significant SNP-CpG associations that are biologically classified either through co-localization with known regulatory regions (promoters and enhancers), or through known correlations with the expression levels of nearby genes. The SNP rs727563 can be used to illustrate the usefulness of this approach. This SNP has been associated with inflammatory bowel disease through GWAS, but it is not located near any gene related to this phenotype. However, geno⁵mC reveals that rs727563 is associated with the methylation state of several CpGs located in promoter regions of genes reported to be involved in inflammatory processes. This case exemplifies how geno⁵mC can be used to infer relevant and previously unknown interactions between described disease-associated SNPs and their functional targets.

Learning from Fifteen Years of Genome-Wide Association Studies in Age-Related Macular Degeneration

Over the last 15 years, genome-wide association studies (GWAS) have greatly advanced our understanding of the genetic landscape of complex phenotypes. Nevertheless, causal interpretations of GWAS data are challenging but crucial to understand underlying mechanisms and pathologies. In this review, we explore to what extend the research community follows up on GWAS data. We have traced the scientific activities responding to the two largest GWAS conducted on age-related macular degeneration (AMD) so far. Altogether 703 articles were manually categorized according to their study type. This demonstrates that follow-up studies mainly involve "Review articles" (33%) or "Genetic association studies" (33%), while 19% of publications report on findings from experimental work. It is striking to note that only three of 16 AMD-associated loci described de novo in 2016 were examined in the four-year follow-up period after publication. A comparative analysis of five studies on gene expression regulation in AMD-associated loci revealed consistent gene candidates for 15 of these loci. Our random survey highlights the fact that functional follow-up studies on GWAS results are still in its early stages hampering a significant refinement of the vast association data and thus a more accurate insight into mechanisms and pathways.

Multi-layered epigenetic regulation of IRS2 expression in the liver of obese individuals with type 2 diabetes

AIMS/HYPOTHESIS

IRS2 is an important molecular switch that mediates insulin signalling in the liver. IRS2 dysregulation is responsible for the phenomenon of selective insulin resistance that is observed in type 2 diabetes. We hypothesise that epigenetic mechanisms are involved in the regulation of IRS2 in the liver of obese and type 2 diabetic individuals.

METHODS

DNA methylation of seven CpG sites was studied by bisulphite pyrosequencing and mRNA and microRNA (miRNA) expression was assessed by quantitative real-time PCR in liver biopsies of 50 obese non-diabetic and 31 obese type 2 diabetic participants, in a cross-sectional setting. Methylation-sensitive luciferase assays and electrophoretic mobility shift assays were performed. Furthermore, HepG2 cells were treated with insulin and high glucose concentrations to induce miRNA expression and IRS2 downregulation.

RESULTS

We found a significant downregulation of IRS2 expression in the liver of obese individuals with type 2 diabetes (0.84 ± 0.08-fold change; p = 0.0833; adjusted p value [p_a] = 0.0417; n = 31) in comparison with non-diabetic obese participants (n = 50). This downregulation correlated with hepatic IRS2 DNA methylation at CpG5. Additionally, CpG6, which is located in intron 1 of IRS2, was hypomethylated in type 2 diabetes; this site spans the sterol regulatory element binding transcription factor 1 (SREBF1) recognition motif, which likely acts as transcriptional repressor. The adjacent polymorphism rs4547213 (G>A) was significantly associated with DNA methylation at a specificity-protein-1 (SP1) binding site (CpG3). Moreover, DNA methylation of cg25924746, a CpG site located in the shore region of the IRS2 promoter-associated CpG island, was increased in the liver of individuals with type 2 diabetes, as compared with those without diabetes. A second epigenetic mechanism, upregulation of hepatic miRNA hsa-let-7e-5p (let-7e-5p) in obese individuals with type 2 diabetes (n = 29) vs non-diabetic obese individuals (n = 49) (1.2 ± 0.08-fold change; p = 0.0332; p_a = 0.0450), is likely to act synergistically with altered IRS2 DNA methylation to decrease IRS2 expression. Mechanistic in vitro experiments demonstrated an acute upregulation of let-7e-5p expression and simultaneous IRS2 downregulation in a liver (HepG2) cell line upon hyperinsulinaemic and hyperglycaemic conditions.

CONCLUSIONS/INTERPRETATION

Our study highlights a new multi-layered epigenetic network that could be involved in subtle dysregulation of IRS2 in the liver of individuals with type 2 diabetes. This might lead to fine-tuning of IRS2 expression and is likely to be supplementary to the already known factors regulating IRS2 expression. Thereby, our findings could support the discovery of new diagnostic and therapeutic strategies for type 2 diabetes. Graphical abstract.

The impact of sex on gene expression across human tissues

Many complex human phenotypes exhibit sex-differentiated characteristics. However, the molecular mechanisms underlying these differences remain largely unknown. We generated a catalog of sex differences in gene expression and in the genetic regulation of gene expression across 44 human tissue sources surveyed by the Genotype-Tissue Expression project (GTEx, v8 release). We demonstrate that sex influences gene expression levels and cellular composition of tissue samples across the human body. A total of 37% of all genes exhibit sex-biased expression in at least one tissue. We identify cis expression quantitative trait loci (eQTLs) with sex-differentiated effects and characterize their cellular origin. By integrating sex-biased eQTLs with genome-wide association study data, we identify 58 gene-trait associations that are driven by genetic regulation of gene expression in a single sex. These findings provide an extensive characterization of sex differences in the human transcriptome and its genetic regulation.

A mega-analysis of expression quantitative trait loci in retinal tissue

Significant association signals from genome-wide association studies (GWAS) point to genomic regions of interest. However, for most loci the causative genetic variant remains undefined. Determining expression quantitative trait loci (eQTL) in a disease relevant tissue is an excellent approach to zoom in on disease- or trait-associated association signals and hitherto on relevant disease mechanisms. To this end, we explored regulation of gene expression in healthy retina (n = 311) and generated the largest cis-eQTL data set available to date. Genotype- and RNA-Seq data underwent rigorous quality control protocols before FastQTL was applied to assess the influence of genetic markers on local (cis) gene expression. Our analysis identified 403,151 significant eQTL variants (eVariants) that regulate 3,007 genes (eGenes) (Q-Value < 0.05). A conditional analysis revealed 744 independent secondary eQTL signals for 598 of the 3,007 eGenes. Interestingly, 99,165 (24.71%) of all unique eVariants regulate the expression of more than one eGene. Filtering the dataset for eVariants regulating three or more eGenes revealed 96 potential regulatory clusters. Of these, 31 harbour 130 genes which are partially regulated by the same genetic signal. To correlate eQTL and association signals, GWAS data from twelve complex eye diseases or traits were included and resulted in identification of 80 eGenes with potential association. Remarkably, expression of 10 genes is regulated by eVariants associated with multiple eye diseases or traits. In conclusion, we generated a unique catalogue of gene expression regulation in healthy retinal tissue and applied this resource to identify potentially pleiotropic effects in highly prevalent human eye diseases. Our study provides an excellent basis to further explore mechanisms of various retinal disease etiologies.

The retina is a multilayered and highly specified neural tissue crucial for high-resolution visual perception and spatial orientation. Environmental and genetic insults to the retina result in many blinding diseases, such as age-related macular degeneration or glaucoma. Commonly, many of these diseases are age-related suggesting that minor changes are accumulating over a life-time, with little or no contribution of strong individual effects. Specifically, this is true for genetic factors known to underlie the etiology of complex diseases including the prevalent eye diseases. In our study, we searched for effects on gene expression due to genetic variation using 311 healthy post-mortem retinal tissue samples. We show that 3,007 of the 16,766 genes investigated are regulated in the retina by genetic variations. Of these, 80 genes are potentially associated to one or more of twelve complex eye diseases or retinal traits tested. Interestingly, 10 genes appear to be involved in the development of several eye traits suggesting that cellular mechanisms may act at a common point in the disease process. Consequently, our study provides the basis to further explore retinal disease pathways and is likely to highlight target molecules for future therapeutic applications.

Genome-wide pQTL analysis of protein expression regulatory networks in the human liver

BACKGROUND

Previous expression quantitative trait loci (eQTL) studies have identified thousands of genetic variants to be associated with gene expression at the mRNA level in the human liver. However, protein expression often correlates poorly with mRNA levels. Thus, protein quantitative trait loci (pQTL) study is required to identify genetic variants that regulate protein expression in human livers.

RESULTS

We conducted a genome-wide pQTL study in 287 normal human liver samples and identified 900 local pQTL variants and 4026 distant pQTL variants. We further discovered 53 genome hotspots of pQTL variants. Transcriptional region mapping analysis showed that 1133 pQTL variants are in transcriptional regulatory regions. Genomic region enrichment analysis of the identified pQTL variants revealed 804 potential regulatory interactions among 595 predicted regulators (e.g., non-coding RNAs) and 394 proteins. Moreover, pQTL variants and trait-variant integration analysis implied several novel mechanisms underlying the relationships between protein expression and liver diseases, such as alcohol dependence. Notably, over 2000 of the identified pQTL variants have not been reported in previous eQTL studies, suggesting extensive involvement of genetic polymorphisms in post-transcriptional regulation of protein expression in human livers.

CONCLUSIONS

We have partially established protein expression regulation networks in human livers and generated a wealth of pQTL data that could serve as a valuable resource for the scientific community.

Expression Quantitative Trait Methylation Analysis Reveals Methylomic Associations With Gene Expression in Childhood Asthma

BACKGROUND

Nasal (airway) epithelial methylation profiles have been associated with asthma, but the effects of such profiles on expression of distant cis-genes are largely unknown.

RESEARCH QUESTION

To identify genes whose expression is associated with proximal and distal CpG probes (within 1 Mb), and to assess whether and how such genes are differentially expressed in atopic asthma.

STUDY DESIGN AND METHODS

Genome-wide expression quantitative trait methylation (eQTM) analysis in nasal epithelium from Puerto Rican subjects (aged 9-20 years) with (n = 219) and without (n = 236) asthma. After the eQTM analysis, a Gene Ontology Enrichment analysis was conducted for the top 500 eQTM genes, and mediation analyses were performed to identify paths from DNA methylation to atopic asthma through gene expression. Asthma was defined as physician-diagnosed asthma and wheeze in the previous year, and atopy was defined as at least one positive IgE to allergens. Atopic asthma was defined as the presence of both atopy and asthma.

RESULTS

We identified 16,867 significant methylation-gene expression pairs (false-discovery rate-adjusted P < .01) in nasal epithelium from study participants. Most eQTM methylation probes were distant (average distance, ∼378 kb) from their target genes, and also more likely to be located in enhancer regions of their target genes in lung tissue than control probes. The top 500 eQTM genes were enriched in pathways for immune processes and epithelial integrity and were more likely to have been previously identified as differentially expressed in atopic asthma. In a mediation analysis, we identified 5,934 paths through which methylation markers could affect atopic asthma through gene expression in nasal epithelium.

INTERPRETATION

Previous epigenome-wide association studies of asthma have estimated the effects of DNA methylation markers on expression of nearby genes in airway epithelium. Our findings suggest that distant epigenetic regulation of gene expression in airway epithelium plays a role in atopic asthma.

Discovery of novel hepatocyte eQTLs in African Americans

African Americans (AAs) are disproportionately affected by metabolic diseases and adverse drug events, with limited publicly available genomic and transcriptomic data to advance the knowledge of the molecular underpinnings or genetic associations to these diseases or drug response phenotypes. To fill this gap, we obtained 60 primary hepatocyte cultures from AA liver donors for genome-wide mapping of expression quantitative trait loci (eQTL) using LAMatrix. We identified 277 eGenes and 19,770 eQTLs, of which 67 eGenes and 7,415 eQTLs are not observed in the Genotype-Tissue Expression Project (GTEx) liver eQTL analysis. Of the eGenes found in GTEx only 25 share the same lead eQTL. These AA-specific eQTLs are less correlated to GTEx eQTLs. in effect sizes and have larger Fst values compared to eQTLs found in both cohorts (overlapping eQTLs). We assessed the overlap between GWAS variants and their tagging variants with AA hepatocyte eQTLs and demonstrated that AA hepatocyte eQTLs can decrease the number of potential causal variants at GWAS loci. Additionally, we identified 75,002 exon QTLs of which 48.8% are not eQTLs in AA hepatocytes. Our analysis provides the first comprehensive characterization of AA hepatocyte eQTLs and highlights the unique discoveries that are made possible due to the increased genetic diversity within the African ancestry genome.

Resources for functional genomic studies of health and development in nonhuman primates

Primates display a wide range of phenotypic variation underlaid by complex genetically regulated mechanisms. The links among DNA sequence, gene function, and phenotype have been of interest from an evolutionary perspective, to understand functional genome evolution and its phenotypic consequences, and from a biomedical perspective to understand the shared and human-specific roots of health and disease. Progress in methods for characterizing genetic, transcriptomic, and DNA methylation (DNAm) variation is driving the rapid development of extensive omics resources, which are now increasingly available from humans as well as a growing number of nonhuman primates (NHPs). The fast growth of large-scale genomic data is driving the emergence of integrated tools and databases, thus facilitating studies of gene functionality across primates. This review describes NHP genomic resources that can aid in exploration of how genes shape primate phenotypes. It focuses on the gene expression trajectories across development in different tissues, the identification of functional genetic variation (including variants deleterious for protein function and regulatory variants modulating gene expression), and DNAm profiles as an emerging tool to understand the process of aging. These resources enable comparative functional genomics approaches to identify species-specific and primate-shared gene functionalities associated with health and development.

A genome-wide multiphenotypic association analysis identified candidate genes and gene ontology shared by four common risky behaviors

BACKGROUND

Risky behaviors can lead to huge economic and health losses. However, limited efforts are paid to explore the genetic mechanisms of risky behaviors.

RESULT

MASH analysis identified a group of target genes for risky behaviors, such as APBB2, MAPT and DCC. For GO enrichment analysis, FUMA detected multiple risky behaviors related GO terms and brain related diseases, such as regulation of neuron differentiation (adjusted P value = 2.84×10^-5), autism spectrum disorder (adjusted P value =1.81×10^-27) and intelligence (adjusted P value =5.89×10^-15).

CONCLUSION

We reported multiple candidate genes and GO terms shared by the four risky behaviors, providing novel clues for understanding the genetic mechanism of risky behaviors.

METHODS

Multivariate Adaptive Shrinkage (MASH) analysis was first applied to the GWAS data of four specific risky behaviors (automobile speeding, drinks per week, ever-smoker, number of sexual partners) to detect the common genetic variants shared by the four risky behaviors. Utilizing genomic functional annotation data of SNPs, the SNPs detected by MASH were then mapped to target genes. Finally, gene set enrichment analysis of the identified candidate genes were conducted by the FUMA platform to obtain risky behaviors related gene ontology (GO) terms as well as diseases and traits, respectively.

Increased circulating levels of Factor H-Related Protein 4 are strongly associated with age-related macular degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness. Genetic variants at the chromosome 1q31.3 encompassing the complement factor H (CFH, FH) and CFH related genes (CFHR1-5) are major determinants of AMD susceptibility, but their molecular consequences remain unclear. Here we demonstrate that FHR-4 plays a prominent role in AMD pathogenesis. We show that systemic FHR-4 levels are elevated in AMD (P-value = 7.1 × 10-6), whereas no difference is seen for FH. Furthermore, FHR-4 accumulates in the choriocapillaris, Bruch's membrane and drusen, and can compete with FH/FHL-1 for C3b binding, preventing FI-mediated C3b cleavage. Critically, the protective allele of the strongest AMD-associated CFH locus variant rs10922109 has the highest association with reduced FHR-4 levels (P-value = 2.2 × 10-56), independently of the AMD-protective CFHR1-3 deletion, and even in those individuals that carry the high-risk allele of rs1061170 (Y402H). Our findings identify FHR-4 as a key molecular player contributing to complement dysregulation in AMD.

A transcriptome-wide association study based on 27 tissues identifies 106 genes potentially relevant for disease pathology in age-related macular degeneration

Genome-wide association studies (GWAS) for late stage age-related macular degeneration (AMD) have identified 52 independent genetic variants with genome-wide significance at 34 genomic loci. Typically, such an approach rarely results in the identification of functional variants implicating a defined gene in the disease process. We now performed a transcriptome-wide association study (TWAS) allowing the prediction of effects of AMD-associated genetic variants on gene expression. The TWAS was based on the genotypes of 16,144 late-stage AMD cases and 17,832 healthy controls, and gene expression was imputed for 27 different human tissues which were obtained from 134 to 421 individuals. A linear regression model including each individuals imputed gene expression data and the respective AMD status identified 106 genes significantly associated to AMD variants in at least one tissue (Q-value < 0.001). Gene enrichment analysis highlighted rather systemic than tissue- or cell-specific processes. Remarkably, 31 of the 106 genes overlapped with significant GWAS signals of other complex traits and diseases, such as neurological or autoimmune conditions. Taken together, our study highlights the fact that expression of genes associated with AMD is not restricted to retinal tissue as could be expected for an eye disease of the posterior pole, but instead is rather ubiquitous suggesting processes underlying AMD pathology to be of systemic nature.

A New Liver Expression Quantitative Trait Locus Map From 1,183 Individuals Provides Evidence for Novel Expression Quantitative Trait Loci of Drug Response, Metabolic, and Sex-Biased Phenotypes

Expression quantitative trait locus (eQTL) studies in human liver are crucial for elucidating how genetic variation influences variability in disease risk and therapeutic outcomes and may help guide strategies to obtain maximal efficacy and safety of clinical interventions. Associations between expression microarray and genome-wide genotype data from four human liver eQTL studies (n = 1,183) were analyzed. More than 2.3 million cis-eQTLs for 15,668 genes were identified. When eQTLs were filtered against a list of 1,496 drug response genes, 187,829 cis-eQTLs for 1,191 genes were identified. Additionally, 1,683 sex-biased cis-eQTLs were identified, as well as 49 and 73 cis-eQTLs that colocalized with genome-wide association study signals for blood metabolite or lipid levels, respectively. Translational relevance of these results is evidenced by linking DPYD eQTLs to differences in safety of chemotherapy, linking the sex-biased regulation of PCSK9 expression to anti-lipid therapy, and identifying the G-protein coupled receptor GPR180 as a novel drug target for hypertriglyceridemia.

Allelic Heterogeneity at the CRP Locus Identified by Whole-Genome Sequencing in Multi-ancestry Cohorts

Whole-genome sequencing (WGS) can improve assessment of low-frequency and rare variants, particularly in non-European populations that have been underrepresented in existing genomic studies. The genetic determinants of C-reactive protein (CRP), a biomarker of chronic inflammation, have been extensively studied, with existing genome-wide association studies (GWASs) conducted in >200,000 individuals of European ancestry. In order to discover novel loci associated with CRP levels, we examined a multi-ancestry population (n = 23,279) with WGS (∼38× coverage) from the Trans-Omics for Precision Medicine (TOPMed) program. We found evidence for eight distinct associations at the CRP locus, including two variants that have not been identified previously (rs11265259 and rs181704186), both of which are non-coding and more common in individuals of African ancestry (∼10% and ∼1% minor allele frequency, respectively, and rare or monomorphic in 1000 Genomes populations of East Asian, South Asian, and European ancestry). We show that the minor (G) allele of rs181704186 is associated with lower CRP levels and decreased transcriptional activity and protein binding in vitro, providing a plausible molecular mechanism for this African ancestry-specific signal. The individuals homozygous for rs181704186-G have a mean CRP level of 0.23 mg/L, in contrast to individuals heterozygous for rs181704186 with mean CRP of 2.97 mg/L and major allele homozygotes with mean CRP of 4.11 mg/L. This study demonstrates the utility of WGS in multi-ethnic populations to drive discovery of complex trait associations of large effect and to identify functional alleles in noncoding regulatory regions.

Investigating skewness to understand gene expression heterogeneity in large patient cohorts

BACKGROUND

Skewness is an under-utilized statistical measure that captures the degree of asymmetry in the distribution of any dataset. This study applied a new metric based on skewness to identify regulators or genes that have outlier expression in large patient cohorts.

RESULTS

We investigated whether specific patterns of skewed expression were related to the enrichment of biological pathways or genomic properties like DNA methylation status. Our study used publicly available datasets that were generated using both RNA-sequencing and microarray technology platforms. For comparison, the datasets selected for this study also included different samples derived from control donors and cancer patients. When comparing the shift in expression skewness between cancer and control datasets, we observed an enrichment of pathways related to the immune function that reflects an increase towards positive skewness in the cancer relative to control datasets. A significant correlation was also detected between expression skewness and the top 500 genes corresponding to the most significant differential DNA methylation occurring in the promotor regions for four Cancer Genome Atlas cancer cohorts.

CONCLUSIONS

Our results indicate that expression skewness can reveal new insights into transcription based on outlier and asymmetrical behaviour in large patient cohorts.

Genetic Control of Expression and Splicing in Developing Human Brain Informs Disease Mechanisms

Tissue-specific regulatory regions harbor substantial genetic risk for disease. Because brain development is a critical epoch for neuropsychiatric disease susceptibility, we characterized the genetic control of the transcriptome in 201 mid-gestational human brains, identifying 7,962 expression quantitative trait loci (eQTL) and 4,635 spliceQTL (sQTL), including several thousand prenatal-specific regulatory regions. We show that significant genetic liability for neuropsychiatric disease lies within prenatal eQTL and sQTL. Integration of eQTL and sQTL with genome-wide association studies (GWAS) via transcriptome-wide association identified dozens of novel candidate risk genes, highlighting shared and stage-specific mechanisms in schizophrenia (SCZ). Gene network analysis revealed that SCZ and autism spectrum disorder (ASD) affect distinct developmental gene co-expression modules. Yet, in each disorder, common and rare genetic variation converges within modules, which in ASD implicates superficial cortical neurons. More broadly, these data, available as a web browser and our analyses, demonstrate the genetic mechanisms by which developmental events have a widespread influence on adult anatomical and behavioral phenotypes.

Biological characterization of expression quantitative trait loci (eQTLs) showing tissue-specific opposite directional effects

Interpreting the susceptible loci documented by genome-wide association studies (GWASs) is of utmost importance in the post-GWAS era. Since most complex traits are contributed by multiple tissues, analyzing tissue-specific effects of expression quantitative trait loci (eQTLs) is a promising approach. Here we describe “opposite eQTL effects”, i.e., gene expression effects of eQTLs that are in the opposite direction between different tissues, as the biologically meaningful annotations of genes and genetic variants for understanding the GWAS loci. The genes and single-nucleotide polymorphisms (SNPs) associated with the opposite eQTL effects (opp-multi-eQTL-Genes and opp-multi-eQTL-SNPs) were extracted from the largest eQTL database provided by the Genotype-Tissue Expression (GTEx) project (release version 7). The opposite eQTL effects were detected even between closely related tissues such as cerebellum and brain cortex, and a significant proportion of the genes having eQTLs were annotated as the opp-multi-eQTL-Genes (2,323 out of 31,212; 7.4%). The opp-multi-eQTL-SNPs showed locational enrichment at the transcription start site and also possible involvement of epigenetic regulation. The biological importance of the opposite eQTL effects was also assessed using the SNPs reported in GWASs (GWAS-SNPs), which demonstrated that a high proportion of the opp-multi-eQTL-SNPs are in linkage disequilibrium with the GWAS-SNPs (2,498 out of 9,290; 26.9%). Based on the results, the opposite eQTL effects can be a common phenomenon in the tissue-specific gene regulation with a possible contribution to the development of complex traits.

Retinal transcriptome and eQTL analyses identify genes associated with age-related macular degeneration

Genome-wide association studies (GWAS) have identified genetic variants at 34 loci contributing to age-related macular degeneration (AMD)1-3. We generated transcriptional profiles of postmortem retinas from 453 controls and cases at distinct stages of AMD and integrated retinal transcriptomes, covering 13,662 protein-coding and 1,462 noncoding genes, with genotypes at more than 9 million common SNPs for expression quantitative trait loci (eQTL) analysis of a tissue not included in Genotype-Tissue Expression (GTEx) and other large datasets4,5. Cis-eQTL analysis identified 10,474 genes under genetic regulation, including 4,541 eQTLs detected only in the retina. Integrated analysis of AMD-GWAS with eQTLs ascertained likely target genes at six reported loci. Using transcriptome-wide association analysis (TWAS), we identified three additional genes, RLBP1, HIC1 and PARP12, after Bonferroni correction. Our studies expand the genetic landscape of AMD and establish the Eye Genotype Expression (EyeGEx) database as a resource for post-GWAS interpretation of multifactorial ocular traits.