Dissecting the genetics of the human transcriptome identifies novel trait-related trans-eQTLs and corroborates the regulatory relevance of non-protein coding loci†

Sex and statin-related genetic associations at the PCSK9 gene locus: results of genome-wide association meta-analysis

BACKGROUND

Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a key player of lipid metabolism with higher plasma levels in women throughout their life. Statin treatment affects PCSK9 levels also showing evidence of sex-differential effects. It remains unclear whether these differences can be explained by genetics.

METHODS

We performed genome-wide association meta-analyses (GWAS) of PCSK9 levels stratified for sex and statin treatment in six independent studies of Europeans (8936 women/11,080 men respectively 14,825 statin-free/5191 statin-treated individuals). Loci associated in one of the strata were tested for statin- and sex-interactions considering all independent signals per locus. Independent variants at the PCSK9 gene locus were then used in a stratified Mendelian Randomization analysis (cis-MR) of PCSK9 effects on low-density lipoprotein cholesterol (LDL-C) levels to detect differences of causal effects between the subgroups.

RESULTS

We identified 11 loci associated with PCSK9 in at least one stratified subgroup (p < 1.0 × 10-6), including the PCSK9 gene locus and five other lipid loci: APOB, TM6SF2, FADS1/FADS2, JMJD1C, and HP/HPR. The interaction analysis revealed eight loci with sex- and/or statin-interactions. At the PCSK9 gene locus, there were four independent signals, one with a significant sex-interaction showing stronger effects in men (rs693668). Regarding statin treatment, there were two significant interactions in PCSK9 missense mutations: rs11591147 had stronger effects in statin-free individuals, and rs11583680 had stronger effects in statin-treated individuals. Besides replicating known loci, we detected two novel genome-wide significant associations: one for statin-treated individuals at 6q11.1 (within KHDRBS2) and one for males at 12q24.22 (near KSR2/NOS1), both with significant interactions. In the MR of PCSK9 on LDL-C, we observed significant causal estimates within all subgroups, but significantly stronger causal effects in statin-free subjects compared to statin-treated individuals.

CONCLUSIONS

We performed the first double-stratified GWAS of PCSK9 levels and identified multiple biologically plausible loci with genetic interaction effects. Our results indicate that the observed sexual dimorphism of PCSK9 and its statin-related interactions have a genetic basis. Significant differences in the causal relationship between PCSK9 and LDL-C suggest sex-specific dosages of PCSK9 inhibitors.

An atlas of genome-wide gene expression and metabolite associations and possible mediation effects towards body mass index

Investigating the cross talk of different omics layers is crucial to understand molecular pathomechanisms of metabolic diseases like obesity. Here, we present a large-scale association meta-analysis of genome-wide whole blood and peripheral blood mononuclear cell (PBMC) gene expressions profiled with Illumina HT12v4 microarrays and metabolite measurements from dried blood spots (DBS) characterized by targeted liquid chromatography tandem mass spectrometry (LC-MS/MS) in three large German cohort studies with up to 7706 samples. We found 37,295 associations comprising 72 amino acids (AA) and acylcarnitine (AC) metabolites (including ratios) and 8579 transcripts. We applied this catalogue of associations to investigate the impact of associating transcript-metabolite pairs on body mass index (BMI) as an example metabolic trait. This is achieved by conducting a comprehensive mediation analysis considering metabolites as mediators of gene expression effects and vice versa. We discovered large mediation networks comprising 27,023 potential mediation effects within 20,507 transcript-metabolite pairs. Resulting networks of highly connected (hub) transcripts and metabolites were leveraged to gain mechanistic insights into metabolic signaling pathways. In conclusion, here, we present the largest available multi-omics integration of genome-wide transcriptome data and metabolite data of amino acid and fatty acid metabolism and further leverage these findings to characterize potential mediation effects towards BMI proposing candidate mechanisms of obesity and related metabolic diseases. KEY MESSAGES: Thousands of associations of 72 amino acid and acylcarnitine metabolites and 8579 genes expand the knowledge of metabolome-transcriptome associations. A mediation analysis of effects on body mass index revealed large mediation networks of thousands of obesity-related gene-metabolite pairs. Highly connected, potentially mediating hub genes and metabolites enabled insight into obesity and related metabolic disease pathomechanisms.

Genetic pathways regulating the longitudinal acquisition of cocaine self-administration in a panel of inbred and recombinant inbred mice

To identify addiction genes, we evaluate intravenous self-administration of cocaine or saline in 84 inbred and recombinant inbred mouse strains over 10 days. We integrate the behavior data with brain RNA-seq data from 41 strains. The self-administration of cocaine and that of saline are genetically distinct. We maximize power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci are identified in the genome-wide association study. A transcriptome-wide association study highlights the Trpv2 ion channel as a key locus for cocaine self-administration as well as identifying 17 additional genes, including Arhgef26, Slc18b1, and Slco5a1. We find numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance of Trpv2, an ionotropic cannabinoid receptor, for the response to cocaine.

Genome-wide aggregated trans-effects on risk of type 1 diabetes: A test of the "omnigenic" sparse effector hypothesis of complex trait genetics

The "omnigenic" hypothesis postulates that the polygenic effects of common SNPs on a typical complex trait are mediated through trans-effects on expression of a relatively sparse set of effector ("core") genes. We tested this hypothesis in a study of 4,964 cases of type 1 diabetes (T1D) and 7,497 controls by using summary statistics to calculate aggregated (excluding the HLA region) trans-scores for gene expression in blood. From associations of T1D with aggregated trans-scores, nine putative core genes were identified, of which three-STAT1, CTLA4 and FOXP3-are genes in which variants cause monogenic forms of autoimmune diabetes. Seven of these genes affect the activity of regulatory T cells, and two are involved in immune responses to microbial lipids. Four T1D-associated genomic regions could be identified as master regulators via trans-effects on gene expression. These results support the sparse effector hypothesis and reshape our understanding of the genetic architecture of T1D.

Whole genome DNA and RNA sequencing of whole blood elucidates the genetic architecture of gene expression underlying a wide range of diseases

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis-eQTL variant-gene transcript (eGene) pairs at p < 5 × 10^-8 (2,855,111 unique cis-eQTL variants and 15,982 unique eGenes) and 1,469,754 trans-eQTL variant-eGene pairs at p < 1e-12 (526,056 unique trans-eQTL variants and 7233 unique eGenes). In addition, 442,379 cis-eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis-eGenes are enriched for immune functions (FDR < 0.05). The cis-eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

Genotype-phenotype correlation of T-cell subtypes reveals senescent and cytotoxic genes in Alzheimer's disease

Recent studies identifying expression quantitative trait loci (eQTLs) in immune cells have uncovered important links between disease risk alleles and gene expression trends in monocytes, T cells and other cell types. However, these studies are generally done with young, healthy subjects, limiting the utility of their findings for age-related conditions such as Alzheimer's disease (AD). We have performed RNA sequencing on four T-cell subsets in genome-wide genotyped and well-characterized AD subjects and age- and sex-matched controls from the Religious Orders Study/Memory and Aging Project. We correlated gene expression data with AD neuropathological traits and with single-nucleotide polymorphisms to detect eQTLs. We identified several significant genes involved in T-cell senescence and cytotoxicity, consistent with T-cell RNA sequencing studies in aged/AD cohorts. We identified unexpected eQTLs previously associated with neuropsychiatric disease traits. Finally, we discovered that pathways related to axon guidance and synaptic function were enriched among trans-eQTLs in coding regions of the genome. Our data strengthen the potential link between T-cell senescence and age-related neurodegenerative disease. In addition, our eQTL data suggest that T-cell phenotypes may influence neuropsychiatric disorders and can be influenced by genes involved in neurodevelopmental processes.

DNA methylation patterns reflect individual's lifestyle independent of obesity

OBJECTIVE

Obesity is driven by modifiable lifestyle factors whose effects may be mediated by epigenetics. Therefore, we investigated lifestyle effects on blood DNA methylation in participants of the LIFE-Adult study, a well-characterised population-based cohort from Germany.

RESEARCH DESIGN AND METHODS

