Genetic study of complex diseases in the post-GWAS era

The relationship between innate/adaptive immunity and gastrointestinal cancer : a multi-omics Mendelian randomization study

BACKGROUND

Innate/adaptive immunity is the key to anti-tumor therapy. However, its causal relationship to Gastrointestinal (GI) cancer remains unclear.

METHODS

Immunity genes were extracted from the MSigDB database. The Genome-wide association studies (GWAS) summary data of GI cancer were integrated with expression quantitative trait loci (eQTL) and DNA methylation quantitative trait loci (mQTL) associated with genes. Summary-data-based Mendelian randomization (SMR) and co-localization analysis were used to reveal causal relationships between genes and GI cancer. Two-sample MR analysis was used for sensitivity analysis. Single cell analysis clarified the enrichment of genes.

RESULTS

Three-step SMR analysis showed that a putative mechanism, cg17294865 CpG site regulating HLA-DRA expression was negatively associated with gastric cancer risk. HLA-DRA was significantly differentially expressed in monocyte/macrophage and myeloid cells in gastric cancer.

CONCLUSION

This study provides evidence that upregulating the expression level of HLA-DRA can reduce the risk of gastric cancer.

Artificial Intelligence and Computational Biology in Gene Therapy: A Review

One of the trending fields in almost all areas of science and technology is artificial intelligence. Computational biology and artificial intelligence can help gene therapy in many steps including: gene identification, gene editing, vector design, development of new macromolecules and modeling of gene delivery. There are various tools used by computational biology and artificial intelligence in this field, such as genomics, transcriptomic and proteomics data analysis, machine learning algorithms and molecular interaction studies. These tools can introduce new gene targets, novel vectors, optimized experiment conditions, predict the outcomes and suggest the best solutions to avoid undesired immune responses following gene therapy treatment.

The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.

Potentially functional genetic variants in RPS6KA4 and MAP2K5 in the MAPK signaling pathway predict HBV-related hepatocellular carcinoma survival

Hepatocellular carcinoma (HCC) ranks the third leading cause of cancer deaths with a dismal 5-year survival rate. The mitogen-activated protein kinase (MAPK) signaling pathway is abnormally activated in HCC to promote growth and aggressive metastatic potential of cancer cells. Therefore, genetic variants in the MAPK signaling pathway may serve as potential predictors of Hepatitis B virus (HBV)-related HCC survival. In the present study, we performed a two-stage survival analysis to evaluate the associations between 10,912 single nucleotide polymorphisms (SNPs) in 79 MAPK signaling pathway genes and the overall survival (OS) of 866 HBV-related HCC patients, followed by functional annotation. In combined datasets, we identified two novel and potential functional SNPs (RPS6KA4 rs600377 T>G and MAP2K5 rs17300363 A>C) as prognostic factors for HBV-related HCC, with adjusted allelic hazards ratios of 1.24 (95% confidence interval [CI] = 1.05-1.46, p = 0.010) and 1.48 (1.15-1.91, p = 0.001), respectively. Furthermore, their combined risk genotypes also predicted a poor survival in a dose-response manner in the combined data set (Ptrend  < 0.001). Additional functional analysis showed that RPS6KA4 rs600377 G and MAP2K5 rs17300363 C alleles were associated with elevated mRNA expression levels of the corresponding genes in normal tissues. These results provide new insights into the role of genetic variants in the MAPK signaling pathway genes in HBV-related HCC survival.

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.

Novel single nucleotide polymorphisms in gestational diabetes mellitus

The association between gestational diabetes mellitus (GDM) and single nucleotide polymorphisms (SNPs) has attracted global research attention. Exploring SNPs can help us further understand the pathogenesis of GDM, predict the risk of GDM, and guide the management of GDM patients. In this review, we summarized the studies on the association between SNPs and GDM, focusing on novel SNPs published in the last 10 years. The SNPs identified to be associated with GDM included HMG20A (rs7178572), CDKAL1 (rs7756992, rs7754840, and rs7747752), ADIPOQ (rs266729 and rs17300539), MTHFR (rs1801133), IL10 (rs3021094), CDKN2B (rs1063192), and TRPM5 (rs35197079). However, the role of SNPs in the prediction, diagnosis, treatment, and prognosis of GDM, as a polygenic disease, needs to be further explored in multiple ethnic populations.

SNPs in 3'UTR miRNA Target Sequences Associated with Individual Drug Susceptibility

The complementary interaction of microRNAs (miRNAs) with their binding sites in the 3'untranslated regions (3'UTRs) of target gene mRNAs represses translation, playing a leading role in gene expression control. MiRNA recognition elements (MREs) in the 3'UTRs of genes often contain single nucleotide polymorphisms (SNPs), which can change the binding affinity for target miRNAs leading to dysregulated gene expression. Accumulated data suggest that these SNPs can be associated with various human pathologies (cancer, diabetes, neuropsychiatric disorders, and cardiovascular diseases) by disturbing the interaction of miRNAs with their MREs located in mRNA 3'UTRs. Numerous data show the role of SNPs in 3'UTR MREs in individual drug susceptibility and drug resistance mechanisms. In this review, we brief the data on such SNPs focusing on the most rigorously proven cases. Some SNPs belong to conventional genes from the drug-metabolizing system (in particular, the genes coding for cytochromes P450 (CYP 450), phase II enzymes (SULT1A1 and UGT1A), and ABCB3 transporter and their expression regulators (PXR and GATA4)). Other examples of SNPs are related to the genes involved in DNA repair, RNA editing, and specific drug metabolisms. We discuss the gene-by-gene studies and genome-wide approaches utilized or potentially utilizable to detect the MRE SNPs associated with individual response to drugs.

Progress Toward a Multiomic Understanding of Traumatic Brain Injury: A Review

Traumatic brain injury (TBI) is not a single disease state but describes an array of conditions associated with insult or injury to the brain. While some individuals with TBI recover within a few days or months, others present with persistent symptoms that can cause disability, neuropsychological trauma, and even death. Understanding, diagnosing, and treating TBI is extremely complex for many reasons, including the variable biomechanics of head impact, differences in severity and location of injury, and individual patient characteristics. Because of these confounding factors, the development of reliable diagnostics and targeted treatments for brain injury remains elusive. We argue that the development of effective diagnostic and therapeutic strategies for TBI requires a deep understanding of human neurophysiology at the molecular level and that the framework of multiomics may provide some effective solutions for the diagnosis and treatment of this challenging condition. To this end, we present here a comprehensive review of TBI biomarker candidates from across the multiomic disciplines and compare them with known signatures associated with other neuropsychological conditions, including Alzheimer’s disease and Parkinson’s disease. We believe that this integrated view will facilitate a deeper understanding of the pathophysiology of TBI and its potential links to other neurological diseases.

Genetic screening in patients of Meige syndrome and blepharospasm

OBJECTIVE

Meige syndrome (MS) is cranial dystonia, including bilateral eyelid spasms (blepharospasm; BSP) and involuntary movements of the jaw muscles (oromandibular dystonia; OMD). Up to now, the pathogenic genes of MS and BSP are still unclear.

METHODS

We performed Sanger sequencing of GNAL, TOR1A, TOR2A, THAP1, and REEP4 exons on 78 patients, including 53 BSP and 25 MS and 96 healthy controls.

RESULTS

c.845G > C[R282P] of TOR1A, c.629delC[p.Gly210AlafsTer60] of TOR2A, c.1322A > G[N441S] of GNAL, c.446G > A[R149Q], and c.649C > T[R217C] of REEP4 were identified and predicated as deleterious probably damaging variants. Three potential alterations of splicing variants of TOR1A and TOR2A were identified in patients. The frequencies of TOR1A rs1435566780 and THAP1 rs545930392 were higher in patients than in controls.

