QTLbase: an integrative resource for quantitative trait loci across multiple human molecular phenotypes

Potential association between PSCA rs2976395 functional variant and pancreatic cancer risk

Correlated regions of systemic interindividual variation (CoRSIV) represent a small proportion of the human genome showing DNA methylation patterns that are the same in all human tissues, are different among individuals, and are partially regulated by genetic variants in cis. In this study we aimed at investigating single-nucleotide polymorphisms (SNPs) within CoRSIVs and their involvement with pancreatic ductal adenocarcinoma (PDAC) risk. We analyzed 29,099 CoRSIV-SNPs and 133,615 CoRSIV-mQTLs in 14,394 cases and 247,022 controls of European and Asian descent. We observed that the A allele of the rs2976395 SNP was associated with increased PDAC risk in Europeans (p = 2.81 × 10-5). This SNP lies in the prostate stem cell antigen gene and is in perfect linkage disequilibrium with a variant (rs2294008) that has been reported to be associated with risk of many other cancer types. The A allele is associated with the DNA methylation level of the gene according to the PanCan-meQTL database and with overexpression according to QTLbase. The expression of the gene has been observed to be deregulated in many tumors of the gastrointestinal tract including pancreatic cancer; however, functional studies are needed to elucidate the function relevance of the association.

A multivariate to multivariate approach for voxel-wise genome-wide association analysis

The joint analysis of imaging-genetics data facilitates the systematic investigation of genetic effects on brain structures and functions with spatial specificity. We focus on voxel-wise genome-wide association analysis, which may involve trillions of single nucleotide polymorphism (SNP)-voxel pairs. We attempt to identify underlying organized association patterns of SNP-voxel pairs and understand the polygenic and pleiotropic networks on brain imaging traits. We propose a bi-clique graph structure (ie, a set of SNPs highly correlated with a cluster of voxels) for the systematic association pattern. Next, we develop computational strategies to detect latent SNP-voxel bi-cliques and an inference model for statistical testing. We further provide theoretical results to guarantee the accuracy of our computational algorithms and statistical inference. We validate our method by extensive simulation studies, and then apply it to the whole genome genetic and voxel-level white matter integrity data collected from 1052 participants of the human connectome project. The results demonstrate multiple genetic loci influencing white matter integrity measures on splenium and genu of the corpus callosum.

Crosstalk between epitranscriptomic and epigenomic modifications and its implication in human diseases

Crosstalk between N6-methyladenosine (m6A) and epigenomes is crucial for gene regulation, but its regulatory directionality and disease significance remain unclear. Here, we utilize quantitative trait loci (QTLs) as genetic instruments to delineate directional maps of crosstalk between m6A and two epigenomic traits, DNA methylation (DNAme) and H3K27ac. We identify 47 m6A-to-H3K27ac and 4,733 m6A-to-DNAme and, in the reverse direction, 106 H3K27ac-to-m6A and 61,775 DNAme-to-m6A regulatory loci, with differential genomic location preference observed for different regulatory directions. Integrating these maps with complex diseases, we prioritize 20 genome-wide association study (GWAS) loci for neuroticism, depression, and narcolepsy in brain; 1,767 variants for asthma and expiratory flow traits in lung; and 249 for coronary artery disease, blood pressure, and pulse rate in muscle. This study establishes disease regulatory paths, such as rs3768410-DNAme-m6A-asthma and rs56104944-m6A-DNAme-hypertension, uncovering locus-specific crosstalk between m6A and epigenomic layers and offering insights into regulatory circuits underlying human diseases.

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

BACKGROUND

Genome-wide association studies (GWAS) have revealed a multitude of candidate genetic variants affecting the risk of developing complex traits and diseases. However, the highlighted regions are typically in the non-coding genome, and uncovering the functional causative single nucleotide variants (SNVs) is challenging. Prioritization of variants is commonly based on genomic annotation with markers of active regulatory elements, but current approaches still poorly predict functional variants. To address this, we systematically analyze six markers of active regulatory elements for their ability to identify functional variants.

RESULTS

We benchmark against molecular quantitative trait loci (molQTL) from assays of regulatory element activity that identify allelic effects on DNA-binding factor occupancy, reporter assay expression, and chromatin accessibility. We identify the combination of DNase footprints and divergent enhancer RNA (eRNA) as markers for functional variants. This signature provides high precision, but with a trade-off of low recall, thus substantially reducing candidate variant sets to prioritize variants for functional validation. We present this as a framework called FINDER-Functional SNV IdeNtification using DNase footprints and eRNA.

CONCLUSIONS

We demonstrate the utility to prioritize variants using leukocyte count trait and analyze variants in linkage disequilibrium with a lead variant to predict a functional variant in asthma. Our findings have implications for prioritizing variants from GWAS, in development of predictive scoring algorithms, and for functionally informed fine mapping approaches.

Genetic regulation of m6A RNA methylation and its contribution in human complex diseases

N6-methyladenosine (m6A) has been established as the most prevalent chemical modification in message RNA (mRNA), playing an essential role in determining the fate of RNA molecules. Dysregulation of m6A has been revealed to lead to abnormal physiological conditions and cause various types of human diseases. Recent studies have delineated the genetic regulatory maps for m6A methylation by mapping the quantitative trait loci of m6A (m6A-QTLs), thereby building up the regulatory circuits linking genetic variants, m6A, and human complex traits. Here, we review the recent discoveries concerning the genetic regulatory maps of m6A, describing the methodological and technical details of m6A-QTL identification, and introducing the key findings of the cis- and trans-acting drivers of m6A. We further delve into the tissue- and ethnicity-specificity of m6A-QTL, the association with other molecular phenotypes in light of genetic regulation, the regulators underlying m6A genetics, and importantly, the functional roles of m6A in mediating human complex diseases. Lastly, we discuss potential research avenues that can accelerate the translation of m6A genetics studies toward the development of therapies for human genetic diseases.

Admixture mapping of cognitive function in diverse Hispanic and Latino adults: Results from the Hispanic Community Health Study/Study of Latinos

INTRODUCTION

We conducted admixture mapping and fine-mapping analyses to identify ancestry-of-origin loci influencing cognitive abilities.

METHODS

We estimated the association of local ancestry intervals across the genome with five neurocognitive measures in 7140 diverse Hispanic and Latino adults (mean age 55 years). We prioritized genetic variants in associated loci and tested them for replication in four independent cohorts.

RESULTS

We identified nine local ancestry-associated regions for the five neurocognitive measures. There was strong biological support for the observed associations to cognitive function at all loci and there was statistical evidence of independent replication at 4q12, 9p22.1, and 13q12.13.

DISCUSSION

Our study identified multiple novel loci harboring genes implicated in cognitive functioning and dementia, and uncovered ancestry-relevant genetic variants. It adds to our understanding of the genetic architecture of cognitive function in Hispanic and Latino adults and demonstrates the power of admixture mapping to discover unique haplotypes influencing cognitive function, complementing genome-wide association studies.

HIGHLIGHTS

We identified nine ancestry-of-origin chromosomal regions associated with five neurocognitive traits. In each associated region, we identified single nucleotide polymorphisms (SNPs) that explained, at least in part, the admixture signal and were tested for replication in independent samples of Black, non-Hispanic White, and Hispanic/Latino adults with the same or similar neurocognitive tests. Statistical evidence of independent replication of the prioritized SNPs was observed for three of the nine associations, at chr4q12, chr9p22.1, and chr13q12.13. At all loci, there was strong biological support for the observed associations to cognitive function and dementia, prioritizing genes such as KIT, implicated in autophagic clearance of neurotoxic proteins and on mast cell and microglial-mediated inflammation; SLC24A2, implicated in synaptic plasticity associated with learning and memory; and MTMR6, implicated in phosphoinositide lipids metabolism.

Genetically proxied HTRA1 protease activity and circulating levels independently predict risk of ischemic stroke and coronary artery disease

Genetic variants in HTRA1 are associated with stroke risk. However, the mechanisms mediating this remain largely unknown, as does the full spectrum of phenotypes associated with genetic variation in HTRA1. Here we show that rare HTRA1 variants are linked to ischemic stroke in the UK Biobank and BioBank Japan. Integrating data from biochemical experiments, we next show that variants causing loss of protease function associated with ischemic stroke, coronary artery disease and skeletal traits in the UK Biobank and MyCode cohorts. Moreover, a common variant modulating circulating HTRA1 mRNA and protein levels enhances the risk of ischemic stroke and coronary artery disease while lowering the risk of migraine and macular dystrophy in genome-wide association study, UK Biobank, MyCode and BioBank Japan data. We found no interaction between proxied HTRA1 activity and levels. Our findings demonstrate the role of HTRA1 for cardiovascular diseases and identify two mechanisms as potential targets for therapeutic interventions.

