PHF20 is crucial for epigenetic control of starvation-induced autophagy through enhancer activation.

Autophagy is a catabolic pathway that maintains cellular homeostasis under various stress conditions, including conditions of nutrient deprivation. To elevate autophagic flux to a sufficient level under stress conditions, transcriptional activation of autophagy genes occurs to replenish autophagy components. Thus, the transcriptional and epigenetic control of the genes regulating autophagy is essential for cellular homeostasis. Here, we applied integrated transcriptomic and epigenomic profiling to reveal the roles of plant homeodomain finger protein 20 (PHF20), which is an epigenetic reader possessing methyl binding activity, in controlling the expression of autophagy genes. Phf20 deficiency led to impaired autophagic flux and autophagy gene expression under glucose starvation. Interestingly, the genome-wide characterization of chromatin states by Assay for Transposase-Accessible Chromatin (ATAC)-sequencing revealed that the PHF20-dependent chromatin remodelling occurs in enhancers that are co-occupied by dimethylated lysine 36 on histone H3 (H3K36me2). Importantly, the recognition of H3K36me2 by PHF20 was found to be highly correlated with increased levels of H3K4me1/2 at the enhancer regions. Collectively, these results indicate that PHF20 regulates autophagy genes through enhancer activation via H3K36me2 recognition as an epigenetic reader. Our findings emphasize the importance of nuclear events in the regulation of autophagy.

Substrate-binding loop interactions with pseudouridine trigger conformational changes that promote catalytic efficiency of pseudouridine kinase PUKI.

Pseudouridine, one major RNA modification, is catabolized into uracil and ribose-5'-phosphate by two sequential enzymatic reactions. In the first step, pseudouridine kinase (PUKI) phosphorylates pseudouridine to pseudouridine 5'-monophosphate. High-fidelity catalysis of pseudouridine by PUKI prevents possible disturbance of in vivo pyrimidine homeostasis. However, the molecular basis of how PUKI selectively phosphorylates pseudouridine over uridine with >100-fold greater efficiency despite minor differences in their Km values has not been elucidated. To investigate this selectivity, in this study we determined the structures of PUKI from Escherichia coli strain B (EcPUKI) in various ligation states. The structure of EcPUKI was determined to be similar to PUKI from Arabidopsis thaliana, including an α/ß core domain and ß-stranded small domain, with dimerization occurring via the ß-stranded small domain. In a binary complex, we show that Ser30 in the substrate-binding loop of the small domain mediates interactions with the hallmark N1 atom of pseudouridine nucleobase, causing conformational changes in its quaternary structure. Kinetic and fluorescence spectroscopic analyses also showed that the Ser30-mediated interaction is a prerequisite for conformational changes and subsequent catalysis by EcPUKI. Furthermore, S30A mutation or EcPUKI complexed with other nucleosides homologous to pseudouridine but lacking the pseudouridine-specific N1 atom did not induce such conformational changes, demonstrating the catalytic significance of the proposed Ser30-mediated interaction. These analyses provide structural and functional evidence for a pseudouridine-dependent conformational change of EcPUKI and its functional linkage to catalysis.

Genetic assessment of pathogenic germline alterations in lysosomal genes among Asian patients with pancreatic ductal adenocarcinoma.