CONCLUSIONS

TOR1A rs1435566780 (c.*16G > C(G > A)) and THAP1 rs545930392 (c.192G > A[K64K]) may contribute to the etiology of MS and BSP. Other identified rare mutations predicted as deleterious probably damaging need further confirmation. Larger MS and BSP cohorts and functional studies will need to be performed further to elucidate the association between these genes and the diseases.

The osteoporosis risk variant rs9820407 at 3p22.1 acts as an allele-specific enhancer to regulate CTNNB1 expression by long-range chromatin loop formation

Previous powerful genome-wide association studies (GWASs) and whole-genome sequencing have identified multiple single-nucleotide polymorphisms (SNPs) located over 69 kb upstream of CTNNB1 at 3p22.1 locus associated with osteoporosis. The CTNNB1 gene encodes β-catenin that is an integral part of adherens junctions and the primary mediator of the canonical Wnt signaling pathway. The causal variants and underlying molecular mechanisms of the osteoporosis susceptibility locus 3p22.1 remains unknown. Through comprehensive computational analyses, including expression quantitative trait locus (eQTL), high-throughput chromatin interaction (Hi-C), epigenomic and functional annotation, four enhancer SNPs (rs9820407, rs9878224, rs454690 and rs9832204) were prioritized as potential causal SNPs at 3p22.1 for osteoporosis. Rs9820407 displayed the strongest enhancer activity in dual-luciferase assays. Specifically, the minor rs9820407-A can preferentially bind transcription factor FOXC1, elevate the enhancer activity and increase CTNNB1 expression. The architectural protein CTCF was presumably involved in long-range chromatin interaction between rs9820407 and CTNNB1. Our study provided a mechanistic insight into how noncoding enhancer SNP rs9820407 distally regulates CTNNB1 expression and modulates osteoporosis risk.

Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications

Osteoporosis is a common skeletal disease, affecting ~200 million people around the world. As a complex disease, osteoporosis is influenced by many factors, including diet (e.g. calcium and protein intake), physical activity, endocrine status, coexisting diseases and genetic factors. In this review, we first summarize the discovery from genome-wide association studies (GWASs) in the bone field in the last 12 years. To date, GWASs and meta-analyses have discovered hundreds of loci that are associated with bone mineral density (BMD), osteoporosis, and osteoporotic fractures. However, the GWAS approach has sometimes been criticized because of the small effect size of the discovered variants and the mystery of missing heritability, these two questions could be partially explained by the newly raised conceptual models, such as omnigenic model and natural selection. Finally, we introduce the clinical use of GWAS findings in the bone field, such as the identification of causal clinical risk factors, the development of drug targets and disease prediction. Despite the fruitful GWAS discoveries in the bone field, most of these GWAS participants were of European descent, and more genetic studies should be carried out in other ethnic populations to benefit disease prediction in the corresponding population.

Transcriptome prediction performance across machine learning models and diverse ancestries

Transcriptome prediction methods such as PrediXcan and FUSION have become popular in complex trait mapping. Most transcriptome prediction models have been trained in European populations using methods that make parametric linear assumptions like the elastic net (EN). To potentially further optimize imputation performance of gene expression across global populations, we built transcriptome prediction models using both linear and non-linear machine learning (ML) algorithms and evaluated their performance in comparison to EN. We trained models using genotype and blood monocyte transcriptome data from the Multi-Ethnic Study of Atherosclerosis (MESA) comprising individuals of African, Hispanic, and European ancestries and tested them using genotype and whole-blood transcriptome data from the Modeling the Epidemiology Transition Study (METS) comprising individuals of African ancestries. We show that the prediction performance is highest when the training and the testing population share similar ancestries regardless of the prediction algorithm used. While EN generally outperformed random forest (RF), support vector regression (SVR), and K nearest neighbor (KNN), we found that RF outperformed EN for some genes, particularly between disparate ancestries, suggesting potential robustness and reduced variability of RF imputation performance across global populations. When applied to a high-density lipoprotein (HDL) phenotype, we show including RF prediction models in PrediXcan revealed potential gene associations missed by EN models. Therefore, by integrating other ML modeling into PrediXcan and diversifying our training populations to include more global ancestries, we may uncover new genes associated with complex traits.

The Polymorphism at PLCB4 Promoter (rs6086746) Changes the Binding Affinity of RUNX2 and Affects Osteoporosis Susceptibility: An Analysis of Bioinformatics-Based Case-Control Study and Functional Validation

PURPOSE

Genome-wide association studies have identified numerous genetic variants that are associated with osteoporosis risk; however, most of them are present in the non-coding regions of the genome and the functional mechanisms are unknown. In this study, we aimed to investigate the potential variation in runt domain transcription factor 2 (RUNX2), which is an osteoblast-specific transcription factor that normally stimulates bone formation and osteoblast differentiation, regarding variants within RUNX2 binding sites and risk of osteoporosis in postmenopausal osteoporosis (PMOP).

METHODS

We performed bioinformatics-based prediction by combining whole genome sequencing and chromatin immunoprecipitation sequencing to screen functional SNPs in the RUNX2 binding site using data from the database of Taiwan Biobank; Case-control studies with 651 postmenopausal women comprising 107 osteoporosis patients, 290 osteopenia patients, and 254 controls at Tri-Service General Hospital between 2015 and 2019 were included. The subjects were examined for bone mass density and classified into normal and those with osteopenia or osteoporosis by T-scoring with dual-energy X-ray absorptiometry. Furthermore, mRNA expression and luciferase reporter assay were used to provide additional evidence regarding the associations identified in the association analyses. Chi-square tests and logistic regression were mainly used for statistical assessment.

RESULTS