Combined Genetic Association and Differed Expression Analysis of UBE2L3 Uncovers a Genetic Regulatory Role of (Immuno)proteasome in IgA Nephropathy

Introduction

IgA nephropathy (IgAN) is a leading cause of end-stage renal disease. The exact pathogenesis of IgAN is not well defined, but some genetic studies have led to a novel discovery that the (immuno)proteasome probably plays an important role in IgAN.

Methods

We firstly analyzed the association of variants in the UBE2L3 region with susceptibility to IgAN in 3,495 patients and 9,101 controls, and then analyzed the association between lead variant and clinical phenotypes in 1,803 patients with regular follow-up data. The blood mRNA levels of members of the ubiquitin-proteasome system including UBE2L3 were analyzed in peripheral blood mononuclear cells from 53 patients and 28 healthy controls. The associations between UBE2L3 and the expression levels of genes involved in Gd-IgA1 production were also explored.

Results

The rs131654 showed the most significant association signal in UBE2L3 region (OR: 1.10, 95% CI: 1.04-1.16, p = 2.29 × 10-3), whose genotypes were also associated with the levels of Gd-IgA1 (p = 0.04). The rs131654 was observed to exert cis-eQTL effects on UBE2L3 in various tissues and cell types, particularly in immune cell types in multiple databases. The UBE2L3, LUBAC, and proteasome subunits were highly expressed in patients compared with healthy controls. High expression levels of UBE2L3 were not only associated with higher proteinuria (r = 0.34, p = 0.01) and lower eGFR (r = -0.28, p = 0.04), but also positively correlated with the gene expression of LUBAC and other proteasome subunits. Additionally, mRNA expression levels of UBE2L3 were also positively correlated with IL-6 and RELA, but negatively correlated with the expression levels of the key enzyme in the process of glycosylation including C1GALT1 and C1GALT1C1.

Conclusion

In conclusion, by combined genetic association and differed expression analysis of UBE2L3, our data support a role of genetically conferred dysregulation of the (immuno)proteasome in regulating galactose-deficient IgA1 in the development of IgAN.

From GWASs toward Mechanistic Understanding with Case Studies in Dermatogenetics

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

Obesity-Dependent Association of the rs10454142 PPP1R21 with Breast Cancer

The purpose of this work was to find a link between the breast cancer (BC)-risk effects of sex hormone-binding globulin (SHBG)-associated polymorphisms and obesity. The study was conducted on a sample of 1498 women (358 BC; 1140 controls) who, depending on the presence/absence of obesity, were divided into two groups: obese (119 BC; 253 controls) and non-obese (239 BC; 887 controls). Genotyping of nine SHBG-associated single nucleotide polymorphisms (SNP)-rs17496332 PRMT6, rs780093 GCKR, rs10454142 PPP1R21, rs3779195 BAIAP2L1, rs440837 ZBTB10, rs7910927 JMJD1C, rs4149056 SLCO1B1, rs8023580 NR2F2, and rs12150660 SHBG-was executed, and the BC-risk impact of these loci was analyzed by logistic regression separately in each group of obese/non-obese women. We found that the BC-risk effect correlated by GWAS with the SHBG-level polymorphism rs10454142 PPP1R21 depends on the presence/absence of obesity. The SHBG-lowering allele C rs10454142 PPP1R21 has a risk value for BC in obese women (allelic model: CvsT, OR = 1.52, 95%CI = 1.10-2.11, and pperm = 0.013; additive model: CCvsTCvsTT, OR = 1.71, 95%CI = 1.15-2.62, and pperm = 0.011; dominant model: CC + TCvsTT, OR = 1.95, 95%CI = 1.13-3.37, and pperm = 0.017) and is not associated with the disease in women without obesity. SNP rs10454142 PPP1R21 and 10 proxy SNPs have adipose-specific regulatory effects (epigenetic modifications of promoters/enhancers, DNA interaction with 51 transcription factors, eQTL/sQTL effects on five genes (PPP1R21, RP11-460M2.1, GTF2A1L, STON1-GTF2A1L, and STON1), etc.), can be "likely cancer driver" SNPs, and are involved in cancer-significant pathways. In conclusion, our study detected an obesity-dependent association of the rs10454142 PPP1R21 with BC in women.

The Impact of Inherited Genetic Variation on DNA Methylation in Prostate Cancer and Benign Tissues of African American and European American Men

BACKGROUND

American men of African ancestry (AA) have higher prostate cancer incidence and mortality rates compared with American men of European ancestry (EA). Differences in genetic susceptibility mechanisms may contribute to this disparity.

METHODS

To gain insights into the regulatory mechanisms of prostate cancer susceptibility variants, we tested the association between SNPs and DNA methylation (DNAm) at nearby CpG sites across the genome in benign and cancer prostate tissue from 74 AA and 74 EA men. Genome-wide SNP data (from benign tissue) and DNAm were generated using Illumina arrays.

RESULTS

Among AA men, we identified 6,298 and 2,641 cis-methylation QTLs (meQTL; FDR of 0.05) in benign and tumor tissue, respectively, with 6,960 and 1,700 detected in EA men. We leveraged genome-wide association study (GWAS) summary statistics to identify previously reported prostate cancer GWAS signals likely to share a common causal variant with a detected meQTL. We identified nine GWAS-meQTL pairs with strong evidence of colocalization (four in EA benign, three in EA tumor, two in AA benign, and three in AA tumor). Among these colocalized GWAS-meQTL pairs, we identified colocalizing expression quantitative trait loci (eQTL) impacting four eGenes with known roles in tumorigenesis.

CONCLUSIONS

These findings highlight epigenetic regulatory mechanisms by which prostate cancer-risk SNPs can modify local DNAm and/or gene expression in prostate tissue.

IMPACT

Overall, our findings showed general consistency in the meQTL landscape of AA and EA men, but meQTLs often differ by tissue type (normal vs. cancer). Ancestry-based linkage disequilibrium differences and lack of AA representation in GWAS decrease statistical power to detect colocalization for some regions.

Novel genetic association of the Furin gene polymorphism rs1981458 with COVID-19 severity among Indian populations

SARS CoV-2, the causative agent for the ongoing COVID-19 pandemic, it enters the host cell by activating the ACE2 receptor with the help of two proteasesi.e., Furin and TMPRSS2. Therefore, variations in these genes may account for differential susceptibility and severity between populations. Previous studies have shown that the role of ACE2 and TMPRSS2 gene variants in understanding COVID-19 susceptibility among Indian populations. Nevertheless, a knowledge gap exists concerning the COVID-19 susceptibility of Furin gene variants among diverse South Asian ethnic groups. Investigating the role of Furin gene variants and their global phylogeographic structure is essential to comprehensively understanding COVID-19 susceptibility in these populations. We have used 450 samples from diverse Indian states and performed linear regression to analyse the Furin gene variant's with COVID-19 Case Fatality Rate (CFR) that could be epidemiologically associated with disease severity outcomes. Associated genetic variants were further evaluated for their expression and regulatory potential through various Insilco analyses. Additionally, we examined the Furin gene using next-generation sequencing (NGS) data from 393 diverse global samples, with a particular emphasis on South Asia, to investigate its Phylogeographic structure among diverse world populations. We found a significant positive association for the SNP rs1981458 with COVID-19 CFR (p < 0.05) among diverse Indian populations at different timelines of the first and second waves. Further, QTL and other regulatory analyses showed various significant associations for positive regulatory roles of rs1981458 and Furin gene, mainly in Immune cells and virus infection process, highlighting their role in host immunity and viral assembly and processing. The Furin protein-protein interaction suggested that COVID-19 may contribute to Pulmonary arterial hypertension via a typical inflammation mechanism. The phylogeographic architecture of the Furin gene demonstrated a closer genetic affinity of South Asia with West Eurasian populations. Therefore, it is worth proposing that for the Furin gene, the COVID-19 susceptibility of South Asians will be more similar to the West Eurasian population. Our previous studies on the ACE2 and TMPRSS2 genes showed genetic affinity of South Asian with East Eurasians and West Eurasians, respectively. Therefore, with the collective information from these three important genes (ACE2, TMPRSS2 and Furin) we modelled COVID-19 susceptibilityof South Asia in between these two major ancestries with an inclination towards West Eurasia. In conclusion, this study, for the first time, concluded the role of rs1981458 in COVID-19 severity among the Indian population and outlined its regulatory potential.This study also highlights that the genetic structure for COVID-19 susceptibilityof South Asia is distinct, however, inclined to the West Eurasian population. We believe this insight may be utilised as a genetic biomarker to identify vulnerable populations, which might be directly relevant for developing policies and allocating resources more effectively during an epidemic.

