Center for Life Sciences

Tremendous amount of RNA sequencing data have been produced by large consortium projects such as TCGA and GTEx, creating new opportunities for data mining and deeper understanding of gene functions. While certain existing web servers are valuable and widely used, many expression analysis functions needed by experimental biologists are still not adequately addressed by these tools. We introduce GEPIA (Gene Expression Profiling Interactive Analysis), a web-based tool to deliver fast and customizable functionalities based on TCGA and GTEx data. GEPIA provides key interactive and customizable functions including differential expression analysis, profiling plotting, correlation analysis, patient survival analysis, similar gene detection and dimensionality reduction analysis. The comprehensive expression analyses with simple clicking through GEPIA greatly facilitate data mining in wide research areas, scientific discussion and the therapeutic discovery process. GEPIA fills in the gap between cancer genomics big data and the delivery of integrated information to end users, thus helping unleash the value of the current data resources. GEPIA is available at http://gepia.cancer-pku.cn/.

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. For example, a key point that needs to be emphasized is thatthere is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process versus those that measure flux through the autophagy pathway (i.e., the completeprocess including the amount and rate of cargo sequestered and degraded). In particular, a block in macroautophagy that results in autophagosome accumulation must be differentiated from stimuli that increase autophagic activity, defined as increasedautophagy induction coupled with increased delivery to, and degradation within, lysosomes (inmost higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in manycases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. It is worth emphasizing here that lysosomal digestion is a stage of autophagy and evaluating its competence is a crucial part of the evaluation of autophagic flux, or complete autophagy. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as forreviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multipleassays to monitor autophagy. Along these lines, because of the potential for pleiotropic effects due to blocking autophagy through genetic manipulation, it is imperative to target by gene knockout or RNA interference more than one autophagyrelated protein. In addition, some individual Atg proteins, or groups of proteins, are involved in other cellular pathways implying that not all Atg proteins can be used as a specific marker for an autophagic process. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular assays, we hope to encourage technical innovation in the field.

How SARS-CoV-2 binds to human cells Scientists are racing to learn the secrets of severe acute respiratory syndrome–coronavirus 2 (SARS-CoV-2), which is the cause of the pandemic disease COVID-19. The first step in viral entry is the binding of the viral trimeric spike protein to the human receptor angiotensin-converting enzyme 2 (ACE2). Yan et al. present the structure of human ACE2 in complex with a membrane protein that it chaperones, B 0 AT1. In the context of this complex, ACE2 is a dimer. A further structure shows how the receptor binding domain of SARS-CoV-2 interacts with ACE2 and suggests that it is possible that two trimeric spike proteins bind to an ACE2 dimer. The structures provide a basis for the development of therapeutics targeting this crucial interaction. Science , this issue p. 1444

Introduced in 2017, the GEPIA (Gene Expression Profiling Interactive Analysis) web server has been a valuable and highly cited resource for gene expression analysis based on tumor and normal samples from the TCGA and the GTEx databases. Here, we present GEPIA2, an updated and enhanced version to provide insights with higher resolution and more functionalities. Featuring 198 619 isoforms and 84 cancer subtypes, GEPIA2 has extended gene expression quantification from the gene level to the transcript level, and supports analysis of a specific cancer subtype, and comparison between subtypes. In addition, GEPIA2 has adopted new analysis techniques of gene signature quantification inspired by single-cell sequencing studies, and provides customized analysis where users can upload their own RNA-seq data and compare them with TCGA and GTEx samples. We also offer an API for batch process and easy retrieval of the analysis results. The updated web server is publicly accessible at http://gepia2.cancer-pku.cn/.

Abstract Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature 1 . Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium 2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses 3–15 , enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated—but distinct—DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Abstract Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale 1–3 . Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4–5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter 4 ; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation 5,6 ; analyses timings and patterns of tumour evolution 7 ; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity 8,9 ; and evaluates a range of more-specialized features of cancer genomes 8,10–18 .

Immunotherapy has revolutionized cancer treatment and rejuvenated the field of tumor immunology. Several types of immunotherapy, including adoptive cell transfer (ACT) and immune checkpoint inhibitors (ICIs), have obtained durable clinical responses, but their efficacies vary, and only subsets of cancer patients can benefit from them. Immune infiltrates in the tumor microenvironment (TME) have been shown to play a key role in tumor development and will affect the clinical outcomes of cancer patients. Comprehensive profiling of tumor-infiltrating immune cells would shed light on the mechanisms of cancer-immune evasion, thus providing opportunities for the development of novel therapeutic strategies. However, the highly heterogeneous and dynamic nature of the TME impedes the precise dissection of intratumoral immune cells. With recent advances in single-cell technologies such as single-cell RNA sequencing (scRNA-seq) and mass cytometry, systematic interrogation of the TME is feasible and will provide insights into the functional diversities of tumor-infiltrating immune cells. In this review, we outline the recent progress in cancer immunotherapy, particularly by focusing on landmark studies and the recent single-cell characterization of tumor-associated immune cells, and we summarize the phenotypic diversities of intratumoral immune cells and their connections with cancer immunotherapy. We believe such a review could strengthen our understanding of the progress in cancer immunotherapy, facilitate the elucidation of immune cell modulation in tumor progression, and thus guide the development of novel immunotherapies for cancer treatment.

autophagic responses. Here, we critically discuss current methods of assessing autophagy and the information they can, or cannot, provide. Our ultimate goal is to encourage intellectual and technical innovation in the field.

With the goal of providing a comprehensive, high-quality resource for both plant transcription factors (TFs) and their regulatory interactions with target genes, we upgraded plant TF database PlantTFDB to version 4.0 (http://planttfdb.cbi.pku.edu.cn/). In the new version, we identified 320 370 TFs from 165 species, presenting a more comprehensive genomic TF repertoires of green plants. Besides updating the pre-existing abundant functional and evolutionary annotation for identified TFs, we generated three new types of annotation which provide more directly clues to investigate functional mechanisms underlying: (i) a set of high-quality, non-redundant TF binding motifs derived from experiments; (ii) multiple types of regulatory elements identified from high-throughput sequencing data; (iii) regulatory interactions curated from literature and inferred by combining TF binding motifs and regulatory elements. In addition, we upgraded previous TF prediction server, and set up four novel tools for regulation prediction and functional enrichment analyses. Finally, we set up a novel companion portal PlantRegMap (http://plantregmap.cbi.pku.edu.cn) for users to access the regulation resource and analysis tools conveniently.

Dysregulation of eating behavior can lead to obesity, which affects 10% of the adult population worldwide and accounts for nearly 3 million deaths every year. Despite this burden on society, we currently lack effective pharmacological treatment options to regulate appetite. We used Drosophila melanogaster larvae to develop a high-throughput whole organism screen for drugs that modulate food intake. In a screen of 3630 small molecules, we identified the serotonin (5-hydroxytryptamine or 5-HT) receptor antagonist metitepine as a potent anorectic drug. Using cell-based assays we show that metitepine is an antagonist of all five Drosophila 5-HT receptors. We screened fly mutants for each of these receptors and found that serotonin receptor 5-HT2A is the sole molecular target for feeding inhibition by metitepine. These results highlight the conservation of molecular mechanisms controlling appetite and provide a method for unbiased whole-organism drug screens to identify novel drugs and molecular pathways modulating food intake.

Exosomes are 40-100 nm nano-sized vesicles that are released from many cell types into the extracellular space. Such vesicles are widely distributed in various body fluids. Recently, mRNAs and microRNAs (miRNAs) have been identified in exosomes, which can be taken up by neighboring or distant cells and subsequently modulate recipient cells. This suggests an active sorting mechanism of exosomal miRNAs, since the miRNA profiles of exosomes may differ from those of the parent cells. Exosomal miRNAs play an important role in disease progression, and can stimulate angiogenesis and facilitate metastasis in cancers. In this review, we will introduce the origin and the trafficking of exosomes between cells, display current research on the sorting mechanism of exosomal miRNAs, and briefly describe how exosomes and their miRNAs function in recipient cells. Finally, we will discuss the potential applications of these miRNA-containing vesicles in clinical settings.

Abstract The SARS-CoV-2 B.1.1.529 (Omicron) variant contains 15 mutations of the receptor-binding domain (RBD). How Omicron evades RBD-targeted neutralizing antibodies requires immediate investigation. Here we use high-throughput yeast display screening 1,2 to determine the profiles of RBD escaping mutations for 247 human anti-RBD neutralizing antibodies and show that the neutralizing antibodies can be classified by unsupervised clustering into six epitope groups (A–F)—a grouping that is highly concordant with knowledge-based structural classifications 3–5 . Various single mutations of Omicron can impair neutralizing antibodies of different epitope groups. Specifically, neutralizing antibodies in groups A–D, the epitopes of which overlap with the ACE2-binding motif, are largely escaped by K417N, G446S, E484A and Q493R. Antibodies in group E (for example, S309) 6 and group F (for example, CR3022) 7 , which often exhibit broad sarbecovirus neutralizing activity, are less affected by Omicron, but a subset of neutralizing antibodies are still escaped by G339D, N440K and S371L. Furthermore, Omicron pseudovirus neutralization showed that neutralizing antibodies that sustained single mutations could also be escaped, owing to multiple synergetic mutations on their epitopes. In total, over 85% of the tested neutralizing antibodies were escaped by Omicron. With regard to neutralizing-antibody-based drugs, the neutralization potency of LY-CoV016, LY-CoV555, REGN10933, REGN10987, AZD1061, AZD8895 and BRII-196 was greatly undermined by Omicron, whereas VIR-7831 and DXP-604 still functioned at a reduced efficacy. Together, our data suggest that infection with Omicron would result in considerable humoral immune evasion, and that neutralizing antibodies targeting the sarbecovirus conserved region will remain most effective. Our results inform the development of antibody-based drugs and vaccines against Omicron and future variants.

Higher-order chromatin structure is emerging as an important regulator of gene expression. Although dynamic chromatin structures have been identified in the genome, the full scope of chromatin dynamics during mammalian development and lineage specification remains to be determined. By mapping genome-wide chromatin interactions in human embryonic stem (ES) cells and four human ES-cell-derived lineages, we uncover extensive chromatin reorganization during lineage specification. We observe that although self-associating chromatin domains are stable during differentiation, chromatin interactions both within and between domains change in a striking manner, altering 36% of active and inactive chromosomal compartments throughout the genome. By integrating chromatin interaction maps with haplotype-resolved epigenome and transcriptome data sets, we find widespread allelic bias in gene expression correlated with allele-biased chromatin states of linked promoters and distal enhancers. Our results therefore provide a global view of chromatin dynamics and a resource for studying long-range control of gene expression in distinct human cell lineages. An analysis of genome-wide chromatin interactions during human embryonic stem cell differentiation reveals changes in chromatic organization and simultaneously identifies allele-resolved chromatin structure and differences in gene expression during differentiation. Higher-order chromatin structures are among the factors influencing gene expression, although how these structures evolve during differentiation and lineage specification in mammalian systems is still unclear. Bing Ren and colleagues have mapped the differences in genome-wide chromatin interactions between human embryonic stem cells and their differentiated progeny. They delineate biases in allelic gene expression that correlate with allele-biased chromatin interactions between distal enhancers and proximal promoters.

A novel coronavirus [severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2)] outbreak has caused a global coronavirus disease 2019 (COVID-19) pandemic, resulting in tens of thousands of infections and thousands of deaths worldwide. The RNA-dependent RNA polymerase [(RdRp), also named nsp12] is the central component of coronaviral replication and transcription machinery, and it appears to be a primary target for the antiviral drug remdesivir. We report the cryo-electron microscopy structure of COVID-19 virus full-length nsp12 in complex with cofactors nsp7 and nsp8 at 2.9-angstrom resolution. In addition to the conserved architecture of the polymerase core of the viral polymerase family, nsp12 possesses a newly identified β-hairpin domain at its N terminus. A comparative analysis model shows how remdesivir binds to this polymerase. The structure provides a basis for the design of new antiviral therapeutics that target viral RdRp.

Developing therapeutics against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) could be guided by the distribution of epitopes, not only on the receptor binding domain (RBD) of the Spike (S) protein but also across the full Spike (S) protein. We isolated and characterized monoclonal antibodies (mAbs) from 10 convalescent COVID-19 patients. Three mAbs showed neutralizing activities against authentic SARS-CoV-2. One mAb, named 4A8, exhibits high neutralization potency against both authentic and pseudotyped SARS-CoV-2 but does not bind the RBD. We defined the epitope of 4A8 as the N-terminal domain (NTD) of the S protein by determining with cryo-eletron microscopy its structure in complex with the S protein to an overall resolution of 3.1 angstroms and local resolution of 3.3 angstroms for the 4A8-NTD interface. This points to the NTD as a promising target for therapeutic mAbs against COVID-19.

