Background:Lung adenocarcinoma (LUAD) is clinically and molecularly defined by oncogenic driver mutations, identification of which has led to the development of driver-targeted therapies and substantial improvements in prognosis for subsets of LUAD patients. Recent studies assessing clinical outcomes in the context of multigenic alterations have identified secondary mutations that might explain differential responses to targeted therapies, chemotherapies and immunotherapies. Genetic inactivation or loss of SMARCA4 , which frequently co-occurs with loss-of-function mutations in STK11 and KEAP1 , is especially predictive of poor prognosis and shorter overall survival in LUAD patients, regardless of driver status. We sought to examine the clinical and functional associations of SMARCA4 deficiency in LUAD, with or without co-associated STK11/KEAP1 loss-of-function. Methods:We examined correlation between SMARCA4 loss, gene expression and prognosis through genomic and transcriptomic profiling of clinically annotated LUAD samples. We generated isogenic cell line models with genetic knockouts of SMARCA4 with or without concomitant STK11 and KEAP1 to profile mutationally-defined genotypes of interest in vitro and in vivo . Lastly, we interrogated the functional dependency of SMARCA4/STK11/KEAP1 triple mutant models on TGF-β signaling to assess its potential as a therapeutic target. Results:SMARCA4/STK11/KEAP1 triple mutant LUAD is associated with poor survival and high frequency of multisite metastasis. SMARCA4/STK11/KEAP1 triple knockout models showed enhanced migration and invasion in vitro , and diversified organotropism in an in vivo intracardiac xenograft metastasis assay. RNA-Seq and DNaseI-Seq of these in vitro models and clinical samples identified upregulation of TGF-β signaling and EMT gene expression signatures, and corresponding changes in chromatin accessibility, in SMARCA4/STK11/KEAP1 triple mutant LUAD. Conclusions:We identify SMARCA4/STK11/KEAP1 triple mutant LUAD as a prognostically significant disease subset and nominate TGF-β signaling as a potential therapeutic target.

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.

Abstract 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.

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 traits 1 , but it remains challenging to distinguish variants that affect regulatory function 2 . Genomic DNase I footprinting enables the quantitative, nucleotide-resolution delineation of sites of transcription factor occupancy within native chromatin 3 – 6 . However, only a small fraction of such sites have been precisely resolved on the human genome sequence 6 . 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 regions 1 , 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.

The human genome encodes millions of regulatory elements, of which only a small fraction are active within a given cell type. Little is known about the global impact of chromatin remodelers on regulatory DNA landscapes and how this translates to gene expression. We use precision genome engineering to reawaken homozygously inactivated SMARCA4, a central ATPase of the human SWI/SNF chromatin remodeling complex, in lung adenocarcinoma cells. Here, we combine DNase I hypersensitivity, histone modification, and transcriptional profiling to show that SMARCA4 dramatically increases both the number and magnitude of accessible chromatin sites genome-wide, chiefly by unmasking sites of low regulatory factor occupancy. By contrast, transcriptional changes are concentrated within well-demarcated remodeling domains wherein expression of specific genes is gated by both distal element activation and promoter chromatin configuration. Our results provide a perspective on how global chromatin remodeling activity is translated to gene expression via regulatory DNA.