Genetic Variations in Elements of the Oxytocinergic Pathway are Associated with Attention/Hyperactivity Problems and Anxiety Problems in Childhood

Genetic alterations related to oxytocin system seem to influence the neurobiology of attention-deficit hyperactivity disorder and anxiety problems leading to greater functional, social and emotional impairment. Here, we analyzed the association of OXTR rs2254298 and CD38 rs6449182 variants with attention/hyperactivity problems and anxiety problems in children. The study enrolled 292 children and adjusted regression model revealed OXTR rs2254298 AA genotype as a risk factor for attention deficit/hyperactivity problems (PR: 2.37; PadjFDR = 0.006), attention problems (PR: 2.71; PadjFDR = 0.003) and anxiety problems (PR: 1.92; PadjFDR = 0.018). CD38 rs6449182 G allele showed as a risk factor for attention deficit/hyperactivity problems (PR: 1.56; PadjFDR = 0.028). Moreover, in silico approach for regulatory roles found markers that influence chromatin accessibility and transcription capacity. Together, these data provide genetic information of oxytocin in developmental and behavioral disorders opening a range of opportunities for future studies that clarify their neurobiology in childhood.

Sex-Hormone-Binding Globulin Gene Polymorphisms and Breast Cancer Risk in Caucasian Women of Russia

In our work, the associations of GWAS (genome-wide associative studies) impact for sex-hormone-binding globulin (SHBG)-level SNPs with the risk of breast cancer (BC) in the cohort of Caucasian women of Russia were assessed. The work was performed on a sample of 1498 women (358 BC patients and 1140 control (non BC) subjects). SHBG correlated in previously GWAS nine polymorphisms such as rs780093 GCKR, rs17496332 PRMT6, rs3779195 BAIAP2L1, rs10454142 PPP1R21, rs7910927 JMJD1C, rs4149056 SLCO1B1, rs440837 ZBTB10, rs12150660 SHBG, and rs8023580 NR2F2 have been genotyped. BC risk effects of allelic and non-allelic SHBG-linked gene SNPs interactions were detected by regression analysis. The risk genetic factor for BC developing is an SHBG-lowering allele variant C rs10454142 PPP1R21 ([additive genetic model] OR = 1.31; 95%CI = 1.08-1.65; pperm = 0.024; power = 85.26%), which determines 0.32% of the cancer variance. Eight of the nine studied SHBG-related SNPs have been involved in cancer susceptibility as part of nine different non-allelic gene interaction models, the greatest contribution to which is made by rs10454142 PPP1R21 (included in all nine models, 100%) and four more SNPs-rs7910927 JMJD1C (five models, 55.56%), rs17496332 PRMT6 (four models, 44.44%), rs780093 GCKR (four models, 44.44%), and rs440837 ZBTB10 (four models, 44.44%). For SHBG-related loci, pronounced functionality in the organism (including breast, liver, fibroblasts, etc.) was predicted in silico, having a direct relationship through many pathways with cancer pathophysiology. In conclusion, our results demonstrated the involvement of SHBG-correlated genes polymorphisms in BC risk in Caucasian women in Russia.

Genetic predisposition of BMP7 polymorphisms to lumbar disk herniation in the Chinese Han population

Bone morphogenetic protein 7 (BMP7) can induce skeletal formation, promote the differentiation of chondrocytes and osteoblasts, and ameliorate intervertebral disc degeneration. The study was designed to evaluate the relationship of BMP7 variants to LDH risk in the Chinese Han population. BMP7 variants were genotyped with the Agena MassARRAY system among 690 LDH patients and 690 healthy controls. The odds ratio (OR) and 95% confidence interval (CI) were calculated by logistic regression. Multi-factor dimension reduction (MDR) (version 3.0.2) software was used to evaluate the effect of BMP7 variant-variant interaction on the susceptibility to LDH. Here, the risk-reducing association between rs230189 and LDH occurrence was found (p = 0.005, OR = 0.79). Specially, rs230189 was associated with decreased LDH risk in females (p = 0.001, OR = 0.60), elder group (p = 0.025, OR = 0.76), subjects with BMI < 24 kg/m2 (p = 0.027, OR = 0.48), nonsmokers (p = 0.001, OR = 0.66), and nondrinkers (p = 0.011, OR = 0.72). Moreover, rs1321862 might be the risk factor for LDH susceptibility among the participants with BMI < 24 kg/m2 (p = 0.024, OR = 1.84). MDR results displayed that rs230189 was the greatest attribution factor on LDH risk in the single-locus model, with an information gain of 0.44%. The present study demonstrated that BMP7 rs230189 g.55771443A>C may play a protective role against LDH risk. Our findings may help to understand the potential mechanism of BMP7 in LDH susceptibility.

COLOCdb: a comprehensive resource for multi-model colocalization of complex traits

