Extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation.

Higher-order chromosomal organization for transcription regulation is poorly understood in eukaryotes. Using genome-wide Chromatin Interaction Analysis with Paired-End-Tag sequencing (ChIA-PET), we mapped long-range chromatin interactions associated with RNA polymerase II in human cells and uncovered widespread promoter-centered intragenic, extragenic, and intergenic interactions. These interactions further aggregated into higher-order clusters, wherein proximal and distal genes were engaged through promoter-promoter interactions. Most genes with promoter-promoter interactions were active and transcribed cooperatively, and some interacting promoters could influence each other implying combinatorial complexity of transcriptional controls. Comparative analyses of different cell lines showed that cell-specific chromatin interactions could provide structural frameworks for cell-specific transcription, and suggested significant enrichment of enhancer-promoter interactions for cell-specific functions. Furthermore, genetically-identified disease-associated noncoding elements were found to be spatially engaged with corresponding genes through long-range interactions. Overall, our study provides insights into transcription regulation by three-dimensional chromatin interactions for both housekeeping and cell-specific genes in human cells.

MYC overexpression leads to increased chromatin interactions at super-enhancers and MYC binding sites.

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter-enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein-protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Core-binding factor fusion downregulation of ADAR2 RNA editing contributes to AML leukemogenesis.

Adenosine-to-inosine RNA editing, which is catalyzed by adenosine deaminases acting on RNA (ADAR) family of enzymes, ADAR1 and ADAR2, has been shown to contribute to multiple cancers. However, other than the chronic myeloid leukemia blast crisis, relatively little is known about its role in other types of hematological malignancies. Here, we found that ADAR2, but not ADAR1 and ADAR3, was specifically downregulated in the core-binding factor (CBF) acute myeloid leukemia (AML) with t(8;21) or inv(16) translocations. In t(8;21) AML, RUNX1-driven transcription of ADAR2 was repressed by the RUNX1-ETO additional exon 9a fusion protein in a dominant-negative manner. Further functional studies confirmed that ADAR2 could suppress leukemogenesis specifically in t(8;21) and inv16 AML cells dependent on its RNA editing capability. Expression of 2 exemplary ADAR2-regulated RNA editing targets coatomer subunit α and component of oligomeric Golgi complex 3 inhibits the clonogenic growth of human t(8;21) AML cells. Our findings support a hitherto, unappreciated mechanism leading to ADAR2 dysregulation in CBF AML and highlight the functional relevance of loss of ADAR2-mediated RNA editing to CBF AML.

MCIBox: a toolkit for single-molecule multi-way chromatin interaction visualization and micro-domains identification.

The emerging ligation-free three-dimensional (3D) genome mapping technologies can identify multiplex chromatin interactions with single-molecule precision. These technologies not only offer new insight into high-dimensional chromatin organization and gene regulation, but also introduce new challenges in data visualization and analysis. To overcome these challenges, we developed MCIBox, a toolkit for multi-way chromatin interaction (MCI) analysis, including a visualization tool and a platform for identifying micro-domains with clustered single-molecule chromatin complexes. MCIBox is based on various clustering algorithms integrated with dimensionality reduction methods that can display multiplex chromatin interactions at single-molecule level, allowing users to explore chromatin extrusion patterns and super-enhancers regulation modes in transcription, and to identify single-molecule chromatin complexes that are clustered into micro-domains. Furthermore, MCIBox incorporates a two-dimensional kernel density estimation algorithm to identify micro-domains boundaries automatically. These micro-domains were stratified with distinctive signatures of transcription activity and contained different cell-cycle-associated genes. Taken together, MCIBox represents an invaluable tool for the study of multiple chromatin interactions and inaugurates a previously unappreciated view of 3D genome structure.

Chromatin Interactions and Regulatory Elements in Cancer: From Bench to Bedside.

