Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles

The rheumatic disease-associated FAM167A-BLK locus encodes DIORA-1, a novel disordered protein expressed highly in bronchial epithelium and alveolar macrophages

Triggering of autoimmunity that leads to rheumatic disease has been suggested to depend upon gene-environment interactions occurring in epithelial barriers and associated immune cells. Genetic studies have identified associations of the FAM167A-BLK locus with rheumatoid arthritis, systemic lupus erythematosus (SLE) and Sjögren's syndrome. While BLK (B lymphocyte kinase) has a well-established role in B cells, family with sequence similarity to 167 member A (FAM167A) and its gene family remain uncharacterized. To begin to understand the role of FAM167A in rheumatic disease pathogenesis, we explored this gene family and cloned and investigated the gene products. Expression of quantitative trait locus analysis was performed in immune cells. FAM167A and FAM167B were cloned from human peripheral blood mononuclear cells (PBMC). Gene conservation and protein properties were analysed by online tools, mRNA expression measured in mouse organs by quantitative polymerase chain reaction (qPCR) and protein expression investigated in human tissues by immunohistochemistry. We found that autoimmune risk genotypes within the FAM167A-BLK locus lead to increased expression of FAM167A. The FAM167 gene family includes two members, FAM167A and FAM167B, which are not homologous to any other annotated gene but are evolutionarily conserved. The encoded proteins, which we denote 'disordered autoimmunity' (DIORA)-1 and DIORA-2, respectively, are characterized by a high content of intrinsic disorder. Notably, DIORA-1 has its highest expression in the lung, detectable in both bronchial epithelium and alveolar macrophages with an endosomal localization pattern. In summary, the FAM167A gene is associated with several rheumatic diseases and encodes a novel disordered protein, DIORA-1, which is expressed highly in the lung, consistent with a potential role in disease pathogenesis.

What has GWAS done for HLA and disease associations?

The major histocompatibility complex (MHC) is located in chromosome 6p21 and contains crucial regulators of immune response, including human leucocyte antigen (HLA) genes, alongside other genes with nonimmunological roles. More recently, a repertoire of noncoding RNA genes, including expressed pseudogenes, has also been identified. The MHC is the most gene dense and most polymorphic part of the human genome. The region exhibits haplotype-specific linkage disequilibrium patterns, contains the strongest cis- and trans-eQTLs/meQTLs in the genome and is known as a hot spot for disease associations. Another layer of complexity is provided to the region by the extreme structural variation and copy number variations. While the HLA-B gene has the highest number of alleles, the HLA-DR/DQ subregion is structurally most variable and shows the highest number of disease associations. Reliance on a single reference sequence has complicated the design, execution and analysis of GWAS for the MHC region and not infrequently, the MHC region has even been excluded from the analysis of GWAS data. Here, we contrast features of the MHC region with the rest of the genome and highlight its complexities, including its functional polymorphisms beyond those determined by single nucleotide polymorphisms or single amino acid residues. One of the several issues with customary GWAS analysis is that it does not address this additional layer of polymorphisms unique to the MHC region. We highlight alternative approaches that may assist with the analysis of GWAS data from the MHC region and unravel associations with all functional polymorphisms beyond single SNPs. We suggest that despite already showing the highest number of disease associations, the true extent of the involvement of the MHC region in disease genetics may not have been uncovered.

Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition

Individual risk of type 2 diabetes (T2D) is modified by perturbations to the mass, distribution and function of adipose tissue. To investigate the mechanisms underlying these associations, we explored the molecular, cellular and whole-body effects of T2D-associated alleles near KLF14. We show that KLF14 diabetes-risk alleles act in adipose tissue to reduce KLF14 expression and modulate, in trans, the expression of 385 genes. We demonstrate, in human cellular studies, that reduced KLF14 expression increases pre-adipocyte proliferation but disrupts lipogenesis, and in mice, that adipose tissue-specific deletion of Klf14 partially recapitulates the human phenotype of insulin resistance, dyslipidemia and T2D. We show that carriers of the KLF14 T2D risk allele shift body fat from gynoid stores to abdominal stores and display a marked increase in adipocyte cell size, and that these effects on fat distribution, and the T2D association, are female specific. The metabolic risk associated with variation at this imprinted locus depends on the sex both of the subject and of the parent from whom the risk allele derives.

A Meta-Analysis of Genome-Wide Association Studies of Growth Differentiation Factor-15 Concentration in Blood

Blood levels of growth differentiation factor-15 (GDF-15), also known as macrophage inhibitory cytokine-1 (MIC-1), have been associated with various pathological processes and diseases, including cardiovascular disease and cancer. Prior studies suggest genetic factors play a role in regulating blood MIC-1/GDF-15 concentration. In the current study, we conducted the largest genome-wide association study (GWAS) to date using a sample of ∼5,400 community-based Caucasian participants, to determine the genetic variants associated with MIC-1/GDF-15 blood concentration. Conditional and joint (COJO), gene-based association, and gene-set enrichment analyses were also carried out to identify novel loci, genes, and pathways. Consistent with prior results, a locus on chromosome 19, which includes nine single nucleotide polymorphisms (SNPs) (top SNP, rs888663, p = 1.690 × 10^-35), was significantly associated with blood MIC-1/GDF-15 concentration, and explained 21.47% of its variance. COJO analysis showed evidence for two independent signals within this locus. Gene-based analysis confirmed the chromosome 19 locus association and in addition, a putative locus on chromosome 1. Gene-set enrichment analyses showed that the“COPI-mediated anterograde transport” gene-set was associated with MIC-1/GDF15 blood concentration with marginal significance after FDR correction (p = 0.067). In conclusion, a locus on chromosome 19 was associated with MIC-1/GDF-15 blood concentration with genome-wide significance, with evidence for a new locus (chromosome 1). Future studies using independent cohorts are needed to confirm the observed associations especially for the chromosomes 1 locus, and to further investigate and identify the causal SNPs that contribute to MIC-1/GDF-15 levels.

Risk of nontyphoidal Salmonella bacteraemia in African children is modified by STAT4

Nontyphoidal Salmonella (NTS) is a major cause of bacteraemia in Africa. The disease typically affects HIV-infected individuals and young children, causing substantial morbidity and mortality. Here we present a genome-wide association study (180 cases, 2677 controls) and replication analysis of NTS bacteraemia in Kenyan and Malawian children. We identify a locus in STAT4, rs13390936, associated with NTS bacteraemia. rs13390936 is a context-specific expression quantitative trait locus for STAT4 RNA expression, and individuals carrying the NTS-risk genotype demonstrate decreased interferon-γ (IFNγ) production in stimulated natural killer cells, and decreased circulating IFNγ concentrations during acute NTS bacteraemia. The NTS-risk allele at rs13390936 is associated with protection against a range of autoimmune diseases. These data implicate interleukin-12-dependent IFNγ-mediated immunity as a determinant of invasive NTS disease in African children, and highlight the shared genetic architecture of infectious and autoimmune disease.

RNAseq analysis of bronchial epithelial cells to identify COPD-associated genes and SNPs

Background

There is a need for more powerful methods to identify low-effect SNPs that contribute to hereditary COPD pathogenesis. We hypothesized that SNPs contributing to COPD risk through cis-regulatory effects are enriched in genes comprised by bronchial epithelial cell (BEC) expression patterns associated with COPD.

Methods

To test this hypothesis, normal BEC specimens were obtained by bronchoscopy from 60 subjects: 30 subjects with COPD defined by spirometry (FEV1/FVC < 0.7, FEV1% < 80%), and 30 non-COPD controls. Targeted next generation sequencing was used to measure total and allele-specific expression of 35 genes in genome maintenance (GM) genes pathways linked to COPD pathogenesis, including seven TP53 and CEBP transcription factor family members. Shrinkage linear discriminant analysis (SLDA) was used to identify COPD-classification models. COPD GWAS were queried for putative cis-regulatory SNPs in the targeted genes.

