Predicting cell types and genetic variations contributing to disease by combining GWAS and epigenetic data

Genome-Wide Association Study of Growth and Sex Traits Provides Insight into Heritable Mechanisms Underlying Growth Development of Macrobrachium nipponense (Oriental River Prawn)

Male hybrid oriental river prawns grow significantly faster than hybrid females. In this study, the growth and sex traits of 181 individuals of Macrobrachium nipponense were recorded, and each individual genotype was evaluated using the 2b-RAD sequencing method. The genetic parameters for growth and sex traits were estimated. A genome-wide association analysis (GWAS) of these traits was performed. In total, 18 growth-related SNPs were detected from 12 chromosomes using a mixed linear model. The most significant loci of weight are located on the position of the SNP (102638935, chromosome 13), which can explain 11.87% of the phenotypic variation. A total of 11 significant SNPs were detected on four chromosomes associated with sex trait (three on chromosome 4, one on chromosome 7 and seven on chromosome 17). The heritability of this trait is 0.8998 and belongs to the range of ultra-high heritability. Genetic correlations were prevalent among the 11 traits examined, the genetic coefficient between sex and body weight reached a significant level of −0.23. This study is the first GWAS for sex of binary and growth traits in oriental river prawn. Our results provide a set of markers for the genetic selection of growth traits and help us to further understand the genetic mechanisms of growth in Macrobrachium nipponense.

Genome-wide DNA methylation profiling after Ayurveda intervention to bronchial asthmatics identifies differential methylation in several transcription factors with immune process related function

BACKGROUND

The Indian traditional medicinal system, Ayurveda, describes several lifestyle practices, processes and medicines as an intervention to treat asthma. Rasayana therapy is one of them and although these treatment modules show improvement in bronchial asthma, their mechanism of action, particularly the effect on DNA methylation, is largely understudied.

OBJECTIVES

Our study aimed at identifying the contribution of DNA methylation changes in modulating bronchial asthma phenotype upon Ayurveda intervention.

MATERIALS AND METHODS

In this study, genome-wide methylation profiling in peripheral blood DNA of healthy controls and bronchial asthmatics before (BT) and after (AT) Ayurveda treatment was performed using array-based profiling of reference-independent methylation status (aPRIMES) coupled to microarray technique.

RESULTS

We identified 4820 treatment-associated DNA methylation signatures (TADS) and 11,643 asthma-associated DNA methylation signatures (AADS), differentially methylated [FDR (≤0.1) adjusted p-values] in AT and HC groups respectively, compared to BT group. Neurotrophin TRK receptor signaling pathway was significantly enriched for differentially methylated genes in bronchial asthmatics, compared to AT and HC subjects. Additionally, we identified over 100 differentially methylated immune-related genes located in the promoter/5'-UTR regions of TADS and AADS. Various immediate-early response and immune regulatory genes with functions such as transcription factor activity (FOXD1, FOXD2, GATA6, HOXA3, HOXA5, MZF1, NFATC1, NKX2-2, NKX2-3, RUNX1, KLF11), G-protein coupled receptor activity (CXCR4, PTGER4), G-protein coupled receptor binding (UCN), DNA binding (JARID2, EBF2, SOX9), SNARE binding (CAPN10), transmembrane signaling receptor activity (GP1BB), integrin binding (ITGA6), calcium ion binding (PCDHGA12), actin binding (TRPM7, PANX1, TPM1), receptor tyrosine kinase binding (PIK3R2), receptor activity (GDNF), histone methyltransferase activity (MLL5), and catalytic activity (TSTA3) were found to show consistent methylation status between AT and HC group in microarray data.

CONCLUSIONS

Our study reports the DNA methylation-regulated genes in bronchial asthmatics showing improvement in symptoms after Ayurveda intervention. DNA methylation regulation in the identified genes and pathways represents the Ayurveda intervention responsive genes and may be further explored as diagnostic, prognostic, and therapeutic biomarkers for bronchial asthma in peripheral blood.

chromMAGMA: regulatory element-centric interrogation of risk variants