Early mammalian development is orchestrated by genome-encoded regulatory elements populated by a changing complement of regulatory factors, creating a dynamic chromatin landscape. To define the spatiotemporal organization of regulatory DNA landscapes during mouse development and maturation, we generated nucleotide-resolution DNA accessibility maps from 15 tissues sampled at 9 intervals spanning post-conception day 9.5 through early adult, and integrated these with 41 adult-stage DNase-seq profiles to create a global atlas of mouse regulatory DNA. Collectively, we delineated >1.8 million DNase I hypersensitive sites (DHSs), with the vast majority displaying temporal and tissue-selective patterning. Here we show that tissue regulatory DNA compartments show sharp embryonic-to-fetal transitions characterized by wholesale turnover of DHSs and progressive domination by a diminishing number of transcription factors. We show further that aligning mouse and human fetal development on a regulatory axis exposes disease-associated variation enriched in early intervals lacking human samples. Our results provide an expansive new resource for decoding mammalian developmental regulatory programs. ### Competing Interest Statement The authors have declared no competing interest.

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 traits, yet it remains challenging to distinguish variants that impact regulatory function. Genomic DNase I footprinting enables quantitative, nucleotide-resolution delineation of sites of transcription factor occupancy within native chromatin. However, to date only a small fraction of such sites have been precisely resolved on the human genome sequence. 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 at nucleotide resolution ~4.5 million compact genomic elements encoding transcription factor occupancy. We map the fine-scale structure of ~1.6 million DHS and show that the overwhelming majority is populated by well-spaced sites of single transcription factor:DNA interaction. Cell context-dependent cis-regulation is chiefly executed by wholesale actuation of accessibility at regulatory DNA versus by differential transcription factor occupancy within accessible elements. We show further that the well-described enrichment of disease- and phenotypic trait-associated genetic variants in regulatory regions is almost entirely attributable to variants localizing within footprints, and that functional variants impacting transcription factor occupancy are nearly evenly partitioned between loss- and gain-of-function alleles. Unexpectedly, we find that the global density of human genetic variation is markedly increased within transcription factor footprints, revealing an unappreciated driver of cis-regulatory evolution. Our results provide a new framework for both global and nucleotide-precision analyses of gene regulatory mechanisms and functional genetic variation.

Non-invasive epigenome editing is a promising strategy for engineering gene expression programs, yet potency, specificity, and persistence remain challenging. Here we show that effective epigenome editing is gated at single-base precision via 'keyhole' sites in endogenous regulatory DNA. Synthetic repressors targeting promoter keyholes can ablate gene expression in up to 99% of primary cells with single-gene specificity and can seamlessly repress multiple genes in combination. Transient exposure of primary T cells to keyhole repressors confers mitotically heritable silencing that persists to the limit of primary cultures in vitro and for at least 4 weeks in vivo, enabling manufacturing of cell products with enhanced therapeutic efficacy. DNA recognition and effector domains can be encoded as separate proteins that reassemble at keyhole sites and function with the same efficiency as single chain effectors, enabling gated control and rapid screening for novel functional domains that modulate endogenous gene expression patterns. Our results provide a powerful and exponentially flexible system for programming gene expression and therapeutic cell products.

DNase I hypersensitive sites (DHSs) are generic markers of regulatory DNA 1 – 5 and contain genetic variations associated with diseases and phenotypic traits 6 – 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.

DNase I hypersensitive sites (DHSs) are generic markers of regulatory DNA and harbor disease- and phenotypic trait-associated genetic variation. We established high-precision maps of DNase I hypersensitive sites from 733 human biosamples encompassing 439 cell and tissue types and states, and integrated these to precisely delineate and numerically index ~3.6 million DHSs encoded within the human genome, providing a common coordinate system for regulatory DNA. Here we show that the expansive scale of cell and tissue states sampled exposes an unprecedented degree of stereotyped actuation of large sets of elements, signaling the operation of distinct genome-scale regulatory programs. We show further that the complex actuation patterns of individual elements can be captured comprehensively by a simple regulatory vocabulary reflecting their dominant cellular manifestation. This vocabulary, in turn, enables comprehensive and quantitative regulatory annotation of both protein-coding genes and the vast array of well-defined but poorly-characterized non-coding RNA genes. Finally, we show that the combination of high-precision DHSs and regulatory vocabularies markedly concentrate disease- and trait-associated non-coding genetic signals both along the genome and across cellular compartments. Taken together, our results provide a common and extensible coordinate system and vocabulary for human regulatory DNA, and a new global perspective on the architecture of human gene regulation.

