Genome-wide eQTLs and heritability for gene expression traits in unrelated individuals

Background

While the possible sources underlying the so-called ‘missing heritability’ evident in current genome-wide association studies (GWAS) of complex traits have been actively pursued in recent years, resolving this mystery remains a challenging task. Studying heritability of genome-wide gene expression traits can shed light on the goal of understanding the relationship between phenotype and genotype. Here we used microarray gene expression measurements of lymphoblastoid cell lines and genome-wide SNP genotype data from 210 HapMap individuals to examine the heritability of gene expression traits.

Results

Heritability levels for expression of 10,720 genes were estimated by applying variance component model analyses and 1,043 expression quantitative loci (eQTLs) were detected. Our results indicate that gene expression traits display a bimodal distribution of heritability, one peak close to 0% and the other summit approaching 100%. Such a pattern of the within-population variability of gene expression heritability is common among different HapMap populations of unrelated individuals but different from that obtained in the CEU and YRI trio samples. Higher heritability levels are shown by housekeeping genes and genes associated with cis eQTLs. Both cis and trans eQTLs make comparable cumulative contributions to the heritability. Finally, we modelled gene-gene interactions (epistasis) for genes with multiple eQTLs and revealed that epistasis was not prevailing in all genes but made a substantial contribution in explaining total heritability for some genes analysed.

Conclusions

We utilised a mixed effect model analysis for estimating genetic components from population based samples. On basis of analyses of genome-wide gene expression from four HapMap populations, we demonstrated detailed exploitation of the distribution of genetic heritabilities for expression traits from different populations, and highlighted the importance of studying interaction at the gene expression level as an important source of variation underlying missing heritability.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-13) contains supplementary material, which is available to authorized users.     

LncRNA <Emphasis Type="Italic">ANRIL</Emphasis> Expression and <Emphasis Type="Italic">ANRIL</Emphasis> Gene Polymorphisms Contribute to the Risk of Ischemic Stroke in the Chinese Han Population

The aim of the present study was to explore the role of lncRNA ANRIL in the pathogenesis of ischemic stroke (IS) and coronary artery disease (CAD) and to determine the association between ANRIL variants and the genetic susceptibility of IS and CAD in the Chinese Han population. A genetic association study including 550 IS patients, 550 CAD patients, and 550 healthy controls was conducted. The expression levels of lncRNA ANRIL, CDKN2A, and CDKN2B were detected using qRT-PCR. Genotyping was performed by Sequenom MassARRAY on an Agena platform. Our study showed that IS patients had an increased lncRNA ANRIL expression (P?=?0.002) and a decreased CDKN2A expression (P?<?0.001) compared with normal controls. A significant difference with regard to the genotype distribution of rs2383207 was found between male IS patients and controls (P?=?0.011). The minor allele of rs2383207 significantly increased the IS risk under a recessive model (OR?=?1.52, 95% CI?=?1.05–2.21, P?=?0.027). The minor allele of rs1333049 was significantly associated with the risk of IS among the male patients under a recessive model (OR?=?1.56, 95% CI?=?1.04–2.35, P?=?0.031). However, no significant association was found between the ANRIL variants and the risk of CAD (all P?>?0.050). In addition, we found a decreased lncRNA ANRIL expression in IS patients who carried the GG genotype of rs1333049 compared with IS patients who carried the CC or CG genotype (P?=?0.041). In summary, we found that IS patients had an increased lncRNA ANRIL expression and a decreased CDKN2A expression compared with the controls, which might play an impellent role in pathological processes of IS. The ANRIL variants rs2383207 and rs1333049 were significantly associated with the risk of IS among males but not females in the Chinese Han population.     

Clustered gene expression changes flank targeted gene loci in knockout mice

Background

Gene expression profiling using microarrays is a powerful technology widely used to study regulatory networks. Profiling of mRNA levels in mutant organisms has the potential to identify genes regulated by the mutated protein.

Methodology/Principle Findings

