Hundreds Of Genomic Loci Research Articles

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm.

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Vitamin D: An Overview of Gene Regulation, Ranging from Metabolism to Genomic Effects.

Vitamin D is a pro-hormone characterized by an intricate metabolism and regulation. It is well known for its role in calcium and phosphate metabolism, and in bone health. However, several studies have assessed a huge number of extra-skeletal functions, ranging from cell proliferation in some oncogenic pathways to antioxidant and immunomodulatory functions. Vitamin D exerts its role by binding to VDRs (vitamin D receptors), which are located in many different tissues. Moreover, VDRs are able to bind hundreds of genomic loci, modulating the expression of various primary target genes. Interestingly, plenty of gene polymorphisms regarding VDRs are described, each one carrying a potential influence against gene expression, with relapses in several chronic diseases and metabolic complications. In this review, we provide an overview of the genetic aspects of vitamin D and VDR, emphasizing the gene regulation of vitamin D, and the genetic modulation of VDR target genes. In addition, we briefly summarize the rare genetic disease linked to vitamin D metabolism.

The DNA methylome of human vascular endothelium and its use in liquid biopsies

Vascular endothelial cells (VECs) are an essential component of each tissue, contribute to multiple pathologies, and are targeted by important drugs. Yet, there is a shortage of biomarkers to assess VEC turnover. To develop DNA methylation-based liquid biopsies for VECs, we determined the methylome of VECs isolated from freshly dissociated human tissues. A comparison with a human cell-type methylome atlas yielded thousands of loci that are uniquely unmethylated in VECs. These sites are typically gene enhancers, often residing adjacent to VEC-specific genes. We also identified hundreds of genomic loci that are differentially methylated in organotypic VECs, indicating that VECs feeding specific organs are distinct cell types with a stable epigenetic identity. We established universal and lung-specific VEC markers and evaluated their presence in circulating cell-free DNA (cfDNA). Nearly 2.5% of cfDNA in the plasma of healthy individuals originates from VECs. Sepsis, graft versus host disease, and cardiac catheterization are associated with elevated levels of VEC-derived cfDNA, indicative of vascular damage. Lung-specific VEC cfDNA is selectively elevated in patients with chronic obstructive pulmonary disease (COPD) or lung cancer, revealing tissue-specific vascular turnover. VEC cfDNA biomarkers inform vascular dynamics in health and disease, potentially contributing to early diagnosis and monitoring of pathologies, and assessment of drug activity. This work was supported by the Beutler Research Program, Helmsley Charitable Trust, JDRF, Grail and the DON Foundation (to Y.D.). Y.D holds the Walter& Greta Stiel Chair in heart studies. B.G., R.S., J.M., D.N., T.K., and Y.D. filed patents on cfDNA analysis.

Identifying susceptibility genes for essential hypertension by transcriptome-wide association study

Hypertension is a leading risk factor of cardiovascular disease and mortality in the population worldwide. Recently, hundreds of genomic loci were reported for hypertension by GWAS, however, the most SNPs are located in intergenic regions of genome, where a functional cause is difficult to determine. In the current study, a TWAS of hypertension was conducted using 452,264 individuals including 84,640 patients. KEGG and GO enrichment analyses were performed for the hypertension-related genes identified via TWAS. PPI network analysis based on the STRING database was also performed to detect TWAS-identified genes in hypertension. We have identified 18,420 genes from the GWAS summary data, and of those 1010 non-overlapping genes expression were significantly associated with hypertension after FDR correction (PFDR <0.05) in four tissues (left heart ventricle, aorta, whole blood, and peripheral blood). The KEGG and GO terms were mostly related to autoimmune mechanisms, and the autoimmune-related pathways have also been enriched using GO analysis for PPI genes. We further performed Mendelian randomization analysis, and the results supported a significant association between autoimmunity and hypertension. Moreover, 15 novel hypertension-susceptible genes were identified in all tissues, and five of the genes (RBM6, HLA-DRB5, UHRF1BP1, LYZ, and TMEM116) were associated with autoimmune system, which provide further evidence supporting an autoimmune mechanism in hypertension. In summary, our study supports that an autoimmune mechanism plays an important role in the development of hypertension, and these findings will provide new biological insights that will assist in deciphering the molecular etiology of hypertension.

