CRISPR-free, programmable RNA pseudouridylation to suppress premature termination codons.

Nonsense mutations, accounting for >20% of disease-associated mutations, lead to premature translation termination. Replacing uridine with pseudouridine in stop codons suppresses translation termination, which could be harnessed to mediate readthrough of premature termination codons (PTCs). Here, we present RESTART, a programmable RNA base editor, to revert PTC-induced translation termination in mammalian cells. RESTART utilizes an engineered guide snoRNA (gsnoRNA) and the endogenous H/ACA box snoRNP machinery to achieve precise pseudouridylation. We also identified and optimized gsnoRNA scaffolds to increase the editing efficiency. Unexpectedly, we found that a minor isoform of pseudouridine synthase DKC1, lacking a C-terminal nuclear localization signal, greatly improved the PTC-readthrough efficiency. Although RESTART induced restricted off-target pseudouridylation, they did not change the coding information nor the expression level of off-targets. Finally, RESTART enables robust pseudouridylation in primary cells and achieves functional PTC readthrough in disease-relevant contexts. Collectively, RESTART is a promising RNA-editing tool for research and therapeutics.

Deep whole-genome analysis of 494 hepatocellular carcinomas.

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

An interactive resource to identify cancer genetic and lineage dependencies targeted by small molecules.

