"Calling cards" for DNA-binding proteins in mammalian cells

Transcriptional precision in photoreceptor development and diseases - Lessons from 25 years of CRX research

The vertebrate retina is made up of six specialized neuronal cell types and one glia that are generated from a common retinal progenitor. The development of these distinct cell types is programmed by transcription factors that regulate the expression of specific genes essential for cell fate specification and differentiation. Because of the complex nature of transcriptional regulation, understanding transcription factor functions in development and disease is challenging. Research on the Cone-rod homeobox transcription factor CRX provides an excellent model to address these challenges. In this review, we reflect on 25 years of mammalian CRX research and discuss recent progress in elucidating the distinct pathogenic mechanisms of four CRX coding variant classes. We highlight how in vitro biochemical studies of CRX protein functions facilitate understanding CRX regulatory principles in animal models. We conclude with a brief discussion of the emerging systems biology approaches that could accelerate precision medicine for CRX-linked diseases and beyond.

Pycallingcards: an integrated environment for visualizing, analyzing, and interpreting Calling Cards data

MOTIVATION

Unraveling the transcriptional programs that control how cells divide, differentiate, and respond to their environments requires a precise understanding of transcription factors' (TFs) DNA-binding activities. Calling cards (CC) technology uses transposons to capture transient TF binding events at one instant in time and then read them out at a later time. This methodology can also be used to simultaneously measure TF binding and mRNA expression from single-cell CC and to record and integrate TF binding events across time in any cell type of interest without the need for purification. Despite these advantages, there has been a lack of dedicated bioinformatics tools for the detailed analysis of CC data.

RESULTS

We introduce Pycallingcards, a comprehensive Python module specifically designed for the analysis of single-cell and bulk CC data across multiple species. Pycallingcards introduces two innovative peak callers, CCcaller and MACCs, enhancing the accuracy and speed of pinpointing TF binding sites from CC data. Pycallingcards offers a fully integrated environment for data visualization, motif finding, and comparative analysis with RNA-seq and ChIP-seq datasets. To illustrate its practical application, we have reanalyzed previously published mouse cortex and glioblastoma datasets. This analysis revealed novel cell-type-specific binding sites and potential sex-linked TF regulators, furthering our understanding of TF binding and gene expression relationships. Thus, Pycallingcards, with its user-friendly design and seamless interface with the Python data science ecosystem, stands as a critical tool for advancing the analysis of TF functions via CC data.

AVAILABILITY AND IMPLEMENTATION

Pycallingcards can be accessed on the GitHub repository: https://github.com/The-Mitra-Lab/pycallingcards.

Calling Cards: A Customizable Platform to Longitudinally Record Protein-DNA Interactions Over Time in Cells and Tissues

Calling Cards is a platform technology to record a cumulative history of transient protein-DNA interactions in the genome of genetically targeted cell types. The record of these interactions is recovered by next-generation sequencing. Compared with other genomic assays, readouts of which provide a snapshot at the time of harvest, Calling Cards enables correlation of historical molecular states to eventual outcomes or phenotypes. To achieve this, Calling Cards uses the piggyBac transposase to insert self-reporting transposon "Calling Cards" into the genome, leaving permanent marks at interaction sites. Calling Cards can be deployed in a variety of in vitro and in vivo biological systems to study gene regulatory networks involved in development, aging, and disease. Out of the box, it assesses enhancer usage but can be adapted to profile-specific transcription factor (TF) binding with custom TF-piggyBac fusion proteins. The Calling Cards workflow has five main stages: delivery of Calling Cards reagents, sample preparation, library preparation, sequencing, and data analysis. Here, we first present a comprehensive guide for experimental design, reagent selection, and optional customization of the platform to study additional TFs. Then, we provide an updated protocol for the five steps, using reagents that improve throughput and decrease costs, including an overview of a newly deployed computational pipeline. This protocol is designed for users with basic molecular biology experience to process samples into sequencing libraries in 2 days. Familiarity with bioinformatic analysis and command line tools is required to set up the pipeline in a high-performance computing environment and to conduct downstream analyses. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Preparation and delivery of Calling Cards reagents Support Protocol 1: Next-generation sequencing quantification of barcode distribution within self-reporting transposon plasmid pool and adeno-associated virus genome Basic Protocol 2: Sample collection and RNA purification Support Protocol 2: Library density quantitative PCR Basic Protocol 3: Sequencing library preparation Basic Protocol 4: Library pooling and sequencing Basic Protocol 5: Data analysis.

Calling Cards: a customizable platform to longitudinally record protein-DNA interactions over time in cells and tissues

Calling Cards is a platform technology to record a cumulative history of transient protein-DNA interactions in the genome of genetically targeted cell types. The record of these interactions is recovered by next generation sequencing. Compared to other genomic assays, whose readout provides a snapshot at the time of harvest, Calling Cards enables correlation of historical molecular states to eventual outcomes or phenotypes. To achieve this, Calling Cards uses the piggyBac transposase to insert self-reporting transposon (SRT) "Calling Cards" into the genome, leaving permanent marks at interaction sites. Calling Cards can be deployed in a variety of in vitro and in vivo biological systems to study gene regulatory networks involved in development, aging, and disease. Out of the box, it assesses enhancer usage but can be adapted to profile specific transcription factor binding with custom transcription factor (TF)-piggyBac fusion proteins. The Calling Cards workflow has five main stages: delivery of Calling Card reagents, sample preparation, library preparation, sequencing, and data analysis. Here, we first present a comprehensive guide for experimental design, reagent selection, and optional customization of the platform to study additional TFs. Then, we provide an updated protocol for the five steps, using reagents that improve throughput and decrease costs, including an overview of a newly deployed computational pipeline. This protocol is designed for users with basic molecular biology experience to process samples into sequencing libraries in 1-2 days. Familiarity with bioinformatic analysis and command line tools is required to set up the pipeline in a high-performance computing environment and to conduct downstream analyses. Basic Protocol 1: Preparation and delivery of Calling Cards reagentsBasic Protocol 2: Sample preparationBasic Protocol 3: Sequencing library preparationBasic Protocol 4: Library pooling and sequencingBasic Protocol 5: Data analysis.

NetProphet 3: a machine learning framework for transcription factor network mapping and multi-omics integration

MOTIVATION

Many methods have been proposed for mapping the targets of transcription factors (TFs) from gene expression data. It is known that combining outputs from multiple methods can improve performance. To date, outputs have been combined by using either simplistic formulae, such as geometric mean, or carefully hand-tuned formulae that may not generalize well to new inputs. Finally, the evaluation of accuracy has been challenging due to the lack of genome-scale, ground-truth networks.

RESULTS

We developed NetProphet3, which combines scores from multiple analyses automatically, using a tree boosting algorithm trained on TF binding location data. We also developed three independent, genome-scale evaluation metrics. By these metrics, NetProphet3 is more accurate than other commonly used packages, including NetProphet 2.0, when gene expression data from direct TF perturbations are available. Furthermore, its integration mode can forge a consensus network from gene expression data and TF binding location data.