BACKGROUND: Lysosomes are closely linked to autophagic activity, which plays a vital role in pancreatic ductal adenocarcinoma (PDAC) biology. The survival of PDAC patients is still poor, and the identification of novel genetic factors for prognosis and treatment is highly required to prevent PDAC-related deaths. This study investigated the germline variants related to lysosomal dysfunction in patients with PDAC and to analyze whether they contribute to the development of PDAC. METHODS: The germline putative pathogenic variants (PPV) in genes involved in lysosomal storage disease (LSD) was compared between patients with PDAC (n = 418) and healthy controls (n = 845) using targeted panel and whole-exome sequencing. Furthermore, pancreatic organoids from wild-type and KrasG12D mice were used to evaluate the effect of lysosomal dysfunction on PDAC development. RNA sequencing (RNA-seq) analysis was performed with established PDAC patient-derived organoids (PDOs) according to the PPV status. RESULTS: The PPV in LSD-related genes was higher in patients with PDAC than in healthy controls (8.13 vs. 4.26%, Log2 OR = 1.65, P = 3.08 × 10-3). The PPV carriers of LSD-related genes with PDAC were significantly younger than the non-carriers (mean age 61.5 vs. 65.3 years, P = 0.031). We further studied a variant of the lysosomal enzyme, galactosylceramidase (GALC), which was the most frequently detected LSD variant in our cohort. Autophagolysosomal activity was hampered when GALC was downregulated, which was accompanied by paradoxically elevated autophagic flux. Furthermore, the number of proliferating Ki-67+ cells increased significantly in pancreatic organoids derived from Galc knockout KrasG12D mice. Moreover, GALC PPV carriers tended to show drug resistance in both PDAC cell line and PDAC PDO, and RNA-seq analysis revealed that various metabolism and gene repair pathways were upregulated in PDAC PDOs harboring a GALC variant. CONCLUSIONS: Genetically defined lysosomal dysfunction is frequently observed in patients with young-onset PDAC. This might contribute to PDAC development by altering metabolism and impairing autophagolysosomal activity, which could be potentially implicated in therapeutic applications for PDAC.

Assessment of attenuation of varicella-zoster virus vaccines based on genomic comparison.

Live attenuated varicella-zoster virus (VZV) vaccines are used to prevent chickenpox and shingles. Single nucleotide polymorphisms (SNPs) that occur during the attenuation of parental strains are critical indicators of vaccine safety. To assess the attenuation of commercial VZV vaccines, genetic variants were comprehensively examined through high-throughput sequencing of viral DNA isolated from four VZV vaccines (Barycela, VarilRix, VariVax, and SKY Varicella). Whole-genome comparison of the four vaccines with the wild-type strain (Dumas) revealed that the sequences are highly conserved on a genome-wide scale. Among the 196 common variants across the four vaccines, 195 were already present in the genome of the parental strain (pOka), indicating that the variants occurred during the generation of the parental strain from the Dumas strain. Compared to the pOka genome, the vaccines exhibited distinct variant frequencies on a genome-wide and within an attenuation-related open reading frame. In particular, attenuation-associated 42 SNPs showed that Barycela, VarilRix, VariVax, and SKY Varicella are in ascending order regarding similarity with pOka-like genotypes, which in turn, might provide genomic evidence for the levels of attenuation. Finally, the phylogenetic network analysis demonstrated that genetic distances from the parental strain correlated with the attenuation levels of the vaccines.

Determinants governing T cell receptor α/ß-chain pairing in repertoire formation of identical twins.

The T cell repertoire in each individual includes T cell receptors (TCRs) of enormous sequence diversity through the pairing of diverse TCR α- and ß-chains, each generated by somatic recombination of paralogous gene segments. Whether the TCR repertoire contributes to susceptibility to infectious or autoimmune diseases in concert with disease-associated major histocompatibility complex (MHC) polymorphisms is unknown. Due to a lack in high-throughput technologies to sequence TCR α-ß pairs, current studies on whether the TCR repertoire is shaped by host genetics have so far relied only on single-chain analysis. Using a high-throughput single T cell sequencing technology, we obtained the largest paired TCRαß dataset so far, comprising 965,523 clonotypes from 15 healthy individuals including 6 monozygotic twin pairs. Public TCR α- and, to a lesser extent, TCR ß-chain sequences were common in all individuals. In contrast, sharing of entirely identical TCRαß amino acid sequences was very infrequent in unrelated individuals, but highly increased in twins, in particular in CD4 memory T cells. Based on nucleotide sequence identity, a subset of these shared clonotypes appeared to be the progeny of T cells that had been generated during fetal development and had persisted for more than 50 y. Additional shared TCRαß in twins were encoded by different nucleotide sequences, implying that genetic determinants impose structural constraints on thymic selection that favor the selection of TCR α-ß pairs with entire sequence identities.

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.