Predicting causal genes from psychiatric genome-wide association studies using high-level etiological knowledge.

Genome-wide association studies have discovered hundreds of genomic loci associated with psychiatric traits, but the causal genes underlying these associations are often unclear, a research gap that has hindered clinical translation. Here, we present a Psychiatric Omnilocus Prioritization Score (PsyOPS) derived from just three binary features encapsulating high-level assumptions about psychiatric disease etiology - namely, that causal psychiatric disease genes are likely to be mutationally constrained, be specifically expressed in the brain, and overlap with known neurodevelopmental disease genes. To our knowledge, PsyOPS is the first method specifically tailored to prioritizing causal genes at psychiatric GWAS loci. We show that, despite its extreme simplicity, PsyOPS achieves state-of-the-art performance at this task, comparable to a prior domain-agnostic approach relying on tens of thousands of features. Genes prioritized by PsyOPS are substantially more likely than other genes at the same loci to have convergent evidence of direct regulation by the GWAS variant according to both DNA looping assays and expression or splicing quantitative trait locus (QTL) maps. We provide examples of genes hundreds of kilobases away from the lead variant, like GABBR1 for schizophrenia, that are prioritized by all three of PsyOPS, DNA looping and QTLs. Our results underscore the power of incorporating high-level knowledge of trait etiology into causal gene prediction at GWAS loci, and comprise a resource for researchers interested in experimentally characterizing psychiatric gene candidates.

High-Throughput DNA FISH (hiFISH).

High-throughput DNA fluorescence in situ hybridization (hiFISH) combines multicolor combinatorial DNA FISH staining with automated image acquisition and analysis to visualize and localize tens to hundreds of genomic loci in up to millions of cells. hiFISH can be used to measure physical distances between pairs of genomic loci, radial distances from genomic loci to the nuclear edge or center, and distances between genomic loci and nuclear structures defined by protein or RNA markers. The resulting large datasets of 3D spatial distances can be used to study cellular heterogeneity in genome architecture and the molecular mechanisms underlying this phenomenon in a variety of cellular systems. In this chapter we provide detailed protocols for hiFISH to measure distances between genomic loci, including all steps involved in DNA FISH probe design and preparation, cell culture, DNA FISH staining in 384-well imaging plates, automated image acquisition and analysis, and, finally, statistical analysis.

PRICKLE1 × FOCAD Interaction Revealed by Genome-Wide vQTL Analysis of Human Facial Traits.

The human face is a highly complex and variable structure resulting from the intricate coordination of numerous genetic and non-genetic factors. Hundreds of genomic loci impacting quantitative facial features have been identified. While these associations have been shown to influence morphology by altering the mean size and shape of facial measures, their effect on trait variance remains unclear. We conducted a genome-wide association analysis for the variance of 20 quantitative facial measurements in 2,447 European individuals and identified several suggestive variance quantitative trait loci (vQTLs). These vQTLs guided us to conduct an efficient search for gene-by-gene (G × G) interactions, which uncovered an interaction between PRICKLE1 and FOCAD affecting cranial base width. We replicated this G × G interaction signal at the locus level in an additional 5,128 Korean individuals. We used the hypomorphic Prickle1Beetlejuice (Prickle1Bj) mouse line to directly test the function of Prickle1 on the cranial base and observed wider cranial bases in Prickle1Bj/Bj. Importantly, we observed that the Prickle1 and Focadhesin proteins co-localize in murine cranial base chondrocytes, and this co-localization is abnormal in the Prickle1Bj/Bj mutants. Taken together, our findings uncovered a novel G × G interaction effect in humans with strong support from both epidemiological and molecular studies. These results highlight the potential of studying measures of phenotypic variability in gene mapping studies of facial morphology.

