eQTL Analysis Research Articles

Identification of Genetic Markers Linked to The Activity of Indoleamine 2,3-Dioxygenase and Kidney Function.

Indoleamine 2,3-dioxygenase (IDO) is a tryptophan-degrading enzyme belonging to the kynurenine pathway. IDO activity has been suggested as a potential biomarker for early diagnosis of chronic kidney disease (CKD). The aim of this study was to perform coincident association analysis to gain genetic insights into the correlation between IDO activity and CKD. This study evaluated the association between IDO activity and CKD using the Korea Association REsource (KARE) cohort. Logistic and linear regression were used to analyze CKD and quantitative phenotypes such as IDO and estimated glomerular filtration rate (eGFR). Our results identified 10 single nucleotide polymorphisms (SNPs) that were coincidently associated with both IDO and CKD (p < 0.001). Among them, rs6550842, rs77624055, and rs35651150 were selected as potential candidates after excluding SNPs with insufficient evidence for having an association with IDO or CKD. Further expression quantitative trait loci (eQTL) analysis for variants at selected loci showed that rs6550842 and rs35651150 significantly affected the expression of NKIRAS1 and SH2D4A genes in human tissues, respectively. Additionally, we highlighted that the NKIRAS1 and BMP6 genes were correlated with IDO activity and CKD through signaling pathways associated with inflammation. Our data suggest that NKIRAS1, SH2D4A, and BMP6 were potential causative genes affecting IDO activity and CKD through integrated analysis. Identifying these genes could aid in early detection and treatment by predicting the risk of CKD associated with IDO activity.

Multi-Ancestry eQTL Analysis of Human Brain Identifies Candidate Causal Variants for Psychiatric Disorders

Multi-Ancestry eQTL Analysis of Human Brain Identifies Candidate Causal Variants for Psychiatric Disorders

Genetic variants of SOS2, MAP2K1 and RASGRF2 in the RAS pathway genes predict survival of HBV-related hepatocellular carcinoma patients.

The RAS pathway participates in the cascade of proliferation and cell division process, and the activated RAS pathway can lead to tumorigenesis including hepatocellular carcinoma (HCC). However, few studies have explored the effects of genetic variants in the RAS pathway-related genes on the survival of patients with HBV-related HCC. In the present study, we assessed the associations between 11,658 single-nucleotide polymorphisms (SNPs) in 62 RAS pathway genes and the overall survival (OS) of 866 HBV-related HCC individuals, which were randomly split (1:1) into discovery and validation datasets. As a result, three potentially functional SNPs were identified, based on multivariable cox proportional hazards regression analyses, in SOS Ras/Rho guanine nucleotide exchange factor 2 (SOS2, rs4632055 A > G), Ras protein-specific guanine nucleotide releasing factor 2 (RASGRF2, rs26418A > G) and mitogen-activated protein kinase 1 (MAP2K1,rs57120695 C > T), which were significantly and independently associated with OS of HBV-related HCC patients [adjusted hazards ratios (HRs) of 1.42, 1.32 and 1.50, respectively; 95% confidence intervals (CI), 1.14 to 1.76, 1.15 to 1.53 and 1.15 to 1.97, respectively; P = 0.001, < 0.001 and 0.003, respectively]. Additionally, the joint effects as the unfavorable genotypes of these three SNPs showed a significant association with the poor survival of HCC (trend test P < 0.001). The expression quantitative trait loci (eQTL) analysis further revealed that the rs4632055 G allele and the rs26418 A allele were associated with lower mRNA expression levels of SOS2 and RASGRF2, respectively. Collectively, these potentially functional SNPs of RASGRF2, SOS2 and M2PAK1 may become potential prognostic biomarkers for HBV-related HCC after hepatectomy.

Abstract 1404: A functional genetic polymorphism of mitochondrion-associated DNAJA3 in nonalcoholic steatohepatitis-related hepatocellular carcinoma