Chromatin interactions regulate gene expression by bringing distal regulatory elements, such as super-enhancers, to promoters in close spatial proximity. It has been recognized that in cancer, chromatin interactions can be dysregulated, leading to aberrant oncogene expression. Chromatin interactions may potentially serve as biomarkers, or be modulated via CRISPR therapy and small molecule inhibitors against transcription. However, these methods face challenges that must be resolved and raise questions for further research. Understanding chromatin interactions is essential for safety aspects of anticancer therapies, such as the mechanism of action of epigenetic regulators and transcription factors in cancer, and potential off-target effects arising from targeting super-enhancers and promoters. In this review article, we discuss how chromatin interactions and regulatory elements may become dysregulated in cancer, potential methods to target them for clinical therapy, and outline outstanding questions that require addressing before epigenetic therapies can translate to the clinic safely and effectively.

TP63, SOX2, and KLF5 Establish a Core Regulatory Circuitry That Controls Epigenetic and Transcription Patterns in Esophageal Squamous Cell Carcinoma Cell Lines.

BACKGROUND & AIMS: We investigated the transcriptome of esophageal squamous cell carcinoma (ESCC) cells, activity of gene regulatory (enhancer and promoter regions), and the effects of blocking epigenetic regulatory proteins. METHODS: We performed chromatin immunoprecipitation sequencing with antibodies against H3K4me1, H3K4me3, and H3K27ac and an assay for transposase-accessible chromatin to map the enhancer regions and accessible chromatin in 8 ESCC cell lines. We used the CRC_Mapper algorithm to identify core regulatory circuitry transcription factors in ESCC cell lines, and determined genome occupancy profiles for 3 of these factors. In ESCC cell lines, expression of transcription factors was knocked down with small hairpin RNAs, promoter and enhancer regions were disrupted by CRISPR/Cas9 genome editing, or bromodomains and extraterminal (BET) family proteins and histone deacetylases (HDACs) were inhibited with ARV-771 and romidepsin, respectively. ESCC cell lines were then analyzed by whole-transcriptome sequencing, immunoprecipitation, immunoblots, immunohistochemistry, and viability assays. Interactions between distal enhancers and promoters were identified and verified with circular chromosome conformation capture sequencing. NOD-SCID mice were given injections of modified ESCC cells, some mice where given injections of HDAC or BET inhibitors, and growth of xenograft tumors was measured. RESULTS: We identified super-enhancer-regulated circuits and transcription factors TP63, SOX2, and KLF5 as core regulatory factors in ESCC cells. Super-enhancer regulation of ALDH3A1 mediated by core regulatory factors was required for ESCC viability. We observed direct interactions between the promoter region of TP63 and functional enhancers, mediated by the core regulatory circuitry transcription factors. Deletion of enhancer regions from ESCC cells decreased expression of the core regulatory circuitry transcription factors and reduced cell viability; these same results were observed with knockdown of each core regulatory circuitry transcription factor. Incubation of ESCC cells with BET and HDAC disrupted the core regulatory circuitry program and the epigenetic modifications observed in these cells; mice given injections of HDAC or BET inhibitors developed smaller xenograft tumors from the ESCC cell lines. Xenograft tumors grew more slowly in mice given the combination of ARV-771 and romidepsin than mice given either agent alone. CONCLUSIONS: In epigenetic and transcriptional analyses of ESCC cell lines, we found the transcription factors TP63, SOX2, and KLF5 to be part of a core regulatory network that determines chromatin accessibility, epigenetic modifications, and gene expression patterns in these cells. A combination of epigenetic inhibitors slowed growth of xenograft tumors derived from ESCC cells in mice.

The enhancer RNA ARIEL activates the oncogenic transcriptional program in T-cell acute lymphoblastic leukemia.

