The Identity of Human Tissue-Emigrant CD8+ T Cells.

Lymphocyte migration is essential for adaptive immune surveillance. However, our current understanding of this process is rudimentary, because most human studies have been restricted to immunological analyses of blood and various tissues. To address this knowledge gap, we used an integrated approach to characterize tissue-emigrant lineages in thoracic duct lymph (TDL). The most prevalent immune cells in human and non-human primate efferent lymph were T cells. Cytolytic CD8+ T cell subsets with effector-like epigenetic and transcriptional signatures were clonotypically skewed and selectively confined to the intravascular circulation, whereas non-cytolytic CD8+ T cell subsets with stem-like epigenetic and transcriptional signatures predominated in tissues and TDL. Moreover, these anatomically distinct gene expression profiles were recapitulated within individual clonotypes, suggesting parallel differentiation programs independent of the expressed antigen receptor. Our collective dataset provides an atlas of the migratory immune system and defines the nature of tissue-emigrant CD8+ T cells that recirculate via TDL.

Distinct Regulation of Th17 and Th1 Cell Differentiation by Glutaminase-Dependent Metabolism.

Activated T cells differentiate into functional subsets with distinct metabolic programs. Glutaminase (GLS) converts glutamine to glutamate to support the tricarboxylic acid cycle and redox and epigenetic reactions. Here, we identify a key role for GLS in T cell activation and specification. Though GLS deficiency diminished initial T cell activation and proliferation and impaired differentiation of Th17 cells, loss of GLS also increased Tbet to promote differentiation and effector function of CD4 Th1 and CD8 CTL cells. This was associated with altered chromatin accessibility and gene expression, including decreased PIK3IP1 in Th1 cells that sensitized to IL-2-mediated mTORC1 signaling. In vivo, GLS null T cells failed to drive Th17-inflammatory diseases, and Th1 cells had initially elevated function but exhausted over time. Transient GLS inhibition, however, led to increased Th1 and CTL T cell numbers. Glutamine metabolism thus has distinct roles to promote Th17 but constrain Th1 and CTL effector cell differentiation.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR.

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

Genetic and Epigenetic Fine Mapping of Complex Trait Associated Loci in the Human Liver.

Deciphering the impact of genetic variation on gene regulation is fundamental to understanding common, complex human diseases. Although histone modifications are important markers of gene regulatory elements of the genome, any specific histone modification has not been assayed in more than a few individuals in the human liver. As a result, the effects of genetic variation on histone modification states in the liver are poorly understood. Here, we generate the most comprehensive genome-wide dataset of two epigenetic marks, H3K4me3 and H3K27ac, and annotate thousands of putative regulatory elements in the human liver. We integrate these findings with genome-wide gene expression data collected from the same human liver tissues and high-resolution promoter-focused chromatin interaction maps collected from human liver-derived HepG2 cells. We demonstrate widespread functional consequences of natural genetic variation on putative regulatory element activity and gene expression levels. Leveraging these extensive datasets, we fine-map a total of 74 GWAS loci that have been associated with at least one complex phenotype. Our results reveal a repertoire of genes and regulatory mechanisms governing complex disease development and further the basic understanding of genetic and epigenetic regulation of gene expression in the human liver tissue.

Canonical Notch signaling is required for bone morphogenetic protein-mediated human osteoblast differentiation.