Large-scale genome-wide association studies (GWAS) have provided profound insights into complex traits and diseases. Yet, deciphering the fine-scale molecular mechanisms of how genetic variants manifest to cause the phenotypes remains a daunting task. Here, we present COLOCdb (https://ngdc.cncb.ac.cn/colocdb), a comprehensive genetic colocalization database by integrating more than 3000 GWAS summary statistics and 13 types of xQTL to date. By employing two representative approaches for the colocalization analysis, COLOCdb deposits results from three key components: (i) GWAS-xQTL, pair-wise colocalization between GWAS loci and different types of xQTL, (ii) GWAS-GWAS, pair-wise colocalization between the trait-associated genetic loci from GWASs and (iii) xQTL-xQTL, pair-wise colocalization between the genetic loci associated with molecular phenotypes in xQTLs. These results together represent the most comprehensive colocalization analysis, which also greatly expands the list of shared variants with genetic pleiotropy. We expect that COLOCdb can serve as a unique and useful resource in advancing the discovery of new biological mechanisms and benefit future functional studies.

Disentangling Accelerated Cognitive Decline from the Normal Aging Process and Unraveling Its Genetic Components: A Neuroimaging-Based Deep Learning Approach

Background

The progressive cognitive decline, an integral component of Alzheimer's disease (AD), unfolds in tandem with the natural aging process. Neuroimaging features have demonstrated the capacity to distinguish cognitive decline changes stemming from typical brain aging and AD between different chronological points.

Objective

To disentangle the normal aging effect from the AD-related accelerated cognitive decline and unravel its genetic components using a neuroimaging-based deep learning approach.

Methods

We developed a deep-learning framework based on a dual-loss Siamese ResNet network to extract fine-grained information from the longitudinal structural magnetic resonance imaging (MRI) data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. We then conducted genome-wide association studies (GWAS) and post-GWAS analyses to reveal the genetic basis of AD-related accelerated cognitive decline.

Results

We used our model to process data from 1,313 individuals, training it on 414 cognitively normal people and predicting cognitive assessment for all participants. In our analysis of accelerated cognitive decline GWAS, we identified two genome-wide significant loci: APOE locus (chromosome 19 p13.32) and rs144614292 (chromosome 11 p15.1). Variant rs144614292 (G > T) has not been reported in previous AD GWA studies. It is within the intronic region of NELL1, which is expressed in neurons and plays a role in controlling cell growth and differentiation. The cell-type-specific enrichment analysis and functional enrichment of GWAS signals highlighted the microglia and immune-response pathways.

Conclusions

Our deep learning model effectively extracted relevant neuroimaging features and predicted individual cognitive decline. We reported a novel variant (rs144614292) within the NELL1 gene.

Fast multiple-trait genome-wide association analysis for correlated longitudinal measurements

Large-scale longitudinal biobank data can be leveraged to identify genetic variation contributing to human diseases progression and traits trajectories. While methods for genome-wide association studies (GWAS) of multiple correlated traits have been proposed, an efficient multiple-trait approach to model longitudinal phenotypes is not currently available. We developed GAMUT, a genome-wide association approach for multiple longitudinal traits. GAMUT employs a mixed-effects model to fit longitudinal outcomes where a fast algorithm for inversion by recursive partitioning of the random effects submatrix is introduced. To evaluate performance of the algorithms introduced and assess their statistical power and type I error, stochastic simulation was conducted. Consistent with our expectation, power was greater for cross-sectional (CS) than longitudinal (LT) effects, particularly with a diminishing LT/CS ratio. With a minimum minor allele count of 3 within genotype by time categories, observed type I error was roughly equal to theoretical genome-wide significance. Additionally, 28 blood-based biomarkers measured at 2 time points on participants of the UK Biobank were used to compare GAMUT against single-trait standard and longitudinal GWAS (including rate of change). Across all biomarkers, we observed 539 (CS) and 248 (LT) significant independent variants for the GAMUT method, and 513 (CS) and 30 (LT) for single-trait longitudinal GWAS, respectively. Only 37 variants were identified by modeling rates of change using standard GWAS.

xQTLbiolinks: a comprehensive and scalable tool for integrative analysis of molecular QTLs

Genome-wide association studies (GWAS) have identified thousands of disease-associated non-coding variants, posing urgent needs for functional interpretation. Molecular Quantitative Trait Loci (xQTLs) such as eQTLs serve as an essential intermediate link between these non-coding variants and disease phenotypes and have been widely used to discover disease-risk genes from many population-scale studies. However, mining and analyzing the xQTLs data presents several significant bioinformatics challenges, particularly when it comes to integration with GWAS data. Here, we developed xQTLbiolinks as the first comprehensive and scalable tool for bulk and single-cell xQTLs data retrieval, quality control and pre-processing from public repositories and our integrated resource. In addition, xQTLbiolinks provided a robust colocalization module through integration with GWAS summary statistics. The result generated by xQTLbiolinks can be flexibly visualized or stored in standard R objects that can easily be integrated with other R packages and custom pipelines. We applied xQTLbiolinks to cancer GWAS summary statistics as case studies and demonstrated its robust utility and reproducibility. xQTLbiolinks will profoundly accelerate the interpretation of disease-associated variants, thus promoting a better understanding of disease etiologies. xQTLbiolinks is available at https://github.com/lilab-bioinfo/xQTLbiolinks.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Meta-analysis of age-related cognitive decline reveals a novel locus for the attention domain and implicates a COVID-19-related gene for global cognitive function

INTRODUCTION

Cognitive abilities have substantial heritability throughout life, as shown by twin- and population-based studies. However, there is limited understanding of the genetic factors related to cognitive decline in aging across neurocognitive domains.

METHODS

We conducted a meta-analysis on 3045 individuals aged ≥65, derived from three population-based cohorts, to identify genetic variants associated with the decline of five neurocognitive domains (attention, memory, executive function, language, visuospatial function) and global cognitive decline. We also conducted gene-based and functional bioinformatics analyses.

RESULTS

Apolipoprotein E (APOE)4 was significantly associated with decline of memory (p = 5.58E-09) and global cognitive function (p = 1.84E-08). We identified a novel association with attention decline on chromosome 9, rs6559700 (p = 2.69E-08), near RASEF. Gene-based analysis also identified a novel gene, TMPRSS11D, involved in the activation of SARS-CoV-2, to be associated with the decline in global cognitive function (p = 4.28E-07).

DISCUSSION

Domain-specific genetic studies can aid in the identification of novel genes and pathways associated with decline across neurocognitive domains.

HIGHLIGHTS

rs6559700 was associated with decline of attention. APOE4 was associated with decline of memory and global cognitive decline. TMPRSS11D, a gene involved in the activation of SARS-CoV-2, was implicated in global cognitive decline. Cognitive domain abilities had both unique and shared molecular pathways across the domains.

Genome-wide analyses identify NEAT1 as genetic modifier of age at onset of amyotrophic lateral sclerosis

BACKGROUND

Patients with amyotrophic lateral sclerosis (ALS) demonstrate great heterogeneity in the age at onset (AAO), which is closely related to the course of disease. However, most genetic studies focused on the risk of ALS, while the genetic background underlying AAO of ALS is still unknown.

METHODS

To identify genetic determinants influencing AAO of ALS, we performed genome-wide association analysis using a Cox proportional hazards model in 2,841 patients with ALS (Ndiscovery = 2,272, Nreplication = 569) in the Chinese population. We further conducted colocalization analysis using public cis-eQTL dataset, and Mendelian randomization analysis to identify risk factors for AAO of ALS. Finally, functional experiments including dual-luciferase reporter assay and RT-qPCR were performed to explore the regulatory effect of the target variant.

RESULTS

The total heritability of AAO of ALS was ~ 0.24. One novel locus rs10128627 (FRMD8) was significantly associated with earlier AAO by ~ 3.15 years (P = 1.54E-08, beta = 0.31, SE = 0.05). This locus was cis-eQTL of NEAT1 in multiple brain tissues and blood. Colocalization analysis detected association signals at this locus between AAO of ALS and expression of NEAT1. Furthermore, functional exploration supported the variant rs10128627 was associated with upregulated expression of NEAT1 in cell models and patients with ALS. Causal inference suggested higher total cholesterol, low-density lipoprotein, and eosinophil were nominally associated with earlier AAO of ALS, while monocyte might delay the AAO.

CONCLUSIONS

Collective evidence from genetic, bioinformatic, and functional results suggested NEAT1 as a key player in the disease progression of ALS. These findings improve the current understanding of the genetic role in AAO of ALS, and provide a novel target for further research on the pathogenesis and therapeutic options to delay the disease onset.

Genome-wide association study of hippocampal blood-oxygen-level-dependent-cerebral blood flow correlation in Chinese Han population

Correlation between blood-oxygen-level-dependent (BOLD) and cerebral blood flow (CBF) has been used as an index of neurovascular coupling. Hippocampal BOLD-CBF correlation is associated with neurocognition, and the reduced correlation is associated with neuropsychiatric disorders. We conducted the first genome-wide association study of the hippocampal BOLD-CBF correlation in 4,832 Chinese Han subjects. The hippocampal BOLD-CBF correlation had an estimated heritability of 16.2-23.9% and showed reliable genome-wide significant association with a locus at 3q28, in which many variants have been linked to neuroimaging and cerebrospinal fluid markers of Alzheimer's disease. Gene-based association analyses showed four significant genes (GMNC, CRTC2, DENND4B, and GATAD2B) and revealed enrichment for mast cell calcium mobilization, microglial cell proliferation, and ubiquitin-related proteolysis pathways that regulate different cellular components of the neurovascular unit. This is the first unbiased identification of the association of hippocampal BOLD-CBF correlation, providing fresh insights into the genetic architecture of hippocampal neurovascular coupling.

Genome-wide analysis identifies novel loci influencing plasma apolipoprotein E concentration and Alzheimer's disease risk

The APOE 2/3/4 polymorphism is the greatest genetic risk factor for Alzheimer's disease (AD). This polymorphism is also associated with variation in plasma ApoE level; while APOE*4 lowers, APOE*2 increases ApoE level. Lower plasma ApoE level has also been suggested to be a risk factor for incident dementia. To our knowledge, no large genome-wide association study (GWAS) has been reported on plasma ApoE level. This study aimed to identify new genetic variants affecting plasma ApoE level as well as to test if baseline ApoE level is associated with cognitive function and incident dementia in a longitudinally followed cohort of the Ginkgo Evaluation of Memory (GEM) study. Baseline plasma ApoE concentration was measured in 3031 participants (95.4% European Americans (EAs)). GWAS analysis was performed on 2580 self-identified EAs where both genotype and plasma ApoE data were available. Lower ApoE concentration was associated with worse cognitive function, but not with incident dementia. As expected, the risk for AD increased from E2/2 through to E4/4 genotypes (P for trend = 4.8E-75). In addition to confirming the expected and opposite associations of APOE*2 (P = 4.73E-79) and APOE*4 (P = 8.73E-12) with ApoE level, GWAS analysis revealed nine additional independent signals in the APOE region, and together they explained about 22% of the variance in plasma ApoE level. We also identified seven new loci on chromosomes 1, 4, 5, 7, 11, 12 and 20 (P range = 5.49E-08 to 5.36E-10) that explained about 9% of the variance in ApoE level. Plasma ApoE level-associated independent variants, especially in the APOE region, were also associated with AD risk and amyloid deposition in the brain, indicating that genetically determined ApoE level variation may be a risk factor for developing AD. These results improve our understanding of the genetic determinants of plasma ApoE level and their potential value in affecting AD risk.

The effects of unsaturated fatty acids on psoriasis: A two-sample Mendelian randomization study

Unsaturated fatty acids have been reported to be associated with the risk of psoriasis. However, the causal relationship between them remains unclear This study aimed to explore the causal relationship between unsaturated FAs and psoriasis. Firstly, we obtained genome-wide association study (GWAS) data for psoriasis from the FINNGEN database (number of cases = 4510, number of controls = 212,242) and different FA levels (number of samples = 114,999) from the IEU OpenGWAS Project. Secondly, the genetic correlation coefficient was calculated using linkage disequilibrium fractional regression. Thirdly, a two-sample Mendelian randomization (MR) analysis was performed using independent instrumental variables (p < 5 × 10-8) to determine the direction of randomization. Finally, expression quantitative trait loci (eQTL)-related analyses of common single nucleotide polymorphisms (SNPs) were carried out to explore the potential molecular mechanisms of unsaturated FAs affecting psoriasis. We found that an increase in the ratio of monounsaturated fatty acids (MUFAs) to total fatty acids could increase the risk of psoriasis (inverse-variance weighted [IVW], adjusted odds ratio [OR] = 1.175; adjusted 95% confidence interval [CI] = 1.045-1.321; adjusted p = .007). However, an increase in the ratio of polyunsaturated fatty acids (PUFAa) to total fatty acids could decrease the risk of psoriasis (IVW, adjusted OR = 0.754; adjusted 95% CI = 0.631-0.901; adjusted p = .002). Moreover, an increase in the ratio of PUFAs to MUFAs could decrease the risk of psoriasis (IVW, adjusted OR = 0.823; adjusted 95% CI = 0.715-0.948; adjusted p = .007). The heterogeneity of data was eliminated, and pleiotropy was not detected. There was no statistical difference in the MR analysis of other fatty acids indices with psoriasis. Further, no statistically significant evidence was found to verify a causal relationship between psoriasis and fatty acid levels in reverse MR. Functional enrichment analysis showed that these eQTL related to common SNPs were mainly involved in organic ion transport, choline metabolism, and the expression of key metabolic factors mediated by PKA, ChREBP, and PP2A. Our study indicated that the ratio of MUFAs to total fatty acids had a positive causal effect on psoriasis, while the ratio of PUFAs to total fatty acids and the ratio of PUFAs to MUFAs had a negative causal effect on psoriasis. Moreover, PKA-, PP2A-, and ChREBP-mediated activation of metabolic factors may play an important role in this process.

Genomics of Treatable Traits in Asthma

The astounding number of genetic variants revealed in the 15 years of genome-wide association studies of asthma has not kept pace with the goals of translational genomics. Moving asthma diagnosis from a nonspecific umbrella term to specific phenotypes/endotypes and related traits may provide insights into features that may be prevented or alleviated by therapeutical intervention. This review provides an overview of the different asthma endotypes and phenotypes and the genomic findings from asthma studies using patient stratification strategies and asthma-related traits. Asthma genomic research for treatable traits has uncovered novel and previously reported asthma loci, primarily through studies in Europeans. Novel genomic findings for asthma phenotypes and related traits may arise from multi-trait and specific phenotyping strategies in diverse populations.

Investigating Neuron Degeneration in Huntington's Disease Using RNA-Seq Based Transcriptome Study

Huntington's disease (HD) is a progressive neurodegenerative disorder caused due to a CAG repeat expansion in the huntingtin (HTT) gene. The primary symptoms of HD include motor dysfunction such as chorea, dystonia, and involuntary movements. The primary motor cortex (BA4) is the key brain region responsible for executing motor/movement activities. Investigating patient and control samples from the BA4 region will provide a deeper understanding of the genes responsible for neuron degeneration and help to identify potential markers. Previous studies have focused on overall differential gene expression and associated biological functions. In this study, we illustrate the relationship between variants and differentially expressed genes/transcripts. We identified variants and their associated genes along with the quantification of genes and transcripts. We also predicted the effect of variants on various regulatory activities and found that many variants are regulating gene expression. Variants affecting miRNA and its targets are also highlighted in our study. Co-expression network studies revealed the role of novel genes. Function interaction network analysis unveiled the importance of genes involved in vesicle-mediated transport. From this unified approach, we propose that genes expressed in immune cells are crucial for reducing neuron death in HD.

Disentangling accelerated cognitive decline from the normal aging process and unraveling its genetic components: A neuroimaging-based deep learning approach

Background

The progressive cognitive decline that is an integral component of AD unfolds in tandem with the natural aging process. Neuroimaging features have demonstrated the capacity to distinguish cognitive decline changes stemming from typical brain aging and Alzheimer's disease between different chronological points.

Methods

We developed a deep-learning framework based on dual-loss Siamese ResNet network to extract fine-grained information from the longitudinal structural magnetic resonance imaging (MRI) data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. We then conducted genome-wide association studies (GWAS) and post-GWAS analyses to reveal the genetic basis of AD-related accelerated cognitive decline.

Results

We used our model to process data from 1,313 individuals, training it on 414 cognitively normal people and predicting cognitive assessment for all participants. In our analysis of accelerated cognitive decline GWAS, we identified two genome-wide significant loci: APOE locus (chromosome 19 p13.32) and rs144614292 (chromosome 11 p15.1). Variant rs144614292 (G>T) has not been reported in previous AD GWA studies. It is within the intronic region of NELL1, which is expressed in neuron and plays a role in controlling cell growth and differentiation. In addition, MUC7 and PROL1/OPRPNon chromosome 4 were significant at the gene level. The cell-type-specific enrichment analysis and functional enrichment of GWAS signals highlighted the microglia and immune-response pathways. Furthermore, we found that the cognitive decline slope GWAS was positively correlated with previous AD GWAS.

Conclusion

Our deep learning model was demonstrated effective on extracting relevant neuroimaging features and predicting individual cognitive decline. We reported a novel variant (rs144614292) within the NELL1 gene. Our approach has the potential to disentangle accelerated cognitive decline from the normal aging process and to determine its related genetic factors, leveraging opportunities for early intervention.

MTMR3 risk alleles enhance Toll Like Receptor 9-induced IgA immunity in IgA nephropathy

Multiple genome-wide association studies (GWASs) have reproducibly identified the MTMR3/HORMAD2/LIF/OSM locus to be associated with IgA nephropathy (IgAN). However, the causal variant(s), implicated gene(s), and altered mechanisms remain poorly understood. Here, we performed fine-mapping analyses based on GWAS datasets encompassing 2762 IgAN cases and 5803 control individuals, and identified rs4823074 as the candidate causal variant that intersects the MTMR3 promoter in B-lymphoblastoid cells. Mendelian randomization studies suggested the risk allele may modulate disease susceptibility by affecting serum IgA levels through increased MTMR3 expression. Consistently, elevated MTMR3 expression in peripheral blood mononuclear cells was observed in patients with IgAN. Further mechanistic studies in vitro demonstrated that MTMR3 increased IgA production dependent upon its phosphatidylinositol 3-phosphate binding domain. Moreover, our study provided the in vivo functional evidence that Mtmr3-/- mice exhibited defective Toll Like Receptor 9-induced IgA production, glomerular IgA deposition, as well as mesangial cell proliferation. RNA-seq and pathway analyses showed that MTMR3 deficiency resulted in an impaired intestinal immune network for IgA production. Thus, our results support the role of MTMR3 in IgAN pathogenesis by enhancing Toll Like Receptor 9-induced IgA immunity.

The Impact of CC16 on Pulmonary Epithelial-Driven Host Responses during Mycoplasma pneumoniae Infection in Mouse Tracheal Epithelial Cells

Club Cell Secretory Protein (CC16) plays many protective roles within the lung; however, the complete biological functions, especially regarding the pulmonary epithelium during infection, remain undefined. We have previously shown that CC16-deficient (CC16-/-) mouse tracheal epithelial cells (MTECs) have enhanced Mp burden compared to CC16-sufficient (WT) MTECs; therefore, in this study, we wanted to further define how the pulmonary epithelium responds to infection in the context of CC16 deficiency. Using mass spectrometry and quantitative proteomics to analyze proteins secreted apically from MTECs grown at an air-liquid interface, we investigated the protective effects that CC16 elicits within the pulmonary epithelium during Mycoplasma pneumoniae (Mp) infection. When challenged with Mp, WT MTECs have an overall reduction in apical protein secretion, whereas CC16-/- MTECs have increased apical protein secretion compared to their unchallenged controls. Following Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) assessment, many of the proteins upregulated from CC16-/- MTECS (unchallenged and during Mp infection) were related to airway remodeling, which were not observed by WT MTECs. These findings suggest that CC16 may be important in providing protection within the pulmonary epithelium during respiratory infection with Mp, which is the major causative agent of community-acquired pneumoniae.

