Faculty Bibliography

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'^1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

BACKGROUND:Xenotransplantation of genetically engineered porcine organs has the potential to address the challenge of organ donor shortage. Two cases of porcine-to-human kidney xenotransplantation were performed, yet the physiological effects on the xenografts and the recipients' immune responses remain largely uncharacterized. METHODS:We performed single-cell RNA sequencing (scRNA-seq) and longitudinal RNA-seq analyses of the porcine kidneys to dissect xenotransplantation-associated cellular dynamics and xenograft-recipient interactions. We additionally performed longitudinal scRNA-seq of the peripheral blood mononuclear cells (PBMCs) to detect recipient immune responses across time. FINDINGS/RESULTS:Although no hyperacute rejection signals were detected, scRNA-seq analyses of the xenografts found evidence of endothelial cell and immune response activation, indicating early signs of antibody-mediated rejection. Tracing the cells' species origin, we found human immune cell infiltration in both xenografts. Human transcripts in the longitudinal bulk RNA-seq revealed that human immune cell infiltration and the activation of interferon-gamma-induced chemokine expression occurred by 12 and 48 h post-xenotransplantation, respectively. Concordantly, longitudinal scRNA-seq of PBMCs also revealed two phases of the recipients' immune responses at 12 and 48-53 h. Lastly, we observed global expression signatures of xenotransplantation-associated kidney tissue damage in the xenografts. Surprisingly, we detected a rapid increase of proliferative cells in both xenografts, indicating the activation of the porcine tissue repair program. CONCLUSIONS:Longitudinal and single-cell transcriptomic analyses of porcine kidneys and the recipient's PBMCs revealed time-resolved cellular dynamics of xenograft-recipient interactions during xenotransplantation. These cues can be leveraged for designing gene edits and immunosuppression regimens to optimize xenotransplantation outcomes. FUNDING/BACKGROUND:This work was supported by NIH RM1HG009491 and DP5OD033430.

In a previous study, heart xenografts from 10-gene-edited pigs transplanted into two human decedents did not show evidence of acute-onset cellular- or antibody-mediated rejection. Here, to better understand the detailed molecular landscape following xenotransplantation, we carried out bulk and single-cell transcriptomics, lipidomics, proteomics and metabolomics on blood samples obtained from the transplanted decedents every 6 h, as well as histological and transcriptomic tissue profiling. We observed substantial early immune responses in peripheral blood mononuclear cells and xenograft tissue obtained from decedent 1 (male), associated with downstream T cell and natural killer cell activity. Longitudinal analyses indicated the presence of ischemia reperfusion injury, exacerbated by inadequate immunosuppression of T cells, consistent with previous findings of perioperative cardiac xenograft dysfunction in pig-to-nonhuman primate studies. Moreover, at 42 h after transplantation, substantial alterations in cellular metabolism and liver-damage pathways occurred, correlating with profound organ-wide physiological dysfunction. By contrast, relatively minor changes in RNA, protein, lipid and metabolism profiles were observed in decedent 2 (female) as compared to decedent 1. Overall, these multi-omics analyses delineate distinct responses to cardiac xenotransplantation in the two human decedents and reveal new insights into early molecular and immune responses after xenotransplantation. These findings may aid in the development of targeted therapeutic approaches to limit ischemia reperfusion injury-related phenotypes and improve outcomes.

Core histone genes display a remarkable diversity of cis-regulatory mechanisms despite their protein sequence conservation. However, the dynamics and significance of this regulatory turnover are not well understood. Here we describe the evolutionary history of core histone gene regulation across 400 million years in budding yeasts. We find that canonical mode of core histone regulation - mediated by the trans-regulator Spt10 - is ancient, likely emerging between 320-380 million years ago and is fixed in the majority of extant species. Unexpectedly, we uncovered the emergence of a novel core histone regulatory mode in the Hanseniaspora genus, from its fast-evolving lineage (FEL), which coincided with the loss of one copy of its paralogous core histones genes. We show that the ancestral Spt10 histone regulatory mode was replaced, via cis-regulatory changes in the histone control regions, by a derived Mcm1 histone regulatory mode and that this rewiring event occurred with no changes to the trans-regulator, Mcm1, itself. Finally, we studied the growth dynamics of the cell cycle and histone synthesis in genetically modified Hanseniaspora uvarum. We find that H. uvarum divides rapidly, with most cells completing a cell cycle within 60 minutes. Interestingly, we observed that the regulatory coupling between histone and DNA synthesis was lost in H. uvarum. Our results demonstrate that core histone gene regulation was fixed anciently in budding yeasts, however it has greatly diverged in the Hanseniaspora FEL.

Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) is a promising tool to study genomic rearrangements. However, the potential of SCRaMbLE to study genomic rearrangements is currently hindered, because a strain containing all 16 synthetic chromosomes is not yet available. Here, we construct SparLox83R, a yeast strain containing 83 loxPsym sites distributed across all 16 chromosomes. SCRaMbLE of SparLox83R produces versatile genome-wide genomic rearrangements, including inter-chromosomal events. Moreover, when combined with synthetic chromosomes, SCRaMbLE of hetero-diploids with SparLox83R leads to increased diversity of genomic rearrangements and relatively faster evolution of traits compared to hetero-diploids only with wild-type chromosomes. Analysis of the SCRaMbLEd strain with increased tolerance to nocodazole demonstrates that genomic rearrangements can perturb the transcriptome and 3D genome structure and consequently impact phenotypes. In summary, a genome with sparsely distributed loxPsym sites can serve as a powerful tool for studying the consequence of genomic rearrangements and accelerating strain engineering in Saccharomyces cerevisiae.