Abstract Nonalcoholic fatty liver disease (NAFLD) is the most common form of liver disease worldwide and may progress to nonalcoholic steatohepatitis (NASH). Although emerging evidence suggests that NASH would be an important driver of hepatocellular carcinoma (HCC), molecularly what drives the transition of NAFLD to NASH and then HCC has not been clearly defined. Since NASH is a mitochondrial disease, we aimed to determine genetic variants associated with mitochondrial dysfunction linked to NASH and NASH-related HCC. We performed an integrative genomics analysis by analyzing a large-scale genome-wide association study of the body fat distribution-related variants in 344,369 individuals from multiple ancestries. And we linked these results to the top 24 significant variants to NASH and NASH-related HCC in 1,009 individuals from the NCI-UMD cohort. We also performed the expression quantitative trait loci (eQTL) analysis on candidate variants. We analyzed the 3D genome architecture especially chromatin topology based on UCSC Genome Browser that influences the expression of candidate variants. As a result, we identified rs3747579-TT, an eQTL to a mitochondrial chaperone Hsp40, also known as DNAJA3 rather than a coding variant, as a risk allele linked to NASH-related HCC. HCC patients with rs3747579-TT variant had a poor prognosis compared to HCC patients with rs3747579-CT or CC variant. We confirmed that HCC with rs3747579-TT had a reduced expression of DNAJA3, encoding a protein that maintains mitochondrial homeostasis. Consistently, HCC patients with low DNAJA3 expression had a poor prognosis. We found a specific site based on chromatin topology that may influence rs3747579 to regulate the expression of DNAJA3. In summary, we identified rs3747579-TT as a functional polymorphic variant of DNAJA3 that may be responsible for the development of NASH and NASH-related HCC. Citation Format: Yuto Shiode, Ching-Wen Chang, Xin Wei Wang. A functional genetic polymorphism of mitochondrion-associated DNAJA3 in nonalcoholic steatohepatitis-related hepatocellular carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1404.

An Association Study of Germline Variants in Bladder Cancer-Related Genes with the Prognosis of Non-Muscle Invasive Bladder Cancer.

Various germline genetic variants are associated with the prognosis of non-muscle invasive bladder cancer (NMIBC). Germline variants in genes frequently somatically mutated in bladder cancer have not been studied thoroughly in relation to risk of recurrence or progression in NMIBC. To identify germline DNA variants in bladder carcinogenesis-related genes associated with recurrence or progression in NMIBC. We analysed associations between single-nucleotide polymorphisms (SNPs) and NMIBC recurrence and progression using data from the Nijmegen Bladder Cancer Study (NBCS, 1,443 patients). We included 5,053 SNPs within 46 genes known to have mutation, overexpression or amplification in bladder cancer. We included all recurrences in the statistical analysis and performed both single variant analysis and gene-based analysis. SNPs and genes that showed significant or suggestive association (false discovery rate P value < 20%) were followed-up in independent cohorts for replication analysis, through eQTL analysis and tests for association of tumour expression levels with NMIBC recurrence and progression. Single variant analysis showed no statistically significant associations with recurrence or progression. In gene-based analysis, the aggregate effect of the 25 SNPs in the Cyclin D1 gene (CCND1) was statistically significantly associated with NMIBC recurrence (Punadj = 0.001, PFDR = 0.046), but not with progression (Punadj = 0.17, PFDR = 0.54). Validation analysis in independent cohorts did not confirm the association of CCND1 with NMIBC recurrence. We could not identify reproducible associations between common germline variants in bladder carcinogenesis-related genes and NMIBC recurrence or progression.

Cross-phenotype association analysis of gastric cancer: in-silico functional annotation based on the disease-gene network.

A gene or variant has pleiotropic effects, and genetic variant identification across multiple phenotypes can provide a comprehensive understanding of biological pathways shared among different diseases or phenotypes. Discovery of genetic loci associated with multiple diseases can simultaneously support general interventions. Several meta-analyses have shown genetic associations with gastric cancer (GC); however, no study has identified associations with other phenotypes using this approach. Here, we applied disease network analysis and gene-based analysis (GBA) to examine genetic variants linked to GC and simultaneously associated with other phenotypes. We conducted a single-nucleotide polymorphism (SNP) level meta-analysis and GBA through a systematic genome-wide association study (GWAS) linked to GC, to integrate published results for the SNP variants and group them into major GC-associated genes. We then performed disease network and expression quantitative trait loci (eQTL) analyses to evaluate cross-phenotype associations and expression levels of GC-related genes. Seven genes (MTX1, GBAP1, MUC1, TRIM46, THBS3, PSCA, and ABO) were associated with GC as well as blood urea nitrogen (BUN), glomerular filtration rate (GFR), and uric acid (UA). In addition, 17 SNPs regulated the expression of genes located on 1q22, 24 SNPs regulated the expression of PSCA on 8q24.3, and rs7849820 regulated the expression of ABO on 9q34.2. Furthermore, rs1057941 and rs2294008 had the highest posterior causal probabilities of being a causal candidate SNP in 1q22, and 8q24.3, respectively. These findings identified seven GC-associated genes exhibiting a cross-association with GFR, BUN, and UA.