The chromatin-binding protein PHF6 functions as an E3 ubiquitin ligase of H2BK120 via H2BK12Ac recognition for activation of trophectodermal genes.

Epigenetic regulation is important for establishing lineage-specific gene expression during early development. Although signaling pathways have been well-studied for regulation of trophectoderm reprogramming, epigenetic regulation of trophectodermal genes with histone modification dynamics have been poorly understood. Here, we identify that plant homeodomain finger protein 6 (PHF6) is a key epigenetic regulator for activation of trophectodermal genes using RNA-sequencing and ChIP assays. PHF6 acts as an E3 ubiquitin ligase for ubiquitination of H2BK120 (H2BK120ub) via its extended plant homeodomain 1 (PHD1), while the extended PHD2 of PHF6 recognizes acetylation of H2BK12 (H2BK12Ac). Intriguingly, the recognition of H2BK12Ac by PHF6 is important for exerting its E3 ubiquitin ligase activity for H2BK120ub. Together, our data provide evidence that PHF6 is crucial for epigenetic regulation of trophectodermal gene expression by linking H2BK12Ac to H2BK120ub modification.

RORα is crucial for attenuated inflammatory response to maintain intestinal homeostasis.

Retinoic acid-related orphan receptor α (RORα) functions as a transcription factor for various biological processes, including circadian rhythm, cancer, and metabolism. Here, we generate intestinal epithelial cell (IEC)-specific RORα-deficient (RORαΔIEC) mice and find that RORα is crucial for maintaining intestinal homeostasis by attenuating nuclear factor κB (NF-κB) transcriptional activity. RORαΔIEC mice exhibit excessive intestinal inflammation and highly activated inflammatory responses in the dextran sulfate sodium (DSS) mouse colitis model. Transcriptome analysis reveals that deletion of RORα leads to up-regulation of NF-κB target genes in IECs. Chromatin immunoprecipitation analysis reveals corecruitment of RORα and histone deacetylase 3 (HDAC3) on NF-κB target promoters and subsequent dismissal of CREB binding protein (CBP) and bromodomain-containing protein 4 (BRD4) for transcriptional repression. Together, we demonstrate that RORα/HDAC3-mediated attenuation of NF-κB signaling controls the balance of inflammatory responses, and therapeutic strategies targeting this epigenetic regulation could be beneficial to the treatment of chronic inflammatory diseases, including inflammatory bowel disease (IBD).

Poly(A)-specific ribonuclease sculpts the 3' ends of microRNAs.

The 3' ends of metazoan microRNAs (miRNAs) are initially defined by the RNase III enzymes during maturation, but subsequently experience extensive modifications by several enzymatic activities. For example, terminal nucleotidyltransferases (TENTs) elongate miRNAs by adding one or a few nucleotides to their 3' ends, which occasionally leads to differential regulation of miRNA stability or function. However, the catalytic entities that shorten miRNAs and the molecular consequences of such shortening are less well understood, especially in vertebrates. Here, we report that poly(A)-specific ribonuclease (PARN) sculpts the 3' ends of miRNAs in human cells. By generating PARN knockout cells and characterizing their miRNAome, we demonstrate that PARN digests the 3' extensions of miRNAs that are derived from the genome or attached by TENTs, thereby effectively reducing the length of miRNAs. Surprisingly, PARN-mediated shortening has little impact on miRNA stability, suggesting that this process likely operates to finalize miRNA maturation, rather than to initiate miRNA decay. PARN-mediated shortening is pervasive across most miRNAs and appears to be a conserved mechanism contributing to the 3' end formation of vertebrate miRNAs. Our findings add miRNAs to the expanding list of noncoding RNAs whose 3' end formation depends on PARN.

Selection of self-priming molecular replicators.