Candidate causal risk variants from genome-wide association studies reside almost exclusively in noncoding regions of the genome and innovative approaches are necessary to understand their biological function. Multi-marker analysis of genomic annotation (MAGMA) is a widely used program that nominates candidate risk genes by mapping single-nucleotide polymorphism summary statistics from genome-wide association studies to gene bodies. We augmented MAGMA to create chromatin-MAGMA (chromMAGMA), a method to nominate candidate risk genes based on the presence of risk variants within noncoding regulatory elements (REs). We applied chromMAGMA to a genetic susceptibility dataset for epithelial ovarian cancer (EOC), a rare gynecologic malignancy characterized by high mortality. This identified 155 unique candidate EOC risk genes across five EOC histotypes; 83% (105/127) of high-grade serous ovarian cancer risk genes had not previously been implicated in this EOC histotype. Risk genes nominated by chromMAGMA converged on mRNA splicing and transcriptional dysregulation pathways. chromMAGMA is a pipeline that nominates candidate risk genes through a gene regulation-focused approach and helps interpret the biological mechanism of noncoding risk variants for complex diseases.

Multi-omics colocalization with genome-wide association studies reveals a context-specific genetic mechanism at a childhood onset asthma risk locus

BACKGROUND

Genome-wide association studies (GWASs) have identified thousands of variants associated with asthma and other complex diseases. However, the functional effects of most of these variants are unknown. Moreover, GWASs do not provide context-specific information on cell types or environmental factors that affect specific disease risks and outcomes. To address these limitations, we used an upper airway epithelial cell (AEC) culture model to assess transcriptional and epigenetic responses to rhinovirus (RV), an asthma-promoting pathogen, and provide context-specific functional annotations to variants discovered in GWASs of asthma.

METHODS

Genome-wide genetic, gene expression, and DNA methylation data in vehicle- and RV-treated upper AECs were collected from 104 individuals who had a diagnosis of airway disease (n=66) or were healthy participants (n=38). We mapped cis expression and methylation quantitative trait loci (cis-eQTLs and cis-meQTLs, respectively) in each treatment condition (RV and vehicle) in AECs from these individuals. A Bayesian test for colocalization between AEC molecular QTLs and adult onset asthma and childhood onset asthma GWAS SNPs, and a multi-ethnic GWAS of asthma, was used to assign the function to variants associated with asthma. We used Mendelian randomization to demonstrate DNA methylation effects on gene expression at asthma colocalized loci.

RESULTS

Asthma and allergic disease-associated GWAS SNPs were specifically enriched among molecular QTLs in AECs, but not in GWASs from non-immune diseases, and in AEC eQTLs, but not among eQTLs from other tissues. Colocalization analyses of AEC QTLs with asthma GWAS variants revealed potential molecular mechanisms of asthma, including QTLs at the TSLP locus that were common to both the RV and vehicle treatments and to both childhood onset and adult onset asthma, as well as QTLs at the 17q12-21 asthma locus that were specific to RV exposure and childhood onset asthma, consistent with clinical and epidemiological studies of these loci.

CONCLUSIONS

This study provides evidence of functional effects for asthma risk variants in AECs and insight into RV-mediated transcriptional and epigenetic response mechanisms that modulate genetic effects in the airway and risk for asthma.

Variant to Gene Mapping to Discover New Targets for Immune Tolerance

The breakdown of immunological tolerance leads to autoimmune disease, and the mechanisms that maintain self-tolerance, especially in humans, are not fully understood. Genome-wide association studies (GWAS) have identified hundreds of human genetic loci statistically linked to autoimmune disease risk, and epigenetic modifications of DNA and chromatin at these loci have been associated with autoimmune disease risk. Because the vast majority of these signals are located far from genes, identifying causal variants, and their functional consequences on the correct effector genes, has been challenging. These limitations have hampered the translation of GWAS findings into novel drug targets and clinical interventions, but recent advances in understanding the spatial organization of the genome in the nucleus have offered mechanistic insights into gene regulation and answers to questions left open by GWAS. Here we discuss the potential for 'variant-to-gene mapping' approaches that integrate GWAS with 3D functional genomic data to identify human genes involved in the maintenance of tolerance.

Comprehensive functional annotation of susceptibility variants associated with asthma