Results

On a network basis, TP53 and CEBP transcription factor pathway gene pair network connections, including key DNA repair gene ERCC5, were significantly different in COPD subjects (e.g., Wilcoxon rank sum test for closeness, p-value = 5.0E-11). ERCC5 SNP rs4150275 association with chronic bronchitis was identified in a set of Lung Health Study (LHS) COPD GWAS SNPs restricted to those in putative regulatory regions within the targeted genes, and this association was validated in the COPDgene non-hispanic white (NHW) GWAS. ERCC5 SNP rs4150275 is linked (D’ = 1) to ERCC5 SNP rs17655 which displayed differential allelic expression (DAE) in BEC and is an expression quantitative trait locus (eQTL) in lung tissue (p = 3.2E-7). SNPs in linkage (D’ = 1) with rs17655 were predicted to alter miRNA binding (rs873601). A classifier model that comprised gene features CAT, CEBPG, GPX1, KEAP1, TP73, and XPA had pooled 10-fold cross-validation receiver operator characteristic area under the curve of 75.4% (95% CI: 66.3%–89.3%). The prevalence of DAE was higher than expected (p = 0.0023) in the classifier genes.

Conclusions

GM genes comprised by COPD-associated BEC expression patterns were enriched for SNPs with cis-regulatory function, including a putative cis-rSNP in ERCC5 that was associated with COPD risk. These findings support additional total and allele-specific expression analysis of gene pathways with high prior likelihood for involvement in COPD pathogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12890-018-0603-y) contains supplementary material, which is available to authorized users.

CD4+ and B Lymphocyte Expression Quantitative Traits at Rheumatoid Arthritis Risk Loci in Patients With Untreated Early Arthritis: Implications for Causal Gene Identification

OBJECTIVE

Rheumatoid arthritis (RA) is a genetically complex disease of immune dysregulation. This study sought to gain further insight into the genetic risk mechanisms of RA by conducting an expression quantitative trait locus (eQTL) analysis of confirmed genetic risk loci in CD4+ T cells and B cells from carefully phenotyped patients with early arthritis who were naive to therapeutic immunomodulation.

METHODS

RNA and DNA were isolated from purified B and/or CD4+ T cells obtained from the peripheral blood of 344 patients with early arthritis. Genotyping and global gene expression measurements were carried out using Illumina BeadChip microarrays. Variants in linkage disequilibrium (LD) with non-HLA RA single-nucleotide polymorphisms (defined as r2 ≥ 0.8) were analyzed, seeking evidence of cis- or trans-eQTLs according to whether the associated probes were or were not within 4 Mb of these LD blocks.

RESULTS

Genes subject to cis-eQTL effects that were common to both CD4+ and B lymphocytes at RA risk loci were FADS1, FADS2, BLK, FCRL3, ORMDL3, PPIL3, and GSDMB. In contrast, those acting on METTL21B, JAZF1, IKZF3, and PADI4 were unique to CD4+ lymphocytes, with the latter candidate risk gene being identified for the first time in this cell subset. B lymphocyte-specific eQTLs for SYNGR1 and CD83 were also found. At the 8p23 BLK-FAM167A locus, adjacent genes were subject to eQTLs whose activity differed markedly between cell types; in particular, the FAM167A effect displayed striking B lymphocyte specificity. No trans-eQTLs approached experiment-wide significance, and linear modeling did not identify a significant influence of biologic covariates on cis-eQTL effect sizes.

CONCLUSION

These findings further refine the understanding of candidate causal genes in RA pathogenesis, thus providing an important platform from which downstream functional studies, directed toward particular cell types, may be prioritized.

Exome-wide association study identifies four novel loci for systemic lupus erythematosus in Han Chinese population

OBJECTIVES

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease of considerable genetic predisposition. Genome-wide association studies have identified tens of common variants for SLE. However, the majority of them reside in non-coding sequences. The contributions of coding variants have not yet been systematically evaluated.

METHODS

We performed a large-scale exome-wide study in 5004 SLE cases and 8179 healthy controls in a Han Chinese population using a custom exome array, and then genotyped 32 variants with suggestive evidence in an independent cohort of 13 246 samples. We further explored the regulatory effect of one novel non-coding single nucleotide polymorphism (SNP) in ex vivo experiments.

RESULTS

We discovered four novel SLE gene regions (LCT, TPCN2, AHNAK2 and TNFRSF13B) encompassing three novel missense variants (XP_016859577.1:p.Asn1639Ser, XP_016859577.1:p.Val219Phe and XP_005267356.1:p.Thr4664Ala) and two non-coding variants (rs10750836 and rs4792801) with genome-wide significance (pmeta <5.00×10-8). These variants are enriched in several chromatin states of primary B cells. The novel intergenic variant rs10750836 exhibited an expression quantitative trait locus effect on the TPCN2 gene in immune cells. Clones containing this novel SNP exhibited gene promoter activity for TPCN2 (P=1.38×10-3) whose expression level was reduced significantly in patients with SLE (P<2.53×10-2) and was suggested to be further modulated by rs10750836 in CD19+ B cells (P=7.57×10-5) in ex vivo experiments.

CONCLUSIONS

This study identified three novel coding variants and four new susceptibility gene regions for SLE. The results provide insights into the biological mechanism of SLE.

Population-based analysis of ocular Chlamydia trachomatis in trachoma-endemic West African communities identifies genomic markers of disease severity

BACKGROUND

Chlamydia trachomatis (Ct) is the most common infectious cause of blindness and bacterial sexually transmitted infection worldwide. Ct strain-specific differences in clinical trachoma suggest that genetic polymorphisms in Ct may contribute to the observed variability in severity of clinical disease.

METHODS

Using Ct whole genome sequences obtained directly from conjunctival swabs, we studied Ct genomic diversity and associations between Ct genetic polymorphisms with ocular localization and disease severity in a treatment-naïve trachoma-endemic population in Guinea-Bissau, West Africa.

RESULTS

All Ct sequences fall within the T2 ocular clade phylogenetically. This is consistent with the presence of the characteristic deletion in trpA resulting in a truncated non-functional protein and the ocular tyrosine repeat regions present in tarP associated with ocular tissue localization. We have identified 21 Ct non-synonymous single nucleotide polymorphisms (SNPs) associated with ocular localization, including SNPs within pmpD (odds ratio, OR = 4.07, p* = 0.001) and tarP (OR = 0.34, p* = 0.009). Eight synonymous SNPs associated with disease severity were found in yjfH (rlmB) (OR = 0.13, p* = 0.037), CTA0273 (OR = 0.12, p* = 0.027), trmD (OR = 0.12, p* = 0.032), CTA0744 (OR = 0.12, p* = 0.041), glgA (OR = 0.10, p* = 0.026), alaS (OR = 0.10, p* = 0.032), pmpE (OR = 0.08, p* = 0.001) and the intergenic region CTA0744-CTA0745 (OR = 0.13, p* = 0.043).

CONCLUSIONS

This study demonstrates the extent of genomic diversity within a naturally circulating population of ocular Ct and is the first to describe novel genomic associations with disease severity. These findings direct investigation of host-pathogen interactions that may be important in ocular Ct pathogenesis and disease transmission.

T cells are influenced by a long non-coding RNA in the autoimmune associated PTPN2 locus