The oncogenic transcription factor TAL1 regulates the transcriptional program in T-ALL. ARID5B is one of the critical downstream targets of TAL1, which further activates the oncogenic regulatory circuit in T-ALL cells. Here, we elucidated the molecular functions of the noncoding RNA, ARID5B-inducing enhancer associated long noncoding RNA (ARIEL), in T-ALL pathogenesis. We demonstrated that ARIEL is specifically activated in TAL1 + T-ALL cases, and its expression is associated with ARID5B enhancer activity. ARIEL recruits mediator proteins to the ARID5B enhancer, promotes enhancer-promoter interactions, and activates the expression of ARID5B, thereby positively regulating the TAL1-induced transcriptional program and the MYC oncogene. The TAL1 complex coordinately regulates the expression of ARIEL Knockdown of ARIEL inhibits cell growth and survival of T-ALL cells in culture and blocks disease progression in a murine xenograft model. Our results indicate that ARIEL plays an oncogenic role as an enhancer RNA in T-ALL.

Aberrant hyperediting of the myeloma transcriptome by ADAR1 confers oncogenicity and is a marker of poor prognosis.

DNA alterations have been extensively reported in multiple myeloma (MM); however, they cannot yet fully explain all the biological and molecular abnormalities in MM, which remains to this day an incurable disease with eventual emergence of refractory disease. Recent years have seen abnormalities at the RNA levels being reported to possess potential biological relevance in cancers. ADAR1-mediated A-to-I editing is an important posttranscriptional mechanism in human physiology, and the biological implication of its abnormality, especially at the global level, is underexplored in MM. In this study, we define the biological implications of A-to-I editing and how it contributes to MM pathogenesis. Here, we identified that the MM transcriptome is aberrantly hyperedited because of the overexpression of ADAR1. These events were associated with patients' survival independent of 1q21 amplifications and could affect patients' responsiveness to different treatment regimes. Our functional assays established ADAR1 to be oncogenic, driving cellular growth and proliferation in an editing-dependent manner. In addition, we identified NEIL1 (base-excision repair gene) as an essential and a ubiquitously edited ADAR1 target in MM. The recoded NEIL1 protein showed defective oxidative damage repair capacity and loss-of-function properties. Collectively, our data demonstrated that ADAR1-mediated A-to-I editing is both clinically and biologically relevant in MM. These data unraveled novel insights into MM molecular pathogenesis at the global RNA level.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

An oestrogen-receptor-alpha-bound human chromatin interactome.

Genomes are organized into high-level three-dimensional structures, and DNA elements separated by long genomic distances can in principle interact functionally. Many transcription factors bind to regulatory DNA elements distant from gene promoters. Although distal binding sites have been shown to regulate transcription by long-range chromatin interactions at a few loci, chromatin interactions and their impact on transcription regulation have not been investigated in a genome-wide manner. Here we describe the development of a new strategy, chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) for the de novo detection of global chromatin interactions, with which we have comprehensively mapped the chromatin interaction network bound by oestrogen receptor alpha (ER-alpha) in the human genome. We found that most high-confidence remote ER-alpha-binding sites are anchored at gene promoters through long-range chromatin interactions, suggesting that ER-alpha functions by extensive chromatin looping to bring genes together for coordinated transcriptional regulation. We propose that chromatin interactions constitute a primary mechanism for regulating transcription in mammalian genomes.

ChIA-PET analysis of transcriptional chromatin interactions.

Long-range chromatin contacts between specific DNA regulatory elements play a pivotal role in gene expression regulation, and a global characterization of these interactions in the 3-dimensional (3D) chromatin structure is imperative in understanding signaling networks and cell states. Chromatin Interaction Analysis using Paired-End Tag sequencing (ChIA-PET) is a method which converts functional chromatin structure into millions of short tag sequences. Combining Chromatin Immunoprecipitation (ChIP), proximity ligation and high-throughput sequencing, ChIA-PET provides a global and unbiased interrogation of higher-order chromatin structures associated with specific protein factors. Here, we describe the detailed procedures of the ChIA-PET methodology, unraveling transcription-associated chromatin contacts in a model human cell line.