Genome-wide association studies (GWAS) have identified hundreds of primarily non-coding disease-susceptibility variants that further need functional interpretation to prioritize and discriminate the disease-relevant variants. We present a comprehensive genome-wide non-coding variant prioritization scheme followed by validation using Pyrosequencing and TaqMan assays in asthma. We implemented a composite Functional Annotation Score (cFAS) to investigate over 32,000 variants consisting of 1525 GWAS-lead asthma-susceptibility variants and their LD proxies (r² ≥ 0.80). Functional annotation pipeline in cFAS revealed 274 variants with significant score at 1% false discovery rate. This study implicates a novel locus 4p16 (SLC26A1) with eQTL variant (rs11936407) and known loci in 17q12-21 and 5q22 which encode ORM1-like protein 3 (ORMDL3, rs406527, and rs12936231) and thymic stromal lymphopoietin (TSLP, rs3806932 and rs10073816) epithelial gene, respectively. Follow-up validation analysis through pyrosequencing of CpG sites in and nearby rs4065275 and rs11936407 showed genotype-dependent hypomethylation on asthma cases compared with healthy controls. Prioritized variants are enriched for asthma-specific histone modification associated with active chromatin (H3K4me1 and H3K27ac) in T cells, B cells, lung, and immune-related interferon gamma signaling pathways. Our findings, together with those from prior studies, suggest that SNPs can affect asthma by regulating enhancer activity, and our comprehensive bioinformatics and functional analysis could lead to biological insights into asthma pathogenesis.Graphic abstract.

Interactions between environmental pollutants and genetic susceptibility in asthma risk

Exposure to air pollution is associated with enhanced risk of developing asthma, notably in the presence of genetic risk factors. Interaction analyses have shown that both outdoor and indoor air pollution interact with genetic variability to increase the incidence of asthma. In this review, we summarize recent progress in candidate gene-based studies, as well as genome-wide gene-air pollution interaction studies. Advances in epigenetics have provided evidence for DNA methylation as a mediator in gene-air pollution interactions. Emerging strategies for study design and statistical analyses may improve power in future studies. Improved air pollution exposure assessment methods and asthma endo-typing can also be expected to increase the ability to detect biologically driven gene-air pollution interaction effects.

Optimising the identification of causal variants across varying genetic architectures in crops

Association studies use statistical links between genetic markers and the phenotype variation across many individuals to identify genes controlling variation in the target phenotype. However, this approach, particularly conducted on a genome-wide scale (GWAS), has limited power to identify the genes responsible for variation in traits controlled by complex genetic architectures. In this study, we employ real-world genotype datasets from four crop species with distinct minor allele frequency distributions, population structures and linkage disequilibrium patterns. We demonstrate that different GWAS statistical approaches provide favourable trade-offs between power and accuracy for traits controlled by different types of genetic architectures. FarmCPU provides the most favourable outcomes for moderately complex traits while a Bayesian approach adopted from genomic prediction provides the most favourable outcomes for extremely complex traits. We assert that by estimating the complexity of genetic architectures for target traits and selecting an appropriate statistical approach for the degree of complexity detected, researchers can substantially improve the ability to dissect the genetic factors controlling complex traits such as flowering time, plant height and yield component.

Integrated Functional Genomic Analysis Enables Annotation of Kidney Genome-Wide Association Study Loci

BACKGROUND

Linking genetic risk loci identified by genome-wide association studies (GWAS) to their causal genes remains a major challenge. Disease-associated genetic variants are concentrated in regions containing regulatory DNA elements, such as promoters and enhancers. Although researchers have previously published DNA maps of these regulatory regions for kidney tubule cells and glomerular endothelial cells, maps for podocytes and mesangial cells have not been available.

METHODS

We generated regulatory DNA maps (DNase-seq) and paired gene expression profiles (RNA-seq) from primary outgrowth cultures of human glomeruli that were composed mainly of podocytes and mesangial cells. We generated similar datasets from renal cortex cultures, to compare with those of the glomerular cultures. Because regulatory DNA elements can act on target genes across large genomic distances, we also generated a chromatin conformation map from freshly isolated human glomeruli.

RESULTS

We identified thousands of unique regulatory DNA elements, many located close to transcription factor genes, which the glomerular and cortex samples expressed at different levels. We found that genetic variants associated with kidney diseases (GWAS) and kidney expression quantitative trait loci were enriched in regulatory DNA regions. By combining GWAS, epigenomic, and chromatin conformation data, we functionally annotated 46 kidney disease genes.

CONCLUSIONS

We demonstrate a powerful approach to functionally connect kidney disease-/trait-associated loci to their target genes by leveraging unique regulatory DNA maps and integrated epigenomic and genetic analysis. This process can be applied to other kidney cell types and will enhance our understanding of genome regulation and its effects on gene expression in kidney disease.