Osteoblast differentiation of bone marrow-derived human mesenchymal stem cells (hMSC) can be induced by stimulation with canonical Notch ligand, Jagged1, or bone morphogenetic proteins (BMPs). However, it remains elusive how these two pathways lead to the same phenotypic outcome. Since Runx2 is regarded as a master regulator of osteoblastic differentiation, we targeted Runx2 with siRNA in hMSC. This abrogated both Jagged1 and BMP2 mediated osteoblastic differentiation, confirming the fundamental role for Runx2. However, while BMP stimulation increased Runx2 and downstream Osterix protein expression, Jagged1 treatment failed to upregulate either, suggesting that canonical Notch signals require basal Runx2 expression. To fully understand the transcriptomic profile of differentiating osteoblasts, RNA sequencing was performed in cells stimulated with BMP2 or Jagged1. There was common upregulation of ALPL and extracellular matrix genes, such as ACAN, HAS3, MCAM, and OLFML2B. Intriguingly, genes encoding components of Notch signaling (JAG1, HEY2, and HES4) were among the top 10 genes upregulated by both stimuli. Indeed, ALPL expression occurred concurrently with Notch activation and inhibiting Notch activity for up to 24 hours after BMP administration with DAPT (a gamma secretase inhibitor) completely abrogated hMSC osteoblastogenesis. Concordantly, RBPJ (recombination signal binding protein for immunoglobulin kappa J region, a critical downstream modulator of Notch signals) binding could be demonstrated within the ALPL and SP7 promoters. As such, siRNA-mediated ablation of RBPJ decreased BMP-mediated osteoblastogenesis. Finally, systemic Notch inhibition using diabenzazepine (DBZ) reduced BMP2-induced calvarial bone healing in mice supporting the critical regulatory role of Notch signaling in BMP-induced osteoblastogenesis.

The Loss of TET2 Promotes CD8+ T Cell Memory Differentiation.

T cell differentiation requires appropriate regulation of DNA methylation. In this article, we demonstrate that the methylcytosine dioxygenase ten-eleven translocation (TET)2 regulates CD8+ T cell differentiation. In a murine model of acute viral infection, TET2 loss promotes early acquisition of a memory CD8+ T cell fate in a cell-intrinsic manner without disrupting Ag-driven cell expansion or effector function. Upon secondary recall, TET2-deficient memory CD8+ T cells demonstrate superior pathogen control. Genome-wide methylation analysis identified a number of differentially methylated regions in TET2-deficient versus wild-type CD8+ T cells. These differentially methylated regions did not occur at the loci of differentially expressed memory markers; rather, several hypermethylated regions were identified in known transcriptional regulators of CD8+ T cell memory fate. Together, these data demonstrate that TET2 is an important regulator of CD8+ T cell fate decisions.

Leveraging epigenomics and contactomics data to investigate SNP pairs in GWAS.

Although Genome Wide Association Studies (GWAS) have led to many valuable insights into the genetic bases of common diseases over the past decade, the issue of missing heritability has surfaced, as the discovered main effect genetic variants found to date do not account for much of a trait's predicted genetic component. We present a workflow, integrating epigenomics and topologically associating domain data, aimed at discovering trait-associated SNP pairs from GWAS where neither SNP achieved independent genome-wide significance. Each analyzed SNP pair consists of one SNP in a putative active enhancer and another SNP in a putative physically interacting gene promoter in a trait-relevant tissue. As a proof-of-principle case study, we used this approach to identify focused collections of SNP pairs that we analyzed in three independent Type 2 diabetes (T2D) GWAS. This approach led us to discover 35 significant SNP pairs, encompassing both novel signals and signals for which we have found orthogonal support from other sources. Nine of these pairs are consistent with eQTL results, two are consistent with our own capture C experiments, and seven involve signals supported by recent T2D literature.

DNA binding by FOXP3 domain-swapped dimer suggests mechanisms of long-range chromosomal interactions.

FOXP3 is a lineage-specific transcription factor that is required for regulatory T cell development and function. In this study, we determined the crystal structure of the FOXP3 forkhead domain bound to DNA. The structure reveals that FOXP3 can form a stable domain-swapped dimer to bridge DNA in the absence of cofactors, suggesting that FOXP3 may play a role in long-range gene interactions. To test this hypothesis, we used circular chromosome conformation capture coupled with high throughput sequencing (4C-seq) to analyze FOXP3-dependent genomic contacts around a known FOXP3-bound locus, Ptpn22. Our studies reveal that FOXP3 induces significant changes in the chromatin contacts between the Ptpn22 locus and other Foxp3-regulated genes, reflecting a mechanism by which FOXP3 reorganizes the genome architecture to coordinate the expression of its target genes. Our results suggest that FOXP3 mediates long-range chromatin interactions as part of its mechanisms to regulate specific gene expression in regulatory T cells.