Integrative Genomic Analysis of m6a-SNPs Identifies Potential Functional Variants Associated with Alzheimer's Disease.

Alzheimer's disease (AD) is a neurodegenerative disorder that affects 35 million people worldwide. However, no potential therapeutics currently are available for AD because of the multiple factors involved in it, such as regulatory factors with their candidate genes, factors associated with the expression levels of its corresponding genes, and many others. To date, 29 novel loci from GWAS have been reported for AD by the Psychiatric Genomics Consortium (PGC2). Nevertheless, the main challenge of the post-GWAS era, namely to detect significant variants of the target disease, has not been conducted for AD. N6-methyladenosine (m6a) is reported as the most prevalent mRNA modification that exists in eukaryotes and that influences mRNA nuclear export, translation, splicing, and the stability of mRNA. Furthermore, studies have also reported m6a's association with neurogenesis and brain development. We carried out an integrative genomic analysis of AD variants from GWAS and m6a-SNPs from m6AVAR to identify the effects of m6a-SNPs on AD and identified the significant variants using the statistically significance value (p-value <0.05). The cis-regularity variants with their corresponding genes and their influence on gene expression in the gene expression profiles of AD patients were determined, and showed 1458 potential m6a-SNPs (based on p-value <0.05) associated with AD. eQTL analysis showed that 258 m6a-SNPs had cis-eQTL signals that overlapped with six significant differentially expressed genes based on p-value <0.05 in two datasets of AD gene expression profiles. A follow-up study to elucidate the impact of our identified m6a-SNPs in the experimental study would validate our findings for AD, which would contribute to the etiology of AD.

Distinct gene dysregulation patterns herald precision medicine potentiality in systemic lupus erythematosus

ObjectivesWe aimed at investigating the whole-blood transcriptome, expression quantitative trait loci (eQTLs), and levels of selected serological markers in patients with SLE versus healthy controls (HC) to gain insight into pathogenesis and identify drug targets. MethodsWe analyzed differentially expressed genes (DEGs) and dysregulated gene modules in a cohort of 350 SLE patients and 497 HC from the European PRECISESADS project (NTC02890121), split into a discovery (60%) and a replication (40%) set. Replicated DEGs qualified for eQTL, pathway enrichment, regulatory network, and druggability analysis. For validation purposes, a separate gene module analysis was performed in an independent cohort (GSE88887). ResultsAnalysis of 521 replicated DEGs identified multiple enriched interferon signaling pathways through Reactome. Gene module analysis yielded 18 replicated gene modules in SLE patients, including 11 gene modules that were validated in GSE88887. Three distinct gene module clusters were defined i.e., “interferon/plasma cells”, “inflammation”, and “lymphocyte signaling”. Predominant downregulation of the lymphocyte signaling cluster denoted renal activity. By contrast, upregulation of interferon-related genes indicated hematological activity and vasculitis. Druggability analysis revealed several potential drugs interfering with dysregulated genes within the “interferon” and “PLK1 signaling events” modules. STAT1 was identified as the chief regulator in the most enriched signaling molecule network. Drugs annotated to 15 DEGs associated with cis-eQTLs included bortezomib for its ability to modulate CTSL activity. Belimumab was annotated to TNFSF13B (BAFF) and daratumumab was annotated to CD38 among the remaining replicated DEGs. ConclusionsModulation of interferon, STAT1, PLK1, B and plasma cell signatures showed promise as viable approaches to treat SLE, pointing to their importance in SLE pathogenesis.