Histone modifications and their role in epigenetics of atopy and allergic diseases

This review covers basic aspects of histone modification and the role of posttranslational histone modifications in the development of allergic diseases, including the immune mechanisms underlying this development. Together with DNA methylation, histone modifications (including histone acetylation, methylation, phosphorylation, ubiquitination, etc.) represent the classical epigenetic mechanisms. However, much less attention has been given to histone modifications than to DNA methylation in the context of allergy. A systematic review of the literature was undertaken to provide an unbiased and comprehensive update on the involvement of histone modifications in allergy and the mechanisms underlying this development. In addition to covering the growing interest in the contribution of histone modifications in regulating the development of allergic diseases, this review summarizes some of the evidence supporting this contribution. There are at least two levels at which the role of histone modifications is manifested. One is the regulation of cells that contribute to the allergic inflammation (T cells and macrophages) and those that participate in airway remodeling [(myo-) fibroblasts]. The other is the direct association between histone modifications and allergic phenotypes. Inhibitors of histone-modifying enzymes may potentially be used as anti-allergic drugs. Furthermore, epigenetic patterns may provide novel tools in the diagnosis of allergic disorders.

Adult Neural Stem Cells: Basic Research and Production Strategies for Neurorestorative Therapy

Over many decades, constructing genetically and phenotypically stable lines of neural stem cells (NSC) for clinical purposes with the aim of restoring irreversibly lost functions of nervous tissue has been one of the major goals for multiple research groups. The unique ability of stem cells to maintain their own pluripotent state even in the adult body has made them into the choice object of study. With the development of the technology for induced pluripotent stem cells (iPSCs) and direct transdifferentiation of somatic cells into the desired cell type, the initial research approaches based on the use of allogeneic NSCs from embryonic or fetal nervous tissue are gradually becoming a thing of the past. This review deals with basic molecular mechanisms for maintaining the pluripotent state of embryonic/induced stem and reprogrammed somatic cells, as well as with currently existing reprogramming strategies. The focus is on performing direct reprogramming while bypassing the stage of iPSCs which is known for genetic instability and an increased risk of tumorigenesis. A detailed description of various protocols for obtaining reprogrammed neural cells used in the therapy of the nervous system pathology is also provided.

Epigenetic regulation of T-helper cell differentiation, memory, and plasticity in allergic asthma

An estimated 300 million people currently suffer from asthma, which causes approximately 250 000 deaths a year. Allergen-specific T-helper (Th) cells produce cytokines that induce many of the hallmark features of asthma including airways hyperreactivity, eosinophilic and neutrophilic inflammation, mucus hypersecretion, and airway remodeling. Cytokine-producing Th subsets including Th1 (IFN-γ), Th2 (IL-4, IL-5, IL-13), Th9 (IL-9), Th17 (IL-17), Th22 (IL-22), and T regulatory (IL-10) cells have all been suggested to play a role in the development of asthma. Th differentiation involves genetic regulation of gene expression through the concerted action of cytokines, transcription factors, and epigenetic regulators. We describe how Th differentiation and plasticity is regulated by epigenetic histone and DNA modifications, with a focus on the regulation of histone methylation by members of the polycomb and trithorax complexes. In addition, we outline environmental influences that could influence epigenetic regulation of Th cells and discuss the potential to regulate Th plasticity and function through drugs targeting the epigenetic machinery. It is also becoming apparent that epigenetic regulation of allergen-specific memory Th cells may be important in the development and persistence of chronic allergies. Finally, we describe how epigenetic modifiers regulate cytokine memory in Th cells and describe recently identified hybrid, plastic, and pathogenic memory Th subsets the context of allergic asthma.

Pathogenic Th cell subsets in chronic inflammatory diseases