Through candidate gene approaches, 3 SNPs in the RUNX2 binding site were selected. A novel SNP rs6086746 in the PLCB4 promoter was identified to be associated with osteoporosis in Chinese populations. Patients with AA allele had higher risk of osteoporosis than those with GG+AG (adjusted OR = 6.89; 95% confidence intervals: 2.23-21.31, p = 0.001). Moreover, the AA genotype exhibited lower bone mass density (p < 0.05). Regarding mRNA expression, there were large differences in the correlation between PLCB4 and different RUNX2 alleles (Cohen's q = 0.91). Functionally, the rs6086746 A allele reduces the RUNX2 binding affinity, thus enhancing the suppression of PLCB4 expression (p < 0.05).

CONCLUSIONS

Our results provide further evidence to support the important role of the SNP rs6086746 in the etiology of osteopenia/osteoporosis, thereby enhancing the current understanding of the susceptibility to osteoporosis. We further studied the mechanism underlying osteoporosis regulation by PLCB4.

Genetics of Progressive Supranuclear Palsy: A Review

Progressive supranuclear palsy (PSP) is an atypical parkinsonism with prominent 4R-tau neuropathology, and the classical clinical phenotype is characterized by vertical supranuclear gaze palsy, unprovoked falls, akinetic-rigid syndrome and cognitive decline. Though PSP is generally regarded as sporadic, there is increasing evidence suggesting that a series of common and rare genetic variants impact on sporadic and familial forms of PSP. To date, more than 10 genes have been reported to show a potential association with PSP. Among these genes, the microtubule-associated protein tau (MAPT) is the risk locus with the strongest effect size on sporadic PSP in the case-control genome-wide association studies (GWAS). Additionally, MAPT mutations are the most common cause of familial PSP while the leucine-rich repeat kinase 2 (LRRK2) is a rare monogenic cause of PSP, and several other gene mutations may mimic the PSP phenotype, like the dynactin subunit 1 (DCTN1). In total, 15 MAPT mutations have been identified in cases with PSP, and the mean age at onset is much earlier than in cases carrying LRRK2 or DCTN1 mutations. GWAS have further identified several risk loci of PSP, proposing molecular pathways related to PSP. The present review focused on genetic studies on PSP and summarized genetic factors of PSP, which may help to elucidate the underlying pathogenesis and provide new perspectives for therapeutic strategies.

USF3 modulates osteoporosis risk by targeting WNT16, RANKL, RUNX2, and two GWAS lead SNPs rs2908007 and rs4531631

Osteoporotic fractures cause major morbidity and mortality in the aging population. Genome-wide association studies (GWAS) have identified USF3 as the novel susceptibility gene of osteoporosis. However, the functional role in bone metabolism and the target gene of the basic helix-loop-helix transcription factor USF3 are unclear. Here, we show that USF3 enhances osteoblast differentiation and suppresses osteoclastogenesis in cultured human osteoblast-like U-2OS cells. Mechanistic studies revealed that transcription factor USF3 antagonistically interacts with anti-osteogenic TWIST1/TCF12 heterodimer in the WNT16 and RUNX2 promoter, and counteracts CREB1 and JUN/FOS in the RANKL promoter. Importantly, the osteoporosis GWAS variant g.1744A>G (rs2908007A>G) located in the WNT16 promoter confers G-allele-specific transcriptional modulation by USF3, TWIST1/TCF12 and TBX5/TBX15, and USF3 transactivates the osteoclastogenesis suppressor WNT16 promoter activity and antagonizes the repression of WNT16 by TWIST1 and TCF12. The risk G allele of osteoporosis GWAS variant g.3260A>G (rs4531631A>G) in the RANKL promoter facilitates the binding of CREB1 and JUN/FOS and enhances transactivation of the osteoclastogenesis contributor RANKL that is inhibited by USF3. Our findings uncovered the functional mechanisms of osteoporosis novel GWAS-associated gene USF3 and lead single nucleotide polymorphisms rs2908007 and rs4531631 in the regulation of bone formation and resorption.

METTL3 polymorphisms and Wilms tumor susceptibility in Chinese children: A five-center case-control study

BACKGROUND

Wilms tumor is a common pediatric tumor worldwide. Methyltransferase like 3 (METTL3) is a core gene of the N⁶ -methyladenosine (m⁶ A) modification that widely affects the transcription of tumor-related genes in eukaryotes. METTL3 has been extensively investigated in various tumors but not Wilms tumor.

METHODS

We describe a five-center case-control study with 414 patients and 1199 controls aiming to explore the associations between METTL3 polymorphisms (rs1061026 T>G, rs1061027 C>A, rs1139130 A>G and rs1263801 G>C) and Wilms tumor susceptibility. A TaqMan real-time polymerase chain reaction was performed for genotyping. Odds ratios (ORs) and 95% confidence intervals (CIs) were reported as evaluation indicators to determine any associations.

RESULTS

Referring to the preliminary analysis results, protective genotypes were identified as rs1061026 TG/GG, rs1061027 CA/AA, rs1139130 GG and rs1263801 GC/CC. The children with three protective genotypes were less likely to develop Wilms tumor than children without protective genotypes (adjusted OR = 0.68, 95% CI = 0.46-0.999, p = 0.0496). Similarly, stratified analysis of the subgroup aged > 18 months, carrying 3 or 4 protective genotypes, was a protective factor for Wilms tumor compared to carrying 0-2 protective genotypes (adjusted OR = 0.59 95% CI = 0.39-0.91, p = 0.016). However, we did not observe any other significant results.

CONCLUSIONS

The combined effect of METTL3 polymorphisms reduce Wilms tumor susceptibility in Chinese children. This conclusion requires further verification.

The open targets post-GWAS analysis pipeline

Motivation

Genome-wide association studies (GWAS) are a powerful method to detect even weak associations between variants and phenotypes; however, many of the identified associated variants are in non-coding regions, and presumably influence gene expression regulation. Identifying potential drug targets, i.e. causal protein-coding genes, therefore, requires crossing the genetics results with functional data.

Results

We present a novel data integration pipeline that analyses GWAS results in the light of experimental epigenetic and cis-regulatory datasets, such as ChIP-Seq, Promoter-Capture Hi-C or eQTL, and presents them in a single report, which can be used for inferring likely causal genes. This pipeline was then fed into an interactive data resource.

Availability and implementation

The analysis code is available at www.github.com/Ensembl/postgap and the interactive data browser at postgwas.opentargets.io.

Functional polymorphisms and transcriptional analysis in the 5' region of the human serotonin receptor 1B gene (HTR1B) and their associations with psychiatric disorders

BACKGROUND

The 5-hydroxytryptamine 1B receptor (5-HT1B) plays an essential role in the serotonin (5-HT) system and is widely involved in a variety of brain activities. HTR1B is the gene encoding 5-HT1B. Genome-wide association studies have shown that HTR1B polymorphisms are closely related to multiple mental and behavioral disorders; however, the functional mechanisms underlying these associations are unknown. This study investigated the effect of several HTR1B haplotypes on regulation of gene expression in vitro and the functional sequences in the 5' regulatory region of HTR1B to determine their potential association with mental and behavioral disorders.

METHODS

Six haplotypes consisting of rs4140535, rs1778258, rs17273700, rs1228814, rs11568817, and rs130058 and several truncated fragments of the 5' regulatory region of HTR1B were transfected into SK-N-SH and HEK-293 cells. The relative fluorescence intensities of the different haplotypes and truncated fragments were detected using a dual-luciferase reporter assay system.

RESULTS

Compared to the major haplotype T-G-T-C-T-A, the relative fluorescence intensities of haplotypes C-A-T-C-T-A, C-G-T-C-T-A, C-G-C-A-G-T, and C-G-T-A-T-A were significantly lower, and that of haplotype C-G-C-A-G-A was significantly higher. Furthermore, the effects of the rs4140535T allele, the rs17273700C-rs11568817G linkage combination, and the rs1228814A allele made their relative fluorescence intensities significantly higher than their counterparts at each locus. Conversely, the rs1778258A and rs130058T alleles decreased the relative fluorescence intensities. In addition, we found that regions from - 1587 to - 1371 bp (TSS, + 1), - 1149 to - 894 bp, - 39 to + 130 bp, + 130 to + 341 bp, and + 341 to + 505 bp upregulated gene expression. In contrast, regions - 603 to - 316 bp and + 130 to + 341 bp downregulated gene expression. Region + 341 to + 505 bp played a decisive role in gene transcription.

CONCLUSIONS

HTR1B 5' regulatory region polymorphisms have regulatory effects on gene expression and potential correlate with several pathology and physiology conditions. This study suggests that a crucial sequence for transcription is located in region + 341 ~ + 505 bp. Regions - 1587 to - 1371 bp, - 1149 to - 894 bp, - 603 to - 316 bp, - 39 to + 130 bp, and + 130 to + 341 bp contain functional sequences that can promote or suppress the HTR1B gene expression.

Bayesian Genome-wide TWAS Method to Leverage both cis- and trans-eQTL Information through Summary Statistics

Transcriptome-wide association studies (TWASs) have been widely used to integrate gene expression and genetic data for studying complex traits. Due to the computational burden, existing TWAS methods do not assess distant trans-expression quantitative trait loci (eQTL) that are known to explain important expression variation for most genes. We propose a Bayesian genome-wide TWAS (BGW-TWAS) method that leverages both cis- and trans-eQTL information for a TWAS. Our BGW-TWAS method is based on Bayesian variable selection regression, which not only accounts for cis- and trans-eQTL of the target gene but also enables efficient computation by using summary statistics from standard eQTL analyses. Our simulation studies illustrated that BGW-TWASs achieved higher power compared to existing TWAS methods that do not assess trans-eQTL information. We further applied BWG-TWAS to individual-level GWAS data (N = ∼3.3K), which identified significant associations between the genetically regulated gene expression (GReX) of ZC3H12B and Alzheimer dementia (AD) (p value = 5.42 × 10-13), neurofibrillary tangle density (p value = 1.89 × 10-6), and global measure of AD pathology (p value = 9.59 × 10-7). These associations for ZC3H12B were completely driven by trans-eQTL. Additionally, the GReX of KCTD12 was found to be significantly associated with β-amyloid (p value = 3.44 × 10-8) which was driven by both cis- and trans-eQTL. Four of the top driven trans-eQTL of ZC3H12B are located within APOC1, a known major risk gene of AD and blood lipids. Additionally, by applying BGW-TWAS with summary-level GWAS data of AD (N = ∼54K), we identified 13 significant genes including known GWAS risk genes HLA-DRB1 and APOC1, as well as ZC3H12B.

The Role of ARID5B in Acute Lymphoblastic Leukemia and Beyond

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children with distinct characteristics among different subtypes. Although the etiology of ALL has not been fully unveiled, initiation of ALL has been demonstrated to partly depend on genetic factors. As indicated by several genome wide association studies (GWASs) and candidate gene analyses, ARID5B, a member of AT-rich interactive domain (ARID) protein family, is associated with the occurrence and prognosis of ALL. However, the mechanisms by which ARID5B genotype impact on the susceptibility and treatment outcome remain vague. In this review, we outline developments in the understanding of ARID5B in the susceptibility of ALL and its therapeutic perspectives, and summarize the underlying mechanisms based on the limited functional studies, hoping to illustrate the possible mechanisms of ARID5B impact and highlight the potential treatment regimens.

Investigation of the association of endothelial nitric oxide synthase (eNOS)-T786C gene polymorphism with the risk of male infertility in an Iranian population

Recent studies have demonstrated that making the sorts of oxygen reactive, such as nitric oxide, can cause oxidative lipid damage, protein damage, and damage to the DNA of cells. Sperm DNA damage effect on the reduction of sperm mobility and damage of acrosome membrane lead to the inability of sperm to fertilize the oocyte. Increasing expression of endothelial nitric oxide synthase (eNOS) gene is seen in various diseases such as cardiovascular and infertility diseases. This study aimed to assess the association between eNOS gene single nucleotide polymorphism/SNP (rs2070744, T786C) and risk of male infertility and the quality of sperm parameters in an Iranian population. In this case-control study, 100 infertile men were enrolled as a patient group. Control groups consisted of 100 fertile men. T786C genotypes were determined using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Results showed that T786C SNP, contained frequent genotype TC (p = 0.000; OR = 0.000; 95% CI = 0.000-0.015), TC + CC genotypes (p = 0.000; OR = 0.000; 95% CI = 0.000-0.015), and C allele (p = 0.000; OR = 0.00; 95% CI = 0.000-0.007), revealed a significant with male infertility. Based on the findings of this study suggested that although T786C SNP could not be applied as an appropriate genetic risk factor for male infertility, it probably may be considered a protective marker for other researchers. However, more comprehensive studies in different populations are required to confirm our data.

Atypical chemokine receptor ACKR2-V41A has decreased CCL2 binding, scavenging, and activation, supporting sustained inflammation and increased Alzheimer's disease risk

A recent genome-wide association study (GWAS) of 59 cerebrospinal fluid (CSF) proteins with a connection to Alzheimer's disease (AD) demonstrated an association between increased levels of chemokine ligand 2 (CCL2) with an atypical chemokine receptor chemokine-binding protein 2 variant V41A (ACKR2-V41A; rs2228467). High levels of CCL2 are associated with increased risk of AD development as well as other inflammatory diseases. In this study we characterized the biological function of the ACKR2-V41A receptor compared to the wild type allele by measuring its ligand binding affinity, CCL2 scavenging efficiency, and cell activation sensitivity. We transfected Chinese hamster ovary cells with plasmids carrying wild type ACKR2 (ACKR2-WT) or the mutant ACKR2-V41A receptor. Binding affinity assays showed that ACKR2-V41A has a lower binding affinity for CCL2 and CCL4 than ACKR2-WT. CCL2 scavenging results aligned with binding affinity assays, with ACKR2-V41A cells scavenging CCL2 with a lower efficiency than ACKR2-WT. Cell activation assays also showed that ACKR2-V41A cells had significantly lower receptor upregulation (β-Arrestin-dependent signaling pathway) upon stimulation compared to ACKR2-WT cells. These findings provide molecular and biological mechanistic insights into the GWAS association of ACKR2-V41A with increased levels of CCL2 in CSF and possibly other chemokine ligands. Increased CCL2 levels are associated with accelerated cognitive decline and increased risk of AD. Understanding how this atypical chemokine receptor allele increases serum markers of inflammation could lead to novel therapeutic solutions for AD.

Computational and functional analyses of T2D GWAS SNPs for transcription factor binding

Genome-wide association studies (GWASs) have successfully identified numerous non-coding genetic variants for type 2 diabetes (T2D), but the functional roles underlying these non-coding variants remain largely unknown. The effects of T2D GWAS lead SNPs on transcriptional factors binding motifs were firstly analyzed via JASPAR, followed by functional validations including dual-luciferase reporter assays, biotin-based DNA pull-down assays, real-time quantitative PCR, and western blotting. The results showed that GWAS SNP rs4430796 conferred T allele specific transcriptional enhancer activity via a PAX6 binding element, and upregulated the expression of HNF1B. GWAS SNP rs4607103 showed a bidirectional modulation of ADAMTS9-AS2 and ADAMTS9 by TCF7L2 in a T allele-specific manner. GWAS SNP rs849135 conferred C allele-specific bidirectional transcriptional enhancer activity via a CREB1 binding element. Our findings have uncovered the functional mechanisms of three T2D GWAS SNPs via affecting the binding of transcription factors, providing new insights into the genetics and molecular pathogenesis of T2D.

Osteoporosis genome-wide association study variant c.3781 C>A is regulated by a novel anti-osteogenic factor miR-345-5p

Upstream transcription factor family member 3 (USF3) c.3781C>A (rs1026364) in the 3'-untranslated region (3'-UTR) has been firmly associated with bone mineral density (BMD) in genome-wide association study (GWAS). However, the molecular mechanism by which it influences BMD and osteoporosis is unknown. Bioinformatics analyses suggested that the risk c.3781A allele creates a target site for hsa-miR-345-5p binding. Luciferase assay validated that the c.3781A allele displayed significantly lower luciferase activities than the c.3781C allele in the human osteoblast cell line hFOB1.19, osteosarcoma cell lines U-2OS and Saos-2, and embryonic kidney cell line 293T. Furthermore, hsa-miR-345-5p regulated USF3 expression on both messenger RNA and protein levels in hFOB1.19 and U937 cells with heterozygous A/C genotype. Transfection of hsa-miR-345-5p antagomiR in heterozygous hFOB1.19 cells significantly increased the expression of osteogenic marker genes RUNX2, OSTERIX, COL1A1, ALP, OPN, OCN, and alkaline phosphatase activity and matrix mineralization level. Importantly, we found that hsa-miR-345-5p also inhibits osteoblast maturation in cell lines U-2OS with hsa-miR-345-5p nonbinding C/C genotype by targeting RUNX3 and SMAD1. Our findings uncovered a novel pathogenetic mechanism of osteoporosis by GWAS variant c.3781C>A-mediated disruption of hsa-miR-345-5p binding at the 3'-UTR of USF3 and the functional role of hsa-miR-345-5p in osteogenic differentiation.

