Multiplex generation and single cell analysis of structural variants in a mammalian genome.

The functional consequences of structural variants (SVs) in mammalian genomes are challenging to study. This is due to several factors, including: 1) their numerical paucity relative to other forms of standing genetic variation such as single nucleotide variants (SNVs) and short insertions or deletions (indels); 2) the fact that a single SV can involve and potentially impact the function of more than one gene and/or cis regulatory element; and 3) the relative immaturity of methods to generate and map SVs, either randomly or in targeted fashion, in in vitro or in vivo model systems. Towards addressing these challenges, we developed Genome-Shuffle-seq, a straightforward method that enables the multiplex generation and mapping of several major forms of SVs (deletions, inversions, translocations) throughout a mammalian genome. Genome-Shuffle-seq is based on the integration of "shuffle cassettes" to the genome, wherein each shuffle cassette contains components that facilitate its site-specific recombination (SSR) with other integrated shuffle cassettes (via Cre-loxP), its mapping to a specific genomic location (via T7-mediated in vitro transcription or IVT), and its identification in single-cell RNA-seq (scRNA-seq) data (via T7-mediated in situ transcription or IST). In this proof-of-concept, we apply Genome-Shuffle-seq to induce and map thousands of genomic SVs in mouse embryonic stem cells (mESCs) in a single experiment. Induced SVs are rapidly depleted from the cellular population over time, possibly due to Cre-mediated toxicity and/or negative selection on the rearrangements themselves. Leveraging T7 IST of barcodes whose positions are already mapped, we further demonstrate that we can efficiently genotype which SVs are present in association with each of many single cell transcriptomes in scRNA-seq data. Finally, preliminary evidence suggests our method may be a powerful means of generating extrachromosomal circular DNAs (ecDNAs). Looking forward, we anticipate that Genome-Shuffle-seq may be broadly useful for the systematic exploration of the functional consequences of SVs on gene expression, the chromatin landscape, and 3D nuclear architecture. We further anticipate potential uses for in vitro modeling of ecDNAs, as well as in paving the path to a minimal mammalian genome.

Resurrecting essential amino acid biosynthesis in mammalian cells.

Major genomic deletions in independent eukaryotic lineages have led to repeated ancestral loss of biosynthesis pathways for nine of the twenty canonical amino acids. While the evolutionary forces driving these polyphyletic deletion events are not well understood, the consequence is that extant metazoans are unable to produce nine essential amino acids (EAAs). Previous studies have highlighted that EAA biosynthesis tends to be more energetically costly, raising the possibility that these pathways were lost from organisms with access to abundant EAAs. It is unclear whether present-day metazoans can reaccept these pathways to resurrect biosynthetic capabilities that were lost long ago or whether evolution has rendered EAA pathways incompatible with metazoan metabolism. Here, we report progress on a large-scale synthetic genomics effort to reestablish EAA biosynthetic functionality in mammalian cells. We designed codon-optimized biosynthesis pathways based on genes mined from Escherichia coli. These pathways were de novo synthesized in 3 kilobase chunks, assembled in yeasto and genomically integrated into a Chinese hamster ovary (CHO) cell line. One synthetic pathway produced valine at a sufficient level for cell viability and proliferation. 13C-tracing verified de novo biosynthesis of valine and further revealed build-up of pathway intermediate 2,3-dihydroxy-3-isovalerate. Increasing the dosage of downstream ilvD boosted pathway performance and allowed for long-term propagation of second-generation cells in valine-free medium at 3.2 days per doubling. This work demonstrates that mammalian metabolism is amenable to restoration of ancient core pathways, paving a path for genome-scale efforts to synthetically restore metabolic functions to the metazoan lineage.

Synthetic regulatory reconstitution reveals principles of mammalian Hox cluster regulation.

Precise Hox gene expression is crucial for embryonic patterning. Intra-Hox transcription factor binding and distal enhancer elements have emerged as the major regulatory modules controlling Hox gene expression. However, quantifying their relative contributions has remained elusive. Here, we introduce "synthetic regulatory reconstitution," a conceptual framework for studying gene regulation, and apply it to the HoxA cluster. We synthesized and delivered variant rat HoxA clusters (130 to 170 kilobases) to an ectopic location in the mouse genome. We found that a minimal HoxA cluster recapitulated correct patterns of chromatin remodeling and transcription in response to patterning signals, whereas the addition of distal enhancers was needed for full transcriptional output. Synthetic regulatory reconstitution could provide a generalizable strategy for deciphering the regulatory logic of gene expression in complex genomes.

A versatile platform for locus-scale genome rewriting and verification.

Routine rewriting of loci associated with human traits and diseases would facilitate their functional analysis. However, existing DNA integration approaches are limited in terms of scalability and portability across genomic loci and cellular contexts. We describe Big-IN, a versatile platform for targeted integration of large DNAs into mammalian cells. CRISPR/Cas9-mediated targeting of a landing pad enables subsequent recombinase-mediated delivery of variant payloads and efficient positive/negative selection for correct clones in mammalian stem cells. We demonstrate integration of constructs up to 143 kb, and an approach for one-step scarless delivery. We developed a staged pipeline combining PCR genotyping and targeted capture sequencing for economical and comprehensive verification of engineered stem cells. Our approach should enable combinatorial interrogation of genomic functional elements and systematic locus-scale analysis of genome function.