CD4(+) T cells play central roles to appropriate protection against pathogens. While, they can also be pathogenic driving inflammatory diseases. Besides the classical model of differentiation of T helper 1 (Th1) and Th2 cells, various CD4(+) T cell subsets, including Th17, Th9, T follicular helper (Tfh) and T regulatory (Treg) cells, have been recognized recently. In this review, we will focus on how these various CD4(+) T cell subsets contribute to the pathogenesis of immune-mediated inflammatory diseases. We will also discuss various unique subpopulations of T helper cells that have been identified. Recent advancement of the basic immunological research revealed that T helper cells are plastic than we imagined. So, we will focus on the molecular mechanisms underlying the generation of the plasticity and heterogeneity of T helper cell subsets. These latest finding regarding T helper cell subsets has pushed us to reconsider the etiology of immune-mediated inflammatory diseases beyond the model based on the conventional Th1/Th2 balance. Toward this end, we put forward another model, "the pathogenic Th population disease induction model", as a possible mechanism for the induction and/or persistence of immune-mediated inflammatory diseases.

Extensive Association of Common Disease Variants with Regulatory Sequence

Overlap between non-coding DNA regulatory sequences and common variant associations can help to identify specific cell and tissue types that are relevant for particular diseases. In a systematic manner, we analyzed variants from 94 genome-wide association studies (reporting at least 12 loci at p<5x10^-8) by projecting them onto 466 epigenetic datasets (characterizing DNase I hypersensitive sites; DHSs) derived from various adult and fetal tissue samples and cell lines including many biological replicates. We were able to confirm many expected associations, such as the involvement of specific immune cell types in immune-related diseases and tissue types in diseases that affect specific organs, for example, inflammatory bowel disease and coronary artery disease. Other notable associations include adrenal glands in coronary artery disease, the immune system in Alzheimer’s disease, and the kidney for bone marrow density. The association signals for some GWAS (for example, myopia or age at menarche) did not show a clear pattern with any of the cell or tissue types studied. In general, the identified variants from GWAS tend to be located outside coding regions. Altogether, we have performed an extensive characterization of GWAS signals in relation to cell and tissue-specific DHSs, demonstrating a key role for regulatory mechanisms in common diseases and complex traits.

17q21 asthma-risk variants switch CTCF binding and regulate IL-2 production by T cells

Asthma and autoimmune disease susceptibility has been strongly linked to genetic variants in the 17q21 haploblock that alter the expression of ORMDL3; however, the molecular mechanisms by which these variants perturb gene expression and the cell types in which this effect is most prominent are unclear. We found several 17q21 variants overlapped enhancers present mainly in primary immune cell types. CD4⁺ T cells showed the greatest increase (threefold) in ORMDL3 expression in individuals carrying the asthma-risk alleles, where ORMDL3 negatively regulated interleukin-2 production. The asthma-risk variants rs4065275 and rs12936231 switched CTCF-binding sites in the 17q21 locus, and 4C-Seq assays showed that several distal cis-regulatory elements upstream of the disrupted ZPBP2 CTCF-binding site interacted with the ORMDL3 promoter region in CD4⁺ T cells exclusively from subjects carrying asthma-risk alleles. Overall, our results suggested that T cells are one of the most prominent cell types affected by 17q21 variants.

Variations in the 17q21 locus are linked to asthma susceptibility and other autoimmune diseases. Here, the authors perform cell type-specific functional genomic analyses of asthma-risk SNPs, and show a genotype specific mechanism of differential gene regulation relevant to immune function.

CD4+ T-cell subsets in inflammatory diseases: beyond the Th1/Th2 paradigm

CD4(+)T cells are crucial for directing appropriate immune responses during host defense and for the pathogenesis of inflammatory diseases. In addition to the classical biphasic model of differentiation of T-helper 1 (Th1) and Th2 cells, unexpected increases in the numbers of CD4(+)T-cell subsets, including Th17, Th9, T follicular-helper (Tfh) and T-regulatory (Treg) cells, have been recognized. In the present review, we focus on how these various T-helper cell subsets contribute to the pathogenesis of immune-mediated inflammatory diseases. In particular, we focus on multiple sclerosis, psoriasis and asthma as typical model diseases in which multiple T-helper cell subsets have recently been suggested to play a role. We will also discuss various unique sub-populations of T-helper cells that have been identified. First, we will introduce the heterogeneous T-helper cell subsets, which are classified by their simultaneous expression of multiple key transcription factors. We will also introduce different kinds of memory-type Th2 cells, which are involved in the pathogenesis of chronic type-2 immune-related diseases. Finally, we will discuss the molecular mechanisms underlying the generation of the plasticity and heterogeneity of T-helper cell subsets. The latest progress in the study of T-helper cell subsets has forced us to reconsider the etiology of immune-mediated inflammatory diseases beyond the model based on the Th1/Th2 balance. To this end, we propose another model--the pathogenic T-helper population disease-induction model--as a possible mechanism for the induction and/or persistence of immune-mediated inflammatory diseases.

