The BET Protein BRD2 Cooperates with CTCF to Enforce Transcriptional and Architectural Boundaries.

Bromodomain and extraterminal motif (BET) proteins are pharmacologic targets for the treatment of diverse diseases, yet the roles of individual BET family members remain unclear. We find that BRD2, but not BRD4, co-localizes with the architectural/insulator protein CCCTC-binding factor (CTCF) genome-wide. CTCF recruits BRD2 to co-bound sites whereas BRD2 is dispensable for CTCF occupancy. Disruption of a CTCF/BRD2-occupied element positioned between two unrelated genes enables regulatory influence to spread from one gene to another, suggesting that CTCF and BRD2 form a transcriptional boundary. Accordingly, single-molecule mRNA fluorescence in situ hybridization (FISH) reveals that, upon site-specific CTCF disruption or BRD2 depletion, expression of the two genes becomes increasingly correlated. HiC shows that BRD2 depletion weakens boundaries co-occupied by CTCF and BRD2, but not those that lack BRD2. These findings indicate that BRD2 supports boundary activity, and they raise the possibility that pharmacologic BET inhibitors can influence gene expression in part by perturbing domain boundary function.

A hyperactive transcriptional state marks genome reactivation at the mitosis-G1 transition.

During mitosis, RNA polymerase II (Pol II) and many transcription factors dissociate from chromatin, and transcription ceases globally. Transcription is known to restart in bulk by telophase, but whether de novo transcription at the mitosis-G1 transition is in any way distinct from later in interphase remains unknown. We tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. Unexpectedly, during the earliest rounds of transcription at the mitosis-G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. Enhancer-promoter chromatin contacts are depleted during mitosis and restored rapidly upon G1 entry but do not spike. Of the chromatin-associated features examined, histone H3 Lys27 acetylation levels at individual loci in mitosis best predict the mitosis-G1 transcriptional spike. Single-molecule RNA imaging supports that the mitosis-G1 transcriptional spike can constitute the maximum transcriptional activity per DNA copy throughout the cell division cycle. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression. Our results raise the possibility that passage through the mitosis-G1 transition might predispose cells to diverge in gene expression states.

Genetic variant pathogenicity prediction trained using disease-specific clinical sequencing data sets.

Recent advances in DNA sequencing have expanded our understanding of the molecular basis of genetic disorders and increased the utilization of clinical genomic tests. Given the paucity of evidence to accurately classify each variant and the difficulty of experimentally evaluating its clinical significance, a large number of variants generated by clinical tests are reported as variants of unknown clinical significance. Population-scale variant databases can improve clinical interpretation. Specifically, pathogenicity prediction for novel missense variants can use features describing regional variant constraint. Constrained genomic regions are those that have an unusually low variant count in the general population. Computational methods have been introduced to capture these regions and incorporate them into pathogenicity classifiers, but these methods have yet to be compared on an independent clinical variant data set. Here, we introduce one variant data set derived from clinical sequencing panels and use it to compare the ability of different genomic constraint metrics to determine missense variant pathogenicity. This data set is compiled from 17,071 patients surveyed with clinical genomic sequencing for cardiomyopathy, epilepsy, or RASopathies. We further use this data set to demonstrate the necessity of disease-specific classifiers and to train PathoPredictor, a disease-specific ensemble classifier of pathogenicity based on regional constraint and variant-level features. PathoPredictor achieves an average precision >90% for variants from all 99 tested disease genes while approaching 100% accuracy for some genes. The accumulation of larger clinical variant training data sets can significantly enhance their performance in a disease- and gene-specific manner.

RH genotype matching for transfusion support in sickle cell disease.

