Systematic bromodomain protein screens identify homologous recombination and R-loop suppression pathways involved in genome integrity.

Bromodomain proteins (BRD) are key chromatin regulators of genome function and stability as well as therapeutic targets in cancer. Here, we systematically delineate the contribution of human BRD proteins for genome stability and DNA double-strand break (DSB) repair using several cell-based assays and proteomic interaction network analysis. Applying these approaches, we identify 24 of the 42 BRD proteins as promoters of DNA repair and/or genome integrity. We identified a BRD-reader function of PCAF that bound TIP60-mediated histone acetylations at DSBs to recruit a DUB complex to deubiquitylate histone H2BK120, to allowing direct acetylation by PCAF, and repair of DSBs by homologous recombination. We also discovered the bromo-and-extra-terminal (BET) BRD proteins, BRD2 and BRD4, as negative regulators of transcription-associated RNA-DNA hybrids (R-loops) as inhibition of BRD2 or BRD4 increased R-loop formation, which generated DSBs. These breaks were reliant on topoisomerase II, and BRD2 directly bound and activated topoisomerase I, a known restrainer of R-loops. Thus, comprehensive interactome and functional profiling of BRD proteins revealed new homologous recombination and genome stability pathways, providing a framework to understand genome maintenance by BRD proteins and the effects of their pharmacological inhibition.

Spatiotemporal transcriptional dynamics of the cycling mouse oviduct.

Female fertility in mammals requires iterative remodeling of the entire adult female reproductive tract across the menstrual/estrous cycle. However, while transcriptome dynamics across the estrous cycle have been reported in human and bovine models, no global analysis of gene expression across the estrous cycle has yet been reported for the mouse. Here, we examined the cellular composition and global transcriptional dynamics of the mouse oviduct along the anteroposterior axis and across the estrous cycle. We observed robust patterns of differential gene expression along the anteroposterior axis, but we found surprisingly few changes in gene expression across the estrous cycle. Notable gene expression differences along the anteroposterior axis included a surprising enrichment for genes related to embryonic development, such as Hox and Wnt genes. The relatively stable transcriptional dynamics across the estrous cycle differ markedly from other mammals, leading us to speculate that this is an evolutionarily derived state that may reflect the extremely rapid five-day mouse estrous cycle. This dataset fills a critical gap by providing an important genomic resource for a highly tractable genetic model of mammalian female reproduction.

The histone variant H2A.Z in yeast is almost exclusively incorporated into the +1 nucleosome in the direction of transcription.

Transcription start sites (TSS) in eukaryotes are characterized by a nucleosome-depleted region (NDR), which appears to be flanked upstream and downstream by strongly positioned nucleosomes incorporating the histone variant H2A.Z. H2A.Z associates with both active and repressed TSS and is important for priming genes for rapid transcriptional activation. However, the determinants of H2A.Z occupancy at specific nucleosomes and its relationship to transcription initiation remain unclear. To further elucidate the specificity of H2A.Z, we determined its genomic localization at single nucleosome resolution, as well as the localization of its chromatin remodelers Swr1 and Ino80. By analyzing H2A.Z occupancy in conjunction with RNA expression data that captures promoter-derived antisense initiation, we find that H2A.Z's bimodal incorporation on either side of the NDR is not a general feature of TSS, but is specifically a marker for bidirectional transcription, such that the upstream flanking -1 H2A.Z-containing nucleosome is more appropriately considered as a +1 H2A.Z nucleosome for antisense transcription. The localization of H2A.Z almost exclusively at the +1 nucleosome suggests that a transcription-initiation dependent process could contribute to its specific incorporation.

Ancient eukaryotic protein interactions illuminate modern genetic traits and disorders.

All eukaryotes share a common ancestor from roughly 1.5 - 1.8 billion years ago, a single-celled, swimming microbe known as LECA, the Last Eukaryotic Common Ancestor. Nearly half of the genes in modern eukaryotes were present in LECA, and many current genetic diseases and traits stem from these ancient molecular systems. To better understand these systems, we compared genes across modern organisms and identified a core set of 10,092 shared protein-coding gene families likely present in LECA, a quarter of which are uncharacterized. We then integrated >26,000 mass spectrometry proteomics analyses from 31 species to infer how these proteins interact in higher-order complexes. The resulting interactome describes the biochemical organization of LECA, revealing both known and new assemblies. We analyzed these ancient protein interactions to find new human gene-disease relationships for bone density and congenital birth defects, demonstrating the value of ancestral protein interactions for guiding functional genetics today.

