Xuya Wang

“Perfect” designer chromosome V and behavior of a ring derivative

INTRODUCTION The Saccharomyces cerevisiae 2.0 project (Sc2.0) aims to modify the yeast genome with a series of densely spaced designer changes. Both a synthetic yeast chromosome arm (synIXR) and the entirely synthetic chromosome (synIII) function with high fitness in yeast. For designer genome synthesis projects, precise engineering of the physical sequence to match the specified design is important for the systematic evaluation of underlying design principles. Yeast can maintain nuclear chromosomes as rings, occurring by chance at repeated sequences, although the cyclized format is unfavorable in meiosis given the possibility of dicentric chromosome formation from meiotic recombination. Here, we describe the de novo synthesis of synthetic yeast chromosome V (synV) in the “Build-A-Genome China” course, perfectly matching the designer sequence and bearing loxPsym sites, distinguishable watermarks, and all the other features of the synthetic genome. We generated a ring synV derivative with user-specified cyclization coordinates and characterized its performance in mitosis and meiosis. RATIONALE Systematic evaluation of underlying Sc2.0 design principles requires that the final assembled synthetic genome perfectly match the designed sequence. Given the size of yeast chromosomes, synthetic chromosome construction is performed iteratively, and new mutations and unpredictable events may occur during synthesis; even a very small number of unintentional nucleotide changes across the genome could have substantial effects on phenotype. Therefore, precisely matching the physical sequence to the designed sequence is crucial for verification of the design principles in genome synthesis. Ring chromosomes can extend those design principles to provide a model for genomic rearrangement, ring chromosome evolution, and human ring chromosome disorders. RESULTS We chemically synthesized, assembled, and incorporated designer chromosome synV (536,024 base pairs) of S. cerevisiae according to Sc2.0 principles, based on the complete nucleotide sequence of native yeast chromosome V (576,874 base pairs). This work was performed as part of the “Build-A-Genome China” course in Tianjin University. We corrected all mutations found—including duplications, substitutions, and indels—in the initial synV strain by using integrative cotransformation of the precise desired changes and by means of a clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)–based method. Altogether, 3331 corrected base pairs were required to match to the designed sequence. We generated a strain that exactly matches all designer sequence changes that displays high fitness under a variety of culture conditions. All corrections were verified with whole-genome sequencing; RNA sequencing revealed only minor changes in gene expression—most notably, decreases in expression of genes relocated near synthetic telomeres as a result of design. We constructed a functional circular synV (ring_synV) derivative in yeast by precisely joining both chromosome ends (telomeres) at specified coordinates. The ring chromosome showed restoration of subtelomeric gene expression levels. The ring_synV strain exhibited fitness comparable with that of the linear synV strain, revealed no change in sporulation frequency, but notably reduced spore viability. In meiosis, heterozygous or homozygous diploid ring_wtV and ring_synV chromosomes behaved similarly, exhibiting substantially higher frequency of the formation of zero-spore tetrads, a type that was not seen in the rod chromosome diploids. Rod synV chromosomes went through meiosis with high spore viability, despite no effort having been made to preserve meiotic competency in the design of synV. CONCLUSION The perfect designer-matched synthetic chromosome V provides strategies to edit sequence variants and correct unpredictable events, such as off-target integration of extra copies of synthetic DNA elsewhere in the genome. We also constructed a ring synthetic chromosome derivative and evaluated its fitness and stability in yeast. Both synV and synVI can be circularized and can power yeast cell growth without affecting fitness when gene content is maintained. These fitness and stability phenotypes of the ring synthetic chromosome in yeast provide a model system with which to probe the mechanism of human ring chromosome disorders. Synthesis, cyclization, and characterization of synV. (A) Synthetic chromosome V (synV, 536,024 base pairs) was designed in silico from native chromosome V (wtV, 576,874 base pairs), with extensive genotype modification designed to be phenotypically neutral. (B) CRISPR/Cas9 strategy for multiplex repair

Synthesis, debugging, and effects of synthetic chromosome consolidation: synVI and beyond