Self-priming amplification of oligonucleotides is possible based on foldback of 3' ends, self-priming, and concatemerization, especially in the presence of phosphorothioate linkages. Such a simple replicative mechanism may have led to the accumulation of specific replicators at or near the origin of life. To determine how early replicators may have competed with one another, we have carried out selections with phosphorothiolated hairpins appended to a short random sequence library (N10). Upon the addition of deoxynucleoside triphosphates and a polymerase, concatemers quickly formed, and those random sequences that templated the insertion of purines, especially during initiation, quickly predominated. Over several serial transfers, particular sequences accumulated, and in isolation these were shown to outcompete less efficient replicators.

CTCF cooperates with CtIP to drive homologous recombination repair of double-strand breaks.

The pleiotropic CCCTC-binding factor (CTCF) plays a role in homologous recombination (HR) repair of DNA double-strand breaks (DSBs). However, the precise mechanistic role of CTCF in HR remains largely unclear. Here, we show that CTCF engages in DNA end resection, which is the initial, crucial step in HR, through its interactions with MRE11 and CtIP. Depletion of CTCF profoundly impairs HR and attenuates CtIP recruitment at DSBs. CTCF physically interacts with MRE11 and CtIP and promotes CtIP recruitment to sites of DNA damage. Subsequently, CTCF facilitates DNA end resection to allow HR, in conjunction with MRE11-CtIP. Notably, the zinc finger domain of CTCF binds to both MRE11 and CtIP and enables proficient CtIP recruitment, DNA end resection and HR. The N-terminus of CTCF is able to bind to only MRE11 and its C-terminus is incapable of binding to MRE11 and CtIP, thereby resulting in compromised CtIP recruitment, DSB resection and HR. Overall, this suggests an important function of CTCF in DNA end resection through the recruitment of CtIP at DSBs. Collectively, our findings identify a critical role of CTCF at the first control point in selecting the HR repair pathway.

iNrich, Rapid and Robust Method to Enrich N-Terminal Proteome in a Highly Multiplexed Platform.

The field of terminal proteomics is limited in that it is optimized for large-scale analysis via multistep processes involving liquid chromatography. Here, we present an integrated N-terminal peptide enrichment method (iNrich) that can handle as little as 25 µg of cell lysate via a single-stage encapsulated solid-phase extraction column. iNrich enables simple, rapid, and reproducible sample processing, treatment of a wide range of protein amounts (25 µg ∼ 1 mg), multiplexed parallel sample preparation, and in-stage sample prefractionation using a mixed-anion-exchange filter. We identified ∼5000 N-terminal peptides (Nt-peptides) from only 100 µg of human cell lysate including Nt-formyl peptides. Multiplexed sample preparation facilitated quantitative and robust enrichment of N-terminome with dozens of samples simultaneously. We further developed the method to incorporate isobaric tags such as a tandem mass tag (TMT) and used it to discover novel peptides during ER stress analysis. The iNrich facilitated high-throughput N-terminomics and degradomics at a low cost using commercially available reagents and apparatus, without requiring arduous procedures.

The ATP-dependent chromatin remodeler Chd1 is recruited by transcription elongation factors and maintains H3K4me3/H3K36me3 domains at actively transcribed and spliced genes.

Chd1 (Chromodomain Helicase DNA Binding Protein 1) is a conserved ATP-dependent chromatin remodeler that maintains the nucleosomal structure of chromatin, but the determinants of its specificity and its impact on gene expression are not well defined. To identify the determinants of Chd1 binding specificity in the yeast genome, we investigated Chd1 occupancy in mutants of several candidate factors. We found that several components of the PAF1 transcription elongation complex contribute to Chd1 recruitment to highly transcribed genes and identified Spt4 as a factor that appears to negatively modulate Chd1 binding to chromatin. We discovered that CHD1 loss alters H3K4me3 and H3K36me3 patterns throughout the yeast genome. Interestingly, the aberrant histone H3 methylation patterns were predominantly observed within 1 kb from the transcription start site, where both histone H3 methylation marks co-occur. A reciprocal change between the two marks was obvious in the absence of Chd1, suggesting a role for CHD1 in establishing or maintaining the boundaries of these largely mutually exclusive histone marks. Strikingly, intron-containing genes were most susceptible to CHD1 loss and exhibited a high degree of histone H3 methylation changes. Intron retention was significantly lower in the absence of CHD1, suggesting that CHD1 function as a chromatin remodeler could indirectly affect RNA splicing.