Non-coding SNPs in the protein tyrosine phosphatase non-receptor type 2 (PTPN2) locus have been linked with several autoimmune diseases, including rheumatoid arthritis, type I diabetes, and inflammatory bowel disease. However, the functional consequences of these SNPs are poorly characterized. Herein, we show in blood cells that SNPs in the PTPN2 locus are highly correlated with DNA methylation levels at four CpG sites downstream of PTPN2 and expression levels of the long non-coding RNA (lncRNA) LINC01882 downstream of these CpG sites. We observed that LINC01882 is mainly expressed in T cells and that anti-CD3/CD28 activated naïve CD4⁺ T cells downregulate the expression of LINC01882. RNA sequencing analysis of LINC01882 knockdown in Jurkat T cells, using a combination of antisense oligonucleotides and RNA interference, revealed the upregulation of the transcription factor ZEB1 and kinase MAP2K4, both involved in IL-2 regulation. Overall, our data suggests the involvement of LINC01882 in T cell activation and hints towards an auxiliary role of these non-coding SNPs in autoimmunity associated with the PTPN2 locus.

Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis

Background

There is a pressing need in rheumatoid arthritis (RA) to identify patients who will not respond to first-line disease-modifying anti-rheumatic drugs (DMARD). We explored whether differences in genomic architecture represented by a chromosome conformation signature (CCS) in blood taken from early RA patients before methotrexate (MTX) treatment could assist in identifying non-response to DMARD and, whether there is an association between such a signature and RA specific expression quantitative trait loci (eQTL).

Methods

We looked for the presence of a CCS in blood from early RA patients commencing MTX using chromosome conformation capture by EpiSwitch™. Using blood samples from MTX responders, non-responders and healthy controls, a custom designed biomarker discovery array was refined to a 5-marker CCS that could discriminate between responders and non-responders to MTX. We cross-validated the predictive power of the CCS by generating 150 randomized groups of 59 early RA patients (30 responders and 29 non-responders) before MTX treatment. The CCS was validated using a blinded, independent cohort of 19 early RA patients (9 responders and 10 non-responders). Last, the loci of the CCS markers were mapped to RA-specific eQTL.

Results

We identified a 5-marker CCS that could identify, at baseline, responders and non-responders to MTX. The CCS consisted of binary chromosome conformations in the genomic regions of IFNAR1, IL-21R, IL-23, CXCL13 and IL-17A. When tested on a cohort of 59 RA patients, the CCS provided a negative predictive value of 90.0% for MTX response. When tested on a blinded independent validation cohort of 19 early RA patients, the signature demonstrated a true negative response rate of 86 and a 90% sensitivity for detection of non-responders to MTX. Only conformations in responders mapped to RA-specific eQTL.

Conclusions

Here we demonstrate that detection of a CCS in blood in early RA is able to predict inadequate response to MTX with a high degree of accuracy. Our results provide a proof of principle that a priori stratification of response to MTX is possible, offering a mechanism to provide alternative treatments for non-responders to MTX earlier in the course of the disease.

Electronic supplementary material

The online version of this article (10.1186/s12967-018-1387-9) contains supplementary material, which is available to authorized users.

Expression level of risk genes of MHC class II is a susceptibility factor for autoimmunity: New insights

To date, the study of the impact of major hystocompatibility complex on autoimmunity has been prevalently focused on structural diversity of MHC molecules in binding and presentation of (auto)antigens to cognate T cells. Recently, a number of experimental evidences suggested new points of view to investigate the complex relationships between MHC gene expression and the individual predisposition to autoimmune diseases. Irrespective of the nature of the antigen, a threshold of MHC-peptide complexes needs to be reached, as well as a threshold of T cell receptors engaged is required, for the activation and proliferation of autoantigen-reactive T cells. Moreover, it is well known that increased expression of MHC class II molecules may alter the T cell receptor repertoire during thymic development, and affect the survival and expansion of mature T cells. Many evidences confirmed that the level of both transcriptional and post-transcriptional regulation are involved in the modulation of the expression of MHC class II genes and that both contribute to the predisposition to autoimmune diseases. Here, we aim to focus some of these regulative aspects to better clarify the role of MHC class II genes in predisposition and development of autoimmunity.

Novel Transcriptional Activity and Extensive Allelic Imbalance in the Human MHC Region

The MHC region encodes HLA genes and is the most complex region in the human genome. The extensively polymorphic nature of the HLA hinders accurate localization and functional assessment of disease risk loci within this region. Using targeted capture sequencing and constructing individualized genomes for transcriptome alignment, we identified 908 novel transcripts within the human MHC region. These include 593 novel isoforms of known genes, 137 antisense strand RNAs, 119 novel long intergenic noncoding RNAs, and 5 transcripts of 3 novel putative protein-coding human endogenous retrovirus genes. We revealed allele-dependent expression imbalance involving 88% of all heterozygous transcribed single nucleotide polymorphisms throughout the MHC transcriptome. Among these variants, the genetic variant associated with Behçet's disease in the HLA-B/MICA region, which tags HLA-B*51, is within novel long intergenic noncoding RNA transcripts that are exclusively expressed from the haplotype with the protective but not the disease risk allele. Further, the transcriptome within the MHC region can be defined by 14 distinct coexpression clusters, with evidence of coregulation by unique transcription factors in at least 9 of these clusters. Our data suggest a very complex regulatory map of the human MHC, and can help uncover functional consequences of disease risk loci in this region.

Resolving TYK2 locus genotype-to-phenotype differences in autoimmunity

Thousands of genetic variants have been identified, which contribute to the development of complex diseases, but determining how to elucidate their biological consequences for translation into clinical benefit is challenging. Conflicting evidence regarding the functional impact of genetic variants in the tyrosine kinase 2 (TYK2) gene, which is differentially associated with common autoimmune diseases, currently obscures the potential of TYK2 as a therapeutic target. We aimed to resolve this conflict by performing genetic meta-analysis across disorders; subsequent molecular, cellular, in vivo, and structural functional follow-up; and epidemiological studies. Our data revealed a protective homozygous effect that defined a signaling optimum between autoimmunity and immunodeficiency and identified TYK2 as a potential drug target for certain common autoimmune disorders.

Deciphering the genetic regulation of peripheral blood transcriptome in pigs through expression genome-wide association study and allele-specific expression analysis

BACKGROUND

Efforts to improve sustainability in livestock production systems have focused on two objectives: investigating the genetic control of immune function as it pertains to robustness and disease resistance, and finding predictive markers for use in breeding programs. In this context, the peripheral blood transcriptome represents an important source of biological information about an individual's health and immunological status, and has been proposed for use as an intermediate phenotype to measure immune capacity. The objective of this work was to study the genetic architecture of variation in gene expression in the blood of healthy young pigs using two approaches: an expression genome-wide association study (eGWAS) and allele-specific expression (ASE) analysis.

RESULTS

The blood transcriptomes of 60-day-old Large White pigs were analyzed by expression microarrays for eGWAS (242 animals) and by RNA-Seq for ASE analysis (38 animals). Using eGWAS, the expression levels of 1901 genes were found to be associated with expression quantitative trait loci (eQTLs). We recovered 2839 local and 1752 distant associations (Single Nucleotide Polymorphism or SNP located less or more than 1 Mb from expression probe, respectively). ASE analyses confirmed the extensive cis-regulation of gene transcription in blood, and revealed allelic imbalance in 2286 SNPs, which affected 763 genes. eQTLs and ASE-genes were widely distributed on all chromosomes. By analyzing mutually overlapping eGWAS results, we were able to describe putative regulatory networks, which were further refined using ASE data. At the functional level, genes with genetically controlled expression that were detected by eGWAS and/or ASE analyses were significantly enriched in biological processes related to RNA processing and immune function. Indeed, numerous distant and local regulatory relationships were detected within the major histocompatibility complex region on chromosome 7, revealing ASE for most class I and II genes.

