diffloop: a computational framework for identifying and analyzing differential DNA loops from sequencing data

Abstract Summary The 3D architecture of DNA within the nucleus is a key determinant of interactions between genes, regulatory elements, and transcriptional machinery. As a result, differences in DNA looping structure are associated with variation in gene expression and cell state. To systematically assess changes in DNA looping architecture between samples, we introduce diffloop, an R/Bioconductor package that provides a suite of functions for the quality control, statistical testing, annotation, and visualization of DNA loops. We demonstrate this functionality by detecting differences between ENCODE ChIA-PET samples and relate looping to variability in epigenetic state. Availability and implementation Diffloop is implemented as an R/Bioconductor package available at https://bioconductor.org/packages/release/bioc/html/diffloop.html Contact [email protected] Supplementary information Supplementary data are available at Bioinformatics online

Heritability of fetal hemoglobin, white cell count, and other clinical traits from a sickle cell disease family cohort

Sickle cell disease (SCD) is the most common monogenic disorder in the world. Notably, there is extensive clinical heterogeneity in SCD that cannot be fully accounted for by known factors, and in particular, the extent to which the phenotypic diversity of SCD can be explained by genetic variation has not been reliably quantified. Here, in a family-based cohort of 449 patients with SCD and 755 relatives, we first show that 5 known modifiers affect 11 adverse outcomes in SCD to varying degrees. We then utilize a restricted maximum likelihood procedure to estimate the heritability of 20 hematologic traits, including fetal hemoglobin (HbF) and white blood cell count (WBC), in the clinically relevant context of inheritance from healthy carriers to SCD patients. We report novel estimations of heritability for HbF at 31.6% (±5.4%) and WBC at 41.2% (±6.8%) in our cohort. Finally, we demonstrate shared genetic bases between HbF, WBC, and other hematologic traits, but surprisingly little overlap between HbF and WBC themselves. In total, our analyses show that HbF and WBC have significant heritable components among individuals with SCD and their relatives, demonstrating the value of using family-based studies to better understand modifiers of SCD

Single Cell Transcriptomics Implicate Novel Monocyte and T Cell Immune Dysregulation in Sarcoidosis

Sarcoidosis is a systemic inflammatory disease characterized by infiltration of immune cells into granulomas. Previous gene expression studies using heterogeneous cell mixtures lack insight into cell-type-specific immune dysregulation. We performed the first single-cell RNA-sequencing study of sarcoidosis in peripheral immune cells in 48 patients and controls. Following unbiased clustering, differentially expressed genes were identified for 18 cell types and bioinformatically assessed for function and pathway enrichment. Our results reveal persistent activation of circulating classical monocytes with subsequent upregulation of trafficking molecules. Specifically, classical monocytes upregulated distinct markers of activation including adhesion molecules, pattern recognition receptors, and chemokine receptors, as well as enrichment of immunoregulatory pathways HMGB1, mTOR, and ephrin receptor signaling. Predictive modeling implicated TGFβ and mTOR signaling as drivers of persistent monocyte activation. Additionally, sarcoidosis T cell subsets displayed patterns of dysregulation. CD4 naïve T cells were enriched for markers of apoptosis and Th17/T(reg) differentiation, while effector T cells showed enrichment of anergy-related pathways. Differentially expressed genes in regulatory T cells suggested dysfunctional p53, cell death, and TNFR2 signaling. Using more sensitive technology and more precise units of measure, we identify cell-type specific, novel inflammatory and regulatory pathways. Based on our findings, we suggest a novel model involving four convergent arms of dysregulation: persistent hyperactivation of innate and adaptive immunity via classical monocytes and CD4 naïve T cells, regulatory T cell dysfunction, and effector T cell anergy. We further our understanding of the immunopathology of sarcoidosis and point to novel therapeutic targets

Transcript-indexed ATAC-seq for precision immune profiling.

T cells create vast amounts of diversity in the genes that encode their T cell receptors (TCRs), which enables individual clones to recognize specific peptide-major histocompatibility complex (MHC) ligands. Here we combined sequencing of the TCR-encoding genes with assay for transposase-accessible chromatin with sequencing (ATAC-seq) analysis at the single-cell level to provide information on the TCR specificity and epigenomic state of individual T cells. By using this approach, termed transcript-indexed ATAC-seq (T-ATAC-seq), we identified epigenomic signatures in immortalized leukemic T cells, primary human T cells from healthy volunteers and primary leukemic T cells from patient samples. In peripheral blood CD4+ T cells from healthy individuals, we identified cis and trans regulators of naive and memory T cell states and found substantial heterogeneity in surface-marker-defined T cell populations. In patients with a leukemic form of cutaneous T cell lymphoma, T-ATAC-seq enabled identification of leukemic and nonleukemic regulatory pathways in T cells from the same individual by allowing separation of the signals that arose from the malignant clone from the background T cell noise. Thus, T-ATAC-seq is a new tool that enables analysis of epigenomic landscapes in clonal T cells and should be valuable for studies of T cell malignancy, immunity and immunotherapy