INTRODUCTION Total synthesis of designer chromosomes and genomes is a new paradigm for the study of genetics and biological systems. The Sc2.0 project is building a designer yeast genome from scratch to test and extend the limits of our biological knowledge. Here we describe the design, rapid assembly, and characterization of synthetic chromosome VI (synVI). Further, we investigate the phenotypic, transcriptomic, and proteomic consequences associated with consolidation of three synthetic chromosomes–synVI, synIII, and synIXR—into a single poly-synthetic strain. RATIONALE A host of Sc2.0 chromosomes, including synVI, have now been constructed in discrete strains. With debugging steps, where the number of bugs scales with chromosome length, all individual synthetic chromosomes have been shown to power yeast cells to near wild-type (WT) fitness. Testing the effects of Sc2.0 chromosome consolidation to uncover possible synthetic genetic interactions and/or perturbations of native cellular networks as the number of designer changes increases is the next major step for the Sc2.0 project. RESULTS SynVI was rapidly assembled using nine sequential steps of SwAP-In (switching auxotrophies progressively by integration), yielding a ~240-kb synthetic chromosome designed to Sc2.0 specifications. We observed partial silencing of the left- and rightmost genes on synVI, each newly positioned subtelomerically relative to their locations on native VI. This result suggests that consensus core X elements of Sc2.0 universal telomere caps are insufficient to fully buffer telomere position effects. The synVI strain displayed a growth defect characterized by an increased frequency of glycerol-negative colonies. The defect mapped to a synVI design feature in the essential PRE4 gene (YFR050C), encoding the β7 subunit of the 20S proteasome. Recoding 10 codons near the 3′ end of the PRE4 open reading frame (ORF) caused a ~twofold reduction in Pre4 protein level without affecting RNA abundance. Reverting the codons to the WT sequence corrected both the Pre4 protein level and the phenotype. We hypothesize that the formation of a stem loop involving recoded codons underlies reduced Pre4 protein level. Sc2.0 chromosomes (synI to synXVI) are constructed individually in discrete strains and consolidated into poly-synthetic (poly-syn) strains by “endoreduplication intercross.” Consolidation of synVI with synthetic chromosomes III (synIII) and IXR (synIXR) yields a triple-synthetic (triple-syn) strain that is ~6% synthetic overall—with almost 70 kb deleted, including 20 tRNAs, and more than 12 kb recoded. Genome sequencing of double-synthetic (synIII synVI, synIII synIXR, synVI synIXR) and triple-syn (synIII synVI synIXR) cells indicates that suppressor mutations are not required to enable coexistence of Sc2.0 chromosomes. Phenotypic analysis revealed a slightly slower growth rate for the triple-syn strain only; the combined effect of tRNA deletions on different chromosomes might underlie this result. Transcriptome and proteome analyses indicate that cellular networks are largely unperturbed by the existence of multiple synthetic chromosomes in a single cell. However, a second bug on synVI was discovered through proteomic analysis and is associated with alteration of the HIS2 transcription start as a consequence of tRNA deletion and loxPsym site insertion. Despite extensive genetic alterations across 6% of the genome, no major global changes were detected in the poly-syn strain “omics” analyses. CONCLUSION Analyses of phenotypes, transcriptomics, and proteomics of synVI and poly-syn strains reveal, in general, WT cell properties and the existence of rare bugs resulting from genome editing. Deletion of subtelomeres can lead to gene silencing, recoding deep within an ORF can yield a translational defect, and deletion of elements such as tRNA genes can lead to a complex transcriptional output. These results underscore the complementarity of transcriptomics and proteomics to identify bugs, the consequences of designer changes in Sc2.0 chromosomes. The consolidation of Sc2.0 designer chromosomes into a single strain appears to be exceptionally well tolerated by yeast. A predictable exception to this is the deletion of tRNAs, which will be restored on a separate neochromosome to avoid synthetic lethal genetic interactions between deleted tRNA genes as additional synthetic chromosomes are introduced. Debugging synVI and characterization of poly-synthetic yeast cells. (A) The second Sc2.0 chromosome to be constructed, synVI, encodes a “bug” that causes a variable colony size, dubbed a “glycerol-negative growth-suppression defect.” (B) Synonymous changes in the essential PRE4 ORF lead to a reduced protein level, which underlies the growth defect

Bug mapping and fitness testing of chemically synthesized chromosome X