CONCLUSIONS

This study represents, to the best of our knowledge, the first genome-wide map of the genetic control of gene expression in porcine peripheral blood. These results represent an interesting resource for the identification of genetic markers and blood biomarkers associated with variations in immunity traits in pigs, as well as any other complex traits for which blood is an appropriate surrogate tissue.

Type 1 diabetes genome-wide association studies: not to be lost in translation

Genetic studies have identified >60 loci associated with the risk of developing type 1 diabetes (T1D). The vast majority of these are identified by genome-wide association studies (GWAS) using large case-control cohorts of European ancestry. More than 80% of the heritability of T1D can be explained by GWAS data in this population group. However, with few exceptions, their individual contribution to T1D risk is low and understanding their function in disease biology remains a huge challenge. GWAS on its own does not inform us in detail on disease mechanisms, but the combination of GWAS data with other omics-data is beginning to advance our understanding of T1D etiology and pathogenesis. Current knowledge supports the notion that genetic variation in both pancreatic β cells and in immune cells is central in mediating T1D risk. Advances, perspectives and limitations of GWAS are discussed in this review.

cis-Acting Complex-Trait-Associated lincRNA Expression Correlates with Modulation of Chromosomal Architecture

Intergenic long noncoding RNAs (lincRNAs) are the largest class of transcripts in the human genome. Although many have recently been linked to complex human traits, the underlying mechanisms for most of these transcripts remain undetermined. We investigated the regulatory roles of a high-confidence and reproducible set of 69 trait-relevant lincRNAs (TR-lincRNAs) in human lymphoblastoid cells whose biological relevance is supported by their evolutionary conservation during recent human history and genetic interactions with other trait-associated loci. Their enrichment in enhancer-like chromatin signatures, interactions with nearby trait-relevant protein-coding loci, and preferential location at topologically associated domain (TAD) boundaries provide evidence that TR-lincRNAs likely regulate proximal trait-relevant gene expression in cis by modulating local chromosomal architecture. This is consistent with the positive and significant correlation found between TR-lincRNA abundance and intra-TAD DNA-DNA contacts. Our results provide insights into the molecular mode of action by which TR-lincRNAs contribute to complex human traits.

Shared genetic origin of asthma, hay fever and eczema elucidates allergic disease biology

Asthma, hay fever (or allergic rhinitis) and eczema (or atopic dermatitis) often coexist in the same individuals1, partly because of a shared genetic origin2–4. To identify shared risk variants, we performed a genome-wide association study (GWAS, n=360,838) of a broad allergic disease phenotype that considers the presence of any one of these three diseases. We identified 136 independent risk variants (P<3x10^-8), including 73 not previously reported, which implicate 132 nearby genes in allergic disease pathophysiology. Disease-specific effects were detected for only six variants, confirming that most represent shared risk factors. Tissue-specific heritability and biological process enrichment analyses suggest that shared risk variants influence lymphocyte-mediated immunity. Six target genes provide an opportunity for drug repositioning, while for 36 genes CpG methylation was found to influence transcription independently of genetic effects. Asthma, hay fever and eczema partly coexist because they share many genetic risk variants that dysregulate the expression of immune-related genes.

Sex influences eQTL effects of SLE and Sjögren's syndrome-associated genetic polymorphisms

BACKGROUND

Systemic lupus erythematosus (SLE) and primary Sjögren's syndrome (pSS) are autoimmune disorders characterized by autoantibodies, dysregulated B cells, and notably high female-to-male incidence ratios. Genome-wide association studies have identified several susceptibility SNPs for both diseases. Many SNPs in the genome are expression quantitative trait loci (eQTLs), with context-dependent effects. Assuming that sex is a biological context, we investigated whether SLE/pSS SNPs act as eQTLs in B cells and used a disease-targeted approach to understand if they display sex-specific effects.

METHODS

We used genome-wide genotype and gene expression data from primary B cells from 125 males and 162 females. The MatrixEQTL R package was used to identify eQTLs within a genomic window of 2 Mb centered on each of 22 established SLE and/or pSS susceptibility SNPs. To find sex-specific eQTLs, we used a linear model with a SNP * sex interaction term.

RESULTS

We found ten SNPs affecting the expression of 16 different genes (FDR < 0.05). rs7574865-INPP1, rs7574865-MYO1B, rs4938573-CD3D, rs11755393-SNRPC, and rs4963128-PHRF1 were novel observations for the immune compartment and B cells. By analyzing the SNP * sex interaction terms, we identified six genes with differentially regulated expression in females compared to males, depending on the genotype of SLE/pSS-associated SNPs: SLC39A8 (BANK1 locus), CD74 (TNIP1 locus), PXK, CTSB (BLK/FAM167A locus), ARCN1 (CXCR5 locus), and DHX9 (NCF2 locus).

CONCLUSIONS

We identified several unknown sex-specific eQTL effects of SLE/pSS-associated genetic polymorphisms and provide novel insight into how gene-sex interactions may contribute to the sex bias in systemic autoimmune diseases.

Interplay of cis and trans mechanisms driving transcription factor binding and gene expression evolution

Noncoding regulatory variants play a central role in the genetics of human diseases and in evolution. Here we measure allele-specific transcription factor binding occupancy of three liver-specific transcription factors between crosses of two inbred mouse strains to elucidate the regulatory mechanisms underlying transcription factor binding variations in mammals. Our results highlight the pre-eminence of cis-acting variants on transcription factor occupancy divergence. Transcription factor binding differences linked to cis-acting variants generally exhibit additive inheritance, while those linked to trans-acting variants are most often dominantly inherited. Cis-acting variants lead to local coordination of transcription factor occupancies that decay with distance; distal coordination is also observed and may be modulated by long-range chromatin contacts. Our results reveal the regulatory mechanisms that interplay to drive transcription factor occupancy, chromatin state, and gene expression in complex mammalian cell states.

Consensus Genome-Wide Expression Quantitative Trait Loci and Their Relationship with Human Complex Trait Disease

Most of the risk loci identified from genome-wide association (GWA) studies do not provide direct information on the biological basis of a disease or on the underlying mechanisms. Recent expression quantitative trait locus (eQTL) association studies have provided information on genetic factors associated with gene expression variation. These eQTLs might contribute to phenotype diversity and disease susceptibility, but interpretation is handicapped by low reproducibility of the expression results. To address this issue, we have generated a set of consensus eQTLs by integrating publicly available data for specific human populations and cell types. Overall, we find over 4000 genes that are involved in high-confidence eQTL relationships. To elucidate the role that eQTLs play in human common diseases, we matched the high-confidence eQTLs to a set of 335 disease risk loci identified from the Wellcome Trust Case Control Consortium GWA study and follow-up studies for 7 human complex trait diseases-bipolar disorder (BD), coronary artery disease (CAD), Crohn's disease (CD), hypertension (HT), rheumatoid arthritis (RA), type 1 diabetes (T1D), and type 2 diabetes (T2D). The results show that the data are consistent with ∼50% of these disease loci arising from an underlying expression change mechanism.

Context-specific effects of genetic variants associated with autoimmune disease