Accurate Identification of DNA Replication Origin by Fusing Epigenomics and Chromatin Interaction Information.

DNA replication initiation is a complex process involving various genetic and epigenomic signatures. The correct identification of replication origins (ORIs) could provide important clues for the study of a variety of diseases caused by replication. Here, we design a computational approach named iORI-Epi to recognize ORIs by incorporating epigenome-based features, sequence-based features, and 3D genome-based features. The iORI-Epi displays excellent robustness and generalization ability on both training datasets and independent datasets of K562 cell line. Further experiments confirm that iORI-Epi is highly scalable in other cell lines (MCF7 and HCT116). We also analyze and clarify the regulatory role of epigenomic marks, DNA motifs, and chromatin interaction in DNA replication initiation of eukaryotic genomes. Finally, we discuss gene enrichment pathways from the perspective of ORIs in different replication timing states and heuristically dissect the effect of promoters on replication initiation. Our computational methodology is worth extending to ORI identification in other eukaryotic species.

MCI-frcnn: A deep learning method for topological micro-domain boundary detection.

Chromatin structural domains, or topologically associated domains (TADs), are a general organizing principle in chromatin biology. RNA polymerase II (RNAPII) mediates multiple chromatin interactive loops, tethering together as RNAPII-associated chromatin interaction domains (RAIDs) to offer a framework for gene regulation. RAID and TAD alterations have been found to be associated with diseases. They can be further dissected as micro-domains (micro-TADs and micro-RAIDs) by clustering single-molecule chromatin-interactive complexes from next-generation three-dimensional (3D) genome techniques, such as ChIA-Drop. Currently, there are few tools available for micro-domain boundary identification. In this work, we developed the MCI-frcnn deep learning method to train a Faster Region-based Convolutional Neural Network (Faster R-CNN) for micro-domain boundary detection. At the training phase in MCI-frcnn, 50 images of RAIDs from Drosophila RNAPII ChIA-Drop data, containing 261 micro-RAIDs with ground truth boundaries, were trained for 7 days. Using this well-trained MCI-frcnn, we detected micro-RAID boundaries for the input new images, with a fast speed (5.26 fps), high recognition accuracy (AUROC = 0.85, mAP = 0.69), and high boundary region quantification (genomic IoU = 76%). We further applied MCI-frcnn to detect human micro-TADs boundaries using human GM12878 SPRITE data and obtained a high region quantification score (mean gIoU = 85%). In all, the MCI-frcnn deep learning method which we developed in this work is a general tool for micro-domain boundary detection.

p53-NEIL1 co-abnormalities induce genomic instability and promote synthetic lethality with Chk1 inhibition in multiple myeloma having concomitant 17p13(del) and 1q21(gain).

Recurrent cytogenetic abnormalities are the main hallmark of multiple myeloma (MM) and patients having 2 or more high-risk prognostic events are associated with extremely poor outcome. 17p13(del) and 1q21(gain) are critical and independent high-risk cytogenetic markers, however, the biological significance underlying the poor outcome in MM patients having co-occurrence of both these chromosomal aberrations has never been interrogated. Herein, we identified that patients harbouring concomitant 17p13(del) with 1q21(gain) demonstrated the worst prognosis as compared to patients with single- (either 17p13(del) or 1q21(gain)) and with no chromosomal events (WT for both chromosomal loci); and they are highly enriched for genomic instability (GI) signature. We discovered that the GI feature in the patients with concomitant 17p13(del)-1q21(gain) was recapitulating the biological properties of myeloma cells with co-existing p53-deficiency and NEIL1 mRNA-hyper-editing (associated with chromosome 17p and 1q, respectively) that have inherent DNA damage response (DDR) and persistent activation of Chk1 pathway. Importantly, this became a vulnerable point for therapeutic targeting whereby the cells with this co-abnormalities demonstrated hyper-sensitivity to siRNA- and pharmacological-mediated-Chk1 inhibition, as observed at both the in vitro and in vivo levels. Mechanistically, this was attributable to the synthetic lethal relationship between p53-NEIL1-Chk1 abnormalities. The Chk1 inhibitor (AZD7762) tested showed good synergism with standard-of-care myeloma drugs, velcade and melphalan, thus further reinforcing the translational potential of this therapeutic approach. In summary, combination of NEIL1-p53 abnormalities with an ensuing Chk1 activation could serve as an Achilles heel and predispose MM cells with co-existing 1q21(gain) and 17p13(del) to therapeutic vulnerability for Chk1 inhibition.