LINC01149 variant modulates MICA expression that facilitates hepatitis B virus spontaneous recovery but increases hepatocellular carcinoma risk

Interpreting disease-causing variants, especially in noncoding regions by genome-wide association studies (GWAS), has become one of the most challenging and demanding tasks. We hypothesized that functional lncRNAs variants in GWAS-identified loci might alter expression level of genes associated with persistent HBV infection and hepatocellular carcinoma (HCC). Integrated bioinformatics approaches were used to prioritize potentially functional variants and a two-stage case-control study (2473 HBV positive HCC patients, 2248 persistent HBV carriers and 2294 spontaneously recovered subjects) was performed to assess the roles of these variants. The rs2844512 G > C variant in LINC01149 was identified to facilitate HBV spontaneous recovery (OR = 0.84, 95% CI = 0.77-0.92) but increase the risk of HCC (OR = 1.21, 95% CI = 1.11-1.32) in combined samples. Subsequent biological assays indicated this variant created a binding site for miR-128-3p and upregulated MICA expression by serving as a miRNA sponge, which might recruit NK-cells to lyse infected cells, but release highly soluble MICA by shedding to induce NK-cells exhaustion and tumor immune evasion. These findings highlight a regulatory circuit between LINC01149 and MICA, mediating by miR-128-3p, and the important role of upregulated MICA in conferring susceptibility to persistent HBV infection and HCC.