Rh alloimmunization remains a challenge for patients with sickle cell disease (SCD) despite transfusion of serologic Rh C, E, and K antigen-matched red cells. Inheritance of altered RH alleles contributes to the prevalence of Rh antibodies after blood transfusion in patients with SCD and explains approximately one-third of cases. The remainder seem to be stimulated by altered Rh proteins on African American donor red cells. Matching patients with donors on the basis of RH genotype may mitigate Rh alloimmunization, but the feasibility and resources required are not known. We compared RH allele frequencies between patients with SCD (n = 857) and African American donors (n = 587) and showed that RH allele frequencies are similar. Overall, 29% of RHD and 53% of RHCE alleles are altered in patients and African American donors. We modeled RH genotype matching compared with serologic Rh D, C, and E, along with K antigen matching, and found that approximately twice the number of African American donors would be required for RH genotype vs Rh serologic matching at our institution. We demonstrated that African American donor recruitment is necessary to maintain an adequate supply of C-, E-, and K-negative donor units to avoid depleting the Rh-negative (RhD-) blood supply. Our results suggest that prophylactic RH genetic matching for patients with SCD is feasible with a donor pool comprised primarily of African-Americans and would optimize the use of our existing minority donor inventory. The current cost of RH genotyping all minority donors and management of the data remain limiting factors.

Functions of BET proteins in erythroid gene expression.

Inhibitors of bromodomain and extraterminal motif proteins (BETs) are being evaluated for the treatment of cancer and other diseases, yet much remains to be learned about how BET proteins function during normal physiology. We used genomic and genetic approaches to examine BET function in a hematopoietic maturation system driven by GATA1, an acetylated transcription factor previously shown to interact with BETs. We found that BRD2, BRD3, and BRD4 were variably recruited to GATA1-regulated genes, with BRD3 binding the greatest number of GATA1-occupied sites. Pharmacologic BET inhibition impaired GATA1-mediated transcriptional activation, but not repression, genome-wide. Mechanistically, BETs promoted chromatin occupancy of GATA1 and subsequently supported transcriptional activation. Using a combination of CRISPR-Cas9-mediated genomic engineering and shRNA approaches, we observed that depletion of either BRD2 or BRD4 alone blunted erythroid gene activation. Surprisingly, depletion of BRD3 only affected erythroid transcription in the context of BRD2 deficiency. Consistent with functional overlap among BET proteins, forced BRD3 expression substantially rescued defects caused by BRD2 deficiency. These results suggest that pharmacologic BET inhibition should be interpreted in the context of distinct steps in transcriptional activation and overlapping functions among BET family members.

Red blood cell alloimmunization in transfused patients with bone marrow failure syndromes.

BACKGROUND: Red blood cell (RBC) alloimmunization is a concern for patients who receive multiple or chronic transfusions. Alloimmunization prevalence in transfused patients with bone marrow failure syndrome (BMFS) is unknown. This study aimed to determine physician practice for RBC antigen matching, immunization rates, and antibody specificities in patients with BMFS. STUDY DESIGN AND METHODS: The clinical records of all patients with BMFS seen at the Children's Hospital of Philadelphia between 2001 and 2015 were reviewed. Immunization rate was determined per 100 units transfused. RESULTS: ABO/D, C, E, and K (CEK) RBC matching was requested for 21.8% of patients. A total of 3782 RBC units were transfused to 87 patients, of which 2551 (67.5%) were CEK matched and 1231 (32.5%) were ABO/D only matched. The majority of units transfused to patients on a chronic transfusion regimen were CEK matched (89.6% of 2728 units). No anti-C, -E, or -K antibodies formed in any patient during the 14-year study period. Two alloantibodies and two autoantibodies formed, resulting in a rate of 0.05 alloantibodies and 0.05 autoantibodies per 100 units transfused. The prevalence of alloimmunization was 2.3%. CONCLUSION: The rate and prevalence of RBC alloimmunization were low in patients with BMFS. CEK matching avoided alloimmunization to these antigens in chronically transfused patients.

Application of simultaneous selective pressures slows adaptation.