AVAILABILITY AND IMPLEMENTATION

All data and code are available at https://zenodo.org/record/7504131#.Y7Wu3i-B2x8.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Measuring transcription factor binding and gene expression using barcoded self-reporting transposon calling cards and transcriptomes

Calling cards technology using self-reporting transposons enables the identification of DNA-protein interactions through RNA sequencing. Although immensely powerful, current implementations of calling cards in bulk experiments on populations of cells are technically cumbersome and require many replicates to identify independent insertions into the same genomic locus. Here, we have drastically reduced the cost and labor requirements of calling card experiments in bulk populations of cells by introducing a DNA barcode into the calling card itself. An additional barcode incorporated during reverse transcription enables simultaneous transcriptome measurement in a facile and affordable protocol. We demonstrate that barcoded self-reporting transposons recover in vitro binding sites for four basic helix-loop-helix transcription factors with important roles in cell fate specification: ASCL1, MYOD1, NEUROD2 and NGN1. Further, simultaneous calling cards and transcriptional profiling during transcription factor overexpression identified both binding sites and gene expression changes for two of these factors. Lastly, we demonstrated barcoded calling cards can record binding in vivo in the mouse brain. In sum, RNA-based identification of transcription factor binding sites and gene expression through barcoded self-reporting transposon calling cards and transcriptomes is an efficient and powerful method to infer gene regulatory networks in a population of cells.

Genome-wide protein-DNA interaction site mapping in bacteria using a double-stranded DNA-specific cytosine deaminase

DNA-protein interactions are central to fundamental cellular processes, yet widely implemented technologies for measuring these interactions on a genome scale in bacteria are laborious and capture only a snapshot of binding events. We devised a facile method for mapping DNA-protein interaction sites in vivo using the double-stranded DNA-specific cytosine deaminase toxin DddA. In 3D-seq (DddA-sequencing), strains containing DddA fused to a DNA-binding protein of interest accumulate characteristic mutations in DNA sequence adjacent to sites occupied by the DNA-bound fusion protein. High-depth sequencing enables detection of sites of increased mutation frequency in these strains, yielding genome-wide maps of DNA-protein interaction sites. We validated 3D-seq for four transcription regulators in two bacterial species, Pseudomonas aeruginosa and Escherichia coli. We show that 3D-seq offers ease of implementation, the ability to record binding event signatures over time and the capacity for single-cell resolution.

Here and there: the double-side transgene localization

Random transgene integration is a powerful tool for developing new genome-wide screening approaches. These techniques have already been used for functional gene annotation by transposon-insertion sequencing, for identif ication of transcription factor binding sites and regulatory sequences, and for dissecting chromatin position effects. Precise localization of transgenes and accurate artifact f iltration are essential for this type of method. To date, many mapping assays have been developed, including Inverse-PCR, TLA, LAM-PCR, and splinkerette PCR. However, none of them is able to ensure localization of both transgene’s f lanking regions simultaneously, which would be necessary for some applications. Here we proposed a cheap and simple NGS-based approach that overcomes this limitation. The developed assay requires using intentionally designed vectors that lack recognition sites of one or a set of restriction enzymes used for DNA fragmentation. By looping and sequencing these DNA fragments, we obtain special data that allows us to link the two f lanking regions of the transposon. This can be useful for precise insertion mapping and for screening approaches in the f ield of chromosome engineering, where chromosomal recombination events between transgenes occur in a cell population. To demonstrate the method’s feasibility, we applied it for mapping SB transposon integration in the human HAP1 cell line. Our technique allowed us to eff iciently localize genomic transposon integrations, which was conf irmed via PCR analysis. For practical application of this approach, we proposed a set of recommendations and a normalization strategy. The developed method can be used for multiplex transgene localization and detection of rearrangements between them.

The qBED track: a novel genome browser visualization for point processes

SUMMARY

Transposon calling cards is a genomic assay for identifying transcription factor binding sites in both bulk and single cell experiments. Here, we describe the qBED format, an open, text-based standard for encoding and analyzing calling card data. In parallel, we introduce the qBED track on the WashU Epigenome Browser, a novel visualization that enables researchers to inspect calling card data in their genomic context. Finally, through examples, we demonstrate that qBED files can be used to visualize non-calling card datasets, such as Combined Annotation-Dependent Depletion scores and GWAS/eQTL hits, and thus may have broad utility to the genomics community.

AVAILABILITY AND IMPLEMENTATION