Genetic analyses of chr11p15.5 region identify MUC5AC-MUC5B associated with asthma-related phenotypes

Objective Genome-wide association studies (GWASs) have identified single nucleotide polymorphisms (SNPs) in chr11p15.5 region associated with asthma and idiopathic interstitial pneumonias (IIPs). We sought to identify functional genes for asthma by combining SNPs and mRNA expression in bronchial epithelial cells (BEC) in the Severe Asthma Research Program (SARP). Methods Correlation analyses of mRNA expression of six candidate genes (AP2A2, MUC6, MUC2, MUC5AC, MUC5B, and TOLLIP) and asthma phenotypes were performed in the longitudinal cohort (n = 156) with RNAseq in BEC, and replicated in the cross-sectional cohort (n = 155). eQTL (n = 114) and genetic association analysis of asthma severity (426 severe vs. 531 non-severe asthma) were performed, and compared with previously published GWASs of IIPs and asthma. Results Higher expression of AP2A2 and MUC5AC and lower expression of MUC5B in BEC were correlated with asthma, asthma exacerbations, and T2 biomarkers (P < 0.01). SNPs associated with asthma and IIPs in previous GWASs were eQTL SNPs for MUC5AC, MUC5B, or TOLLIP, however, they were not in strong linkage disequilibrium. The risk alleles for asthma or protective alleles for IIPs were associated with higher expression of MUC5AC and lower expression of MUC5B. rs11603634, rs12788104, and rs28415845 associated with moderate-to-severe asthma or adult onset asthma in previous GWASs were not associated with asthma severity (P > 0.8). Conclusions SNPs associated with asthma in chr11p15.5 region are not associated with asthma severity neither with IIPs. Higher expression of MUC5AC and lower expression of MUC5B are risk for asthma but protective for IIPs.

Genome-wide identification of RNA modification-related single nucleotide polymorphisms associated with rheumatoid arthritis

BackgroundRNA modification plays important roles in many biological processes, such as gene expression control. The aim of this study was to identify single nucleotide polymorphisms related to RNA modification (RNAm-SNPs) for rheumatoid arthritis (RA) as putative functional variants.MethodsWe examined the association of RNAm-SNPs with RA in summary data from a genome-wide association study of 19,234 RA cases and 61,565 controls. We performed eQTL and pQTL analyses for the RNAm-SNPs to find associated gene expression and protein levels. Furthermore, we examined the associations of gene expression and circulating protein levels with RA using two-sample Mendelian randomization analysis methods.ResultsA total of 160 RNAm-SNPs related to m6A, m1A, A-to-I, m7G, m5C, m5U and m6Am modifications were identified to be significantly associated with RA. These RNAm-SNPs were located in 62 protein-coding genes, which were significantly enriched in immune-related pathways. RNAm-SNPs in important RA susceptibility genes, such as PADI2, SPRED2, PLCL2, HLA-A, HLA-B, HLA-DRB1, HLA-DPB1, TRAF1 and TXNDC11, were identified. Most of these RNAm-SNPs showed eQTL effects, and the expression levels of 26 of the modifiable genes (e.g., PADI2, TRAF1, HLA-A, HLA-DRB1, HLA-DPB1 and HLA-B) in blood cells were associated with RA. Circulating protein levels, such as CFB, GZMA, HLA-DQA2, IL21, LRPAP1 and TFF3, were affected by RNAm-SNPs and were associated with RA.ConclusionThe present study identified RNAm-SNPs in the reported RA susceptibility genes and suggested that RNAm-SNPs may affect RA risk by affecting the expression levels of corresponding genes and proteins.

Tuberculosis severity associates with variants and eQTLs related to vascular biology and infection-induced inflammation.