Using tissues from multiple lines of knockout mice we have examined genome-wide changes in gene expression. We report that a significant proportion of changed genes were found near the targeted gene.

Conclusions/Significance

The apparent clustering of these genes was explained by the presence of flanking DNA from the parental ES cell. We provide recommendations for the analysis and reporting of microarray data from knockout mice     

Impact of the genome wide supported NRGN gene on anterior cingulate morphology in schizophrenia

Background

The rs12807809 single-nucleotide polymorphism in NRGN is a genetic risk variant with genome-wide significance for schizophrenia. The frequency of the T allele of rs12807809 is higher in individuals with schizophrenia than in those without the disorder. Reduced immunoreactivity of NRGN, which is expressed exclusively in the brain, has been observed in Brodmann areas (BA) 9 and 32 of the prefrontal cortex in postmortem brains from patients with schizophrenia compared with those in controls.

Methods

Genotype effects of rs12807809 were investigated on gray matter (GM) and white matter (WM) volumes using magnetic resonance imaging (MRI) with a voxel-based morphometry (VBM) technique in a sample of 99 Japanese patients with schizophrenia and 263 healthy controls.

Results

Although significant genotype-diagnosis interaction either on GM or WM volume was not observed, there was a trend of genotype-diagnosis interaction on GM volume in the left anterior cingulate cortex (ACC). Thus, the effects of NRGN genotype on GM volume of patients with schizophrenia and healthy controls were separately investigated. In patients with schizophrenia, carriers of the risk T allele had a smaller GM volume in the left ACC (BA32) than did carriers of the non-risk C allele. Significant genotype effect on other regions of the GM or WM was not observed for either the patients or controls.

Conclusions

Our findings suggest that the genome-wide associated genetic risk variant in the NRGN gene may be related to a small GM volume in the ACC in the left hemisphere in patients with schizophrenia.     

Genetic Variation in the Epidermal Transglutaminase Genes Is Not Associated with Atopic Dermatitis

Background

Atopic dermatitis (AD) is a common chronic inflammatory skin disorder where epidermal barrier dysfunction is a major factor in the pathogenesis. The identification of AD susceptibility genes related to barrier dysfunction is therefore of importance. The epidermal transglutaminases (TGM1, TGM3 and TGM5) encodes essential cross-linking enzymes in the epidermis.

Objective

To determine whether genetic variability in the epidermal transglutaminases contributes to AD susceptibility.

Methods

Forty-seven single nucleotide polymorphisms (SNPs) in the TGM1, TGM3 and TGM5 gene region were tested for genetic association with AD, independently and in relation to FLG genotype, using a pedigree disequilibrium test (PDT) in a Swedish material consisting of 1753 individuals from 539 families. In addition, a German case-control material, consisting of 533 AD cases and 1996 controls, was used for in silico analysis of the epidermal TGM regions. Gene expression of the TGM1, TGM3 and TGM5 gene was investigated by relative quantification with Real Time PCR (qRT-PCR). Immunohistochemical (IHC) analysis was performed to detect TG1, TG3 and TG5 protein expression in the skin of patients and healthy controls.

Results

PDT analysis identified a significant association between the TGM1 SNP rs941505 and AD with allergen-specific IgE in the Swedish AD family material. However, the association was not replicated in the German case-control material. No significant association was detected for analyzed SNPs in relation to FLG genotype. TG1, TG3 and TG5 protein expression was detected in AD skin and a significantly increased TGM3 mRNA expression was observed in lesional skin by qRT-PCR.

Conclusion

Although TGM1 and TGM3 may be differentially expressed in AD skin, the results from the genetic analysis suggest that genetic variation in the epidermal transglutaminases is not an important factor in AD susceptibility.     

The Effect of Single Nucleotide Polymorphisms from Genome Wide Association Studies in Multiple Sclerosis on Gene Expression

Background