TIGAR: An Improved Bayesian Tool for Transcriptomic Data Imputation Enhances Gene Mapping of Complex Traits

The transcriptome-wide association studies (TWASs) that test for association between the study trait and the imputed gene expression levels from cis-acting expression quantitative trait loci (cis-eQTL) genotypes have successfully enhanced the discovery of genetic risk loci for complex traits. By using the gene expression imputation models fitted from reference datasets that have both genetic and transcriptomic data, TWASs facilitate gene-based tests with GWAS data while accounting for the reference transcriptomic data. The existing TWAS tools like PrediXcan and FUSION use parametric imputation models that have limitations for modeling the complex genetic architecture of transcriptomic data. Therefore, to improve on this, we employ a nonparametric Bayesian method that was originally proposed for genetic prediction of complex traits, which assumes a data-driven nonparametric prior for cis-eQTL effect sizes. The nonparametric Bayesian method is flexible and general because it includes both of the parametric imputation models used by PrediXcan and FUSION as special cases. Our simulation studies showed that the nonparametric Bayesian model improved both imputation R2 for transcriptomic data and the TWAS power over PrediXcan when ≥1% cis-SNPs co-regulate gene expression and gene expression heritability ≤0.2. In real applications, the nonparametric Bayesian method fitted transcriptomic imputation models for 57.8% more genes over PrediXcan, thus improving the power of follow-up TWASs. We implement both parametric PrediXcan and nonparametric Bayesian methods in a convenient software tool "TIGAR" (Transcriptome-Integrated Genetic Association Resource), which imputes transcriptomic data and performs subsequent TWASs using individual-level or summary-level GWAS data.

Genetic Influences on Patient-Oriented Outcomes in Traumatic Brain Injury: A Living Systematic Review of Non-Apolipoprotein E Single-Nucleotide Polymorphisms

There is a growing literature on the impact of genetic variation on outcome in traumatic brain injury (TBI). Whereas a substantial proportion of these publications have focused on the apolipoprotein E (APOE) gene, several have explored the influence of other polymorphisms. We undertook a systematic review of the impact of single-nucleotide polymorphisms (SNPs) in non–apolipoprotein E (non-APOE) genes associated with patient outcomes in adult TBI). We searched EMBASE, MEDLINE, CINAHL, and gray literature from inception to the beginning of August 2017 for studies of genetic variance in relation to patient outcomes in adult TBI. Sixty-eight articles were deemed eligible for inclusion into the systematic review. The SNPs described were in the following categories: neurotransmitter (NT) in 23, cytokine in nine, brain-derived neurotrophic factor (BDNF) in 12, mitochondrial genes in three, and miscellaneous SNPs in 21. All studies were based on small patient cohorts and suffered from potential bias. A range of SNPs associated with genes coding for monoamine NTs, BDNF, cytokines, and mitochondrial proteins have been reported to be associated with variation in global, neuropsychiatric, and behavioral outcomes. An analysis of the tissue, cellular, and subcellular location of the genes that harbored the SNPs studied showed that they could be clustered into blood–brain barrier associated, neuroprotective/regulatory, and neuropsychiatric/degenerative groups. Several small studies report that various NT, cytokine, and BDNF-related SNPs are associated with variations in global outcome at 6–12 months post-TBI. The association of these SNPs with neuropsychiatric and behavioral outcomes is less clear. A definitive assessment of role and effect size of genetic variation in these genes on outcome remains uncertain, but could be clarified by an adequately powered genome-wide association study with appropriate recording of outcomes.

Chromatin Conformation Links Putative Enhancers in Intracranial Aneurysm-Associated Regions to Potential Candidate Genes

Background

We previously showed that intracranial aneurysm (IA)–associated single‐nucleotide polymorphisms are enriched in promoters and putative enhancers identified in the human circle of Willis, on which IAs develop, suggesting a role for promoters and enhancers in IAs. We further investigated the role of putative enhancers in the pathogenesis of IA by identifying their potential target genes and validating their regulatory activity.

Methods and Results

Using our previously published circle of Willis chromatin immunoprecipitation and sequencing data, we selected 34 putative enhancers in IA‐associated regions from genome‐wide association studies. We then used a chromatin conformation capture technique to prioritize target genes and found that 15 putative enhancers interact with the promoters of 6 target genes: SOX17,CDKN2B,MTAP,CNNM2,RPEL1, and GATA6. Subsequently, we assessed the activity of these putative enhancers in vivo in zebrafish embryos and confirmed activity for 8 putative enhancers. Last, we found that all 6 target genes are expressed in the circle of Willis, on the basis of RNA sequencing data and in situ hybridization. Furthermore, in situ hybridization showed that these genes are expressed in multiple cell types in the circle of Willis.

Conclusions

In 4 of 6 IA‐associated genome‐wide association study regions, we identified 8 putative enhancers that are active in vivo and interact with 6 nearby genes, suggesting that these genes are regulated by the identified putative enhancers. These genes, SOX17,CDKN2B,MTAP,CNNM2,RPEL1, and GATA6, are therefore potential candidate genes involved in IA pathogenesis and should be studied using animal models in the future.

A Positive Causal Influence of IL-18 Levels on the Risk of T2DM: A Mendelian Randomization Study

A large number of clinical studies have shown that interleukin-18 (IL-18) plasma levels are positively correlated with the pathogenesis and development of type 2 diabetes mellitus (T2DM), but it remains unclear whether IL-18 causes T2DM, primarily due to the influence of reverse causality and residual confounding factors. Genome-wide association studies have led to the discovery of numerous common variants associated with IL-18 and T2DM and opened unprecedented opportunities for investigating possible associations between genetic traits and diseases. In this study, we employed a two-sample Mendelian randomization (MR) method to analyze the causal relationships between IL-18 plasma levels and T2DM using IL18-related SNPs as genetic instrumental variables (IVs). We first selected eight SNPs that were significantly associated with IL-18 but independent of T2DM. We then used these SNPs as IVs to evaluate their effects on T2DM using the inverse-variance weighted (IVW) method. Finally, we conducted sensitivity analysis and MR-Egger regression analysis to evaluate the heterogeneity and pleiotropic effects of each variant. The results based on the IVW method demonstrate that high IL-18 plasma levels significantly increase the risk of T2DM, and no heterogeneity or pleiotropic effects appeared after the sensitivity and MR-Egger analyses.