DNA targeting Class 2 CRISPR-Cas effector nucleases, including the well-studied Cas9 proteins, evolved protospacer-adjacent motif (PAM) and guide RNA interactions that sequentially license their binding and cleavage activities at protospacer target sites. Both interactions are nucleic acid sequence specific but function constitutively; thus, they provide intrinsic spatial control over DNA targeting activities but naturally lack temporal control. Here we show that engineered Cas9 fusion proteins which bind to nascent RNAs near a protospacer can facilitate spatiotemporal coupling between transcription and DNA targeting at that protospacer: Transcription-associated Cas9 Targeting (TraCT). Engineered TraCT is enabled when suboptimal PAM interactions limit basal activity in vivo and when one or more nascent RNA substrates are still tethered to the actively transcribing target DNA in cis. We further show that this phenomenon can be exploited for selective editing at one of two identical targets in distinct gene loci, or, in diploid allelic loci that are differentially transcribed. Our work demonstrates that temporal control over Cas9's targeting activity at specific DNA sites may be engineered without modifying Cas9's core domains and guide RNA components or their expression levels. More broadly, it establishes RNA binding in cis as a mechanism that can conditionally stimulate CRISPR-Cas DNA targeting in eukaryotes.

Enhancer function is frequently investigated piecemeal using truncated reporter assays or single deletion analysis. Thus it remains unclear to what extent enhancer function at native loci relies on surrounding genomic context. Using the Big-IN technology for targeted integration of large DNAs, we analyzed the regulatory architecture of the murine Igf2/H19 locus, a paradigmatic model of enhancer selectivity. We assembled payloads containing a 157-kb functional Igf2/H19 locus and engineered mutations to genetically direct CTCF occupancy at the imprinting control region (ICR) that switches the target gene of the H19 enhancer cluster. Contrasting the activity of payloads delivered to the endogenous locus or to a safe harbor locus (Hprt) revealed that the Igf2/H19 locus includes additional, previously unknown long-range regulatory elements. Exchanging components of the Igf2/H19 locus with the well-studied Sox2 locus showed that the H19 enhancer cluster functioned poorly out of context, and required its native surroundings to activate Sox2 expression. Conversely, the Sox2 locus control region (LCR) could activate both Igf2 and H19 outside its native context, but its activity was only partially modulated by CTCF occupancy at the ICR. Analysis of regulatory DNA actuation across different cell types revealed that, while the H19 enhancers are tightly coordinated within their native locus, the Sox2 LCR acts more independently. We show that these enhancer clusters typify broader classes of loci genome-wide. Our results show that unexpected dependencies may influence even the most studied functional elements, and our synthetic regulatory genomics approach permits large-scale manipulation of complete loci to investigate the relationship between locus architecture and function.

We describe construction of the synthetic yeast chromosome XI (synXI) and reveal the effects of redesign at non-coding DNA elements. The 660-kb synthetic yeast genome project (Sc2.0) chromosome was assembled from synthesized DNA fragments before CRISPR-based methods were used in a process of bug discovery, redesign, and chromosome repair, including precise compaction of 200 kb of repeat sequence. Repaired defects were related to poor centromere function and mitochondrial health and were associated with modifications to non-coding regions. As part of the Sc2.0 design, loxPsym sequences for Cre-mediated recombination are inserted between most genes. Using the GAP1 locus from chromosome XI, we show that these sites can facilitate induced extrachromosomal circular DNA (eccDNA) formation, allowing direct study of the effects and propagation of these important molecules. Construction and characterization of synXI contributes to our understanding of non-coding DNA elements, provides a useful tool for eccDNA study, and will inform future synthetic genome design.