Autoimmune diseases such as rheumatoid arthritis and coeliac disease are typical examples of complex genetic diseases caused by a combination of genetic and non-genetic risk factors. Insight into the genetic risk factors (single nucleotide polymorphisms (SNPs)) has increased since genome-wide association studies (GWAS) became possible in 2007 and, for individual diseases, SNPs can now explain some 15-50% of genetic risk. GWAS have also shown that some 50% of the genetic risk factors for individual autoimmune diseases overlap between different diseases. Thus, shared risk factors may converge to pathways that, when perturbed by genetic variation, predispose to autoimmunity in general. This raises the question of what determines disease specificity, and suggests that identical risk factors may have different effects in various autoimmune diseases. Addressing this question requires translation of genetic risk factors to causal genes and then to molecular and cellular pathways. Since >90% of the genetic risk factors are found in the non-coding part of the genome (i.e. outside the exons of protein-coding genes) and can have an impact on gene regulation, there is an urgent need to better understand the non-coding part of the genome. Here, we will outline the methods being used to unravel the gene regulatory networks perturbed in autoimmune diseases and the importance of doing this in the relevant cell types. We will highlight findings in coeliac disease, which manifests in the small intestine, to demonstrate how cell type and disease context can impact on the consequences of genetic risk factors.

Identification of expression quantitative trait loci (eQTLs) in human peripheral blood mononuclear cells (PBMCs) and shared with liver and brain

PBMCs are essential for immunity and involved in various diseases. To identify genetic variations contributing to PBMCs transcriptome-wide gene expression, we performed a genome-wide eQTL analysis by using genome-wide SNPs data and transcriptome-wide mRNA expression data. To assess whether there are common regulation patterns shared among different tissues/organs, public datasets were utilized to identify common eQTLs shared with PBMCs in lymphoblastoid, monocytes, liver, and brain. Allelic expression imbalance (AEI) assay was employed to validate representative eQTLs identified. We identified 443 cis- and 2386 trans-eSNPs (FDR <0.05), which regulated 128 and 635 target genes, respectively. A transcriptome-wide expression regulation network was constructed, highlighting the importance of 28 pleiotropic eSNPs and 18 dually (cis- and trans-) regulated genes. Three genes, that is, TIPRL, HSPB8, and EGLN3, were commonly regulated by hundreds of eSNPs and constituted a very complex interaction network. Strikingly, the missense SNP rs371513 trans- regulated 25 target genes, which were functionally related to poly(A) RNA binding. Among 8904 eQTLs (P < 0.001) identified herein in PBMCs, a minority (163) was overlapped with lymphoblastoid, monocytes, liver, and/or brain. Besides, two cis-eSNPs in PBMC were confirmed by AEI. The present results demonstrated a comprehensive expression regulation network for human PBMCs and may provide novel insights into the pathogenesis of immunological diseases related to PBMCs.

Genetic correlations reveal the shared genetic architecture of transcription in human peripheral blood

Transcript co-expression is regulated by a combination of shared genetic and environmental factors. Here, we estimate the proportion of co-expression that is due to shared genetic variance. To do so, we estimated the genetic correlations between each pairwise combination of 2469 transcripts that are highly heritable and expressed in whole blood in 1748 unrelated individuals of European ancestry. We identify 556 pairs with a significant genetic correlation of which 77% are located on different chromosomes, and report 934 expression quantitative trait loci, identified in an independent cohort, with significant effects on both transcripts in a genetically correlated pair. We show significant enrichment for transcription factor control and physical proximity through chromatin interactions as possible mechanisms of shared genetic control. Finally, we construct networks of interconnected transcripts and identify their underlying biological functions. Using genetic correlations to investigate transcriptional co-regulation provides valuable insight into the nature of the underlying genetic architecture of gene regulation.

Covariance of gene expression pairs is due to a combination of shared genetic and environmental factors. Here the authors estimate the genetic correlation between highly heritable pairs and identify transcription factor control and chromatin interactions as possible mechanisms of correlation.

Unique Allelic eQTL Clusters in Human MHC Haplotypes

The control of gene regulation within the major histocompatibility complex (MHC) remains poorly understood, despite several expression quantitative trait loci (eQTL) studies revealing an association of MHC gene expression with independent tag-single nucleotide polymorphisms (SNPs). MHC haplotype variation may exert a greater effect on gene expression phenotype than specific single variants. To explore the effect of MHC haplotype sequence diversity on gene expression phenotypes across the MHC, we examined the MHC transcriptomic landscape at haplotype-specific resolution for three prominent MHC haplotypes (A2-B46-DR9, A33-B58-DR3, and A1-B8-DR3) derived from MHC-homozygous B-lymphoblastoid cell lines (B-LCLs). We demonstrate that MHC-wide gene expression patterns are dictated by underlying haplotypes, and identify 36 differentially expressed genes. By mapping these haplotype sequence variations to known eQTL, we provide evidence that unique allelic combinations of eQTL, embedded within haplotypes, are correlated with the level of expression of 17 genes. Interestingly, the influence of haplotype sequence on gene expression is not homogenous across the MHC. We show that haplotype sequence polymorphisms within or proximate to HLA-A, HLA-C, C4A, and HLA-DRB regions exert haplotype-specific gene regulatory effects, whereas the expression of genes in other parts of the MHC region are not affected by the haplotype sequence. Overall, we demonstrate that MHC haplotype sequence diversity can impact phenotypic outcome via the alteration of transcriptional variability, indicating that a haplotype-based approach is fundamental for the assessment of trait associations in the MHC.

Effects of Type 1 Diabetes Risk Alleles on Immune Cell Gene Expression