The association of nucleotide-binding oligomerization domain 2 gene polymorphisms with the risk of asthma in the Chinese Han population

BACKGROUND

Genetic background is one of the important risk factors for development of asthma. The nucleotide-binding oligomerization domain 2 (NOD2) has been involved in the pathogenesis of asthma. The purpose of this study was to explore the relationship between NOD2 gene polymorphisms and asthma susceptibility in the Chinese Han population.

METHODS

Children with asthma (n = 309) and Healthy children (n = 163) were recruited from Yancheng Third People's Hospital, Yancheng, China, between January 2016 and December 2017. The NOD2 gene polymorphisms were measured by the Snapshot SNP genotyping assays. Genotyping was performed for 4 tag SNPs of NOD2. Serum IFN-β levels were measured by ELISA.

RESULTS

The serum IFN-β levels were significantly lower in Asthmatic children than those in the controls (p < 0.001). Low levels of IFN-β may be related to the susceptibility to severe asthma. The rs3135499 C allele was associated with a significantly increased risk of asthma as compared with the rs3135499 A allele.

CONCLUSION

The rs3135499 polymorphism of NOD2 gene and IFN-β may play a role in the pathogenesis of asthma.

A generalized model for combining dependent SNP-level summary statistics and its extensions to statistics of other levels

Here I propose a fundamentally new flexible model to reveal the association between a trait and a set of genetic variants in a genomic region/gene. This model was developed for the situation when original individual-level phenotype and genotype data are not available, but the researcher possesses the results of statistical analyses conducted on these data (namely, SNP-level summary Z score statistics and SNP-by-SNP correlations). The new model was analytically derived from the classical multiple linear regression model applied for the region-based association analysis of individual-level phenotype and genotype data by using the linear compression of data, where the SNP-by-SNP correlations are among the explanatory variables, and the summary Z score statistics are categorized as the response variables. I analytically show that the regional association analysis methods developed within the framework of the classical multiple linear regression model with additive effects of genetic variants can be reformulated in terms of the new model without the loss of information. The results obtained from the regional association analysis utilizing the classical model and those derived using the proposed model are identical when SNP-by-SNP correlations and SNP-level statistics are estimated from the same genetic data.

Combined study on the single nucleotide polymorphisms in the follicle-stimulating hormone receptor (Ser680Asn) and anti-Müllerian hormone receptor type II (-482A>G) as genetic markers in assisted reproduction

Background Infertile women may have underlying genetic abnormalities. There is, at present, a significant number of studies on the relation between the follicle stimulating hormone receptor (FSHR) or anti-Müllerian hormone type II receptor (AMHRII) polymorphisms and response to in-vitro fertilisation (IVF) treatment. However, it is not yet clear which genotype or combination of genotypes is favourable towards a better ovarian stimulation and pregnancy outcome. Materials and methods In this study we assessed the distribution of the genotypes of FSHR Ser680Asn and of AMHRII -482A>G gene polymorphisms in a group of 126 infertile women and a control group of 100 fertile women by using real-time polymerase chain reaction (RT-PCR). Results Statistical analysis showed that the frequency of the genotypes is similar in both control and IVF/ intracytoplasmic sperm injection (ICSI) groups. Further investigation of the frequency of the nine possible combinations of these polymorphisms in the groups revealed no correlation between infertility and combination of the polymorphisms. Women with one polymorphism have on average 5.5 units higher levels of AMH compared to women carrying no polymorphism. In women with no polymorphisms, for each unit of FSH increase, the average concentration of blood AMH is expected to be 72% lower. Conclusion The distribution of the FSHR Ser680Asn and of the AMHRII -482A>G gene polymorphisms, in the Greek population is similar in fertile and infertile women. The study showed that FSH and AMH correlated levels in certain cases could be used to estimate a patient's ovarian reserve.

Insights of Acute Lymphoblastic Leukemia with Development of Genomic Investigation

Treatment outcomes for acute lymphoblastic leukemia (ALL), especially pediatric ALL, have greatly improved due to the risk-adapted therapy. Combination of drug development, clinical practice, as well as basic genetic researches has brought the survival rate of ALL from less than 10% to more than 90% today, not only increasing the treatment efficacy but also limiting adverse drug reactions (ADRs). In this review, we summarized the landscape identification of ALL genetic alterations, which provided the opportunity to increase the survival rate and especially minimize the relapse risk of ALL, and highlighted the importance of the development of new technologies of genomic investigation for translational medicine.

Timing during translation matters: synonymous mutations in human pathologies influence protein folding and function

Ribosomes translate mRNAs with non-uniform speed. Translation velocity patterns are a conserved feature of mRNA and have evolved to fine-tune protein folding, expression and function. Synonymous single-nucleotide polymorphisms (sSNPs) that alter programmed translational speed affect expression and function of the encoded protein. Synergistic advances in next-generation sequencing have led to the identification of sSNPs associated with disease penetrance. Here, we draw on studies with disease-related proteins to enhance our understanding of mechanistic contributions of sSNPs to functional alterations of the encoded protein. We emphasize the importance of identification of sSNPs along with disease-causing mutations to understand genotype-phenotype relationships.

TrioMDR: Detecting SNP interactions in trio families with model-based multifactor dimensionality reduction

Single nucleotide polymorphism (SNP) interactions can explain the missing heritability of common complex diseases. Many interaction detection methods have been proposed in genome-wide association studies, and they can be divided into two types: population-based and family-based. Compared with population-based methods, family-based methods are robust vs. population stratification. Several family-based methods have been proposed, among which Multifactor Dimensionality Reduction (MDR)-based methods are popular and powerful. However, current MDR-based methods suffer from heavy computational burden. Furthermore, they do not allow for main effect adjustment. In this work we develop a two-stage model-based MDR approach (TrioMDR) to detect multi-locus interaction in trio families (i.e., two parents and one affected child). TrioMDR combines the MDR framework with logistic regression models to check interactions, so TrioMDR can adjust main effects. In addition, unlike consuming permutation procedures used in traditional MDR-based methods, TrioMDR utilizes a simple semi-parameter P-values correction procedure to control type I error rate, this procedure only uses a few permutations to achieve the significance of a multi-locus model and significantly speeds up TrioMDR. We performed extensive experiments on simulated data to compare the type I error and power of TrioMDR under different scenarios. The results demonstrate that TrioMDR is fast and more powerful in general than some recently proposed methods for interaction detection in trios. The R codes of TrioMDR are available at: https://github.com/TrioMDR/TrioMDR.

From lipid locus to drug target through human genomics

In the last decade, over 175 genetic loci have robustly been associated to levels of major circulating blood lipids. Most loci are specific to one or two lipids, whereas some (SUGP1, ZPR1, TRIB1, HERPUD1, and FADS1) are associated to all. While exposing the polygenic architecture of circulating lipids and the underpinnings of dyslipidaemia, these genome-wide association studies (GWAS) have provided further evidence of the critical role that lipids play in coronary heart disease (CHD) risk, as indicated by the 2.7-fold enrichment for macrophage gene expression in atherosclerotic plaques and the association of 25 loci (such as PCSK9, APOB, ABCG5-G8, KCNK5, LPL, HMGCR, NPC1L1, CETP, TRIB1, ABO, PMAIP1-MC4R, and LDLR) with CHD. These GWAS also confirmed known and commonly used therapeutic targets, including HMGCR (statins), PCSK9 (antibodies), and NPC1L1 (ezetimibe). As we head into the post-GWAS era, we offer suggestions for how to move forward beyond genetic risk loci, towards refining the biology behind the associations and identifying causal genes and therapeutic targets. Deep phenotyping through lipidomics and metabolomics will refine and increase the resolution to find causal and druggable targets, and studies aimed at demonstrating gene transcriptional and regulatory effects of lipid associated loci will further aid in identifying these targets. Thus, we argue the need for deeply phenotyped, large genetic association studies to reduce costs and failures and increase the efficiency of the drug discovery pipeline. We conjecture that in the next decade a paradigm shift will tip the balance towards a data-driven approach to therapeutic target development and the application of precision medicine where human genomics takes centre stage.