The qBED track is available on the WashU Epigenome Browser (http://epigenomegateway.wustl.edu/browser), beginning with version 46. Source code for the WashU Epigenome Browser with qBED support is available on GitHub (http://github.com/arnavm/eg-react and http://github.com/lidaof/eg-react). A complete definition of the qBED format is available as part of the WashU Epigenome Browser documentation (https://eg.readthedocs.io/en/latest/tracks.html#qbed-track). We have also released a tutorial on how to upload qBED data to the browser (http://dx.doi.org/10.17504/protocols.io.bca8ishw).

Recording of DNA-binding events reveals the importance of a repurposed Candida albicans regulatory network for gut commensalism

Candida albicans is a fungal component of the human gut microbiota and an opportunistic pathogen. C. albicans transcription factors (TFs), Wor1 and Efg1, are master regulators of an epigenetic switch required for fungal mating that also control colonization of the mammalian gut. We show that additional mating regulators, WOR2, WOR3, WOR4, AHR1, CZF1, and SSN6, also influence gut commensalism. Using Calling Card-seq to record Candida TF DNA-binding events in the host, we examine the role and relationships of these regulators during murine gut colonization. By comparing in-host transcriptomes of regulatory mutants with enhanced versus diminished commensal fitness, we also identify a set of candidate commensalism effectors. These include Cht2, a GPI-linked chitinase whose gene is bound by Wor1, Czf1, and Efg1 in vivo, that we show promotes commensalism. Thus, the network required for a C. albicans sexual switch is biochemically active in the host intestine and repurposed to direct commensalism.

Brd4-bound enhancers drive cell-intrinsic sex differences in glioblastoma

Consistent sex differences in incidence and outcome have been reported in numerous cancers including brain tumors. GBM, the most common and aggressive primary brain tumor, occurs with higher incidence and shorter survival in males compared to females. Brd4 is essential for regulating transcriptome-wide gene expression and specifying cell identity, including that of GBM. We report that sex-biased Brd4 activity drives sex differences in GBM and renders male and female tumor cells differentially sensitive to BET inhibitors. The observed sex differences in BETi treatment strongly indicate that sex differences in disease biology translate into sex differences in therapeutic responses. This has critical implications for clinical use of BET inhibitors further affirming the importance of inclusion of sex as a biological variable.

Sex can be an important determinant of cancer phenotype, and exploring sex-biased tumor biology holds promise for identifying novel therapeutic targets and new approaches to cancer treatment. In an established isogenic murine model of glioblastoma (GBM), we discovered correlated transcriptome-wide sex differences in gene expression, H3K27ac marks, large Brd4-bound enhancer usage, and Brd4 localization to Myc and p53 genomic binding sites. These sex-biased gene expression patterns were also evident in human glioblastoma stem cells (GSCs). These observations led us to hypothesize that Brd4-bound enhancers might underlie sex differences in stem cell function and tumorigenicity in GBM. We found that male and female GBM cells exhibited sex-specific responses to pharmacological or genetic inhibition of Brd4. Brd4 knockdown or pharmacologic inhibition decreased male GBM cell clonogenicity and in vivo tumorigenesis while increasing both in female GBM cells. These results were validated in male and female patient-derived GBM cell lines. Furthermore, analysis of the Cancer Therapeutic Response Portal of human GBM samples segregated by sex revealed that male GBM cells are significantly more sensitive to BET (bromodomain and extraterminal) inhibitors than are female cells. Thus, Brd4 activity is revealed to drive sex differences in stem cell and tumorigenic phenotypes, which can be abrogated by sex-specific responses to BET inhibition. This has important implications for the clinical evaluation and use of BET inhibitors.

The piggyBac-derived protein 5 (PGBD5) transposes both the closely and the distantly related piggyBac-like elements Tcr-pble and Ifp2

The vertebrate piggyBac derived transposase 5 (PGBD5) encodes a domesticated transposase, which is active and able to transpose its distantly related piggyBac-like element (pble), Ifp2. This raised the question whether PGBD5 would be more effective at mobilizing a phylogenetically closely related pble element. We aimed to identify the pble most closely related to the pgbd5 gene. We updated the landscape of vertebrate pgbd genes to develop efficient filters and identify the most closely related pble to each of these genes. We found that Tcr-pble is phylogenetically the closest pble to the pgbd5 gene. Furthermore, we evaluated the capacity of two murine and human PGBD5 isoforms, Mm523 and Hs524, to transpose both Tcr-pble and Ifp2 elements. We found that both pbles could be transposed by Mm523 with similar efficiency. However, integrations of both pbles occurred through both proper transposition and improper PGBD5-dependent recombination. This suggested that the ability of PGBD5 to bind both pbles may not be based on the primary sequence of element ends, but may involve recognition of inner DNA motifs, possibly related to palindromic repeats. In agreement with this hypothesis, we identified internal palindromic repeats near the end of 24 pble sequences, which display distinct sequences.

The C-terminal Domain of piggyBac Transposase Is Not Required for DNA Transposition

PiggyBac(PB)-like elements (pble) are members of a eukaryotic DNA transposon family. This family is of interest to evolutionary genomics because pble transposases have been domesticated at least 9 times in vertebrates. The amino acid sequence of pble transposases can be split into three regions: an acidic N-terminal domain (~100 aa), a central domain (~400 aa) containing a DD[D/E] catalytic triad, and a cysteine-rich domain (CRD; ~90 aa). Two recent reports suggested that a functional CRD is required for pble transposase activity. Here we found that two CRD-deficient pble transposases, a PB variant and an isoform encoded by the domesticated PB-derived vertebrate transposase gene 5 (pgbd5) trigger transposition of the Ifp2 pble. When overexpressed in HeLa cells, these CRD-deficient transposases can insert Ifp2 elements with proper and improper transposon ends, associated with deleterious effects on cells. Finally, we found that mouse CRD-deficient transposase Pgbd5, as well as PB, do not insert pbles at random into chromosomes. Transposition events occurred more often in genic regions, in the neighbourhood of the transcription start sites and were often found in genes predominantly expressed in the human central nervous system.

High-Resolution Genome-Wide Occupancy in Candida spp. Using ChEC-seq

To persist in their dynamic human host environments, fungal pathogens must sense and adapt by modulating their gene expression to fulfill their cellular needs. Understanding transcriptional regulation on a global scale would uncover cellular processes linked to persistence and virulence mechanisms that could be targeted for antifungal therapeutics. Infections associated with the yeast Candida albicans, a highly prevalent fungal pathogen, and the multiresistant related species Candida auris are becoming a serious public health threat. To define the set of a gene regulated by a transcriptional regulator in C. albicans, chromatin immunoprecipitation (ChIP)-based techniques, including ChIP with microarray technology (ChIP-chip) or ChIP-DNA sequencing (ChIP-seq), have been widely used. Here, we describe a new set of PCR-based micrococcal nuclease (MNase)-tagging plasmids for C. albicans and other Candida spp. to determine the genome-wide location of any transcriptional regulator of interest using chromatin endogenous cleavage (ChEC) coupled to high-throughput sequencing (ChEC-seq). The ChEC procedure does not require protein-DNA cross-linking or sonication, thus avoiding artifacts related to epitope masking or the hyper-ChIPable euchromatic phenomenon. In a proof-of-concept application of ChEC-seq, we provided a high-resolution binding map of the SWI/SNF chromatin remodeling complex, a master regulator of fungal fitness in C. albicans, in addition to the transcription factor Nsi1 that is an ortholog of the DNA-binding protein Reb1 for which genome-wide occupancy was previously established in Saccharomyces cerevisiae The ChEC-seq procedure described here will allow a high-resolution genomic location definition which will enable a better understanding of transcriptional regulatory circuits that govern fungal fitness and drug resistance in these medically important fungi.IMPORTANCE Systemic fungal infections caused by Candida albicans and the "superbug" Candida auris are becoming a serious public health threat. The ability of these yeasts to cause disease is linked to their faculty to modulate the expression of genes that mediate their escape from the immune surveillance and their persistence in the different unfavorable niches within the host. Comprehensive knowledge on gene expression control of fungal fitness is consequently an interesting framework for the identification of essential infection processes that could be hindered by chemicals as potential therapeutics. Here, we expanded the use of ChEC-seq, a technique that was initially developed in the yeast model Saccharomyces cerevisiae to identify genes that are modulated by a transcriptional regulator, in pathogenic yeasts from the genus Candida This robust technique will allow a better characterization of key gene expression regulators and their contribution to virulence and antifungal resistance in these pathogenic yeasts.

Quantitative analysis of transcription factor binding and expression using calling cards reporter arrays

We report a tool, Calling Cards Reporter Arrays (CCRA), that measures transcription factor (TF) binding and the consequences on gene expression for hundreds of synthetic promoters in yeast. Using Cbf1p and MAX, we demonstrate that the CCRA method is able to detect small changes in binding free energy with a sensitivity comparable to in vitro methods, enabling the measurement of energy landscapes in vivo. We then demonstrate the quantitative analysis of cooperative interactions by measuring Cbf1p binding at synthetic promoters with multiple sites. We find that the cooperativity between Cbf1p dimers varies sinusoidally with a period of 10.65 bp and energetic cost of 1.37 K_BT for sites that are positioned ‘out of phase’. Finally, we characterize the binding and expression of a group of TFs, Tye7p, Gcr1p and Gcr2p, that act together as a ‘TF collective’, an important but poorly characterized model of TF cooperativity. We demonstrate that Tye7p often binds promoters without its recognition site because it is recruited by other collective members, whereas these other members require their recognition sites, suggesting a hierarchy where these factors recruit Tye7p but not vice versa. Our experiments establish CCRA as a useful tool for quantitative investigations into TF binding and function.

Self-Reporting Transposons Enable Simultaneous Readout of Gene Expression and Transcription Factor Binding in Single Cells

Cellular heterogeneity confounds in situ assays of transcription factor (TF) binding. Single-cell RNA sequencing (scRNA-seq) deconvolves cell types from gene expression, but no technology links cell identity to TF binding sites (TFBS) in those cell types. We present self-reporting transposons (SRTs) and use them in single-cell calling cards (scCC), a novel assay for simultaneously measuring gene expression and mapping TFBS in single cells. The genomic locations of SRTs are recovered from mRNA, and SRTs deposited by exogenous, TF-transposase fusions can be used to map TFBS. We then present scCC, which map SRTs from scRNA-seq libraries, simultaneously identifying cell types and TFBS in those same cells. We benchmark multiple TFs with this technique. Next, we use scCC to discover BRD4-mediated cell-state transitions in K562 cells. Finally, we map BRD4 binding sites in the mouse cortex at single-cell resolution, establishing a new method for studying TF biology in situ.

A viral toolkit for recording transcription factor-DNA interactions in live mouse tissues

Transcription factors (TFs) enact precise regulation of gene expression through site-specific, genome-wide binding. Common methods for TF-occupancy profiling, such as chromatin immunoprecipitation, are limited by requirement of TF-specific antibodies and provide only end-point snapshots of TF binding. Alternatively, TF-tagging techniques, in which a TF is fused to a DNA-modifying enzyme that marks TF-binding events across the genome as they occur, do not require TF-specific antibodies and offer the potential for unique applications, such as recording of TF occupancy over time and cell type specificity through conditional expression of the TF-enzyme fusion. Here, we create a viral toolkit for one such method, calling cards, and demonstrate that these reagents can be delivered to the live mouse brain and used to report TF occupancy. Further, we establish a Cre-dependent calling cards system and, in proof-of-principle experiments, show utility in defining cell type-specific TF profiles and recording and integrating TF-binding events across time. This versatile approach will enable unique studies of TF-mediated gene regulation in live animal models.

Dual threshold optimization and network inference reveal convergent evidence from TF binding locations and TF perturbation responses

A high-confidence map of the direct, functional targets of each transcription factor (TF) requires convergent evidence from independent sources. Two significant sources of evidence are TF binding locations and the transcriptional responses to direct TF perturbations. Systematic data sets of both types exist for yeast and human, but they rarely converge on a common set of direct, functional targets for a TF. Even the few genes that are both bound and responsive may not be direct functional targets. Our analysis shows that when there are many nonfunctional binding sites and many indirect targets, nonfunctional sites are expected to occur in the cis-regulatory DNA of indirect targets by chance. To address this problem, we introduce dual threshold optimization (DTO), a new method for setting significance thresholds on binding and perturbation-response data, and show that it improves convergence. It also enables comparison of binding data to perturbation-response data that have been processed by network inference algorithms, which further improves convergence. The combination of dual threshold optimization and network inference greatly expands the high-confidence TF network map in both yeast and human. Next, we analyze a comprehensive new data set measuring the transcriptional response shortly after inducing overexpression of a yeast TF. We also present a new yeast binding location data set obtained by transposon calling cards and compare it to recent ChIP-exo data. These new data sets improve convergence and expand the high-confidence network synergistically.

RNA-guided piggyBac transposition in human cells

Safer and more efficient methods for directing therapeutic genes to specific sequences could increase the repertoire of treatable conditions. Many current approaches act passively, first initiating a double-stranded break, then relying on host repair to uptake donor DNA. Alternatively, we delivered an actively integrating transposase to the target sequence to initiate gene insertion. We fused the hyperactive piggyBac transposase to the highly specific, catalytically dead SpCas9-HF1 (dCas9) and designed guide RNAs (gRNAs) to the CCR5 safe harbor sequence. We introduced mutations to the native DNA-binding domain of piggyBac to reduce non-specific binding of the transposase and cause the fusion protein to favor binding by dCas9. This strategy enabled us, for the first time, to direct transposition to the genome using RNA. We showed that increasing the number of gRNAs improved targeting efficiency. Interestingly, over half of the recovered insertions were found at a single TTAA hotspot. We also found that the fusion increased the error rate at the genome-transposon junction. We isolated clonal cell lines containing a single insertion at CCR5 and demonstrated long-term expression from this locus. These vectors expand the utility of the piggyBac system for applications in targeted gene addition for biomedical research and gene therapy.

In Vivo Piggybac-Based Gene Delivery towards Murine Pancreatic Parenchyma Confers Sustained Expression of Gene of Interest

The pancreas is a glandular organ that functions in the digestive system and endocrine system of vertebrates. The most common disorders involving the pancreas are diabetes, pancreatitis, and pancreatic cancer. In vivo gene delivery targeting the pancreas is important for preventing or curing such diseases and for exploring the biological function of genes involved in the pathogenesis of these diseases. Our previous experiments demonstrated that adult murine pancreatic cells can be efficiently transfected by exogenous plasmid DNA following intraparenchymal injection and subsequent in vivo electroporation using tweezer-type electrodes. Unfortunately, the induced gene expression was transient. Transposon-based gene delivery, such as that facilitated by piggyBac (PB), is known to confer stable integration of a gene of interest (GOI) into host chromosomes, resulting in sustained expression of the GOI. In this study, we investigated the use of the PB transposon system to achieve stable gene expression when transferred into murine pancreatic cells using the above-mentioned technique. Expression of the GOI (coding for fluorescent protein) continued for at least 1.5 months post-gene delivery. Splinkerette-PCR-based analysis revealed the presence of the consensus sequence TTAA at the junctional portion between host chromosomes and the transgenes; however, this was not observed in all samples. This plasmid-based PB transposon system enables constitutive expression of the GOI in pancreas for potential therapeutic and biological applications.

The impact of transposable element activity on therapeutically relevant human stem cells

Human stem cells harbor significant potential for basic and clinical translational research as well as regenerative medicine. Currently ~ 3000 adult and ~ 30 pluripotent stem cell-based, interventional clinical trials are ongoing worldwide, and numbers are increasing continuously. Although stem cells are promising cell sources to treat a wide range of human diseases, there are also concerns regarding potential risks associated with their clinical use, including genomic instability and tumorigenesis concerns. Thus, a deeper understanding of the factors and molecular mechanisms contributing to stem cell genome stability are a prerequisite to harnessing their therapeutic potential for degenerative diseases. Chemical and physical factors are known to influence the stability of stem cell genomes, together with random mutations and Copy Number Variants (CNVs) that accumulated in cultured human stem cells. Here we review the activity of endogenous transposable elements (TEs) in human multipotent and pluripotent stem cells, and the consequences of their mobility for genomic integrity and host gene expression. We describe transcriptional and post-transcriptional mechanisms antagonizing the spread of TEs in the human genome, and highlight those that are more prevalent in multipotent and pluripotent stem cells. Notably, TEs do not only represent a source of mutations/CNVs in genomes, but are also often harnessed as tools to engineer the stem cell genome; thus, we also describe and discuss the most widely applied transposon-based tools and highlight the most relevant areas of their biomedical applications in stem cells. Taken together, this review will contribute to the assessment of the risk that endogenous TE activity and the application of genetically engineered TEs constitute for the biosafety of stem cells to be used for substitutive and regenerative cell therapies.

Enzymatic methods for genome-wide profiling of protein binding sites

Genome-wide mapping of protein-DNA interactions is a staple approach in many areas of modern molecular biology. Genome-wide profiles of protein-binding sites are most commonly generated by chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq). Although ChIP-seq has played a central role in studying genome-wide protein binding, recent work has highlighted systematic biases in the technique that warrant technical and interpretive caution and underscore the need for orthogonal techniques to both confirm the results of ChIP-seq studies and uncover new insights not accessible to ChIP. Several such techniques, based on genetic or immunological targeting of enzymatic activity to specific genomic loci, have been developed. Here, we review the development, applications and future prospects of these methods as complements to ChIP-based approaches and as powerful techniques in their own right.

