Search

Abstract The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Higher prevalence and worst outcome have been reported among people with systemic lupus erythematosus with non-European ancestries, with both genetic and socioeconomic variables as contributing factors. In Mexico, studies assessing the inequities related to quality of life for Systemic Lupus Erythematosus patients remain sparse. This study aims to assess the inequities related to quality of life in a cohort of Mexican people with SLE. This study included 942 individuals with SLE from the Mexican Lupus Registry (LupusRGMX) and two healthy control groups. Self-answered surveys were collected via the Research Electronic Data Capture platform between May 2021 and January 2023. Data was analyzed as a cross-sectional study. A random forest model was implemented to assess potential predictive variables. Permutation tests were performed to analyze the effect health providers had on diagnosis lag and quality of life's differences among socioeconomic levels. Partial correlation analysis between the number of patients and rheumatologists registered was also performed. Systemic Lupus Erythematosus participants had significantly lower quality of life than healthy people (p-values < 0.0001). Socioeconomic status, delay in diagnosis, and corticosteroid consumption were the factors that influenced QoL the most (RMSE = 9.53 with the importance variable validated); lower quality of life was associated with lower socioeconomic status (p-value < 0.0001). Disparities in health services were reflected in longer diagnosis time among people with public health providers (p-value = 0.0419). A significant association between diagnosed patients and available rheumatologists by geographical state was observed (ρ = 0.4, p-value = 0.0259), which can be translated into restricted access to specialists. Since most of our cohort exhibited low socioeconomic status, it is important to consider them as a vulnerable population; this study settles the necessity to deepen the effects of the socioeconomic disparities, allowing to design public policies and strategies aimed to reduce Systemic Lupus Erythematosus disparities, therefore improving quality of life of Mexican people with Systemic Lupus Erythematosus.