Computational analyses of obesity associated loci generated by genome-wide association studies

OBJECTIVES

Genome-wide association studies (GWASs) have discovered associations of numerous SNPs and genes with obesity. However, the underlying molecular mechanisms through which these SNPs and genes affect the predisposition to obesity remain not fully understood. Aims of our study are to comprehensively characterize obesity GWAS SNPs and genes through computational approaches.

METHODS

For obesity GWAS identified SNPs, functional annotation, effects on miRNAs binding and impact on protein phosphorylation were performed via RegulomeDB and 3DSNP, miRNASNP, and the PhosSNP 1.0 database, respectively. For obesity associated genes, protein-protein interaction network construction, gene ontology and pathway enrichment analyses were performed by STRING, PANTHER and STRING, respectively.

RESULTS

A total of 445 SNPs are significantly associated with obesity related phenotypes at threshold P < 5×10-8. A number of SNPs were eQTLs for obesity associated genes, some SNPs located at binding sites of obesity related transcription factors. SNPs that might affect miRNAs binding and protein phosphorylation were identified. Protein-protein interaction network analysis identified the highly-interconnected "hub" genes. Obesity associated genes mainly involved in metabolic process and catalytic activity, and significantly enriched in 15 signal pathways.

CONCLUSIONS

Our results provided the targets for follow-up experimental testing and further shed new light on obesity pathophysiology.

Computational and functional characterization of four SNPs in the SOST locus associated with osteoporosis

The SOST gene encodes sclerostin, a C-terminal cysteine knot-like domain containing key negative regulator of osteoblastic bone formation that inhibits LRP5/6-mediated canonical Wnt signaling. Numerous single nucleotide polymorphisms (SNPs) in the SOST locus are firmly associated with bone mineral density (BMD) and fracture in genome-wide association studies (GWAS) and candidate gene association studies. However, the validation and mechanistic elucidation of causal genetic variants, especially for SNPs located beyond the promoter-proximal region, remain largely unresolved. By employing computational and experimental approaches, here we identify four SNPs rs1230399, rs7220711, rs1107748 and rs75901553 as functional variants which display allelic variation in SOST gene expression. The osteoporosis associated SNP rs1230399 in the SOST distal upstream regulatory region shows FOXA1 binding activity with subsequent transinactivation in a T allele-specific manner. The BMD GWAS lead SNPs rs7220711 and rs1107748 both reside in the 52-kb regulatory element deletion 35-kb downstream of the SOST gene which leads to Van Buchem disease. The rs7220711-A has a higher affinity for the transcriptional repressors MAFF or MAFK homodimers than rs7220711-G, while rs1107748 confers C allele specific transcriptional enhancer activity via a CTCF binding element. The variant rs75901553 C>T located in a conserved site of the SOST 3' UTR abolishes a target binding site for miR-98-5p which is negatively responsive to parathyroid hormone or 17β-estradiol in osteoblastic cell lines. Our findings uncover the biological consequences of four independent genetic variants in the SOST region and their important roles in SOST expression via diverse mechanisms, providing new insights into the genetics and molecular pathogenesis of osteoporosis.

Association of Cx43 rs2071166 polymorphism with an increased risk for atrial septal defect

Atrial septal defect is one of the most common CHD. The pathogenesis of atrial septal defect still remains unknown. Cx43 is the most prevalent connexin in the mammalian heart during development. Its genetic variants can cause several CHD. The aim of our study was to investigate the association of genetic variations of the Cx43 with sporadic atrial septal defect. A total of 450 paediatric patients were recruited, including 150 cases with atrial septal defect and 300 healthy controls. The promoter region of Cx43 was analysed by sequencing after polymerase chain reaction. All data were analysed by using the Statistic Package for Social Science 19.0 software. The frequency of the single nucleotide polymorphism rs2071166 was significantly higher in atrial septal defect cases than in healthy controls. The CC genotype at rs2071166 site in Cx43 was correlated with an increased risk for atrial septal defect (p<0.0001, odds ratio=3.891, 95% confidence interval 1.948-7.772) and the C allele was positively correlated with atrial septal defect (p=0.007, odds ratio=1.567, 95% confidence interval 1.129-2.175). In conclusion, our results confirmed that rs2071166 in Cx43 may be relevant with an increased atrial septal defect risk.

COPD GWAS variant at 19q13.2 in relation with DNA methylation and gene expression

Chronic obstructive pulmonary disease (COPD) is among the major health burdens in adults. While cigarette smoking is the leading risk factor, a growing number of genetic variations have been discovered to influence disease susceptibility. Epigenetic modifications may mediate the response of the genome to smoking and regulate gene expression. Chromosome 19q13.2 region is associated with both smoking and COPD, yet its functional role is unclear. Our study aimed to determine whether rs7937 (RAB4B, EGLN2), a top genetic variant in 19q13.2 region identified in genome-wide association studies of COPD, is associated with differential DNA methylation in blood (N = 1490) and gene expression in blood (N = 721) and lungs (N = 1087). We combined genetic and epigenetic data from the Rotterdam Study (RS) to perform the epigenome-wide association analysis of rs7937. Further, we used genetic and transcriptomic data from blood (RS) and from lung tissue (Lung expression quantitative trait loci mapping study), to perform the transcriptome-wide association study of rs7937. Rs7937 was significantly (FDR < 0.05) and consistently associated with differential DNA methylation in blood at 4 CpG sites in cis, independent of smoking. One methylation site (cg11298343-EGLN2) was also associated with COPD (P = 0.001). Additionally, rs7937 was associated with gene expression levels in blood in cis (EGLN2), 42% mediated through cg11298343, and in lung tissue, in cis and trans (NUMBL, EGLN2, DNMT3A, LOC101929709 and PAK2). Our results suggest that changes of DNA methylation and gene expression may be intermediate steps between genetic variants and COPD, but further causal studies in lung tissue should confirm this hypothesis.

The role of histone modification and a regulatory single-nucleotide polymorphism (rs2071166) in the Cx43 promoter in patients with TOF

Abnormal level of Cx43 expression could result in CHD. Epigenetic modification and disease-associated, non-coding SNPs might influence gene transcription and expression. Our study aimed to determine the role of histone modification and an rSNP (rs2071166) in the Cx43 promoter in patients with TOF. Our results indicate that H3K18ac bind to Cx43 promoter and that their levels are reduced in TOF patients relative to controls. The relationship between the non-coding SNP in the Cx43 gene and TOF patients was evaluated in 158 patients and 300 controls. The C allele of rs2071166 was confirmed to result in an increased risk of TOF (OR = 1.586, 95%CI 1.149–2.189). Individuals with the CC genotype at rs2071166 also showed a significant susceptibility to TOF (OR = 2.961, 95%CI 1.452–6.038). The mRNA level in TOF who were CC genotype was lower than that in patients with the AA/AC genotype. Functional analysis in cells and transgenic zebrafish models showed that rs2071166 decreased the activity of the promoter and could block the interaction between RXRα and RARE. This is the first study to illustrate that epigenetic modification and an rSNP in the Cx43 promoter region play a critical role in TOF by impacting the transcriptional activity and expression level of Cx43.