Beneficial mutations that arise in an evolving asexual population may compete or interact in ways that alter the overall rate of adaptation through mechanisms such as clonal or functional interference. The application of multiple selective pressures simultaneously may allow for a greater number of adaptive mutations, increasing the opportunities for competition between selectively advantageous alterations, and thereby reducing the rate of adaptation. We evolved a strain of Saccharomyces cerevisiae that could not produce its own histidine or uracil for ~500 generations under one or three selective pressures: limitation of the concentration of glucose, histidine, and/or uracil in the media. The rate of adaptation was obtained by measuring evolved relative fitness using competition assays. Populations evolved under a single selective pressure showed a statistically significant increase in fitness on those pressures relative to the ancestral strain, but the populations evolved on all three pressures did not show a statistically significant increase in fitness over the ancestral strain on any single pressure. Simultaneously limiting three essential nutrients for a population of S. cerevisiae effectively slows the rate of evolution on any one of the three selective pressures applied, relative to the single selective pressure cases. We identify possible mechanisms for fitness changes seen between populations evolved on one or three limiting nutrient pressures by high-throughput sequencing. Adding multiple selective pressures to evolving disease like cancer and infectious diseases could reduce the rate of adaptation and thereby may slow disease progression, prolong drug efficacy and prevent deaths.

Rapid and accurate interpretation of clinical exomes using Phenoxome: a computational phenotype-driven approach.

Clinical exome sequencing (CES) has become the preferred diagnostic platform for complex pediatric disorders with suspected monogenic etiologies. Despite rapid advancements, the major challenge still resides in identifying the casual variants among the thousands of variants detected during CES testing, and thus establishing a molecular diagnosis. To improve the clinical exome diagnostic efficiency, we developed Phenoxome, a robust phenotype-driven model that adopts a network-based approach to facilitate automated variant prioritization. Phenoxome dissects the phenotypic manifestation of a patient in concert with their genomic profile to filter and then prioritize variants that are likely to affect the function of the gene (potentially pathogenic variants). To validate our method, we have compiled a clinical cohort of 105 positive patient samples that represent a wide range of genetic heterogeneity. Phenoxome identifies the causative variants within the top 5, 10, or 25 candidates in more than 50%, 71%, or 88% of these exomes, respectively. Furthermore, we show that our method is optimized for clinical testing by outperforming the current state-of-art method. We have demonstrated the performance of Phenoxome using a clinical cohort and showed that it enables rapid and accurate interpretation of clinical exomes. Phenoxome is available at https://phenoxome.chop.edu/ .

The Development and Validation of Clinical Exome-Based Panels Using ExomeSlicer: Considerations and Proof of Concept Using an Epilepsy Panel.

Exome-based panels are becoming the preferred diagnostic strategy in clinical laboratories. This approach enables dynamic gene content update and, if needed, cost-effective reflex to whole-exome sequencing. Currently, no guidelines or appropriate resources are available to support the clinical implementation of exome-based panels. Here, we highlight principles and important considerations for the clinical development and validation of exome-based panels. In addition, we developed ExomeSlicer, a novel, web-based resource, which uses empirical exon-level next-generation sequencing quality metrics to predict and visualize technically challenging exome-wide regions in any gene or genes of interest. Exome sequencing data from 100 clinical epilepsy cases were used to illustrate the clinical utility of ExomeSlicer in predicting poor-quality regions and its impact on streamlining the ad hoc Sanger sequencing fill in burden. With the use of ExomeSlicer, >2100 low complexity and/or high-homology regions affecting >1615 genes across the exome were also characterized. These regions can be a source of false-positive or false-negative variant calls, which can lead to misdiagnoses in tested patients and/or inaccurate functional annotations. We provide important considerations and a novel resource for the clinical development of exome-based panels.

Exploiting genetic variation to uncover rules of transcription factor binding and chromatin accessibility.

Single-nucleotide variants that underlie phenotypic variation can affect chromatin occupancy of transcription factors (TFs). To delineate determinants of in vivo TF binding and chromatin accessibility, we introduce an approach that compares ChIP-seq and DNase-seq data sets from genetically divergent murine erythroid cell lines. The impact of discriminatory single-nucleotide variants on TF ChIP signal enables definition at single base resolution of in vivo binding characteristics of nuclear factors GATA1, TAL1, and CTCF. We further develop a facile complementary approach to more deeply test the requirements of critical nucleotide positions for TF binding by combining CRISPR-Cas9-mediated mutagenesis with ChIP and targeted deep sequencing. Finally, we extend our analytical pipeline to identify nearby contextual DNA elements that modulate chromatin binding by these three TFs, and to define sequences that impact kb-scale chromatin accessibility. Combined, our approaches reveal insights into the genetic basis of TF occupancy and their interplay with chromatin features.