De novo assembly and delivery to mouse cells of a 101 kb functional human gene.

Design and large-scale synthesis of DNA has been applied to the functional study of viral and microbial genomes. New and expanded technology development is required to unlock the transformative potential of such bottom-up approaches to the study of larger mammalian genomes. Two major challenges include assembling and delivering long DNA sequences. Here, we describe a workflow for de novo DNA assembly and delivery that enables functional evaluation of mammalian genes on the length scale of 100 kilobase pairs (kb). The DNA assembly step is supported by an integrated robotic workcell. We demonstrate assembly of the 101 kb human HPRT1 gene in yeast from 3 kb building blocks, precision delivery of the resulting construct to mouse embryonic stem cells, and subsequent expression of the human protein from its full-length human gene in mouse cells. This workflow provides a framework for mammalian genome writing. We envision utility in producing designer variants of human genes linked to disease and their delivery and functional analysis in cell culture or animal models.

Ancestral reconstruction reveals mechanisms of ERK regulatory evolution.

Protein kinases are crucial to coordinate cellular decisions and therefore their activities are strictly regulated. Previously we used ancestral reconstruction to determine how CMGC group kinase specificity evolved (Howard et al., 2014). In the present study, we reconstructed ancestral kinases to study the evolution of regulation, from the inferred ancestor of CDKs and MAPKs, to modern ERKs. Kinases switched from high to low autophosphorylation activity at the transition to the inferred ancestor of ERKs 1 and 2. Two synergistic amino acid changes were sufficient to induce this change: shortening of the ß3-αC loop and mutation of the gatekeeper residue. Restoring these two mutations to their inferred ancestral state led to a loss of dependence of modern ERKs 1 and 2 on the upstream activating kinase MEK in human cells. Our results shed light on the evolutionary mechanisms that led to the tight regulation of a kinase that is central in development and disease.

Msn2/4 regulate expression of glycolytic enzymes and control transition from quiescence to growth.

Nutrient availability and stresses impact a cell's decision to enter a growth state or a quiescent state. Acetyl-CoA stimulates cell growth under nutrient-limiting conditions, but how cells generate acetyl-CoA under starvation stress is less understood. Here, we show that general stress response factors, Msn2 and Msn4, function as master transcriptional regulators of yeast glycolysis via directly binding and activating genes encoding glycolytic enzymes. Yeast cells lacking Msn2 and Msn4 exhibit prevalent repression of glycolytic genes and a significant delay of acetyl-CoA accumulation and reentry into growth from quiescence. Thus Msn2/4 exhibit a dual role in activating carbohydrate metabolism genes and stress response genes. These results suggest a possible mechanism by which starvation-induced stress response factors may prime quiescent cells to reenter growth through glycolysis when nutrients are limited.

Functions of Peptidoglycan Recognition Proteins (Pglyrps) at the Ocular Surface: Bacterial Keratitis in Gene-Targeted Mice Deficient in Pglyrp-2, -3 and -4.

PURPOSE: Functions of antimicrobial peptidoglycan recognition proteins (Pglyrp1-4) at the ocular surface are poorly understood. Earlier, we reported an antibacterial role for Pglyrp-1 in Pseudomonas aeruginosa keratitis. Here we investigated functions of three other related genes Pglyrp-2, -3 and -4 in a mouse model of P. aeruginosa keratitis. METHODS: Wild type (WT) and each of the Pglyrp-null genotypes were challenged with P. aeruginosa keratitis. The eyes were scored in a blinded manner 24 and 48h post infection. Viable bacterial counts and inflammatory factors (IL-12, TNF-α, IFN-γ, CCL2, IL-6 and IL-10) were measured in whole eye homogenates using cytometric bead arrays. Expressions of Pglyrp-1-4, mouse beta defensins (mBD)-2,-3, cathelicidin-related antimicrobial peptide (CRAMP) were determined by qRTPCR in total RNA extracts of uninfected and infected eyes of WT and each of the Pglyrp-null mouse types. RESULTS: The Pglyrp-2-/- mice showed reduced disease and lower induction of pro-inflammatory TNF-α (p = 0.02) than WT or the other Pglyrp null mice. Viable bacterial yield was significantly lower in the Pglyrp-2-/- (p = 0.0007) and the Pglyrp-4-/- (p = 0.098) mice. With regards to expression of these antimicrobial genes, Pglyrp-2 expression was induced after infection in WT mice. Pglyrp-3 expression was low before and after infection in WT mice, while Pglyrp-4 expression was slightly elevated after infection in WT, Pglyrp-2 and -3 null mice. Pglyrp-1 expression was slightly elevated after infection in all genotypes without statistical significance. Transcripts for antimicrobial peptides mBD2, mBD3 and CRAMP were elevated in infected Pglyrp-2-/- males without statistical significance. CONCLUSIONS: Efficient resolution of keratitis in the Pglyrp-2-/- mice may be due to a reduced pro-inflammatory microenvironment and synergistic antibacterial activities of defensins, CRAMP and Pglyrp-1. Therefore, in ocular infections the pro-inflammatory functions of Pglyrp-2 must be regulated to benefit the host.