Genetic studies have identified 61 variants associated with the risk of developing Type 1 Diabetes (T1D). The functions of most of the non-HLA (Human Leukocyte Antigen) genetic variants remain unknown. We found that only 16 of these risk variants could potentially be linked to a protein-coding change. Therefore, we investigated whether these variants affected susceptibility by regulating changes in gene expression. To do so, we examined whole transcriptome profiles of 600 samples from the Type 1 Diabetes Genetics Consortium (T1DGC). These comprised four different immune cell types (Epstein-Barr virus (EBV)-transformed B cells, either basal or after stimulation; and cluster of differentiation (CD)4+ and CD8+ T cells). Many of the T1D-associated risk variants regulated expression of either neighboring (cis-) or distant (trans-) genes. In brief, 24 of the non-HLA T1D variants affected the expression of 31 nearby genes (cis) while 25 affected 38 distant genes (trans). The effects were highly significant (False Discovery Rate p < 0.001). In addition, we searched in public databases for expression effects of T1D single nucleotide polymorphisms (SNPs) in other immune cell types such as CD14+ monocytes, lipopolysaccharide (LPS) stimulated monocytes, and CD19+ B cells. In this paper, we review the (expression quantitative trait loci (eQTLs) associated with each of the 60 T1D variants and provide a summary of the genes impacted by T1D risk alleles in various immune cells. We then review the methodological steps involved in analyzing the function of genome wide association studies (GWAS)-identified variants, with emphasis on those affecting gene expression. We also discuss recent advancements in the methodologies and their advantages. We conclude by suggesting future study designs that will aid in the study of T1D risk variants.

Methods to analyze big data in pharmacogenomics research

The scale and scope of pharmacogenomics research continues to expand as the cost and efficiency of molecular data generation techniques advance. These new technologies give rise to enormous opportunity for the identification of important genetic and genomic factors important for drug treatment response. With this opportunity come significant challenges. Most of these can be categorized as 'big data' issues, facing not only pharmacogenomics, but other fields in the life sciences as well. In this review, we describe some of the analysis techniques and tools being implemented for genetic/genomic discovery in pharmacogenomics.

The Role of Osteopontin and Its Gene on Glucocorticoid Response in Myasthenia Gravis

Biomarkers that assess treatment response for patients with the autoimmune disorder, myasthenia gravis (MG), have not been evaluated to a significant extent. We hypothesized the pro-inflammatory cytokine, osteopontin (OPN), may be associated with variability of response to glucocorticoids (GCs) in patients with MG. A cohort of 250 MG patients treated with standardized protocol of GCs was recruited, and plasma OPN and polymorphisms of its gene, secreted phosphoprotein 1 (SPP1), were evaluated. Mean OPN levels were higher in patients compared to healthy controls. Carriers of rs11728697*T allele (allele definition: one of two or more alternative forms of a gene) were more frequent in the poorly GC responsive group compared to the GC responsive group indicating an association of rs11728697*T allele with GC non-responsiveness. One risk haplotype (AGTACT) was identified associated with GC non-responsiveness compared with GC responsive MG group. Genetic variations of SPP1 were found associated with the response to GC among MG patients.

Polygenic burdens on cell-specific pathways underlie the risk of rheumatoid arthritis

Recent evidence suggests that a substantial portion of complex disease risk alleles modify gene expression in a cell-specific manner. To identify candidate causal genes and biological pathways of immune-related complex diseases, we conducted expression quantitative trait loci (eQTL) analysis on five subsets of immune cells (CD4⁺ T cells, CD8⁺ T cells, B cells, natural killer (NK) cells and monocytes) and unfractionated peripheral blood from 105 healthy Japanese volunteers. We developed a three-step analytical pipeline comprising (i) prediction of individual gene expression using our eQTL database and public epigenomic data, (ii) gene-level association analysis and (iii) prediction of cell-specific pathway activity by integrating the direction of eQTL effects. By applying this pipeline to rheumatoid arthritis data sets, we identified candidate causal genes and a cytokine pathway (upregulation of tumor necrosis factor (TNF) in CD4⁺ T cells). Our approach is an efficient way to characterize the polygenic contributions and potential biological mechanisms of complex diseases.

Mapping eQTLs with RNA-seq reveals novel susceptibility genes, non-coding RNAs and alternative-splicing events in systemic lupus erythematosus

Studies attempting to functionally interpret complex-disease susceptibility loci by GWAS and eQTL integration have predominantly employed microarrays to quantify gene-expression. RNA-Seq has the potential to discover a more comprehensive set of eQTLs and illuminate the underlying molecular consequence. We examine the functional outcome of 39 variants associated with Systemic Lupus Erythematosus (SLE) through the integration of GWAS and eQTL data from the TwinsUK microarray and RNA-Seq cohort in lymphoblastoid cell lines. We use conditional analysis and a Bayesian colocalisation method to provide evidence of a shared causal-variant, then compare the ability of each quantification type to detect disease relevant eQTLs and eGenes. We discovered the greatest frequency of candidate-causal eQTLs using exon-level RNA-Seq, and identified novel SLE susceptibility genes (e.g. NADSYN1 and TCF7) that were concealed using microarrays, including four non-coding RNAs. Many of these eQTLs were found to influence the expression of several genes, supporting the notion that risk haplotypes may harbour multiple functional effects. Novel SLE associated splicing events were identified in the T-reg restricted transcription factor, IKZF2, and other candidate genes (e.g. WDFY4) through asQTL mapping using the Geuvadis cohort. We have significantly increased our understanding of the genetic control of gene-expression in SLE by maximising the leverage of RNA-Seq and performing integrative GWAS-eQTL analysis against gene, exon, and splice-junction quantifications. We conclude that to better understand the true functional consequence of regulatory variants, quantification by RNA-Seq should be performed at the exon-level as a minimum, and run in parallel with gene and splice-junction level quantification.

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.

Genetic regulation of gene expression in the epileptic human hippocampus

Epilepsy is a serious and common neurological disorder. Expression quantitative loci (eQTL) analysis is a vital aid for the identification and interpretation of disease-risk loci. Many eQTLs operate in a tissue- and condition-specific manner. We have performed the first genome-wide cis-eQTL analysis of human hippocampal tissue to include not only normal (n = 22) but also epileptic (n = 22) samples. We demonstrate that disease-associated variants from an epilepsy GWAS meta-analysis and a febrile seizures (FS) GWAS are significantly more enriched with epilepsy-eQTLs than with normal hippocampal eQTLs from two larger independent published studies. In contrast, GWAS meta-analyses of two other brain diseases associated with hippocampal pathology (Alzheimer's disease and schizophrenia) are more enriched with normal hippocampal eQTLs than with epilepsy-eQTLs. These observations suggest that an eQTL analysis that includes disease-affected brain tissue is advantageous for detecting additional risk SNPs for the afflicting and closely related disorders, but not for distinct diseases affecting the same brain regions. We also show that epilepsy eQTLs are enriched within epilepsy-causing genes: an epilepsy cis-gene is significantly more likely to be a causal gene for a Mendelian epilepsy syndrome than to be a causal gene for another Mendelian disorder. Epilepsy cis-genes, compared to normal hippocampal cis-genes, are more enriched within epilepsy-causing genes. Hence, we utilize the epilepsy eQTL data for the functional interpretation of epilepsy disease-risk variants and, thereby, highlight novel potential causal genes for sporadic epilepsy. In conclusion, an epilepsy-eQTL analysis is superior to normal hippocampal tissue eQTL analyses for identifying the variants and genes underlying epilepsy.

Understanding Human Autoimmunity and Autoinflammation Through Transcriptomics

Transcriptomics, the high-throughput characterization of RNAs, has been instrumental in defining pathogenic signatures in human autoimmunity and autoinflammation. It enabled the identification of new therapeutic targets in IFN-, IL-1- and IL-17-mediated diseases. Applied to immunomonitoring, transcriptomics is starting to unravel diagnostic and prognostic signatures that stratify patients, track molecular changes associated with disease activity, define personalized treatment strategies, and generally inform clinical practice. Herein, we review the use of transcriptomics to define mechanistic, diagnostic, and predictive signatures in human autoimmunity and autoinflammation. We discuss some of the analytical approaches applied to extract biological knowledge from high-dimensional data sets. Finally, we touch upon emerging applications of transcriptomics to study eQTLs, B and T cell repertoire diversity, and isoform usage.

Network-based analysis reveals novel gene signatures in peripheral blood of patients with chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is currently the third leading cause of death and there is a huge unmet clinical need to identify disease biomarkers in peripheral blood. Compared to gene level differential expression approaches to identify gene signatures, network analyses provide a biologically intuitive approach which leverages the co-expression patterns in the transcriptome to identify modules of co-expressed genes.

METHODS

A weighted gene co-expression network analysis (WGCNA) was applied to peripheral blood transcriptome from 238 COPD subjects to discover co-expressed gene modules. We then determined the relationship between these modules and forced expiratory volume in 1 s (FEV1). In a second, independent cohort of 381 subjects, we determined the preservation of these modules and their relationship with FEV1. For those modules that were significantly related to FEV1, we determined the biological processes as well as the blood cell-specific gene expression that were over-represented using additional external datasets.

RESULTS

Using WGCNA, we identified 17 modules of co-expressed genes in the discovery cohort. Three of these modules were significantly correlated with FEV1 (FDR < 0.1). In the replication cohort, these modules were highly preserved and their FEV1 associations were reproducible (P < 0.05). Two of the three modules were negatively related to FEV1 and were enriched in IL8 and IL10 pathways and correlated with neutrophil-specific gene expression. The positively related module, on the other hand, was enriched in DNA transcription and translation and was strongly correlated to CD4+, CD8+ T cell-specific gene expression.

CONCLUSIONS

Network based approaches are promising tools to identify potential biomarkers for COPD.

TRIAL REGISTRATION

The ECLIPSE study was funded by GlaxoSmithKline, under ClinicalTrials.gov identifier NCT00292552 and GSK No. SCO104960.

Gene Regulatory Elements, Major Drivers of Human Disease

Gene expression changes, the driving forces for cellular diversity in multicellular organisms, are regulated by a diverse set of gene regulatory elements that direct transcription in specific cells. Mutations in these elements, ranging from chromosomal aberrations to single-nucleotide polymorphisms, are a major cause of human disease. However, we currently have a very limited understanding of how regulatory element genotypes lead to specific phenotypes. In this review, we discuss the various methods of regulatory element identification, the different types of mutations they harbor, and their impact on human disease. We highlight how these variations can affect transcription of multiple genes in gene regulatory networks. In addition, we describe how novel technologies, such as massively parallel reporter assays and CRISPR/Cas9 genome editing, are beginning to provide a better understanding of the functional roles that these elements have and how their alteration can lead to specific phenotypes.

Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease

Recent advances in highly multiplexed immunoassays have allowed systematic large-scale measurement of hundreds of plasma proteins in large cohort studies. In combination with genotyping, such studies offer the prospect to 1) identify mechanisms involved with regulation of protein expression in plasma, and 2) determine whether the plasma proteins are likely to be causally implicated in disease. We report here the results of genome-wide association (GWA) studies of 83 proteins considered relevant to cardiovascular disease (CVD), measured in 3,394 individuals with multiple CVD risk factors. We identified 79 genome-wide significant (p<5e-8) association signals, 55 of which replicated at P<0.0007 in separate validation studies (n = 2,639 individuals). Using automated text mining, manual curation, and network-based methods incorporating information on expression quantitative trait loci (eQTL), we propose plausible causal mechanisms for 25 trans-acting loci, including a potential post-translational regulation of stem cell factor by matrix metalloproteinase 9 and receptor-ligand pairs such as RANK-RANK ligand. Using public GWA study data, we further evaluate all 79 loci for their causal effect on coronary artery disease, and highlight several potentially causal associations. Overall, a majority of the plasma proteins studied showed evidence of regulation at the genetic level. Our results enable future studies of the causal architecture of human disease, which in turn should aid discovery of new drug targets.