Induction of the Immunoproteasome Subunit Lmp7 Links Proteostasis and Immunity in α-Synuclein Aggregation Disorders.

Accumulation of aggregated α-synuclein into Lewy bodies is thought to contribute to the onset and progression of dopaminergic neuron degeneration in Parkinson's disease (PD) and related disorders. Although protein aggregation is associated with perturbation of proteostasis, how α-synuclein aggregation affects the brain proteome and signaling remains uncertain. In a mouse model of α-synuclein aggregation, 6% of 6215 proteins and 1.6% of 8183 phosphopeptides changed in abundance, indicating conservation of proteostasis and phosphorylation signaling. The proteomic analysis confirmed changes in abundance of proteins that regulate dopamine synthesis and transport, synaptic activity and integrity, and unearthed changes in mRNA binding, processing and protein translation. Phosphorylation signaling changes centered on axonal and synaptic cytoskeletal organization and structural integrity. Proteostatic responses included a significant increase in the levels of Lmp7, a component of the immunoproteasome. Increased Lmp7 levels and activity were also quantified in postmortem human brains with PD and dementia with Lewy bodies. Functionally, the immunoproteasome degrades α-synuclein aggregates and generates potentially antigenic peptides. Expression and activity of the immunoproteasome may represent testable targets to induce adaptive responses that maintain proteome integrity and modulate immune responses in protein aggregation disorders.

Oral nitrite restores age-dependent phenotypes in eNOS-null mice.

Alterations in the synthesis and bioavailability of NO are central to the pathogenesis of cardiovascular and metabolic disorders. Although endothelial NO synthase-derived (eNOS-derived) NO affects mitochondrial long-chain fatty acid ß-oxidation, the pathophysiological significance of this regulation remains unclear. Accordingly, we determined the contributions of eNOS/NO signaling in the adaptive metabolic responses to fasting and in age-induced metabolic dysfunction. Four-month-old eNOS-/- mice are glucose intolerant and exhibit serum dyslipidemia and decreased capacity to oxidize fatty acids. However, during fasting, eNOS-/- mice redirect acetyl-CoA to ketogenesis to elevate circulating levels of ß-hydroxybutyrate similar to wild-type mice. Treatment of 4-month-old eNOS-/- mice with nitrite for 10 days corrected the hypertension and serum hyperlipidemia and normalized the rate of fatty acid oxidation. Fourteen-month-old eNOS-/- mice exhibited metabolic derangements, resulting in reduced utilization of fat to generate energy, lower resting metabolic activity, and diminished physical activity. Seven-month administration of nitrite to eNOS-/- mice reversed the age-dependent metabolic derangements and restored physical activity. While the eNOS/NO signaling is not essential for the metabolic adaptation to fasting, it is critical for regulating systemic metabolic homeostasis in aging. The development of age-dependent metabolic disorder is prevented by low-dose replenishment of bioactive NO.

Common variants in MMP20 at 11q22.2 predispose to 11q deletion and neuroblastoma risk.

MYCN amplification and 11q deletion are two inversely correlated prognostic factors of poor outcome in neuroblastoma. Here we identify common variants at 11q22.2 within MMP20 that associate with neuroblastoma cases harboring 11q deletion (rs10895322), using GWAS in 113 European-American cases and 5109 ancestry-matched controls. The association is replicated in 44 independent cases and 1902 controls. Our study yields novel insights into the genetic underpinnings of neuroblastoma, demonstrating that the inherited common variants reported contribute to the origin of intra-tumor genetic heterogeneity in neuroblastoma.Chromosomal abnormalities such as 11q deletion are associated with poor prognosis in neuroblastoma. Here, the authors perform a genome-wide association study and identify an association between a variant within a Matrix metalloproteinase (MMP) gene member, MMP20, and 11q-deletion subtype neuroblastoma.