Large-scale sequence and structural comparisons of human naive and antigen-experienced antibody repertoires.

Elucidating how antigen exposure and selection shape the human antibody repertoire is fundamental to our understanding of B-cell immunity. We sequenced the paired heavy- and light-chain variable regions (VH and VL, respectively) from large populations of single B cells combined with computational modeling of antibody structures to evaluate sequence and structural features of human antibody repertoires at unprecedented depth. Analysis of a dataset comprising 55,000 antibody clusters from CD19(+)CD20(+)CD27(-) IgM-naive B cells, >120,000 antibody clusters from CD19(+)CD20(+)CD27(+) antigen-experienced B cells, and >2,000 RosettaAntibody-predicted structural models across three healthy donors led to a number of key findings: (i) VH and VL gene sequences pair in a combinatorial fashion without detectable pairing restrictions at the population level; (ii) certain VH:VL gene pairs were significantly enriched or depleted in the antigen-experienced repertoire relative to the naive repertoire; (iii) antigen selection increased antibody paratope net charge and solvent-accessible surface area; and (iv) public heavy-chain third complementarity-determining region (CDR-H3) antibodies in the antigen-experienced repertoire showed signs of convergent paired light-chain genetic signatures, including shared light-chain third complementarity-determining region (CDR-L3) amino acid sequences and/or Vκ,λ-Jκ,λ genes. The data reported here address several longstanding questions regarding antibody repertoire selection and development and provide a benchmark for future repertoire-scale analyses of antibody responses to vaccination and disease.

Subtype-specific addiction of the activated B-cell subset of diffuse large B-cell lymphoma to FOXP1.

High expression of the forkhead box P1 (FOXP1) transcription factor distinguishes the aggressive activated B cell (ABC) diffuse large B-cell lymphoma (DLBCL) subtype from the better prognosis germinal center B-cell (GCB)-DLBCL subtype and is highly correlated with poor outcomes. A genetic or functional role for FOXP1 in lymphomagenesis, however, remains unknown. Here, we report that sustained FOXP1 expression is vital for ABC-DLBCL cell-line survival. Genome-wide analyses revealed direct and indirect FOXP1 transcriptional enforcement of ABC-DLBCL hallmarks, including the classical NF-κB and MYD88 (myeloid differentiation primary response gene 88) pathways. FOXP1 promoted gene expression underlying transition of the GCB cell to the plasmablast--the transient B-cell stage targeted in ABC-DLBCL transformation--by antagonizing pathways distinctive of GCB-DLBCL, including that of the GCB "master regulator," BCL6 (B-cell lymphoma 6). Cell-line derived FOXP1 target genes that were highly correlated with FOXP1 expression in primary DLBCL accurately segregated the corresponding clinical subtypes of a large cohort of primary DLBCL isolates and identified conserved pathways associated with ABC-DLBCL pathology.

Simultaneous mapping of transcript ends at single-nucleotide resolution and identification of widespread promoter-associated non-coding RNA governed by TATA elements.