mGWAS-Explorer 2.0: Causal Analysis and Interpretation of Metabolite-Phenotype Associations

Metabolomics-based genome-wide association studies (mGWAS) are key to understanding the genetic regulations of metabolites in complex phenotypes. We previously developed mGWAS-Explorer 1.0 to link single-nucleotide polymorphisms (SNPs), metabolites, genes and phenotypes for hypothesis generation. It has become clear that identifying potential causal relationships between metabolites and phenotypes, as well as providing deep functional insights, are crucial for further downstream applications. Here, we introduce mGWAS-Explorer 2.0 to support the causal analysis between >4000 metabolites and various phenotypes. The results can be interpreted within the context of semantic triples and molecular quantitative trait loci (QTL) data. The underlying R package is released for reproducible analysis. Using two case studies, we demonstrate that mGWAS-Explorer 2.0 is able to detect potential causal relationships between arachidonic acid and Crohn's disease, as well as between glycine and coronary heart disease.

Identification and validation of six acute myocardial infarction-associated variants, including a novel prognostic marker for cardiac mortality

Background

Acute myocardial infarction (AMI) is one of the leading causes of death worldwide, and approximately half of AMI-related deaths occur before the affected individual reaches the hospital. The present study aimed to identify and validate genetic variants associated with AMI and their role as prognostic markers.