NRSF and BDNF polymorphisms as biomarkers of cognitive dysfunction in adults with newly diagnosed epilepsy

Cognitive dysfunction is a common comorbidity in people with epilepsy, but its causes remain unclear. It may be related to the etiology of the disorder, the consequences of seizures, or the effects of antiepileptic drug treatment. Genetics may also play a contributory role. We investigated the influence of variants in the genes encoding neuron-restrictive silencer factor (NRSF) and brain-derived neurotrophic factor (BDNF), proteins previously associated with cognition and epilepsy, on cognitive function in people with newly diagnosed epilepsy. A total of 82 patients who had previously undergone detailed neuropsychological assessment were genotyped for single nucleotide polymorphisms (SNPs) across the NRSF and BDNF genes. Putatively functional SNPs were included in a genetic association analysis with specific cognitive domains, including memory, psychomotor speed, and information processing. Cross-sectional and longitudinal designs were used to explore genetic influences on baseline cognition at diagnosis and change from baseline over the first year since diagnosis, respectively. We found a statistically significant association between genotypic variation and memory function at both baseline (NRSF: rs1105434, rs2227902 and BDNF: rs1491850, rs2030324, rs11030094) and in our longitudinal analysis (NRSF: rs2227902 and BDNF: rs12273363). Psychomotor speed was also associated with genotype (NRSF rs3796529) in the longitudinal assessment. In line with our previous work on general cognitive function in the healthy aging population, we observed an additive interaction between risk alleles for the NRSF rs2227902 (G) and BDNF rs6265 (A) polymorphisms which was again consistent with a significantly greater decline in delayed recall over the first year since diagnosis. These findings support a role for the NRSF-BDNF pathway in the modulation of cognitive function in patients with newly diagnosed epilepsy.

Decoding the non-coding genome: elucidating genetic risk outside the coding genome

Current evidence emerging from genome-wide association studies indicates that the genetic underpinnings of complex traits are likely attributable to genetic variation that changes gene expression, rather than (or in combination with) variation that changes protein-coding sequences. This is particularly compelling with respect to psychiatric disorders, as genetic changes in regulatory regions may result in differential transcriptional responses to developmental cues and environmental/psychosocial stressors. Until recently, however, the link between transcriptional regulation and psychiatric genetic risk has been understudied. Multiple obstacles have contributed to the paucity of research in this area, including challenges in identifying the positions of remote (distal from the promoter) regulatory elements (e.g. enhancers) and their target genes and the underrepresentation of neural cell types and brain tissues in epigenome projects - the availability of high-quality brain tissues for epigenetic and transcriptome profiling, particularly for the adolescent and developing brain, has been limited. Further challenges have arisen in the prediction and testing of the functional impact of DNA variation with respect to multiple aspects of transcriptional control, including regulatory-element interaction (e.g. between enhancers and promoters), transcription factor binding and DNA methylation. Further, the brain has uncommon DNA-methylation marks with unique genomic distributions not found in other tissues - current evidence suggests the involvement of non-CG methylation and 5-hydroxymethylation in neurodevelopmental processes but much remains unknown. We review here knowledge gaps as well as both technological and resource obstacles that will need to be overcome in order to elucidate the involvement of brain-relevant gene-regulatory variants in genetic risk for psychiatric disorders.

Epigenome-modifying tools in asthma

Asthma is a chronic disease which causes recurrent breathlessness affecting 300 million people worldwide of whom 250,000 die annually. The epigenome is a set of heritable modifications and tags that affect the genome without changing the intrinsic DNA sequence. These marks include DNA methylation, modifications to histone proteins around which DNA is wrapped and expression of noncoding RNA. Alterations in all of these processes have been reported in patients with asthma. In some cases these differences are linked to disease severity and susceptibility and may account for the limited value of genetic studies in asthma. Animal models of asthma suggest that epigenetic modifications and processes are linked to asthma and may be tractable targets for therapeutic intervention.

Disentangling the Effects of Colocalizing Genomic Annotations to Functionally Prioritize Non-coding Variants within Complex-Trait Loci