Super Enhancer-Mediated Upregulation of HJURP Promotes Growth and Survival of t(4;14)-Positive Multiple Myeloma.

Multiple myeloma is an incurable malignancy with marked clinical and genetic heterogeneity. The cytogenetic abnormality t(4;14) (p16.3;q32.3) confers aggressive behavior in multiple myeloma. Recently, essential oncogenic drivers in a wide range of cancers have been shown to be controlled by super-enhancers (SE). We used chromatin immunoprecipitation sequencing of the active enhancer marker histone H3 lysine 27 acetylation (H3K27ac) to profile unique SEs in t(4;14)-translocated multiple myeloma. The histone chaperone HJURP was aberrantly overexpressed in t(4;14)-positive multiple myeloma due to transcriptional activation by a distal SE induced by the histone lysine methyltransferase NSD2. Silencing of HJURP with short hairpin RNA or CRISPR interference of SE function impaired cell viability and led to apoptosis. Conversely, HJURP overexpression promoted cell proliferation and abrogated apoptosis. Mechanistically, the NSD2/BRD4 complex positively coregulated HJURP transcription by binding the promoter and active elements of its SE. In summary, this study introduces SE profiling as an efficient approach to identify new targets and understand molecular pathogenesis in specific subtypes of cancer. Moreover, HJURP could be a valuable therapeutic target in patients with t(4;14)-positive myeloma. SIGNIFICANCE: A super-enhancer screen in t(4;14) multiple myeloma serves to identify genes that promote growth and survival of myeloma cells, which may be evaluated in future studies as therapeutic targets.

Chromatin Rewiring by Mismatch Repair Protein MSH2 Alters Cell Adhesion Pathways and Sensitivity to BET Inhibition in Gastric Cancer.

Mutations in the DNA mismatch repair gene MSH2 are causative of microsatellite instability (MSI) in multiple cancers. Here, we discovered that besides its well-established role in DNA repair, MSH2 exerts a novel epigenomic function in gastric cancer. Unbiased CRISPR-based mass spectrometry combined with genome-wide CRISPR functional screening revealed that in early-stage gastric cancer MSH2 genomic binding is not randomly distributed but rather is associated specifically with tumor-associated super-enhancers controlling the expression of cell adhesion genes. At these loci, MSH2 genomic binding was required for chromatin rewiring, de novo enhancer-promoter interactions, maintenance of histone acetylation levels, and regulation of cell adhesion pathway expression. The chromatin function of MSH2 was independent of its DNA repair catalytic activity but required MSH6, another DNA repair gene, and recruitment to gene loci by the SWI/SNF chromatin remodeler SMARCA4/BRG1. Loss of MSH2 in advanced gastric cancers was accompanied by deficient cell adhesion pathway expression, epithelial-mesenchymal transition, and enhanced tumorigenesis in vitro and in vivo. However, MSH2-deficient gastric cancers also displayed addiction to BAZ1B, a bromodomain-containing family member, and consequent synthetic lethality to bromodomain and extraterminal motif (BET) inhibition. Our results reveal a role for MSH2 in gastric cancer epigenomic regulation and identify BET inhibition as a potential therapy in MSH2-deficient gastric malignancies. SIGNIFICANCE: DNA repair protein MSH2 binds and regulates cell adhesion genes by enabling enhancer-promoter interactions, and loss of MSH2 causes deficient cell adhesion and bromodomain and extraterminal motif inhibitor synthetic lethality in gastric cancer.