Background Linking genetic risk loci identified by genome-wide association studies (GWAS)to their causal genes remains a major challenge. Disease-associated genetic variants are concentrated in regions containing regulatory DNA elements, such as promoters and enhancers. Although researchers have previously published DNA maps of these regulatory regions for kidney tubule cells and glomerular endothelial cells, maps for podocytes and mesangial cells have not been available. Methods We generated regulatory DNA maps (DNase-seq) and paired gene expression profiles (RNA-seq) from primary outgrowth cultures of human glomeruli that were composed mainly of podocytes and mesangial cells. We generated similar datasets from renal cortex cultures, to compare with those of the glomerular cultures. Because regulatory DNA elements can act on target genes across large genomic distances, we also generated a chromatin conformation map from freshly isolated human glomeruli. Results We identified thousands of unique regulatory DNA elements, many located close to transcription factor genes, which the glomerular and cortex samples expressed at different levels. We found that genetic variants associated with kidney diseases (GWAS) and kidney expression quantitative trait loci were enriched in regulatory DNA regions. By combining GWAS, epigenomic, and chromatin conformation data, we functionally annotated 46 kidney disease genes. Conclusions We demonstrate a powerful approach to functionally connect kidney disease-/trait-associated loci to their target genes by leveraging unique regulatory DNA maps and integrated epigenomic and genetic analysis. This process can be applied to other kidney cell types and will enhance our understanding of genome regulation and its effects on gene expression in kidney disease.

BackgroundTranscriptional 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.MethodsWe 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).FindingsThough 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.InterpretationIntegrated 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.

Transcriptional dysregulation drives cancer formation but the underlying mechanisms are still poorly understood. As a model system, we used renal cell carcinoma (RCC), the most common malignant kidney tumor which canonically activates the hypoxia-inducible transcription factor (HIF) pathway. We performed genome-wide chromatin accessibility and transcriptome profiling on paired tumor/normal samples and found that numerous transcription factors with a RCC-selective expression pattern also demonstrated evidence of HIF binding in the vicinity of their gene body. Some of these transcription factors influenced the tumor’s regulatory landscape, notably the stem cell transcription factor POU5F1 ( OCT4 ). Unexpectedly, we discovered a HIF-pathway-responsive cryptic promoter embedded within a human-specific retroviral repeat element that drives POU5F1 expression in RCC via a novel transcript. Elevat POU5F1 expression levels were correlated with advanced tumor stage and poorer overall survival in RCC patients. Thus, integrated transcriptomic and epigenomic analysis of even a small number of primary patient samples revealed remarkably convergent shared regulatory landscapes and a novel mechanism for dysregulated expression of POU5F1 in RCC.

Context: Articular cartilage possesses poor natural healing mechanisms, and a variety of non-cell-based and cell-based treatments aim to promote regeneration of hyaline cartilage. Data Sources: A review of the literature to December 2013 using PubMed with search criteria including the keywords stem cell, cell therapy, cell transplantation, cartilage, chondral, and chondrogenic. Study Selection: Forty-five articles were identified that employed local mesenchymal stem cell (MSC) therapy for joint disorders in humans. Nine comparative studies were identified, consisting of 3 randomized trials, 5 cohort studies, and 1 case-control study. Study Type: Clinical review. Level of Evidence: Level 4. Data Extraction: Studies were assessed for stem cell source, method of implantation, comparison groups, and concurrent surgical techniques. Results: Two studies comparing MSC treatment to autologous chondrocyte implantation found similar efficacy. Three studies reported clinical benefits with intra-articular MSC injection over non-MSC controls for cases undergoing debridement with or without marrow stimulation, although a randomized study found no significant clinical difference at 2-year follow-up but reported better 18-month magnetic resonance imaging and histologic scores in the MSC group. No human studies have compared intra-articular MSC therapy to non-MSC techniques for osteoarthritis in the absence of surgery. Conclusion: Mesenchymal stem cell–based therapies appear safe and effective for joint disorders in large animal preclinical models. Evidence for use in humans, particularly, comparison with more established treatments such as autologous chondrocyte implantation and microfracture, is limited.

The laboratory mouse shares the majority of its protein-coding genes with humans, making it the premier model organism in biomedical research, yet the two mammals differ in significant ways. To gain greater insights into both shared and species-specific transcriptional and cellular regulatory programs in the mouse, the Mouse ENCODE Consortium has mapped transcription, DNase I hypersensitivity, transcription factor binding, chromatin modifications and replication domains throughout the mouse genome in diverse cell and tissue types. By comparing with the human genome, we not only confirm substantial conservation in the newly annotated potential functional sequences, but also find a large degree of divergence of sequences involved in transcriptional regulation, chromatin state and higher order chromatin organization. Our results illuminate the wide range of evolutionary forces acting on genes and their regulatory regions, and provide a general resource for research into mammalian biology and mechanisms of human diseases.