Macroautophagy/autophagy is an essential, conserved self-eating process that cells perform to allow degradation of intracellular components, including soluble proteins, aggregated proteins, organelles, macromolecular complexes, and foreign bodies. The process requires formation of a double-membrane structure containing the sequestered cytoplasmic material, the autophagosome, that ultimately fuses with the lysosome. This review will define this process and the cellular pathways required, from the formation of the double membrane to the fusion with lysosomes in molecular terms, and in particular highlight the recent progress in our understanding of this complex process.

Abstract Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron sublineages BA.2.12.1, BA.4 and BA.5 exhibit higher transmissibility than the BA.2 lineage 1 . The receptor binding and immune-evasion capability of these recently emerged variants require immediate investigation. Here, coupled with structural comparisons of the spike proteins, we show that BA.2.12.1, BA.4 and BA.5 (BA.4 and BA.5 are hereafter referred collectively to as BA.4/BA.5) exhibit similar binding affinities to BA.2 for the angiotensin-converting enzyme 2 (ACE2) receptor. Of note, BA.2.12.1 and BA.4/BA.5 display increased evasion of neutralizing antibodies compared with BA.2 against plasma from triple-vaccinated individuals or from individuals who developed a BA.1 infection after vaccination. To delineate the underlying antibody-evasion mechanism, we determined the escape mutation profiles 2 , epitope distribution 3 and Omicron-neutralization efficiency of 1,640 neutralizing antibodies directed against the receptor-binding domain of the viral spike protein, including 614 antibodies isolated from people who had recovered from BA.1 infection. BA.1 infection after vaccination predominantly recalls humoral immune memory directed against ancestral (hereafter referred to as wild-type (WT)) SARS-CoV-2 spike protein. The resulting elicited antibodies could neutralize both WT SARS-CoV-2 and BA.1 and are enriched on epitopes on spike that do not bind ACE2. However, most of these cross-reactive neutralizing antibodies are evaded by spike mutants L452Q, L452R and F486V. BA.1 infection can also induce new clones of BA.1-specific antibodies that potently neutralize BA.1. Nevertheless, these neutralizing antibodies are largely evaded by BA.2 and BA.4/BA.5 owing to D405N and F486V mutations, and react weakly to pre-Omicron variants, exhibiting narrow neutralization breadths. The therapeutic neutralizing antibodies bebtelovimab 4 and cilgavimab 5 can effectively neutralize BA.2.12.1 and BA.4/BA.5, whereas the S371F, D405N and R408S mutations undermine most broadly sarbecovirus-neutralizing antibodies. Together, our results indicate that Omicron may evolve mutations to evade the humoral immunity elicited by BA.1 infection, suggesting that BA.1-derived vaccine boosters may not achieve broad-spectrum protection against new Omicron variants.