Identifying genomic annotations that differentiate causal from trait-associated variants is essential to fine mapping disease loci. Although many studies have identified non-coding functional annotations that overlap disease-associated variants, these annotations often colocalize, complicating the ability to use these annotations for fine mapping causal variation. We developed a statistical approach (Genomic Annotation Shifter [GoShifter]) to assess whether enriched annotations are able to prioritize causal variation. GoShifter defines the null distribution of an annotation overlapping an allele by locally shifting annotations; this approach is less sensitive to biases arising from local genomic structure than commonly used enrichment methods that depend on SNP matching. Local shifting also allows GoShifter to identify independent causal effects from colocalizing annotations. Using GoShifter, we confirmed that variants in expression quantitative trail loci drive gene-expression changes though DNase-I hypersensitive sites (DHSs) near transcription start sites and independently through 3' UTR regulation. We also showed that (1) 15%-36% of trait-associated loci map to DHSs independently of other annotations; (2) loci associated with breast cancer and rheumatoid arthritis harbor potentially causal variants near the summits of histone marks rather than full peak bodies; (3) variants associated with height are highly enriched in embryonic stem cell DHSs; and (4) we can effectively prioritize causal variation at specific loci.

Genetic-epigenetic dysregulation of thymic TSH receptor gene expression triggers thyroid autoimmunity

Graves disease (GD) is an autoimmune condition caused by interacting genetic and environmental factors. Genetic studies have mapped several single-nucleotide polymorphisms (SNPs) that are strongly associated with GD, but the mechanisms by which they trigger disease are unknown. We hypothesized that epigenetic modifications induced by microenvironmental influences of cytokines can reveal the functionality of GD-associated SNPs. We analyzed genome-wide histone H3 lysine 4 methylation and gene expression in thyroid cells induced by IFNα, a key cytokine secreted during viral infections, and overlapped them with known GD-associated SNPs. We mapped an open chromatin region overlapping two adjacent GD-associated SNPs (rs12101255 and rs12101261) in intron 1 of the thyroid stimulating hormone receptor (TSHR) gene. We then demonstrated that this region functions as a regulatory element through binding of the transcriptional repressor promyelocytic leukemia zinc finger protein (PLZF) at the rs12101261 site. Repression by PLZF depended on the rs12101261 disease susceptibility allele and was increased by IFNα. Intrathymic TSHR expression was decreased in individuals homozygous for the rs12101261 disease-associated genotype compared with carriers of the disease-protective allele. Our studies discovered a genetic-epigenetic interaction involving a noncoding SNP in the TSHR gene that regulates thymic TSHR gene expression and facilitates escape of TSHR-reactive T cells from central tolerance, triggering GD.

Epigenomic analysis of primary human T cells reveals enhancers associated with TH2 memory cell differentiation and asthma susceptibility

A characteristic feature of asthma is the aberrant accumulation, differentiation or function of memory CD4(+) T cells that produce type 2 cytokines (TH2 cells). By mapping genome-wide histone modification profiles for subsets of T cells isolated from peripheral blood of healthy and asthmatic individuals, we identified enhancers with known and potential roles in the normal differentiation of human TH1 cells and TH2 cells. We discovered disease-specific enhancers in T cells that differ between healthy and asthmatic individuals. Enhancers that gained the histone H3 Lys4 dimethyl (H3K4me2) mark during TH2 cell development showed the highest enrichment for asthma-associated single nucleotide polymorphisms (SNPs), which supported a pathogenic role for TH2 cells in asthma. In silico analysis of cell-specific enhancers revealed transcription factors, microRNAs and genes potentially linked to human TH2 cell differentiation. Our results establish the feasibility and utility of enhancer profiling in well-defined populations of specialized cell types involved in disease pathogenesis.

On the identification of potential regulatory variants within genome wide association candidate SNP sets

BACKGROUND

Genome wide association studies (GWAS) are a population-scale approach to the identification of segments of the genome in which genetic variations may contribute to disease risk. Current methods focus on the discovery of single nucleotide polymorphisms (SNPs) associated with disease traits. As there are many SNPs within identified risk loci, and the majority of these are situated within non-coding regions, a key challenge is to identify and prioritize variants affecting regulatory sequences that are likely to contribute to the phenotype assessed.

METHODS