INTRODUCTION Design and construction of an extensively modified yeast genome is a direct means to interrogate the integrity, comprehensiveness, and accuracy of the knowledge amassed by the yeast community to date. The international synthetic yeast genome project (Sc2.0) aims to build an entirely designer, synthetic Saccharomyces cerevisiae genome. The synthetic genome is designed to increase genome stability and genetic flexibility while maintaining cell fitness near that of the wild type. A major challenge for a genome synthesis lies in identifying and eliminating fitness-reducing sequence variants referred to as “bugs.” RATIONALE Debugging is imperative for successfully building a fit strain encoding a synthetic genome. However, it is time-consuming and laborious to replace wild-type genes and measure strain fitness systematically. The Sc2.0 PCRTag system, which specifies recoded sequences within open reading frames (ORFs), is designed to distinguish synthetic from wild-type DNA in a simple polymerase chain reaction (PCR) assay. This system provides an opportunity to efficiently map bugs to the related genes by using a pooling strategy and subsequently correct them. Further, as we identify bugs in designer sequences, we will identify gaps in our knowledge and gain a deeper understanding of genome biology, allowing refinement of future design strategies. RESULTS We chemically synthesized yeast chromosome X, synX, designed to be 707,459 base pairs. A high-throughput mapping strategy called pooled PCRTag mapping (PoPM) was developed to identify unexpected bugs during chromosome assembly. With this method, the genotypes of pools of colonies with normal or defective fitness are assessed by PCRTag analysis. The PoPM method exploits the patchwork structure of synthetic and wild-type sequences observed in the majority of putative synthetic DNA integrants or meiotic progeny derived from synthetic/wild-type strain backcross. PCRTag analysis with both synthetic and wild-type specific primers, carried out with genomic DNA extracted from the two pools of clones (normal fitness versus a specific growth defect), can be used to identify regions of synthetic DNA missing from the normal fitness pool and, analogously, sections of wild-type DNA absent from the specific growth-defect pool. In this way, the defect can be efficiently mapped to a very small overlapping region, and subsequent systematic analysis of designed changes in that region can be used to identify the bug. Several bugs were identified and corrected, including a growth defect mapping to a specific synonymously recoded PCRTag sequence in the essential FIP1 ORF and the effect of introducing a loxPsym site that unexpectedly altered the the promoter function of a nearby gene, ATP2. In addition, meiotic crossover was employed to repair the massive duplications and rearrangements in the synthetic chromosome. The debugged synX strain exhibited high fitness under a variety of conditions tested and in competitive growth with the wild-type strain. CONCLUSION Synthetic yeast chromosome X was chemically synthesized from scratch, a rigorous, incremental step toward complete synthesis of the whole yeast genome. Thousands of designer modifications in synX revealed extensive flexibility of the yeast genome. We developed an efficient mapping method, PoPM, to identify bugs during genome synthesis, generalizable to any watermarked synthetic chromosome, and several details of yeast biology were uncovered by debugging. Considering the numerous gene-associated PCRTags available in the synthetic chromosomes, PoPM may represent a powerful tool to map interesting phenotypes of mutated synthetic strains or even mutated wild-type strains to the relevant genes. It may also be useful to study yeast genetic interactions when an unexpected phenotype is generated by alterations in two or more genes, substantially expanding understanding of yeast genomic and cellular functions. The PoPM method is also likely to be useful for mapping phenotype(s) resulting from the genome SCRaMbLE system. Characterization of synX and debugging by pooled PCRTag mapping. (Top) Design overview of synthetic chromosome X. (Bottom) Flow diagram of pooled PCRTag mapping (PoPM). Debugging a genome sequence is imperative for successfully building a synthetic genome. As part of the effort to build a designer eukaryotic genome, yeast synthetic chromosome X (synX), designed as 707,459 base pairs, was synthesized chemically. SynX exhibited good fitness under a wide variety of conditions. A highly efficient mapping strategy called pooled PCRTag mapping (PoPM), which can be generalized to any watermarked synthetic chromosome, was developed to identify genetic alterations that affect cell fitness (“bugs”). A series of bugs were corrected that included a large region bearing complex amplifications, a growth defect mapping to a recoded sequence in FIP1, and a loxPsym site affecting promoter function of ATP2. PoPM is a powerful tool for synthetic yeast genome debugging and an efficient strategy for phenotype-genotype mapping.

An Analysis of the Sensitivity of Proteogenomic Mapping of Somatic Mutations and Novel Splicing Events in Cancer