Several proteins that circulate in blood have been linked to cardiovascular disease through the use of classic epidemiology and correlation studies. If individuals with higher risk of disease have higher levels of a protein, the protein may be associated with disease. However, this does not necessarily mean that the protein causes disease; it may merely be an innocent bystander or a consequence of the disease process. To establish whether a protein causes disease, a genetic approach, insensitive to reverse causation, can be used. Instead of correlating the levels of the protein itself, gene variants that regulate the protein levels are used in the analysis. This approach requires prior knowledge of which genetic variants are linked to individual proteins. Therefore we completed a map of how common genetic variants affect the blood concentration levels of 83 proteins that have been implicated in cardiovascular disease. By using this map of cause-to-effect findings, we gained insights into the regulation of a majority of the proteins under study and how they relate to risk of coronary artery disease. This study provides a map of genetic regulation of important cardiovascular plasma proteins, insights into their upstream regulatory environment, as well as novel leads for cardiovascular drug development.

Living in an adaptive world: Genomic dissection of the genus Homo and its immune response

More than a decade after the sequencing of the human genome, a deluge of genome-wide population data are generating a portrait of human genetic diversity at an unprecedented level of resolution. Genomic studies have provided new insight into the demographic and adaptive history of our species, Homo sapiens, including its interbreeding with other hominins, such as Neanderthals, and the ways in which natural selection, in its various guises, has shaped genome diversity. These studies, combined with functional genomic approaches, such as the mapping of expression quantitative trait loci, have helped to identify genes, functions, and mechanisms of prime importance for host survival and involved in phenotypic variation and differences in disease risk. This review summarizes new findings in this rapidly developing field, focusing on the human immune response. We discuss the importance of defining the genetic and evolutionary determinants driving immune response variation, and highlight the added value of population genomic approaches in settings relevant to immunity and infection.

Functional variation in allelic methylomes underscores a strong genetic contribution and reveals novel epigenetic alterations in the human epigenome

Background

The functional impact of genetic variation has been extensively surveyed, revealing that genetic changes correlated to phenotypes lie mostly in non-coding genomic regions. Studies have linked allele-specific genetic changes to gene expression, DNA methylation, and histone marks but these investigations have only been carried out in a limited set of samples.

Results

We describe a large-scale coordinated study of allelic and non-allelic effects on DNA methylation, histone mark deposition, and gene expression, detecting the interrelations between epigenetic and functional features at unprecedented resolution. We use information from whole genome and targeted bisulfite sequencing from 910 samples to perform genotype-dependent analyses of allele-specific methylation (ASM) and non-allelic methylation (mQTL). In addition, we introduce a novel genotype-independent test to detect methylation imbalance between chromosomes. Of the ~2.2 million CpGs tested for ASM, mQTL, and genotype-independent effects, we identify ~32% as being genetically regulated (ASM or mQTL) and ~14% as being putatively epigenetically regulated. We also show that epigenetically driven effects are strongly enriched in repressed regions and near transcription start sites, whereas the genetically regulated CpGs are enriched in enhancers. Known imprinted regions are enriched among epigenetically regulated loci, but we also observe several novel genomic regions (e.g., HOX genes) as being epigenetically regulated. Finally, we use our ASM datasets for functional interpretation of disease-associated loci and show the advantage of utilizing naïve T cells for understanding autoimmune diseases.

Conclusions

Our rich catalogue of haploid methylomes across multiple tissues will allow validation of epigenome association studies and exploration of new biological models for allelic exclusion in the human genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1173-7) contains supplementary material, which is available to authorized users.

Pathogenic implications for autoimmune mechanisms derived by comparative eQTL analysis of CD4+ versus CD8+ T cells

Inappropriate activation or inadequate regulation of CD4+ and CD8+ T cells may contribute to the initiation and progression of multiple autoimmune and inflammatory diseases. Studies on disease-associated genetic polymorphisms have highlighted the importance of biological context for many regulatory variants, which is particularly relevant in understanding the genetic regulation of the immune system and its cellular phenotypes. Here we show cell type-specific regulation of transcript levels of genes associated with several autoimmune diseases in CD4+ and CD8+ T cells including a trans-acting regulatory locus at chr12q13.2 containing the rs1131017 SNP in the RPS26 gene. Most remarkably, we identify a common missense variant in IL27, associated with type 1 diabetes that results in decreased functional activity of the protein and reduced expression levels of downstream IRF1 and STAT1 in CD4+ T cells only. Altogether, our results indicate that eQTL mapping in purified T cells provides novel functional insights into polymorphisms and pathways associated with autoimmune diseases.

A functional SNP associated with atopic dermatitis controls cell type-specific methylation of the VSTM1 gene locus

BACKGROUND

Expression quantitative trait loci (eQTL) databases represent a valuable resource to link disease-associated SNPs to specific candidate genes whose gene expression is significantly modulated by the SNP under investigation. We previously identified signal inhibitory receptor on leukocytes-1 (SIRL-1) as a powerful regulator of human innate immune cell function. While it is constitutively high expressed on neutrophils, on monocytes the SIRL-1 surface expression varies strongly between individuals. The underlying mechanism of regulation, its genetic control as well as potential clinical implications had not been explored yet.

METHODS

Whole blood eQTL data of a Chinese cohort was used to identify SNPs regulating the expression of VSTM1, the gene encoding SIRL-1. The genotype effect was validated by flow cytometry (cell surface expression), correlated with electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP) and bisulfite sequencing (C-methylation) and its functional impact studied the inhibition of reactive oxygen species (ROS).

RESULTS