The type 2 diabetes presumed causal variant within TCF7L2 resides in an element that controls the expression of ACSL5.

AIMS/HYPOTHESIS: One of the most strongly associated type 2 diabetes loci reported to date resides within the TCF7L2 gene. Previous studies point to the T allele of rs7903146 in intron 3 as the causal variant at this locus. We aimed to identify the actual gene(s) under the influence of this variant. METHODS: Using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein-9 nuclease, we generated a 1.4 kb deletion of the genomic region harbouring rs7903146 in the HCT116 cell line, followed by global gene expression analysis. We then carried out a combination of circularised chromosome conformation capture (4C) and Capture C in cell lines, HCT116 and NCM460 in order to ascertain which promoters of these perturbed genes made consistent physical contact with this genomic region. RESULTS: We observed 99 genes with significant differential expression (false discovery rate [FDR] cut-off:10%) and an effect size of at least twofold. The subsequent promoter contact analyses revealed just one gene, ACSL5, which resides in the same topologically associating domain as TCF7L2. The generation of additional, smaller deletions (66 bp and 104 bp) comprising rs7903146 showed consistently reduced ACSL5 mRNA levels across all three deletions of up to 30-fold, with commensurate loss of acyl-CoA synthetase long-chain family member 5 (ACSL5) protein. Notably, the deletion of this single-nucleotide polymorphism region abolished significantly detectable chromatin contacts with the ACSL5 promoter. We went on to confirm that contacts between rs7903146 and the ACSL5 promoter regions were conserved in human colon tissue. ACSL5 encodes ACSL5, an enzyme with known roles in fatty acid metabolism. CONCLUSIONS/INTERPRETATION: This 'variant to gene mapping' effort implicates the genomic location harbouring rs7903146 as a regulatory region for ACSL5.

Patient reported outcomes among treatment modalities for prostate cancer.

INTRODUCTION: To characterize patient reported outcomes for urinary and sexual function using International Prostate Symptom Score (IPSS) and Sexual Health Inventory for Men (SHIM) comparing intensity modulated radiation therapy (IMRT), low dose rate brachytherapy (LDR), post-prostatectomy IMRT (PPRT), and radical prostatectomy (RP). MATERIALS AND METHODS: Patients treated for prostate cancer from 2001-2012 completed self-reported SHIM and IPSS surveys. Subgroups were created by baseline score. Mean change from baseline was determined at each time point for the cohort and subgroups. Statistical analysis was performed with generalized estimating equation method. Incontinence was not captured in the questionnaires. RESULTS: A total of 14,523 IPSS surveys from 3,515 men were evaluated. Patients treated with IMRT experienced a minimal decrease in IPSS score from baseline. PPRT scores did not differ from IMRT at any time point (range: +/- 3 points from baseline in IPSS score over 50 months). LDR had an initial IPSS rise (between 5-10 points on the IPSS over 1-9 months) versus IMRT but returned to comparable levels at 34 months. RP was associated with a lower IPSS versus IMRT. LDR had the largest rise from baseline, with return toward baseline. A total of 2,624 SHIM surveys from 857 men were evaluated. LDR and PPRT did not differ from IMRT at any time point (range: +/- 5 points from baseline in SHIM score for 36 months). RP experienced the largest decline from baseline (up to -7 points on SHIM score), at 3 to 7 months; RP had a larger early decrease in SHIM score versus IMRT between 3 and 22 months, after which there was no difference. CONCLUSIONS: IPSS and SHIM score patterns differed among treatment modalities. These data can be used to predict changes in urinary and sexual function over time based on modality and baseline score.

Use of Postprostatectomy Radiation Therapy at an NCI-Designated Comprehensive Cancer Center.