Improvements in mass spectrometry (MS)-based peptide sequencing provide a new opportunity to determine whether polymorphisms, mutations, and splice variants identified in cancer cells are translated. Herein, we apply a proteogenomic data integration tool (QUILTS) to illustrate protein variant discovery using whole genome, whole transcriptome, and global proteome datasets generated from a pair of luminal and basal-like breast-cancer-patient-derived xenografts (PDX). The sensitivity of proteogenomic analysis for singe nucleotide variant (SNV) expression and novel splice junction (NSJ) detection was probed using multiple MS/MS sample process replicates defined here as an independent tandem MS experiment using identical sample material. Despite analysis of over 30 sample process replicates, only about 10% of SNVs (somatic and germline) detected by both DNA and RNA sequencing were observed as peptides. An even smaller proportion of peptides corresponding to NSJ observed by RNA sequencing were detected (<0.1%). Peptides mapping to DNA-detected SNVs without a detectable mRNA transcript were also observed, suggesting that transcriptome coverage was incomplete (∼80%). In contrast to germline variants, somatic variants were less likely to be detected at the peptide level in the basal-like tumor than in the luminal tumor, raising the possibility of differential translation or protein degradation effects. In conclusion, this large-scale proteogenomic integration allowed us to determine the degree to which mutations are translated and identify gaps in sequence coverage, thereby benchmarking current technology and progress toward whole cancer proteome and transcriptome analysis.

Proteogenomic integration reveals therapeutic targets in breast cancer xenografts.

Recent advances in mass spectrometry (MS) have enabled extensive analysis of cancer proteomes. Here, we employed quantitative proteomics to profile protein expression across 24 breast cancer patient-derived xenograft (PDX) models. Integrated proteogenomic analysis shows positive correlation between expression measurements from transcriptomic and proteomic analyses; further, gene expression-based intrinsic subtypes are largely re-capitulated using non-stromal protein markers. Proteogenomic analysis also validates a number of predicted genomic targets in multiple receptor tyrosine kinases. However, several protein/phosphoprotein events such as overexpression of AKT proteins and ARAF, BRAF, HSP90AB1 phosphosites are not readily explainable by genomic analysis, suggesting that druggable translational and/or post-translational regulatory events may be uniquely diagnosed by MS. Drug treatment experiments targeting HER2 and components of the PI3K pathway supported proteogenomic response predictions in seven xenograft models. Our study demonstrates that MS-based proteomics can identify therapeutic targets and highlights the potential of PDX drug response evaluation to annotate MS-based pathway activities.

Breast tumors educate the proteome of stromal tissue in an individualized but coordinated manner

Proteomic analysis of the tumor-associated stroma reveals extensive and coordinated regulation by breast cancers. Profiling the tumor stroma proteome Communication between a tumor and cells in the surrounding stroma contributes to tumor growth, progression, and drug resistance. Thus, targeting this communication, in the primary tumor and especially in metastatic niches, may be an effective way to treat cancer. Wang et al. grew patient breast tumors subcutaneously in mice and obtained species-distinguished proteomic profiles of the tumors (human) and tumor-associated stroma (mouse). The authors found that all breast tumors consistently altered clustered subsets of the stromal proteome, particularly proteins involved in immune signaling, but that these varied in a subtype- and stage-specific manner. These findings may have future implications for treatment stratification and provide a platform from which to understand this experimental model and tumor-stroma interactions on a large-scale protein level. Cancer forms specialized microenvironmental niches that promote local invasion and colonization. Engrafted patient-derived xenografts (PDXs) locally invade and colonize naïve stroma in mice while enabling unambiguous molecular discrimination of human proteins in the tumor from mouse proteins in the microenvironment. To characterize how patient breast tumors form a niche and educate naïve stroma, subcutaneous breast cancer PDXs were globally profiled by species-specific quantitative proteomics. Regulation of PDX stromal proteins by breast tumors was extensive, with 35% of the stromal proteome altered by tumors consistently across different animals and passages. Differentially regulated proteins in the stroma clustered into six signatures, which included both known and previously unappreciated contributors to tumor invasion and colonization. Stromal proteomes were coordinately regulated; however, the sets of proteins altered by each tumor were highly distinct. Integrated analysis of tumor and stromal proteins, a comparison made possible in these xenograft models, indicated that the known hallmarks of cancer contribute pleiotropically to establishing and maintaining the microenvironmental niche of the tumor. Education of the stroma by the tumor is therefore an intrinsic property of breast tumors that is highly individualized, yet proceeds by consistent, nonrandom, and defined tumor-promoting molecular alterations.

Mass spectrometry-based proteomics reveals potential roles of NEK9 and MAP2K4 in resistance to PI3K inhibitors in triple negative breast cancers