Functional Genomics Identify a Regulatory Risk Variation rs4420550 in the 16p11.2 Schizophrenia-Associated Locus

BackgroundGenome-wide association studies (GWASs) have reported hundreds of genomic loci associated with schizophrenia, yet identifying the functional risk variations is a key step in elucidating the underlying mechanisms. MethodsWe applied multiple bioinformatics and molecular approaches, including expression quantitative trait loci analyses, epigenome signature identification, luciferase reporter assay, chromatin conformation capture, homology-directed genome editing by CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/Cas9), RNA sequencing, and ATAC-Seq (assay for transposase-accessible chromatin using sequencing). ResultsWe found that the schizophrenia GWAS risk variations at 16p11.2 were significantly associated with messenger RNA levels of multiple genes in human brain, and one of the leading expression quantitative trait loci genes, MAPK3, is located ∼200 kb away from these risk variations in the genome. Further analyses based on the epigenome marks in human brain and cell lines suggested that a noncoding single nucleotide polymorphism, rs4420550 (p = 2.36 × 10−9 in schizophrenia GWAS), was within a DNA enhancer region, which was validated via in vitro luciferase reporter assays. The chromatin conformation capture experiment showed that the rs4420550 region physically interacted with the MAPK3 promoter and TAOK2 promoter. Precise CRISPR/Cas9 editing of a single base pair in cells followed by RNA sequencing further confirmed the regulatory effects of rs4420550 on the transcription of 16p11.2 genes, and ATAC-Seq demonstrated that rs4420550 affected chromatin accessibility at the 16p11.2 region. The rs4420550-[A/A] cells showed significantly higher proliferation rates compared with rs4420550-[G/G] cells. ConclusionsThese results together suggest that rs4420550 is a functional risk variation, and this study illustrates an example of comprehensive functional characterization of schizophrenia GWAS risk loci.

Genome-Scale Imaging of the 3D Organization and Transcriptional Activity of Chromatin

The 3D organization of chromatin regulates many genome functions. Our understanding of 3D genome organization requires tools to directly visualize chromatin conformation in its native context. Here we report an imaging technology for visualizing chromatin organization across multiple scales in single cells with high genomic throughput. First we demonstrate multiplexed imaging of hundreds of genomic loci by sequential hybridization, which allows high-resolution conformation tracing of whole chromosomes. Next we report a multiplexed error-robust fluorescence in situ hybridization (MERFISH)-based method for genome-scale chromatin tracing and demonstrate simultaneous imaging of more than 1,000 genomic loci and nascent transcripts of more than 1,000 genes together with landmark nuclear structures. Using this technology, we characterize chromatin domains, compartments, and trans-chromosomal interactions and their relationship to transcription in single cells. We envision broad application of this high-throughput, multi-scale, and multi-modal imaging technology, which provides an integrated view of chromatin organization in its native structural and functional context.

Complexities of Understanding Function from CKD-Associated DNA Variants.

Genome-wide association studies (GWASs) have facilitated the unbiased discovery of hundreds of genomic loci associated with CKD and kidney function. The vast majority of disease-associated DNA variants are noncoding. Those that are causal in CKD pathogenesis likely modulate transcription of target genes in a cell type-specific manner. To gain novel biological insights into mechanisms driving the development of CKD, the causal variants (which are usually not the most significant variant reported in a GWAS), their target genes, and causal cell types need to be identified. This functional validation requires a large number of new data sets, complex bioinformatics analyses, and experimental cellular and in vivo studies. Here, we review the basic principles and some of the current approaches being leveraged to assign functional significance to a genotype-phenotype association.

Mapping Alzheimer's Disease Variants to Their Target Genes Using Computational Analysis of Chromatin Configuration.