A massively parallel reporter assay dissects the influence of chromatin structure on cis-regulatory activity

A gene's position in the genome can profoundly affect its expression because regional differences in chromatin modulate the activity of locally acting cis-regulatory sequences (CRSs). Here we study how CRSs and regional chromatin act in concert on a genome-wide scale. We present a massively parallel reporter gene assay that measures the activities of hundreds of different CRSs, each integrated at many specific genomic locations. Although genome location strongly affected CRS activity, the relative strengths of CRSs were maintained at all chromosomal locations. The intrinsic activities of CRSs also correlated with their activities in plasmid-based assays. We explain our data with a quantitative model in which expression levels are set by independent contributions from local CRSs and the regional chromatin environment, rather than by more complex sequence- or protein-specific interactions between these two factors. The methods we present will help investigators determine when regulatory information is integrated in a modular fashion and when regulatory sequences interact in more complex ways.

Transposase mapping identifies the genomic targets of BAP1 in uveal melanoma

Background

BAP1 is a histone deubiquitinase that acts as a tumor and metastasis suppressor associated with disease progression in human cancer. We have used the “Calling Card System” of transposase-directed transposon insertion mapping to identify the genomic targets of BAP1 in uveal melanoma (UM). This system was developed to identify the genomic loci visited by transcription factors that bind directly to DNA; our study is the first use of the system with a chromatin-remodeling factor that binds to histones but does not interact directly with DNA.