Activation of PI3K signaling is frequently observed in triple-negative breast cancer (TNBC), yet PI3K inhibitors have shown limited clinical activity. To investigate intrinsic and adaptive mechanisms of resistance, we analyzed a panel of patient-derived xenograft models of TNBC with varying responsiveness to buparlisib, a pan-PI3K inhibitor. In a subset of patient-derived xenografts, resistance was associated with incomplete inhibition of PI3K signaling and upregulated MAPK/MEK signaling in response to buparlisib. Outlier phosphoproteome and kinome analyses identified novel candidates functionally important to buparlisib resistance, including NEK9 and MAP2K4. Knockdown of NEK9 or MAP2K4 reduced both baseline and feedback MAPK/MEK signaling and showed synthetic lethality with buparlisib in vitro A complex in/del frameshift in PIK3CA decreased sensitivity to buparlisib via NEK9/MAP2K4-dependent mechanisms. In summary, our study supports a role for NEK9 and MAP2K4 in mediating buparlisib resistance and demonstrates the value of unbiased omic analyses in uncovering resistance mechanisms to targeted therapy.Significance: Integrative phosphoproteogenomic analysis is used to determine intrinsic resistance mechanisms of triple-negative breast tumors to PI3K inhibition. Cancer Res; 78(10); 2732-46. ©2018 AACR.

Transcription factor profiling reveals molecular choreography and key regulators of human retrotransposon expression

Significance Retrotransposons replicate through RNA intermediates that are reverse transcribed and inserted at new genomic locations. LINE-1 (L1) elements constitute ∼17% of the human genome, making them the most successful retrotransposons in the human genome by mass. The activity of L1s was shown first in the germline or during early embryogenesis. More recent studies demonstrate a wider prevalence of L1 expression in somatic cells including neurons, aging cells, and different types of cancer. In this study, we developed the MapRRCon pipeline and performed a comprehensive computational analysis of L1 transcriptional regulators using ENCODE ChIP-seq datasets. We revealed the binding of various transcription factors, including Myc and CTCF, to the 5′ UTR promoter of the youngest human L1 family (L1HS) and their potential functional impact on L1HS expression. Transposable elements (TEs) represent a substantial fraction of many eukaryotic genomes, and transcriptional regulation of these factors is important to determine TE activities in human cells. However, due to the repetitive nature of TEs, identifying transcription factor (TF)-binding sites from ChIP-sequencing (ChIP-seq) datasets is challenging. Current algorithms are focused on subtle differences between TE copies and thus bias the analysis to relatively old and inactive TEs. Here we describe an approach termed “MapRRCon” (mapping repeat reads to a consensus) which allows us to identify proteins binding to TE DNA sequences by mapping ChIP-seq reads to the TE consensus sequence after whole-genome alignment. Although this method does not assign binding sites to individual insertions in the genome, it provides a landscape of interacting TFs by capturing factors that bind to TEs under various conditions. We applied this method to screen TFs’ interaction with L1 in human cells/tissues using ENCODE ChIP-seq datasets and identified 178 of the 512 TFs tested as bound to L1 in at least one biological condition with most of them (138) localized to the promoter. Among these L1-binding factors, we focused on Myc and CTCF, as they play important roles in cancer progression and 3D chromatin structure formation. Furthermore, we explored the transcriptomes of The Cancer Genome Atlas breast and ovarian tumor samples in which a consistent anti-/correlation between L1 and Myc/CTCF expression was observed, suggesting that these two factors may play roles in regulating L1 transcription during the development of such tumors.

The largest SWI/SNF polyglutamine domain is a pH sensor

Polyglutamines are known to form aggregates in pathogenic contexts, such as in Huntington’s disease, however little is known about their role in normal biological processes. We found that a polyglutamine domain in the SNF5 subunit of the yeast SWI/SNF complex, histidines within this sequence, and transient intracellular acidification are required for efficient transcriptional regulation during carbon starvation. We hypothesized that a pH-driven oligomerization of the SNF5 polyglutamine region is required for transcriptional reprogramming. In support of this idea, we found that a synthetic spidroin domain from spider silk, which is soluble at pH 7 but oligomerizes at pH ~ 6.3, could partially complement the function of the SNF5 polyglutamine. These results suggest that the SNF5 polyglutamine domain acts as a pH-driven transcriptional regulator.

Corrigendum: Proteogenomic integration reveals therapeutic targets in breast cancer xenografts.

Nature Communications 8: Article number: 14864 (2017)); Published: 28 March 2017; Updated: 25 April 2017 The original version of this Article contained a typographical error in the spelling of the author Beifang Niu, which was incorrectly given as Beifung Niu. This has now been corrected in both thePDF and HTML versions of the Article.