The single-cell eQTLGen consortium

In recent years, functional genomics approaches combining genetic information with bulk RNA-sequencing data have identified the downstream expression effects of disease-associated genetic risk factors through so-called expression quantitative trait locus (eQTL) analysis. Single-cell RNA-sequencing creates enormous opportunities for mapping eQTLs across different cell types and in dynamic processes, many of which are obscured when using bulk methods. Rapid increase in throughput and reduction in cost per cell now allow this technology to be applied to large-scale population genetics studies. To fully leverage these emerging data resources, we have founded the single-cell eQTLGen consortium (sc-eQTLGen), aimed at pinpointing the cellular contexts in which disease-causing genetic variants affect gene expression. Here, we outline the goals, approach and potential utility of the sc-eQTLGen consortium. We also provide a set of study design considerations for future single-cell eQTL studies.

Systematic genetic analysis of the MHC region reveals mechanistic underpinnings of HLA type associations with disease

The MHC region is highly associated with autoimmune and infectious diseases. Here we conduct an in-depth interrogation of associations between genetic variation, gene expression and disease. We create a comprehensive map of regulatory variation in the MHC region using WGS from 419 individuals to call eight-digit HLA types and RNA-seq data from matched iPSCs. Building on this regulatory map, we explored GWAS signals for 4083 traits, detecting colocalization for 180 disease loci with eQTLs. We show that eQTL analyses taking HLA type haplotypes into account have substantially greater power compared with only using single variants. We examined the association between the 8.1 ancestral haplotype and delayed colonization in Cystic Fibrosis, postulating that downregulation of RNF5 expression is the likely causal mechanism. Our study provides insights into the genetic architecture of the MHC region and pinpoints disease associations that are due to differential expression of HLA genes and non-HLA genes.

Hematopoietic Npc1 mutation shifts gut microbiota composition in Ldlr-/- mice on a high-fat, high-cholesterol diet

While the link between diet-induced changes in gut microbiota and lipid metabolism in metabolic syndrome (MetS) has been established, the contribution of host genetics is rather unexplored. As several findings suggested a role for the lysosomal lipid transporter Niemann-Pick type C1 (NPC1) in macrophages during MetS, we here explored whether a hematopoietic Npc1 mutation, induced via bone marrow transplantation, influences gut microbiota composition in low-density lipoprotein receptor knockout (Ldlr^-/-) mice fed a high-fat, high-cholesterol (HFC) diet for 12 weeks. Ldlr^-/- mice fed a HFC diet mimic a human plasma lipoprotein profile and show features of MetS, providing a model to explore the role of host genetics on gut microbiota under MetS conditions. Fecal samples were used to profile the microbial composition by 16 s ribosomal RNA gene sequencing. The hematopoietic Npc1 mutation shifted the gut microbiota composition and increased microbial richness and diversity. Variations in plasma lipid levels correlated with microbial diversity and richness as well as with several bacterial genera. This study suggests that host genetic influences on lipid metabolism affect the gut microbiome under MetS conditions. Future research investigating the role of host genetics on gut microbiota might therefore lead to identification of diagnostic and therapeutic targets for MetS.

Epigenome-wide association study of body mass index, and the adverse outcomes of adiposity

Overweight and obesity affect ~1.5 billion people worldwide, and are major risk factors for type-2 diabetes (T2D), cardiovascular disease and related metabolic and inflammatory disturbances.1,2 Although the mechanisms linking adiposity to its clinical sequelae are poorly understood, recent studies suggest that adiposity may influence DNA methylation,3–6 a key regulator of gene expression and molecular phenotype.7 Here we use epigenome-wide association to show that body mass index (BMI, a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci at P<1x10^-7, range P=9.2x10^-8 to 6.0x10^-46; N=10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find the methylation loci are enriched for functional genomic features in multiple tissues (P<0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P<9.0x10^-6, range P=5.5x10^-6 to 6.1x10^-35, N=1,785 samples). The methylation loci identified highlight genes involved in lipid and lipoprotein metabolism, substrate transport, and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future type-2 diabetes (relative risk per 1SD increase in Methylation Risk Score: 2.3 [2.07-2.56]; P=1.1x10^-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type-2 diabetes and other adverse clinical consequences of obesity.