SCAN: SNP and copy number annotation

sc2GWAS: a comprehensive platform linking single cell and GWAS traits of human

Identifying cell populations associated with risk variants is essential for uncovering cell-specific mechanisms that drive disease development and progression. Integrating genome-wide association studies (GWAS) with single-cell RNA sequencing (scRNA-seq) has become an effective strategy for detecting trait-cell relationships. The accumulation of trait-related single cell data has led to an urgent need for its comprehensively processing. To address this, we developed sc2GWAS (https://bio.liclab.net/sc2GWAS/), which aims to document large-scale GWAS trait-cell regulatory pairs at single-cell resolution and provide comprehensive annotations and enrichment analyses for these related pairs. The current version of sc2GWAS curates a total of 15 078 310 candidate trait-cell pairs from > 6 300 000 individual cells, offering a valuable resource for exploring complex regulatory relationships between traits and cells. We applied strict quality control measures on both scRNA-seq data and GWAS data, ensuring the reliability and accuracy of the datasets for the identification of trait-relevant cells and genes. In addition, sc2GWAS provides ranked lists of trait-relevant genes and extensive (epi) genetic annotations, making it a valuable resource for downstream analyses. We demonstrate the utility of the platform by investigating Alzheimer's disease, where we identified significant associations between the disease and microglial cells, with the APOE gene emerging as particularly significant. This platform facilitates detailed research into complex trait-cell and trait-gene interactions, we anticipate that sc2GWAS will become a comprehensive and valuable platform for exploring GWAS trait-cell regulatory mechanisms.

Genetic Risk of Ankylosing Spondylitis and Second-Line Therapy Need in Crohn's Disease: A Mendelian Randomization Study

Background: Crohn's disease (CD) and Ankylosing Spondylitis (AS) are chronic conditions with overlapping inflammatory pathways. This research investigates the genetic association between AS and the requirement for more aggressive therapeutic interventions in CD, suggesting a likelihood of increased severity in CD progression among individuals diagnosed with AS. Methods: This study utilized two-sample Mendelian randomization (TSMR) to analyze GWAS datasets for AS and CD requiring second-line treatment. Instrumental variables were selected based on single-nucleotide polymorphisms of genome-wide significance. Analytical methods included inverse-variance weighted (IVW), MR Egger, and other MR approaches, alongside sensitivity analysis, to validate the findings. Results: Our results indicated a significant association between AS genetic predisposition and the increased need for second-line treatments in CD. The IVW method showed an Odds Ratio (OR) of 2.16, and MR Egger provided an OR of 2.71, both were statistically significant. This association persisted even after the exclusion of influential outlier SNP rs2517655, confirming the robustness of our findings. Conclusions: This study suggests that genetic factors contributing to AS may influence the progression of CD, potentially necessitating more intensive treatment strategies. These findings underscore the importance of early screening in patients with co-existing AS and CD for tailoring treatment approaches, thus advancing personalized medicine in the management of these complex conditions.

Deciphering the genetics and mechanisms of predisposition to multiple myeloma

Multiple myeloma (MM) is an incurable malignancy of plasma cells. Epidemiological studies indicate a substantial heritable component, but the underlying mechanisms remain unclear. Here, in a genome-wide association study totaling 10,906 cases and 366,221 controls, we identify 35 MM risk loci, 12 of which are novel. Through functional fine-mapping and Mendelian randomization, we uncover two causal mechanisms for inherited MM risk: longer telomeres; and elevated levels of B-cell maturation antigen (BCMA) and interleukin-5 receptor alpha (IL5RA) in plasma. The largest increase in BCMA and IL5RA levels is mediated by the risk variant rs34562254-A at TNFRSF13B. While individuals with loss-of-function variants in TNFRSF13B develop B-cell immunodeficiency, rs34562254-A exerts a gain-of-function effect, increasing MM risk through amplified B-cell responses. Our results represent an analysis of genetic MM predisposition, highlighting causal mechanisms contributing to MM development.

KnockTF 2.0: a comprehensive gene expression profile database with knockdown/knockout of transcription (co-)factors in multiple species

Transcription factors (TFs), transcription co-factors (TcoFs) and their target genes perform essential functions in diseases and biological processes. KnockTF 2.0 (http://www.licpathway.net/KnockTF/index.html) aims to provide comprehensive gene expression profile datasets before/after T(co)F knockdown/knockout across multiple tissue/cell types of different species. Compared with KnockTF 1.0, KnockTF 2.0 has the following improvements: (i) Newly added T(co)F knockdown/knockout datasets in mice, Arabidopsis thaliana and Zea mays and also an expanded scale of datasets in humans. Currently, KnockTF 2.0 stores 1468 manually curated RNA-seq and microarray datasets associated with 612 TFs and 172 TcoFs disrupted by different knockdown/knockout techniques, which are 2.5 times larger than those of KnockTF 1.0. (ii) Newly added (epi)genetic annotations for T(co)F target genes in humans and mice, such as super-enhancers, common SNPs, methylation sites and chromatin interactions. (iii) Newly embedded and updated search and analysis tools, including T(co)F Enrichment (GSEA), Pathway Downstream Analysis and Search by Target Gene (BLAST). KnockTF 2.0 is a comprehensive update of KnockTF 1.0, which provides more T(co)F knockdown/knockout datasets and (epi)genetic annotations across multiple species than KnockTF 1.0. KnockTF 2.0 facilitates not only the identification of functional T(co)Fs and target genes but also the investigation of their roles in the physiological and pathological processes.

eRNAbase: a comprehensive database for decoding the regulatory eRNAs in human and mouse

Enhancer RNAs (eRNAs) transcribed from distal active enhancers serve as key regulators in gene transcriptional regulation. The accumulation of eRNAs from multiple sequencing assays has led to an urgent need to comprehensively collect and process these data to illustrate the regulatory landscape of eRNAs. To address this need, we developed the eRNAbase (http://bio.liclab.net/eRNAbase/index.php) to store the massive available resources of human and mouse eRNAs and provide comprehensive annotation and analyses for eRNAs. The current version of eRNAbase cataloged 10 399 928 eRNAs from 1012 samples, including 858 human samples and 154 mouse samples. These eRNAs were first identified and uniformly processed from 14 eRNA-related experiment types manually collected from GEO/SRA and ENCODE. Importantly, the eRNAbase provides detailed and abundant (epi)genetic annotations in eRNA regions, such as super enhancers, enhancers, common single nucleotide polymorphisms, expression quantitative trait loci, transcription factor binding sites, CRISPR/Cas9 target sites, DNase I hypersensitivity sites, chromatin accessibility regions, methylation sites, chromatin interactions regions, topologically associating domains and RNA spatial interactions. Furthermore, the eRNAbase provides users with three novel analyses including eRNA-mediated pathway regulatory analysis, eRNA-based variation interpretation analysis and eRNA-mediated TF-target gene analysis. Hence, eRNAbase is a powerful platform to query, browse and visualize regulatory cues associated with eRNAs.

CRdb: a comprehensive resource for deciphering chromatin regulators in human

Chromatin regulators (CRs) regulate epigenetic patterns on a partial or global scale, playing a critical role in affecting multi-target gene expression. As chromatin immunoprecipitation sequencing (ChIP-seq) data associated with CRs are rapidly accumulating, a comprehensive resource of CRs needs to be built urgently for collecting, integrating, and processing these data, which can provide abundant annotated information on CR upstream and downstream regulatory analyses as well as CR-related analysis functions. This study established an integrative CR resource, named CRdb (http://cr.liclab.net/crdb/), with the aim of curating a large number of available resources for CRs and providing extensive annotations and analyses of CRs to help biological researchers clarify the regulation mechanism and function of CRs. The CRdb database comprised a total of 647 CRs and 2,591 ChIP-seq samples from more than 300 human tissues and cell types. These samples have been manually curated from NCBI GEO/SRA and ENCODE. Importantly, CRdb provided the abundant and detailed genetic annotations in CR-binding regions based on ChIP-seq. Furthermore, CRdb supported various functional annotations and upstream regulatory information on CRs. In particular, it embedded four types of CR regulatory analyses: CR gene set enrichment, CR-binding genomic region annotation, CR-TF co-occupancy analysis, and CR regulatory axis analysis. CRdb is a useful and powerful resource that can help in exploring the potential functions of CRs and their regulatory mechanism in diseases and biological processes.

Genome-wide association study for systemic lupus erythematosus in an egyptian population

Systemic lupus erythematosus (SLE) susceptibility has a strong genetic component. Genome-wide association studies (GWAS) across trans-ancestral populations show both common and distinct genetic variants of susceptibility across European and Asian ancestries, while many other ethnic populations remain underexplored. We conducted the first SLE GWAS on Egyptians-an admixed North African/Middle Eastern population-using 537 patients and 883 controls. To identify novel susceptibility loci and replicate previously known loci, we performed imputation-based association analysis with 6,382,276 SNPs while accounting for individual admixture. We validated the association analysis using adaptive permutation tests (n = 109). We identified a novel genome-wide significant locus near IRS1/miR-5702 (Pcorrected = 1.98 × 10-8) and eight novel suggestive loci (Pcorrected < 1.0 × 10-5). We also replicated (Pperm < 0.01) 97 previously known loci with at least one associated nearby SNP, with ITGAM, DEF6-PPARD and IRF5 the top three replicated loci. SNPs correlated (r 2 > 0.8) with lead SNPs from four suggestive loci (ARMC9, DIAPH3, IFLDT1, and ENTPD3) were associated with differential gene expression (3.5 × 10-95 < p < 1.0 × 10-2) across diverse tissues. These loci are involved in cellular proliferation and invasion-pathways prominent in lupus and nephritis. Our study highlights the utility of GWAS in an admixed Egyptian population for delineating new genetic associations and for understanding SLE pathogenesis.

TcoFBase: a comprehensive database for decoding the regulatory transcription co-factors in human and mouse

Transcription co-factors (TcoFs) play crucial roles in gene expression regulation by communicating regulatory cues from enhancers to promoters. With the rapid accumulation of TcoF associated chromatin immunoprecipitation sequencing (ChIP-seq) data, the comprehensive collection and integrative analyses of these data are urgently required. Here, we developed the TcoFBase database (http://tcof.liclab.net/TcoFbase), which aimed to document a large number of available resources for mammalian TcoFs and provided annotations and enrichment analyses of TcoFs. TcoFBase curated 2322 TcoFs and 6759 TcoFs associated ChIP-seq data from over 500 tissues/cell types in human and mouse. Importantly, TcoFBase provided detailed and abundant (epi) genetic annotations of ChIP-seq based TcoF binding regions. Furthermore, TcoFBase supported regulatory annotation information and various functional annotations for TcoFs. Meanwhile, TcoFBase embedded five types of TcoF regulatory analyses for users, including TcoF gene set enrichment, TcoF binding genomic region annotation, TcoF regulatory network analysis, TcoF-TF co-occupancy analysis and TcoF regulatory axis analysis. TcoFBase was designed to be a useful resource that will help reveal the potential biological effects of TcoFs and elucidate TcoF-related regulatory mechanisms.

From GWAS to Gene: Transcriptome-Wide Association Studies and Other Methods to Functionally Understand GWAS Discoveries

Since their inception, genome-wide association studies (GWAS) have identified more than a hundred thousand single nucleotide polymorphism (SNP) loci that are associated with various complex human diseases or traits. The majority of GWAS discoveries are located in non-coding regions of the human genome and have unknown functions. The valley between non-coding GWAS discoveries and downstream affected genes hinders the investigation of complex disease mechanism and the utilization of human genetics for the improvement of clinical care. Meanwhile, advances in high-throughput sequencing technologies reveal important genomic regulatory roles that non-coding regions play in the transcriptional activities of genes. In this review, we focus on data integrative bioinformatics methods that combine GWAS with functional genomics knowledge to identify genetically regulated genes. We categorize and describe two types of data integrative methods. First, we describe fine-mapping methods. Fine-mapping is an exploratory approach that calibrates likely causal variants underneath GWAS signals. Fine-mapping methods connect GWAS signals to potentially causal genes through statistical methods and/or functional annotations. Second, we discuss gene-prioritization methods. These are hypothesis generating approaches that evaluate whether genetic variants regulate genes via certain genetic regulatory mechanisms to influence complex traits, including colocalization, mendelian randomization, and the transcriptome-wide association study (TWAS). TWAS is a gene-based association approach that investigates associations between genetically regulated gene expression and complex diseases or traits. TWAS has gained popularity over the years due to its ability to reduce multiple testing burden in comparison to other variant-based analytic approaches. Multiple types of TWAS methods have been developed with varied methodological designs and biological hypotheses over the past 5 years. We dive into discussions of how TWAS methods differ in many aspects and the challenges that different TWAS methods face. Overall, TWAS is a powerful tool for identifying complex trait-associated genes. With the advent of single-cell sequencing, chromosome conformation capture, gene editing technologies, and multiplexing reporter assays, we are expecting a more comprehensive understanding of genomic regulation and genetically regulated genes underlying complex human diseases and traits in the future.

Variation in genes implicated in B-cell development and antibody production affects susceptibility to pemphigus

Pemphigus foliaceus (PF) is an autoimmune blistering skin disease characterized by the presence of pathogenic autoantibodies against desmoglein 1, a component of intercellular desmosome junctions. PF occurs sporadically across the globe and is endemic in some Brazilian regions. Because PF is a B-cell-mediated disease, we aimed to study the impact of variants within genes encoding molecules involved in the different steps of B-cell development and antibody production on the susceptibility of endemic PF. We analysed 3,336 single nucleotide polymorphisms (SNPs) from 167 candidate genes genotyped with Illumina microarray in a cohort of 227 PF patients and 193 controls. After quality control and exclusion of non-informative and redundant SNPs, 607 variants in 149 genes remained in the logistic regression analysis, in which sex and ancestry were included as covariates. Our results revealed 10 SNPs within or nearby 11 genes that were associated with susceptibility to endemic PF (OR >1.56; p < 0.005): rs6657275*G (TGFB2); rs1818545*A (RAG1/RAG2/IFTAP);rs10781530*A (PAXX), rs10870140*G and rs10781522*A (TRAF2); rs535068*A (TNFRSF1B); rs324011*A (STAT6);rs6432018*C (YWHAQ); rs17149161*C (YWHAG); and rs2070729*C (IRF1). Interestingly, these SNPs have been previously associated with differential gene expression, mostly in peripheral blood, in publicly available databases. For the first time, we show that polymorphisms in genes involved in B-cell development and antibody production confer differential susceptibility to endemic PF, and therefore are candidates for possible functional studies to understand immunoglobulin gene rearrangement and its impact on diseases.

Qtlizer: comprehensive QTL annotation of GWAS results

Exploration of genetic variant-to-gene relationships by quantitative trait loci such as expression QTLs is a frequently used tool in genome-wide association studies. However, the wide range of public QTL databases and the lack of batch annotation features complicate a comprehensive annotation of GWAS results. In this work, we introduce the tool “Qtlizer” for annotating lists of variants in human with associated changes in gene expression and protein abundance using an integrated database of published QTLs. Features include incorporation of variants in linkage disequilibrium and reverse search by gene names. Analyzing the database for base pair distances between best significant eQTLs and their affected genes suggests that the commonly used cis-distance limit of 1,000,000 base pairs might be too restrictive, implicating a substantial amount of wrongly and yet undetected eQTLs. We also ranked genes with respect to the maximum number of tissue-specific eQTL studies in which a most significant eQTL signal was consistent. For the top 100 genes we observed the strongest enrichment with housekeeping genes (P = 2 × 10^–6) and with the 10% highest expressed genes (P = 0.005) after grouping eQTLs by r² > 0.95, underlining the relevance of LD information in eQTL analyses. Qtlizer can be accessed via https://genehopper.de/qtlizer or by using the respective Bioconductor R-package (https://doi.org/10.18129/B9.bioc.Qtlizer).

Pleiotropic loci underlying bone mineral density and bone size identified by a bivariate genome-wide association analysis

Aiming to identify pleiotropic genomic loci for bone mineral density and bone size, we performed a bivariate GWAS in five discovery samples and replicated in two large-scale samples. We identified 2 novel loci at 2q37.1 and 6q26. Our findings provide insight into common genetic architecture underlying both traits.

INTRODUCTION

Bone mineral density (BMD) and bone size (BS) are two important factors that contribute to the development of osteoporosis and osteoporotic fracture. Both BMD and BS are highly heritable and they are genetically correlated. In this study, we aim to identify pleiotropic loci associated with BMD and BS.

METHODS

We conducted a bivariate genome-wide association (GWA) analysis of hip BMD and hip BS in 6180 participants from 5 samples, followed by in silico replication in the UK Biobank study of BMD (N = 426,824) and the deCODE study of BS (N = 28,954), respectively.

RESULTS

SNPs from 2 genomic loci were significant at the genome-wide significance (GWS) level (p lt; 5 × 10-8) in the discovery samples and were successfully replicated in the replication samples (2q37.1, lead SNP rs7575512, discovery p = 1.49 × 10-10, replication p = 0.05; 6q26, lead SNP rs1040724, discovery p = 1.95 × 10-8, replication p = 0.03). Functional annotations suggested functional relevance of the identified variants to bone development.

CONCLUSION

Our findings provide insight into the common genetic architecture underlying BMD and BS, and enhance our understanding of the potential mechanism of osteoporosis fracture.

High carbohydrate and noodle/meat-rich dietary patterns interact with the minor haplotype in the 22q13 loci to increase its association with non-alcoholic fatty liver disease risk in Koreans

Non-alcoholic fatty liver(NAFLD) is prevalent in Asians despite the low obesity rate. We hypothesized that the haplotype of genetic variants in the 22q13 loci has a strong association with non-alcoholic fatty liver disease (NAFLD) that can be identified by genome-wide association study and that lifestyles may interact with the haplotype. We tested the hypothesis in middle-aged and elderly adults in a large city hospital-based cohort from the KoGES study. Men and women diagnosed with fatty liver, but who respectively consumed over 40 and 30 g ethanol per day were excluded. The haplotype of the selected SNPs from the 22q13 loci that influences NAFLD risk was generated. Among the 27374 participants, 1486 (5.4%) were diagnosed with NAFLD. LARGE_rs240072, RBFOX2_rs11089778, TRIOBP_rs12628603, PNPLA3_rs738409, and PARVB_rs2073080 in the 22q13 loci were included in the haplotype. Participants with the minor haplotype had 1.8, 2.3, and 1.8 times higher in the risk for NAFLD and serum AST and ALT activities, respectively, than those with the major haplotype. BMI, waist circumferences, serum glucose concentrations, and blood pressure interacted with the haplotype for NAFLD risk. We also found that a high carbohydrate intake and a dietary pattern characterized by high noodle and meat consumption significantly interacted with the minor haplotype to increase the risk of NAFLD. We hypothesized that the high incidence of NAFLD among Koreans, despite a relatively low incidence of obesity, might be due to genetic factors and perhaps their interactions with dietary patterns. The hypothesis was accepted since this study confirmed that participants with the minor allele of the haplotype in the 22q13 loci had a higher NAFLD risk that was exacerbated by high intakes of carbohydrates and a dietary pattern characterized by high noodle and meat consumption.

OncoBase: a platform for decoding regulatory somatic mutations in human cancers

Whole-exome and whole-genome sequencing have revealed millions of somatic mutations associated with different human cancers, and the vast majority of them are located outside of coding sequences, making it challenging to directly interpret their functional effects. With the rapid advances in high-throughput sequencing technologies, genome-scale long-range chromatin interactions were detected, and distal target genes of regulatory elements were determined using three-dimensional (3D) chromatin looping. Herein, we present OncoBase (http://www.oncobase.biols.ac.cn/), an integrated database for annotating 81 385 242 somatic mutations in 68 cancer types from more than 120 cancer projects by exploring their roles in distal interactions between target genes and regulatory elements. OncoBase integrates local chromatin signatures, 3D chromatin interactions in different cell types and reconstruction of enhancer-target networks using state-of-the-art algorithms. It employs informative visualization tools to display the integrated local and 3D chromatin signatures and effects of somatic mutations on regulatory elements. Enhancer-promoter interactions estimated from chromatin interactions are integrated into a network diffusion system that quantitatively prioritizes somatic mutations and target genes from a large pool. Thus, OncoBase is a useful resource for the functional annotation of regulatory noncoding regions and systematically benchmarking the regulatory effects of embedded noncoding somatic mutations in human carcinogenesis.

Genetic networks in Parkinson's and Alzheimer's disease

Parkinson’s disease (PD) and Alzheimer’s disease (AD) are the most common neurodegenerative diseases and there is increasing evidence that they share common physiological and pathological links. Here we have conducted the largest network analysis of PD and AD based on their gene expressions in blood to date. We identified modules that were not preserved between disease and healthy control (HC) networks, and important hub genes and transcription factors (TFs) in these modules. We highlighted that the PD module not preserved in HCs was associated with insulin resistance, and HDAC6 was identified as a hub gene in this module which may have the role of influencing tau phosphorylation and autophagic flux in neurodegenerative disease. The AD module associated with regulation of lipolysis in adipocytes and neuroactive ligand-receptor interaction was not preserved in healthy and mild cognitive impairment networks and the key hubs TRPC5 and BRAP identified as potential targets for therapeutic treatments of AD. Our study demonstrated that PD and AD share common disrupted genetics and identified novel pathways, hub genes and TFs that may be new areas for mechanistic study and important targets in both diseases.

Identification of key genes associated with multiple sclerosis based on gene expression data from peripheral blood mononuclear cells

The aim of this study was to identify the potential key candidate genes of multiple sclerosis (MS) and uncover mechanisms in MS. We combined data from the microarray expression profile of three MS stages and performed bioinformatics analysis. Differentially expressed genes (DEGs) were identified among the distinct stages of MS and healthy controls, and a total of 349 shared DEGs were identified. Gene ontology (GO) and pathway enrichment analyses showed that the DEGs were significantly enriched in the biological processes (BPs) of purine-related metabolic processes and signaling, especially the common DEGs, which were enriched in some immunological processes. Most of the DEGs were enriched in signaling pathways associated with the immune system, some immune diseases and infectious disease pathways. Through a protein-protein interaction (PPI) network analysis and a gene expression regulatory network constructed with MS-related miRNAs, we confirmed FOS, TP53, VEGFA, JUN, HIF1A, RB1, PTGS2, CXCL8, OAS2, NFKBIA and OAS1 as candidate genes of MS. Furthermore , we explored the potential SNPs associated with MS by database mining. In conclusion, this study provides the identified genes, SNPs, biological processes, and cellular pathways associated with MS. The uncovered candidate genes may be potential biomarkers involved in the diagnosis and therapy of MS.

Genomics of periodontal disease and tooth morbidity

In this review we critically summarize the evidence base and the progress to date regarding the genomic basis of periodontal disease and tooth morbidity (ie, dental caries and tooth loss), and discuss future applications and research directions in the context of precision oral health and care. Evidence for these oral/dental traits from genome-wide association studies first emerged less than a decade ago. Basic and translational research activities in this domain are now under way by multiple groups around the world. Key departure points in the oral health genomics discourse are: (a) some heritable variation exists for periodontal and dental diseases; (b) the environmental component (eg, social determinants of health and behavioral risk factors) has a major influence on the population distribution but probably interacts with factors of innate susceptibility at the person-level; (c) sizeable, multi-ethnic, well-characterized samples or cohorts with high-quality measures on oral health outcomes and genomics information are required to make decisive discoveries; (d) challenges remain in the measurement of oral health and disease, with current periodontitis and dental caries traits capturing only a part of the health-disease continuum, and are little or not informed by the underlying biology; (e) the substantial individual heterogeneity that exists in the clinical presentation and lifetime trajectory of oral disease can be identified and leveraged in a precision medicine framework or, if unappreciated, can hamper translational efforts. In this review we discuss how composite or biologically informed traits may offer improvements over clinically defined ones for the genomic interrogation of oral diseases. We demonstrate the utility of the results of genome-wide association studies for the development and testing of a genetic risk score for severe periodontitis. We conclude that exciting opportunities lie ahead for improvements in the oral health of individual patients and populations via advances in our understanding of the genomic basis of oral health and disease. The pace of new discoveries and their equitable translation to practice will largely depend on investments in the education and training of the oral health care workforce, basic and population research, and sustained collaborative efforts..

A phenome-wide association study (PheWAS) in the Population Architecture using Genomics and Epidemiology (PAGE) study reveals potential pleiotropy in African Americans

We performed a hypothesis-generating phenome-wide association study (PheWAS) to identify and characterize cross-phenotype associations, where one SNP is associated with two or more phenotypes, between thousands of genetic variants assayed on the Metabochip and hundreds of phenotypes in 5,897 African Americans as part of the Population Architecture using Genomics and Epidemiology (PAGE) I study. The PAGE I study was a National Human Genome Research Institute-funded collaboration of four study sites accessing diverse epidemiologic studies genotyped on the Metabochip, a custom genotyping chip that has dense coverage of regions in the genome previously associated with cardio-metabolic traits and outcomes in mostly European-descent populations. Here we focus on identifying novel phenome-genome relationships, where SNPs are associated with more than one phenotype. To do this, we performed a PheWAS, testing each SNP on the Metabochip for an association with up to 273 phenotypes in the participating PAGE I study sites. We identified 133 putative pleiotropic variants, defined as SNPs associated at an empirically derived p-value threshold of p<0.01 in two or more PAGE study sites for two or more phenotype classes. We further annotated these PheWAS-identified variants using publicly available functional data and local genetic ancestry. Amongst our novel findings is SPARC rs4958487, associated with increased glucose levels and hypertension. SPARC has been implicated in the pathogenesis of diabetes and is also known to have a potential role in fibrosis, a common consequence of multiple conditions including hypertension. The SPARC example and others highlight the potential that PheWAS approaches have in improving our understanding of complex disease architecture by identifying novel relationships between genetic variants and an array of common human phenotypes.

Exonic Variants in Aging-Related Genes Are Predictive of Phenotypic Aging Status

Background: Recent studies investigating longevity have revealed very few convincing genetic associations with increased lifespan. This is, in part, due to the complexity of biological aging, as well as the limited power of genome-wide association studies, which assay common single nucleotide polymorphisms (SNPs) and require several thousand subjects to achieve statistical significance. To overcome such barriers, we performed comprehensive DNA sequencing of a panel of 20 genes previously associated with phenotypic aging in a cohort of 200 individuals, half of whom were clinically defined by an "early aging" phenotype, and half of whom were clinically defined by a "late aging" phenotype based on age (65-75 years) and the ability to walk up a flight of stairs or walk for 15 min without resting. A validation cohort of 511 late agers was used to verify our results. Results: We found early agers were not enriched for more total variants in these 20 aging-related genes than late agers. Using machine learning methods, we identified the most predictive model of aging status, both in our discovery and validation cohorts, to be a random forest model incorporating damaging exon variants [Combined Annotation-Dependent Depletion (CADD) > 15]. The most heavily weighted variants in the model were within poly(ADP-ribose) polymerase 1 (PARP1) and excision repair cross complementation group 5 (ERCC5), both of which are involved in a canonical aging pathway, DNA damage repair. Conclusion: Overall, this study implemented a framework to apply machine learning to identify sequencing variants associated with complex phenotypes such as aging. While the small sample size making up our cohort inhibits our ability to make definitive conclusions about the ability of these genes to accurately predict aging, this study offers a unique method for exploring polygenic associations with complex phenotypes.

Systems Approach to Identify Common Genes and Pathways Associated with Response to Selective Serotonin Reuptake Inhibitors and Major Depression Risk

Despite numerous studies on major depressive disorder (MDD) susceptibility, the precise underlying molecular mechanism has not been elucidated which restricts the development of etiology-based disease-modifying drug. Major depressive disorder treatment is still symptomatic and is the leading cause of (~30%) failure of the current antidepressant therapy. Here we comprehended the probable genes and pathways commonly associated with antidepressant response and MDD. A systematic review was conducted, and candidate genes/pathways associated with antidepressant response and MDD were identified using an integrative genetics approach. Initially, single nucleotide polymorphisms (SNPs)/genes found to be significantly associated with antidepressant response were systematically reviewed and retrieved from the candidate studies and genome-wide association studies (GWAS). Also, significant variations concerning MDD susceptibility were extracted from GWAS only. We found 245 (Set A) and 800 (Set B) significantly associated genes with antidepressant response and MDD, respectively. Further, gene set enrichment analysis revealed the top five co-occurring molecular pathways (p ≤ 0.05) among the two sets of genes: Cushing syndrome, Axon guidance, cAMP signaling pathway, Insulin secretion, and Glutamatergic synapse, wherein all show a very close relation to synaptic plasticity. Integrative analyses of candidate gene and genome-wide association studies would enable us to investigate the putative targets for the development of disease etiology-based antidepressant that might be more promising than current ones.

Genomic landscape analyses of reprogrammed cells using integrative and non-integrative methods reveal variable cancer-associated alterations

Recent development of cell reprogramming technologies brought a major hope for future cell therapy applications by the use of these cells or their derivatives. For this purpose, one of the major requirements is the absence of genomic alterations generating a risk of cell transformation. Here we analyzed by microarray-based comparative genomic hybridization human iPSC generated by two non-integrative and one integrative method at pluripotent stage as well as in corresponding teratomas. We show that all iPSC lines exhibit copy number variations (CNV) of several genes deregulated in oncogenesis. These cancer-associated genomic alterations were more pronounced in virally programmed hiPSCs and their derivative teratoma as compared to those found in iPSC generated by mRNA-mediated reprogramming. Bioinformatics analysis showed the involvement of these genes in human leukemia and carcinoma. We conclude that genetic screening should become a standard procedure to ensure that hiPSCs are free from cancer-associated genomic alterations before clinical use.

Protocols, Methods, and Tools for Genome-Wide Association Studies (GWAS) of Dental Traits

Oral health and disease are known to be influenced by complex interactions between environmental (e.g., social and behavioral) factors and innate susceptibility. Although the exact contribution of genomics and other layers of "omics" to oral health is an area of active research, it is well established that the susceptibility to dental caries, periodontal disease, and other oral and craniofacial traits is substantially influenced by the human genome. A comprehensive understanding of these genomic factors is necessary for the realization of precision medicine in the oral health domain. To aid in this direction, the advent and increasing affordability of high-throughput genotyping has enabled the simultaneous interrogation of millions of genetic polymorphisms for association with oral and craniofacial traits. Specifically, genome-wide association studies (GWAS) of dental caries and periodontal disease have provided initial insights into novel loci and biological processes plausibly implicated in these two common, complex, biofilm-mediated diseases. This paper presents a summary of protocols, methods, tools, and pipelines for the conduct of GWAS of dental caries, periodontal disease, and related traits. The protocol begins with the consideration of different traits for both diseases and outlines procedures for genotyping, quality control, adjustment for population stratification, heritability and association analyses, annotation, reporting, and interpretation. Methods and tools available for GWAS are being constantly updated and improved; with this in mind, the presented approaches have been successfully applied in numerous GWAS and meta-analyses among tens of thousands of individuals, including dental traits such as dental caries and periodontal disease. As such, they can serve as a guide or template for future genomic investigations of these and other traits.

Meta-analysis of genome-wide association studies of aggressive and chronic periodontitis identifies two novel risk loci

Periodontitis is one of the most common inflammatory diseases, with a prevalence of 11% worldwide for the severe forms and an estimated heritability of 50%. It is classified into the widespread moderate form chronic periodontitis (CP) and the rare early-onset and severe phenotype aggressive periodontitis (AgP). These different disease manifestations are thought to share risk alleles and predisposing environmental factors. To obtain novel insights into the shared genetic etiology and the underlying molecular mechanisms of both forms, we performed a two step-wise meta-analysis approach using genome-wide association studies of both phenotypes. Genotypes from imputed genome-wide association studies (GWAS) of AgP and CP comprising 5,095 cases and 9,908 controls of North-West European genetic background were included. Two loci were associated with periodontitis at a genome-wide significance level. They located within the pseudogene MTND1P5 on chromosome 8 (rs16870060-G, P = 3.69 × 10^-9, OR = 1.36, 95% CI = [1.23-1.51]) and intronic of the long intergenic non-coding RNA LOC107984137 on chromosome 16, downstream of the gene SHISA9 (rs729876-T, P = 9.77 × 10^-9, OR = 1.24, 95% CI = [1.15-1.34]). This study identified novel risk loci of periodontitis, adding to the genetic basis of AgP and CP.

Genome-wide association meta-analysis of coronary artery disease and periodontitis reveals a novel shared risk locus

Evidence for a shared genetic basis of association between coronary artery disease (CAD) and periodontitis (PD) exists. To explore the joint genetic basis, we performed a GWAS meta-analysis. In the discovery stage, we used a German aggressive periodontitis sample (AgP-Ger; 680 cases vs 3,973 controls) and the CARDIoGRAMplusC4D CAD meta-analysis dataset (60,801 cases vs 123,504 controls). Two SNPs at the known CAD risk loci ADAMTS7 (rs11634042) and VAMP8 (rs1561198) passed the pre-assigned selection criteria (P_AgP-Ger < 0.05; P_CAD < 5 × 10⁻⁸; concordant effect direction) and were replicated in an independent GWAS meta-analysis dataset of PD (4,415 cases vs 5,935 controls). SNP rs1561198 showed significant association (PD[Replication]: P = 0.008 OR = 1.09, 95% CI = [1.02–1.16]; PD [Discovery + Replication]: P = 0.0002, OR = 1.11, 95% CI = [1.05–1.17]). For the associated haplotype block, allele specific cis-effects on VAMP8 expression were reported. Our data adds to the shared genetic basis of CAD and PD and indicate that the observed association of the two disease conditions cannot be solely explained by shared environmental risk factors. We conclude that the molecular pathway shared by CAD and PD involves VAMP8 function, which has a role in membrane vesicular trafficking, and is manipulated by pathogens to corrupt host immune defense.

Variation within voltage-gated calcium channel genes and antipsychotic treatment response in a South African first episode schizophrenia cohort

Voltage-gated calcium channels have been implicated in schizophrenia aetiology; however, little is known about their involvement in antipsychotic treatment response. This study investigated variants within the calcium channel subunit genes for association with antipsychotic treatment response in a first episode schizophrenia cohort. Twelve regulatory variants within seven genes were shown to be significantly associated with treatment outcome. Most notably, the CACNA1B rs2229949 CC genotype was associated with improved negative symptomology, where the C allele was predicted to abolish a miRNA-binding site (has-mir-5002-3p), suggesting a possible mechanism of action through which this variant may have an effect. These results implicate the calcium channel subunits in antipsychotic treatment response and suggest that increased activation of these channels may be explored to enhance or predict antipsychotic treatment outcome.

Genetic variant for behavioral regulation factor of executive function and its possible brain mechanism in attention deficit hyperactivity disorder

As a childhood-onset psychiatric disorder, attention deficit hyperactivity disorder (ADHD) is complicated by phenotypic and genetic heterogeneity. Lifelong executive function deficits in ADHD are described in many literatures and have been proposed as endophenotypes of ADHD. However, its genetic basis is still elusive. In this study, we performed a genome-wide association study of executive function, rated with Behavioral Rating Inventory of Executive Function (BRIEF), in ADHD children. We identified one significant variant (rs852004, P = 2.51e-08) for the overall score of BRIEF. The association analyses for each component of executive function found this locus was more associated with inhibit and monitor components. Further principle component analysis and confirmatory factor analysis provided an ADHD-specific executive function pattern including inhibit and monitor factors. SNP rs852004 was mainly associated with the Behavioral Regulation factor. Meanwhile, we found the significant locus was associated with ADHD symptom. The Behavioral Regulation factor mediated its effect on ADHD symptom. Functional magnetic resonance imaging (fMRI) analyses further showed evidence that this variant affected the activity of inhibition control related brain regions. It provided new insights for the genetic basis of executive function in ADHD.

Exon Array Biomarkers for the Differential Diagnosis of Schizophrenia and Bipolar Disorder

This study developed potential blood-based biomarker tests for diagnosing and differentiating schizophrenia (SZ), bipolar disorder type I (BD), and normal control (NC) subjects using mRNA gene expression signatures. A total of 90 subjects (n = 30 each for the three groups of subjects) provided blood samples at two visits. The Affymetrix exon microarray was used to profile the expression of over 1.4 million probesets. We selected potential biomarker panels using the temporal stability of the probesets and also back-tested them at two different visits for each subject. The 18-gene biomarker panels, using logistic regression modeling, correctly differentiated the three groups of subjects with high accuracy across the two different clinical visits (83-88% accuracy). The results are also consistent with the actual data and the "leave-one-out" analyses, indicating that the models should be predictive when applied to independent data cohorts. Many of the SZ and BD subjects were taking antipsychotic and mood stabilizer medications at the time of blood draw, raising the possibility that these drugs could have affected some of the differential transcription signatures. Using an independent Illumina data set of gene expression data from antipsychotic medication-free SZ subjects, the 18-gene biomarker panels produced a receiver operating characteristic curve accuracy greater than 0.866 in patients that were less than 30 years of age and medication free. We confirmed select transcripts by quantitative PCR and the nCounter® System. The episodic nature of psychiatric disorders might lead to highly variable results depending on when blood is collected in relation to the severity of the disease/symptoms. We have found stable trait gene panel markers for lifelong psychiatric disorders that may have diagnostic utility in younger undiagnosed subjects where there is a critical unmet need. The study requires replication in subjects for ultimate proof of the utility of the differential diagnosis.

Genetic Polymorphisms in Cytokine Genes in Colombian Patients with Ocular Toxoplasmosis

Toxoplasmosis is caused by infection with the protozoan parasite Toxoplasma gondii, which has the capacity to infect all warm-blooded animals worldwide. Toxoplasmosis is a major cause of visual defects in the Colombian population; however, the association between genetic polymorphisms in cytokine genes and susceptibility to ocular toxoplasmosis has not been studied in this population. This work evaluates the associations between polymorphisms in genes coding for the cytokines tumor necrosis factor alpha (TNF-α) (rs1799964, rs1800629, rs1799724, rs1800630, and rs361525), interleukin 1β (IL-1β) (rs16944, rs1143634, and rs1143627), IL-1α (rs1800587), gamma interferon (IFN-γ) (rs2430561), and IL-10 (rs1800896 and rs1800871) and the presence of ocular toxoplasmosis (OT) in a sample of a Colombian population (61 patients with OT and 116 healthy controls). Genotyping was performed with the "dideoxynucleotide (ddNTP) primer extension" technique. Functional-effect predictions of single nucleotide polymorphisms (SNPs) were done by using FuncPred. A polymorphism in the IL-10 gene promoter (-1082G/A) was significantly more prevalent in OT patients than in controls (P = 1.93e-08; odds ratio [OR] = 5.27e+03; 95% confidence interval [CI] = 3.18 to 8.739; Bonferroni correction [BONF] = 3.48e-07). In contrast, haplotype "AG" of the IL-10 gene promoter polymorphisms (rs1800896 and rs1800871) was present at a lower frequency in OT patients (P = 7e-04; OR = 0.10; 95% CI = 0.03 to 0.35). The +874A/T polymorphism of IFN-γ was associated with OT (P = 3.37e-05; OR = 4.2; 95% CI = 2.478 to 7.12; BONF = 6.07e-04). Haplotype "GAG" of the IL-1β gene promoter polymorphisms (rs1143634, rs1143627, and rs16944) appeared to be significantly associated with OT (P = 0.0494). The IL-10, IFN-γ, and IL-1β polymorphisms influence the development of OT in the Colombian population.

A novel variant associated with HDL-C levels by modifying DAGLB expression levels: An annotation-based genome-wide association study

Although numbers of genome-wide association studies (GWAS) have been performed for serum lipid levels, limited heritability has been explained. Studies showed that combining data from GWAS and expression quantitative trait loci (eQTLs) signals can both enhance the discovery of trait-associated SNPs and gain a better understanding of the mechanism. We performed an annotation-based, multistage genome-wide screening for serum-lipid-level-associated loci in totally 6863 Han Chinese. A serum high-density lipoprotein cholesterol (HDL-C) associated variant rs1880118 (hg19 chr7:g. 6435220G>C) was replicated (P_combined = 1.4E-10). rs1880118 was associated with DAGLB (diacylglycerol lipase, beta) expression levels in subcutaneous adipose tissue (P = 5.9E-42) and explained 47.7% of the expression variance. After the replication, an active segment covering variants tagged by rs1880118 near 5' of DAGLB was annotated using histone modification and transcription factor binding signals. The luciferase report assay revealed that the segment containing the minor alleles showed increased transcriptional activity compared with segment contains the major alleles, which was consistent with the eQTL analyses. The expression-trait association tests indicated the association between the DAGLB and serum HDL-C levels using gene-based approaches called "TWAS" (P = 3.0E-8), "SMR" (P = 1.1E-4), and "Sherlock" (P = 1.6E-6). To summarize, we identified a novel HDL-C-associated variant which explained nearly half of the expression variance of DAGLB. Integrated analyses established a genotype-gene-phenotype three-way association and expanded our knowledge of DAGLB in lipid metabolism.

Convergent downstream candidate mechanisms of independent intergenic polymorphisms between co-classified diseases implicate epistasis among noncoding elements

Eighty percent of DNA outside protein coding regions was shown biochemically functional by the ENCODE project, enabling studies of their interactions. Studies have since explored how convergent downstream mechanisms arise from independent genetic risks of one complex disease. However, the cross-talk and epistasis between intergenic risks associated with distinct complex diseases have not been comprehensively characterized. Our recent integrative genomic analysis unveiled downstream biological effectors of disease-specific polymorphisms buried in intergenic regions, and we then validated their genetic synergy and antagonism in distinct GWAS. We extend this approach to characterize convergent downstream candidate mechanisms of distinct intergenic SNPs across distinct diseases within the same clinical classification. We construct a multipartite network consisting of 467 diseases organized in 15 classes, 2,358 disease-associated SNPs, 6,301 SNPassociated mRNAs by eQTL, and mRNA annotations to 4,538 Gene Ontology mechanisms. Functional similarity between two SNPs (similar SNP pairs) is imputed using a nested information theoretic distance model for which p-values are assigned by conservative scale-free permutation of network edges without replacement (node degrees constant). At FDR≤5%, we prioritized 3,870 intergenic SNP pairs associated, among which 755 are associated with distinct diseases sharing the same disease class, implicating 167 intergenic SNPs, 14 classes, 230 mRNAs, and 134 GO terms. Co-classified SNP pairs were more likely to be prioritized as compared to those of distinct classes confirming a noncoding genetic underpinning to clinical classification (odds ratio ∼3.8; p≤10-25). The prioritized pairs were also enriched in regions bound to the same/interacting transcription factors and/or interacting in long-range chromatin interactions suggestive of epistasis (odds ratio ∼ 2,500; p≤10-25). This prioritized network implicates complex epistasis between intergenic polymorphisms of co-classified diseases and offers a roadmap for a novel therapeutic paradigm: repositioning medications that target proteins within downstream mechanisms of intergenic disease-associated SNPs. Supplementary information and software: http://lussiergroup.org/publications/disease_class.

Functional mapping and annotation of genetic associations with FUMA

A main challenge in genome-wide association studies (GWAS) is to pinpoint possible causal variants. Results from GWAS typically do not directly translate into causal variants because the majority of hits are in non-coding or intergenic regions, and the presence of linkage disequilibrium leads to effects being statistically spread out across multiple variants. Post-GWAS annotation facilitates the selection of most likely causal variant(s). Multiple resources are available for post-GWAS annotation, yet these can be time consuming and do not provide integrated visual aids for data interpretation. We, therefore, develop FUMA: an integrative web-based platform using information from multiple biological resources to facilitate functional annotation of GWAS results, gene prioritization and interactive visualization. FUMA accommodates positional, expression quantitative trait loci (eQTL) and chromatin interaction mappings, and provides gene-based, pathway and tissue enrichment results. FUMA results directly aid in generating hypotheses that are testable in functional experiments aimed at proving causal relations.

Prioritizing genetic variants is a major challenge in genome-wide association studies. Here, the authors develop FUMA, a web-based bioinformatics tool that uses a combination of positional, eQTL and chromatin interaction mapping to prioritize likely causal variants and genes.

Genome-wide quantitative trait loci mapping of the human cerebrospinal fluid proteome

Cerebrospinal fluid (CSF) is virtually the only one accessible source of proteins derived from the central nervous system (CNS) of living humans and possibly reflects the pathophysiology of a variety of neuropsychiatric diseases. However, little is known regarding the genetic basis of variation in protein levels of human CSF. We examined CSF levels of 1,126 proteins in 133 subjects and performed a genome-wide association analysis of 514,227 single nucleotide polymorphisms (SNPs) to detect protein quantitative trait loci (pQTLs). To be conservative, Spearman's correlation was used to identify an association between genotypes of SNPs and protein levels. A total of 421 cis and 25 trans SNP-protein pairs were significantly correlated at a false discovery rate (FDR) of less than 0.01 (nominal P < 7.66 × 10-9). Cis-only analysis revealed additional 580 SNP-protein pairs with FDR < 0.01 (nominal P < 2.13 × 10-5). pQTL SNPs were more likely, compared to non-pQTL SNPs, to be a disease/trait-associated variants identified by previous genome-wide association studies. The present findings suggest that genetic variations play an important role in the regulation of protein expression in the CNS. The obtained database may serve as a valuable resource to understand the genetic bases for CNS protein expression pattern in humans.

Biological networks in Parkinson's disease: an insight into the epigenetic mechanisms associated with this disease

BACKGROUND

Parkinson's disease (PD) is the second most prevalent neurodegenerative disorders in the world. Studying PD from systems biology perspective involving genes and their regulators might provide deeper insights into the complex molecular interactions associated with this disease.

RESULT

We have studied gene co-expression network obtained from a PD-specific microarray data. The co-expression network identified 11 hub genes, of which eight genes are not previously known to be associated with PD. Further study on the functionality of these eight novel hub genes revealed that these genes play important roles in several neurodegenerative diseases. Furthermore, we have studied the tissue-specific expression and histone modification patterns of the novel hub genes. Most of these genes possess several histone modification sites those are already known to be associated with neurodegenerative diseases. Regulatory network namely mTF-miRNA-gene-gTF involves microRNA Transcription Factor (mTF), microRNA (miRNA), gene and gene Transcription Factor (gTF). Whereas long noncoding RNA (lncRNA) mediated regulatory network involves miRNA, gene, mTF and lncRNA. mTF-miRNA-gene-gTF regulatory network identified a novel feed-forward loop. lncRNA-mediated regulatory network identified novel lncRNAs of PD and revealed the two-way regulatory pattern of PD-specific miRNAs where miRNAs can be regulated by both the TFs and lncRNAs. SNP analysis of the most significant genes of the co-expression network identified 20 SNPs. These SNPs are present in the 3' UTR of known PD genes and are controlled by those miRNAs which are also involved in PD.

CONCLUSION

Our study identified eight novel hub genes which can be considered as possible candidates for future biomarker identification studies for PD. The two regulatory networks studied in our work provide a detailed overview of the cellular regulatory mechanisms where the non-coding RNAs namely miRNA and lncRNA, can act as epigenetic regulators of PD. SNPs identified in our study can be helpful for identifying PD at an earlier stage. Overall, this study may impart a better comprehension of the complex molecular interactions associated with PD from systems biology perspective.

Multiple phenotype association tests using summary statistics in genome-wide association studies

We study in this article jointly testing the associations of a genetic variant with correlated multiple phenotypes using the summary statistics of individual phenotype analysis from Genome-Wide Association Studies (GWASs). We estimated the between-phenotype correlation matrix using the summary statistics of individual phenotype GWAS analyses, and developed genetic association tests for multiple phenotypes by accounting for between-phenotype correlation without the need to access individual-level data. Since genetic variants often affect multiple phenotypes differently across the genome and the between-phenotype correlation can be arbitrary, we proposed robust and powerful multiple phenotype testing procedures by jointly testing a common mean and a variance component in linear mixed models for summary statistics. We computed the p-values of the proposed tests analytically. This computational advantage makes our methods practically appealing in large-scale GWASs. We performed simulation studies to show that the proposed tests maintained correct type I error rates, and to compare their powers in various settings with the existing methods. We applied the proposed tests to a GWAS Global Lipids Genetics Consortium summary statistics data set and identified additional genetic variants that were missed by the original single-trait analysis.

An independent component analysis confounding factor correction framework for identifying broad impact expression quantitative trait loci

Genome-wide expression Quantitative Trait Loci (eQTL) studies in humans have provided numerous insights into the genetics of both gene expression and complex diseases. While the majority of eQTL identified in genome-wide analyses impact a single gene, eQTL that impact many genes are particularly valuable for network modeling and disease analysis. To enable the identification of such broad impact eQTL, we introduce CONFETI: Confounding Factor Estimation Through Independent component analysis. CONFETI is designed to address two conflicting issues when searching for broad impact eQTL: the need to account for non-genetic confounding factors that can lower the power of the analysis or produce broad impact eQTL false positives, and the tendency of methods that account for confounding factors to model broad impact eQTL as non-genetic variation. The key advance of the CONFETI framework is the use of Independent Component Analysis (ICA) to identify variation likely caused by broad impact eQTL when constructing the sample covariance matrix used for the random effect in a mixed model. We show that CONFETI has better performance than other mixed model confounding factor methods when considering broad impact eQTL recovery from synthetic data. We also used the CONFETI framework and these same confounding factor methods to identify eQTL that replicate between matched twin pair datasets in the Multiple Tissue Human Expression Resource (MuTHER), the Depression Genes Networks study (DGN), the Netherlands Study of Depression and Anxiety (NESDA), and multiple tissue types in the Genotype-Tissue Expression (GTEx) consortium. These analyses identified both cis-eQTL and trans-eQTL impacting individual genes, and CONFETI had better or comparable performance to other mixed model confounding factor analysis methods when identifying such eQTL. In these analyses, we were able to identify and replicate a few broad impact eQTL although the overall number was small even when applying CONFETI. In light of these results, we discuss the broad impact eQTL that have been previously reported from the analysis of human data and suggest that considerable caution should be exercised when making biological inferences based on these reported eQTL.

Genome-wide association screens for Achilles tendon and ACL tears and tendinopathy

Achilles tendinopathy or rupture and anterior cruciate ligament (ACL) rupture are substantial injuries affecting athletes, associated with delayed recovery or inability to return to competition. To identify genetic markers that might be used to predict risk for these injuries, we performed genome-wide association screens for these injuries using data from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort consisting of 102,979 individuals. We did not find any single nucleotide polymorphisms (SNPs) associated with either of these injuries with a p-value that was genome-wide significant (p<5x10-8). We found, however, four and three polymorphisms with p-values that were borderline significant (p<10-6) for Achilles tendon injury and ACL rupture, respectively. We then tested SNPs previously reported to be associated with either Achilles tendon injury or ACL rupture. None showed an association in our cohort with a false discovery rate of less than 5%. We obtained, however, moderate to weak evidence for replication in one case; specifically, rs4919510 in MIR608 had a p-value of 5.1x10-3 for association with Achilles tendon injury, corresponding to a 7% chance of false replication. Finally, we tested 2855 SNPs in 90 candidate genes for musculoskeletal injury, but did not find any that showed a significant association below a false discovery rate of 5%. We provide data containing summary statistics for the entire genome, which will be useful for future genetic studies on these injuries.

A candidate transacting modulator of fetal hemoglobin gene expression in the Arab-Indian haplotype of sickle cell anemia

Fetal hemoglobin (HbF) levels are higher in the Arab-Indian (AI) β-globin gene haplotype of sickle cell anemia compared with African-origin haplotypes. To study genetic elements that effect HbF expression in the AI haplotype we completed whole genome sequencing in 14 Saudi AI haplotype sickle hemoglobin homozygotes-seven selected for low HbF (8.2% ± 1.3%) and seven selected for high HbF (23.5% ± 2.6%). An intronic single nucleotide polymorphism (SNP) in ANTXR1, an anthrax toxin receptor (chromosome 2p13), was associated with HbF. These results were replicated in two independent Saudi AI haplotype cohorts of 120 and 139 patients, but not in 76 Saudi Benin haplotype, 894 African origin haplotype and 44 AI haplotype patients of Indian origin, suggesting that this association is effective only in the Saudi AI haplotype background. ANTXR1 variants explained 10% of the HbF variability compared with 8% for BCL11A. These two genes had independent, additive effects on HbF and together explained about 15% of HbF variability in Saudi AI sickle cell anemia patients. ANTXR1 was expressed at mRNA and protein levels in erythroid progenitors derived from induced pluripotent stem cells (iPSCs) and CD34⁺ cells. As CD34⁺ cells matured and their HbF decreased ANTXR1 expression increased; as iPSCs differentiated and their HbF increased, ANTXR1 expression decreased. Along with elements in cis to the HbF genes, ANTXR1 contributes to the variation in HbF in Saudi AI haplotype sickle cell anemia and is the first gene in trans to HBB that is associated with HbF only in carriers of the Saudi AI haplotype. Am. J. Hematol. 91:1118-1122, 2016. © 2016 Wiley Periodicals, Inc.

A combined analysis of genome-wide expression profiling of bipolar disorder in human prefrontal cortex

Numbers of gene expression profiling studies of bipolar disorder have been published. Besides different array chips and tissues, variety of the data processes in different cohorts aggravated the inconsistency of results of these genome-wide gene expression profiling studies. By searching the gene expression databases, we obtained six data sets for prefrontal cortex (PFC) of bipolar disorder with raw data and combinable platforms. We used standardized pre-processing and quality control procedures to analyze each data set separately and then combined them into a large gene expression matrix with 101 bipolar disorder subjects and 106 controls. A standard linear mixed-effects model was used to calculate the differentially expressed genes (DEGs). Multiple levels of sensitivity analyses and cross validation with genetic data were conducted. Functional and network analyses were carried out on basis of the DEGs. In the result, we identified 198 unique differentially expressed genes in the PFC of bipolar disorder and control. Among them, 115 DEGs were robust to at least three leave-one-out tests or different pre-processing methods; 51 DEGs were validated with genetic association signals. Pathway enrichment analysis showed these DEGs were related with regulation of neurological system, cell death and apoptosis, and several basic binding processes. Protein-protein interaction network further identified one key hub gene. We have contributed the most comprehensive integrated analysis of bipolar disorder expression profiling studies in PFC to date. The DEGs, especially those with multiple validations, may denote a common signature of bipolar disorder and contribute to the pathogenesis of disease.

A childhood acute lymphoblastic leukemia genome-wide association study identifies novel sex-specific risk variants

Childhood acute lymphoblastic leukemia (ALL) occurs more frequently in males. Reasons behind sex differences in childhood ALL risk are unknown. In the present genome-wide association study (GWAS), we explored the genetic basis of sex differences by comparing genotype frequencies between male and female cases in a case-only study to assess effect-modification by sex.The case-only design included 236 incident cases of childhood ALL consecutively recruited at the Texas Children's Cancer Center in Houston, Texas from 2007 to 2012. All cases were non-Hispanic whites, aged 1 to 10 years, and diagnosed with confirmed B-cell precursor ALL. Genotyping was performed using the Illumina HumanCoreExome BeadChip on the Illumina Infinium platform. Besides the top 100 statistically most significant results, results were also analyzed by the top 100 highest effect size with a nominal statistical significance (P <0.05).The statistically most significant sex-specific association (P = 4 × 10) was with the single nucleotide polymorphism (SNP) rs4813720 (RASSF2), an expression quantitative trait locus (eQTL) for RASSF2 in peripheral blood. rs4813720 is also a strong methylation QTL (meQTL) for a CpG site (cg22485289) within RASSF2 in pregnancy, at birth, childhood, and adolescence. cg22485289 is one of the hypomethylated CpG sites in ALL compared with pre-B cells. Two missense SNPs, rs12722042 and 12722039, in the HLA-DQA1 gene yielded the highest effect sizes (odds ratio [OR] ∼ 14; P <0.01) for sex-specific results. The HLA-DQA1 SNPs belong to DQA1*01 and confirmed the previously reported male-specific association with DQA1*01. This finding supports the proposed infection-related etiology in childhood ALL risk for males. Further analyses revealed that most SNPs (either direct effect or through linkage disequilibrium) were within active enhancers or active promoter regions and had regulatory effects on gene expression levels.Cumulative data suggested that RASSF2 rs4813720, which correlates with increased RASSF2 expression, may counteract the suppressor effect of estrogen-regulated miR-17-92 on RASSF2 resulting in protection in males. Given the amount of sex hormone-related mechanisms suggested by our findings, future studies should examine prenatal or early postnatal programming by sex hormones when hormone levels show a large variation.

Polymorphisms involving gain or loss of CpG sites are significantly enriched in trait-associated SNPs

Some single nucleotide polymorphisms (SNPs) influence the existence of CpG sites, the basis of DNA modification such as methylation and hydroxymethylation. These polymorphisms can lead to gain or loss of CpG sites and were defined as CpG site related SNPs (cgSNPs) in this study. The cgSNPs change DNA sequence and might potentially affect DNA modification such as methylation. However, the functional consequence of cgSNPs is poorly understood. We observed that a considerable proportion (23.0%) of common variants were cgSNPs in human genome. Mutations involving loss of CpG sites were associated with reduced levels of methylation (~20.2%) using The Cancer Genome Atlas (TCGA) data. Using public databases (SCAN and seeQTL) of expression quantitative trait loci (eQTLs), we found that the cgSNPs were significantly enriched in eQTLs via logistic regression and simulation test. Furthermore, we observed that cgSNPs were more likely to be trait-associated loci especially cancers using a catalog of published genome-wide association studies (GWAS) recorded by National Human Genome Research Institute (NHGRI). Our results indicated that cgSNP might be meaningful as annotation either in SNP functional prediction or in screening for trait-associated SNPs.

Identification of genetic modifiers of age-at-onset for familial Parkinson's disease

Parkinson's disease (PD) is the most common cause of neurodegenerative movement disorder and the second most common cause of dementia. Genes are thought to have a stronger effect on age-at-onset of PD than on risk, yet there has been a phenomenal success in identifying risk loci but not age-at-onset modifiers. We conducted a genome-wide study for age-at-onset. We analysed familial and non-familial PD separately, per prior evidence for strong genetic effect on age-at-onset in familial PD. GWAS was conducted in 431 unrelated PD individuals with at least one affected relative (familial PD) and 1544 non-familial PD from the NeuroGenetics Research Consortium (NGRC); an additional 737 familial PD and 2363 non-familial PD were used for replication. In familial PD, two signals were detected and replicated robustly: one mapped to LHFPL2 on 5q14.1 (PNGRC = 3E-8, PReplication = 2E-5, PNGRC + Replication = 1E-11), the second mapped to TPM1 on 15q22.2 (PNGRC = 8E-9, PReplication = 2E-4, PNGRC + Replication = 9E-11). The variants that were associated with accelerated onset had low frequencies (<0.02). The LHFPL2 variant was associated with earlier onset by 12.33 [95% CI: 6.2; 18.45] years in NGRC, 8.03 [2.95; 13.11] years in replication, and 9.79 [5.88; 13.70] years in the combined data. The TPM1 variant was associated with earlier onset by 15.30 [8.10; 22.49] years in NGRC, 9.29 [1.79; 16.79] years in replication, and 12.42 [7.23; 17.61] years in the combined data. Neither LHFPL2 nor TPM1 was associated with age-at-onset in non-familial PD. LHFPL2 (function unknown) is overexpressed in brain tumours. TPM1 encodes a highly conserved protein that regulates muscle contraction, and is a tumour-suppressor gene.

Phenome-Wide Association Study to Explore Relationships between Immune System Related Genetic Loci and Complex Traits and Diseases

We performed a Phenome-Wide Association Study (PheWAS) to identify interrelationships between the immune system genetic architecture and a wide array of phenotypes from two de-identified electronic health record (EHR) biorepositories. We selected variants within genes encoding critical factors in the immune system and variants with known associations with autoimmunity. To define case/control status for EHR diagnoses, we used International Classification of Diseases, Ninth Revision (ICD-9) diagnosis codes from 3,024 Geisinger Clinic MyCode^® subjects (470 diagnoses) and 2,899 Vanderbilt University Medical Center BioVU biorepository subjects (380 diagnoses). A pooled-analysis was also carried out for the replicating results of the two data sets. We identified new associations with potential biological relevance including SNPs in tumor necrosis factor (TNF) and ankyrin-related genes associated with acute and chronic sinusitis and acute respiratory tract infection. The two most significant associations identified were for the C6orf10 SNP rs6910071 and “rheumatoid arthritis” (ICD-9 code category 714) (p_METAL = 2.58 x 10⁻⁹) and the ATN1 SNP rs2239167 and “diabetes mellitus, type 2” (ICD-9 code category 250) (p_METAL = 6.39 x 10⁻⁹). This study highlights the utility of using PheWAS in conjunction with EHRs to discover new genotypic-phenotypic associations for immune-system related genetic loci.

Genomic approaches for understanding the genetics of complex disease

There are thousands of known associations between genetic variants and complex human phenotypes, and the rate of novel discoveries is rapidly increasing. Translating those associations into knowledge of disease mechanisms remains a fundamental challenge because the associated variants are overwhelmingly in noncoding regions of the genome where we have few guiding principles to predict their function. Intersecting the compendium of identified genetic associations with maps of regulatory activity across the human genome has revealed that phenotype-associated variants are highly enriched in candidate regulatory elements. Allele-specific analyses of gene regulation can further prioritize variants that likely have a functional effect on disease mechanisms; and emerging high-throughput assays to quantify the activity of candidate regulatory elements are a promising next step in that direction. Together, these technologies have created the ability to systematically and empirically test hypotheses about the function of noncoding variants and haplotypes at the scale needed for comprehensive and systematic follow-up of genetic association studies. Major coordinated efforts to quantify regulatory mechanisms across genetically diverse populations in increasingly realistic cell models would be highly beneficial to realize that potential.

Variants in Adjacent Oxytocin/Vasopressin Gene Region and Associations with ASD Diagnosis and Other Autism Related Endophenotypes

BACKGROUND

There has been increasing interest in oxytocin (peptide: OT, gene: OXT) as a treatment pathway for neurodevelopmental disorders such as Autism Spectrum Disorder (ASD). Neurodevelopmental disorders affect functional, social, and intellectual abilities. With advances in molecular biology, research has connected multiple gene regions to the clinical presentation of ASD. Studies have also shown that the neuropeptide hormones OT and arginine vasopressin (AVP) influence mammalian social and territorial behaviors and may have treatment potential for neurodevelopmental disorders. Published data examining molecular and phenotypic variation in ASD, such as cognitive abilities, are limited. Since most studies have focused on the receptors in the OT-AVP system, we investigated genetic variation within peptide genes for association with phenotypic ASD features that help identify subgroups within the spectrum.

METHODS

In this study, TDT analysis was carried out utilizing FBAT in 207 probands (156 trios) and a European Ancestry (EA) subsample (108 trios).The evolutionarily related and adjacent genes of OXT and AVP were studied for associations between the tagged single nucleotide polymorphisms and ASD diagnosis, social abilities, restrictive and repetitive behaviors, and IQ for cognitive abilities. Additionally, relationships with whole blood serotonin (WB5HT) were explored because of the developmental relationships connecting plasma levels of OT and WB5HT within ASD.

RESULTS

RESULTS indicate significant association between OXT rs6084258 (p = 0.001) and ASD. Associations with several endophenotypes were also noted: OXT rs6133010 was associated with IQ (full scale IQ, p = 0.008; nonverbal IQ, p = 0.010, verbal IQ, p = 0.006); and OXT rs4813625 and OXT rs877172 were associated with WB5HT levels (EA, p = 0.027 and p = 0.033, respectively). Additionally, we measured plasma OT (pOT) levels in a subsample (N = 54). RESULTS show the three polymorphisms, OXT rs6084258, OXT rs11697250, and OXT rs877172, have significant association with pOT (EA, p = 0.011, p = 0.010, and p = 0.002, respectively).

CONCLUSIONS

These findings suggest that SNPs near OXT and AVP are associated with diagnosis of ASD, social behaviors, restricted and repetitive behaviors, IQ, pOT, and WB5HT. Future studies need to replicate these findings and examine gene-interactions in other neurodevelopmental disorders. Mechanisms of action may influence early social and cognitive development that may or may not be limited to ASD diagnosis.

Integrative genomics analyses unveil downstream biological effectors of disease-specific polymorphisms buried in intergenic regions

Functionally altered biological mechanisms arising from disease-associated polymorphisms, remain difficult to characterise when those variants are intergenic, or, fall between genes. We sought to identify shared downstream mechanisms by which inter- and intragenic single-nucleotide polymorphisms (SNPs) contribute to a specific physiopathology. Using computational modelling of 2 million pairs of disease-associated SNPs drawn from genome-wide association studies (GWAS), integrated with expression Quantitative Trait Loci (eQTL) and Gene Ontology functional annotations, we predicted 3,870 inter–intra and inter–intra SNP pairs with convergent biological mechanisms (FDR<0.05). These prioritised SNP pairs with overlapping messenger RNA targets or similar functional annotations were more likely to be associated with the same disease than unrelated pathologies (OR>12). We additionally confirmed synergistic and antagonistic genetic interactions for a subset of prioritised SNP pairs in independent studies of Alzheimer’s disease (entropy P=0.046), bladder cancer (entropy P=0.039), and rheumatoid arthritis (PheWAS case–control P<10⁻⁴). Using ENCODE data sets, we further statistically validated that the biological mechanisms shared within prioritised SNP pairs are frequently governed by matching transcription factor binding sites and long-range chromatin interactions. These results provide a ‘roadmap’ of disease mechanisms emerging from GWAS and further identify candidate therapeutic targets among downstream effectors of intergenic SNPs.