Methods

The transposase piggyBac (PBase) was fused to BAP1 and expressed in OCM-1A UM cells. The insertion of transposons near BAP1 binding sites in UM cells were identified by genomic sequencing. We also examined RNA expression in the same OCM-1A UM cells after BAP1 depletion to identify BAP1 binding sites associated with BAP1-responsive genes. Sets of significant genes were analyzed for common pathways, transcription factor binding sites, and ability to identify molecular tumor classes.

Results

We found a strong correlation between multiple calling-card transposon insertions targeted by BAP1-PBase and BAP1-responsive expression of adjacent genes. BAP1-bound genomic loci showed narrow distributions of insertions and were near transcription start sites, consistent with recruitment of BAP1 to these sites by specific DNA-binding proteins. Sequence consensus analysis of BAP1-bound sites showed enrichment of motifs specific for YY1, NRF1 and Ets transcription factors, which have been shown to interact with BAP1 in other cell types. Further, a subset of the BAP1 genomic target genes was able to discriminate aggressive tumors in published gene expression data from primary UM tumors.

Conclusions

The calling card methodology works equally well for chromatin regulatory factors that do not interact directly with DNA as for transcription factors. This technique has generated a new and expanded list of BAP1 targets in UM that provides important insight into metastasis pathways and identifies novel potential therapeutic targets.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0424-0) contains supplementary material, which is available to authorized users.

MINCE-Seq: Mapping In Vivo Nascent Chromatin with EdU and Sequencing

The epigenome has been mapped in different cell types to understand the relationship between the chromatin landscape and the control of gene expression. Most mapping studies profile a large population of cells in various stages of the cell cycle, which results in an average snapshot of the chromatin landscape. However, chromatin is highly dynamic, undergoing rapid changes during active processes such as replication, transcription, repair, and remodeling. Hence, we need methods to map chromatin as a function of time. To address this problem in the context of replication, we developed the method MINCE-seq (Mapping In vivo Nascent Chromatin with EdU and sequencing). MINCE-seq is a genome-wide method that uses the passage of replication fork as a starting point to map the chromatin landscape as a function of time. MINCE-seq can measure chromatin dynamics in a time scale of minutes and at the resolution of individual nucleosome positions and transcription factor-binding sites genome-wide.

O Cdc7 kinase where art thou?

Although Cdc7 protein kinase is important for regulating DNA replication in all eukaryotes and is a target for cancer therapy, it has never been localized in cells. Recently, a novel molecular genomic method used by our laboratory to localize Cdc7 to regions of chromosomes. Originally, mutations in the CDC7 gene were found in the classic cdc mutant collection of Hartwell et al. (Genetics 74:267-286, 1973). The CDC7 gene was found to encode a protein kinase called DDK that has been studied for many years, establishing its precise role in the initiation of DNA replication at origins. Recently, clinical studies are underway with DDK inhibitors against DDK in cancer patients. However, the conundrum is that Cdc7 has never been detected at origins of replication even though many studies have suggested it should be there. We used "Calling Card" system in which DNA binding proteins are localized to the genome via retrotransposon insertion and deep-sequencing methods. We have shown that Cdc7 localizes at many regions of the genome and was enriched at functional origins of replication. These results are consistent with DDK's role in many additional genomic processes including mutagenesis, chromatid cohesion, and meiotic recombination. Thus, the main conclusion from our studies is that Cdc7 kinase is found at many locations in the genome, but is enriched at functional origins of replication. Furthermore, we propose that application of the Calling Card system to other eukaryotes should be useful in identification of functional origins in other eukaryotic cells.

Transposons As Tools for Functional Genomics in Vertebrate Models

Genetic tools and mutagenesis strategies based on transposable elements are currently under development with a vision to link primary DNA sequence information to gene functions in vertebrate models. By virtue of their inherent capacity to insert into DNA, transposons can be developed into powerful tools for chromosomal manipulations. Transposon-based forward mutagenesis screens have numerous advantages including high throughput, easy identification of mutated alleles, and providing insight into genetic networks and pathways based on phenotypes. For example, the Sleeping Beauty transposon has become highly instrumental to induce tumors in experimental animals in a tissue-specific manner with the aim of uncovering the genetic basis of diverse cancers. Here, we describe a battery of mutagenic cassettes that can be applied in conjunction with transposon vectors to mutagenize genes, and highlight versatile experimental strategies for the generation of engineered chromosomes for loss-of-function as well as gain-of-function mutagenesis for functional gene annotation in vertebrate models, including zebrafish, mice, and rats.