Pan-cancer pervasive upregulation of 3' UTR splicing drives tumourigenesis.

Most mammalian genes generate messenger RNAs with variable untranslated regions (UTRs) that are important post-transcriptional regulators. In cancer, shortening at 3' UTR ends via alternative polyadenylation can activate oncogenes. However, internal 3' UTR splicing remains poorly understood as splicing studies have traditionally focused on protein-coding alterations. Here we systematically map the pan-cancer landscape of 3' UTR splicing and present this in SpUR ( http://www.cbrc.kaust.edu.sa/spur/home/ ). 3' UTR splicing is widespread, upregulated in cancers, correlated with poor prognosis and more prevalent in oncogenes. We show that antisense oligonucleotide-mediated inhibition of 3' UTR splicing efficiently reduces oncogene expression and impedes tumour progression. Notably, CTNNB1 3' UTR splicing is the most consistently dysregulated event across cancers. We validate its upregulation in hepatocellular carcinoma and colon adenocarcinoma, and show that the spliced 3' UTR variant is the predominant contributor to its oncogenic functions. Overall, our study highlights the importance of 3' UTR splicing in cancer and may launch new avenues for RNA-based anti-cancer therapeutics.

Chromatin interaction neural network (ChINN): a machine learning-based method for predicting chromatin interactions from DNA sequences.

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples.

Interplay and cooperation between SREBF1 and master transcription factors regulate lipid metabolism and tumor-promoting pathways in squamous cancer.

Squamous cell carcinomas (SCCs) comprise one of the most common histologic types of human cancer. Transcriptional dysregulation of SCC cells is orchestrated by tumor protein p63 (TP63), a master transcription factor (TF) and a well-researched SCC-specific oncogene. In the present study, both Gene Set Enrichment Analysis (GSEA) of SCC patient samples and in vitro loss-of-function assays establish fatty-acid metabolism as a key pathway downstream of TP63. Further studies identify sterol regulatory element binding transcription factor 1 (SREBF1) as a central mediator linking TP63 with fatty-acid metabolism, which regulates the biosynthesis of fatty-acids, sphingolipids (SL), and glycerophospholipids (GPL), as revealed by liquid chromatography tandem mass spectrometry (LC-MS/MS)-based lipidomics. Moreover, a feedback co-regulatory loop consisting of SREBF1/TP63/Kruppel like factor 5 (KLF5) is identified, which promotes overexpression of all three TFs in SCCs. Downstream of SREBF1, a non-canonical, SCC-specific function is elucidated: SREBF1 cooperates with TP63/KLF5 to regulate hundreds of cis-regulatory elements across the SCC epigenome, which converge on activating cancer-promoting pathways. Indeed, SREBF1 is essential for SCC viability and migration, and its overexpression is associated with poor survival in SCC patients. Taken together, these data shed light on mechanisms of transcriptional dysregulation in cancer, identify specific epigenetic regulators of lipid metabolism, and uncover SREBF1 as a potential therapeutic target and prognostic marker in SCC.

The use of multiple displacement amplification to amplify complex DNA libraries.

Complex libraries for genomic DNA and cDNA sequencing analyses are typically amplified using bacterial propagation. To reduce biases, large numbers of colonies are plated and scraped from solid-surface agar. This process is time consuming, tedious and limits scaling up. At the same time, multiple displacement amplification (MDA) has been recently developed as a method for in vitro amplification of DNA. However, MDA has no selection function for the removal of ligation multimers. We developed a novel method of briefly introducing ligation reactions into bacteria to select single insert DNA clones followed by MDA to amplify. We applied these methods to a Gene Identification Signatures with Paired-End diTags (GIS-PET) library, which is a complex transcriptome library created by pairing short tags from the 5' and 3' ends of cDNA fragments together, and demonstrated that this selection and amplification strategy is unbiased and efficient.