To study the evolutionary dynamics of regulatory DNA, we mapped >1.3 million deoxyribonuclease I–hypersensitive sites (DHSs) in 45 mouse cell and tissue types, and systematically compared these with human DHS maps from orthologous compartments. We found that the mouse and human genomes have undergone extensive cis-regulatory rewiring that combines branch-specific evolutionary innovation and loss with widespread repurposing of conserved DHSs to alternative cell fates, and that this process is mediated by turnover of transcription factor (TF) recognition elements. Despite pervasive evolutionary remodeling of the location and content of individual cis-regulatory regions, within orthologous mouse and human cell types the global fraction of regulatory DNA bases encoding recognition sites for each TF has been strictly conserved. Our findings provide new insights into the evolutionary forces shaping mammalian regulatory DNA landscapes.

The basic body plan and major physiological axes have been highly conserved during mammalian evolution, yet only a small fraction of the human genome sequence appears to be subject to evolutionary constraint. To quantify cis- versus trans-acting contributions to mammalian regulatory evolution, we performed genomic DNase I footprinting of the mouse genome across 25 cell and tissue types, collectively defining ∼8.6 million transcription factor (TF) occupancy sites at nucleotide resolution. Here we show that mouse TF footprints conjointly encode a regulatory lexicon that is ∼95% similar with that derived from human TF footprints. However, only ∼20% of mouse TF footprints have human orthologues. Despite substantial turnover of the cis-regulatory landscape, nearly half of all pairwise regulatory interactions connecting mouse TF genes have been maintained in orthologous human cell types through evolutionary innovation of TF recognition sequences. Furthermore, the higher-level organization of mouse TF-to-TF connections into cellular network architectures is nearly identical with human. Our results indicate that evolutionary selection on mammalian gene regulation is targeted chiefly at the level of trans-regulatory circuitry, enabling and potentiating cis-regulatory plasticity.

DNase I hypersensitive sites (DHSs) are markers of regulatory DNA and have underpinned the discovery of all classes of cis-regulatory elements including enhancers, promoters, insulators, silencers and locus control regions. Here we present the first extensive map of human DHSs identified through genome-wide profiling in 125 diverse cell and tissue types. We identify ∼2.9 million DHSs that encompass virtually all known experimentally validated cis-regulatory sequences and expose a vast trove of novel elements, most with highly cell-selective regulation. Annotating these elements using ENCODE data reveals novel relationships between chromatin accessibility, transcription, DNA methylation and regulatory factor occupancy patterns. We connect ∼580,000 distal DHSs with their target promoters, revealing systematic pairing of different classes of distal DHSs and specific promoter types. Patterning of chromatin accessibility at many regulatory regions is organized with dozens to hundreds of co-activated elements, and the transcellular DNase I sensitivity pattern at a given region can predict cell-type-specific functional behaviours. The DHS landscape shows signatures of recent functional evolutionary constraint. However, the DHS compartment in pluripotent and immortalized cells exhibits higher mutation rates than that in highly differentiated cells, exposing an unexpected link between chromatin accessibility, proliferative potential and patterns of human variation. An extensive map of human DNase I hypersensitive sites, markers of regulatory DNA, in 125 diverse cell and tissue types is described; integration of this information with other ENCODE-generated data sets identifies new relationships between chromatin accessibility, transcription, DNA methylation and regulatory factor occupancy patterns. This paper describes the first extensive map of human DNaseI hypersensitive sites — markers of regulatory DNA — in 125 diverse cell and tissue types. Integration of this information with other data sets generated by ENCODE (Encyclopedia of DNA Elements) identified new relationships between chromatin accessibility, transcription, DNA methylation and regulatory-factor occupancy patterns. Evolutionary-conservation analysis revealed signatures of recent functional constraint within DNaseI hypersensitive sites.