Neuregulin signaling pathway in smoking behavior

Understanding molecular processes that link comorbid traits such as addictions and mental disorders can provide novel therapeutic targets. Neuregulin signaling pathway (NSP) has previously been implicated in schizophrenia, a neurodevelopmental disorder with high comorbidity to smoking. Using a Finnish twin family sample, we have previously detected association between nicotine dependence and ERBB4 (a neuregulin receptor), and linkage for smoking initiation at the ERBB4 locus on 2q33. Further, Neuregulin3 has recently been shown to associate with nicotine withdrawal in a behavioral mouse model. In this study, we scrutinized association and linkage between 15 036 common, low frequency and rare genetic variants in 10 NSP genes and phenotypes encompassing smoking and alcohol use. Using the Finnish twin family sample (N=1998 from 740 families), we detected 66 variants (representing 23 LD blocks) significantly associated (false discovery rate P<0.05) with smoking initiation, nicotine dependence and nicotine withdrawal. We comprehensively annotated the associated variants using expression (eQTL) and methylation quantitative trait loci (meQTL) analyses in a Finnish population sample. Among the 66 variants, we identified 25 eQTLs (in NRG1 and ERBB4), 22 meQTLs (in NRG3, ERBB4 and PSENEN), a missense variant in NRG1 (rs113317778) and a splicing disruption variant in ERBB4 (rs13385826). Majority of the QTLs in blood were replicated in silico using publicly available databases, with additional QTLs observed in brain. In conclusion, our results support the involvement of NSP in smoking behavior but not in alcohol use and abuse, and disclose functional potential for 56 of the 66 associated single-nucleotide polymorphism.

Diversity and inclusion in genomic research: why the uneven progress?

Conducting genomic research in diverse populations has led to numerous advances in our understanding of human history, biology, and health disparities, in addition to discoveries of vital clinical significance. Conducting genomic research in diverse populations is also important in ensuring that the genomic revolution does not exacerbate health disparities by facilitating discoveries that will disproportionately benefit well-represented populations. Despite the general agreement on the need for genomic research in diverse populations in terms of equity and scientific progress, genomic research remains largely focused on populations of European descent. In this article, we describe the rationale for conducting genomic research in diverse populations by reviewing examples of advances facilitated by their inclusion. We also explore some of the factors that perpetuate the disproportionate attention on well-represented populations. Finally, we discuss ongoing efforts to ameliorate this continuing bias. Collaborative and intensive efforts at all levels of research, from the funding of studies to the publication of their findings, will be necessary to ensure that genomic research does not conserve historical inequalities or curtail the contribution that genomics could make to the health of all humanity.

Association of regulatory TPH2 polymorphisms with higher reduction in depressive symptoms in children and adolescents treated with fluoxetine

Genetic variability related to the brain serotonergic system has a significant impact on both the susceptibility to psychiatric disorders, such as major depressive disorder (MDD), and the response to antidepressant drugs, such as fluoxetine. TPH2 is one of the most important serotonergic candidate genes in selective serotonin reuptake inhibitors (SSRIs) pharmacogenetic studies. The aim of the present study was to evaluate the influence of regulatory polymorphisms that are specifically located in human TPH2 transcription factor binding sites (TFBSs), and therefore could be functional by altering gene expression, on clinical improvement in children and adolescents treated with fluoxetine. The selection of SNPs was also based on their linkage disequilibrium with TPH2 rs4570625, a genetic variant with questionable functionality, which was previously associated with clinical response in our pediatric population. A total of 83 children and adolescents were clinically evaluated 12weeks after initiating antidepressant treatment with fluoxetine for the first time. Clinical improvement was assessed by reductions in depressive symptoms measured using the Children's Depression Inventory (CDI) scale. The polymorphisms rs11179002, rs60032326 and rs34517220 were, for the first time in the literature, significantly associated with higher clinical improvement. The strongest association was found for rs34517220. In particular, minor allele homozygotes showed higher score reductions on the CDI scale compared with the major allele carriers. Interestingly, this polymorphism is located in a human TPH2 TFBS for two relevant transcription factors in the serotoninergic neurons, Foxa1 and Foxa2, which together with the high level of significance found for this SNP, could indicate that rs34517220 is in fact the crucial functional genetic variant related to the fluoxetine response. These results provide new evidence for the role of regulatory genetic variants that could modulate human TPH2 expression in the SSRI antidepressant response.

Functional annotation of Alzheimer's disease associated loci revealed by GWASs

Genome-wide association studies (GWASs) discovered a number of SNPs and genes associated with Alzheimer's disease (AD). However, how these SNPs and genes influence the liability to AD is not fully understood. We deployed computational approaches to explore the function and action mechanisms of AD -related SNPs and genes identified by GWASs, including the effects of 195 GWAS lead SNPs and 338 proxy SNPs on miRNAs binding and protein phosphorylation, their RegulomeDB and 3DSNP scores, and gene ontology, pathway enrichment and protein-protein interaction network of 126 AD-associated genes. Our computational analysis identified 6 lead SNPs (rs10119, rs1048699, rs148763909, rs610932, rs6857 and rs714948) and 2 proxy SNPs (rs12539172 and rs2847655) that potentially impacted the miRNA binding. Lead SNP rs2296160 and proxy SNPs rs679620 and rs2228145 were identified as PhosSNPs potentially influencing protein phosphorylation. AD-associated genes showed enrichment of “regulation of beta-amyloid formation”, “regulation of neurofibrillary tangle assembly”, “leukocyte mediated immunity” and “protein-lipid complex assembly” signaling pathway. Protein-protein interaction network and functional module analyses identified highly-interconnected “hub” genes (APOE, PICALM, BIN1, ABCA7, CD2AP, CLU, CR1, MS4A4E and MS4A6A) and bottleneck genes (APOE, TOMM40, NME8, PICALM, CD2AP, ZCWPW1, FAM180B, GAB2 and PTK2B) that created three tight subnetworks. Our results provided the targets for further experimental assessment and further insight on AD pathophysiology.

Whole-genome sequencing of monozygotic twins discordant for schizophrenia indicates multiple genetic risk factors for schizophrenia

Schizophrenia is a common disorder with a high heritability, but its genetic architecture is still elusive. We implemented whole-genome sequencing (WGS) analysis of 8 families with monozygotic (MZ) twin pairs discordant for schizophrenia to assess potential association of de novo mutations (DNMs) or inherited variants with susceptibility to schizophrenia. Eight non-synonymous DNMs (including one splicing site) were identified and shared by twins, which were either located in previously reported schizophrenia risk genes (p.V24689I mutation in TTN, p.S2506T mutation in GCN1L1, IVS3+1G > T in DOCK1) or had a benign to damaging effect according to in silico prediction analysis. By searching the inherited rare damaging or loss-of-function (LOF) variants and common susceptible alleles from three classes of schizophrenia candidate genes, we were able to distill genetic alterations in several schizophrenia risk genes, including GAD1, PLXNA2, RELN and FEZ1. Four inherited copy number variations (CNVs; including a large deletion at 16p13.11) implicated for schizophrenia were identified in four families, respectively. Most of families carried both missense DNMs and inherited risk variants, which might suggest that DNMs, inherited rare damaging variants and common risk alleles together conferred to schizophrenia susceptibility. Our results support that schizophrenia is caused by a combination of multiple genetic factors, with each DNM/variant showing a relatively small effect size.