Population genomics in a disease targeted primary cell model
Ontology highlight
ABSTRACT: The common genetic variants associated with complex traits typically lie in non-coding DNA and may alter gene regulation in a cell-type specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an eQTL study of primary human osteoblasts (HOb) derived from unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top BMD-associated SNPs were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment of converging cis-eQTLs as compared to LCLs. The top BMD loci with SNPs showing strong cis-effects on gene expression in HObs were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All variants were tested in the Swedish MrOs Cohort (n=3014), providing evidence for two novel BMD loci. These variants were then tested in the Rotterdam Study (n=2100), yielding converging evidence for BMD association at one locus. The cis-regulatory effect was further fine-mapped to the proximal promoter of the gene. Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.
Project description:<p>The common genetic variants associated with complex traits typically lie in non-coding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus. The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.</p> <p>Reprinted from Grundberg E, Kwan T, Ge B et al., Genome Res Nov 2009, with permission from Genome Research, Cold Spring Harbor Press, PMID:<a href="http://www.ncbi.nlm.nih.gov/pubmed/19654370" target="_blank">19654370</a>.</p>
Project description:<p>The common genetic variants associated with complex traits typically lie in non-coding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus. The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.</p> <p>Reprinted from Grundberg E, Kwan T, Ge B et al., Genome Res Nov 2009, with permission from Genome Research, Cold Spring Harbor Press, PMID:<a href="http://www.ncbi.nlm.nih.gov/pubmed/19654370" target="_blank">19654370</a>.</p>
Project description:We hypothesized that mapping genetic variants associated with promoter and enhancer functions can provide novel insights into the mechanism through which eQTLs (expression quantitative trait loci) exert their effects on gene expression. To this end, we quantified genome-wide promoter usage and enhancer activity and tested the resulting molecular phenotypes for association with nearby genetic variants to discover cis-QTLs. To perform such analyzes, we have produced deep transcriptome profiling of 154 unrelated central European individuals, applying deepCAGE (Cap Analysis of Gene Expression) to nuclear enriched total RNA extracted from EBV transformed lymphoblastoid cell lines (LCLs). The LCLs were either purchased from Coriell Cell Repository (CEU, n=86) or from the GenCord collection (n=68, Gutierrez-Arcelus M et al. Elife 2013).
Project description:Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long non-coding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor SNPs demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (FDR correction<0.05) with BMD through TWAS and had a colocalizing eQTL (regional colocalization probability (RCP)>0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis.
Project description:Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long non-coding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor SNPs demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (FDR correction<0.05) with BMD through TWAS and had a colocalizing eQTL (regional colocalization probability (RCP)>0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis.
Project description:Genome-wide association studies (GWAS) have identified 19 risk variants associated with colorectal cancer. As most of these risk variants reside outside the coding regions of genes, we conducted cis-expression quantitative trait loci (cis-eQTL) analyses to investigate possible regulatory functions on the expression of neighboring genes. Forty microsatellite stable and CpG island methylator phenotype-negative colorectal tumors and paired adjacent-normal colon tissues were used for genome-wide SNP and gene expression profiling in our cis-eQTL analyses. This submission represents transcriptome component of study.
Project description:Genome-wide association studies (GWAS) have identified 19 risk variants associated with colorectal cancer. As most of these risk variants reside outside the coding regions of genes, we conducted cis-expression quantitative trait loci (cis-eQTL) analyses to investigate possible regulatory functions on the expression of neighboring genes.
Project description:Genome-wide association studies (GWAS) have boosted our knowledge of genetic risk variants in autoimmune diseases (AIDs). Most of the risk variants are located within or near genes with immunological functions, and the majority is found to be non-coding, pointing towards a regulatory role. We have performed a cis expression quantitative trait locus (eQTL) screen to investigate whether single nucleotide polymorphisms (SNPs) associated with AIDs influence gene expression in thymus. Genotyping was performed using the Immunochip and 353 AID associated SNPs were tested against expression of surrounding genes (+/- 1 Mb) from human thymic tissue (N=42). We identified eight genes where the expression was associated with AID risk SNPs at a study-wide level of significance (P < 2.57x10-5). Five genes (FCRL3, RNASET2, C2orf74, SIRPG and SYS1) displayed cis eQTL signals also in other tissues, while for two loci (NPIPB8 and LOC388814), the eQTL signal appear to be thymus-specific. Since many AID risk variants from GWAS have been subsequently fine-mapped in recent Immunochip projects, we explored the overlap between these novel AID risk variants and the thymic eQTL regions. Moreover, we examined the functional annotation of the seven expression altering SNPs (eSNPs). Our study reveals autoimmune risk variants that act as eQTLs in thymus. We have highlighted functional variants within these genetic regions that potentially can represent causal autoimmune risk variants. Total RNA from 42 human thymic samples were obtained from children undergoing cardiac surgery.
Project description:The majority of variants associated with complex traits and common diseases identified by genome-wide association studies (GWAS) map to noncoding regions of the genome with unknown regulatory effects. By leveraging ancestrally diverse biobank-scale GWAS data, massively parallel CRISPR screens and single cell transcriptomic and proteomic sequencing, we discovered target genes of noncoding variants for blood trait loci. For 91 GWAS loci, we identified 124 target genes in cis, which were often — but not always — the closest genes to the fine-mapped variant. Using precise variant insertion via base editing, we connect specific variants with gene expression changes. We also identified trans-effect networks of noncoding loci when cis target genes encoded transcription factors or microRNAs, such as GFI1B and miR-142. Trans-regulatory networks were themselves enriched for fine-mapped GWAS variants, demonstrating polygenic contributions to complex traits. Co-expression clustering of GFI1B trans-target genes identifies gene networks specific to different blood cell fates and differentiation stages. This platform will enable massively parallel assays to characterize the target genes and mechanisms of human noncoding variants in both cis and trans.
Project description:Recent genome-wide association studies (GWAS) have identified a number of novel genetic associations with complex human diseases. In spite of these successes, results from GWAS generally explain only a small proportion of disease heritability, an observation termed the M-bM-^@M-^\missing heritability problem.M-bM-^@M-^] Several sources for the missing heritability have been proposed, including the contribution of many common variants with small individual effect sizes, which cannot be reliably found using the standard GWAS approach. The goal of our study was to explore a complementary approach, which combines GWAS results with functional data in order to identify novel genetic associations with small effect sizes. To do so, we conducted a GWAS for lymphocyte count, a physiologic quantitative trait associated with asthma, in 462 Hutterites. In parallel, we performed a genome-wide gene expression study in lymphoblastoid cell lines (LCLs) from 96 Hutterites. We found significant support for genetic associations using the GWAS data when we considered variants near the 193 genes whose expression levels across individuals were most correlated with lymphocyte counts. Interestingly, these variants are also enriched with signatures of an association with asthma susceptibility, an observation we were able to replicate. The associated loci include genes previously implicated in asthma susceptibility, as well as novel candidate genes enriched for functions related to T cell receptor signaling and ATP synthesis. Our results, therefore, establish a new set of asthma susceptibility candidate genes. More generally, our observations support the notion that many loci of small effects influence variation in lymphocyte count and asthma susceptibility. 96 RNA samples were collected (1 subsequently excluded) from lymphoblastoid cell lines derived from Hutterite subjects chosen to represent the extremes of absolute lymphocyte count. The high and low absolute lymphocyte count groups were balanced with respect to age, gender, and relatedness.