Genome-wide association studies (GWAS) have successfully identified hundreds of genomic loci that are associated with human traits and disease. However, because the majority of the genome-wide significant (GWS) loci fall onto the non-coding genome, the functional impact of many remain unknown. Three-dimensional chromatin interactions identified by Hi-C or its derivatives can provide useful tools to annotate these loci by linking non-coding variants to their actionable genes. Here, we outline a protocol to map GWAS non-coding variants to their putative genes using Alzheimer's disease (AD) GWAS and Hi-C datasets from human adult brain tissue. Putative causal single-nucleotide polymorphisms (SNPs) are identified by application of fine-mapping algorithms. SNPs are then mapped to their putative target genes using enhancer-promoter interactions based on Hi-C. The resulting gene set represents AD risk genes, as they are potentially regulated by AD risk variants. To garner further biological insights into molecular mechanisms underlying AD, we characterize AD risk genes using developmental brain expression data and brain single-cell expression profiles. This protocol can be expanded to any GWAS and Hi-C datasets to identify putative target genes and molecular mechanisms underlying various human traits and diseases.

Trouble With Tribbles-1.

Genome-wide association studies (GWAS) have identified hundreds of genomic loci in humans that are significantly associated with plasma cholesterol, triglycerides, and coronary artery disease. Although some loci contain genes with known regulatory roles in lipid metabolism and atherosclerosis, the majority were being implicated for the first time. The 8q24 locus, containing the gene TRIB1 ( Tribbles-1), is the only novel GWAS locus that associates with all 4 plasma lipid traits and coronary artery disease, an observation that has spurred immense interest in this locus. Subsequent in vivo loss and gain of function studies confirmed that Trib1 plays a role in hepatic lipid metabolism, validating the initial genetic observation. Yet, many challenges remain in discerning the nature of the association between the TRIB1 locus and cardiometabolic phenotypes. Is TRIB1 the causal gene at the 8q24 locus and what is the functional consequence of the associated noncoding variation? Is the relationship between TRIB1 and the transcription factor C/EBPα (CCAAT/enhancer-binding protein alpha) the primary molecular mechanism governing the genetic association or could it be an as yet unknown function for this interesting pseudokinase? Is hepatic TRIB1 the sole regulator of lipid metabolism or could extrahepatic TRIB1 play a role as well? The following review summarizes key findings related to these questions and highlights both the challenges and excitement in pursuing translational research of a novel gene in the post-GWAS era.

Common and Rare Genetic Risk Factors Converge in Protein Interaction Networks Underlying Schizophrenia.

Hundreds of genomic loci have been identified with the recent advances of schizophrenia in genome-wide association studies (GWAS) and sequencing studies. However, the functional interactions among those genes remain largely unknown. We developed a network-based approach to integrate multiple genetic risk factors, which lead to the discovery of new susceptibility genes and causal sub-networks, or pathways in schizophrenia. We identified significantly and consistently over-represented pathways in the largest schizophrenia GWA studies, which are highly relevant to synaptic plasticity, neural development and signaling transduction, such as long-term potentiation, neurotrophin signaling pathway, and the ERBB signaling pathway. We also demonstrated that genes targeted by common SNPs are more likely to interact with genes harboring de novo mutations (DNMs) in the protein-protein interaction (PPI) network, suggesting a mutual interplay of both common and rare variants in schizophrenia. We further developed an edge-based search algorithm to identify the top-ranked gene modules associated with schizophrenia risk. Our results suggest that the N-methyl-D-aspartate receptor (NMDAR) interactome may play a leading role in the pathology of schizophrenia, as it is highly targeted by multiple types of genetic risk factors.

Abstract 5331: Repression of TET genes and enhancement of DNA methyltransferase activity are critical for induction of aberrant DNA methylation