Characterization and anti-inflammation role of swine IFITM3 gene

IFITM3 is involved in cell adhesion, apoptosis, immune, and antivirus activity. Furthermore, IFITM3 gene has been considered as a preferential marker for inflammatory diseases, and positive correlation to pathological grades. Therefore, we assumed that IFITM3 was regulated by different signal pathways. To better understand IFITM3 function in inflammatory response, we cloned swine IFITM3 gene, and detected IFITM3 distribution in tissues, as well as characterized this gene. Results indicated that the length of swine IFITM3 gene was 438 bp, encoding 145 amino acids. IFITM3 gene expression abundance was higher in spleen and lungs. Moreover, we next constructed the eukaryotic expression vector PBIFM3 and transfected into PK15 cells, finally obtained swine IFITM3 gene stable expression cell line. Meanwhile, we explored the effects of LPS on swine IFITM3 expression. Results showed that LPS increased IFITM3 mRNA abundance and exhibited time-dependent effect for LPS treatment. To further demonstrate the mechanism that IFITM3 regulated type I IFNs production, we also detected the important molecules expression of TLR4 signaling pathway. In transfected and non-transfected IFITM3 PK15 cells, LPS exacerbated the relative expression of TLR4-NFκB signaling molecules. However, the IFITM3 overexpression suppressed the inflammatory development of PK15 cells. In conclusion, these data indicated that the overexpression of swine IFITM3 could decrease the inflammatory response through TLR4 signaling pathway, and participate in type I interferon production. These findings may lead to an improved understanding of the biological function of IFITM3 in inflammation.

An optimized, broadly applicable piggyBac transposon induction system

The piggyBac (PB) transposon has been used in a number of biological applications. The insertion of PB transposons into the genome can disrupt genes or regulatory regions, impacting cellular function, so for many experiments it is important that PB transposition is tightly controlled. Here, we systematically characterize three methods for the post-translational control of the PB transposon in four cell lines. We investigated fusions of the PB transposase with ERT2 and two degradation domains (FKBP-DD, DHFR-DD), in multiple orientations, and determined (i) the fold-induction achieved, (ii) the absolute transposition efficiency of the activated construct and (iii) the effects of two inducer molecules on cellular transcription and function. We found that the FKBP-DD confers the PB transposase with a higher transposition activity and better dynamic range than can be achieved with the other systems. In addition, we found that the FKBP-DD regulates transposon activity in a reversible and dose-dependent manner. Finally, we showed that Shld1, the chemical inducer of FKBP-DD, does not interfere with stem cell differentiation, whereas tamoxifen has significant effects. We believe the FKBP-based PB transposon induction will be useful for transposon-mediated genome engineering, insertional mutagenesis and the genome-wide mapping of transcription factor binding.

Swine IRF3/IRF7 attenuates inflammatory responses through TLR4 signaling pathway

To explore the role of IRF3/IRF7 during inflammatory responses, we investigated the effects of swine IRF3/IRF7 on TLR4 signaling pathway and inflammatory factors expression in porcine kidney epithelial PK15 cell lines. We successfully constructed eukaryotic vectors PB-IRF3 and PB-IRF7, transfected these vectors into PK15 cells and observed GFP under a fluorescence microscope. In addition, RT-PCR was also used to detect transfection efficiency. We found that IRF3/IRF7 was efficiently overexpressed in PK15 cells. Moreover, we evaluated the effects of IRF3/IRF7 on the TLR4 signaling pathway and inflammatory factors by RT-PCR. Transfected cells were treated with lipopolysaccharide (LPS) alone, or in combination with a TBK1 inhibitor (LiCl). We revealed that IRF3/IRF7 enhanced IFNα production, and decreased IL-6 mRNA expression. Blocking the TBK1 pathway, inhibited the changes in IFNα, but not IL-6 mRNA. This illustrated that IRF3/IRF7 enhanced IFNα production through TLR4/TBK1 signaling pathway and played an anti-inflammatory role, while IRF3/IRF7 decreased IL-6 expression independent of the TBK1 pathway. Trends in MyD88, TRAF6, TBK1 and NFκB mRNA variation were similar in all treatments. LPS increased MyD88, TRAF6, TBK1 and NFκB mRNA abundance in PBR3/PBR7 and PBv cells, while LiCl blocked the LPS-mediated effects. The levels of these four factors in PBR3/PBR7 cells were higher than those in PBv. These results demonstrated that IRF3/IRF7 regulated the inflammatory response through the TLR4 signaling pathway. Overexpression of swine IRF3/IRF7 in PK15 cells induced type I interferons production, and attenuated inflammatory responses through TLR4 signaling pathway.

Model-based transcriptome engineering promotes a fermentative transcriptional state in yeast

The ability to rationally manipulate the transcriptional states of cells would be of great use in medicine and bioengineering. We have developed an algorithm, NetSurgeon, which uses genome-wide gene-regulatory networks to identify interventions that force a cell toward a desired expression state. We first validated NetSurgeon extensively on existing datasets. Next, we used NetSurgeon to select transcription factor deletions aimed at improving ethanol production in Saccharomyces cerevisiae cultures that are catabolizing xylose. We reasoned that interventions that move the transcriptional state of cells using xylose toward that of cells producing large amounts of ethanol from glucose might improve xylose fermentation. Some of the interventions selected by NetSurgeon successfully promoted a fermentative transcriptional state in the absence of glucose, resulting in strains with a 2.7-fold increase in xylose import rates, a 4-fold improvement in xylose integration into central carbon metabolism, or a 1.3-fold increase in ethanol production rate. We conclude by presenting an integrated model of transcriptional regulation and metabolic flux that will enable future efforts aimed at improving xylose fermentation to prioritize functional regulators of central carbon metabolism.

Identifying Cancer Driver Genes Using Replication-Incompetent Retroviral Vectors

Identifying novel genes that drive tumor metastasis and drug resistance has significant potential to improve patient outcomes. High-throughput sequencing approaches have identified cancer genes, but distinguishing driver genes from passengers remains challenging. Insertional mutagenesis screens using replication-incompetent retroviral vectors have emerged as a powerful tool to identify cancer genes. Unlike replicating retroviruses and transposons, replication-incompetent retroviral vectors lack additional mutagenesis events that can complicate the identification of driver mutations from passenger mutations. They can also be used for almost any human cancer due to the broad tropism of the vectors. Replication-incompetent retroviral vectors have the ability to dysregulate nearby cancer genes via several mechanisms including enhancer-mediated activation of gene promoters. The integrated provirus acts as a unique molecular tag for nearby candidate driver genes which can be rapidly identified using well established methods that utilize next generation sequencing and bioinformatics programs. Recently, retroviral vector screens have been used to efficiently identify candidate driver genes in prostate, breast, liver and pancreatic cancers. Validated driver genes can be potential therapeutic targets and biomarkers. In this review, we describe the emergence of retroviral insertional mutagenesis screens using replication-incompetent retroviral vectors as a novel tool to identify cancer driver genes in different cancer types.

Past Roadblocks and New Opportunities in Transcription Factor Network Mapping