Tuberculosis (TB) remains a major public health problem globally, even compared to COVID-19. Genome-wide studies have failed to discover genes that explain a large proportion of genetic risk for adult pulmonary TB, and even fewer have examined genetic factors underlying TB severity, an intermediate trait impacting disease experience, quality of life, and risk of mortality. No prior severity analyses used a genome-wide approach. As part of our ongoing household contact study in Kampala, Uganda, we conducted a genome-wide association study (GWAS) of TB severity measured by TBScore, in two independent cohorts of culture-confirmed adult TB cases (n = 149 and n = 179). We identified 3 SNPs (P<1.0 x 10-7) including one on chromosome 5, rs1848553, that was GWAS significant (meta-analysis p = 2.97x10-8). All three SNPs are in introns of RGS7BP and have effect sizes corresponding to clinically meaningful reductions in disease severity. RGS7BP is highly expressed in blood vessels and plays a role in infectious disease pathogenesis. Other genes with suggestive associations defined gene sets involved in platelet homeostasis and transport of organic anions. To explore functional implications of the TB severity-associated variants, we conducted eQTL analyses using expression data from Mtb-stimulated monocyte-derived macrophages. A single variant (rs2976562) associated with monocyte SLA expression (p = 0.03) and subsequent analyses indicated that SLA downregulation following MTB stimulation associated with increased TB severity. Src Like Adaptor (SLAP-1), encoded by SLA, is highly expressed in immune cells and negatively regulates T cell receptor signaling, providing a potential mechanistic link to TB severity. These analyses reveal new insights into the genetics of TB severity with regulation of platelet homeostasis and vascular biology being central to consequences for active TB patients. This analysis also reveals genes that regulate inflammation can lead to differences in severity. Our findings provide an important step in improving TB patient outcomes.

Meta-analysis of erosive hand osteoarthritis identifies four common variants that associate with relatively large effect

ObjectivesErosive hand osteoarthritis (EHOA) is a severe subset of hand osteoarthritis (OA). It is unclear if EHOA is genetically different from other forms of OA. Sequence variants at ten loci...

Abstract P550: Genome-Wide Association Analysis Identified Genetic Variants Associated With Cardiac Structure and Function in the Hispanic Community Health Study/Study of Latinos

Introduction: The role of genetic variants on cardiac structure and function among the Hispanic/Latino population is unclear. Hypothesis: We hypothesize that common genetic variants are associated with cardiac structure and function measures in the Hispanic/Latino population. Methods: We performed genome-wide association analyses of nine echo measures among individuals from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL), including: left ventricular end-systolic volume (LVESV), LV end-diastolic volume, mean wall thickness, LV mass index, LV ejection fraction (LVEF), stroke volume (SV), septal early mitral annular velocity, ratio of early ventricular filling peak velocity to early mitral annular velocity, and left atrial volume index. We further conducted enrichment analysis to locate candidate gene pathways and performed replication using data from UK BioBank (UKB) and BioBank of Japan (BBJ). Results: We included 4471 Hispanic/Latino participants from two sub-cohorts of HCHS/SOL (mean age 55 and 59, 68% female). Fourteen previously unreported variants at 3 loci were associated with LVEF, SV, or LVESV (p-value &lt; 5х10 -8 , Figure ). Twelve of the 14 variants were associated with LVEF, located on or nearby ncRNA gene LOC105374813 , as well as 1 SNP on gene IPO5 with SV and 1 SNP on gene PGBD5 with LVESV. There were 199 variants that reached suggestive significance (p&lt;5х10 -7 ) for nine measures, with most SNPs had effect directions consistent compared to UKB and BBJ findings. FUMA enrichment analysis showed SNP associated with SV on IPO5 had higher expression in kidney, liver, and heart tissue compared to other tissues. For variants associated with LVEF, eQTL analysis with heart tissue showed that gene RP11-342K6.2 and GCFC2 had chromatin interaction in the region. Conclusions: We identified 14 loci associated with cardiac structure and function among Hispanics/Latinos, providing supporting evidence of genetic contribution to cardiac dysfunction in Hispanic/Latino populations.

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Pitfalls and opportunities for applying latent variables in single-cell eQTL analyses