Polygenic Risk Assessment Reveals Pleiotropy between Sarcoidosis and Inflammatory Disorders in the Context of Genetic Ancestry

Sarcoidosis is a complex disease of unknown etiology characterized by the presence of granulomatous inflammation. Though various immune system pathways have been implicated in disease, the relationship between the genetic determinants of sarcoidosis and other inflammatory disorders has not been characterized. Herein, we examined the degree of genetic pleiotropy common to sarcoidosis and other inflammatory disorders to identify shared pathways and disease systems pertinent to sarcoidosis onset. To achieve this, we quantify the association of common variant polygenic risk scores from nine complex inflammatory disorders with sarcoidosis risk. Enrichment analyses of genes implicated in pleiotropic associations were further used to elucidate candidate pathways. In European-Americans, we identify significant pleiotropy between risk of sarcoidosis and risk of asthma (R2=2.03%; p=8.89×10−9), celiac disease (R2=2.03%; p=8.21×10−9), primary biliary cirrhosis (R2=2.43%; p=2.01×10−10), and rheumatoid arthritis (R2=4.32%; p=2.50×10−17). These associations validate in African Americans only after accounting for the proportion of genome-wide European ancestry, where we demonstrate similar effects of polygenic risk for African-Americans with the highest levels of European ancestry. Variants and genes implicated in European-American pleiotropic associations were enriched for pathways involving interleukin-12, interleukin-27, and cell adhesion molecules, corroborating the hypothesized immunopathogenesis of disease

Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types.

We introduce an approach to identify disease-relevant tissues and cell types by analyzing gene expression data together with genome-wide association study (GWAS) summary statistics. Our approach uses stratified linkage disequilibrium (LD) score regression to test whether disease heritability is enriched in regions surrounding genes with the highest specific expression in a given tissue. We applied our approach to gene expression data from several sources together with GWAS summary statistics for 48 diseases and traits (average N = 169,331) and found significant tissue-specific enrichments (false discovery rate (FDR) < 5%) for 34 traits. In our analysis of multiple tissues, we detected a broad range of enrichments that recapitulated known biology. In our brain-specific analysis, significant enrichments included an enrichment of inhibitory over excitatory neurons for bipolar disorder, and excitatory over inhibitory neurons for schizophrenia and body mass index. Our results demonstrate that our polygenic approach is a powerful way to leverage gene expression data for interpreting GWAS signals

Functional dissection of inherited non-coding variation influencing multiple myeloma risk

Funding Information: This work was supported by grants from the Knut and Alice Wallenberg Foundation (2012.0193 and 2017.0436), the Swedish Research Council (2017-02023 and 2018-00424), the Swedish Cancer Society (2017/265), the Nordic Cancer Union (R217-A13329-18-S65), Arne and Inga-Britt Lundberg’s Stiftelse (2017-0055), European Research Council (EU-MSCA-COFUND 754299 CanFaster), Myeloma UK and Cancer Research UK (C1298/A8362), The National Institute of Health (R01 DK103794 and R01HL146500), the New York Stem Cell Foundation, a gift from the Lodish Family to Boston Children’s Hospital, and Mr. Ralph Stockwell. We thank Ellinor Johnsson for her assistance between 2011 and 2020. We are indebted to the patients who participated in the study. Publisher Copyright: © 2022, The Author(s).Thousands of non-coding variants have been associated with increased risk of human diseases, yet the causal variants and their mechanisms-of-action remain obscure. In an integrative study combining massively parallel reporter assays (MPRA), expression analyses (eQTL, meQTL, PCHiC) and chromatin accessibility analyses in primary cells (caQTL), we investigate 1,039 variants associated with multiple myeloma (MM). We demonstrate that MM susceptibility is mediated by gene-regulatory changes in plasma cells and B-cells, and identify putative causal variants at six risk loci (SMARCD3, WAC, ELL2, CDCA7L, CEP120, and PREX1). Notably, three of these variants co-localize with significant plasma cell caQTLs, signaling the presence of causal activity at these precise genomic positions in an endogenous chromosomal context in vivo. Our results provide a systematic functional dissection of risk loci for a hematologic malignancy.Peer reviewe

The Polygenic and Monogenic Basis of Blood Traits and Diseases

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation. Analysis of blood cell traits in the UK Biobank and other cohorts illuminates the full genetic architecture of hematopoietic phenotypes, with evidence supporting the omnigenic model for complex traits and linking polygenic burden with monogenic blood diseases

The Polygenic and Monogenic Basis of Blood Traits and Diseases

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.</p

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p &lt; 5 × 10−9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies. Delineation of the genetic architecture of hematological traits in a multi-ethnic dataset allows identification of rare variants with strong effects specific to non-European populations and improved fine mapping of GWAS variants using the trans-ethnic approach