FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases.

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.

Spin Coupling Effect on Geometry-Dependent X-Ray Absorption of Diradicals.

We theoretically investigate the influence of diradical electron spin coupling on the time-resolved X-ray absorption spectra of the photochemical ring opening of furanone. We predict geometry-dependent carbon K-edge signals involving transitions from core orbitals to both singly and unoccupied molecular orbitals. The most obvious features of the ring opening come from the carbon atom directly involved in the bond breaking through its transition to both the newly formed singly occupied and the available lowest unoccupied molecular orbitals (SOMO and LUMO, respectively). In addition to this primary feature, the singlet spin coupling of four unpaired electrons that arises in the core-to-LUMO states creates additional geometry dependence in some spectral features with both oscillator strengths and relative excitation energies varying observably as a function of the ring opening. We attribute this behavior to a spin-occupancy-induced selection rule, which occurs when singlet spin coupling is enforced in the diradical state. Notably, one of these geometry-sensitive core-to-LUMO transitions excites core electrons from a backbone carbon not involved in the bond breaking, providing a novel nonlocal X-ray probe of chemical dynamics arising from electron spin coupling.

Femtosecond symmetry breaking and coherent relaxation of methane cations via x-ray spectroscopy.

Understanding the relaxation pathways of photoexcited molecules is essential to gain atomistic-level insight into photochemistry. We performed a time-resolved study of ultrafast molecular symmetry breaking through geometric relaxation (Jahn-Teller distortion) on the methane cation. Attosecond transient absorption spectroscopy with soft x-rays at the carbon K-edge revealed that the distortion occurred within 10 ± 2 femtoseconds after few-femtosecond strong-field ionization of methane. The distortion activated coherent oscillations in the asymmetric scissoring vibrational mode of the symmetry-broken cation, which were detected in the x-ray signal. These oscillations were damped within 58 ± 13 femtoseconds because vibrational coherence was lost with the energy redistributing into lower-frequency vibrational modes. This study completely reconstructs the molecular relaxation dynamics of this prototypical example and opens avenues for exploring complex systems.

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron-Withdrawing Groups.

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the 1ππ* (S2) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron-withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the 1nπ* (S1) and 3ππ* (T1) states on a 1.6 ± 0.4 ps time scale, which is similar to the 1.5 ps time scale earlier observed for acetylacetone [ J. Am. Chem. Soc. 2017, 139, 16576-16583, DOI: 10.1021/jacs.7b07532]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

Correction: Jahn-Teller distortion and dissociation of CCl4 + by transient X-ray spectroscopy simultaneously at the carbon K- and chlorine L-edge.

[This corrects the article DOI: 10.1039/D2SC02402K.].

Jahn-Teller distortion and dissociation of CCl4 + by transient X-ray spectroscopy simultaneously at the carbon K- and chlorine L-edge.

X-ray Transient Absorption Spectroscopy (XTAS) and theoretical calculations are used to study CCl4 + prepared by 800 nm strong-field ionization. XTAS simultaneously probes atoms at the carbon K-edge (280-300 eV) and chlorine L-edge (195-220 eV). Comparison of experiment to X-ray spectra computed by orbital-optimized density functional theory (OO-DFT) indicates that after ionization, CCl4 + undergoes symmetry breaking driven by Jahn-Teller distortion away from the initial tetrahedral structure (Td) in 6 ± 2 fs. The resultant symmetry-broken covalently bonded form subsequently separates to a noncovalently bound complex between CCl3 + and Cl over 90 ± 10 fs, which is again predicted by theory. Finally, after more than 800 fs, L-edge signals for atomic Cl are observed, indicating dissociation to free CCl3 + and Cl. The results for Jahn-Teller distortion to the symmetry-broken form of CCl4 + and formation of the Cl-CCl+ 3 complex characterize previously unobserved new species along the route to dissociation.

Associations Between Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants and Risk of Coronavirus Disease 2019 (COVID-19) Hospitalization Among Confirmed Cases in Washington State: A Retrospective Cohort Study.