Using latent variables in gene expression data can help correct unobserved confounders and increase statistical power for expression quantitative trait Loci (eQTL) detection. The probabilistic estimation of expression residuals (PEER) and principal component analysis (PCA) are widely used methods that can remove unwanted variation and improve eQTL discovery power in bulk RNA-seq analysis. However, their performance has not been evaluated extensively in single-cell eQTL analysis, especially for different cell types. Potential challenges arise due to the structure of single-cell RNA-seq data, including sparsity, skewness, and mean-variance relationship. Here, we show by a series of analyses that PEER and PCA require additional quality control and data transformation steps on the pseudo-bulk matrix to obtain valid latent variables; otherwise, it can result in highly correlated factors (Pearson's correlation r = 0.63 ~ 0.99). Incorporating valid PFs/PCs in the eQTL association model would identify 1.7 ~ 13.3% more eGenes. Sensitivity analysis showed that the pattern of change between the number of eGenes detected and fitted PFs/PCs varied significantly in different cell types. In addition, using highly variable genes to generate latent variables could achieve similar eGenes discovery power as using all genes but save considerable computational resources (~ 6.2-fold faster).

Association of Plant-Based and High-Protein Diets with a Lower Obesity Risk Defined by Fat Mass in Middle-Aged and Elderly Persons with a High Genetic Risk of Obesity

Obesity has become a severe public health challenge globally. The present study aimed to identify separate and interactive dietary, genetic, and other factors that increase the risk of obesity as measured by body fat (BF) mass. We utilized a genome-wide association study to identify genetic variants associated with high fat mass (obesity; n = 10,502) and combined them to generate polygenic risk scores (PRS) of genetic variants interacting with each other in adults aged over 40 while excluding body-fat-related diseases in a city-hospital-based cohort (n = 53,828). It was validated in Ansan/Ansung plus rural cohorts (n = 13,007). We then evaluated dietary and lifestyle factors in subjects to assess what factors might help overcome a genetic propensity for higher BF. The three-SNP model included brain-derived neurotrophic factor (BDNF)_rs6265, fat-mass- and obesity-associated protein (FTO)_rs1421085, and SEC16B_rs509325. The genes with the minor alleles of ADCY3_rs6545790 and BAIAP2_rs35867081 increased their gene expression in the visceral and subcutaneous adipocytes, but their gene expression decreased in the hypothalamus in eQTL analysis. In the three-SNP model, the PRS was associated with BF mass by 1.408 and 1.396 times after adjusting covariates 1 (age, gender, survey year, residence area, education, and income) and 2 (covariates in model 1 plus energy intake, alcohol intake, regular exercise, and smoking status), respectively. However, when separating subjects by PRS of the three-SNP model, a plant-based diet was the most significant factor associated with low BF, followed by high-protein diets and lower energy intakes. They could offset the effects of high genetic risk for high BF. In conclusion, modulating nutrient intakes might overcome a high genetic risk for obesity. Dietary choices favoring more plant-based and higher-protein foods might help prevent increased BF in Asians and potentially people of other ethnicities with high polygenetic risk scores.

Genetic Underpinnings of Brain Structural Connectome for Young Adults

With distinct advantages in power over behavioral phenotypes, brain imaging traits have become emerging endophenotypes to dissect molecular contributions to behaviors and neuropsychiatric illnesses. Among different imaging features, brain structural connectivity (i.e., structural connectome) which summarizes the anatomical connections between different brain regions is one of the most cutting edge while under-investigated traits; and the genetic influence on the structural connectome variation remains highly elusive. Relying on a landmark imaging genetics study for young adults, we develop a biologically plausible brain network response shrinkage model to comprehensively characterize the relationship between high dimensional genetic variants and the structural connectome phenotype. Under a unified Bayesian framework, we accommodate the topology of brain network and biological architecture within the genome; and eventually establish a mechanistic mapping between genetic biomarkers and the associated brain sub-network units. An efficient expectation-maximization algorithm is developed to estimate the model and ensure computing feasibility. In the application to the Human Connectome Project Young Adult (HCP-YA) data, we establish the genetic underpinnings which are highly interpretable under functional annotation and brain tissue eQTL analysis, for the brain white matter tracts connecting the hippocampus and two cerebral hemispheres. We also show the superiority of our method in extensive simulations. Supplementary materials for this article are available online.

Genetic associations between bipolar disorder and brain structural phenotypes.

Patients with bipolar disorder (BD) and their first-degree relatives exhibit alterations in brain volume and cortical structure, whereas the underlying genetic mechanisms remain unclear. In this study, based on the published genome-wide association studies (GWAS), the extent of polygenic overlap between BD and 15 brain structural phenotypes was investigated using linkage disequilibrium score regression and MiXeR tool, and the shared genomic loci were discovered by conjunctional false discovery rate (conjFDR) and expression quantitative trait loci (eQTL) analyses. MiXeR estimated the overall measure of polygenic overlap between BD and brain structural phenotypes as 4-53% on a 0-100% scale (as quantified by the Dice coefficient). Subsequent conjFDR analyses identified 54 independent loci (71 risk single-nucleotide polymorphisms) jointly associated with BD and brain structural phenotypes with a conjFDR < 0.05, among which 33 were novel that had not been reported in the previous BD GWAS. Follow-up eQTL analyses in respective brain regions both confirmed well-known risk genes (e.g. CACNA1C, NEK4, GNL3, MAPK3) and discovered novel risk genes (e.g. LIMK2 and CAMK2N2). This study indicates a substantial shared genetic basis between BD and brain structural phenotypes, and provides novel insights into the developmental origin of BD and related biological mechanisms.

Novel Genetic Variants in TP37, PIK3R1, CALM1, and PLCG2 of the Neurotrophin Signaling Pathway Are Associated with the Progression from Mild Cognitive Impairment to Alzheimer's Disease.

Alzheimer's disease (AD) is a common neurodegenerative disease and mild cognitive impairment (MCI) is considered as the prodromal stage of AD. Previous studies showed that changes in the neurotrophin signaling pathway could lead to cognitive decline in AD. However, the association of single nucleotide polymorphisms (SNPs) in genes that are involved in this pathway with AD progression from MCI remains unclear. We investigated the associations between SNPs involved in the neurotrophin signaling pathway with AD progression. We performed single-locus analysis to identify neurotrophin-signaling-related SNPs associated with the AD progression using 767 patients from the Alzheimer's Disease Neuroimaging Initiative study and 1,373 patients from the National Alzheimer's Coordinating Center study. We constructed polygenic risk scores (PRSs) using the identified independent non-APOE SNPs and evaluated its prediction performance on AD progression. We identified 25 SNPs significantly associated with AD progression with Bayesian false-discovery probability ≤0.8. Based on the linkage disequilibrium clumping and expression quantitative trait loci analysis, we found 6 potentially functional SNPs that were associated with AD progression independently. The PRS analysis quantified the combined effects of these SNPs on longitudinal cognitive assessments and biomarkers from cerebrospinal fluid and neuroimaging. The addition of PRSs to the prediction model for 3-year progression to AD from MCI significantly increased the predictive accuracy. Genetic variants in the specific genes of the neurotrophin signaling pathway are predictors of AD progression. eQTL analysis supports that these SNPs regulate expression of key genes involved in the neurotrophin signaling pathway.

Genome-Wide Association Study on Longitudinal Change in Fasting Plasma Glucose in Korean Population.

Genome-wide association studies (GWAS) on type 2 diabetes mellitus (T2DM) have identified more than 400 distinct genetic loci associated with diabetes and nearly 120 loci for fasting plasma glucose (FPG) and fasting insulin level to date. However, genetic risk factors for the longitudinal deterioration of FPG have not been thoroughly evaluated. We aimed to identify genetic variants associated with longitudinal change of FPG over time. We used two prospective cohorts in Korean population, which included a total of 10,528 individuals without T2DM. GWAS of repeated measure of FPG using linear mixed model was performed to investigate the interaction of genetic variants and time, and meta-analysis was conducted. Genome-wide complex trait analysis was used for heritability calculation. In addition, expression quantitative trait loci (eQTL) analysis was performed using the Genotype-Tissue Expression project. A small portion (4%) of the genome-wide single nucleotide polymorphism (SNP) interaction with time explained the total phenotypic variance of longitudinal change in FPG. A total of four known genetic variants of FPG were associated with repeated measure of FPG levels. One SNP (rs11187850) showed a genome-wide significant association for genetic interaction with time. The variant is an eQTL for NOC3 like DNA replication regulator (NOC3L) gene in pancreas and adipose tissue. Furthermore, NOC3L is also differentially expressed in pancreatic β-cells between subjects with or without T2DM. However, this variant was not associated with increased risk of T2DM nor elevated FPG level. We identified rs11187850, which is an eQTL of NOC3L, to be associated with longitudinal change of FPG in Korean population.