The high rate of clinical response to protein-kinase-targeting drugs matched to cancer patients with specific genomic alterations has prompted efforts to use cancer cell line (CCL) profiling to identify additional biomarkers of small-molecule sensitivities. We have quantitatively measured the sensitivity of 242 genomically characterized CCLs to an Informer Set of 354 small molecules that target many nodes in cell circuitry, uncovering protein dependencies that: (1) associate with specific cancer-genomic alterations and (2) can be targeted by small molecules. We have created the Cancer Therapeutics Response Portal (http://www.broadinstitute.org/ctrp) to enable users to correlate genetic features to sensitivity in individual lineages and control for confounding factors of CCL profiling. We report a candidate dependency, associating activating mutations in the oncogene ß-catenin with sensitivity to the Bcl-2 family antagonist, navitoclax. The resource can be used to develop novel therapeutic hypotheses and to accelerate discovery of drugs matched to patients by their cancer genotype and lineage.

Paradoxical Mitophagy Regulation by PINK1 and TUFm.

Aberrant mitophagy has been implicated in a broad spectrum of disorders. PINK1, Parkin, and ubiquitin have pivotal roles in priming mitophagy. However, the entire regulatory landscape and the precise control mechanisms of mitophagy remain to be elucidated. Here, we uncover fundamental mitophagy regulation involving PINK1 and a non-canonical role of the mitochondrial Tu translation elongation factor (TUFm). The mitochondrion-cytosol dual-localized TUFm interacts with PINK1 biochemically and genetically, which is an evolutionarily conserved Parkin-independent route toward mitophagy. A PINK1-dependent TUFm phosphoswitch at Ser222 determines conversion from activating to suppressing mitophagy. PINK1 modulates differential translocation of TUFm because p-S222-TUFm is restricted predominantly to the cytosol, where it inhibits mitophagy by impeding Atg5-Atg12 formation. The self-antagonizing feature of PINK1/TUFm is critical for the robustness of mitophagy regulation, achieved by the unique kinetic parameters of p-S222-TUFm, p-S65-ubiquitin, and their common kinase PINK1. Our findings provide new mechanistic insights into mitophagy and mitophagy-associated disorders.

Hypothetical chloroplast reading frame 51 encodes a photosystem I assembly factor in cyanobacteria.

Hypothetical chloroplast open reading frames (ycfs) are putative genes in the plastid genomes of photosynthetic eukaryotes. Many ycfs are also conserved in the genomes of cyanobacteria, the presumptive ancestors of present-day chloroplasts. The functions of many ycfs are still unknown. Here, we generated knock-out mutants for ycf51 (sll1702) in the cyanobacterium Synechocystis sp. PCC 6803. The mutants showed reduced photoautotrophic growth due to impaired electron transport between photosystem II (PSII) and PSI. This phenotype results from greatly reduced PSI content in the ycf51 mutant. The ycf51 disruption had little effect on the transcription of genes encoding photosynthetic complex components and the stabilization of the PSI complex. In vitro and in vivo analyses demonstrated that Ycf51 cooperates with PSI assembly factor Ycf3 to mediate PSI assembly. Furthermore, Ycf51 interacts with the PSI subunit PsaC. Together with its specific localization in the thylakoid membrane and the stromal exposure of its hydrophilic region, our data suggest that Ycf51 is involved in PSI complex assembly. Ycf51 is conserved in all sequenced cyanobacteria, including the earliest branching cyanobacteria of the Gloeobacter genus, and is also present in the plastid genomes of glaucophytes. However, Ycf51 has been lost from other photosynthetic eukaryotic lineages. Thus, Ycf51 is a PSI assembly factor that has been functionally replaced during the evolution of oxygenic photosynthetic eukaryotes.

High UV damage and low repair, but not cytosine deamination, stimulate mutation hotspots at ETS binding sites in melanoma.

Noncoding mutation hotspots have been identified in melanoma and many of them occur at the binding sites of E26 transformation-specific (ETS) proteins; however, their formation mechanism and functional impacts are not fully understood. Here, we used UV (Ultraviolet) damage sequencing data and analyzed cyclobutane pyrimidine dimer (CPD) formation, DNA repair, and CPD deamination in human cells at single-nucleotide resolution. Our data show prominent CPD hotspots immediately after UV irradiation at ETS binding sites, particularly at sites with a conserved TTCCGG motif, which correlate with mutation hotspots identified in cutaneous melanoma. Additionally, CPDs are repaired slower at ETS binding sites than in flanking DNA. Cytosine deamination in CPDs to uracil is suggested as an important step for UV mutagenesis. However, we found that CPD deamination is significantly suppressed at ETS binding sites, particularly for the CPD hotspot on the 5' side of the ETS motif, arguing against a role for CPD deamination in promoting ETS-associated UV mutations. Finally, we analyzed a subset of frequently mutated promoters, including the ribosomal protein genes RPL13A and RPS20, and found that mutations in the ETS motif can significantly reduce the promoter activity. Thus, our data identify high UV damage and low repair, but not CPD deamination, as the main mechanism for ETS-associated mutations in melanoma and uncover important roles of often-overlooked mutation hotspots in perturbing gene transcription.

eccDNA-pipe: an integrated pipeline for identification, analysis and visualization of extrachromosomal circular DNA from high-throughput sequencing data.

Extrachromosomal circular DNA (eccDNA) is currently attracting considerable attention from researchers due to its significant impact on tumor biogenesis. High-throughput sequencing (HTS) methods for eccDNA identification are continually evolving. However, an efficient pipeline for the integrative and comprehensive analysis of eccDNA obtained from HTS data is still lacking. Here, we introduce eccDNA-pipe, an accessible software package that offers a user-friendly pipeline for conducting eccDNA analysis starting from raw sequencing data. This dataset includes data from various sequencing techniques such as whole-genome sequencing (WGS), Circle-seq and Circulome-seq, obtained through short-read sequencing or long-read sequencing. eccDNA-pipe presents a comprehensive solution for both upstream and downstream analysis, encompassing quality control and eccDNA identification in upstream analysis and downstream tasks such as eccDNA length distribution analysis, differential analysis of genes enriched with eccDNA and visualization of eccDNA structures. Notably, eccDNA-pipe automatically generates high-quality publication-ready plots. In summary, eccDNA-pipe provides a comprehensive and user-friendly pipeline for customized analysis of eccDNA research.

The AAA-ATPase Yta4/ATAD1 interacts with the mitochondrial divisome to inhibit mitochondrial fission.

Mitochondria are in a constant balance of fusion and fission. Excessive fission or deficient fusion leads to mitochondrial fragmentation, causing mitochondrial dysfunction and physiological disorders. How the cell prevents excessive fission of mitochondria is not well understood. Here, we report that the fission yeast AAA-ATPase Yta4, which is the homolog of budding yeast Msp1 responsible for clearing mistargeted tail-anchored (TA) proteins on mitochondria, plays a critical role in preventing excessive mitochondrial fission. The absence of Yta4 leads to mild mitochondrial fragmentation in a Dnm1-dependent manner but severe mitochondrial fragmentation upon induction of mitochondrial depolarization. Overexpression of Yta4 delocalizes the receptor proteins of Dnm1, i.e., Fis1 (a TA protein) and Mdv1 (the bridging protein between Fis1 and Dnm1), from mitochondria and reduces the localization of Dnm1 to mitochondria. The effect of Yta4 overexpression on Fis1 and Mdv1, but not Dnm1, depends on the ATPase and translocase activities of Yta4. Moreover, Yta4 interacts with Dnm1, Mdv1, and Fis1. In addition, Yta4 competes with Dnm1 for binding Mdv1 and decreases the affinity of Dnm1 for GTP and inhibits Dnm1 assembly in vitro. These findings suggest a model, in which Yta4 inhibits mitochondrial fission by inhibiting the function of the mitochondrial divisome composed of Fis1, Mdv1, and Dnm1. Therefore, the present work reveals an uncharacterized molecular mechanism underlying the inhibition of mitochondrial fission.

Author Correction: The Apostasia genome and the evolution of orchids.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Molecular mechanism of specific DNA sequence recognition by NRF1.

Nuclear respiratory factor 1 (NRF1) regulates the expression of genes that are vital for mitochondrial biogenesis, respiration, and various other cellular processes. While NRF1 has been reported to bind specifically to GC-rich promoters as a homodimer, the precise molecular mechanism governing its recognition of target gene promoters has remained elusive. To unravel the recognition mechanism, we have determined the crystal structure of the NRF1 homodimer bound to an ATGCGCATGCGCAT dsDNA. In this complex, NRF1 utilizes a flexible linker to connect its dimerization domain (DD) and DNA binding domain (DBD). This configuration allows one NRF1 monomer to adopt a U-turn conformation, facilitating the homodimer to specifically bind to the two TGCGC motifs in the GCGCATGCGC consensus sequence from opposite directions. Strikingly, while the NRF1 DBD alone could also bind to the half-site (TGCGC) DNA of the consensus sequence, the cooperativity between DD and DBD is essential for the binding of the intact GCGCATGCGC sequence and the transcriptional activity of NRF1. Taken together, our results elucidate the molecular mechanism by which NRF1 recognizes specific DNA sequences in the promoters to regulate gene expression.

Genomic analysis reveals a cryptic pangolin species.

Eight extant species of pangolins are currently recognized. Recent studies found that two mitochondrial haplotypes identified in confiscations in Hong Kong could not be assigned to any known pangolin species, implying the existence of a species. Here, we report that two additional mitochondrial haplotypes identified in independent confiscations from Yunnan align with the putative species haplotypes supporting the existence of this mysterious species/population. To verify the new species scenario we performed a comprehensive analysis of scale characteristics and 138 whole genomes representing all recognized pangolin species and the cryptic new species, 98 of which were generated here. Our morphometric results clearly attributed this cryptic species to Asian pangolins (Manis sp.) and the genomic data provide robust and compelling evidence that it is a pangolin species distinct from those recognized previously, which separated from the Philippine pangolin and Malayan pangolin over 5 Mya. Our study provides a solid genomic basis for its formal recognition as the ninth pangolin species or the fifth Asian one, supporting a new taxonomic classification of pangolins. The effects of glacial climate changes and recent anthropogenic activities driven by illegal trade are inferred to have caused its population decline with the genomic signatures showing low genetic diversity, a high level of inbreeding, and high genetic load. Our finding greatly expands current knowledge of pangolin diversity and evolution and has vital implications for conservation efforts to prevent the extinction of this enigmatic and endangered species from the wild.

Structural basis for specific DNA sequence recognition by the transcription factor NFIL3.

The CCAAT/enhancer-binding proteins (C/EBPs) constitute a family of pivotal transcription factors involved in tissue development, cellular function, proliferation, and differentiation. NFIL3, as one of them, plays an important role in regulating immune cell differentiation, circadian clock system, and neural regeneration, yet its specific DNA recognition mechanism remains enigmatic. In this study, we showed by the ITC binding experiments that NFIL3 prefers to bind to the TTACGTAA DNA motif. Our structural studies revealed that the α-helical NFIL3 bZIP domain dimerizes through its leucine zipper region, and binds to DNA via its basic region. The two basic regions of the NFIL3 bZIP dimer were pushed apart upon binding to DNA, facilitating the snug accommodation of the two basic regions within the major grooves of the DNA. Remarkably, our binding and structural data also revealed that both NFIL3 and C/EBPα/ß demonstrate a shared preference for the TTACGTAA sequence. Furthermore, our study revealed that disease-associated mutations within the NFIL3 bZIP domain result in either reduction or complete disruption of its DNA binding ability. These discoveries not only provide valuable insights into the DNA binding mechanisms of NFIL3 but also elucidate the causal role of NFIL3 mutations in disease pathogenesis.

Benchmarking spatial and single-cell transcriptomics integration methods for transcript distribution prediction and cell type deconvolution.

Spatial transcriptomics approaches have substantially advanced our capacity to detect the spatial distribution of RNA transcripts in tissues, yet it remains challenging to characterize whole-transcriptome-level data for single cells in space. Addressing this need, researchers have developed integration methods to combine spatial transcriptomic data with single-cell RNA-seq data to predict the spatial distribution of undetected transcripts and/or perform cell type deconvolution of spots in histological sections. However, to date, no independent studies have comparatively analyzed these integration methods to benchmark their performance. Here we present benchmarking of 16 integration methods using 45 paired datasets (comprising both spatial transcriptomics and scRNA-seq data) and 32 simulated datasets. We found that Tangram, gimVI, and SpaGE outperformed other integration methods for predicting the spatial distribution of RNA transcripts, whereas Cell2location, SpatialDWLS, and RCTD are the top-performing methods for the cell type deconvolution of spots. We provide a benchmark pipeline to help researchers select optimal integration methods to process their datasets.

Histone tyrosine sulfation by SULT1B1 regulates H4R3me2a and gene transcription.

Tyrosine sulfation is a common posttranslational modification in mammals. To date, it has been thought to be limited to secreted and transmembrane proteins, but little is known about tyrosine sulfation on nuclear proteins. Here we report that SULT1B1 is a histone sulfotransferase that can sulfate the tyrosine 99 residue of nascent histone H3 in cytosol. The sulfated histone H3 can be transported into the nucleus and majorly deposited in the promoter regions of genes in chromatin. While the H3Y99 residue is buried inside octameric nucleosome, dynamically regulated subnucleosomal structures provide chromatin-H3Y99sulf the opportunity of being recognized and bound by PRMT1, which deposits H4R3me2a in chromatin. Disruption of H3Y99sulf reduces PRMT1 binding to chromatin, H4R3me2a level and gene transcription. These findings reveal the mechanisms underlying H3Y99 sulfation and its cross-talk with H4R3me2a to regulate gene transcription. This study extends the spectrum of tyrosine sulfation on nuclear proteins and the repertoire of histone modifications regulating chromatin functions.

DDX58 variant triggers IFN-ß-induced autophagy in trabecular meshwork and influences intraocular pressure.

Singleton-Merten syndrome (SMS) is a rare immunogenetic disorder affecting multiple systems, characterized by dental dysplasia, aortic calcification, glaucoma, skeletal abnormalities, and psoriasis. Glaucoma, a key feature of both classical and atypical SMS, remains poorly understood in terms of its molecular mechanism caused by DDX58 mutation. This study presented a novel DDX58 variant (c.1649A>C [p.Asp550Ala]) in a family with childhood glaucoma. Functional analysis showed that DDX58 variant caused an increase in IFN-stimulated gene expression and high IFN-ß-based type-I IFN. As the trabecular meshwork (TM) is responsible for controlling intraocular pressure (IOP), we examine the effect of IFN-ß on TM cells. Our study is the first to demonstrate that IFN-ß significantly reduced TM cell viability and function by activating autophagy. In addition, anterior chamber injection of IFN-ß remarkably increased IOP level in mice, which can be attenuated by treatments with autophagy inhibitor chloroquine. To uncover the specific mechanism underlying IFN-ß-induced autophagy in TM cells, we performed microarray analysis in IFN-ß-treated and DDX58 p.Asp550Ala TM cells. It showed that RSAD2 is necessary for IFN-ß-induced autophagy. Knockdown of RSAD2 by siRNA significantly decreased autophagy flux induced by IFN-ß. Our findings suggest that DDX58 mutation leads to the overproduction of IFN-ß, which elevates IOP by modulating autophagy through RSAD2 in TM cells.

ICBP90, an epigenetic regulator, induces DKK3 promoter methylation, promotes glioma progression, and reduces sensitivity to cis-platinum.

Glioma is the most common brain malignancy, characterized by high morbidity, high mortality, and treatment-resistance. Inverted CCAAT box Binding Protein of 90 kDa (ICBP90) has been reported to be involved in tumor progression and the maintenance of DNA methylation. Herein, we constructed ICBP90 over-expression and knockdown glioma cell lines, and found that ICBP90 knockdown inhibited glioma cell proliferation, migration, and invasion. ICBP90 silencing potentially enhanced cellular sensitivity to cis-platinum (DDP) and exacerbated DDP-induced pyroptosis, manifested by the elevated levels of gasdermin D-N-terminal and cleaved caspase 1; whereas, ICBP90 over-expression exhibited the opposite effects. Consistently, ICBP90 knockdown inhibited tumor growth in an in vivo mouse xenograft study using U251 cells stably expressing sh-ICBP90 and oe-ICBP90. Further experiments found that ICBP90 reduced the expression of Dickkopf 3 homolog (DKK3), a negative regulator of ß-catenin, by binding its promoter and inducing DNA methylation. ICBP90 knockdown prevented the nuclear translocation of ß-catenin and suppressed the expression of c-Myc and cyclin D1. Besides, DKK3 over-expression restored the effects of ICBP90 over-expression on cell proliferation, migration, invasion, and DDP sensitivity. Our findings suggest that ICBP90 inhibits the expression of DKK3 in glioma by maintaining DKK3 promoter methylation, thereby conducing to ICBP90-mediated carcinogenesis and drug insensitivity.

Structural basis for specific DNA sequence motif recognition by the TFAP2 transcription factors.

The TFAP2 family regulates gene expression during differentiation, development, and organogenesis, and includes five homologs in humans. They all possess a highly conserved DNA binding domain (DBD) followed by a helix-span-helix (HSH) domain. The DBD-HSH tandem domain specifically binds to a GCC(N3)GGC consensus sequence, but the precise recognition mechanisms remain unclear. Here, we found that TFAP2 preferred binding to the GCC(N3)GGC sequence, and the pseudo-palindromic GCC and GGC motifs and the length of the central spacer between the two motifs determined their binding specificity. Structural studies revealed that the two flat amphipathic α-helical HSH domains of TFAP2A stacked with each other to form a dimer via hydrophobic interactions, while the stabilized loops from both DBD domains inserted into two neighboring major grooves of the DNA duplex to form base-specific interactions. This specific DNA binding mechanism controlled the length of the central spacer and determined the DNA sequence specificity of TFAP2. Mutations of the TFAP2 proteins are implicated in various diseases. We illustrated that reduction or disruption of the DNA binding ability of the TFAP2 proteins is the primary cause of TFAP2 mutation-associated diseases. Thus, our findings also offer valuable insights into the pathogenesis of disease-associated mutations in TFAP2 proteins.

Mechanical cloak via data-driven aperiodic metamaterial design.

SignificanceAn invisibility cloak to conceal objects from an outside observer has long been a subject of interest in metamaterial design. While cloaks have been manufactured for optical, thermal, and electric fields, limited progress has been made for mechanical cloaks. Most existing designs rely on mapping-based methods, which have so far been limited to special base cells and a narrow selection of voids with simple shapes. In this study, we develop a fundamentally different approach by exploiting data-driven designs to offer timely, customized solutions to mechanical cloaking that were previously difficult to obtain. Through simulations and experimental validations, we show that excellent cloaking performance can be achieved for various boundary conditions, shapes of voids, base cells, and even multiple voids.

Structural insights into DNA recognition by the BEN domain of the transcription factor BANP.

The BEN domain-containing transcription factors regulate transcription by recruiting chromatin-modifying factors to specific chromatin regions via their DNA-binding BEN domains. The BEN domain of BANP has been shown to bind to a CGCG DNA sequence or an AAA-containing matrix attachment regions DNA sequence. Consistent with these in vivo observations, we identified an optimal DNA-binding sequence of AAATCTCG by protein binding microarray, which was also confirmed by our isothermal titration calorimetry and mutagenesis results. We then determined crystal structures of the BANP BEN domain in apo form and in complex with a CGCG-containing DNA, respectively, which revealed that the BANP BEN domain mainly used the electrostatic interactions to bind DNA with some base-specific interactions with the TC motifs. Our isothermal titration calorimetry results also showed that BANP bound to unmethylated and methylated DNAs with comparable binding affinities. Our complex structure of BANP-mCGCG revealed that the BANP BEN domain bound to the unmethylated and methylated DNAs in a similar mode and cytosine methylation did not get involved in binding, which is also consistent with our observations from the complex structures of the BEND6 BEN domain with the CGCG or CGmCG DNAs. Taken together, our results further elucidate the elements important for DNA recognition and transcriptional regulation by the BANP BEN domain-containing transcription factor.