Abstract Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (encoded by IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKε. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a regulatory pathway involving rs2297550, Ikaros, and IKKε implicated by human genetics in risk for SLE.

Abstract Fine-mapping and functional studies implicate rs117701653, a non-coding single nucleotide polymorphism in the CD28/CTLA4/ICOS locus, as a risk variant for rheumatoid arthritis and type 1 diabetes. Here, using DNA pulldown, mass spectrometry, genome editing and eQTL analysis, we establish that the disease-associated risk allele is functional, reducing affinity for the inhibitory chromosomal regulator SMCHD1 to enhance expression of inducible T-cell costimulator (ICOS) in memory CD4+ T cells from healthy donors. Higher ICOS expression is paralleled by an increase in circulating T peripheral helper (Tph) cells and, in rheumatoid arthritis patients, of blood and joint fluid Tph cells as well as circulating plasmablasts. Correspondingly, ICOS ligation and carriage of the rs117701653 risk allele accelerate T cell differentiation into CXCR5-PD-1high Tph cells producing IL-21 and CXCL13. Thus, mechanistic dissection of a functional non-coding variant in human autoimmunity discloses a previously undefined pathway through which ICOS regulates Tph development and abundance.

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

IntroductionChronic rhinosinusitis (CRS) is a chronic inflammatory disease of the sinonasal mucosa with distinct endotypes including type 2 (T2) high eosinophilic CRS with nasal polyps (eCRSwNP), T2 low non-eosinophilic CRS with nasal polyps (neCRSwNP), and CRS without nasal polyps (CRSsNP).MethodsGiven the heterogeneity of disease, we hypothesized that assessment of single cell RNA sequencing (scRNA-seq) across this spectrum of disease would reveal connections between infiltrating and activated immune cells and the epithelial and stromal populations that reside in sinonasal tissue.ResultsHere we find increased expression of genes encoding glycolytic enzymes in epithelial cells (EpCs), stromal cells, and memory T-cell subsets from patients with eCRSwNP, as compared to healthy controls. In basal EpCs, this is associated with a program of cell motility and Rho GTPase effector expression. Across both stromal and immune subsets, glycolytic programming was associated with extracellular matrix interactions, proteoglycan generation, and collagen formation. Furthermore, we report increased cell-cell interactions between EpCs and stromal/immune cells in eCRSwNP compared to healthy control tissue, and we nominate candidate receptor-ligand pairs that may drive tissue remodeling.DiscussionThese findings support a role for glycolytic reprograming in T2-elicited tissue remodeling and implicate increased cellular crosstalk in eCRSwNP.

T peripheral helper (Tph) cells, identified in the synovium of adults with seropositive rheumatoid arthritis, drive B cell maturation and antibody production in non-lymphoid tissues. We sought to determine if similarly dysregulated T cell-B cell interactions underlie another form of inflammatory arthritis, juvenile oligoarthritis (oligo JIA). Clonally expanded Tph cells able to promote B cell antibody production preferentially accumulated in the synovial fluid (SF) of oligo JIA patients with antinuclear antibodies (ANA) compared to autoantibody-negative patients. Single-cell transcriptomics enabled further definition of the Tph gene signature in inflamed tissues and showed that Tph cells from ANA-positive patients upregulated genes associated with B cell help to a greater extent than patients without autoantibodies. T cells that co-expressed regulatory T and B cell-help factors were identified. The phenotype of these Tph-like Treg cells suggests an ability to restrain T cell-B cell interactions in tissues. Our findings support the central role of disordered T cell-help to B cells in autoantibody-positive arthritides.

Oligoarticular juvenile idiopathic arthritis (oligo JIA) is the most common form of chronic inflammatory arthritis in children, yet the cause of this disease remains unknown. To understand immune responses in oligo JIA, we immunophenotyped synovial fluid T cells with flow cytometry, bulk RNA-Seq, single-cell RNA-Seq (scRNA-Seq), DNA methylation studies, and Treg suppression assays. In synovial fluid, CD4+, CD8+, and γδ T cells expressed Th1-related markers, whereas Th17 cells were not enriched. Th1 skewing was prominent in CD4+ T cells, including Tregs, and was associated with severe disease. Transcriptomic studies confirmed a Th1 signature in CD4+ T cells from synovial fluid. The regulatory gene expression signature was preserved in Tregs, even those exhibiting Th1 polarization. These Th1-like Tregs maintained Treg-specific methylation patterns and suppressive function, supporting the stability of this Treg population in the joint. Although synovial fluid CD4+ T cells displayed an overall Th1 phenotype, scRNA-Seq uncovered heterogeneous effector and regulatory subpopulations, including IFN-induced Tregs, peripheral helper T cells, and cytotoxic CD4+ T cells. In conclusion, oligo JIA is characterized by Th1 polarization that encompasses Tregs but does not compromise their regulatory identity. Targeting Th1-driven inflammation and augmenting Treg function may represent important therapeutic approaches in oligo JIA.

Between 5 and 15% of latent Mycobacterium tuberculosis infections develop into active tuberculosis (TB). Here, Luo et al. report a genome-wide association study for early TB progression in a total of 4002 active TB cases and their household contacts in Peru and they identify a locus on 3q23 in which ATP1B3 is mapped as the likely causal gene.

Innate T cells (ITC) contain many subsets and are poised to promptly respond to antigens and pathogens, but how this poised state is maintained is still unclear. Here the authors perform single-cell RNA-seq to align the various ITC subsets along an ‘innateness gradient’ that is associated with changes in proliferation and effector functions.

Abstract Background Cytokines are critical to human disease and are attractive therapeutic targets given their widespread influence on gene regulation and transcription. Defining the downstream regulatory mechanisms influenced by cytokines is central to defining drug and disease mechanisms. One promising strategy is to use interactions between expression quantitative trait loci (eQTLs) and cytokine levels to define target genes and mechanisms. Results In a clinical trial for anti-IL-6 in patients with systemic lupus erythematosus, we measure interferon (IFN) status, anti-IL-6 drug exposure, and whole blood genome-wide gene expression at three time points. We show that repeat transcriptomic measurements increases the number of cis eQTLs identified compared to using a single time point. We observe a statistically significant enrichment of in vivo eQTL interactions with IFN status and anti-IL-6 drug exposure and find many novel interactions that have not been previously described. Finally, we find transcription factor binding motifs interrupted by eQTL interaction SNPs, which point to key regulatory mediators of these environmental stimuli and therefore potential therapeutic targets for autoimmune diseases. In particular, genes with IFN interactions are enriched for ISRE binding site motifs, while those with anti-IL-6 interactions are enriched for IRF4 motifs. Conclusions This study highlights the potential to exploit clinical trial data to discover in vivo eQTL interactions with therapeutically relevant environmental variables.

Abstract Background Detailed molecular analyses of cells from rheumatoid arthritis (RA) synovium hold promise in identifying cellular phenotypes that drive tissue pathology and joint damage. The Accelerating Medicines Partnership RA/SLE Network aims to deconstruct autoimmune pathology by examining cells within target tissues through multiple high-dimensional assays. Robust standardized protocols need to be developed before cellular phenotypes at a single cell level can be effectively compared across patient samples. Methods Multiple clinical sites collected cryopreserved synovial tissue fragments from arthroplasty and synovial biopsy in a 10% DMSO solution. Mechanical and enzymatic dissociation parameters were optimized for viable cell extraction and surface protein preservation for cell sorting and mass cytometry, as well as for reproducibility in RNA sequencing (RNA-seq). Cryopreserved synovial samples were collectively analyzed at a central processing site by a custom-designed and validated 35-marker mass cytometry panel. In parallel, each sample was flow sorted into fibroblast, T-cell, B-cell, and macrophage suspensions for bulk population RNA-seq and plate-based single-cell CEL-Seq2 RNA-seq. Results Upon dissociation, cryopreserved synovial tissue fragments yielded a high frequency of viable cells, comparable to samples undergoing immediate processing. Optimization of synovial tissue dissociation across six clinical collection sites with ~ 30 arthroplasty and ~ 20 biopsy samples yielded a consensus digestion protocol using 100 μg/ml of Liberase™ TL enzyme preparation. This protocol yielded immune and stromal cell lineages with preserved surface markers and minimized variability across replicate RNA-seq transcriptomes. Mass cytometry analysis of cells from cryopreserved synovium distinguished diverse fibroblast phenotypes, distinct populations of memory B cells and antibody-secreting cells, and multiple CD4+ and CD8+ T-cell activation states. Bulk RNA-seq of sorted cell populations demonstrated robust separation of synovial lymphocytes, fibroblasts, and macrophages. Single-cell RNA-seq produced transcriptomes of over 1000 genes/cell, including transcripts encoding characteristic lineage markers identified. Conclusions We have established a robust protocol to acquire viable cells from cryopreserved synovial tissue with intact transcriptomes and cell surface phenotypes. A centralized pipeline to generate multiple high-dimensional analyses of synovial tissue samples collected across a collaborative network was developed. Integrated analysis of such datasets from large patient cohorts may help define molecular heterogeneity within RA pathology and identify new therapeutic targets and biomarkers.

The importance of natural gene expression variation for human behavior is undisputed, but its impact on circadian physiology remains mostly unexplored. Using umbilical cord fibroblasts, we have determined by genome-wide association how common genetic variation impacts upon cellular circadian function. Gene set enrichment points to differences in protein catabolism as one major source of clock variation in humans. The two most significant alleles regulated expression of COPS7B, a subunit of the COP9 signalosome. We further show that the signalosome complex is imported into the nucleus in timed fashion to stabilize the essential circadian protein BMAL1, a novel mechanism to oppose its proteasome-mediated degradation. Thus, circadian clock properties depend in part upon a genetically-encoded competition between stabilizing and destabilizing forces, and genetic alterations in these mechanisms provide one explanation for human chronotype.

Understanding how genetic variation affects distinct cellular phenotypes, such as gene expression levels, alternative splicing and DNA methylation levels, is essential for better understanding of complex diseases and traits. Furthermore, how inter-individual variation of DNA methylation is associated to gene expression is just starting to be studied. In this study, we use the GenCord cohort of 204 newborn Europeans' lymphoblastoid cell lines, T-cells and fibroblasts derived from umbilical cords. The samples were previously genotyped for 2.5 million SNPs, mRNA-sequenced, and assayed for methylation levels in 482,421 CpG sites. We observe that methylation sites associated to expression levels are enriched in enhancers, gene bodies and CpG island shores. We show that while the correlation between DNA methylation and gene expression can be positive or negative, it is very consistent across cell-types. However, this epigenetic association to gene expression appears more tissue-specific than the genetic effects on gene expression or DNA methylation (observed in both sharing estimations based on P-values and effect size correlations between cell-types). This predominance of genetic effects can also be reflected by the observation that allele specific expression differences between individuals dominate over tissue-specific effects. Additionally, we discover genetic effects on alternative splicing and interestingly, a large amount of DNA methylation correlating to alternative splicing, both in a tissue-specific manner. The locations of the SNPs and methylation sites involved in these associations highlight the participation of promoter proximal and distant regulatory regions on alternative splicing. Overall, our results provide high-resolution analyses showing how genome sequence variation has a broad effect on cellular phenotypes across cell-types, whereas epigenetic factors provide a secondary layer of variation that is more tissue-specific. Furthermore, the details of how this tissue-specificity may vary across inter-relations of molecular traits, and where these are occurring, can yield further insights into gene regulation and cellular biology as a whole.

DNA methylation is an essential epigenetic mark whose role in gene regulation and its dependency on genomic sequence and environment are not fully understood. In this study we provide novel insights into the mechanistic relationships between genetic variation, DNA methylation and transcriptome sequencing data in three different cell-types of the GenCord human population cohort. We find that the association between DNA methylation and gene expression variation among individuals are likely due to different mechanisms from those establishing methylation-expression patterns during differentiation. Furthermore, cell-type differential DNA methylation may delineate a platform in which local inter-individual changes may respond to or act in gene regulation. We show that unlike genetic regulatory variation, DNA methylation alone does not significantly drive allele specific expression. Finally, inferred mechanistic relationships using genetic variation as well as correlations with TF abundance reveal both a passive and active role of DNA methylation to regulatory interactions influencing gene expression.

The genetic basis of gene expression variation has long been studied with the aim to understand the landscape of regulatory variants, but also more recently to assist in the interpretation and elucidation of disease signals. To date, many studies have looked in specific tissues and population-based samples, but there has been limited assessment of the degree of inter-population variability in regulatory variation. We analyzed genome-wide gene expression in lymphoblastoid cell lines from a total of 726 individuals from 8 global populations from the HapMap3 project and correlated gene expression levels with HapMap3 SNPs located in cis to the genes. We describe the influence of ancestry on gene expression levels within and between these diverse human populations and uncover a non-negligible impact on global patterns of gene expression. We further dissect the specific functional pathways differentiated between populations. We also identify 5,691 expression quantitative trait loci (eQTLs) after controlling for both non-genetic factors and population admixture and observe that half of the cis-eQTLs are replicated in one or more of the populations. We highlight patterns of eQTL-sharing between populations, which are partially determined by population genetic relatedness, and discover significant sharing of eQTL effects between Asians, European-admixed, and African subpopulations. Specifically, we observe that both the effect size and the direction of effect for eQTLs are highly conserved across populations. We observe an increasing proximity of eQTLs toward the transcription start site as sharing of eQTLs among populations increases, highlighting that variants close to TSS have stronger effects and therefore are more likely to be detected across a wider panel of populations. Together these results offer a unique picture and resource of the degree of differentiation among human populations in functional regulatory variation and provide an estimate for the transferability of complex trait variants across populations.

Population-scale genome sequencing allows the characterization of functional effects of a broad spectrum of genetic variants underlying human phenotypic variation. Here, we investigate the influence of rare and common genetic variants on gene expression patterns, using variants identified from sequencing data from the 1000 genomes project in an African and European population sample and gene expression data from lymphoblastoid cell lines. We detect comparable numbers of expression quantitative trait loci (eQTLs) when compared to genotypes obtained from HapMap 3, but as many as 80% of the top expression quantitative trait variants (eQTVs) discovered from 1000 genomes data are novel. The properties of the newly discovered variants suggest that mapping common causal regulatory variants is challenging even with full resequencing data; however, we observe significant enrichment of regulatory effects in splice-site and nonsense variants. Using RNA sequencing data, we show that 46.2% of nonsynonymous variants are differentially expressed in at least one individual in our sample, creating widespread potential for interactions between functional protein-coding and regulatory variants. We also use allele-specific expression to identify putative rare causal regulatory variants. Furthermore, we demonstrate that outlier expression values can be due to rare variant effects, and we approximate the number of such effects harboured in an individual by effect size. Our results demonstrate that integration of genomic and RNA sequencing analyses allows for the joint assessment of genome sequence and genome function.