Abstract Chronic inflammation is deeply involved in the development of human cancers by inducing epigenetic alterations, such as aberrant DNA methylation. Among various inflammatory cytokines, the expression of Il1b, Tnf, and Nos2, is associated with induction of aberrant DNA methylation. However, the molecular mechanisms by which these cytokines induce aberrant DNA methylation remain unclear. Here, we show that the activation of the NF-kB signaling pathway, downstream of IL-1β and TNF-α, caused up-regulation of specific miRNAs, such as miR-20a, miR-26b, and miR-29c, and that these miRNAs caused repression of the Tet methylcytosine dioxygenases (Tet) genes, Tet1, Tet2, and Tet3. TET repression by overexpression of one of these miRNAs in cultured cells appeared to be insufficient for induction of aberrant DNA methylation as detected by an Infinium MethylationEPIC BeadChip array. However, triple knockout of TET genes induced strong DNA methylation at thousands of genomic loci. At the same time, exposure to nitric oxide, produced by Nos2, enhanced enzymatic activity of DNA methyltransferases (DNMTs). Treatment of cultured cells with nitric oxide induced weak DNA methylation at hundreds of genomic loci. These results show that both the repression of TET genes and enhancement of DNMT activity, induced by chronic inflammation, are critical for induction of aberrant DNA methylation. ''Vicious'' combination of such dysregulations might have a synergistic effect on induction of aberrant DNA methylation. Citation Format: Hideyuki Takeshima, Tohru Niwa, Harumi Yamada, Satoshi Yamashita, Mika Wakabayashi, Toshikazu Ushijima. Repression of TET genes and enhancement of DNA methyltransferase activity are critical for induction of aberrant DNA methylation [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5331.

H3.3K27M mutant proteins reprogram epigenome by sequestering the PRC2 complex to poised enhancers

Expression of histone H3.3K27M mutant proteins in human diffuse intrinsic pontine glioma (DIPG) results in a global reduction of tri-methylation of H3K27 (H3K27me3), and paradoxically, H3K27me3 peaks remain at hundreds of genomic loci, a dichotomous change that lacks mechanistic insights. Here, we show that the PRC2 complex is sequestered at poised enhancers, but not at active promoters with high levels of H3.3K27M proteins, thereby contributing to the global reduction of H3K27me3. Moreover, the levels of H3.3K27M proteins are low at the retained H3K27me3 peaks and consequently having minimal effects on the PRC2 activity at these loci. H3K27me3-mediated silencing at specific tumor suppressor genes, including Wilms Tumor 1, promotes proliferation of DIPG cells. These results support a model in which the PRC2 complex is redistributed to poised enhancers in H3.3K27M mutant cells and contributes to tumorigenesis in part by locally enhancing H3K27me3, and hence silencing of tumor suppressor genes.

GENE-10. H3K27ME3-MEDIATED SILENCING TUMOR SUPPRESSORS SUPPORTS THE PROLIFERATION OF PEDIATRIC BRAIN TUMOR CELLS HARBORING THE H3.3K27M MUTATION

The lysine 27 to methionine mutation of histone H3.3 (H3.3K27M) is detected in over 75% of diffuse intrinsic pontine glioma (DIPG). The H3.3K27M mutant proteins inhibit H3K27 methyltransferase complex PRC2, resulting in a global reduction of tri-methylation of H3K27 (H3K27me3). Paradoxically, high levels of H3K27me3 were also detected at hundreds of genomic loci. However, it is not known how and why H3K27me3 is redistributed in DIPG cells. Here we show that lower levels of H3.3K27M mutant proteins at some genomic loci contribute to the retention of H3K27me3 peaks. But more importantly, Jarid2, a PRC2-associated protein, strongly correlates the presence of H3K27me3 and relieves the H3.3K27M-mediated inhibition in vivo and in vitro. Furthermore, we show that H3K27me3-mediated silencing of tumor suppressor gene Wilms Tumor 1 (WT1) supports the proliferation of DIPG cells and reaction of WT1 inhibits DIPG proliferation. Together, these studies reveal mechanisms whereby H3K27me3 is retained in the environment of global loss of this mark, and how persistence of this mark contributes to DIPG tumorigenesis.

SpotLearn: Convolutional Neural Network for Detection of Fluorescence In Situ Hybridization (FISH) Signals in High-Throughput Imaging Approaches.

DNA fluorescence in situ hybridization (FISH) is the technique of choice to map the position of genomic loci in three-dimensional (3D) space at the single allele level in the cell nucleus. High-throughput DNA FISH methods have recently been developed using complex libraries of fluorescently labeled synthetic oligonucleotides and automated fluorescence microscopy, enabling large-scale interrogation of genomic organization. Although the FISH signals generated by high-throughput methods can, in principle, be analyzed by traditional spot-detection algorithms, these approaches require user intervention to optimize each interrogated genomic locus, making analysis of tens or hundreds of genomic loci in a single experiment prohibitive. We report here the design and testing of two separate machine learning-based workflows for FISH signal detection in a high-throughput format. The two methods rely on random forest (RF) classification or convolutional neural networks (CNNs), respectively. Both workflows detect DNA FISH signals with high accuracy in three separate fluorescence microscopy channels for tens of independent genomic loci, without the need for manual parameter value setting on a per locus basis. In particular, the CNN workflow, which we named SpotLearn, is highly efficient and accurate in the detection of DNA FISH signals with low signal-to-noise ratio (SNR). We suggest that SpotLearn will be useful to accurately and robustly detect diverse DNA FISH signals in a high-throughput fashion, enabling the visualization and positioning of hundreds of genomic loci in a single experiment.

Functional implications of disease-specific variants in loci jointly associated with coeliac disease and rheumatoid arthritis.

Hundreds of genomic loci have been associated with a significant number of immune-mediated diseases, and a large proportion of these associated loci are shared among traits. Both the molecular mechanisms by which these loci confer disease susceptibility and the extent to which shared loci are implicated in a common pathogenesis are unknown. We therefore sought to dissect the functional components at loci shared between two autoimmune diseases: coeliac disease (CeD) and rheumatoid arthritis (RA). We used a cohort of 12 381 CeD cases and 7827 controls, and another cohort of 13 819 RA cases and 12 897 controls, all genotyped with the Immunochip platform. In the joint analysis, we replicated 19 previously identified loci shared by CeD and RA and discovered five new non-HLA loci shared by CeD and RA. Our fine-mapping results indicate that in nine of 24 shared loci the associated variants are distinct in the two diseases. Using cell-type-specific histone markers, we observed that loci which pointed to the same variants in both diseases were enriched for marks of promoters active in CD14+ and CD34+ immune cells (P < 0.001), while loci pointing to distinct variants in one of the two diseases showed enrichment for marks of more specialized cell types, like CD4+ regulatory T cells in CeD (P < 0.0001) compared with Th17 and CD15+ in RA (P = 0.0029).

Identification of &lt;em&gt;Sleeping Beauty&lt;/em&gt; Transposon Insertions in Solid Tumors using Linker-mediated PCR

Genomic, proteomic, transcriptomic, and epigenomic analyses of human tumors indicate that there are thousands of anomalies within each cancer genome compared to matched normal tissue. Based on these analyses it is evident that there are many undiscovered genetic drivers of cancer1. Unfortunately these drivers are hidden within a much larger number of passenger anomalies in the genome that do not directly contribute to tumor formation. Another aspect of the cancer genome is that there is considerable genetic heterogeneity within similar tumor types. Each tumor can harbor different mutations that provide a selective advantage for tumor formation2. Performing an unbiased forward genetic screen in mice provides the tools to generate tumors and analyze their genetic composition, while reducing the background of passenger mutations. The Sleeping Beauty (SB) transposon system is one such method3. The SB system utilizes mobile vectors (transposons) that can be inserted throughout the genome by the transposase enzyme. Mutations are limited to a specific cell type through the use of a conditional transposase allele that is activated by Cre Recombinase. Many mouse lines exist that express Cre Recombinase in specific tissues. By crossing one of these lines to the conditional transposase allele (e.g. Lox-stop-Lox-SB11), the SB system is activated only in cells that express Cre Recombinase. The Cre Recombinase will excise a stop cassette that blocks expression of the transposase allele, thereby activating transposon mutagenesis within the designated cell type. An SB screen is initiated by breeding three strains of transgenic mice so that the experimental mice carry a conditional transposase allele, a concatamer of transposons, and a tissue-specific Cre Recombinase allele. These mice are allowed to age until tumors form and they become moribund. The mice are then necropsied and genomic DNA is isolated from the tumors. Next, the genomic DNA is subjected to linker-mediated-PCR (LM-PCR) that results in amplification of genomic loci containing an SB transposon. LM-PCR performed on a single tumor will result in hundreds of distinct amplicons representing the hundreds of genomic loci containing transposon insertions in a single tumor4. The transposon insertions in all tumors are analyzed and common insertion sites (CISs) are identified using an appropriate statistical method5. Genes within the CIS are highly likely to be oncogenes or tumor suppressor genes, and are considered candidate cancer genes. The advantages of using the SB system to identify candidate cancer genes are: 1) the transposon can easily be located in the genome because its sequence is known, 2) transposition can be directed to almost any cell type and 3) the transposon is capable of introducing both gain- and loss-of-function mutations6. The following protocol describes how to devise and execute a forward genetic screen using the SB transposon system to identify candidate cancer genes (Figure 1).