We found a significant association of a single CpG-SNP, rs612529T/C, located in the promoter of VSTM1. Through flow cytometry analysis we confirmed that primarily in the monocytes the protein level of SIRL-1 is strongly associated with genotype of this SNP. In monocytes, the T allele of this SNP facilitates binding of the transcription factors YY1 and PU.1, of which the latter has been recently shown to act as docking site for modifiers of DNA methylation. In line with this notion rs612529T associates with a complete demethylation of the VSTM1 promoter correlating with the allele-specific upregulation of SIRL-1 expression. In monocytes, this upregulation strongly impacts the IgA-induced production of ROS by these cells. Through targeted association analysis we found a significant Meta P value of 1.14 × 10^-6 for rs612529 for association to atopic dermatitis (AD).

CONCLUSION

Low expression of SIRL-1 on monocytes is associated with an increased risk for the manifestation of an inflammatory skin disease. It thus underlines the role of both the cell subset and this inhibitory immune receptor in maintaining immune homeostasis in the skin. Notably, the genetic regulation is achieved by a single CpG-SNP, which controls the overall methylation state of the promoter gene segment.

GWAS for serum galactose-deficient IgA1 implicates critical genes of the O-glycosylation pathway

Aberrant O-glycosylation of serum immunoglobulin A1 (IgA1) represents a heritable pathogenic defect in IgA nephropathy, the most common form of glomerulonephritis worldwide, but specific genetic factors involved in its determination are not known. We performed a quantitative GWAS for serum levels of galactose-deficient IgA1 (Gd-IgA1) in 2,633 subjects of European and East Asian ancestry and discovered two genome-wide significant loci, in C1GALT1 (rs13226913, P = 3.2 x 10⁻¹¹) and C1GALT1C1 (rs5910940, P = 2.7 x 10⁻⁸). These genes encode molecular partners essential for enzymatic O-glycosylation of IgA1. We demonstrated that these two loci explain approximately 7% of variability in circulating Gd-IgA1 in Europeans, but only 2% in East Asians. Notably, the Gd-IgA1-increasing allele of rs13226913 is common in Europeans, but rare in East Asians. Moreover, rs13226913 represents a strong cis-eQTL for C1GALT1 that encodes the key enzyme responsible for the transfer of galactose to O-linked glycans on IgA1. By in vitro siRNA knock-down studies, we confirmed that mRNA levels of both C1GALT1 and C1GALT1C1 determine the rate of secretion of Gd-IgA1 in IgA1-producing cells. Our findings provide novel insights into the genetic regulation of O-glycosylation and are relevant not only to IgA nephropathy, but also to other complex traits associated with O-glycosylation defects, including inflammatory bowel disease, hematologic disease, and cancer.

O-glycosylation is a common type of post-translational modification of proteins; specific abnormalities in the mechanism of O-glycosylation have been implicated in cancer, inflammatory and blood diseases. However, the molecular basis of abnormal O-glycosylation in these complex disorders is not known. We studied the genetic basis of defective O-glycosylation of serum immunoglobulin A1 (IgA1), that represents the key pathogenic defect in IgA nephropathy, the most common form of primary glomerulonephritis worldwide. We report our results of the first genome-wide association study for this trait using serum assays in 2,633 individuals of European and East-Asian ancestry. In our genome scan, we observed two significant signals with large effects, on chromosomes 7p21.3 and Xq24, jointly explaining about 7% of trait variability. These signals implicate two genes that encode molecular partners essential for enzymatic O-glycosylation of IgA1 and mucins, and represent potential new targets for therapy.

Genome-wide association study of primary sclerosing cholangitis identifies new risk loci and quantifies the genetic relationship with inflammatory bowel disease

Primary sclerosing cholangitis (PSC) is a rare progressive disorder leading to bile duct destruction; ∼75% of patients have comorbid inflammatory bowel disease (IBD). We undertook the largest genome-wide association study of PSC (4,796 cases and 19,955 population controls) and identified four new genome-wide significant loci. The most associated SNP at one locus affects splicing and expression of UBASH3A, with the protective allele (C) predicted to cause nonstop-mediated mRNA decay and lower expression of UBASH3A. Further analyses based on common variants suggested that the genome-wide genetic correlation (r_G) between PSC and ulcerative colitis (UC) (r_G = 0.29) was significantly greater than that between PSC and Crohn's disease (CD) (r_G = 0.04) (P = 2.55 × 10^-15). UC and CD were genetically more similar to each other (r_G = 0.56) than either was to PSC (P < 1.0 × 10^-15). Our study represents a substantial advance in understanding of the genetics of PSC.

Systems Genetics as a Tool to Identify Master Genetic Regulators in Complex Disease

Systems genetics stems from systems biology and similarly employs integrative modeling approaches to describe the perturbations and phenotypic effects observed in a complex system. However, in the case of systems genetics the main source of perturbation is naturally occurring genetic variation, which can be analyzed at the systems-level to explain the observed variation in phenotypic traits. In contrast with conventional single-variant association approaches, the success of systems genetics has been in the identification of gene networks and molecular pathways that underlie complex disease. In addition, systems genetics has proven useful in the discovery of master trans-acting genetic regulators of functional networks and pathways, which in many cases revealed unexpected gene targets for disease. Here we detail the central components of a fully integrated systems genetics approach to complex disease, starting from assessment of genetic and gene expression variation, linking DNA sequence variation to mRNA (expression QTL mapping), gene regulatory network analysis and mapping the genetic control of regulatory networks. By summarizing a few illustrative (and successful) examples, we highlight how different data-modeling strategies can be effectively integrated in a systems genetics study.

Identifying Cell Type-Specific Transcription Factors by Integrating ChIP-seq and eQTL Data-Application to Monocyte Gene Regulation

We describe a novel computational approach to identify transcription factors (TFs) that are candidate regulators in a human cell type of interest. Our approach involves integrating cell type-specific expression quantitative trait locus (eQTL) data and TF data from chromatin immunoprecipitation-to-tag-sequencing (ChIP-seq) experiments in cell lines. To test the method, we used eQTL data from human monocytes in order to screen for TFs. Using a list of known monocyte-regulating TFs, we tested the hypothesis that the binding sites of cell type-specific TF regulators would be concentrated in the vicinity of monocyte eQTLs. For each of 397 ChIP-seq data sets, we obtained an enrichment ratio for the number of ChIP-seq peaks that are located within monocyte eQTLs. We ranked ChIP-seq data sets according to their statistical significances for eQTL overlap, and from this ranking, we observed that monocyte-regulating TFs are more highly ranked than would be expected by chance. We identified 27 TFs that had significant monocyte enrichment scores and mapped them into a protein interaction network. Our analysis uncovered two novel candidate monocyte-regulating TFs, BCLAF1 and SIN3A. Our approach is an efficient method to identify candidate TFs that can be used for any cell/tissue type for which eQTL data are available.

Identification of context-dependent expression quantitative trait loci in whole blood

Genetic risk factors often localize to noncoding regions of the genome with unknown effects on disease etiology. Expression quantitative trait loci (eQTLs) help to explain the regulatory mechanisms underlying these genetic associations. Knowledge of the context that determines the nature and strength of eQTLs may help identify cell types relevant to pathophysiology and the regulatory networks underlying disease. Here we generated peripheral blood RNA-seq data from 2,116 unrelated individuals and systematically identified context-dependent eQTLs using a hypothesis-free strategy that does not require previous knowledge of the identity of the modifiers. Of the 23,060 significant cis-regulated genes (false discovery rate (FDR) ≤ 0.05), 2,743 (12%) showed context-dependent eQTL effects. The majority of these effects were influenced by cell type composition. A set of 145 cis-eQTLs depended on type I interferon signaling. Others were modulated by specific transcription factors binding to the eQTL SNPs.