DEEPOMICS FFPE, a deep neural network model, identifies DNA sequencing artifacts from formalin fixed paraffin embedded tissue with high accuracy.

Formalin-fixed, paraffin-embedded (FFPE) tissue specimens are routinely used in pathological diagnosis, but their large number of artifactual mutations complicate the evaluation of companion diagnostics and analysis of next-generation sequencing data. Identification of variants with low allele frequencies is challenging because existing FFPE filtering tools label all low-frequency variants as artifacts. To address this problem, we aimed to develop DEEPOMICS FFPE, an AI model that can classify a true variant from an artifact. Paired whole exome sequencing data from fresh frozen and FFPE samples from 24 tumors were obtained from public sources and used as training and validation sets at a ratio of 7:3. A deep neural network model with three hidden layers was trained with input features using outputs of the MuTect2 caller. Contributing features were identified using the SHapley Additive exPlanations algorithm and optimized based on training results. The performance of the final model (DEEPOMICS FFPE) was compared with those of existing models (MuTect filter, FFPolish, and SOBDetector) by using well-defined test datasets. We found 41 discriminating properties for FFPE artifacts. Optimization of property quantification improved the model performance. DEEPOMICS FFPE removed 99.6% of artifacts while maintaining 87.1% of true variants, with an F1-score of 88.3 in the entire dataset not used for training, which is significantly higher than those of existing tools. Its performance was maintained even for low-allele-fraction variants with a specificity of 0.995, suggesting that it can be used to identify subclonal variants. Different from existing methods, DEEPOMICS FFPE identified most of the sequencing artifacts in the FFPE samples while retaining more of true variants, including those of low allele frequencies. The newly developed tool DEEPOMICS FFPE may be useful in designing capture panels for personalized circulating tumor DNA assay and identifying candidate neoepitopes for personalized vaccine design. DEEPOMICS FFPE is freely available on the web ( http://deepomics.co.kr/ffpe ) for research.

Transcription and chromatin-based surveillance mechanism controls suppression of cryptic antisense transcription.

Phosphorylation of the RNA polymerase II C-terminal domain Y1S2P3T4S5P6S7 consensus sequence coordinates key events during transcription, and its deregulation leads to defects in transcription and RNA processing. Here, we report that the histone deacetylase activity of the fission yeast Hos2/Set3 complex plays an important role in suppressing cryptic initiation of antisense transcription when RNA polymerase II phosphorylation is dysregulated due to the loss of Ssu72 phosphatase. Interestingly, although single Hos2 and Set3 mutants have little effect, loss of Hos2 or Set3 combined with ssu72Δ results in a synergistic increase in antisense transcription globally and correlates with elevated sensitivity to genotoxic agents. We demonstrate a key role for the Ssu72/Hos2/Set3 mechanism in the suppression of cryptic antisense transcription at the 3' end of convergent genes that are most susceptible to these defects, ensuring the fidelity of gene expression within dense genomes of simple eukaryotes.

Profiling RNA Polymerase II Phosphorylation Genome-Wide in Fission Yeast.

The RNA polymerase II carboxyl-terminal domain (CTD) consists of tandem repeats of consensus sequence Tyr1-Ser2-Pro3-Thr4-Ser5-Pro6-Ser7. Dynamic posttranslational modifications of the CTD generate a CTD code crucial for the cotranscriptional recruitment of factors that control transcription, chromatin modification, and RNA processing. Analysis of CTD phosphorylation by ChIP (Chromatin ImmunoPrecipitation) coupled with high-throughput DNA sequencing (ChIP-seq) is a powerful tool to investigate the changes in CTD phosphorylation during the transcription cycle. In this chapter, we describe a ChIP-seq protocol to profile the different CTD phospho-marks in fission yeast. Using this protocol, we have found that Tyr1P, Ser2P, and Thr4P signals are highest at gene 3' ends, whereas Ser5P is enriched across the gene bodies.

Elongation/Termination Factor Exchange Mediated by PP1 Phosphatase Orchestrates Transcription Termination.

Termination of RNA polymerase II (Pol II) transcription is a key step that is important for 3' end formation of functional mRNA, mRNA release, and Pol II recycling. Even so, the underlying termination mechanism is not yet understood. Here, we demonstrate that the conserved and essential termination factor Seb1 is found on Pol II near the end of the RNA exit channel and the Rpb4/7 stalk. Furthermore, the Seb1 interaction surface with Pol II largely overlaps with that of the elongation factor Spt5. Notably, Seb1 co-transcriptional recruitment is dependent on Spt5 dephosphorylation by the conserved PP1 phosphatase Dis2, which also dephosphorylates threonine 4 within the Pol II heptad repeated C-terminal domain. We propose that Dis2 orchestrates the transition from elongation to termination phase during the transcription cycle by mediating elongation to termination factor exchange and dephosphorylation of Pol II C-terminal domain.

Histone deacetylation promotes transcriptional silencing at facultative heterochromatin.

It is important to accurately regulate the expression of genes involved in development and environmental response. In the fission yeast Schizosaccharomyces pombe, meiotic genes are tightly repressed during vegetative growth. Despite being embedded in heterochromatin these genes are transcribed and believed to be repressed primarily at the level of RNA. However, the mechanism of facultative heterochromatin formation and the interplay with transcription regulation is not understood. We show genome-wide that HDAC-dependent histone deacetylation is a major determinant in transcriptional silencing of facultative heterochromatin domains. Indeed, mutation of class I/II HDACs leads to increased transcription of meiotic genes and accumulation of their mRNAs. Mechanistic dissection of the pho1 gene where, in response to phosphate, transient facultative heterochromatin is established by overlapping lncRNA transcription shows that the Clr3 HDAC contributes to silencing independently of SHREC, but in an lncRNA-dependent manner. We propose that HDACs promote facultative heterochromatin by establishing alternative transcriptional silencing.

The conserved protein Seb1 drives transcription termination by binding RNA polymerase II and nascent RNA.

Termination of RNA polymerase II (Pol II) transcription is an important step in the transcription cycle, which involves the dislodgement of polymerase from DNA, leading to release of a functional transcript. Recent studies have identified the key players required for this process and showed that a common feature of these proteins is a conserved domain that interacts with the phosphorylated C-terminus of Pol II (CTD-interacting domain, CID). However, the mechanism by which transcription termination is achieved is not understood. Using genome-wide methods, here we show that the fission yeast CID-protein Seb1 is essential for termination of protein-coding and non-coding genes through interaction with S2-phosphorylated Pol II and nascent RNA. Furthermore, we present the crystal structures of the Seb1 CTD- and RNA-binding modules. Unexpectedly, the latter reveals an intertwined two-domain arrangement of a canonical RRM and second domain. These results provide important insights into the mechanism underlying eukaryotic transcription termination.

Exosome Cofactors Connect Transcription Termination to RNA Processing by Guiding Terminated Transcripts to the Appropriate Exonuclease within the Nuclear Exosome.

The yeast Nrd1 interacts with the C-terminal domain (CTD) of RNA polymerase II (RNApII) through its CTD-interacting domain (CID) and also associates with the nuclear exosome, thereby acting as both a transcription termination and RNA processing factor. Previously, we found that the Nrd1 CID is required to recruit the nuclear exosome to the Nrd1 complex, but it was not clear which exosome subunits were contacted. Here, we show that two nuclear exosome cofactors, Mpp6 and Trf4, directly and competitively interact with the Nrd1 CID and differentially regulate the association of Nrd1 with two catalytic subunits of the exosome. Importantly, Mpp6 promotes the processing of Nrd1-terminated transcripts preferentially by Dis3, whereas Trf4 leads to Rrp6-dependent processing. This suggests that Mpp6 and Trf4 may play a role in choosing a particular RNA processing route for Nrd1-terminated transcripts within the exosome by guiding the transcripts to the appropriate exonuclease.

The RNA polymerase II C-terminal domain-interacting domain of yeast Nrd1 contributes to the choice of termination pathway and couples to RNA processing by the nuclear exosome.

The RNA polymerase II (RNApII) C-terminal domain (CTD)-interacting domain (CID) proteins are involved in two distinct RNApII termination pathways and recognize different phosphorylated forms of CTD. To investigate the role of differential CTD-CID interactions in the choice of termination pathway, we altered the CTD-binding specificity of Nrd1 by domain swapping. Nrd1 with the CID from Rtt103 (Nrd1(CID(Rtt103))) causes read-through transcription at many genes, but can also trigger termination where multiple Nrd1/Nab3-binding sites and the Ser(P)-2 CTD co-exist. Therefore, CTD-CID interactions target specific termination complexes to help choose an RNApII termination pathway. Interactions of Nrd1 with both CTD and nascent transcripts contribute to efficient termination by the Nrd1 complex. Surprisingly, replacing the Nrd1 CID with that from Rtt103 reduces binding to Rrp6/Trf4, and RNA transcripts terminated by Nrd1(CID(Rtt103)) are predominantly processed by core exosome. Thus, the Nrd1 CID couples Ser(P)-5 CTD not only to termination, but also to RNA processing by the nuclear exosome.

Cd2+ binds to Atx1 and affects the physical interaction between Atx1 and Ccc2 in Saccharomyces cerevisiae.

The ATX1 deletion strain of Saccharomyces cerevisiae is more resistant to Cd(2+) than the wild-type. To investigate the function of Atx1 in Cd(2+) toxicity, we used a metal-binding assay to study the interaction between Atx1 and Cd(2+) in vitro. Using circular dichroism and two-hybrid analyses, we found that Atx1 can bind Cd(2+) specifically and that Cd(2+) binding to Atx1 affects the physical interaction between Atx1 and Ccc2. These results imply that Atx1 delivers Cd(2+) to Ccc2 and that this delivery is, at least in part, responsible for Cd(2+) toxicity in S. cerevisiae.

Cadmium regulates copper homoeostasis by inhibiting the activity of Mac1, a transcriptional activator of the copper regulon, in Saccharomyces cerevisiae.

Cadmium is a toxic metal and the mechanism of its toxicity has been studied in various model systems from bacteria to mammals. We employed Saccharomyces cerevisiae as a model system to study cadmium toxicity at the molecular level because it has been used to identify the molecular mechanisms of toxicity found in higher organisms. cDNA microarray and Northern blot analyses revealed that cadmium salts inhibited the expression of genes related to copper metabolism. Western blotting, Northern blotting and chromatin immunoprecipitation experiments indicated that CTR1 expression was inhibited at the transcriptional level through direct inhibition of the Mac1 transcriptional activator. The decreased expression of CTR1 results in cellular copper deficiency and inhibition of Fet3 activity, which eventually impairs iron uptake. In this way, cadmium exhibits a negative effect on both iron and copper homoeostasis.

FgEnd1 is a putative component of the endocytic machinery and mediates ferrichrome uptake in F. graminearum.

The function of endocytic pathway in filamentous fungi has remained elusive. Recently, we have identified that FgEnd1, which has a 27% amino acid homology and shares specific EH3 domain with ScEnd3 of Saccharomyces cerevisiae, is a putative member of the endocytic machinery in Fusarium graminearum. The failure of the scend3 mutant to uptake Lucifer yellow (LY) was recovered by introducing FgEnd1 into S. cerevisiae. The deletion of fgend1 in F. graminearum resulted in a 2-fold decrease in the rate of uptake of the endocytic marker FM4-64 when compared to wild-type cells. The rate of uptake was similar to that seen in latrunculin A (Lat-A)-treated cells. Furthermore, fgend1 deletion strain of F. graminearum showed lower ferrichrome (FC) uptake activity than wild-type F. graminearum, and the same rate as LatA-treated cells. Taken together, these results suggest that FgEnd1 is a putative member of the endocytic machinery, although it acts through a different mechanism from ScEnd3 or ScEnd4 of S. cerevisiae.

A novel function of Aft1 in regulating ferrioxamine B uptake: Aft1 modulates Arn3 ubiquitination in Saccharomyces cerevisiae.

Aft1 is a transcriptional activator in Saccharomyces cerevisiae that responds to iron availability and regulates the expression of genes in the iron regulon, such as FET3, FTR1 and the ARN family. Using a two-hybrid screen, we found that Aft1 physically interacts with the FOB (ferrioxamine B) transporter Arn3. This interaction modulates the ability of Arn3 to take up FOB. The interaction between Arn3 and Aft1 was confirmed by beta-galactosidase, co-immunoprecipitation and SPR (surface plasmon resonance) assays. Truncated Aft1 had a stronger interaction with Arn3 and caused a higher FOB-uptake activity than full-length Aft1. Interestingly, only full-length Aft1 induced the correct localization of Arn3 in response to FOB. Furthermore, we found Aft1 affected Arn3 ubiquitination. These results suggest that Aft1 interacts with Arn3 and may regulate the ubiquitination of Arn3 in the cytosolic compartment.

A novel role for thioredoxin reductase in the iron metabolism of S. cerevisiae.

Intracellular levels of iron are tightly regulated. Saccharomyces cerevisiae uses well-defined pathways to extract iron molecules from the environment. Once inside the cell, the iron molecules must be transferred to target sites via an intracellular iron transporter. Although analogous carriers have been described for other metals, such as copper, an iron transporter has yet to be identified. We used two-dimensional gel electrophoresis and mass spectrometry techniques to attempt to identify the iron transporter from cytosolic fraction of S. cerevisiae. In this study, we identified the iron-binding activity of thioredoxin reductase, and our data suggest a potential role for this enzyme in intracellular iron transport.

Functional characterization of Hsp33 protein from Bacillus psychrosaccharolyticus; additional function of HSP33 on resistance to solvent stress.

Psychrophiles have been known as efficient organism to degrade organic solvent. To investigate the mechanism of solvent stress and identify the factors that affect the solvent stress in psychrophiles, we selected Bacillus psychrosaccharolyticus one of the psychrophiles and two-dimensional gel electrophoresis was performed. Among the protein spots analyzed by 2-DE, five spots induced in 3% IPA stress conditions were identified by MS/MS, and one of these spots was identified as a Hsp33 family. The Hsp33 protein sequence of B. psychrosaccharolyticus exhibited a high similarity with the corresponding proteins of other bacteria. The Hsp33 protein of B. psychrosaccharolyticus has a highly conserved zinc-binding domain (CXCX, CXXC) that includes four cysteine residues in the C-terminus. In addition, the transcriptional induction of the HSP33 of B. psychrosaccharolyticus was confirmed by Northern blot analysis, and formation of free thiol linkage was induced under stress conditions such as exposure to solvents, heat-shock, and oxidative stress. Furthermore, over-expressed strains of HSP33 of B. psychrosaccharolyticus in Escherichia coli improved stress tolerance to the organic solvent when compared with the wild-type. These data suggest that the solvent stress condition was similar to heat-shock or oxidative stress, especially through the triggering of induction and activation of a redox-regulatory chaperone, Hsp33, and Hsp33 plays a critical role in the tolerance to stress.