PURPOSE: Characterize use of postprostatectomy radiation (PPRT) for patients with prostate cancer at an NCI-designated comprehensive cancer center. METHODS: We queried our prospective prostate cancer database for patients treated with 60 to 68 Gy of radiation therapy (RT) to the prostate bed after prostatectomy from 2003 to 2011. Prostatectomy cases were obtained from billing records. Patients with an intact prostate treated with definitive RT served as a control for the change in volume of patients with prostate cancer treated in the department. Chi-square analysis assessed differences between adjuvant and salvage RT cohorts. Spearman correlation assessed yearly trends in prostate-specific antigen (PSA) level at the time of referral for RT. Linear regression models tested trends for number of PPRT cases, prostatectomies, and patients with intact prostate receiving radiation across years. RESULTS: PPRT was used to treat 475 men at Fox Chase Cancer Center from 2003 to 2011 (83 adjuvant and 392 salvage). Over time, an increased proportion of patients receiving RT to the prostate were treated with PPRT. No increase was seen in the proportion of patients treated with adjuvant RT compared with salvage RT (P=.5). Patients receiving adjuvant RT were younger, had higher pathologic Gleason score, pathologic T stage, and rates of positive margins than those receiving salvage RT. Pre-RT PSA values were inversely correlated with year (P=.005). The number of patients referred for salvage RT with a PSA of 0.5 ng/mL or less increased significantly from 7.9% in 2003 to 26.6% in 2011 (P=.002). CONCLUSIONS: A larger proportion of patients treated with RT for localized prostate cancer are now receiving PPRT. No increase was seen in the proportion of patients treated with adjuvant RT. Over time, patients with lower PSAs were referred for salvage RT.

Postmastectomy radiation therapy for T3N0: a SEER analysis.

BACKGROUND: There is conflicting evidence regarding the benefit of postmastectomy radiation therapy (PMRT) for pathologic stage T3N0M0 breast cancers. We analyzed data from the Surveillance, Epidemiology, and End Results (SEER) database to investigate the benefit of PMRT in this patient population. METHODS: We queried the SEER database for T3N0M0 breast cancer patients diagnosed from 2000 to 2010 who underwent modified radical mastectomy. We excluded males, patients with unknown radiation timing/type, other primary tumors, and survival <6 months. A total of 2525 patients were included in the analysis. We performed univariate and multivariate statistical analysis using chi-square tests, log-rank tests, and Cox proportional hazards regression. The primary endpoints were overall survival (OS) and cancer-specific survival (CSS). RESULTS: Of the 2525 patients identified, 1063 received PMRT. The median follow-up was 56 months (range, 6-131 months). On univariate analysis, PMRT improved OS (76.5% vs 61.8%, P<.01) and CSS (85.0% vs 82.4%, P<.01) at 8 years. The use of PMRT remained significant on multivariate analysis: PMRT improved OS (hazard ratio 0.63, P<.001) and CSS (hazard ratio 0.77, P = .045). Low tumor grade (P<.01) and marital status of "married" (P = .01) also was a predictor of improved CSS on multivariate analysis. CONCLUSIONS: PMRT was associated with significant improvements in both CSS and OS in patients with T3N0M0 breast cancers treated with modified radical mastectomy from 2000 to 2010. PMRT should be strongly considered in T3N0M0 patients. Postmastectomy radiation therapy is associated with significant improvements in overall and cause-specific survival in patients with T3N0M0 breast cancers treated with modified radical mastectomy from 2000 to 2010 in the SEER database. Postmastectomy radiation therapy should be strongly considered for patients who have T3N0M0 tumors.

A ChIP-seq-defined genome-wide map of MEF2C binding reveals inflammatory pathways associated with its role in bone density determination.