Understanding the relationships between regulatory factor binding, chromatin structure, cis-regulatory elements and RNA-regulation mechanisms relies on accurate information about transcription start sites (TSS) and polyadenylation sites (PAS). Although several approaches have identified transcript ends in yeast, limitations of resolution and coverage have remained, and definitive identification of TSS and PAS with single-nucleotide resolution has not yet been achieved. We developed SMORE-seq (simultaneous mapping of RNA ends by sequencing) and used it to simultaneously identify the strongest TSS for 5207 (90%) genes and PAS for 5277 (91%) genes. The new transcript annotations identified by SMORE-seq showed improved distance relationships with TATA-like regulatory elements, nucleosome positions and active RNA polymerase. We found 150 genes whose TSS were downstream of the annotated start codon, and additional analysis of evolutionary conservation and ribosome footprinting suggests that these protein-coding sequences are likely to be mis-annotated. SMORE-seq detected short non-coding RNAs transcribed divergently from more than a thousand promoters in wild-type cells under normal conditions. These divergent non-coding RNAs were less evident at promoters containing canonical TATA boxes, suggesting a model where transcription initiation at promoters by RNAPII is bidirectional, with TATA elements serving to constrain the directionality of initiation.

Context-dependent genomic locus effects on antibody production in recombinant Chinese hamster ovary cells generated through random integration.

High-yield production of therapeutic protein using Chinese hamster ovary (CHO) cells requires stable cell line development (CLD). CLD typically uses random integration of transgenes; however, this results in clonal variation and subsequent laborious clone screening. Therefore, site-specific integration of a protein expression cassette into a desired chromosomal locus showing high transcriptional activity and stability, referred to as a hot spot, is emerging. Although positional effects are important for therapeutic protein expression, the sequence-specific mechanisms by which hotspots work are not well understood. In this study, we performed whole-genome sequencing (WGS) to locate randomly inserted vectors in the genome of recombinant CHO cells expressing high levels of monoclonal antibodies (mAbs) and experimentally validated these locations and vector compositions. The integration site was characterized by active histone marks and potential enhancer activities, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) mediated indel mutations in the region upstream of the integration site led to a significant reduction in specific antibody productivity by up to 30%. Notably, the integration site and its core region did not function equivalently outside the native genomic context, showing a minimal effect on the increase in exogenous protein expression in the host cell line. We also observed a superior production capacity of the mAb expressing cell line compared to that of the host cell line. Collectively, this study demonstrates that developing recombinant CHO cell lines to produce therapeutic proteins at high levels requires a balance of factors including transgene configuration, genomic locus landscape, and host cell properties.

Proteogenomic landscape of human pancreatic ductal adenocarcinoma in an Asian population reveals tumor cell-enriched and immune-rich subtypes.

We report a proteogenomic analysis of pancreatic ductal adenocarcinoma (PDAC). Mutation-phosphorylation correlations identified signaling pathways associated with somatic mutations in significantly mutated genes. Messenger RNA-protein abundance correlations revealed potential prognostic biomarkers correlated with patient survival. Integrated clustering of mRNA, protein and phosphorylation data identified six PDAC subtypes. Cellular pathways represented by mRNA and protein signatures, defining the subtypes and compositions of cell types in the subtypes, characterized them as classical progenitor (TS1), squamous (TS2-4), immunogenic progenitor (IS1) and exocrine-like (IS2) subtypes. Compared with the mRNA data, protein and phosphorylation data further classified the squamous subtypes into activated stroma-enriched (TS2), invasive (TS3) and invasive-proliferative (TS4) squamous subtypes. Orthotopic mouse PDAC models revealed a higher number of pro-tumorigenic immune cells in TS4, inhibiting T cell proliferation. Our proteogenomic analysis provides significantly mutated genes/biomarkers, cellular pathways and cell types as potential therapeutic targets to improve stratification of patients with PDAC.

Fast Growing Furious Races for Targeting Fibroblast Growth Factor Receptors.

Targeting fibroblast growth factor receptors (FGFRs) has been slow compared to other targeted cancer therapies for receptor tyrosine kinases, such as epidermal growth factor receptors. The low efficacy and variable response have limited the growth of FGFR inhibitors in clinical use. Nevertheless, recent systematic and genomic approaches have identified the biological conditions for effectively targeting FGFRs and can accelerate the development of targeted drugs. Under clinical and preclinical trials, the inhibitors started fast growing furious races to target FGFRs. Finally, FGFRs will be more actionable and targetable with more precise and effective drugs at the end of the race, passing the finish line.