Synthetic repertoires derived from convalescent COVID-19 patients enable discovery of SARS-CoV-2 neutralizing antibodies and a novel quaternary binding modality.

The ongoing evolution of SARS-CoV-2 into more easily transmissible and infectious variants has sparked concern over the continued effectiveness of existing therapeutic antibodies and vaccines. Hence, together with increased genomic surveillance, methods to rapidly develop and assess effective interventions are critically needed. Here we report the discovery of SARS-CoV-2 neutralizing antibodies isolated from COVID-19 patients using a high-throughput platform. Antibodies were identified from unpaired donor B-cell and serum repertoires using yeast surface display, proteomics, and public light chain screening. Cryo-EM and functional characterization of the antibodies identified N3-1, an antibody that binds avidly (Kd,app = 68 pM) to the receptor binding domain (RBD) of the spike protein and robustly neutralizes the virus in vitro. This antibody likely binds all three RBDs of the trimeric spike protein with a single IgG. Importantly, N3-1 equivalently binds spike proteins from emerging SARS-CoV-2 variants of concern, neutralizes UK variant B.1.1.7, and binds SARS-CoV spike with nanomolar affinity. Taken together, the strategies described herein will prove broadly applicable in interrogating adaptive immunity and developing rapid response biological countermeasures to emerging pathogens.

Detection and benchmarking of somatic mutations in cancer genomes using RNA-seq data.

To detect functional somatic mutations in tumor samples, whole-exome sequencing (WES) is often used for its reliability and relative low cost. RNA-seq, while generally used to measure gene expression, can potentially also be used for identification of somatic mutations. However there has been little systematic evaluation of the utility of RNA-seq for identifying somatic mutations. Here, we develop and evaluate a pipeline for processing RNA-seq data from glioblastoma multiforme (GBM) tumors in order to identify somatic mutations. The pipeline entails the use of the STAR aligner 2-pass procedure jointly with MuTect2 from genome analysis toolkit (GATK) to detect somatic variants. Variants identified from RNA-seq data were evaluated by comparison against the COSMIC and dbSNP databases, and also compared to somatic variants identified by exome sequencing. We also estimated the putative functional impact of coding variants in the most frequently mutated genes in GBM. Interestingly, variants identified by RNA-seq alone showed better representation of GBM-related mutations cataloged by COSMIC. RNA-seq-only data substantially outperformed the ability of WES to reveal potentially new somatic mutations in known GBM-related pathways, and allowed us to build a high-quality set of somatic mutations common to exome and RNA-seq calls. Using RNA-seq data in parallel with WES data to detect somatic mutations in cancer genomes can thus broaden the scope of discoveries and lend additional support to somatic variants identified by exome sequencing alone.

Bivalent Chromatin Domains in Glioblastoma Reveal a Subtype-Specific Signature of Glioma Stem Cells.

Glioblastoma multiforme (GBM) can be clustered by gene expression into four main subtypes associated with prognosis and survival, but enhancers and other gene-regulatory elements have not yet been identified in primary tumors. Here, we profiled six histone modifications and CTCF binding as well as gene expression in primary gliomas and identified chromatin states that define distinct regulatory elements across the tumor genome. Enhancers in mesenchymal and classical tumor subtypes drove gene expression associated with cell migration and invasion, whereas enhancers in proneural tumors controlled genes associated with a less aggressive phenotype in GBM. We identified bivalent domains marked by activating and repressive chromatin modifications. Interestingly, the gene interaction network from common (subtype-independent) bivalent domains was highly enriched for homeobox genes and transcription factors and dominated by SHH and Wnt signaling pathways. This subtype-independent signature of early neural development may be indicative of poised dedifferentiation capacity in glioblastoma and could provide potential targets for therapy.Significance: Enhancers and bivalent domains in glioblastoma are regulated in a subtype-specific manner that resembles gene regulation in glioma stem cells. Cancer Res; 78(10); 2463-74. ©2018 AACR.