Multiple sclerosis (MS) is a complex neurological disorder. Its aetiology involves both environmental and genetic factors. Recent genome-wide association studies have identified a number of single nucleotide polymorphisms (SNPs) associated with susceptibility to (MS). We investigated whether these genetic variations were associated with alteration in gene expression.

Methods/Principal Findings

We used a database of mRNA expression and genetic variation derived from immortalised peripheral lymphocytes to investigate polymorphisms associated with MS for correlation with gene expression. Several SNPs were found to be associated with changes in expression: in particular two with HLA-DQA1, HLA-DQA2, HLA-DQB1, HLA-DRB1, HLA-DRB4 and HLA-DRB5, one with ZFP57, one with CD58, two with IL7 and FAM164A, and one with FAM119B, TSFM and KUB3. We found minimal cross-over with a recent whole genome expression study in MS patients.

Discussion

We have shown that many susceptibility loci in MS are associated with changes in gene expression using an unbiased expression database. Several of these findings suggest novel gene candidates underlying the effects of MS-associated genetic variation.     

An efficient algorithm to explore liquid association on a genome-wide scale

Background

The growing wealth of public available gene expression data has made the systemic studies of how genes interact in a cell become more feasible. Liquid association (LA) describes the extent to which coexpression of two genes may vary based on the expression level of a third gene (the controller gene). However, genome-wide application has been difficult and resource-intensive. We propose a new screening algorithm for more efficient processing of LA estimation on a genome-wide scale and apply its use to a Saccharomyces cerevisiae data set.

Results

On a test subset of the data, the fast screening algorithm achieved >99.8% agreement with the exhaustive search of LA values, while reduced run time by 81–93 %. Using a well-known yeast cell-cycle data set with 6,178 genes, we identified triplet combinations with significantly large LA values. In an exploratory gene set enrichment analysis, the top terms for the controller genes in these triplets with large LA values are involved in some of the most fundamental processes in yeast such as energy regulation, transportation, and sporulation.

Conclusion

In summary, in this paper we propose a novel, efficient algorithm to explore LA on a genome-wide scale and identified triplets of interest in cell cycle pathways using the proposed method in a yeast data set. A software package named fastLiquidAssociation for implementing the algorithm is available through http://www.bioconductor.org.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0371-5) contains supplementary material, which is available to authorized users.     

Performance of Polygenic Scores for Predicting Phobic Anxiety

Context

Anxiety disorders are common, with a lifetime prevalence of 20% in the U.S., and are responsible for substantial burdens of disability, missed work days and health care utilization. To date, no causal genetic variants have been identified for anxiety, anxiety disorders, or related traits.

Objective

To investigate whether a phobic anxiety symptom score was associated with 3 alternative polygenic risk scores, derived from external genome-wide association studies of anxiety, an internally estimated agnostic polygenic score, or previously identified candidate genes.

Design

Longitudinal follow-up study. Using linear and logistic regression we investigated whether phobic anxiety was associated with polygenic risk scores derived from internal, leave-one out genome-wide association studies, from 31 candidate genes, and from out-of-sample genome-wide association weights previously shown to predict depression and anxiety in another cohort.

Setting and Participants

Study participants (n = 11,127) were individuals from the Nurses'' Health Study and Health Professionals Follow-up Study.

Main Outcome Measure

Anxiety symptoms were assessed via the 8-item phobic anxiety scale of the Crown Crisp Index at two time points, from which a continuous phenotype score was derived.

Results

We found no genome-wide significant associations with phobic anxiety. Phobic anxiety was also not associated with a polygenic risk score derived from the genome-wide association study beta weights using liberal p-value thresholds; with a previously published genome-wide polygenic score; or with a candidate gene risk score based on 31 genes previously hypothesized to predict anxiety.

Conclusion

There is a substantial gap between twin-study heritability estimates of anxiety disorders ranging between 20–40% and heritability explained by genome-wide association results. New approaches such as improved genome imputations, application of gene expression and biological pathways information, and incorporating social or environmental modifiers of genetic risks may be necessary to identify significant genetic predictors of anxiety.