Site-Specific Proteomic Mapping Identifies Selectively Modified Regulatory Cysteine Residues in Functionally Distinct Protein Networks.

S-Acylation, S-glutathionylation, S-nitrosylation, and S-sulfenylation are prominent, chemically distinct modifications that regulate protein function, redox sensing, and trafficking. Although the biological significance of these modifications is increasingly appreciated, their integration in the proteome remains unknown. Novel mass spectrometry-based technologies identified 2,596 predominately unique sites in 1,319 mouse liver proteins under physiological conditions. Structural analysis localized the modifications in unique, evolutionary conserved protein segments, outside commonly annotated functional regions. Contrary to expectations, propensity for modification did not correlate with biophysical properties that regulate cysteine reactivity. However, the in vivo chemical reactivity is fine-tuned for specificity, demonstrated by the nominal complementation between the four modifications and quantitative proteomics which showed that a reduction in S-nitrosylation is not correlated with increased S-glutathionylation. A comprehensive survey uncovered clustering of modifications within biologically related protein networks. The data provide the first evidence for the occurrence of distinct, endogenous protein networks that undergo redox signaling through specific cysteine modifications.

Regulation of brain glutamate metabolism by nitric oxide and S-nitrosylation.

Nitric oxide (NO) is a signaling intermediate during glutamatergic neurotransmission in the central nervous system (CNS). NO signaling is in part accomplished through cysteine S-nitrosylation, a posttranslational modification by which NO regulates protein function and signaling. In our investigation of the protein targets and functional impact of S-nitrosylation in the CNS under physiological conditions, we identified 269 S-nitrosocysteine residues in 136 proteins in the wild-type mouse brain. The number of sites was significantly reduced in the brains of mice lacking endothelial nitric oxide synthase (eNOS(-/-)) or neuronal nitric oxide synthase (nNOS(-/-)). In particular, nNOS(-/-) animals showed decreased S-nitrosylation of proteins that participate in the glutamate/glutamine cycle, a metabolic process by which synaptic glutamate is recycled or oxidized to provide energy. (15)N-glutamine-based metabolomic profiling and enzymatic activity assays indicated that brain extracts from nNOS(-/-) mice converted less glutamate to glutamine and oxidized more glutamate than those from mice of the other genotypes. GLT1 [also known as EAAT2 (excitatory amino acid transporter 2)], a glutamate transporter in astrocytes, was S-nitrosylated at Cys(373) and Cys(561) in wild-type and eNOS(-/-) mice, but not in nNOS(-/-) mice. A form of rat GLT1 that could not be S-nitrosylated at the equivalent sites had increased glutamate uptake compared to wild-type GLT1 in cells exposed to an S-nitrosylating agent. Thus, NO modulates glutamatergic neurotransmission through the selective, nNOS-dependent S-nitrosylation of proteins that govern glutamate transport and metabolism.

Exome sequencing identifies recurrent mutations in NF1 and RASopathy genes in sun-exposed melanomas.

We report on whole-exome sequencing (WES) of 213 melanomas. Our analysis established NF1, encoding a negative regulator of RAS, as the third most frequently mutated gene in melanoma, after BRAF and NRAS. Inactivating NF1 mutations were present in 46% of melanomas expressing wild-type BRAF and RAS, occurred in older patients and showed a distinct pattern of co-mutation with other RASopathy genes, particularly RASA2. Functional studies showed that NF1 suppression led to increased RAS activation in most, but not all, melanoma cases. In addition, loss of NF1 did not predict sensitivity to MEK or ERK inhibitors. The rebound pathway, as seen by the induction of phosphorylated MEK, occurred in cells both sensitive and resistant to the studied drugs. We conclude that NF1 is a key tumor suppressor lost in melanomas, and that concurrent RASopathy gene mutations may enhance its role in melanomagenesis.

Computational analysis in cancer exome sequencing.