An atlas of cancer-associated T cells The tumor microenvironment contains many different kinds of immune cells, the composition, function, and roles of which are unclear. Using single-cell RNA sequencing of T cells in 21 cancer types from more than 300 patients, Zheng et al . identified differences in transcript composition that could be used to catalog different T cell types (see the Perspective by van der Leun and Schumacher). These annotations identified the different roles of specific types of CD4 + and CD8 + T cells among the different tumor types. Some of these clusters revealed evidence for two developmental paths for T cells, one of which shows a trajectory toward the “exhausted” T cell state, knowledge of which may be useful in developing future cancer immunotherapies. —LMZ

Pluripotent stem cells can be induced from somatic cells, providing an unlimited cell resource, with potential for studying disease and use in regenerative medicine. However, genetic manipulation and technically challenging strategies such as nuclear transfer used in reprogramming limit their clinical applications. Here, we show that pluripotent stem cells can be generated from mouse somatic cells at a frequency up to 0.2% using a combination of seven small-molecule compounds. The chemically induced pluripotent stem cells resemble embryonic stem cells in terms of their gene expression profiles, epigenetic status, and potential for differentiation and germline transmission. By using small molecules, exogenous "master genes" are dispensable for cell fate reprogramming. This chemical reprogramming strategy has potential use in generating functional desirable cell types for clinical applications.

The pandemic of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has become a global crisis. Replication of SARS-CoV-2 requires the viral RNA-dependent RNA polymerase (RdRp) enzyme, a target of the antiviral drug remdesivir. Here we report the cryo-electron microscopy structure of the SARS-CoV-2 RdRp, both in the apo form at 2.8-angstrom resolution and in complex with a 50-base template-primer RNA and remdesivir at 2.5-angstrom resolution. The complex structure reveals that the partial double-stranded RNA template is inserted into the central channel of the RdRp, where remdesivir is covalently incorporated into the primer strand at the first replicated base pair, and terminates chain elongation. Our structures provide insights into the mechanism of viral RNA replication and a rational template for drug design to combat the viral infection.