We focused investigation on SNPs within lung and breast cancer GWAS loci that reached genome-wide significance for potential roles in gene regulation with a specific focus on SNPs likely to disrupt transcription factor binding sites. Within risk loci, the regulatory potential of sub-regions was classified using relevant open chromatin and epigenetic high throughput sequencing data sets from the ENCODE project in available cancer and normal cell lines. Furthermore, transcription factor affinity altering variants were predicted by comparison of position weight matrix scores between disease and reference alleles. Lastly, ChIP-seq data of transcription associated factors and topological domains were included as binding evidence and potential gene target inference.

RESULTS

The sets of SNPs, including both the disease-associated markers and those in high linkage disequilibrium with them, were significantly over-represented in regulatory sequences of cancer and/or normal cells; however, over-representation was generally not restricted to disease-relevant tissue specific regions. The calculated regulatory potential, allelic binding affinity scores and ChIP-seq binding evidence were the three criteria used to prioritize candidates. Fitting all three criteria, we highlighted breast cancer susceptibility SNPs and a borderline lung cancer relevant SNP located in cancer-specific enhancers overlapping multiple distinct transcription associated factor ChIP-seq binding sites.

CONCLUSION

Incorporating high throughput sequencing epigenetic and transcription factor data sets from both cancer and normal cells into cancer genetic studies reveals potential functional SNPs and informs subsequent characterization efforts.

Asthma genetics and personalised medicine

Unbiased genetic approaches, especially genome-wide association studies, have identified novel genetic targets in the pathogenesis of asthma, but so far these targets account for only a small proportion of the heritability of asthma. Recognition of the importance of disease heterogeneity, the need for improved disease phenotyping, and the fact that genes involved in the inception of asthma are likely to be different from those involved in severity widens the scope of asthma genetics. The identification of genes implicated in several causal pathways suggests that genetic scores could be used to capture the effect of genetic variations on individuals. Gene-environment interaction adds another layer of complexity, which is being successfully explored by epigenetic approaches. Pharmacogenetics is one example of how gene-environment interactions are already being taken into account in the identification of drug responders and non-responders, and patients most susceptible to adverse effects. Such applications represent one component of personalised medicine, an approach that places the individual at the centre of health care.

Genome-wide association studies in asthma; perhaps, the end of the beginning

PURPOSE OF REVIEW

A large number of genetic loci contribute towards an individual's susceptibility to asthma and other complex diseases. Genome-wide association studies (GWASs) have provided us with a wealth of loci associated with asthma susceptibility, asthma endotypes and responsiveness to asthma medications. The reproducibility of these genetic loci across different studies highlights the interplay of general and population-specific risk alleles in asthma. Although GWASs have been successful in identifying disease-associated loci, there is still large potential for such studies to provide further insights into asthma pathogenesis.

RECENT FINDINGS

GWASs over the past year have extended study design well beyond the simple case-control and continuous phenotype association formats, for example, including interactions with environmental factors, integrating GWAS data with epigenetic data and GWASs in animal models, incorporating pathway analyses and utilising emerging sequencing technologies.

SUMMARY

Moving beyond traditional GWAS formats is likely to significantly enhance our understanding of the genetic basis for asthma. This review discusses where we are after half a decade of asthma GWASs, and focuses on advances over the past year that show where the GWAS field is headed in the future.

Joint analysis of functional genomic data and genome-wide association studies of 18 human traits

Annotations of gene structures and regulatory elements can inform genome-wide association studies (GWASs). However, choosing the relevant annotations for interpreting an association study of a given trait remains challenging. I describe a statistical model that uses association statistics computed across the genome to identify classes of genomic elements that are enriched with or depleted of loci influencing a trait. The model naturally incorporates multiple types of annotations. I applied the model to GWASs of 18 human traits, including red blood cell traits, platelet traits, glucose levels, lipid levels, height, body mass index, and Crohn disease. For each trait, I used the model to evaluate the relevance of 450 different genomic annotations, including protein-coding genes, enhancers, and DNase-I hypersensitive sites in over 100 tissues and cell lines. The fraction of phenotype-associated SNPs influencing protein sequence ranged from around 2% (for platelet volume) up to around 20% (for low-density lipoprotein cholesterol), repressed chromatin was significantly depleted for SNPs associated with several traits, and cell-type-specific DNase-I hypersensitive sites were enriched with SNPs associated with several traits (for example, the spleen in platelet volume). Finally, reweighting each GWAS by using information from functional genomics increased the number of loci with high-confidence associations by around 5%.