Lifestyle scores (LS) based on diet, physical activity, smoking and alcohol intake were calculated in 4107 participants of the LIFE-Adult study. Fifty subjects with an extremely healthy lifestyle and 50 with an extremely unhealthy lifestyle (5th and 95th percentiles LS) were selected for genome-wide DNA methylation analysis in blood samples employing Illumina Infinium® Methylation EPIC BeadChip system technology.

RESULTS

Differences in DNA methylation patterns between body mass index groups (<25 vs. >30 kg/m² ) were rather marginal compared to inter-lifestyle differences (0 vs. 145 differentially methylated positions [DMPs]), which identified 4682 differentially methylated regions (DMRs; false discovery rate [FDR <5%) annotated to 4426 unique genes. A DMR annotated to the glutamine-fructose-6-phosphate transaminase 2 (GFPT2) locus showed the strongest hypomethylation (∼6.9%), and one annotated to glutamate rich 1 (ERICH1) showed the strongest hypermethylation (∼5.4%) in healthy compared to unhealthy lifestyle individuals. Intersection analysis showed that diet, physical activity, smoking and alcohol intake equally contributed to the observed differences, which affected, among others, pathways related to glutamatergic synapses (adj. p < .01) and axon guidance (adj. p < .05). We showed that methylation age correlates with chronological age and waist-to-hip ratio with lower DNA methylation age (DNAmAge) acceleration distances in participants with healthy lifestyles. Finally, two identified top DMPs for the alanyl aminopeptidase (ANPEP) locus also showed the strongest expression quantitative trait methylation in blood.

CONCLUSIONS

DNA methylation patterns help discriminate individuals with a healthy versus unhealthy lifestyle, which may mask subtle methylation differences derived from obesity.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis -eQTL variant-gene transcript (eGene) pairs at p  < 5x10 - 8 (2,855,111 unique cis -eQTL variants and 15,982 unique eGenes) and 1,469,754 trans -eQTL variant-eGene pairs at p  < 1e-12 (526,056 unique trans -eQTL variants and 7,233 unique eGenes). In addition, 442,379 cis -eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis- eGenes are enriched for immune functions (FDR < 0.05). The cis -eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

Genetically regulated gene expression and proteins revealed discordant effects

Background

Although gene-expression (GE) and protein levels are typically strongly genetically regulated, their correlation is known to be low. Here we investigate this phenomenon by focusing on the genetic background of this correlation in order to understand the similarities and differences in the genetic regulation of these omics layers.

Methods and results

We performed locus-wide association studies of 92 protein levels measured in whole blood for 2,014 samples of European ancestry and found that 66 are genetically regulated. Three female- and one male-specific effects were detected. We estimated the genetically regulated GE for all significant genes in 49 GTEx v8 tissues. A total of 7 proteins showed negative correlations with their respective GE across multiple tissues. Finally, we tested for causal links of GE on protein expression via Mendelian Randomization, and confirmed a negative causal effect of GE on protein level for five of these genes in a total of 63 gene-tissue pairs: BLMH, CASP3, CXCL16, IL6R, and SFTPD. For IL6R, we replicated the negative causal effect on coronary-artery disease (CAD), while its GE was positively linked to CAD.

Conclusion

While total GE and protein levels are only weakly correlated, we found high correlations between their genetically regulated components across multiple tissues. Of note, strong negative causal effects of tissue-specific GE on five protein levels were detected. Causal network analyses revealed that GE effects on CAD risks was in general mediated by protein levels.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis -eQTL variant-gene transcript (eGene) pairs at p <5×10 -8 (2,855,111 unique cis -eQTL variants and 15,982 unique eGenes) and 1,469,754 trans -eQTL variant-eGene pairs at p <1e-12 (526,056 unique trans -eQTL variants and 7,233 unique eGenes). In addition, 442,379 cis -eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis- eGenes are enriched for immune functions (FDR <0.05). The cis -eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

First genome-wide association study of 99 body measures derived from 3-dimensional body scans

Body height, body mass index, hip and waist circumference are important risk factors or outcome variables in clinical and epidemiological research with complex underlying genetics. However, these classical anthropometric traits represent only a very limited view on the human body and other traits with potentially higher functional specificity are not yet studied to a larger extent. Participants of LIFE-Adult were assessed by three-dimensional body scanner VITUS XXL determining 99 high-quality anthropometric traits in parallel. Genotyping was performed by Axiom Genome-Wide CEU 1 Array Plate microarray technology and imputation was done using 1000 Genomes phase 3 reference panel. Combined phenotype and genetic information are available for a total of 7,562 participants. Largest heritabilities were estimated for height traits (maximum heritability with h² = 44% for neck height) and 61 traits achieved values larger than 20%. By genome-wide analyses, we identified 16 loci associated with at least one of the 99 traits. Ten of these loci were not described for association with classical anthropometric traits so far. The strongest novel association was observed for 7p14.3 (rs11979006, P = 2.12 × 10⁻⁹) for the trait Back Width with ZNRF2 as the most plausible candidate gene. Loci established for association with classical anthropometric traits were subjected to anthropometric phenome-wide association analysis. From the reported 709 loci, 211 are co-associated with body scanner traits (enrichment: OR = 1.96, P = 1.08 × 10⁻⁶¹). We conclude that genetics of 3D laser-based anthropometry is promising to identify novel loci and to improve the functional understanding of established ones.

Transcriptome Analyses of Adipose Tissue Samples Identify EGFL6 as a Candidate Gene Involved in Obesity-Related Adipose Tissue Dysfunction in Children

Obesity develops early in childhood and is accompanied by early signs of adipose tissue (AT) dysfunction and metabolic disease in children. In order to analyse the molecular processes during obesity-related AT accumulation in children, we investigated genome-wide expression profiles in AT samples, isolated adipocytes, and stromal vascular fraction (SVF) cells and assessed their relation to obesity as well as biological and functional AT parameters. We detected alterations in gene expression associated with obesity and related parameters, i.e., BMI SDS, adipocyte size, macrophage infiltration, adiponectin, and/or leptin. While differential gene expression in AT and adipocytes shared an enrichment in metabolic pathways and pathways related to extracellular structural organisation, SVF cells showed an overrepresentation in inflammatory pathways. In adipocytes, we found the strongest positive association for epidermal growth factor-like protein 6 (EGFL6) with adipocyte hypertrophy. EGFL6 was also upregulated during in vitro adipocyte differentiation. In children, EGFL6 expression was positively correlated to parameters of AT dysfunction and metabolic disease such as macrophage infiltration into AT, hs-CRP, leptin levels, and HOMA-IR. In conclusion, we provide evidence for early alterations in AT gene expression related to AT dysfunction in children and identified EGFL6 as potentially being involved in processes underlying the pathogenesis of metabolic disease.

Whole Blood Metabolite Profiles Reflect Changes in Energy Metabolism in Heart Failure

A variety of atherosclerosis and cardiovascular disease (ASCVD) phenotypes are tightly linked to changes in the cardiac energy metabolism that can lead to a loss of metabolic flexibility and to unfavorable clinical outcomes. We conducted an association analysis of 31 ASCVD phenotypes and 97 whole blood amino acids, acylcarnitines and derived ratios in the LIFE-Adult (n = 9646) and LIFE-Heart (n = 5860) studies, respectively. In addition to hundreds of significant associations, a total of 62 associations of six phenotypes were found in both studies. Positive associations of various amino acids and a range of acylcarnitines with decreasing cardiovascular health indicate disruptions in mitochondrial, as well as peroxisomal fatty acid oxidation. We complemented our metabolite association analyses with whole blood and peripheral blood mononuclear cell (PBMC) gene-expression analyses of fatty acid oxidation and ketone-body metabolism related genes. This revealed several differential expressions for the heart failure biomarker N-terminal prohormone of brain natriuretic peptide (NT-proBNP) in peripheral blood mononuclear cell (PBMC) gene expression. Finally, we constructed and compared three prediction models of significant stenosis in the LIFE-Heart study using (1) traditional risk factors only, (2) the metabolite panel only and (3) a combined model. Area under the receiver operating characteristic curve (AUC) comparison of these three models shows an improved prediction accuracy for the combined metabolite and classical risk factor model (AUC = 0.78, 95%-CI: 0.76–0.80). In conclusion, we improved our understanding of metabolic implications of ASCVD phenotypes by observing associations with metabolite concentrations and gene expression of the mitochondrial and peroxisomal fatty acid oxidation. Additionally, we demonstrated the predictive potential of the metabolite profile to improve classification of patients with significant stenosis.

Human genetic diversity regulating the TLR10/TLR1/TLR6 locus confers increased cytokines in response to Chlamydia trachomatis

Human genetic diversity can have profound effects on health outcomes upon exposure to infectious agents. For infections with Chlamydia trachomatis (C. trachomatis), the wide range of genital and ocular disease manifestations are likely influenced by human genetic differences that regulate interactions between C. trachomatis and host cells. We leveraged this diversity in cellular responses to demonstrate the importance of variation at the Toll-like receptor 1 (TLR1), TLR6, and TLR10 locus to cytokine production in response to C. trachomatis. We determined that a single-nucleotide polymorphism (SNP) (rs1057807), located in a region that forms a loop with the TLR6 promoter, is associated with increased expression of TLR1, TLR6, and TLR10 and secreted levels of ten C. trachomatis-induced cytokines. Production of these C. trachomatis-induced cytokines is primarily dependent on MyD88 and TLR6 based on experiments using inhibitors, blocking antibodies, RNAi, and protein overexpression. Population genetic analyses further demonstrated that the mean IL-6 response of cells from two European populations were higher than the mean response of cells from three African populations and that this difference was partially attributable to variation in rs1057807 allele frequency. In contrast, a SNP associated with a different pro-inflammatory cytokine (rs2869462 associated with the chemokine CXCL10) exhibited an opposite response, underscoring the complexity of how different genetic variants contribute to an individual's immune response. This multidisciplinary study has identified a long-range chromatin interaction and genetic variation that regulates TLR6 to broaden our understanding of how human genetic variation affects the C. trachomatis-induced immune response.

Genome-wide meta-analysis of phytosterols reveals five novel loci and a detrimental effect on coronary atherosclerosis

Phytosterol serum concentrations are under tight genetic control. The relationship between phytosterols and coronary artery disease (CAD) is controversially discussed. We perform a genome-wide meta-analysis of 32 phytosterol traits reflecting resorption, cholesterol synthesis and esterification in six studies with up to 9758 subjects and detect ten independent genome-wide significant SNPs at seven genomic loci. We confirm previously established associations at ABCG5/8 and ABO and demonstrate an extended locus heterogeneity at ABCG5/8 with different functional mechanisms. New loci comprise HMGCR, NPC1L1, PNLIPRP2, SCARB1 and APOE. Based on these results, we perform Mendelian Randomization analyses (MR) revealing a risk-increasing causal relationship of sitosterol serum concentrations and CAD, which is partly mediated by cholesterol. Here we report that phytosterols are polygenic traits. MR add evidence of both, direct and indirect causal effects of sitosterol on CAD.

Genetic Regulation of Cytokine Response in Patients with Acute Community-Acquired Pneumonia

Background: Community-acquired pneumonia (CAP) is an acute disease condition with a high risk of rapid deteriorations. We analysed the influence of genetics on cytokine regulation to obtain a better understanding of patient’s heterogeneity. Methods: For up to N = 389 genotyped participants of the PROGRESS study of hospitalised CAP patients, we performed a genome-wide association study of ten cytokines IL-1β, IL-6, IL-8, IL-10, IL-12, MCP-1 (MCAF), MIP-1α (CCL3), VEGF, VCAM-1, and ICAM-1. Consecutive secondary analyses were performed to identify independent hits and corresponding causal variants. Results: 102 SNPs from 14 loci showed genome-wide significant associations with five of the cytokines. The most interesting associations were found at 6p21.1 for VEGF (p = 1.58 × 10⁻²⁰), at 17q21.32 (p = 1.51 × 10⁻⁹) and at 10p12.1 (p = 2.76 × 10⁻⁹) for IL-1β, at 10p13 for MIP-1α (CCL3) (p = 2.28 × 10⁻⁹), and at 9q34.12 for IL-10 (p = 4.52 × 10⁻⁸). Functionally plausible genes could be assigned to the majority of loci including genes involved in cytokine secretion, granulocyte function, and cilial kinetics. Conclusion: This is the first context-specific genetic association study of blood cytokine concentrations in CAP patients revealing numerous biologically plausible candidate genes. Two of the loci were also associated with atherosclerosis with probable common or consecutive pathomechanisms.

Genetic Regulation of Transcription in the Endometrium in Health and Disease

The endometrium is a complex and dynamic tissue essential for fertility and implicated in many reproductive disorders. The tissue consists of glandular epithelium and vascularised stroma and is unique because it is constantly shed and regrown with each menstrual cycle, generating up to 10 mm of new mucosa. Consequently, there are marked changes in cell composition and gene expression across the menstrual cycle. Recent evidence shows expression of many genes is influenced by genetic variation between individuals. We and others have reported evidence for genetic effects on hundreds of genes in endometrium. The genetic factors influencing endometrial gene expression are highly correlated with the genetic effects on expression in other reproductive (e.g., in uterus and ovary) and digestive tissues (e.g., salivary gland and stomach), supporting a shared genetic regulation of gene expression in biologically similar tissues. There is also increasing evidence for cell specific genetic effects for some genes. Sample size for studies in endometrium are modest and results from the larger studies of gene expression in blood report genetic effects for a much higher proportion of genes than currently reported for endometrium. There is also emerging evidence for the importance of genetic variation on RNA splicing. Gene mapping studies for common disease, including diseases associated with endometrium, show most variation maps to intergenic regulatory regions. It is likely that genetic risk factors for disease function through modifying the program of cell specific gene expression. The emerging evidence from our gene mapping studies coupled with tissue specific studies, and the GTEx, eQTLGen and EpiMap projects, show we need to expand our understanding of the complex regulation of gene expression. These data also help to link disease genetic risk factors to specific target genes. Combining our data on genetic regulation of gene expression in endometrium, and cell types within the endometrium with gene mapping data for endometriosis and related diseases is beginning to uncover the specific genes and pathways responsible for increased risk of these diseases.

eQTLs as causal instruments for the reconstruction of hormone linked gene networks

Hormones act within in highly dynamic systems and much of the phenotypic response to variation in hormone levels is mediated by changes in gene expression. The increase in the number and power of large genetic association studies has led to the identification of hormone linked genetic variants. However, the biological mechanisms underpinning the majority of these loci are poorly understood. The advent of affordable, high throughput next generation sequencing and readily available transcriptomic databases has shown that many of these genetic variants also associate with variation in gene expression levels as expression Quantitative Trait Loci (eQTLs). In addition to further dissecting complex genetic variation, eQTLs have been applied as tools for causal inference. Many hormone networks are driven by transcription factors, and many of these genes can be linked to eQTLs. In this mini-review, we demonstrate how causal inference and gene networks can be used to describe the impact of hormone linked genetic variation upon the transcriptome within an endocrinology context.

Identification of Disease-Related Genes That Are Common between Alzheimer's and Cardiovascular Disease Using Blood Genome-Wide Transcriptome Analysis

Accumulating evidence has suggested a shared pathophysiology between Alzheimer’s disease (AD) and cardiovascular disease (CVD). Based on genome-wide transcriptomes, specifically those of blood samples, we identify the shared disease-related signatures between AD and CVD. In addition to gene expressions in blood, the following prior knowledge were utilized to identify several candidate disease-related gene (DRG) sets: protein–protein interactions, transcription factors, disease–gene relationship databases, and single nucleotide polymorphisms. We selected the respective DRG sets for AD and CVD that show a high accuracy for disease prediction in bulk and single-cell gene expression datasets. Then, gene regulatory networks (GRNs) were constructed from each of the AD and CVD DRG sets to identify the upstream regulating genes. Using the GRNs, we identified two common upstream genes (GPBP1 and SETDB2) between the AD and CVD GRNs. In summary, this study has identified the potential AD- and CVD-related genes and common hub genes between these sets, which may help to elucidate the shared mechanisms between these two diseases.

Meta-GWAS of PCSK9 levels detects two novel loci at APOB and TM6SF2

Background

Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a key player in lipid metabolism, as it degrades low-density lipoprotein (LDL) receptors from hepatic cell membranes. So far, only variants of the PCSK9 gene locus were found to be associated with PCSK9 levels. Here we aimed to identify novel genetic loci that regulate PCSK9 levels and how they relate to other lipid traits. Additionally, we investigated to what extend the causal effect of PCSK9 on coronary artery disease (CAD) is mediated by low-density lipoprotein–cholesterol (LDL–C).

Methods and Results

We performed a genome-wide association study meta-analysis of PCSK9 levels in up to 12 721 samples of European ancestry. The estimated heritability was 10.3%, which increased to 12.6% using only samples from patients without statin treatment. We successfully replicated the known PCSK9 hit consisting of three independent signals. Interestingly, in a study of 300 African Americans, we confirmed the locus with a different PCSK9 variant. Beyond PCSK9, our meta-analysis detected three novel loci with genome-wide significance. Co-localization analysis with cis-eQTLs and lipid traits revealed biologically plausible candidate genes at two of them: APOB and TM6SF2. In a bivariate Mendelian Randomization analysis, we detected a strong effect of PCSK9 on LDL-C, but not vice versa. LDL-C mediated 63% of the total causal effect of PCSK9 on CAD.

Conclusion

Our study identified novel genetic loci with plausible candidate genes affecting PCSK9 levels. Ethnic heterogeneity was observed at the PCSK9 locus itself. Although the causal effect of PCSK9 on CAD is mainly mediated by LDL-C, an independent direct effect also occurs.

Joint eQTL mapping and inference of gene regulatory network improves power of detecting both cis- and trans-eQTLs

MOTIVATION

Genetic variations of expression quantitative trait loci (eQTLs) play a critical role in influencing complex traits and diseases development. Two main factors that affect the statistical power of detecting eQTLs are: (i) relatively small size of samples available, and (ii) heavy burden of multiple testing due to a very large number of variants to be tested. The later issue is particularly severe when one tries to identify trans-eQTLs that are far away from the genes they influence. If one can exploit co-expressed genes jointly in eQTL-mapping, effective sample size can be increased. Furthermore, using the structure of the gene regulatory network (GRN) may help to identify trans-eQTLs without increasing multiple testing burden.

RESULTS

In this article, we use the structure equation model (SEM) to model both GRN and effect of eQTLs on gene expression, and then develop a novel algorithm, named sparse SEM for eQTL mapping (SSEMQ), to conduct joint eQTL mapping and GRN inference. The SEM can exploit co-expressed genes jointly in eQTL mapping and also use GRN to determine trans-eQTLs. Computer simulations demonstrate that our SSEMQ significantly outperforms nine existing eQTL mapping methods. SSEMQ is further used to analyze two real datasets of human breast and whole blood tissues, yielding a number of cis- and trans-eQTLs.

AVAILABILITY AND IMPLEMENTATION

R package ssemQr is available at https://github.com/Ivis4ml/ssemQr.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.

From DNA human sequence to the chromatin higher order organisation and its biological meaning: Using biomolecular interaction networks to understand the influence of structural variation on spatial genome organisation and its functional effect

The three-dimensional structure of the human genome has been proven to have a significant functional impact on gene expression. The high-order spatial chromatin is organised first by looping mediated by multiple protein factors, and then it is further formed into larger structures of topologically associated domains (TADs) or chromatin contact domains (CCDs), followed by A/B compartments and finally the chromosomal territories (CTs). The genetic variation observed in human population influences the multi-scale structures, posing a question regarding the functional impact of structural variants reflected by the variability of the genes expression patterns. The current methods of evaluating the functional effect include eQTLs analysis which uses statistical testing of influence of variants on spatially close genes. Rarely, non-coding DNA sequence changes are evaluated by their impact on the biomolecular interaction network (BIN) reflecting the cellular interactome that can be analysed by the classical graph-theoretic algorithms. Therefore, in the second part of the review, we introduce the concept of BIN, i.e. a meta-network model of the complete molecular interactome developed by integrating various biological networks. The BIN meta-network model includes DNA-protein binding by the plethora of protein factors as well as chromatin interactions, therefore allowing connection of genomics with the downstream biomolecular processes present in a cell. As an illustration, we scrutinise the chromatin interactions mediated by the CTCF protein detected in a ChIA-PET experiment in the human lymphoblastoid cell line GM12878. In the corresponding BIN meta-network the DNA spatial proximity is represented as a graph model, combined with the Proteins-Interaction Network (PIN) of human proteome using the Gene Association Network (GAN). Furthermore, we enriched the BIN with the signalling and metabolic pathways and Gene Ontology (GO) terms to assert its functional context. Finally, we mapped the Single Nucleotide Polymorphisms (SNPs) from the GWAS studies and identified the chromatin mutational hot-spots associated with a significant enrichment of SNPs related to autoimmune diseases. Afterwards, we mapped Structural Variants (SVs) from healthy individuals of 1000 Genomes Project and identified an interesting example of the missing protein complex associated with protein Q6GYQ0 due to a deletion on chromosome 14. Such an analysis using the meta-network BIN model is therefore helpful in evaluating the influence of genetic variation on spatial organisation of the genome and its functional effect in a cell.

Low tolerance for transcriptional variation at cohesin genes is accompanied by functional links to disease-relevant pathways

BACKGROUND

The cohesin complex plays an essential role in genome organisation and cell division. A full complement of the cohesin complex and its regulators is important for normal development, since heterozygous mutations in genes encoding these components can be sufficient to produce a disease phenotype. The implication that genes encoding the cohesin subunits or cohesin regulators must be tightly controlled and resistant to variability in expression has not yet been formally tested.

METHODS

Here, we identify spatial-regulatory connections with potential to regulate expression of cohesin loci (Mitotic: SMC1A, SMC3, STAG1, STAG2, RAD21/RAD21-AS; Meiotic: SMC1B, STAG3, REC8, RAD21L1), cohesin-ring support genes (NIPBL, MAU2, WAPL, PDS5A, PDS5B) and CTCF, including linking their expression to that of other genes. We searched the genome-wide association studies (GWAS) catalogue for SNPs mapped or attributed to cohesin genes by GWAS (GWAS-attributed) and the GTEx catalogue for SNPs mapped to cohesin genes by cis-regulatory variants in one or more of 44 tissues across the human body (expression quantitative trail locus-attributed).

RESULTS

Connections that centre on the cohesin ring subunits provide evidence of coordinated regulation that has little tolerance for perturbation. We used the CoDeS3D SNP-gene attribution methodology to identify transcriptional changes across a set of genes coregulated with the cohesin loci that include biological pathways such as extracellular matrix production and proteasome-mediated protein degradation. Remarkably, many of the genes that are coregulated with cohesin loci are themselves intolerant to loss-of-function.

CONCLUSIONS

The results highlight the importance of robust regulation of cohesin genes and implicate novel pathways that may be important in the human cohesinopathy disorders.

Pathway analysis for genome-wide genetic variation data: Analytic principles, latest developments, and new opportunities

Pathway analysis, also known as gene-set enrichment analysis, is a multilocus analytic strategy that integrates a priori, biological knowledge into the statistical analysis of high-throughput genetics data. Originally developed for the studies of gene expression data, it has become a powerful analytic procedure for in-depth mining of genome-wide genetic variation data. Astonishing discoveries were made in the past years, uncovering genes and biological mechanisms underlying common and complex disorders. However, as massive amounts of diverse functional genomics data accrue, there is a pressing need for newer generations of pathway analysis methods that can utilize multiple layers of high-throughput genomics data. In this review, we provide an intellectual foundation of this powerful analytic strategy, as well as an update of the state-of-the-art in recent method developments. The goal of this review is threefold: (1) introduce the motivation and basic steps of pathway analysis for genome-wide genetic variation data; (2) review the merits and the shortcomings of classic and newly emerging integrative pathway analysis tools; and (3) discuss remaining challenges and future directions for further method developments.

RNA Biogenesis Instructs Functional Inter-Chromosomal Genome Architecture

Three-dimensional (3D) genome organization has emerged as an important layer of gene regulation in development and disease. The functional properties of chromatin folding within individual chromosomes (i.e., intra-chromosomal or in cis) have been studied extensively. On the other hand, interactions across different chromosomes (i.e., inter-chromosomal or in trans) have received less attention, being often regarded as background noise or technical artifacts. This viewpoint has been challenged by emerging evidence of functional relationships between specific trans chromatin interactions and epigenetic control, transcription, and splicing. Therefore, it is an intriguing possibility that the key processes involved in the biogenesis of RNAs may both shape and be in turn influenced by inter-chromosomal genome architecture. Here I present the rationale behind this hypothesis, and discuss a potential experimental framework aimed at its formal testing. I present a specific example in the cardiac myocyte, a well-studied post-mitotic cell whose development and response to stress are associated with marked rearrangements of chromatin topology both in cis and in trans. I argue that RNA polymerase II clusters (i.e., transcription factories) and foci of the cardiac-specific splicing regulator RBM20 (i.e., splicing factories) exemplify the existence of trans-interacting chromatin domains (TIDs) with important roles in cellular homeostasis. Overall, I propose that inter-molecular 3D proximity between co-regulated nucleic acids may be a pervasive functional mechanism in biology.

Identification of Reproducible BCL11A Alterations in Schizophrenia Through Individual-Level Prediction of Coexpression

Previous studies have provided evidence for an alteration of genetic coexpression in schizophrenia (SCZ). However, such analyses have thus far lacked biological specificity for individual genes, which may be critical for identifying illness-relevant effects. Therefore, we applied machine learning to identify gene-specific coexpression differences at the individual subject level and compared these between individuals with SCZ, bipolar disorder, major depressive disorder (MDD), autism spectrum disorder (ASD), and healthy controls. Utilizing transcriptome-wide gene expression data from 21 independent datasets, comprising a total of 9509 participants, we identified a reproducible decrease of BCL11A coexpression across 4 SCZ datasets that showed diagnostic specificity for SCZ when compared with ASD and MDD. We further demonstrate that individual-level coexpression differences can be combined in multivariate coexpression scores that show reproducible illness classification across independent datasets in SCZ and ASD. This study demonstrates that machine learning can capture gene-specific coexpression differences at the individual subject level for SCZ and identify novel biomarker candidates.

Cell type-specific genetic regulation of gene expression across human tissues

The Genotype-Tissue Expression (GTEx) project has identified expression and splicing quantitative trait loci in cis (QTLs) for the majority of genes across a wide range of human tissues. However, the functional characterization of these QTLs has been limited by the heterogeneous cellular composition of GTEx tissue samples. We mapped interactions between computational estimates of cell type abundance and genotype to identify cell type-interaction QTLs for seven cell types and show that cell type-interaction expression QTLs (eQTLs) provide finer resolution to tissue specificity than bulk tissue cis-eQTLs. Analyses of genetic associations with 87 complex traits show a contribution from cell type-interaction QTLs and enables the discovery of hundreds of previously unidentified colocalized loci that are masked in bulk tissue.

Genome-wide association analysis of pulse wave velocity traits provide new insights into the causal relationship between arterial stiffness and blood pressure

Background

The pathophysiology of arterial stiffness is not completely understood. Pulse wave velocity (PWV) is an established marker for arterial stiffness. We compare genetics of three PWV modes, namely carotid-femoral PWV (cfPWV), brachial-ankle (baPWV) and brachial-femoral (bfPWV), reflecting different vascular segments to analyse association with genetic variants, heritability and genetic correlation with other biological traits. Furthermore we searched for shared genetic architecture concerning PWV, blood pressure (BP) and coronary artery disease (CAD) and examined the causal relationship between PWV and BP.

Methods and results

We performed a genome-wide association study (GWAS) for cfPWV, baPWV and bfPWV in LIFE-Adult (N = 3,643–6,734). We analysed the overlap of detected genetic loci with those of BP and CAD and performed genetic correlation analyses. By bidirectional Mendelian Randomization, we assessed the causal relationships between PWV and BP. For cfPWV we identified a new locus with genome-wide significance near SLC4A7 on cytoband 3p24.1 (lead SNP rs939834: p = 2.05x10^-8). We replicated a known PWV locus on cytoband 14q32.2 near RP11-61O1.1 (lead SNPs: rs17773233, p = 1.38x10^-4; rs1381289, p = 1.91x10^-4) For baPWV we estimated a heritability of 28% and significant genetic correlation with hypertension (rg = 0.27, p = 6.65x10^-8). We showed a positive causal effect of systolic blood pressure on PWV modes (cfPWV: p = 1.51x10^-4; bfPWV: p = 1.45x10^-3; baPWV: p = 6.82x10^-15).

Conclusions

We identified a new locus for arterial stiffness and successfully replicated an earlier proposed locus. PWV shares common genetic architecture with BP and CAD. BP causally affects PWV. Larger studies are required to further unravel the genetic determinants and effects of PWV.

Investigation of prediction accuracy and the impact of sample size, ancestry, and tissue in transcriptome-wide association studies

In transcriptome-wide association studies (TWAS), gene expression values are predicted using genotype data and tested for association with a phenotype. The power of this approach to detect associations relies, at least in part, on the accuracy of the prediction. Here we compare the prediction accuracy of six different methods-LASSO, Ridge regression, Elastic net, Best Linear Unbiased Predictor, Bayesian Sparse Linear Mixed Model, and Random Forests-by performing cross-validation using data from the Geuvadis Project. We also examine prediction accuracy (a) at different sample sizes, (b) when ancestry of the prediction model training and testing populations is different, and (c) when the tissue used to train the model is different from the tissue to be predicted. We find that, for most genes, the expression cannot be accurately predicted, but in general sparse statistical models tend to outperform polygenic models at prediction. Average prediction accuracy is reduced when the model training set size is reduced or when predicting across ancestries and is marginally reduced when predicting across tissues. We conclude that using sparse statistical models and the development of large reference panels across multiple ethnicities and tissues will lead to better prediction of gene expression, and thus may improve TWAS power.

Genome-wide analysis of carotid plaque burden suggests a role of IL5 in men

BACKGROUND

Carotid artery plaque is an established marker of subclinical atherosclerosis with pronounced sex-dimorphism. Here, we aimed to identify genetic variants associated with carotid plaque burden (CPB) and to examine potential sex-specific genetic effects on plaque sizes.

METHODS AND RESULTS

We defined six operationalizations of CPB considering plaques in common carotid arteries, carotid bulb, and internal carotid arteries. We performed sex-specific genome-wide association analyses for all traits in the LIFE-Adult cohort (n = 727 men and n = 550 women) and tested significantly associated loci for sex-specific effects. In order to identify causal genes, we analyzed candidate gene expression data for correlation with CPB traits and corresponding sex-specific effects. Further, we tested if previously reported SNP associations with CAD and plaque prevalence are also associated with CBP. We found seven loci with suggestive significance for CPB (p<3.33x10-7), explaining together between 6 and 13% of the CPB variance. Sex-specific analysis showed a genome-wide significant hit for men at 5q31.1 (rs201629990, β = -0.401, p = 5.22x10-9), which was not associated in women (β = -0.127, p = 0.093) with a significant difference in effect size (p = 0.008). Analyses of gene expression data suggested IL5 as the most plausible candidate, as it reflected the same sex-specific association with CPBs (p = 0.037). Known plaque prevalence or CAD loci showed no enrichment in the association with CPB.

CONCLUSIONS

We showed that CPB is a complementary trait in analyzing genetics of subclinical atherosclerosis. We detected a novel locus for plaque size in men only suggesting a role of IL5. Several estrogen response elements in this locus point towards a functional explanation of the observed sex-specific effect.

Genomic resources for dissecting the role of non-protein coding variation in gene-environment interactions

The majority of single nucleotide variants (SNVs) identified in Genome Wide Association Studies (GWAS) fall within non-protein coding DNA and have the potential to alter gene expression. Non-protein coding DNA can control gene expression by acting as transcription factor (TF) binding sites or by regulating the organization of DNA into chromatin. SNVs in non-coding DNA sequences can disrupt TF binding and chromatin structure and this can result in pathology. Further, environmental health studies have shown that exposure to xenobiotics can disrupt the ability of TFs to regulate entire gene networks and result in pathology. However, there is a large amount of interindividual variability in exposure-linked health outcomes. One explanation for this heterogeneity is that genetic variation and exposure combine to disrupt gene regulation, and this eventually manifests in disease. Many resources exist that annotate common variants from GWAS and combine them with conservation, functional genomics, and TF binding data. These annotation tools provide clues regarding the biological implications of an SNV, as well as lead to the generation of hypotheses regarding potentially disrupted target genes, epigenetic markers, pathways, and cell types. Collectively this information can be used to predict how SNVs can alter an individual's response to exposure and disease risk. A basic understanding of the regulatory information contained within non-protein coding DNA is needed to predict the biological consequences of SNVs, and to determine how these SNVs impact exposure-related disease. We hope that this review will aid in the characterization of disease-associated genetic variation in the non-protein coding genome.

Tissue specific regulation of transcription in endometrium and association with disease

STUDY QUESTION

Are genetic effects on endometrial gene expression tissue specific and/or associated with reproductive traits and diseases?

SUMMARY ANSWER

Analyses of RNA-sequence data and individual genotype data from the endometrium identified novel and disease associated, genetic mechanisms regulating gene expression in the endometrium and showed evidence that these mechanisms are shared across biologically similar tissues.

WHAT IS KNOWN ALREADY

The endometrium is a complex tissue vital for female reproduction and is a hypothesized source of cells initiating endometriosis. Understanding genetic regulation specific to, and shared between, tissue types can aid the identification of genes involved in complex genetic diseases.

STUDY DESIGN, SIZE, DURATION

RNA-sequence and genotype data from 206 individuals was analysed and results were compared with large publicly available datasets.

PARTICIPANTS/MATERIALS, SETTING, METHODS

RNA-sequencing and genotype data from 206 endometrial samples was used to identify the influence of genetic variants on gene expression, via expression quantitative trait loci (eQTL) analysis and to compare these endometrial eQTLs with those in other tissues. To investigate the association between endometrial gene expression regulation and reproductive traits and diseases, we conducted a tissue enrichment analysis, transcriptome-wide association study (TWAS) and summary data-based Mendelian randomisation (SMR) analyses. Transcriptomic data was used to test differential gene expression between women with and without endometriosis.

MAIN RESULTS AND THE ROLE OF CHANCE

A tissue enrichment analysis with endometriosis genome-wide association study summary statistics showed that genes surrounding endometriosis risk loci were significantly enriched in reproductive tissues. A total of 444 sentinel cis-eQTLs (P < 2.57 × 10-9) and 30 trans-eQTLs (P < 4.65 × 10-13) were detected, including 327 novel cis-eQTLs in endometrium. A large proportion (85%) of endometrial eQTLs are present in other tissues. Genetic effects on endometrial gene expression were highly correlated with the genetic effects on reproductive (e.g. uterus, ovary) and digestive tissues (e.g. salivary gland, stomach), supporting a shared genetic regulation of gene expression in biologically similar tissues. The TWAS analysis indicated that gene expression at 39 loci is associated with endometriosis, including five known endometriosis risk loci. SMR analyses identified potential target genes pleiotropically or causally associated with reproductive traits and diseases including endometriosis. However, without taking account of genetic variants, a direct comparison between women with and without endometriosis showed no significant difference in endometrial gene expression.

LARGE SCALE DATA

The eQTL dataset generated in this study is available at http://reproductivegenomics.com.au/shiny/endo_eqtl_rna/. Additional datasets supporting the conclusions of this article are included within the article and the supplementary information files, or are available on reasonable request.

LIMITATIONS, REASONS FOR CAUTION

Data are derived from fresh tissue samples and expression levels are an average of expression from different cell types within the endometrium. Subtle cell-specifc expression changes may not be detected and differences in cell composition between samples and across the menstrual cycle will contribute to sample variability. Power to detect tissue specific eQTLs and differences between women with and without endometriosis was limited by the sample size in this study. The statistical approaches used in this study identify the likely gene targets for specific genetic risk factors, but not the functional mechanism by which changes in gene expression may influence disease risk.

WIDER IMPLICATIONS OF THE FINDINGS

Our results identify novel genetic variants that regulate gene expression in endometrium and the majority of these are shared across tissues. This allows analysis with large publicly available datasets to identify targets for female reproductive traits and diseases. Much larger studies will be required to identify genetic regulation of gene expression that will be specific to endometrium.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the National Health and Medical Research Council (NHMRC) under project grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321, GNT1083405 and GNT1107258. G.W.M is supported by a NHMRC Fellowship (GNT1078399). J.Y is supported by an ARC Fellowship (FT180100186). There are no competing interests.

Human brain arousal in the resting state: a genome-wide association study

Arousal affects cognition, emotion, and behavior and has been implicated in the etiology of psychiatric disorders. Although environmental conditions substantially contribute to the level of arousal, stable interindividual characteristics are well-established and a genetic basis has been suggested. Here we investigated the molecular genetics of brain arousal in the resting state by conducting a genome-wide association study (GWAS). We selected N = 1877 participants from the population-based LIFE-Adult cohort. Participants underwent a 20-min eyes-closed resting state EEG, which was analyzed using the computerized VIGALL 2.1 (Vigilance Algorithm Leipzig). At the SNP-level, GWAS analyses revealed no genome-wide significant locus (p < 5E-8), although seven loci were suggestive (p < 1E-6). The strongest hit was an expression quantitative trait locus (eQTL) of TMEM159 (lead-SNP: rs79472635, p = 5.49E-8). Importantly, at the gene-level, GWAS analyses revealed significant evidence for TMEM159 (p = 0.013, Bonferroni-corrected). By mapping our SNPs to the GWAS results from the Psychiatric Genomics Consortium, we found that all corresponding markers of TMEM159 showed nominally significant associations with Major Depressive Disorder (MDD; 0.006 ≤ p ≤ 0.011). More specifically, variants associated with high arousal levels have previously been linked to an increased risk for MDD. In line with this, the MetaXcan database suggests increased expression levels of TMEM159 in MDD, as well as Autism Spectrum Disorder, and Alzheimer's Disease. Furthermore, our pathway analyses provided evidence for a role of sodium/calcium exchangers in resting state arousal. In conclusion, the present GWAS identifies TMEM159 as a novel candidate gene which may modulate the risk for psychiatric disorders through arousal mechanisms. Our results also encourage the elaboration of the previously reported interrelations between ion-channel modulators, sleep-wake behavior, and psychiatric disorders.

Circulating Oxytocin Is Genetically Determined and Associated With Obesity and Impaired Glucose Tolerance

CONTEXT

Despite the emerging evidence on the role of oxytocin (OXT) in metabolic diseases, there is a lack of well-powered studies addressing the relationship of circulating OXT with obesity and diabetes.

OBJECTIVES AND DESIGN

Here, we measured OXT in a study cohort (n = 721; 396 women, 325 men; mean age ± SD, 47.7 ± 15.2 years) with subphenotypes related to obesity, including anthropometric traits such as body mass index [BMI (mean ± SD), 26.8 ± 4.6 kg/m2], waist-to-hip ratio (WHR; 0.88 ± 0.09), blood parameters (glucose, 5.32 ± 0.50 mmol/L; insulin, 5.3 ± 3.3 µU/mL), and oral glucose tolerance test to clarify the association with OXT. We also tested in a genome-wide association study (GWAS) whether the interindividual variation in OXT serum levels might be explained by genetic variation.

RESULTS

The OXT concentration was increased in subjects with elevated BMI and positively correlated with WHR, waist circumference, and triglyceride levels. The OXT concentration in subjects with BMI <25 kg/m2 was significantly lower (n = 256; 78.6 pg/mL) than in subjects with a BMI between 25 and 30 kg/m2 (n = 314; 98.5 pg/mL, P = 6 × 10-6) and with BMI >30 kg/m2 (n = 137; 106.4 pg/mL, P = 8 × 10-6). OXT levels were also positively correlated with plasma glucose and insulin and were elevated in subjects with impaired glucose tolerance (P = 4.6 × 10-3). Heritability of OXT was estimated at 12.8%. In a GWAS, two hits in linkage disequilibrium close (19 kb) to the OXT reached genome-wide significant association (top-hit rs12625893, P = 3.1 × 10-8, explained variance 3%).

CONCLUSIONS

Our data show that OXT is genetically affected by a variant near OXT and is associated with obesity and impaired glucose tolerance.

Genetic Association Study of Eight Steroid Hormones and Implications for Sexual Dimorphism of Coronary Artery Disease

CONTEXT

Steroid hormones are important regulators of physiological processes in humans and are under genetic control. A link to coronary artery disease (CAD) is supposed.

OBJECTIVE

Our main objective was to identify genetic loci influencing steroid hormone levels. As a secondary aim, we searched for causal effects of steroid hormones on CAD.

DESIGN

We conducted genome-wide meta-association studies for eight steroid hormones: cortisol, dehydroepiandrosterone sulfate (DHEAS), estradiol, and testosterone in two independent cohorts (LIFE-Adult, LIFE-Heart, maximum n = 7667), and progesterone, 17-hydroxyprogesterone, androstenedione, and aldosterone in LIFE-Heart only (maximum n = 2070). All genome-wide significant loci were tested for sex interactions. Furthermore, we tested whether previously reported CAD single-nucleotide polymorphisms were associated with our steroid hormone panel and investigated causal links between hormone levels and CAD status using Mendelian randomization (MR) approaches.

RESULTS

We discovered 15 novel associated loci for 17-hydroxyprogesterone, progesterone, DHEAS, cortisol, androstenedione, and estradiol. Five of these loci relate to genes directly involved in steroid metabolism, that is, CYP21A1, CYP11B1, CYP17A1, STS, and HSD17B12, almost completing the set of steroidogenic enzymes with genetic associations. Sexual dimorphisms were found for seven of the novel loci. Other loci correspond, for example, to the WNT4/β-catenin pathway. MR revealed that cortisol, androstenedione, 17-hydroxyprogesterone, and DHEA-S had causal effects on CAD. We also observed enrichment of cortisol and testosterone associations among known CAD hits.

CONCLUSION

Our study greatly improves insight into genetic regulation of steroid hormones and their dependency on sex. These results could serve as a basis for analyzing sexual dimorphism in other complex diseases.

Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria.

Imputed gene associations identify replicable trans-acting genes enriched in transcription pathways and complex traits

Regulation of gene expression is an important mechanism through which genetic variation can affect complex traits. A substantial portion of gene expression variation can be explained by both local (cis) and distal (trans) genetic variation. Much progress has been made in uncovering cis-acting expression quantitative trait loci (cis-eQTL), but trans-eQTL have been more difficult to identify and replicate. Here we take advantage of our ability to predict the cis component of gene expression coupled with gene mapping methods such as PrediXcan to identify high confidence candidate trans-acting genes and their targets. That is, we correlate the cis component of gene expression with observed expression of genes in different chromosomes. Leveraging the shared cis-acting regulation across tissues, we combine the evidence of association across all available Genotype-Tissue Expression Project tissues and find 2,356 trans-acting/target gene pairs with high mappability scores. Reassuringly, trans-acting genes are enriched in transcription and nucleic acid binding pathways and target genes are enriched in known transcription factor binding sites. Interestingly, trans-acting genes are more significantly associated with selected complex traits and diseases than target or background genes, consistent with percolating trans effects. Our scripts and summary statistics are publicly available for future studies of trans-acting gene regulation.

Power, false discovery rate and Winner's Curse in eQTL studies

Investigation of the genetic architecture of gene expression traits has aided interpretation of disease and trait-associated genetic variants; however, key aspects of expression quantitative trait loci (eQTL) study design and analysis remain understudied. We used extensive, empirically driven simulations to explore eQTL study design and the performance of various analysis strategies. Across multiple testing correction methods, false discoveries of genes with eQTLs (eGenes) were substantially inflated when false discovery rate (FDR) control was applied to all tests and only appropriately controlled using hierarchical procedures. All multiple testing correction procedures had low power and inflated FDR for eGenes whose causal SNPs had small allele frequencies using small sample sizes (e.g. frequency <10% in 100 samples), indicating that even moderately low frequency eQTL SNPs (eSNPs) in these studies are enriched for false discoveries. In scenarios with ≥80% power, the top eSNP was the true simulated eSNP 90% of the time, but substantially less frequently for very common eSNPs (minor allele frequencies >25%). Overestimation of eQTL effect sizes, so-called ‘Winner’s Curse’, was common in low and moderate power settings. To address this, we developed a bootstrap method (BootstrapQTL) that led to more accurate effect size estimation. These insights provide a foundation for future eQTL studies, especially those with sampling constraints and subtly different conditions.

A catalog of genetic loci associated with kidney function from analyses of a million individuals

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Sodium valproate rescues expression of TRANK1 in iPSC-derived neural cells that carry a genetic variant associated with serious mental illness

Biological characterization of genetic variants identified in genome-wide association studies (GWAS) remains a substantial challenge. Here we used human-induced pluripotent stem cells (iPSC) and their neural derivatives to characterize common variants on chromosome 3p22 that have been associated by GWAS with major mental illnesses. IPSC-derived neural progenitor cells carrying the risk allele of the single nucleotide polymorphism (SNP), rs9834970, displayed lower baseline TRANK1 expression that was rescued by chronic treatment with therapeutic dosages of valproic acid (VPA). VPA had the greatest effects on TRANK1 expression in iPSC, NPC, and astrocytes. Although rs9834970 has no known function, we demonstrated that a nearby SNP, rs906482, strongly affects binding by the transcription factor, CTCF, and that the high-affinity allele usually occurs on haplotypes carrying the rs9834970 risk allele. Decreased expression of TRANK1 perturbed expression of many genes involved in neural development and differentiation. These findings have important implications for the pathophysiology of major mental illnesses and the development of novel therapeutics.

Characterizing human genomic coevolution in locus-gene regulatory interactions

Background

Coevolution has been used to identify and predict interactions and functional relationships between proteins of many different organisms including humans. Current efforts in annotating the human genome increasingly show that non-coding DNA sequence has important functional and regulatory interactions. Furthermore, regulatory elements do not necessarily reside in close proximity of the coding region for their target genes.

Results

We characterize coevolution as it appears in locus-gene interactions in the human genome, focusing on expression Quantitative Trait - Locus (eQTL) interactions. Our results show that in these interactions the conservation status of the loci is predictive of the conservation status of their target genes. Furthermore, comparing the phylogenetic histories of intra-chromosomal pairs of loci and transcription start sites, we find that pairs that appear coevolved are enriched for cis-eQTL interactions. Exploring this property we found that coevolution might be useful in prioritizing association tests in cis-eQTL detection.

Conclusions

The relationship between the conservation status of pairs of loci and protein coding transcription start sites reveal correlations with regulatory interactions. Pairs that appear coevolved are enriched for intra-chromosomal regulatory interactions, thus our results suggest that measures of coevolution can be useful for prediction and detection of new interactions. Measures of coevolution are genome-wide and could potentially be used to prioritize the detection of distant or inter-chromosomal interactions such as trans-eQTL interactions in the human genome.

Analysis of the genomic architecture of a complex trait locus in hypertensive rat models links Tmem63c to kidney damage

Unraveling the genetic susceptibility of complex diseases such as chronic kidney disease remains challenging. Here, we used inbred rat models of kidney damage associated with elevated blood pressure for the comprehensive analysis of a major albuminuria susceptibility locus detected in these models. We characterized its genomic architecture by congenic substitution mapping, targeted next-generation sequencing, and compartment-specific RNA sequencing analysis in isolated glomeruli. This led to prioritization of transmembrane protein Tmem63c as a novel potential target. Tmem63c is differentially expressed in glomeruli of allele-specific rat models during onset of albuminuria. Patients with focal segmental glomerulosclerosis exhibited specific TMEM63C loss in podocytes. Functional analysis in zebrafish revealed a role for tmem63c in mediating the glomerular filtration barrier function. Our data demonstrate that integrative analysis of the genomic architecture of a complex trait locus is a powerful tool for identification of new targets such as Tmem63c for further translational investigation.

The human kidneys filter the entire volume of the blood about 300 times each day. This ability depends on specialized cells, known as podocytes, which wrap around some of the blood vessels in the kidney. These cells control which molecules leave the blood based on their size. Normally large molecules like proteins are blocked, while smaller molecules including waste products, toxins, excess water and salts pass through into the urine.

If this filtration system is damaged, by high blood pressure, for example, it can lead to chronic kidney disease. A hallmark of this disease, often called CKD for short, is high levels of the protein albumin in the urine. Previous studies involving rats with high blood pressure have found several regions of the genome that contribute to high levels of albumin in the urine, including one on chromosome 6. However, this region contains several genes and it was unclear which genes affected the condition.

Schulz et al. set out to narrow down the list and find specific genes that might contribute to elevated albumin in the urine of rats with high blood pressure. This search identified the gene for a protein called TMEM63c as a likely candidate. This protein spans the outer membrane of podocyte cells. Analysis of kidney biopsies showed that patients with chronic kidney disease also had low levels of this protein in their podocytes. Further experiments, this time in zebrafish, showed that reducing the activity of the gene for tmem63c led to damaged podocytes and a leakier filter in the kidneys.

The results suggest that this gene plays an important role in the integrity of the kidneys filtration barrier. It is possible that faulty versions of this gene are behind some cases of chronic kidney disease. If this proves to be the case, a better understanding of the role of this gene may lead to new treatments for the condition.

Shaping the nebulous enhancer in the era of high-throughput assays and genome editing

Since the 1st discovery of transcriptional enhancers in 1981, their textbook definition has remained largely unchanged in the past 37 years. With the emergence of high-throughput assays and genome editing, which are switching the paradigm from bottom-up discovery and testing of individual enhancers to top-down profiling of enhancer activities genome-wide, it has become increasingly evidenced that this classical definition has left substantial gray areas in different aspects. Here we survey a representative set of recent research articles and report the definitions of enhancers they have adopted. The results reveal that a wide spectrum of definitions is used usually without the definition stated explicitly, which could lead to difficulties in data interpretation and downstream analyses. Based on these findings, we discuss the practical implications and suggestions for future studies.

Nuclear genetic regulation of the human mitochondrial transcriptome

Mitochondria play important roles in cellular processes and disease, yet little is known about how the transcriptional regime of the mitochondrial genome varies across individuals and tissues. By analyzing >11,000 RNA-sequencing libraries across 36 tissue/cell types, we find considerable variation in mitochondrial-encoded gene expression along the mitochondrial transcriptome, across tissues and between individuals, highlighting the importance of cell-type specific and post-transcriptional processes in shaping mitochondrial-encoded RNA levels. Using whole-genome genetic data we identify 64 nuclear loci associated with expression levels of 14 genes encoded in the mitochondrial genome, including missense variants within genes involved in mitochondrial function (TBRG4, MTPAP and LONP1), implicating genetic mechanisms that act in trans across the two genomes. We replicate ~21% of associations with independent tissue-matched datasets and find genetic variants linked to these nuclear loci that are associated with cardio-metabolic phenotypes and Vitiligo, supporting a potential role for variable mitochondrial-encoded gene expression in complex disease.

Novel hypotheses emerging from GWAS in migraine?

Recent technical advances in genetics made large-scale genome-wide association studies (GWAS) in migraine feasible and have identified over 40 common DNA sequence variants that affect risk for migraine types. Most of the variants, which are all single nucleotide polymorphisms (SNPs), show robust association with migraine as evidenced by the fact that the vast majority replicate in subsequent independent studies. However, despite thorough bioinformatic efforts aimed at linking the migraine risk SNPs with genes and their molecular pathways, there remains quite some discussion as to how successful this endeavour has been, and their current practical use for the diagnosis and treatment of migraine patients. Although existing genetic information seems to favour involvement of vascular mechanisms, but also neuronal and other mechanisms such as metal ion homeostasis and neuronal migration, the complexity of the underlying genetic pathophysiology presents challenges to advancing genetic knowledge to clinical use. A major issue is to what extent one can rely on bioinformatics to pinpoint the actual disease genes, and from this the linked pathways. In this Commentary, we will provide an overview of findings from GWAS in migraine, current hypotheses of the disease pathways that emerged from these findings, and some of the major drawbacks of the approaches used to identify the genes and pathways. We argue that more functional research is urgently needed to turn the hypotheses that emerge from GWAS in migraine to clinically useful information.

Genome-wide meta-analysis identifies novel determinants of circulating serum progranulin

Progranulin is a secreted protein with important functions in processes including immune and inflammatory response, metabolism and embryonic development. The present study aimed at identification of genetic factors determining progranulin concentrations. We conducted a genome-wide association meta-analysis for serum progranulin in three independent cohorts from Europe: Sorbs (N = 848) and KORA (N = 1628) from Germany and PPP-Botnia (N = 335) from Finland (total N = 2811). Single nucleotide polymorphisms (SNPs) associated with progranulin levels were replicated in two additional German cohorts: LIFE-Heart Study (Leipzig; N = 967) and Metabolic Syndrome Berlin Potsdam (Berlin cohort; N = 833). We measured mRNA expression of genes in peripheral blood mononuclear cells (PBMC) by micro-arrays and performed mRNA expression quantitative trait and expression-progranulin association studies to functionally substantiate identified loci. Finally, we conducted siRNA silencing experiments in vitro to validate potential candidate genes within the associated loci. Heritability of circulating progranulin levels was estimated at 31.8% and 26.1% in the Sorbs and LIFE-Heart cohort, respectively. SNPs at three loci reached study-wide significance (rs660240 in CELSR2-PSRC1-MYBPHL-SORT1, rs4747197 in CDH23-PSAP and rs5848 in GRN) explaining 19.4%/15.0% of the variance and 61%/57% of total heritability in the Sorbs/LIFE-Heart Study. The strongest evidence for association was at rs660240 (P = 5.75 × 10-50), which was also associated with mRNA expression of PSRC1 in PBMC (P = 1.51 × 10-21). Psrc1 knockdown in murine preadipocytes led to a consecutive 30% reduction in progranulin secretion. In conclusion, the present meta-GWAS combined with mRNA expression identified three loci associated with progranulin and supports the role of PSRC1 in the regulation of progranulin secretion.

Most chromatin interactions are not in linkage disequilibrium

Chromatin interactions and linkage disequilibrium (LD) are both pairwise measurements between genomic loci that show block patterns along mammalian chromosomes. Their values are generally high for sites that are nearby in the linear genome but abruptly drop across block boundaries. One function of chromatin boundaries is to insulate regulatory domains from one another. Since recombination is depressed within genes and between distal regulatory elements and their promoters, we hypothesized that LD and chromatin contact frequency might be correlated genome-wide with the boundaries of LD blocks and chromatin domains frequently coinciding. To comprehensively address this question, we compared chromatin contacts in 22 cell types to LD across billions of pairs of loci in the human genome. These computationally intensive analyses revealed that there is no concordance between LD and chromatin interactions, even at genomic distances below 25 kilobases (kb) where both tend to be high. At genomic distances where LD is approximately zero, chromatin interactions are frequent. While LD is somewhat elevated between distal regulatory elements and their promoters, LD block boundaries are depleted—not enriched—at chromatin boundaries. Finally, gene expression and ontology data suggest that chromatin contacts identify regulatory variants more reliably than do LD and genomic proximity. We conclude that the genomic architectures of genetic and physical interactions are independent, with important implications for gene regulatory evolution, interpretation of genetic association studies, and precision medicine.