BACKGROUND: The coronavirus disease 2019 (COVID-19) pandemic is dominated by variant viruses; the resulting impact on disease severity remains unclear. Using a retrospective cohort study, we assessed the hospitalization risk following infection with 7 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants. METHODS: Our study includes individuals with positive SARS-CoV-2 reverse transcription polymerase chain reaction (RT-PCR) in the Washington Disease Reporting System with available viral genome data, from 1 December 2020 to 14 January 2022. The analysis was restricted to cases with specimens collected through sentinel surveillance. Using a Cox proportional hazards model with mixed effects, we estimated hazard ratios (HR) for hospitalization risk following infection with a variant, adjusting for age, sex, calendar week, and vaccination. RESULTS: In total, 58 848 cases were sequenced through sentinel surveillance, of which 1705 (2.9%) were hospitalized due to COVID-19. Higher hospitalization risk was found for infections with Gamma (HR 3.20, 95% confidence interval [CI] 2.40-4.26), Beta (HR 2.85, 95% CI 1.56-5.23), Delta (HR 2.28 95% CI 1.56-3.34), or Alpha (HR 1.64, 95% CI 1.29-2.07) compared to infections with ancestral lineages; Omicron (HR 0.92, 95% CI .56-1.52) showed no significant difference in risk. Following Alpha, Gamma, or Delta infection, unvaccinated patients show higher hospitalization risk, while vaccinated patients show no significant difference in risk, both compared to unvaccinated, ancestral lineage cases. Hospitalization risk following Omicron infection is lower with vaccination. CONCLUSIONS: Infection with Alpha, Gamma, or Delta results in a higher hospitalization risk, with vaccination attenuating that risk. Our findings support hospital preparedness, vaccination, and genomic surveillance.

Differences in nanoscale organization of regulatory active and inactive human chromatin.

Methodological advances in conformation capture techniques have fundamentally changed our understanding of chromatin architecture. However, the nanoscale organization of chromatin and its cell-to-cell variance are less studied. Analyzing genome-wide data from 733 human cell and tissue samples, we identified 2 prototypical regions that exhibit high or absent hypersensitivity to deoxyribonuclease I, respectively. These regulatory active or inactive regions were examined in the lymphoblast cell line K562 by using high-throughput super-resolution microscopy. In both regions, we systematically measured the physical distance of 2 fluorescence in situ hybridization spots spaced by only 5 kb of DNA. Unexpectedly, the resulting distance distributions range from very compact to almost elongated configurations of more than 200-nm length for both the active and inactive regions. Monte Carlo simulations of a coarse-grained model of these chromatin regions based on published data of nucleosome occupancy in K562 cells were performed to understand the underlying mechanisms. There was no parameter set for the simulation model that can explain the microscopically measured distance distributions. Obviously, the chromatin state given by the strength of internucleosomal interaction, nucleosome occupancy, or amount of histone H1 differs from cell to cell, which results in the observed broad distance distributions. This large variability was not expected, especially in inactive regions. The results for the mechanisms for different distance distributions on this scale are important for understanding the contacts that mediate gene regulation. Microscopic measurements show that the inactive region investigated here is expected to be embedded in a more compact chromatin environment. The simulation results of this region require an increase in the strength of internucleosomal interactions. It may be speculated that the higher density of chromatin is caused by the increased internucleosomal interaction strength.

Integrative analysis of 3604 GWAS reveals multiple novel cell type-specific regulatory associations.