One of the principal mechanisms by which cells differentiate and respond to changes in external signals or conditions is by changing the activity levels of transcription factors (TFs). This changes the transcription rates of target genes via the cell's TF network, which ultimately contributes to reconfiguring cellular state. Since microarrays provided our first window into global cellular state, computational biologists have eagerly attacked the problem of mapping TF networks, a key part of the cell's control circuitry. In retrospect, however, steady-state mRNA abundance levels were a poor substitute for TF activity levels and gene transcription rates. Likewise, mapping TF binding through chromatin immunoprecipitation proved less predictive of functional regulation and less amenable to systematic elucidation of complete networks than originally hoped. This review explains these roadblocks and the current, unprecedented blossoming of new experimental techniques built on second-generation sequencing, which hold out the promise of rapid progress in TF network mapping.

Transposon Calling Cards

Identifying the genomic targets of transcription factors is an important step in understanding the regulatory networks of gene transcription in yeast. We have developed a method that utilizes what we refer to as transposon "calling cards," in which a transcription factor directs the Ty5 retrotransposase to insert transposons into the genome adjacent to where the transcription factor binds. This method is designed to be multiplexed with many barcoded transcription factors and has the potential to decrease the labor required for the study of large numbers of transcription factors.

Genome-wide Profiling Reveals Remarkable Parallels Between Insertion Site Selection Properties of the MLV Retrovirus and the piggyBac Transposon in Primary Human CD4(+) T Cells

The inherent risks associated with vector insertion in gene therapy need to be carefully assessed. We analyzed the genome-wide distributions of Sleeping Beauty (SB) and piggyBac (PB) transposon insertions as well as MLV retrovirus and HIV lentivirus insertions in human CD4⁺ T cells with respect to a panel of 40 chromatin states. The distribution of SB transposon insertions displayed the least deviation from random, while the PB transposon and the MLV retrovirus showed unexpected parallels across all chromatin states. Both MLV and PB insertions are enriched at transcriptional start sites (TSSs) and co-localize with BRD4-associated sites. We demonstrate physical interaction between the PB transposase and bromodomain and extraterminal domain proteins (including BRD4), suggesting convergent evolution of a tethering mechanism that directs integrating genetic elements into TSSs. We detect unequal biases across the four systems with respect to targeting genes whose deregulation has been previously linked to serious adverse events in gene therapy clinical trials. The SB transposon has the highest theoretical chance of targeting a safe harbor locus in the human genome. The data underscore the significance of vector choice to reduce the mutagenic load on cells in clinical applications.

ChEC-seq kinetics discriminates transcription factor binding sites by DNA sequence and shape in vivo

Chromatin endogenous cleavage (ChEC) uses fusion of a protein of interest to micrococcal nuclease (MNase) to target calcium-dependent cleavage to specific genomic loci in vivo. Here we report the combination of ChEC with high-throughput sequencing (ChEC-seq) to map budding yeast transcription factor (TF) binding. Temporal analysis of ChEC-seq data reveals two classes of sites for TFs, one displaying rapid cleavage at sites with robust consensus motifs and the second showing slow cleavage at largely unique sites with low-scoring motifs. Sites with high-scoring motifs also display asymmetric cleavage, indicating that ChEC-seq provides information on the directionality of TF-DNA interactions. Strikingly, similar DNA shape patterns are observed regardless of motif strength, indicating that the kinetics of ChEC-seq discriminates DNA recognition through sequence and/or shape. We propose that time-resolved ChEC-seq detects both high-affinity interactions of TFs with consensus motifs and sites preferentially sampled by TFs during diffusion and sliding.

In chromatin endogenous cleavage (ChEC), micrococcal nuclease (MNase) is fused to a protein of interest and its cleavage is thus targeted to specific genomic loci in vivo. Here, the authors show that time-resolved ChEC-seq (high-throughput sequencing after ChEC) can detect DNA shape patterns regardless of motif strength.

Oligonucleotide gap-fill ligation for mutation detection and sequencing in situ

In clinical diagnostics a great need exists for targeted in situ multiplex nucleic acid analysis as the mutational status can offer guidance for effective treatment. One well-established method uses padlock probes for mutation detection and multiplex expression analysis directly in cells and tissues. Here, we use oligonucleotide gap-fill ligation to further increase specificity and to capture molecular substrates for in situ sequencing. Short oligonucleotides are joined at both ends of a padlock gap probe by two ligation events and are then locally amplified by target-primed rolling circle amplification (RCA) preserving spatial information. We demonstrate the specific detection of the A3243G mutation of mitochondrial DNA and we successfully characterize a single nucleotide variant in the ACTB mRNA in cells by in situ sequencing of RCA products generated by padlock gap-fill ligation. To demonstrate the clinical applicability of our assay, we show specific detection of a point mutation in the EGFR gene in fresh frozen and formalin-fixed, paraffin-embedded (FFPE) lung cancer samples and confirm the detected mutation by in situ sequencing. This approach presents several advantages over conventional padlock probes allowing simpler assay design for multiplexed mutation detection to screen for the presence of mutations in clinically relevant mutational hotspots directly in situ.

piggyBac-ing models and new therapeutic strategies

DNA transposons offer an efficient nonviral method of permanently modifying the genomes of mammalian cells. The piggyBac transposon system has proven effective in genomic engineering of mammalian cells for preclinical applications, including gene discovery, simultaneous multiplexed genome modification, animal transgenesis, gene transfer in vivo achieving long-term gene expression in animals, and the genetic modification of clinically relevant cell types, such as induced pluripotent stem cells and human T lymphocytes. piggyBac has many desirable features, including seamless excision of transposons from the genomic DNA and the potential to target integration events to desired DNA sequences. In this review, we explore these recent applications and also highlight the unique advantages of using piggyBac for developing new molecular therapeutic strategies.

Transcription factor regulation and chromosome dynamics during pseudohyphal growth

A multiplexed analysis of the transcriptional regulation of yeast pseudohyphal growth recorded the binding of 28 different transcription factors with barcoded transposons. A core set of target genes is identified, and a process of DNA looping at the FLO11 locus that provides transcriptional memory for expression of the gene is described.

Pseudohyphal growth is a developmental pathway seen in some strains of yeast in which cells form multicellular filaments in response to environmental stresses. We used multiplexed transposon “Calling Cards” to record the genome-wide binding patterns of 28 transcription factors (TFs) in nitrogen-starved yeast. We identified TF targets relevant for pseudohyphal growth, producing a detailed map of its regulatory network. Using tools from graph theory, we identified 14 TFs that lie at the center of this network, including Flo8, Mss11, and Mfg1, which bind as a complex. Surprisingly, the DNA-binding preferences for these key TFs were unknown. Using Calling Card data, we predicted the in vivo DNA-binding motif for the Flo8-Mss11-Mfg1 complex and validated it using a reporter assay. We found that this complex binds several important targets, including FLO11, at both their promoter and termination sequences. We demonstrated that this binding pattern is the result of DNA looping, which regulates the transcription of these targets and is stabilized by an interaction with the nuclear pore complex. This looping provides yeast cells with a transcriptional memory, enabling them more rapidly to execute the filamentous growth program when nitrogen starved if they had been previously exposed to this condition.