Materials and methods

We conducted a replication study of 29 previously identified novel loci containing 85 genetic variants associated with early-onset AMI using a new independent set of 2,920 Koreans [88 patients with early- and 1,085 patients with late-onset AMI, who underwent percutaneous coronary intervention (PCI), and 1,747 healthy controls].

Results

Of the 85 previously reported early-onset variants, six were confirmed in our genome-wide association study with a false discovery rate of less than 0.05. Notably, rs12639023, a cis-eQTL located in the intergenic region between LINC02005 and CNTN3, significantly increased longitudinal cardiac mortality and recurrent AMI. CNTN3 is known to play a role in altering vascular permeability. Another variant, rs78631167, located upstream of PLAUR and known to function in fibrinolysis, was moderately replicated in this study. By surveying the nearby genomic region around rs78631167, we identified a significant novel locus (rs8109584) located 13 bp downstream of rs78631167. The present study showed that six of the early-onset variants of AMI are applicable to both early- and late-onset cases.

Conclusion

Our results confirm markers that can potentially be utilized to predict, screen, prevent, and treat candidate patients with AMI and highlight the potential of rs12639023 as a prognostic marker for cardiac mortality in AMI.

PopTradeOff: A database for exploring population-specificity of adaptive evolution, disease susceptibility, and drug responsiveness