Genome-wide association studies (GWAS) have demonstrated that genetic variation at the MADS box transcription enhancer factor 2, polypeptide C (MEF2C) locus is robustly associated with bone mineral density, primarily at the femoral neck. MEF2C is a transcription factor known to operate via the Wnt signaling pathway. Our hypothesis was that MEF2C regulates the expression of a set of molecular pathways critical to skeletal function. Drawing on our laboratory and bioinformatic experience with ChIP-seq, we analyzed ChIP-seq data for MEF2C available via the ENCODE project to gain insight in to its global genomic binding pattern. We aligned the ChIP-seq data generated for GM12878 (an established lymphoblastoid cell line) and, using the analysis package HOMER, a total of 17,611 binding sites corresponding to 8,118 known genes were observed. We then performed a pathway analysis of the gene list using Ingenuity. At 5 kb, the gene list yielded 'EIF2 Signaling' as the most significant annotation, with a P value of 5.01 × 10(-26). Moving further out, this category remained the top pathway at 50 and 100 kb, then dropped to just second place at 500 kb and beyond by 'Molecular Mechanisms of Cancer'. In addition, at 50 kb and beyond 'RANK Signaling in Osteoclasts' was a consistent feature and resonates with the main general finding from GWAS of bone density. We also observed that MEF2C binding sites were significantly enriched primarily near inflammation associated genes identified from GWAS; indeed, a similar enrichment for inflammation genes has been reported previously using a similar approach for the vitamin D receptor, an established key regulator of bone turnover. Our analyses point to known connective tissue and skeletal processes but also provide novel insights in to networks involved in skeletal regulation. The fact that a specific GWAS category is enriched points to a possible role of inflammation through which it impacts bone mineral density.

Absolute lymphocyte count: a potential prognostic factor for Merkel cell carcinoma.

BACKGROUND: Absolute lymphocyte count (ALC) is a laboratory value commonly obtained during workup of patients with Merkel cell carcinoma (MCC). OBJECTIVE: We report the prognostic impact of ALC as a surrogate of immune status in MCC. METHODS: A complete blood cell count was available for 64 patients with MCC in the month before definitive surgery, chemotherapy, or radiation. Statistical analysis was performed with classification and regression tree analysis, log rank test, and Cox model. RESULTS: Median overall survival (OS) for the cohort was 97 months. Median OS for patients with an ALC less than 1.1 k/mm(3) was 18.8 versus 110.1 months for those with ALC greater than or equal to 1.1 k/mm(3) (P = .002, hazard ratio 0.29). Multivariate analysis of OS controlling for ALC, sex, stage, adjuvant chemotherapy, hematologic malignancy, and immunosuppression demonstrated ALC as a prognostic factor (P = .03). Disease-free survival at 36 months for ALC less than 1.1 k/mm(3) was 26.9% versus 64.4% for those with ALC greater than or equal to 1.1 k/mm(3) (P = .01). ALC was not a significant predictor for disease-free survival on multivariate analysis (P = .12). LIMITATIONS: This is a single-institution retrospective data set. CONCLUSION: ALC is associated with OS but not disease-free survival in MCC using a threshold of less than 1.1 k/mm(3). This test may provide additional prognostic information for patients with MCC.

Projecting genetic associations through gene expression patterns highlights disease etiology and drug mechanisms.

Genes act in concert with each other in specific contexts to perform their functions. Determining how these genes influence complex traits requires a mechanistic understanding of expression regulation across different conditions. It has been shown that this insight is critical for developing new therapies. Transcriptome-wide association studies have helped uncover the role of individual genes in disease-relevant mechanisms. However, modern models of the architecture of complex traits predict that gene-gene interactions play a crucial role in disease origin and progression. Here we introduce PhenoPLIER, a computational approach that maps gene-trait associations and pharmacological perturbation data into a common latent representation for a joint analysis. This representation is based on modules of genes with similar expression patterns across the same conditions. We observe that diseases are significantly associated with gene modules expressed in relevant cell types, and our approach is accurate in predicting known drug-disease pairs and inferring mechanisms of action. Furthermore, using a CRISPR screen to analyze lipid regulation, we find that functionally important players lack associations but are prioritized in trait-associated modules by PhenoPLIER. By incorporating groups of co-expressed genes, PhenoPLIER can contextualize genetic associations and reveal potential targets missed by single-gene strategies.