Targeting piggyBac transposon integrations in the human genome

DNA based transposon systems offer a technology for active and efficient delivery of genes into human cells. An emerging field is directed at manipulating such systems to achieve site-directed integration as compared to un-targeted integration which occurs with native or unmodified transposon systems. The naturally active piggyBac transposon system is derived from insects but has been shown to be very efficient in gene-modifying human cells. Recent efforts have utilized the fusion of DNA binding domains to the piggyBac transposase enzyme with the goal of targeting integration to specific locations in the human genome. In this chapter, we describe methodology for engineering and characterizing chimeric piggyBac transposase enzymes, including experimental approaches for evaluating activity and targeting specificity in the human genome.

A domesticated PiggyBac transposase interacts with heterochromatin and catalyzes reproducible DNA elimination in Tetrahymena

The somatic genome of the ciliated protist Tetrahymena undergoes DNA elimination of defined sequences called internal eliminated sequences (IESs), which account for ~30% of the germline genome. During DNA elimination, IES regions are heterochromatinized and assembled into heterochromatin bodies in the developing somatic nucleus. The domesticated piggyBac transposase Tpb2p is essential for the formation of heterochromatin bodies and DNA elimination. In this study, we demonstrate that the activities of Tpb2p involved in forming heterochromatin bodies and executing DNA elimination are genetically separable. The cysteine-rich domain of Tpb2p, which interacts with the heterochromatin-specific histone modifications, is necessary for both heterochromatin body formation and DNA elimination, whereas the endonuclease activity of Tpb2p is only necessary for DNA elimination. Furthermore, we demonstrate that the endonuclease activity of Tpb2p in vitro and the endonuclease activity that executes DNA elimination in vivo have similar substrate sequence preferences. These results strongly indicate that Tpb2p is the endonuclease that directly catalyzes the excision of IESs and that the boundaries of IESs are at least partially determined by the combination of Tpb2p-heterochromatin interaction and relaxed sequence preference of the endonuclease activity of Tpb2p.

piggyBac transposase tools for genome engineering

The transposon piggyBac is being used increasingly for genetic studies. Here, we describe modified versions of piggyBac transposase that have potentially wide-ranging applications, such as reversible transgenesis and modified targeting of insertions. piggyBac is distinguished by its ability to excise precisely, restoring the donor site to its pretransposon state. This characteristic makes piggyBac useful for reversible transgenesis, a potentially valuable feature when generating induced pluripotent stem cells without permanent alterations to genomic sequence. To avoid further genome modification following piggyBac excision by reintegration, we generated an excision competent/integration defective (Exc(+)Int(-)) transposase. Our findings also suggest the position of a target DNA-transposase interaction. Another goal of genome engineering is to develop reagents that can guide transgenes to preferred genomic regions. Others have shown that piggyBac transposase can be active when fused to a heterologous DNA-binding domain. An Exc(+)Int(-) transposase, the intrinsic targeting of which is defective, might also be a useful intermediate in generating a transposase whose integration activity could be rescued and redirected by fusion to a site-specific DNA-binding domain. We show that fusion to two designed zinc finger proteins rescued the Int(-) phenotype. Successful guided transgene integration into genomic DNA would have broad applications to gene therapy and molecular genetics. Thus, an Exc(+)Int(-) transposase is a potentially useful reagent for genome engineering and provides insight into the mechanism of transposase-target DNA interaction.

Mapping functional transcription factor networks from gene expression data

A critical step in understanding how a genome functions is determining which transcription factors (TFs) regulate each gene. Accordingly, extensive effort has been devoted to mapping TF networks. In Saccharomyces cerevisiae, protein–DNA interactions have been identified for most TFs by ChIP-chip, and expression profiling has been done on strains deleted for most TFs. These studies revealed that there is little overlap between the genes whose promoters are bound by a TF and those whose expression changes when the TF is deleted, leaving us without a definitive TF network for any eukaryote and without an efficient method for mapping functional TF networks. This paper describes NetProphet, a novel algorithm that improves the efficiency of network mapping from gene expression data. NetProphet exploits a fundamental observation about the nature of TF networks: The response to disrupting or overexpressing a TF is strongest on its direct targets and dissipates rapidly as it propagates through the network. Using S. cerevisiae data, we show that NetProphet can predict thousands of direct, functional regulatory interactions, using only gene expression data. The targets that NetProphet predicts for a TF are at least as likely to have sites matching the TF's binding specificity as the targets implicated by ChIP. Unlike most ChIP targets, the NetProphet targets also show evidence of functional regulation. This suggests a surprising conclusion: The best way to begin mapping direct, functional TF-promoter interactions may not be by measuring binding. We also show that NetProphet yields new insights into the functions of several yeast TFs, including a well-studied TF, Cbf1, and a completely unstudied TF, Eds1.

The piggyBac transposon displays local and distant reintegration preferences and can cause mutations at noncanonical integration sites

The DNA transposon piggyBac is widely used as a tool in mammalian experimental systems for transgenesis, mutagenesis, and genome engineering. We have characterized genome-wide insertion site preferences of piggyBac by sequencing a large set of integration sites arising from transposition from two separate genomic loci and a plasmid donor in mouse embryonic stem cells. We found that piggyBac preferentially integrates locally to the excision site when mobilized from a chromosomal location and identified other nonlocal regions of the genome with elevated insertion frequencies. piggyBac insertions were associated with expressed genes and markers of open chromatin structure and were excluded from heterochromatin. At the nucleotide level, piggyBac prefers to insert into TA-rich regions within a broader GC-rich context. We also found that piggyBac can insert into sites other than its known TTAA insertion site at a low frequency (2%). Such insertions introduce mismatches that are repaired with signatures of host cell repair pathways. Transposons could be mobilized from plasmids with the observed noncanonical flanking regions, indicating that piggyBac could generate point mutations in the genome.

Functional characterization of piggyBat from the bat Myotis lucifugus unveils an active mammalian DNA transposon

A revelation of the genomic age has been the contributions of the mobile DNA segments called transposable elements to chromosome structure, function, and evolution in virtually all organisms. Substantial fractions of vertebrate genomes derive from transposable elements, being dominated by retroelements that move via RNA intermediates. Although many of these elements have been inactivated by mutation, several active retroelements remain. Vertebrate genomes also contain substantial quantities and a high diversity of cut-and-paste DNA transposons, but no active representative of this class has been identified in mammals. Here we show that a cut-and-paste element called piggyBat, which has recently invaded the genome of the little brown bat (Myotis lucifugus) and is a member of the piggyBac superfamily, is active in its native form in transposition assays in bat and human cultured cells, as well as in the yeast Saccharomyces cerevisiae. Our study suggests that some DNA transposons are still actively shaping some mammalian genomes and reveals an unprecedented opportunity to study the mechanism, regulation, and genomic impact of cut-and-paste transposition in a natural mammalian host.