Exome sequencing in cancer is a powerful tool for identifying mutational events across the coding region of human genes. Here, we describe computational methods that use exome sequencing reads from cancer samples to identify somatic single nucleotide variants (SNVs), copy number alterations, and short insertions and deletions (InDels). We further describe analytical methods to generate lists of driver genes with more mutational events than expected by chance.

Identification of PLX4032-resistance mechanisms and implications for novel RAF inhibitors.

BRAF inhibitors improve melanoma patient survival, but resistance invariably develops. Here we report the discovery of a novel BRAF mutation that confers resistance to PLX4032 employing whole-exome sequencing of drug-resistant BRAF(V600K) melanoma cells. We further describe a new screening approach, a genome-wide piggyBac mutagenesis screen that revealed clinically relevant aberrations (N-terminal BRAF truncations and CRAF overexpression). The novel BRAF mutation, a Leu505 to His substitution (BRAF(L505H) ), is the first resistance-conferring second-site mutation identified in BRAF mutant cells. The mutation replaces a small nonpolar amino acid at the BRAF-PLX4032 interface with a larger polar residue. Moreover, we show that BRAF(L505H) , found in human prostate cancer, is itself a MAPK-activating, PLX4032-resistant oncogenic mutation. Lastly, we demonstrate that the PLX4032-resistant melanoma cells are sensitive to novel, next-generation BRAF inhibitors, especially the 'paradox-blocker' PLX8394, supporting its use in clinical trials for treatment of melanoma patients with BRAF-mutations.

Correlated evolution of positions within mammalian cis elements.

Transcriptional regulation critically depends on proper interactions between transcription factors (TF) and their cognate DNA binding sites. The widely used model of TF-DNA binding--the positional weight matrix (PWM)--presumes independence between positions within the binding site. However, there is evidence to show that the independence assumption may not always hold, and the extent of interposition dependence is not completely known. We hypothesize that the interposition dependence should partly be manifested as correlated evolution at the positions. We report a maximum-likelihood (ML) approach to infer correlated evolution at any two positions within a PWM, based on a multiple alignment of 5 mammalian genomes. Application to a genome-wide set of putative cis elements in human promoters reveals a prevalence of correlated evolution within cis elements. We found that the interdependence between two positions decreases with increasing distance between the positions. The interdependent positions tend to be evolutionarily more constrained and moreover, the dependence patterns are relatively similar across structurally related transcription factors. Although some of the detected mutational dependencies may be due to context-dependent genomic hyper-mutation, notably CG to TG, the majority is likely due to context-dependent preferences for specific nucleotide combinations within the cis elements. Patterns of evolution at individual nucleotide positions within mammalian TF binding sites are often significantly correlated, suggesting interposition dependence. The proposed methodology is also applicable to other classes of non-coding functional elements. A detailed investigation of mutational dependencies within specific motifs could reveal preferred nucleotide combinations that may help refine the DNA binding models.

Adjusting for background mutation frequency biases improves the identification of cancer driver genes.

A common goal of tumor sequencing projects is finding genes whose mutations are selected for during tumor development. This is accomplished by choosing genes that have more non-synonymous mutations than expected from an estimated background mutation frequency. While this background frequency is unknown, it can be estimated using both the observed synonymous mutation frequency and the non-synonymous to synonymous mutation ratio. The synonymous mutation frequency can be determined across all genes or in a gene-specific manner. This choice introduces an interesting trade-off. A gene-specific frequency adjusts for an underlying mutation bias, but is difficult to estimate given missing synonymous mutation counts. Using a genome-wide synonymous frequency is more robust, but is less suited for adjusting biases. Studying four evaluation criteria for identifying genes with high non-synonymous mutation burden (reflecting preferential selection of expressed genes, genes with mutations in conserved bases, genes with many protein interactions, and genes that show loss of heterozygosity), we find that the gene-specific synonymous frequency is superior in the gene expression and protein interaction tests. In conclusion, the use of the gene-specific synonymous mutation frequency is well suited for assessing a gene's non-synonymous mutation burden.