Integrative Genomic Analyses Identify LncRNA Regulatory Networks across Pediatric Leukemias and Solid Tumors.

Long noncoding RNAs (lncRNA) play an important role in gene regulation and contribute to tumorigenesis. While pan-cancer studies of lncRNA expression have been performed for adult malignancies, the lncRNA landscape across pediatric cancers remains largely uncharted. Here, we curated RNA sequencing data for 1,044 pediatric leukemia and extracranial solid tumors and integrated paired tumor whole genome sequencing and epigenetic data in relevant cell line models to explore lncRNA expression, regulation, and association with cancer. A total of 2,657 lncRNAs were robustly expressed across six pediatric cancers, including 1,142 exhibiting histotype-elevated expression. DNA copy number alterations contributed to lncRNA dysregulation at a proportion comparable to protein coding genes. Application of a multidimensional framework to identify and prioritize lncRNAs impacting gene networks revealed that lncRNAs dysregulated in pediatric cancer are associated with proliferation, metabolism, and DNA damage hallmarks. Analysis of upstream regulation via cell type-specific transcription factors further implicated distinct histotype-elevated and developmental lncRNAs. Integration of these analyses prioritized lncRNAs for experimental validation, and silencing of TBX2-AS1, the top-prioritized neuroblastoma-specific lncRNA, resulted in significant growth inhibition of neuroblastoma cells, confirming the computational predictions. Taken together, these data provide a comprehensive characterization of lncRNA regulation and function in pediatric cancers and pave the way for future mechanistic studies. SIGNIFICANCE: Comprehensive characterization of lncRNAs in pediatric cancer leads to the identification of highly expressed lncRNAs across childhood cancers, annotation of lncRNAs showing histotype-specific elevated expression, and prediction of lncRNA gene regulatory networks.

Variant-to-gene-mapping analyses reveal a role for pancreatic islet cells in conferring genetic susceptibility to sleep-related traits.

We investigated the potential role of sleep-trait associated genetic loci in conferring a degree of their effect via pancreatic α- and ß-cells, given that both sleep disturbances and metabolic disorders, including type 2 diabetes and obesity, involve polygenic contributions and complex interactions. We determined genetic commonalities between sleep and metabolic disorders, conducting linkage disequilibrium genetic correlation analyses with publicly available GWAS summary statistics. Then we investigated possible enrichment of sleep-trait associated SNPs in promoter-interacting open chromatin regions within α- and ß-cells, intersecting public GWAS reports with our own ATAC-seq and high-resolution promoter-focused Capture C data generated from both sorted human α-cells and an established human beta-cell line (EndoC-ßH1). Finally, we identified putative effector genes physically interacting with sleep-trait associated variants in α- and EndoC-ßH1cells running variant-to-gene mapping and establish pathways in which these genes are significantly involved. We observed that insomnia, short and long sleep-but not morningness-were significantly correlated with type 2 diabetes, obesity and other metabolic traits. Both the EndoC-ßH1 and α-cells were enriched for insomnia loci (p = .01; p = .0076), short sleep loci (p = .017; p = .022) and morningness loci (p = 2.2 × 10-7; p = .0016), while the α-cells were also enriched for long sleep loci (p = .034). Utilizing our promoter contact data, we identified 63 putative effector genes in EndoC-ßH1 and 76 putative effector genes in α-cells, with these genes showing significant enrichment for organonitrogen and organophosphate biosynthesis, phosphatidylinositol and phosphorylation, intracellular transport and signaling, stress responses and cell differentiation. Our data suggest that a subset of sleep-related loci confer their effects via cells in pancreatic islets.

Implicating effector genes at COVID-19 GWAS loci using promoter-focused Capture-C in disease-relevant immune cell types.