The influence of adaptive evolution on disease susceptibility has drawn attention; however, the extent of the influence, whether favored mutations also influence drug responses, and whether the associations between the three are population-specific remain unknown. Using a reported deep learning network to integrate seven statistical tests for detecting selection signals, we predicted favored mutations in the genomes of 17 human populations and integrated these favored mutations with reported GWAS sites and drug response-related variants into the database PopTradeOff (http://www.gaemons.net/PopFMIntro). The database also contains genome annotation information on the SNP, sequence, gene, and pathway levels. The preliminary data analyses suggest that substantial associations exist between adaptive evolution, disease susceptibility, and drug responses and that the associations are highly population-specific. The database may be valuable for disease studies, drug development, and personalized medicine.

Association between HSPA8 Gene Variants and Ischemic Stroke: A Pilot Study Providing Additional Evidence for the Role of Heat Shock Proteins in Disease Pathogenesis

HSPA8 is involved in many stroke-associated cellular processes, playing a pivotal role in the protein quality control system. Here we report the results of the pilot study aimed at determining whether HSPA8 SNPs are linked to the risk of ischemic stroke (IS). DNA samples from 2139 Russians (888 IS patients and 1251 healthy controls) were genotyped for tagSNPs (rs1461496, rs10892958, and rs1136141) in the HSPA8 gene using probe-based PCR. SNP rs10892958 of HSPA8 was associated with an increased risk (risk allele G) of IS in smokers (OR = 1.37; 95% CI = 1.07-1.77; p = 0.01) and patients with low fruit and vegetable consumption (OR = 1.36; 95% CI = 1.14-1.63; p = 0.002). SNP rs1136141 of HSPA8 was also associated with an increased risk of IS (risk allele A) exclusively in smokers (OR = 1.68; 95% CI = 1.23-2.28; p = 0.0007) and in patients with a low fruit and vegetable intake (OR = 1.29; 95% CI = 1.05-1.60; p = 0.04). Sex-stratified analysis revealed an association of rs10892958 HSPA8 with an increased risk of IS in males (risk allele G; OR = 1.30; 95% CI = 1.05-1.61; p = 0.01). Thus, SNPs rs10892958 and rs1136141 in the HSPA8 gene represent novel genetic markers of IS.

SERPINE1 mRNA Binding Protein 1 Is Associated with Ischemic Stroke Risk: A Comprehensive Molecular-Genetic and Bioinformatics Analysis of SERBP1 SNPs

The SERBP1 gene is a well-known regulator of SERPINE1 mRNA stability and progesterone signaling. However, the chaperone-like properties of SERBP1 have recently been discovered. The present pilot study investigated whether SERBP1 SNPs are associated with the risk and clinical manifestations of ischemic stroke (IS). DNA samples from 2060 unrelated Russian subjects (869 IS patients and 1191 healthy controls) were genotyped for 5 common SNPs-rs4655707, rs1058074, rs12561767, rs12566098, and rs6702742 SERBP1-using probe-based PCR. The association of SNP rs12566098 with an increased risk of IS (risk allele C; p = 0.001) was observed regardless of gender or physical activity level and was modified by smoking, fruit and vegetable intake, and body mass index. SNP rs1058074 (risk allele C) was associated with an increased risk of IS exclusively in women (p = 0.02), non-smokers (p = 0.003), patients with low physical activity (p = 0.04), patients with low fruit and vegetable consumption (p = 0.04), and BMI ≥25 (p = 0.007). SNPs rs1058074 (p = 0.04), rs12561767 (p = 0.01), rs12566098 (p = 0.02), rs6702742 (p = 0.036), and rs4655707 (p = 0.04) were associated with shortening of activated partial thromboplastin time. Thus, SERBP1 SNPs represent novel genetic markers of IS. Further studies are required to confirm the relationship between SERBP1 polymorphism and IS risk.

Genome-wide gene by environment study of time spent in daylight and chronotype identifies emerging genetic architecture underlying light sensitivity

STUDY OBJECTIVES

Light is the primary stimulus for synchronizing the circadian clock in humans. There are very large interindividual differences in the sensitivity of the circadian clock to light. Little is currently known about the genetic basis for these interindividual differences.

METHODS

We performed a genome-wide gene-by-environment interaction study (GWIS) in 280 897 individuals from the UK Biobank cohort to identify genetic variants that moderate the effect of daytime light exposure on chronotype (individual time of day preference), acting as "light sensitivity" variants for the impact of daylight on the circadian system.

RESULTS

We identified a genome-wide significant SNP mapped to the ARL14EP gene (rs3847634; p < 5 × 10-8), where additional minor alleles were found to enhance the morningness effect of daytime light exposure (βGxE = -.03, SE = 0.005) and were associated with increased gene ARL14EP expression in brain and retinal tissues. Gene-property analysis showed light sensitivity loci were enriched for genes in the G protein-coupled glutamate receptor signaling pathway and genes expressed in Per2+ hypothalamic neurons. Linkage disequilibrium score regression identified Bonferroni significant genetic correlations of greater light sensitivity GWIS with later chronotype and shorter sleep duration. Greater light sensitivity was nominally genetically correlated with insomnia symptoms and risk for post-traumatic stress disorder (PTSD).

CONCLUSIONS

This study is the first to assess light as an important exposure in the genomics of chronotype and is a critical first step in uncovering the genetic architecture of human circadian light sensitivity and its links to sleep and mental health.

An autoimmune pleiotropic SNP modulates IRF5 alternative promoter usage through ZBTB3-mediated chromatin looping

Genetic sharing is extensively observed for autoimmune diseases, but the causal variants and their underlying molecular mechanisms remain largely unknown. Through systematic investigation of autoimmune disease pleiotropic loci, we found most of these shared genetic effects are transmitted from regulatory code. We used an evidence-based strategy to functionally prioritize causal pleiotropic variants and identify their target genes. A top-ranked pleiotropic variant, rs4728142, yielded many lines of evidence as being causal. Mechanistically, the rs4728142-containing region interacts with the IRF5 alternative promoter in an allele-specific manner and orchestrates its upstream enhancer to regulate IRF5 alternative promoter usage through chromatin looping. A putative structural regulator, ZBTB3, mediates the allele-specific loop to promote IRF5-short transcript expression at the rs4728142 risk allele, resulting in IRF5 overactivation and M1 macrophage polarization. Together, our findings establish a causal mechanism between the regulatory variant and fine-scale molecular phenotype underlying the dysfunction of pleiotropic genes in human autoimmunity.

GPGPS: a robust prognostic gene pair signature of glioma ensembling IDH mutation and 1p/19q co-deletion

MOTIVATION

Many studies have shown that IDH mutation and 1p/19q co-deletion can serve as prognostic signatures of glioma. Although these genetic variations affect the expression of one or more genes, the prognostic value of gene expression related to IDH and 1p/19q status is still unclear.

RESULTS

We constructed an ensemble gene pair signature for the risk evaluation and survival prediction of glioma based on the prior knowledge of the IDH and 1p/19q status. First, we separately built two gene pair signatures IDH-GPS and 1p/19q-GPS and elucidated that they were useful transcriptome markers projecting from corresponding genome variations. Then, the gene pairs in these two models were assembled to develop an integrated model named Glioma Prognostic Gene Pair Signature (GPGPS), which demonstrated high area under the curves (AUCs) to predict 1-, 3- and 5-year overall survival (0.92, 0.88 and 0.80) of glioma. GPGPS was superior to the single GPSs and other existing prognostic signatures (avg AUC = 0.70, concordance index = 0.74). In conclusion, the ensemble prognostic signature with 10 gene pairs could serve as an independent predictor for risk stratification and survival prediction in glioma. This study shed light on transferring knowledge from genetic alterations to expression changes to facilitate prognostic studies.

AVAILABILITY AND IMPLEMENTATION

Codes are available at https://github.com/Kimxbzheng/GPGPS.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DNA methylation QTL mapping across diverse human tissues provides molecular links between genetic variation and complex traits

Studies of DNA methylation (DNAm) in solid human tissues are relatively scarce; tissue-specific characterization of DNAm is needed to understand its role in gene regulation and its relevance to complex traits. We generated array-based DNAm profiles for 987 human samples from the Genotype-Tissue Expression (GTEx) project, representing 9 tissue types and 424 subjects. We characterized methylome and transcriptome correlations (eQTMs), genetic regulation in cis (mQTLs and eQTLs) across tissues and e/mQTLs links to complex traits. We identified mQTLs for 286,152 CpG sites, many of which (>5%) show tissue specificity, and mQTL colocalizations with 2,254 distinct GWAS hits across 83 traits. For 91% of these loci, a candidate gene link was identified by integration of functional maps, including eQTMs, and/or eQTL colocalization, but only 33% of loci involved an eQTL and mQTL present in the same tissue type. With this DNAm-focused integrative analysis, we contribute to the understanding of molecular regulatory mechanisms in human tissues and their impact on complex traits.

Integrative Meta-Analysis of Huntington's Disease Transcriptome Landscape

Huntington's disease (HD) is a neurodegenerative disorder with autosomal dominant inheritance caused by glutamine expansion in the Huntingtin gene (HTT). Striatal projection neurons (SPNs) in HD are more vulnerable to cell death. The executive striatal population is directly connected with the Brodmann Area (BA9), which is mainly involved in motor functions. Analyzing the disease samples from BA9 from the SRA database provides insights related to neuron degeneration, which helps to identify a promising therapeutic strategy. Most gene expression studies examine the changes in expression and associated biological functions. In this study, we elucidate the relationship between variants and their effect on gene/downstream transcript expression. We computed gene and transcript abundance and identified variants from RNA-seq data using various pipelines. We predicted the effect of genome-wide association studies (GWAS)/novel variants on regulatory functions. We found that many variants affect the histone acetylation pattern in HD, thereby perturbing the transcription factor networks. Interestingly, some variants affect miRNA binding as well as their downstream gene expression. Tissue-specific network analysis showed that mitochondrial, neuroinflammation, vasculature, and angiogenesis-related genes are disrupted in HD. From this integrative omics analysis, we propose that abnormal neuroinflammation acts as a two-edged sword that indirectly affects the vasculature and associated energy metabolism. Rehabilitation of blood-brain barrier functionality and energy metabolism may secure the neuron from cell death.

Systems genomics in age-related macular degeneration

Genomic studies in age-related macular degeneration (AMD) have identified genetic variants that account for the majority of AMD risk. An important next step is to understand the functional consequences and downstream effects of the identified AMD-associated genetic variants. Instrumental for this next step are 'omics' technologies, which enable high-throughput characterization and quantification of biological molecules, and subsequent integration of genomics with these omics datasets, a field referred to as systems genomics. Single cell sequencing studies of the retina and choroid demonstrated that the majority of candidate AMD genes identified through genomic studies are expressed in non-neuronal cells, such as the retinal pigment epithelium (RPE), glia, myeloid and choroidal cells, highlighting that many different retinal and choroidal cell types contribute to the pathogenesis of AMD. Expression quantitative trait locus (eQTL) studies in retinal tissue have identified putative causal genes by demonstrating a genetic overlap between gene regulation and AMD risk. Linking genetic data to complement measurements in the systemic circulation has aided in understanding the effect of AMD-associated genetic variants in the complement system, and supports that protein QTL (pQTL) studies in plasma or serum samples may aid in understanding the effect of genetic variants and pinpointing causal genes in AMD. A recent epigenomic study fine-mapped AMD causal variants by determing regulatory regions in RPE cells differentiated from induced pluripotent stem cells (iPSC-RPE). Another approach that is being employed to pinpoint causal AMD genes is to produce synthetic DNA assemblons representing risk and protective haplotypes, which are then delivered to cellular or animal model systems. Pinpointing causal genes and understanding disease mechanisms is crucial for the next step towards clinical translation. Clinical trials targeting proteins encoded by the AMD-associated genomic loci C3, CFB, CFI, CFH, and ARMS2/HTRA1 are currently ongoing, and a phase III clinical trial for C3 inhibition recently showed a modest reduction of lesion growth in geographic atrophy. The EYERISK consortium recently developed a genetic test for AMD that allows genotyping of common and rare variants in AMD-associated genes. Polygenic risk scores (PRS) were applied to quantify AMD genetic risk, and may aid in predicting AMD progression. In conclusion, genomic studies represent a turning point in our exploration of AMD. The results of those studies now serve as a driving force for several clinical trials. Expanding to omics and systems genomics will further decipher function and causality from the associations that have been reported, and will enable the development of therapies that will lessen the burden of AMD.

DNA methylation and transcriptomic features are preserved throughout disease recurrence and chemoresistance in high grade serous ovarian cancers

Background

Little is known about the role of global DNA methylation in recurrence and chemoresistance of high grade serous ovarian cancer (HGSOC).

Methods

We performed whole genome bisulfite sequencing and transcriptome sequencing in 62 primary and recurrent tumors from 28 patients with stage III/IV HGSOC, of which 11 patients carried germline, pathogenic BRCA1 and/or BRCA2 mutations.

Results

Landscapes of genome-wide methylation (on average 24.2 million CpGs per tumor) and transcriptomes in primary and recurrent tumors showed extensive heterogeneity between patients but were highly preserved in tumors from the same patient. We identified significant differences in the burden of differentially methylated regions (DMRs) in tumors from BRCA1/2 compared to non-BRCA1/2 carriers (mean 659 DMRs and 388 DMRs in paired comparisons respectively). We identified overexpression of immune pathways in BRCA1/2 carriers compared to non-carriers, implicating an increased immune response in improved survival (P = 0.006) in these BRCA1/2 carriers.

Conclusion

These findings indicate methylome and gene expression programs established in the primary tumor are conserved throughout disease progression, even after extensive chemotherapy treatment, and that changes in methylation and gene expression are unlikely to serve as drivers for chemoresistance in HGSOC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02440-z.

QTLbase2: an enhanced catalog of human quantitative trait loci on extensive molecular phenotypes

Deciphering the fine-scale molecular mechanisms that shape the genetic effects at disease-associated loci from genome-wide association studies (GWAS) remains challenging. The key avenue is to identify the essential molecular phenotypes that mediate the causal variant and disease under particular biological conditions. Therefore, integrating GWAS signals with context-specific quantitative trait loci (QTLs) (such as different tissue/cell types, disease states, and perturbations) from extensive molecular phenotypes would present important strategies for full understanding of disease genetics. Via persistent curation and systematic data processing of large-scale human molecular trait QTLs (xQTLs), we updated our previous QTLbase database (now QTLbase2, http://mulinlab.org/qtlbase) to comprehensively analyze and visualize context-specific QTLs across 22 molecular phenotypes and over 95 tissue/cell types. Overall, the resource features the following major updates and novel functions: (i) 960 more genome-wide QTL summary statistics from 146 independent studies; (ii) new data for 10 previously uncompiled QTL types; (iii) variant query scope expanded to fit 195 QTL datasets based on whole-genome sequencing; (iv) supports filtering and comparison of QTLs for different biological conditions, such as stimulation types and disease states; (v) a new linkage disequilibrium viewer to facilitate variant prioritization across tissue/cell types and QTL types.

The FoxP1 gene regulates lung function, production of matrix metalloproteinases and inflammatory mediators, and viability of lung epithelia

Background

Genes involved in lung development may become dysregulated in adult life and contribute to the pathogenesis of lung diseases. Multiple genes regulate lung development, including Forkhead box protein P1-4 (FoxP1-4).

Methods

We examined the association between variants in the FoxP1-4 genes and lung function using data from a GWAS that included close to 400,000 individuals and 20 million SNPs.

Results

More than 100 variants in the FoxP1 gene, but none in the FoxP2-4 genes, are associated with lung function. The sentinel variant in the FoxP1 gene associated with FEV1 was rs1499894 (C > T), while the sentinel variant in the FoxP1 gene associated with FVC was rs35480566 (A > G). Those with the T allele instead of the C allele for rs1499894, or the G allele instead of the A allele for rs35480566 had increased FoxP1 mRNA levels in transcriptomic data, higher FEV1 and FVC, and reduced odds of being diagnosed with idiopathic pulmonary fibrosis. Further, knockdown of FoxP1 in lung epithelial cells by RNA interference led to increased mRNA levels for matrix metalloproteinases 1, 2, 3 and pro-inflammatory cytokines IL-6 & IL-8, as well as reduced cell viability after exposure to cigarette smoke—all processes implicated in the pathogenesis of COPD and IPF.

Conclusions

Our results suggest that the protein encoded by the FoxP1 gene may protect against the development of COPD and IPF. A causal role for FoxP1 in the pathogenesis of COPD and IPF may warrant further investigation, and FoxP1 may be a novel therapeutic target for these lung disorders.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12931-022-02213-4.

Uncovering the extensive trade-off between adaptive evolution and disease susceptibility

Favored mutations in the human genome may make the carriers adapt to changing environments and lifestyles but also susceptible to specific diseases. The scale and details of the trade-off between adaptive evolution and disease susceptibility are unclear because most favored mutations in different populations remain unidentified. As no statistical test can discriminate favored mutations from nearby hitchhiking neutral ones, we report a deep-learning network (DeepFavored) to integrate multiple statistical tests and divide identifying favored mutations into two subtasks. We identify favored mutations in three human populations and analyzed the correlation between favored/hitchhiking mutations and genome-wide association study (GWAS) sites. Both favored and hitchhiking neutral mutations are enriched in GWAS sites with population-specific features, and the enrichment and population specificity are prominent in genes in specific Gene Ontology (GO) terms. These provide evidence for extensive and population-specific trade-offs between adaptive evolution and disease susceptibility. The unveiled scale helps understand and investigate differences and diseases of humans.

A large-scale polygenic risk score analysis identified candidate proteins associated with anxiety, depression and neuroticism

Psychiatric disorders and neuroticism are closely associated with central nervous system, whose proper functioning depends on efficient protein renewal. This study aims to systematically analyze the association between anxiety / depression / neuroticism and each of the 439 proteins. 47,536 pQTLs of 439 proteins in brain, plasma and cerebrospinal fluid (CSF) were collected from recent genome-wide association study. Polygenic risk scores (PRS) of the 439 proteins were then calculated using the UK Biobank cohort, including 120,729 subjects of neuroticism, 255,354 subjects of anxiety and 316,513 subjects of depression. Pearson correlation analyses were performed to evaluate the correlation between each protein and each of the mental traits by using calculated PRSs as the instrumental variables of protein. In general population, six correlations were identified in plasma and CSF such as plasma protease C1 inhibitor (C1-INH) with neuroticism score (r = - 0.011, P = 2.56 × 10^- 9) in plasma, C1-INH with neuroticism score (r = -0.010, P = 3.09 × 10^- 8) in CSF, and ERBB1 with self-reported depression (r = - 0.012, P = 4.65 × 10^- 5) in CSF. C1-INH and ERBB1 may induce neuroticism and depression by affecting brain function and synaptic development. Gender subgroup analyses found that BST1 was correlated with neuroticism score in male CSF (r = - 0.011, P = 1.80 × 10^- 5), while CNTN2 was correlated with depression score in female brain (r = - 0.013, P = 6.43 × 10^- 4). BST1 and CNTN2 may be involved in nervous system metabolism and brain health. Six common candidate proteins were associated with all three traits (P < 0.05) and were confirmed in relevant proteomic studies, such as C1-INH in plasma, CNTN2 and MSP in the brain. Our results provide novel clues for revealing the roles of proteins in the development of anxiety, depression and neuroticism.

Multi-ancestry genome-wide association study of asthma exacerbations

BACKGROUND

Asthma exacerbations are a serious public health concern due to high healthcare resource utilization, work/school productivity loss, impact on quality of life, and risk of mortality. The genetic basis of asthma exacerbations has been studied in several populations, but no prior study has performed a multi-ancestry meta-analysis of genome-wide association studies (meta-GWAS) for this trait. We aimed to identify common genetic loci associated with asthma exacerbations across diverse populations and to assess their functional role in regulating DNA methylation and gene expression.

METHODS

A meta-GWAS of asthma exacerbations in 4989 Europeans, 2181 Hispanics/Latinos, 1250 Singaporean Chinese, and 972 African Americans analyzed 9.6 million genetic variants. Suggestively associated variants (p ≤ 5 × 10-5 ) were assessed for replication in 36,477 European and 1078 non-European asthma patients. Functional effects on DNA methylation were assessed in 595 Hispanic/Latino and African American asthma patients and in publicly available databases. The effect on gene expression was evaluated in silico.

RESULTS

One hundred and twenty-six independent variants were suggestively associated with asthma exacerbations in the discovery phase. Two variants independently replicated: rs12091010 located at vascular cell adhesion molecule-1/exostosin like glycosyltransferase-2 (VCAM1/EXTL2) (discovery: odds ratio (ORT allele ) = 0.82, p = 9.05 × 10-6 and replication: ORT allele  = 0.89, p = 5.35 × 10-3 ) and rs943126 from pantothenate kinase 1 (PANK1) (discovery: ORC allele  = 0.85, p = 3.10 × 10-5 and replication: ORC allele  = 0.89, p = 1.30 × 10-2 ). Both variants regulate gene expression of genes where they locate and DNA methylation levels of nearby genes in whole blood.

CONCLUSIONS

This multi-ancestry study revealed novel suggestive regulatory loci for asthma exacerbations located in genomic regions participating in inflammation and host defense.

Interaction between COMT Val158 Met polymorphism and childhood trauma predicts risk for depression in men

Depression is a disabling illness with complex etiology. While the Catechol-O-Methyltransferase (COMT) gene, in particular the functional Val¹⁵⁸ Met polymorphism, has been related to depression, the mechanisms underlying this gene-disease association are not completely understood. Therefore, we explore the association of COMT Val¹⁵⁸ Met polymorphism with depression as well as its interaction with childhood trauma in 1,136 young adults from a population-based study carried out in the city of Pelotas-Brazil. The diagnosis was performed through the Mini International Neuropsychiatric Interview 5.0 (MINI 5.0), and trauma was assessed with the childhood trauma questionnaire (CTQ). Total DNA was extracted and genotyped by real-time PCR and the QTLbase dataset was queried to perform large-scale quantitative trait locus (QTL) analysis. Our research showed no direct association between the Val¹⁵⁸ Met polymorphism and the diagnosis of depression (women: χ2=0.10, d=1, p=0.751 and men: χ2=0.003, df=1 p=0.956). However, the Met-allele of the Val¹⁵⁸ Met polymorphism modified the effect of childhood trauma in men [OR=2.58 (95% CI:1.05-6.29); p=0.038] conferring risk for depression only on those who suffer from trauma. The conditional effect from moderation analysis showed that trauma impacts the risk of depression only in men carrying the Met-allele (Effect: 0.9490, Standard Error (SE): 0.2570; p=0.0002). QTLbase and dataset for Val¹⁵⁸ Met polymorphism were consistent for markers that influence chromatin accessibility transcription capacity including histone methylation and acetylation. The changes caused in gene regulation by childhood trauma exposure and polymorphism may serve as evidence of the mechanism whereby the interaction increases susceptibility to this disorder in men.