Identification of &lt;em&gt;Sleeping Beauty&lt;/em&gt; Transposon Insertions in Solid Tumors using Linker-mediated PCR

Genomic, proteomic, transcriptomic, and epigenomic analyses of human tumors indicate that there are thousands of anomalies within each cancer genome compared to matched normal tissue. Based on these analyses it is evident that there are many undiscovered genetic drivers of cancer(1). Unfortunately these drivers are hidden within a much larger number of passenger anomalies in the genome that do not directly contribute to tumor formation. Another aspect of the cancer genome is that there is considerable genetic heterogeneity within similar tumor types. Each tumor can harbor different mutations that provide a selective advantage for tumor formation(2). Performing an unbiased forward genetic screen in mice provides the tools to generate tumors and analyze their genetic composition, while reducing the background of passenger mutations. The Sleeping Beauty (SB) transposon system is one such method(3). The SB system utilizes mobile vectors (transposons) that can be inserted throughout the genome by the transposase enzyme. Mutations are limited to a specific cell type through the use of a conditional transposase allele that is activated by Cre Recombinase. Many mouse lines exist that express Cre Recombinase in specific tissues. By crossing one of these lines to the conditional transposase allele (e.g. Lox-stop-Lox-SB11), the SB system is activated only in cells that express Cre Recombinase. The Cre Recombinase will excise a stop cassette that blocks expression of the transposase allele, thereby activating transposon mutagenesis within the designated cell type. An SB screen is initiated by breeding three strains of transgenic mice so that the experimental mice carry a conditional transposase allele, a concatamer of transposons, and a tissue-specific Cre Recombinase allele. These mice are allowed to age until tumors form and they become moribund. The mice are then necropsied and genomic DNA is isolated from the tumors. Next, the genomic DNA is subjected to linker-mediated-PCR (LM-PCR) that results in amplification of genomic loci containing an SB transposon. LM-PCR performed on a single tumor will result in hundreds of distinct amplicons representing the hundreds of genomic loci containing transposon insertions in a single tumor(4). The transposon insertions in all tumors are analyzed and common insertion sites (CISs) are identified using an appropriate statistical method(5). Genes within the CIS are highly likely to be oncogenes or tumor suppressor genes, and are considered candidate cancer genes. The advantages of using the SB system to identify candidate cancer genes are: 1) the transposon can easily be located in the genome because its sequence is known, 2) transposition can be directed to almost any cell type and 3) the transposon is capable of introducing both gain- and loss-of-function mutations(6). The following protocol describes how to devise and execute a forward genetic screen using the SB transposon system to identify candidate cancer genes (Figure 1).