BACKGROUND: SARS-CoV-2 infection results in a broad spectrum of COVID-19 disease, from mild or no symptoms to hospitalization and death. COVID-19 disease severity has been associated with some pre-existing conditions and the magnitude of the adaptive immune response to SARS-CoV-2, and a recent genome-wide association study (GWAS) of the risk of critical illness revealed a significant genetic component. To gain insight into how human genetic variation attenuates or exacerbates disease following SARS-CoV-2 infection, we implicated putatively functional COVID risk variants in the cis-regulatory landscapes of human immune cell types with established roles in disease severity and used high-resolution chromatin conformation capture to map these disease-associated elements to their effector genes. RESULTS: This functional genomic approach implicates 16 genes involved in viral replication, the interferon response, and inflammation. Several of these genes (PAXBP1, IFNAR2, OAS1, OAS3, TNFAIP8L1, GART) were differentially expressed in immune cells from patients with severe versus moderate COVID-19 disease, and we demonstrate a previously unappreciated role for GART in T cell-dependent antibody-producing B cell differentiation in a human tonsillar organoid model. CONCLUSIONS: This study offers immunogenetic insight into the basis of COVID-19 disease severity and implicates new targets for therapeutics that limit SARS-CoV-2 infection and its resultant life-threatening inflammation.

3D promoter architecture re-organization during iPSC-derived neuronal cell differentiation implicates target genes for neurodevelopmental disorders.

Neurodevelopmental disorders are thought to arise from interrupted development of the brain at an early age. Genome-wide association studies (GWAS) have identified hundreds of loci associated with susceptibility to neurodevelopmental disorders; however, which noncoding variants regulate which genes at these loci is often unclear. To implicate neuronal GWAS effector genes, we performed an integrated analysis of transcriptomics, epigenomics and chromatin conformation changes during the development from Induced pluripotent stem cell-derived neuronal progenitor cells (NPCs) into neurons using a combination of high-resolution promoter-focused Capture-C, ATAC-seq and RNA-seq. We observed that gene expression changes during the NPC-to-neuron transition were highly dependent on both promoter accessibility changes and long-range interactions which connect distal cis-regulatory elements (enhancer or silencers) to developmental-stage-specific genes. These genome-scale promoter-cis-regulatory-element atlases implicated 454 neurodevelopmental disorder-associated, putative causal variants mapping to 600 distal targets. These putative effector genes were significantly enriched for pathways involved in the regulation of neuronal development and chromatin organization, with 27 % expressed in a stage-specific manner. The intersection of open chromatin and chromatin conformation revealed development-stage-specific gene regulatory architectures during neuronal differentiation, providing a rich resource to aid characterization of the genetic and developmental basis of neurodevelopmental disorders.

CRISPR-Cas9-Mediated Genome Editing Confirms EPDR1 as an Effector Gene at the BMD GWAS-Implicated 'STARD3NL' Locus.

Genome-wide-association studies (GWASs) have discovered genetic signals robustly associated with BMD, but typically not the precise localization of effector genes. By intersecting genome-wide promoter-focused Capture C and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data generated in human mesenchymal progenitor cell (hMSC)-derived osteoblasts, consistent contacts were previously predicted between the EPDR1 promoter and multiple BMD-associated candidate causal variants at the 'STARD3NL' locus. RNAi knockdown of EPDR1 expression in hMSC-derived osteoblasts was shown to lead to inhibition of osteoblastogenesis. To fully characterize the physical connection between these putative noncoding causal variants at this locus and the EPDR1 gene, clustered regularly interspaced short-palindromic repeat Cas9 endonuclease (CRISPR-Cas9) genome editing was conducted in hFOB1.19 cells across the single open-chromatin region harboring candidates for the underlying causal variant, rs1524068, rs6975644, and rs940347, all in close proximity to each other. RT-qPCR and immunoblotting revealed dramatic and consistent downregulation of EPDR1 specifically in the edited differentiated osteoblast cells. Consistent with EPDR1 expression changes, alkaline phosphatase staining was also markedly reduced in the edited differentiated cells. Collectively, CRISPR-Cas9 genome editing in the hFOB1.19 cell model supports previous observations, where this regulatory region harboring GWAS-implicated variation operates through direct long-distance physical contact, further implicating a key role for EPDR1 in osteoblastogenesis and BMD determination. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.