BACKGROUND: Genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) are known to preferentially co-locate to active regulatory elements in tissues and cell types relevant to disease aetiology. Further characterisation of associated cell type-specific regulation can broaden our understanding of how GWAS signals may contribute to disease risk. RESULTS: To gain insight into potential functional mechanisms underlying GWAS associations, we developed FORGE2 ( https://forge2.altiusinstitute.org/ ), which is an updated version of the FORGE web tool. FORGE2 uses an expanded atlas of cell type-specific regulatory element annotations, including DNase I hotspots, five histone mark categories and 15 hidden Markov model (HMM) chromatin states, to identify tissue- and cell type-specific signals. An analysis of 3,604 GWAS from the NHGRI-EBI GWAS catalogue yielded at least one significant disease/trait-tissue association for 2,057 GWAS, including > 400 associations specific to epigenomic marks in immune tissues and cell types, > 30 associations specific to heart tissue, and > 60 associations specific to brain tissue, highlighting the key potential of tissue- and cell type-specific regulatory elements. Importantly, we demonstrate that FORGE2 analysis can separate previously observed accessible chromatin enrichments into different chromatin states, such as enhancers or active transcription start sites, providing a greater understanding of underlying regulatory mechanisms. Interestingly, tissue-specific enrichments for repressive chromatin states and histone marks were also detected, suggesting a role for tissue-specific repressed regions in GWAS-mediated disease aetiology. CONCLUSION: In summary, we demonstrate that FORGE2 has the potential to uncover previously unreported disease-tissue associations and identify new candidate mechanisms. FORGE2 is a transparent, user-friendly web tool for the integrative analysis of loci discovered from GWAS.

Associations between SARS-CoV-2 variants and risk of COVID-19 hospitalization among confirmed cases in Washington State: a retrospective cohort study.

BACKGROUND: The COVID-19 pandemic is dominated by variant viruses; the resulting impact on disease severity remains unclear. Using a retrospective cohort study, we assessed the hospitalization risk following infection with seven SARS-CoV-2 variants. METHODS: Our study includes individuals with positive SARS-CoV-2 RT-PCR in the Washington Disease Reporting System with available viral genome data, from December 1, 2020 to January 14, 2022. The analysis was restricted to cases with specimens collected through sentinel surveillance. Using a Cox proportional hazards model with mixed effects, we estimated hazard ratios (HR) for hospitalization risk following infection with a variant, adjusting for age, sex, calendar week, and vaccination. FINDINGS: 58,848 cases were sequenced through sentinel surveillance, of which 1705 (2.9%) were hospitalized due to COVID-19. Higher hospitalization risk was found for infections with Gamma (HR 3.20, 95%CI 2.40-4.26), Beta (HR 2.85, 95%CI 1.56-5.23), Delta (HR 2.28 95%CI 1.56-3.34) or Alpha (HR 1.64, 95%CI 1.29-2.07) compared to infections with ancestral lineages; Omicron (HR 0.92, 95%CI 0.56-1.52) showed no significant difference in risk. Following Alpha, Gamma, or Delta infection, unvaccinated patients show higher hospitalization risk, while vaccinated patients show no significant difference in risk, both compared to unvaccinated, ancestral lineage cases. Hospitalization risk following Omicron infection is lower with vaccination. CONCLUSION: Infection with Alpha, Gamma, or Delta results in a higher hospitalization risk, with vaccination attenuating that risk. Our findings support hospital preparedness, vaccination, and genomic surveillance. SUMMARY: Hospitalization risk following infection with SARS-CoV-2 variant remains unclear. We find a higher hospitalization risk in cases infected with Alpha, Beta, Gamma, and Delta, but not Omicron, with vaccination lowering risk. Our findings support hospital preparedness, vaccination, and genomic surveillance.

X-ray transient absorption reveals the 1Au (nπ*) state of pyrazine in electronic relaxation.

Electronic relaxation in organic chromophores often proceeds via states not directly accessible by photoexcitation. We report on the photoinduced dynamics of pyrazine that involves such states, excited by a 267 nm laser and probed with X-ray transient absorption spectroscopy in a table-top setup. In addition to the previously characterized 1B2u (ππ*) (S2) and 1B3u (nπ*) (S1) states, the participation of the optically dark 1Au (nπ*) state is assigned by a combination of experimental X-ray core-to-valence spectroscopy, electronic structure calculations, nonadiabatic dynamics simulations, and X-ray spectral computations. Despite 1Au (nπ*) and 1B3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The 1Au (nπ*) state is populated in 200 ± 50 femtoseconds after electronic excitation and plays a key role in the relaxation of pyrazine to the ground state.

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations.

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization [D. Hait and M. Head-Gordon, J. Chem. Theory Comput. 16, 1699 (2020)] algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO+, and the N K-edge of NO2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. We also present general guidelines for computing core-excited states from orbital optimized DFT.

Index and biological spectrum of human DNase I hypersensitive sites.

DNase I hypersensitive sites (DHSs) are generic markers of regulatory DNA1-5 and contain genetic variations associated with diseases and phenotypic traits6-8. We created high-resolution maps of DHSs from 733 human biosamples encompassing 438 cell and tissue types and states, and integrated these to delineate and numerically index approximately 3.6 million DHSs within the human genome sequence, providing a common coordinate system for regulatory DNA. Here we show that these maps highly resolve the cis-regulatory compartment of the human genome, which encodes unexpectedly diverse cell- and tissue-selective regulatory programs at very high density. These programs can be captured comprehensively by a simple vocabulary that enables the assignment to each DHS of a regulatory barcode that encapsulates its tissue manifestations, and global annotation of protein-coding and non-coding RNA genes in a manner orthogonal to gene expression. Finally, we show that sharply resolved DHSs markedly enhance the genetic association and heritability signals of diseases and traits. Rather than being confined to a small number of distal elements or promoters, we find that genetic signals converge on congruently regulated sets of DHSs that decorate entire gene bodies. Together, our results create a universal, extensible coordinate system and vocabulary for human regulatory DNA marked by DHSs, and provide a new global perspective on the architecture of human gene regulation.

Global reference mapping of human transcription factor footprints.

Combinatorial binding of transcription factors to regulatory DNA underpins gene regulation in all organisms. Genetic variation in regulatory regions has been connected with diseases and diverse phenotypic traits1, but it remains challenging to distinguish variants that affect regulatory function2. Genomic DNase I footprinting enables the quantitative, nucleotide-resolution delineation of sites of transcription factor occupancy within native chromatin3-6. However, only a small fraction of such sites have been precisely resolved on the human genome sequence6. Here, to enable comprehensive mapping of transcription factor footprints, we produced high-density DNase I cleavage maps from 243 human cell and tissue types and states and integrated these data to delineate about 4.5 million compact genomic elements that encode transcription factor occupancy at nucleotide resolution. We map the fine-scale structure within about 1.6 million DNase I-hypersensitive sites and show that the overwhelming majority are populated by well-spaced sites of single transcription factor-DNA interaction. Cell-context-dependent cis-regulation is chiefly executed by wholesale modulation of accessibility at regulatory DNA rather than by differential transcription factor occupancy within accessible elements. We also show that the enrichment of genetic variants associated with diseases or phenotypic traits in regulatory regions1,7 is almost entirely attributable to variants within footprints, and that functional variants that affect transcription factor occupancy are nearly evenly partitioned between loss- and gain-of-function alleles. Unexpectedly, we find increased density of human genetic variation within transcription factor footprints, revealing an unappreciated driver of cis-regulatory evolution. Our results provide a framework for both global and nucleotide-precision analyses of gene regulatory mechanisms and functional genetic variation.

Single-molecule regulatory architectures captured by chromatin fiber sequencing.

Gene regulation is chiefly determined at the level of individual linear chromatin molecules, yet our current understanding of cis-regulatory architectures derives from fragmented sampling of large numbers of disparate molecules. We developed an approach for precisely stenciling the structure of individual chromatin fibers onto their composite DNA templates using nonspecific DNA N6-adenine methyltransferases. Single-molecule long-read sequencing of chromatin stencils enabled nucleotide-resolution readout of the primary architecture of multikilobase chromatin fibers (Fiber-seq). Fiber-seq exposed widespread plasticity in the linear organization of individual chromatin fibers and illuminated principles guiding regulatory DNA actuation, the coordinated actuation of neighboring regulatory elements, single-molecule nucleosome positioning, and single-molecule transcription factor occupancy. Our approach and results open new vistas on the primary architecture of gene regulation.

Integrated epigenomic profiling reveals endogenous retrovirus reactivation in renal cell carcinoma.

BACKGROUND: Transcriptional dysregulation drives cancer formation but the underlying mechanisms are still poorly understood. Renal cell carcinoma (RCC) is the most common malignant kidney tumor which canonically activates the hypoxia-inducible transcription factor (HIF) pathway. Despite intensive study, novel therapeutic strategies to target RCC have been difficult to develop. Since the RCC epigenome is relatively understudied, we sought to elucidate key mechanisms underpinning the tumor phenotype and its clinical behavior. METHODS: We performed genome-wide chromatin accessibility (DNase-seq) and transcriptome profiling (RNA-seq) on paired tumor/normal samples from 3 patients undergoing nephrectomy for removal of RCC. We incorporated publicly available data on HIF binding (ChIP-seq) in a RCC cell line. We performed integrated analyses of these high-resolution, genome-scale datasets together with larger transcriptomic data available through The Cancer Genome Atlas (TCGA). FINDINGS: Though HIF transcription factors play a cardinal role in RCC oncogenesis, we found that numerous transcription factors with a RCC-selective expression pattern also demonstrated evidence of HIF binding near their gene body. Examination of chromatin accessibility profiles revealed that some of these transcription factors influenced the tumor's regulatory landscape, notably the stem cell transcription factor POU5F1 (OCT4). Elevated POU5F1 transcript levels were correlated with advanced tumor stage and poorer overall survival in RCC patients. Unexpectedly, we discovered a HIF-pathway-responsive promoter embedded within a endogenous retroviral long terminal repeat (LTR) element at the transcriptional start site of the PSOR1C3 long non-coding RNA gene upstream of POU5F1. RNA transcripts are induced from this promoter and read through PSOR1C3 into POU5F1 producing a novel POU5F1 transcript isoform. Rather than being unique to the POU5F1 locus, we found that HIF binds to several other transcriptionally active LTR elements genome-wide correlating with broad gene expression changes in RCC. INTERPRETATION: Integrated transcriptomic and epigenomic analysis of matched tumor and normal tissues from even a small number of primary patient samples revealed remarkably convergent shared regulatory landscapes. Several transcription factors appear to act downstream of HIF including the potent stem cell transcription factor POU5F1. Dysregulated expression of POU5F1 is part of a larger pattern of gene expression changes in RCC that may be induced by HIF-dependent reactivation of dormant promoters embedded within endogenous retroviral LTRs.

Clinical outcomes with a systolic blood pressure lower than 120 mmHg in older patients with high disease burden.

BACKGROUND: The target blood pressure in older patients is controversial. Recent studies provided clinical evidence supporting a target systolic blood pressure <120 mmHg in patients >50 years at high risk of cardiovascular events. METHODS: Retrospective study of 380 consecutive patients ≥60 years with stages 1-5 pre-dialysis chronic kidney disease seen between January 2013 and November 2015. The outcomes of a systolic blood pressure <120 mmHg in older patients with chronic kidney disease and multiple comorbidities were analyzed. RESULTS: Sixty-eight patients had a systolic blood pressure <120 mmHg, 312 patients had a systolic blood pressure ≥120 mmHg. Forty-three patients died during the follow up (11.3%). Patients with a systolic blood pressure <120 mmHg had a higher risk of death: 21 (30.9%) vs 22 (7%). Primary cause of death: Cardiovascular: 11 (25.6%), infectious 9 (20.9%), cancer 5 (11.6%), renal failure 6 (13.9%), COPD/pulmonary fibrosis 2 (4.6%), end stage liver disease 3 (6.9%), traumatic brain injury 1 (2.3%), gastrointestinal hemorrhage 4 (9.3%), complications of diabetes 1 (2.3%), unknown 1 (2.3%). After adjusting for confounding factors, a systolic blood pressure <120 mmHg remained associated with increased mortality. There was a trend to more cardiovascular outcomes in those with a lower blood pressure. CONCLUSIONS: A systolic blood pressure below 120 mmHg in older patients with high disease burden was associated with adverse outcomes. Individualization of blood pressure therapy to each specific patient is warranted.

Rare Variants of the Serotonin Transporter Are Associated With Psychiatric Comorbidity in Irritable Bowel Syndrome.

Alterations in serotonin signaling are suspected in the pathophysiology of irritable bowel syndrome (IBS). By modulating the extracellular reuptake of serotonin, the serotonin transporter (SERT) acts as a key regulator of the bioavailability of serotonin. This study is the first to investigate the impact of rare SERT variants (i.e., those with a minor allele frequency of < 1%) on the risk for IBS, gastrointestinal (GI) symptom level, response to cognitive-behavioral treatment, and psychiatric comorbidity. We sequenced a 0.19 megabase chromosomal stretch containing the SERT gene and surrounding regions in a community sample of 304 IBS patients and 83 controls. We found no significant associations between rare variants in and around the SERT gene and IBS risk, GI symptom profile, or response to treatment. We found preliminary evidence, however, that IBS subjects with a history of either depression or anxiety were significantly more likely to carry multiple rare likely functional variant alleles than IBS patients without psychiatric comorbidity.

Large-scale identification of sequence variants influencing human transcription factor occupancy in vivo.

The function of human regulatory regions depends exquisitely on their local genomic environment and on cellular context, complicating experimental analysis of common disease- and trait-associated variants that localize within regulatory DNA. We use allelically resolved genomic DNase I footprinting data encompassing 166 individuals and 114 cell types to identify >60,000 common variants that directly influence transcription factor occupancy and regulatory DNA accessibility in vivo. The unprecedented scale of these data enables systematic analysis of the impact of sequence variation on transcription factor occupancy in vivo. We leverage this analysis to develop accurate models of variation affecting the recognition sites for diverse transcription factors and apply these models to discriminate nearly 500,000 common regulatory variants likely to affect transcription factor occupancy across the human genome. The approach and results provide a new foundation for the analysis and interpretation of noncoding variation in complete human genomes and for systems-level investigation of disease-associated variants.