Open Access

PGBD5: a neural-specific intron-containing piggyBac transposase domesticated over 500 million years ago and conserved from cephalochordates to humans

  • Thomas Pavelitz1,
  • Lucas T Gray1,
  • Stephanie L Padilla1, 2,
  • Arnold D Bailey1 and
  • Alan M Weiner1Email author
Mobile DNA20134:23

DOI: 10.1186/1759-8753-4-23

Received: 3 June 2013

Accepted: 4 October 2013

Published: 1 November 2013



piggyBac domain (PGBD) transposons are found in organisms ranging from fungi to humans. Three domesticated piggyBac elements have been described. In the ciliates Paramecium tetraurelia and Tetrahymena thermophila, homologs known as piggyMacs excise internal eliminated sequences from germline micronuclear DNA during regeneration of the new somatic macronucleus. In primates, a PGBD3 element inserted into the Cockayne syndrome group B (CSB) gene over 43 Mya serves as an alternative 3′ terminal exon, enabling the CSB gene to generate both full length CSB and a conserved CSB-PGBD3 fusion protein that joins an N-terminal CSB domain to the C-terminal transposase domain.


We describe a fourth domesticated piggyBac element called PGBD5. We show that i) PGBD5 was first domesticated in the common ancestor of the cephalochordate Branchiostoma floridae (aka lancelet or amphioxus) and vertebrates, and is conserved in all vertebrates including lamprey but cannot be found in more basal urochordates, hemichordates, or echinoderms; ii) the lancelet, lamprey, and human PGBD5 genes are syntenic and orthologous; iii) no potentially mobile ancestral PGBD5 elements can be identified in other more deeply rooted organisms; iv) although derived from an IS4-related transposase of the RNase H clan, PGBD5 protein is unlikely to retain enzymatic activity because the catalytic DDD(D) motif is not conserved; v) PGBD5 is preferentially expressed in certain granule cell lineages of the brain and in the central nervous system based on available mouse and human in situ hybridization data, and the tissue-specificity of documented mammalian EST and mRNA clones; vi) the human PGBD5 promoter and gene region is rich in bound regulatory factors including the neuron-restrictive silencer factors NRSF/REST and CoREST, as well as SIN3, KAP1, STAT3, and CTCF; and vii) despite preferential localization within the nucleus, PGBD5 protein is unlikely to bind DNA or chromatin as neither DNase I digestion nor high salt extraction release PGBD5 from fractionated mouse brain nuclei.


We speculate that the neural-specific PGBD5 transposase was domesticated >500 My after cephalochordates and vertebrates split from urochordates, and that PGBD5 may have played a role in the evolution of a primitive deuterostome neural network into a centralized nervous system.


Branchiostoma Central nervous system Cephalochordate Domestication Granule cells PGBD5 PiggyBac


piggyBac family transposons have been identified in fungi, protozoa, cnidarians, plants, insects, crustaceans, echinoderms, urochordates (aka tunicates), hemichordates (acorn worm), fish, amphibia, and mammals suggesting both horizontal transmission and occasional domestication [1]. For example, the domesticated piggyMac transposases are catalytically active and required for programmed genome rearrangements in the ciliates Paramecium tetraurelia[2] and Tetrahymena thermophila[3]. However, piggyBacs are not the only active domesticated transposases: the RAG1/2 recombinases of the human immune system are presumably descended from an ancient Transib transposase [47] and, most recently, the human THAP9 gene has been found to encode a catalytically active P-element DNA transposase of as yet unknown function [8]. Some domesticated transposons such as the centomere protein CENPB retain nuclear localization and specific DNA binding but have lost enzymatic function [9, 10], whereas others such as SETMAR/Metnase [11, 12] and the Cockayne syndrome Group B-PiggyBac domain 3 (CSB-PGBD3) fusion protein [1315] retain site-specific DNA binding but gain new functions by fusion with upstream coding exons.

Most mammalian genomes contain only a handful of decayed piggyBac transposons. piggyBac transposons ceased activity over 35–40 Mya in the anthropoid primate lineage [16] and only somewhat more recently in the mouse lemur (prosimian) lineage [17]. The only exception known is the little brown bat, Myotis lucifugus, which contains thousands of active piggyBat elements [18, 19]. Humans have 5 substantially complete piggyBac elements, designated PGBD1, 2, 3, 4, and 5. PGBD1, 2, and 3 have multiple coding exons, but in each case the piggyBac transposase-related sequence is encoded by a single uninterrupted 3′ terminal exon. Thus, PGBD1 and 2 may resemble the PGBD3 transposon in which the transposase ORF is flanked upstream by a 3′ splice site and downstream by a polyadenylation site [13, 15]. As a result, insertion of PGBD3 into intron 5 of the CSB host gene enables the transposon to take advantage of the CSB promoter, using transposon-encoded alternative mRNA splicing and polyadenylation signals to express transposase as a CSB-PGBD3 fusion protein. The mouse genome lacks PGBD2, 3, and 4 homologs, but contains a modestly conserved PGBD1 suggesting introduction early in the mammalian radiation.

We now find that PGBD5, unlike other vertebrate piggyBacs, has been highly conserved in sequence and synteny from the primitive cephalochordate Branchiostoma floridae and the lamprey Petromyzon marinus (an agnathan or jawless fish) to humans. PGBD5, like other piggyBac elements, belongs to the RNase H clan of Pfam structures (, and exhibits a highly significant match (1.7e-80) to the Hidden Markov Model for eubacterial and archaeal IS4 transposases [20]. Whereas PGBD3 has sustained only a single D to N mutation in the essential catalytic triad DDD(D) [13] and retains the ability to bind the upstream piggyBac terminal inverted repeat [15], PGBD5 lacks 3 of the 4 conserved catalytic piggyBac aspartates [21] and does not appear to bind either DNA or chromatin; moreover, in contrast to all other piggyBac elements except the ciliate piggyMacs [2, 3], the PGBD5 transposase domain is encoded not by a single uninterrupted ORF but by 7 exons separated by long canonical introns most of which are conserved in position although intron loss and gain, or perhaps sliding, may also be involved. Taken together, the genomic data suggest that human PGBD5 is the most highly conserved piggyBac sequence known, and that it dates back over 500 My to the beginning of the chordate lineage.

An ancestral deuterostome is thought to have given rise to hemichordates and echinoderms through one line of descent and to chordates (including urochordates, cephalochordates, and vertebrates) through the other [22]. We show that PGBD5 is conserved from the primitive cephalochordate Branchiostoma floridae (aka lancelet or amphioxus) to humans, but appears to be absent in hemichordates, echinoderms, and urochordates. We also show that mouse PGBD5 is mainly nuclear and, consistent with publically available in situ hybridization and expression data, preferentially expressed in specific areas of the brain and central nervous system (CNS) that are enriched in granule cells. We also show that PGBD5 exon 1 is typically located far upstream from exons 2–7, and is embedded in a CpG island that is rich in bound neural transcriptional and regulatory factors. Finally, the PGBD5 transposase does not retain the catalytic DDD(D) motif found in active piggyBac elements, and is not released from nuclei by DNase I digestion or high salt extraction, suggesting that the transposase is not only inactive but fails to associate with either DNA or chromatin in vivo. Although many other scenarios can be imagined, we suggest, on the basis of these data, that PGBD5 may be expressed mainly, albeit not exclusively, in granule cells — an anatomically and functionally diverse population of small neurons, some of which are capable of adult neurogenesis [2325].

Results and discussion

Assembly of the human PGBD5 gene

An initial BLAST search of the human genome hg18 build in 2007 using the complete PGBD3 transposon from intron 5 of the human CSB gene [13] identified PGBD5 as a distant homolog of PGBD3. Unlike most other piggyBac elements which appear to be previously or currently mobile and to lack introns within the transposase domain, PGBD5 is a single copy gene with multiple introns within the transposase domain and is highly conserved among other vertebrates in intron/exon structure and protein sequence. Intriguingly, the annotated human PGBD5 open reading frame (ORF) in both the UCSC and RefSeq gene tracks of the UCSC Genome Browser continued through the 5′ end of the first annotated exon, which was flanked by canonical 3′ and 5′ splice sites. We therefore searched for additional exons upstream. A single expressed sequence tag (EST) clone CX757968 was annotated as joining a remote 5′ exon in frame to the downstream ORF. A 68 kb intron separating exon 1 from exons 2–7 seemed improbable, especially when represented at the time by only a single EST; however, as illustrated in Figure 1 and discussed in detail below, the evidence for this remote exon 1 quickly became convincing: i) CX757968 was annotated as joining exons 1–3, but upon sequencing we discovered that this EST is in fact a partial 3′ mRNA spanning PGBD5 exons 1–7 (see Methods for details); ii) splicing of CX757968 exon 1 to PGBD5 exons 2–7 generates a continuous ORF in which exon 1 encodes protein motifs that are moderately conserved in the orthologous exon 1 of other vertebrate PGBD5 genes, as well as the more ancient lamprey and lancelet/amphioxus genes; iii) exon 1 lies within the major CpG island located in or near each of the orthologous vertebrate PGBD5 genes; iv) the CpG island is rich in occupied transcription factor (TF) binding sites including the Neural-Restrictive Silencing Factor/RE1-Silencing TF (NRSF/REST); v) the great majority of available vertebrate PGBD5 mRNAs and ESTs are derived from neural or pluripotent tissues; vi) PGBD5 is almost exclusively expressed in the brain, with preferential expression in certain granule cell lineages; and vii) the CpG island and/or nearby sequences exhibit histone modification patterns associated with poised or weak promoters in 9 non-neural cell lines.
Figure 1

Clustal alignment of representative PGBD5 orthologs including Petromyzon marinus (sea lamprey) and Branchiostoma floridae (lancelet or amphioxus). For simplicity, many complete or partial vertebrate PGBD5s have been omitted from the alignment. The N-terminal motifs encoded by exon 1 are moderately conserved in all species including zebra finch, suggesting that chicken exon 1 lies within an unsequenced 770 bp gap located 29 kb upstream of exon 2 and immediately downstream of the sole substantial CpG island in the vicinity of the gene. Although human PGBD5 lacks 3 of the 4 catalytically active aspartates that are often conserved among diverse piggyBac elements, the positions of these four aspartates in a ClustalW alignment of piggyBac proteins most closely related to the active cabbage looper moth (Trichoplusia ni) transposase including human PGBD1, 2, 3, 4, and 5 [21] are highlighted in yellow. The Pfam homology with eubacterial and archeal IS4 transposases of the RNase H clan spans almost all of human PGBD5 exons 2–7 (residues 121–487). Figure key: black carets, vertebrate introns; blue caret, lamprey intron apparently orthologous to lancelet although shifted by 3 residues; red carets, lancelet introns; red underline, lancelet protein sequence derived from genomic tandem repeats (Additional file 1); red dashes, 13 residue deletion resulting from exclusion of predicted lancelet exon 5 which is embedded within the 108 bp genomic tandem repeats and would, if included, result in the 57 residue insertion; magenta underline, predicted nuclear localization signal not conserved in active Trichoplusia ni transposase; yellow highlight, position of four conserved, catalytic aspartates in active piggyBac transposases and homologs including human PGBD1, 2, 3, 4, and 5; gray XXX, regions of known length but undetermined sequence arbitrarily positioned in the clustal alignment; zfish, zfinch, xenopu, coelac, lampre, lancel are zebrafish, zebra finch, Xenopus tropicalis, coelacanth, lamprey, and lancelet respectively. Amino acid residues are colored according to the EBI Clustal convention for side chains (red, AVFPMILW; blue, DE; magenta, RK; green, STYHCNGQ; others, grey). To avoid prejudicial judgments regarding the relationship between highly divergent sequences, we refrained from assigning a similarity or homology score to each residue.

The search for more distantly related PGBD5 orthologs

As shown in Figure 1 (also see Methods and Additional file 1), we readily found orthologs of PGBD5 in the lamprey Petromyzon marinus, a jawless (agnathan) fish with a cartilaginous skeleton that is generally considered one of the most ancestral living vertebrates [26] and also in the lancelet Branchiostoma floridae (aka amphioxus), an even more primitive marine cephalochordate (Figure 2) with a notochord surmounted by a nerve chord leading to a single anterior eye.
Figure 2

Simplified phylogenetic tree of organisms examined. PGBD5 homologs are found in cephalochordates and all vertebrates examined, but nowhere else. This cladogram does not imply either the timing or degree of evolutionary divergence.

Saccoglossus kowalevichii (acorn worm) is a direct developing hemichordate that can provide important clues regarding the origin of chordates (Figure 2). An ancestral deuterostome is thought to have given rise to hemichordates and echinoderms through one line of descent, and to chordates, including the urochordates (aka tunicates), cephalochordates, and vertebrates, through another; moreover, a strong case can be made that similar genetic pathways are used to build the hemichordate and chordate body plans [22] and that hemichordates and chordates use similar genes to define anteroposterior neuroanatomy despite the evolutionary leap from a diffuse nerve network in hemichordates to the centralized chordate nervous system [27, 28].

A tBLASTn search of the NCBI genomic and EST databases for the acorn worm using the complete lancelet PGBD5 protein sequence as query generated only a few weak, fragmentary hits, although ESTs from three different developmental stages were included among the sequence libraries (mixed blastula and gastrula, mixed late gastrula and neurula, and early juveniles) [29]. In contrast, BLASTing the Skow_1.1 draft assembly of the S. kowalevichii genome ( revealed two groups of distant piggyBac homologs, 2 on contigs 47239 and 124986, and 3 on contigs 10046, 73983, and 91916; however, neither group was more closely related to human PGBD5 than PGBD3, and all 5 homologs had continuous intronless ORFs unlike PGBD5. Although there are precedents for functional, developmentally regulated retrogenes [30], exceedingly weak homology and the absence of introns suggest that acorn worm lacks a PGBD5 ortholog.

Similarly, only distant intronless piggyBac ORFs were found in three representative deuterostomes, the echinoderm Strongylocentrotus purpuratus (sea urchin), and the urochordates Ciona intestinalis (sea squirt) and Oikopleura dioica (an appendicularian). We note, however, that an argument can be made on the basis of whole genome comparisons that chordates are not necessarily monophyletic and that Oikopleura, instead of cephalochordates, may be the closest living relative of vertebrates [31]. No PGBD5 orthologs or homologs were found in the nematode C. elegans, the dipteran D. melanogaster, the cnidarian Hydra magnipapillata, the gastropod mollusc Aplysia californica (sea hare), or the poriferan Amphimedon queenslandica (marine sponge). We conclude that PGBD5 is likely to have originated in the chordate lineage and not in the ancestral deuterostome that gave rise to hemichordates and chordates (Figure 2).

Orthology of lancelet and lamprey PGBD5s with vertebrate PGBD5s

Orthologous genes from lancelet scaffolds tend to be concentrated in specific regions of human chromosomes, consistent with conservation of gene linkage between lancelet and vertebrates on a whole chromosome scale (macro-synteny) but limited conservation of local gene order (micro-synteny) [32]. A careful comparison of the micro and macro environment of the lancelet, lamprey, and human PGBD5 genes reveals substantial conservation at both levels, providing strong evidence that the lancelet, lamprey, and vertebrate PGBD5s are orthologous (Figure 3). Specifically, i) the synteny and orientation of 3 genes neighboring lancelet PGBD5 appears to be conserved in humans, and the change in gene order from EFCAB2-HNRNPU-PGBD5-EGLN1 in lancelet to EFCAB2-HNRNPU-EGLN1-PGBD5 in humans is consistent with both macro- and micro-synteny; ii) PGBD5 and COG2 are divergently transcribed in lamprey and human, but not closely linked in lancelet, suggesting that this gene order may have been established at the base of the vertebrate lineage; iii) GALNT2 and HNRNPU are syntenic in lamprey and human, consistent with macro-synteny; and iv) the lancelet, lamprey, and human PGBD5s each appear to be single copy, highly homologous to each other, but only distantly related to other piggyBac families and primate PGBD1, 2, 3, and 4 [1]. We conclude that lancelet, lamprey, and human PGBD5 are orthologs.
Figure 3

Lancelet, lamprey, and human PGBD5 are orthologous. Only shared syntenic genes are shown. The schematic is not drawn to scale, introns are not shown, genomic gaps are unsequenced, and distances are indicated in kb (red). The PGBD5 orthologs are oriented for clarity; the order and orientation of the two lamprey scaffolds is arbitrary. The 5′ and 3′ ends of lamprey GALNT2 are joined in the UCSC browser based on homology to the 5′ and 3′ ends of vertebrate GALNT2. Although >500 families of transposable elements constitute approximately 30% of the lancelet genome [32], no other PGBD5 homologs or fragments are found in the lancelet v2.0 draft genome; very similar PGBD5 sequences are present in both scaffolds 83 and 84 of the v1.0 draft genome (genome.​jgi-psf.​org/​Brafl1/​Brafl1.​home.​html), but PGBD5 appears only once in the v2.0 draft genome (downloadable from the UCSC browser at in a sequence context most closely resembling scaffold 84 (Additional file 1).

We were unable to trace the lancelet PGBD5 gene neighborhood further back in time to a more basal chordate. Homologs of EFCAB2, HNRNPU, EGLN1/2, and GALNT2 in the JGI v2 draft genome of the tunicate Ciona intestinalis are located on different scaffolds, and no COG2 homolog could be found by a BLAT search with human COG2 ( In another tunicate, the appendicularian Oikopleura dioica, only a GALNT2 homolog could be found among the existing scaffolds (http://​www.​genoscope.​cns.​fr/​externe/​GenomeBrowser/​Oikopleura/​). Nor could any PGBD5 homologs be found in either the Ciona or Oikopleura draft genomes. Thus, in the absence of flanking direct repeats or terminal inverted repeats in extant PGBD5s, we do not know whether the ancestral PGBD5 was a transposable element like most other piggyBacs, or even a cellular gene; whether and when it transposed or recombined into the neighborhood of EFCAB2, HNRNPU, and EGLN1/2; whether it had introns before arriving in the neighborhood, gained them once there, or has both gained and lost introns as demonstrated for chordate GIN1 and GIN2 integrase-related proteins [33]; whether it acquired a new, possibly neural-specific promoter and/or 5′ exon upon transposition or recombination, or continued to use an internal promoter as observed for the Trichoplusia ni piggyBac element in both insect [34] and mammalian cells [35]; and whether any ancestral PGBD5s have survived or are now all extinct.

Consistent with PGBD5 orthology, all vertebrates from lamprey to human share the same six introns (Figure 1, black carets). A seventh lamprey intron (blue caret) located in the middle of vertebrate exon 2 is shared with lancelet (red caret) and (despite a 3 amino acid residue shift between lancelet and lamprey) may be an orthologous intron that was lost in the last common ancestor of lamprey and true vertebrates. Although lancelet lacks vertebrate intron 1, it shares vertebrate introns 2 and 3; vertebrate introns 4, 5, and 6 are shifted by 11 to 42 amino acid residues relative to lancelet introns 4, 5, and 6, most probably reflecting a combination of intron loss and gain, rather than intron sliding which appears to be restricted to one or a few nucleotides [36]. We were unable to trace any PGBD5 sequences or intron positions further back than lancelet (Additional file 2).

We wish to stress that PGBD5 domestication (assuming PGBD5 descent from a transposable element) could have been a multistep process in which PGBD5 was initially selected or fixed for different reasons from those for which the gene is retained today. Indeed, PGBD5 could be retained for different reasons in lancelet and/or lamprey than in higher vertebrates which exhibit significantly greater PGBD5 sequence conservation (Figure 1, and Additional file 3). A similar scenario could hold for the SETMAR/Metnase fusion protein which evolved stepwise over a period of 18 My [11, 12]. Orthology can define a line of descent, but not necessarily related gene functions, regulation, or selection.

PGBD5 is primarily if not exclusively expressed in brain and CNS

The majority of available vertebrate PGBD5 mRNAs and ESTs are derived from embryonic, fetal, or adult brain (Additional file 4), consistent with Affymetrix expression profiles for 79 human tissues in the BioGPS gene portal ( [37] indicating that PGBD5 is strongly expressed in brain, very weakly expressed in spinal cord, and also weakly expressed in 4 different leukemias and lymphomas (K562 chronic myelogenous leukemia, HL60 promyelocytic leukemia, and both Daudi and Raji Burkitt’s lymphoma) but not in a variety of normal B- and T-cells (Figure 4).
Figure 4

PGBD5 is primarily if not exclusively expressed in brain and CNS. Affymetrix expression profiles for 79 human tissues and cell lines as described [37] on the BioGPS Gene portal ( indicate that PGBD5 is strongly expressed in brain, and very weakly expressed in spinal cord as well as 4 different leukemia and lymphoma cell lines but not in normal B- and T-cells. Vertical gray dashed lines, levels of expression; vertical magenta solid lines, mean, 4× mean, and 30× mean expression.

To localize PGBD5 expression more precisely within the brain and CNS, we examined public in situ hybridization databases including Max Planck’s Genepaint [38] and the Allen Brain Atlas [39]. Consistent with the BioGPS dataset, the PGBD5 transcript is expressed in restricted regions of the murine CNS. Over the course of development, PGBD5 is expressed in immature cells of the medial pallium and prepontine isthmus which give rise to the hippocampus and cerebellum, respectively (Figure 5, yellow arrowheads). In midline sagittal sections of the embryo, PGBD5 expression can be seen in aspects of the presumptive hypothalamus and medulla (Figure 5, grey arrowheads). In the adult mouse brain, PGBD5 is predominantly expressed in granule cells, including those of the olfactory bulb, hippocampus and cerebellum (Figure 5, white arrowheads). Although embryonic PGBD5 expression in the medial pallium and prepontine isthmus apparently persists in the adult hippocampus and cerebellum, other embryonic sites of PGBD5 expression are lost in the adult including the hypothalamus and medulla.
Figure 5

PGBD5 expression in the adult and developing mouse brain. Panels A-D illustrate expression of PGBD5 in sagittal brain sections (right) and a corresponding view from the Allen Institute Brain Atlas (left). The age of the mouse and the sagittal position relative to midline (0.00 mm) are indicated within each panel. The in situ hybridization stains are from public expression databases including The Allen Institute for Brain Science [39] for the adult mouse and Max Planck’s Genepaint [38] for embryos. The atlas images are color coded: olfactory bulb (green); cerebellum (orange); medial pallium (red); hypothalamus (brown); prepontine hindbrain (purple); and medullary hindbrain (blue). PGBD5 expression is restricted to a subset of cells within each nucleus denoted by a more saturated color. For example, PGBD5 expression in the cerebellar granule cell layer is indicated by a darker orange than the surrounding cerebellar cortex. The Max Planck and Allen Institute in situ hybridization probes both span exons 2–7, and are nearly identical; the absence of exon 1 is unlikely to affect the in situ patterns because there is no evidence in either the ENCODE (Additional file 5) or Chromatin State Segmentation data (Figure 6) of a functional promoter between exons 1 and 2.

Interestingly, PGBD5 is expressed in the granule layer of the olfactory bulb and dentate gyrus in the adult brain. Cells in both these regions are a mixture of mature and newborn granule neurons. In rodents, almost all neurogenesis occurs during embryogenesis and early perinatal development; however, the subventricular zone surrounding the lateral ventricles and the subgranular zone of the dentate gyrus are adult neurogenic niches, containing stem cells that give rise to adult-born granule cells in the olfactory bulb and dentate gyrus, respectively (for review see [24]). In the olfactory bulb, mature granule cells function as inhibitory interneurons that are thought to mediate contrast between odor stimuli [40] whereas adult-born granule cells in this region are induced by olfactory sensory inputs and may be involved in odor memory (for review see [41]). In the hippocampal dentate gyrus, mature granule cells mediate pattern completion during memory formation whereas adult-born granule cells are believed to facilitate pattern separation [42]. Although PGBD5 expression is seen in granule layers of the olfactory bulb and dentate gyrus by in situ hybridization, double labeling would be required to determine whether PGBD5 is expressed by mature or adult-born granule cells. Nevertheless, because PGBD5 continues to be expressed as the medial pallium develops into the mature hippocampus, we suspect that PGBD5-expressing granule neurons of the dentate gyrus are resident cells born during embryogenesis.

Using the Allen Institute’s microarray database, we also screened for PGBD5 expression in the adult human brain. As in murine brain, PGBD5 is expressed throughout the hippocampal formation, and in the pontine and aspects of the medulla including the raphe and arcuate. However, little or no PGBD5 expression was detected in human cerebellum by microarray analysis compared to mouse cerebellum by in situ hybridization.

The human PGBD5 promoter is embedded in a prominent CpG island and binds the neuron-restrictive silencer factors NRSF/REST and CoREST

Human PGBD5 exon 1 is phylogenetically conserved from branchiostomes to vertebrates (Figure 1), co-localizes with the single most prominent CpG island in the entire PGBD5 gene region, and is located 68 kb upstream of PGBD5 exons 2–7 encoding the piggyBac transposase domain (Figure 6 and Additional file 5). This CpG island is likely to harbor the PGBD5 promoter because i) the island overlaps the single strongest cluster of DNase I hypersensitive sites (ENCODE/Duke) that are diagnostic of open chromatin; ii) the island also overlaps the only strong chromatin marks in the Chromatin State Segmentation database (ENCODE/Broad) that are diagnostic of weak and/or poised promoters; and iii) the most highly occupied TF binding sites (ENCODE/Broad/HudsonAlpha/ Stanford/Duke/University of Washington) either overlap the island or cluster 24 kb upstream in the presumptive enhancer.
Figure 6

PGBD5 chromatin structure confirms location of promoter and transcription start site. (Upper panel) Structure of the human PGBD5 gene (top), chromatin state segmentation data from 9 non-neuronal cell lines (middle), and CpG islands (bottom). (Lower panel) An expanded view of the PGBD5 promoter and presumptive enhancer region. Chromatin state segmentation data are from the same 9 human cell types: GM12878 (EBV-transformed B-lymphocytes), H1-hESC (embryonic stem cells), K562 (chronic myelogenous leukemia), HepG2 (hepatocellular carcinoma), HMEC (mammary epithelial cells), HSMM (skeletal muscle myoblasts), HUVEC (umbilical vein endothelial cells), NHEK (epidermal keratinocytes), and NHLF (lung fibroblasts). The DNase I hypersensitive sites are from H1-hESC and GM12878. The TF ChIP-seq data are a composite of 24 neural and non-neural cell lines including H1-hESC, A549 (adenocarcinomic alveolar basal epithelial cells), BE2_C (brain neuroblastoma), HEK293 (embryonic kidney), HeLa (cervical carcinoma), HUVEC, Jurkat (T-lymphocyte), K562, HepG2, NB4 (acute promyelocytic leukemia), PANC-1 (pancreatic carcinoma), PFSK-1, SK-N-MC (neuronal epithelioma), SK-N-SH_RA (neuroblastoma differentiated with retinoic acid), and U87 (primary glioblastoma), and 10 EBV-transformed B-lymphocytes from various ethnic backgrounds. For a color key to the chromatin state segmentation data, see inset (lower panel). For CpG islands, darker green indicates higher CpG density. For DNase I hypersensitive sites, stronger signals indicate higher sensitivity. For TF ChIP-seq data, darker gray reflects higher occupancy. All data is publicly available as UCSC Genome Browser tracks: Open Chromatin by DNase I HS (ENCODE/Duke University); TF ChIP-seq (ENCODE/Broad Institute/HudsonAlpha Institute/Stanford/Duke/University of Washington); and chromatin state segmentation using a Hidden Markov Model to identify genome-wide patterns in ChIP-seq data for histone methylations and acetylations, histone variant H2AZ, RNAP II, and CTCF (ChromHMM from ENCODE/Broad Institute). Only 8 of the 15 distinguishable chromatin states [48, 49] are seen in this particular genomic interval of the 9 cell lines. Figure adapted from visualizations of current ENCODE databases available on the UCSC Genome Browser for hg19.

Consistent with BioGPS analysis (Figure 4), a majority of the 62 non-neural cell lines used for acquisition of the ENCODE TF ChIP-seq data (Figure 6 and Additional file 5) exhibit a single highly occupied site in the PGBD5 promoter for NRSF/REST, a factor which serves as a scaffold for assembly of many other proteins that can repress neural genes in non-neural cells [43]. For example, the N-terminal domain of NRSF/REST binds the SIN3A co-repressor, which in turn binds the repressive HDAC1/2 (histone deacetylases 1 and 2) and STAT3, thereby acting as a context-dependent ISGF3/STAT3 transcriptional switch [44]; the C-terminal domain of NRSF/REST binds the CoREST co-repressor which not only binds HDAC1/2, but also the repressive DNMT1 (DNA methyltransferase 1), histone K4 demethylase, and histone K9 methyltransferase, thereby regulating gene networks that control neural stem cell fate decisions [45]; and the internal domain of NRSF/REST binds a family of DNA motifs that regulate neuronal gene networks [46]. Thus, binding of NRSF/REST to the CpG island in non-neural cells is likely to explain not only colocalization of HDAC2 but also the Ini1, BRG1, and BAF155 components of the SWI/SNF chromatin remodeling complex required for repression of neuronal genes by the C-terminal CoREST corepressor [47].

Similarly, although none of the 9 cell lines used to build the Hidden Markov Model for chromatin state segmentation is neuronal [48, 49], all 9 lines display chromatin marks that are diagnostic of weak and poised promoters overlapping both the CpG island and exon 1 (Figure 6). Indeed, as might be expected, totipotent human H1 embryonic stem cells exhibit by far the most dramatic DNase I hypersensitivity indicative of open chromatin and the broadest chromatin state segmentation profile indicative of potential promoter activity.

NRSF/REST is unlikely to be the only repressor of PGBD5 in non-neural cells. Unexpectedly, loss of functional NRSF/REST during mouse and chicken development de-represses only a few neural genes in non-neural tissues, providing strong evidence for redundant repressors and/or the absence of neural activators in non-neural cells [50]. Consistent with this view, the ability of NRSF/REST to serve as a tumor suppressor in colorectal cancer [51] suggests that inappropriate expression of neural genes in non-neural tissues can have adverse consequences.

PGBD5 is nuclear in mouse brain but not bound to chromatin or DNA

Although a domesticated transposase might be expected to function within the nucleus, domestication and exaptation are full of surprises [52]. Vertebrate PGBD5s do indeed appear to have a strong candidate nuclear localization signal (NLS) with the sequence LRKRKKRKF, and the lamprey (ERRRRKKRKF) and lancelet (DRKRRKKKRF) sequences differ only slightly (Figure 1, magenta underline); however, this candidate NLS sequence is unique to the PGBD5 family and is located in the middle of the conserved transposase core, whereas the experimentally confirmed bipartite NLS in the cabbage looper moth transposase is C-terminally located [53]. Moreover, NLS sequences are typically diverse, whereas the candidate PGBD5 NLS is nearly invariant, suggesting that it may not be an NLS or is highly conserved for unrelated structural or functional reasons.

Given strong if not exclusive expression of PGBD5 in brain and CNS (Figures 4 and 5), as well as unexpected expression in B- and T-cell tumors (Figure 4), we initially decided to examine the subcellular localization of PGBD5 in a panel of pediatric brain tumor cell lines [54]. However, our recombinant PGBD5 (rPGBD5) antibody did not immunoprecipitate any protein of the expected size (about 59 kDa) although several prominent larger bands were visible on western blots (data not shown). In contrast, the rPGBD5 antibody cleanly immunoprecipitated a protein of about 59 kDa from the Daudi Burkitt’s lymphoma B-cell line and also recognized this protein on western blots along with an array of large prominent bands resembling those seen with the pediatric brain tumors (data not shown). These large proteins do not immunoprecipitate with the rPGBD5 antibody (Figures 7 and 8), and cross-react weakly only after denaturation and immobilization on the membrane.
Figure 7

PGBD5 partitions between nucleus and cytoplasm in the Daudi Burkitt’s lymphoma B-cell line. Nuclear and cytoplasmic fractions were prepared as described in Methods. (Upper panel) Whole cell, cytoplasmic, and nuclear fractions of 2×105 cells/lane were resolved by 7% SDS-PAGE, blotted, and probed with anti-rPGBD5 antibody. Mouse whole brain extract and human rPGBD5 served as controls. PGBD5 is indicated by an asterisk. Immunoreactive proteins larger than PGBD5 do not immunoprecipitate with anti-rPGBD5 antibody (compare with Figure 8). These proteins may cross-react weakly with antibody after denaturation and immobilization on the membrane, just as the antibody cross-reacts weakly with the hexahistidine tag but not the HA tag on western blotting (data not shown). (Lower panels) Daudi whole cell, cytoplasmic, and nuclear fractions were also probed for lamin A and β-actin to confirm successful partition.
Figure 8

Immunoprecipitation of endogenous PGBD5 from mouse brain. Mouse whole brain homogenate (5 or 25 μL) was immunoprecipitated with anti-rPGBD5 antibody. The proteins were resolved by SDS-PAGE, blotted, and probed with the same anti-rPGBD5 antibody as in Figure 7. rPGBD5, recombinant PGBD5; E1 and E2, first and second SDS elutions from the beads. Note that proteins larger than PGBD5 are present in whole brain homogenate (rightmost lane) but do not immunoprecipitate (lanes E1 and E2).

Although PGBD5 expression in leukemias and lymphomas (Figure 4) may be a fortuitous consequence of cell transformation, the convenience of working in an established cell line persuaded us to ask whether PGBD5 is nuclear or cytoplasmic in Daudi cells. Growing cells were washed and the plasma membrane disrupted by 0.1% NP40. Cytoplasmic and nuclear fractions were separated by centrifugation, the proteins resolved by SDS-PAGE, and PGBD5 assayed by western blotting with anti-rPGBD5 antibody (Figure 7). PGBD5 partitioned about equally to cytoplasmic and nuclear fractions whereas the large cross-reacting bands partitioned almost exclusively to the cytoplasm, suggesting that these may be abundant cytoskeletal proteins.

We ultimately decided to pursue PGBD5 subcellular localization in vivo using mouse whole brain homogenates because the mouse and human in situ hybridization data indicated expression of PGBD5 in only certain brain cell types (Figures 4 and 5). A BALBc mouse was euthanized, the brain including meninges removed and minced, and the tissue disrupted by Dounce homogenization with the loose and tight pestles in buffer containing protease inhibitors. The resulting homogenate was sonicated to open all cells, and then used without fractionation for immunoprecipitation and western blots. As expected, our polyclonal antibody against rPGBD5 cross-reacted with the nearly identical mouse PGBD5 (Figure 1), immunoprecipitating a single band from mouse brain that was close in size to rPGBD5 (Figure 8). As also observed for Daudi (Figure 7), mouse brain contains proteins that are larger than PGBD5 and cross-react on westerns but do not immunoprecipitate.

To date, all domesticated transposases of known or suspected function appear to be involved in nuclear DNA or chromatin transactions including SETMAR/Metnase [11, 55, 56], RAG1/2 [47], piggyMac [2, 3], piggyBat [18, 19], THAP9 [8], the TF ZBED6 [57], the major mammalian centromere protein CENP-B [9, 10], and the CSB-PGBD3 fusion protein [1315, 58]. We therefore asked whether mouse brain PGBD5 was also nuclear and associated with chromatin or DNA. A frozen whole mouse brain (PelFreez) was crushed in liquid nitrogen, and then homogenized with the loose pestle to disrupt cells while leaving nuclei intact as judged by DAPI staining. Cytoplasmic and nuclear fractions were separated by centrifugation and assayed by western blotting for PGBD5, using histone H3 and glutathione-S-transferase to confirm successful cell fractionation (Figure 9). PGBD5 partitioned about equally into the cytoplasmic and nuclear fractions, but the crude nuclear pellet was heavily contaminated with insoluble brain matter. To separate intact nuclei from lighter membranous material, the crude nuclear pellet was centrifuged through a 1.8-M sucrose cushion. PGBD5 again partitioned with histone H3 (Figure 9, rightmost lane) confirming nuclear localization.
Figure 9

PGBD5 partition in mouse brain. Whole frozen mouse brains were thawed, Dounced with a tight pestle, and partitioned as described in Methods. Crude nuclei were pelleted through a sucrose cushion to remove debris and residual cytoplasm. Aliquots corresponding to equal cell equivalents were assayed by western blotting for PGBD5, glutathione-S-transferase (cytoplasmic marker), and histone H3 (nuclear marker). Nuclear integrity was determined visually by DAPI staining. As in Figure 8, PGBD5 partitioned about equally between nucleus and cytoplasm.

Lastly, we asked whether nuclear PGBD5 could be released by increasing salt concentrations as is the case for non-histone chromosomal proteins and histones. Aliquots of the crude nuclear pellet were extracted with NaCl concentrations from 50 to 800 mM, with and without prior DNase I digestion (Figure 10A). After re-pelleting the nuclei, supernatant and pellet fractions were assayed for PGBD5 by SDS-PAGE and western blotting. PGBD5 was not released by complete DNase I digestion or by NaCl concentrations as high as 800 mM which efficiently extract core histone H3 (Figure 10B, compare right and left panels). These results suggest that PGBD5 is a structural component of the nucleus, perhaps bound to the nuclear matrix or lamina, but do not rule out a weaker or reversible association with chromatin or DNA.
Figure 10

Salt extraction of sonicated crude nuclei from mouse brain. (A) Crude nuclei from mouse whole brain were lysed and chromatin sheared by sonication. Nuclear lysis was assayed by DAPI staining before and after sonication (top left panel). Sonication was monitored by agarose gel electrophoresis with and without DNase I digestion (top right panel). (B) Sonicated nuclei were extracted with NaCl at the indicated concentrations (mM) with or without prior DNase I digestion. Supernatant and pellet fractions were separated by centrifugation, resolved by SDS-PAGE, and assayed by western blotting for PGBD5 (top panels) or histone H3 (bottom panels).


We do not know why the intron-containing PGBD5 element first appears in cephalochordates like the lancelet, which has only minimal sensory awareness of light, touch, vibration, and may sense pheromones, but lacks a neural crest, most placodes, and the organizer dividing midbrain and hindbrain [28]. Nor do we know why the ur-PGBD5 element was lost from (or failed to invade) more primitive urochordate and deuterostome lineages. We also cannot explain why a domesticated transposase that presumably once catalyzed DNA transactions would no longer bind DNA or chromatin yet still localize to the nucleus. Nonetheless, although the potential roles of PGBD5 in chordate brain development, neural cell differentiation, and evolution remain a mystery, it is unlikely to be a coincidence that PGBD5, arguably the most unusual and highly conserved member of the large piggyBac superfamily, first arose in the earliest basal chordate with a primitive CNS.

A major question regarding the origin of our own CNS is whether the chordate dorsal CNS arose from a dorsal centralized, ventral centralized, or uncentralized neural network in the ancestral deuterostome. Yet, the anteroposterior ectodermal expression of 38 orthologous genes involved in chordate CNS patterning is nearly identical in hemichordates and chordates, although hemichordates have a diffuse neural net and chordates a centralized CNS [27]. Admittedly, major advances in neuroanatomy and function are often based on existing genes, but it is possible that the seemingly de novo appearance of PGBD5 in a basal cephalochordate may have been a functional innovation that contributed to subsequent evolution of the rudimentary lancelet CNS [59].


Human PGBD5 cDNA reconstructed

The human CX757968 EST clone, generated from a pluripotent male blastocyst cell line, was obtained from the I.M.A.G.E. Consortium (now distributed by Thermo Scientific, Pittsburgh, PA, USA). This EST is annotated as spanning PGBD5 exons 1–3, but upon resequencing we found it to be a 3′ partial mRNA spanning the last 43 residues of predicted exon 1, all of exons 2–7, and a 1779 nt 3′ UTR between the translation termination codon and poly(A) tail (new GenBank Accession No. KF670820). The missing 67 N-terminal residues downstream of a phylogenetically conserved initiator methionine were embedded in an exceedingly GC-rich CpG island which resisted genomic PCR. Instead, we reconstructed the 5′ end of the PGBD5 ORF de novo by Klenow extension and Pfu polymerase amplification of overlapping synthetic oligodeoxynucleotides [60]. The final PGBD5 cassette was a 2636 bp AscI fragment containing a 1575 bp PGBD5 ORF followed by a 1061 bp 3′ untranslated region; the cassette was cloned into a modified pET45b vector (pET45b-HA-Asc) for bacterial protein production.

Search for additional PGBD5 orthologs or homologs

The complete human PGBD5 protein sequence encoded by exons 1–7 was used as a tBLASTn query to identify PGBD5 orthologs in the rat, mouse, zebra finch, chicken, frog, zebrafish, fugu, coelacanth, sea lamprey, and lancelet genomes. tBLASTn searches using the human, zebrafish, and lancelet PGBD5 proteins as queries were then carried out on the entire non-redundant nucleotide sequence collection (NCBI/NLM/NIH) as well as individual genome databases for 6 organisms: the sea squirt Ciona intestinalis, a urochordate (US Department of Energy, Joint Genome Institute); the sea urchin Strongylocentrotus purpuratus, an echinoderm (Human Genome Sequencing Center, Baylor College of Medicine); the acorn worm Saccoglossus kowalevichii, a hemichordate (Human Genome Sequencing Center, Baylor College of Medicine); Hydra magnipapillata, a Cnidarian (Center for Integrative Genomics, University of California at Berkeley); the sea hare Aplysia californica, a gastropod mollusk (Broad Institute through the UCSC Genome Browser); and the marine sponge Amphimedon queenslandica, a Poriferan (NCBI/NLM/NIH).

Assembly of the sea lamprey and coelacanth PGBD5 sequences

A Petromyzon marinus (sea lamprey) PGBD5 sequence was initially found using the human PGBD5 query for a tBLASTn search of the March 2007 draft genome assembly (WUSTL v.3.0) generated by the Washington University Genome Sequencing Center and currently available on the UCSC Genome Browser. Contig9988 (with PGBD5 on the + strand) and contig34746 (with PGBD5 on the – strand) overlap, and each contig fills in large unsequenced tracts in the other. We found no hint of the conserved exon 1 although contig34746 has 2,565 bp upstream of exon 2; however, exon 1 is 68 kb upstream of exon 2 in humans, and might also be far upstream in lamprey. Sea lamprey has recently been found to undergo extensive DNA diminution during development, significantly complicating genome analysis [6163]. The WUSTL v.3.0 lamprey assembly was subsequently recomputed and made available to us (JJ Smith, University of Kentucky, Lexington) in advance of publication [63]. The ortholog of human PGBD5 exons 2–7 was found in the 442 kb scaffold_316.1-442057 of the recomputed assembly, and annotated manually using the Ensembl predicted protein ENSPMAP00000004323 based on the WUSTL v.3.0 assembly. Although conserved exon 1 could not be found in the 204 kb of gapped sequence upstream of exon 2, it was ultimately located in a separate 4 kb scaffold_19254.1-4118. This small scaffold most likely maps within the large scaffold in 1 of 3 large unsequenced gaps (8, 29, and 63 kb) upstream of exon 2. PGBD5 exons 2–7 in the recomputed assembly differ only trivially from the Ensembl prediction based on the WUSTL assembly. The Latimeria chalumnae (coelacanth) ortholog of PGBD5 was kindly retrieved from contig_02766.1-262,491 of the draft genome assembly (CT Amemiya, Benaroya Research Institute and Department of Biology, University of Washington, Seattle) in advance of publication [64]. As was the case for lamprey, ORFs with homology to PGBD5 were identified and splice sites assigned manually. Coelacanth genome sequences have now been published by both the Amemiya [64] and Okada groups [65].

Reconstruction of the lancelet PGBD5 gene

The complete lancelet PGBD5 gene spanning lancelet exons 1–7 is annotated as Protein and Transcript ID 79338 in the Branchiostoma floridae v1.0 genome assembly and was retrieved from the – strand of scaffold 84 generated by the US Department of Energy Joint Genome Institute (JGI) site ( The resulting lancelet PGBD5 protein sequence is based on predicted splice sites that (with one exception described below) maximize similarity to vertebrate PGBD5 homologs; the sole independent support for the predicted protein sequence is cDNA clone bflv024o04 encoding the N-terminal 83–86 residues of the lancelet protein ( [66]. Curiously, scaffold 84 has 2 × 60 bp imperfect tandem repeats overlapping exon 3, and 3 × 108 bp imperfect tandem repeats overlapping exon 4 (Additional file 1) whereas the equivalent lancelet genomic sequence on the UCSC Genome Browser (chrUn: 288364515–288370215) has 7 instead of 2 × 60 bp repeats and 6 instead of 3 × 108 bp repeats (Additional file 1). More curiously, as discussed below, the first of 2 tandem repeats overlapping exon 3 in scaffold 84 immediately abuts the 5′ splice site of exon 3, and the first of the 3 tandem repeats overlapping exon 4 immediately abuts the 5′ splice site of intron 4 (Figure 1, red underlines). Thus, neither the 2 × 60 bp repeats nor the 3 × 108 bp repeats would substantially affect the homology of lancelet with the lamprey and vertebrate proteins, were it not for the fact that JGI Protein and Transcript ID 79338 annotates almost all of the second and third 108 bp repeats as exon 5, potentially introducing an additional 57 residues into the lancelet protein with no apparent homology to the other PGBD5 proteins (Figure 1). As JGI exon 5 relies on predicted 5′ and 3′ splice sites without confirmation by EST or mRNA data, and the second of the 2 × 60 bp repeats is not translated, we decided to exclude JGI exon 5 from the lancelet PGBD5 protein, instead jump splicing exon 4 directly to exon 6. The resulting 13 residue deletion (Figure 1, red dashes) seems somewhat more plausible, both structurally and evolutionarily, than a 57 residue insertion in the midst of a highly conserved region of the PGBD5 protein.

Our attempts to eliminate the apparent 13 residue deletion, to reduce the predicted 57 residue insertion, or otherwise increase homology of lancelet with other PGBD5 sequences, either by use of alternative 5′ and 3′ splice sites, or by internal deletion of one or more imperfect tandem repeats in various registers, were unsuccessful. Drs Jia-Xing Yue and Nik Putnam (Rice University) graciously reexamined the original Sanger sequencing data [32] as well as new Illumina data generated since the original publication, but the new data independently confirmed the original sequence. Given unforced homology of lancelet PGBD5 with the lamprey and vertebrate proteins except in the vicinity of intron 4 and predicted exon 5, the imperfect 3 × 108 bp tandem triplication could reflect actual duplications, genetic heterozygosity, or even sequence assembly errors, although the re-sequencing data argue strongly against this last possibility.

Lancelet PGBD5 genomic duplications include mRNA splice sites

The 60 and 108 bp tandem repeats that coincide with or overlap 5′ and 3′ splice sites echo earlier observations in lancelet: some lancelet globin genes have short tandem duplications of intron/exon boundaries (dubbed “mirages”) which apparently provide a mechanism for the creation of new introns within the genes [67]. Similarly, both the lancelet alcohol dehydrogenase (Adh) gene [68] and the troponin C (TnC) gene [69] contain “mirages” that duplicate splice sites yet do not generate aberrant mRNAs. The PGBD5 duplications are imperfect, unlike the globin, Adh, and TnC mirages, but all include intron/exon boundaries. These remarkable similarities suggest that tandem repeats in lancelet PGBD5 gene are real, but do not interfere with accurate mRNA splicing.

Protein production

The PGBD5 cassette in the pET45b-HA-Asc vector generates recombinant PGBD5 (rPGBD5) with N-terminal hexahistidine and hemagglutinin tags and a predicted molecular mass of 65.9 kDa, about 7 kDa larger than the predicted 58.5 kDa for endogenous human PGBD5 protein. The rPGBD5 expression construct was transformed into Rosetta cells (Novagen EMD), induced with IPTG, and shaken overnight at room temperature. Cells were opened by sonication and cleared supernatants subjected to consecutive ammonium sulfate cuts of 0–25%, 25–40%, and 40–55% saturation. rPGBD5 was located in the 25–40% pellet by Western blotting, dialyzed, and affinity-purified over a Talon (cobalt) column followed by desalting over G25 Sephadex and further purification by ion exchange chromatography on CM-Sepharose. rPGBD5 protein purity was assessed by SDS-PAGE followed by Coomassie Brilliant Blue or silver staining.

Antibody production

A rabbit polyclonal antibody against rPGBD5 was raised by R&R Research (Stanwood, WA, USA), affinity purified by binding to an NHS-Sepharose column coupled to rPGBD5, and eluted by low pH. The antibody immunoprecipitated human and mouse PGBD5 specifically, but cross-reacted with several abundant larger proteins after denaturation and immobilization on the polyvinylidene difluoride (PVDF) membrane for western blotting (Figures 7 and 8); the antibody also cross-reacted weakly with the hexahistidine tag but not with the HA tag on western blots (data not shown).

Cell culture and transfection

The human HT1080 fibrosarcoma cells were grown in MEM alpha with 5% FBS and 1× penicillin/streptomycin. The human Burkitt’s lymphoma B lymphoblast Daudi cell line was grown in RPMI 1640 with 10% FBS and 1× penicillin/streptomycin. The Trans IT-LT1 reagent was used for transfections as recommended (Mirus Bio, Madison, WI, USA).

Cell extracts

Daudi cells were washed once with PBS, resuspended in cold buffer D (20 mM Hepes pH 7.5, 50 mM NaCl, 2 mM MgCl2, 5% glycerol) and lysed by incubation with 0.1% NP-40 for 10 min on ice. Nuclei were separated from cytoplasm by centrifugation for 5 min at 13,000 rpm in a microfuge. An equal volume of 4× SDS-PAGE loading buffer was added and the sample sonicated to reduce viscosity.

Mouse brain extracts

Mouse whole brain extracts were initially prepared from fresh BALB/c brain by Dounce homogenization in buffer D containing 1 mM PMSF, 0.5X protease inhibitor cocktail (Roche) and sonicated to shear genomic DNA. Subsequently, cell fractions were prepared from flash frozen brains of 8–10 week old Swiss Webster mice (Pel-Freez Biologicals, Rodgers, AR, USA). The brain was thawed on ice, diced with a sterile scalpel, and homogenized by 25 strokes with the loose pestle of a 5 mL Dounce in 5 mL of buffer D. The resulting emulsion was sonicated to generate whole brain homogenates, or centrifuged for 20 min at 10,000 rpm in a Beckman JA-21 rotor to separate cytoplasm and nuclei. The crude nuclear pellet was washed with an equal volume of buffer D + 1 mM PMSF, respun, and then resuspended in 2.5 mL of buffer D + 10% glycerol + 1 mM PMSF for storage at −70°C. The fractionation and integrity of nuclei throughout the preparation were monitored by air drying an aliquot on a glass cover slip. The cover slips were immersed in PBS+DAPI, washed briefly in PBS, and inspected under the microscope for intact fluorescent nuclei. To determine whether PGBD5 binds to chromatin or DNA, crude nuclei were repeatedly sonicated in 10 sec bursts, assaying each time by DAPI staining until all nuclei were broken. Nuclear aliquots were treated with 0.1 U/μL freshly prepared DNAse I. The desired final concentration of NaCl was then added to DNase I treated and untreated nuclei, and the aliquots incubated for 15 min at room temperature before centrifugation to separate supernatant from the salt insoluble pellet. To prepare pure nuclei, crude nuclei were pelleted, resuspended in buffer D + 1 M sucrose, layered over a 1.8M sucrose cushion buffer D, and centrifuged for 1 h at 22,000 rpm (32,000 × g) in an Beckman Optima tabletop ultracentrifuge. The pelleted nuclei were resuspended in buffer D + 0.32 M sucrose and assayed by DAPI staining and western blotting. Mouse brain extracts were immunoprecipitated using anti-rPGBD5 antibody and the Thermo Pierce Direct IP kit for SDS PAGE and western blotting.



Central nervous system


Cockayne syndrome Group B


Expressed sequence tag


Open reading frame


Polyacrylamide gel electrophoresis


PiggyBac domain


Recombinant PGBD5


Transcription factor.



We thank Chris Amemiya (Benaroya Research Institute and Department of Biology, University of Washington, Seattle) and Jeramiah Smith (Department of Biology, University of Kentucky, Lexington) for providing unpublished Petromyzon marinus (sea lamprey) and Latimeria chalumnae (coelacanth) sequence data with homology to PGBD5. We also thank Jia-Xing Yue and Nik Putnam (Rice University, Houston) for reexamining their original Sanger sequencing data in light of newer Illumina data for the lancelet genome. We thank Richard Palmiter (HHMI and Department of Biochemistry, University of Washington, Seattle) for advice, and for fresh BALB/c mouse brain. We are deeply grateful to one of the referees for a trenchant critique stressing the value of synteny as evidence for orthology. This work was supported by NIH awards R01 GM41624 (AMW), the Cell and Molecular Biology Training Program T32 GM007270 (LTG), and a postdoctoral fellowship from the Hilda and Preston Davis Foundation (SLP). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Authors’ Affiliations

Department of Biochemistry, School of Medicine, University of Washington
Howard Hughes Medical Institute, University of Washington


  1. Sarkar A, Sim C, Hong YS, Hogan JR, Fraser MJ, Robertson HM, Collins FH: Molecular evolutionary analysis of the widespread piggyBac transposon family and related “domesticated” sequences. Mol Genet Genomics 2003, 270: 173-180. 10.1007/s00438-003-0909-0View ArticlePubMed
  2. Baudry C, Malinsky S, Restituito M, Kapusta A, Rosa S, Meyer E, Betermier M: PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements in the ciliate Paramecium tetraurelia. Genes Dev 2009, 23: 2478-2483. 10.1101/gad.547309PubMed CentralView ArticlePubMed
  3. Cheng CY, Vogt A, Mochizuki K, Yao MC: A domesticated piggyBac transposase plays key roles in heterochromatin dynamics and DNA cleavage during programmed DNA deletion in Tetrahymena thermophila. Mol Biol Cell 2010, 21: 1753-1762. 10.1091/mbc.E09-12-1079PubMed CentralView ArticlePubMed
  4. Fugmann SD, Messier C, Novack LA, Cameron RA, Rast JP: An ancient evolutionary origin of the Rag1/2 gene locus. Proc Natl Acad Sci USA 2006, 103: 3728-3733. 10.1073/pnas.0509720103PubMed CentralView ArticlePubMed
  5. Bernstein RM, Schluter SF, Bernstein H, Marchalonis JJ: Primordial emergence of the recombination activating gene 1 (RAG1): sequence of the complete shark gene indicates homology to microbial integrases. Proc Natl Acad Sci USA 1996, 93: 9454-9459. 10.1073/pnas.93.18.9454PubMed CentralView ArticlePubMed
  6. Kapitonov VV, Jurka J: RAG1 core and V(D)J recombination signal sequences were derived from Transib transposons. PLoS Biol 2005, 3: e181. 10.1371/journal.pbio.0030181PubMed CentralView ArticlePubMed
  7. Panchin Y, Moroz LL: Molluscan mobile elements similar to the vertebrate recombination-activating genes. Biochem Biophys Res Commun 2008, 369: 818-823. 10.1016/j.bbrc.2008.02.097PubMed CentralView ArticlePubMed
  8. Majumdar S, Singh A, Rio DC: The human THAP9 gene encodes an active P-element DNA transposase. Science 2013, 339: 446-448. 10.1126/science.1231789PubMed CentralView ArticlePubMed
  9. Smit AF, Riggs AD: Tiggers and DNA transposon fossils in the human genome. Proc Natl Acad Sci USA 1996, 93: 1443-1448. 10.1073/pnas.93.4.1443PubMed CentralView ArticlePubMed
  10. Casola C, Hucks D, Feschotte C: Convergent domestication of pogo-like transposases into centromere-binding proteins in fission yeast and mammals. Mol Biol Evol 2008, 25: 29-41.PubMed CentralView ArticlePubMed
  11. Cordaux R, Udit S, Batzer MA, Feschotte C: Birth of a chimeric primate gene by capture of the transposase gene from a mobile element. Proc Natl Acad Sci USA 2006, 103: 8101-8106. 10.1073/pnas.0601161103PubMed CentralView ArticlePubMed
  12. Hromas R, Williamson EA, Fnu S, Lee YJ, Park SJ, Beck BD, You JS, Laitao A, Nickoloff JA, Lee SH: Chk1 phosphorylation of Metnase enhances DNA repair but inhibits replication fork restart. Oncogene 2012, 31: 4245-4254. 10.1038/onc.2011.586PubMed CentralView ArticlePubMed
  13. Newman JC, Bailey AD, Fan HY, Pavelitz T, Weiner AM: An abundant evolutionarily conserved CSB-PiggyBac fusion protein expressed in Cockayne syndrome. PLoS Genet 2008, 4: e1000031. 10.1371/journal.pgen.1000031PubMed CentralView ArticlePubMed
  14. Bailey AD, Gray LT, Pavelitz T, Newman JC, Horibata K, Tanaka K, Weiner AM: The conserved Cockayne syndrome B-piggyBac fusion protein (CSB-PGBD3) affects DNA repair and induces both interferon-like and innate antiviral responses in CSB-null cells. DNA Repair (Amst) 2012, 11: 488-501. 10.1016/j.dnarep.2012.02.004View Article
  15. Gray LT, Fong KK, Pavelitz T, Weiner AM: Tethering of the conserved piggyBac transposase fusion protein CSB-PGBD3 to chromosomal AP-1 proteins regulates expression of nearby genes in humans. PLoS Genet 2012, 8: e1002972. 10.1371/journal.pgen.1002972PubMed CentralView ArticlePubMed
  16. Pace JK 2nd, Feschotte C: The evolutionary history of human DNA transposons: evidence for intense activity in the primate lineage. Genome Res 2007, 17: 422-432. 10.1101/gr.5826307PubMed CentralView ArticlePubMed
  17. Pagan HJ, Smith JD, Hubley RM, Ray DA: PiggyBac-ing on a primate genome: novel elements, recent activity and horizontal transfer. Genome Biol Evol 2010, 2: 293-303. 10.1093/gbe/evq021PubMed CentralView ArticlePubMed
  18. Ray DA, Feschotte C, Pagan HJ, Smith JD, Pritham EJ, Arensburger P, Atkinson PW, Craig NL: Multiple waves of recent DNA transposon activity in the bat, Myotis lucifugus. Genome Res 2008, 18: 717-728. 10.1101/gr.071886.107PubMed CentralView ArticlePubMed
  19. Mitra R, Li X, Kapusta A, Mayhew D, Mitra RD, Feschotte C, Craig NL: Functional characterization of piggyBat from the bat Myotis lucifugus unveils an active mammalian DNA transposon. Proc Natl Acad Sci USA 2013, 110: 234-239. 10.1073/pnas.1217548110PubMed CentralView ArticlePubMed
  20. De Palmenaer D, Siguier P, Mahillon J: IS4 family goes genomic. BMC Evol Biol 2008, 8: 18. 10.1186/1471-2148-8-18PubMed CentralView ArticlePubMed
  21. Keith JH, Schaeper CA, Fraser TS, Fraser MJ Jr: Mutational analysis of highly conserved aspartate residues essential to the catalytic core of the piggyBac transposase. BMC Mol Biol 2008, 9: 73. 10.1186/1471-2199-9-73PubMed CentralView ArticlePubMed
  22. Swalla BJ: Building divergent body plans with similar genetic pathways. Heredity (Edinb) 2006, 97: 235-243. 10.1038/sj.hdy.6800872View Article
  23. Coskun V, Zhao J, Sun YE: Neurons or glia? Can SHP2 know it all? Sci STKE 2007, 2007: pe58.PubMed CentralView ArticlePubMed
  24. Zhao C, Deng W, Gage FH: Mechanisms and functional implications of adult neurogenesis. Cell 2008, 132: 645-660. 10.1016/j.cell.2008.01.033View ArticlePubMed
  25. Islam O, Gong X, Rose-John S, Heese K: Interleukin-6 and neural stem cells: more than gliogenesis. Mol Biol Cell 2009, 20: 188-199. 10.1091/mbc.E08-05-0463PubMed CentralView ArticlePubMed
  26. Wullimann MF: Basal ganglia: insights into origins from lamprey brains. Curr Biol 2011, 21: R497-R500. 10.1016/j.cub.2011.05.052View ArticlePubMed
  27. Gerhart J, Lowe C, Kirschner M: Hemichordates and the origin of chordates. Curr Opin Genet Dev 2005, 15: 461-467. 10.1016/j.gde.2005.06.004View ArticlePubMed
  28. Holland LZ: Chordate roots of the vertebrate nervous system: expanding the molecular toolkit. Nat Rev Neurosci 2009, 10: 736-746. 10.1038/nrn2703View ArticlePubMed
  29. Freeman RM Jr, Wu M, Cordonnier-Pratt MM, Pratt LH, Gruber CE, Smith M, Lander ES, Stange-Thomann N, Lowe CJ, Gerhart J, Kirschner M: cDNA sequences for transcription factors and signaling proteins of the hemichordate Saccoglossus kowalevskii: efficacy of the expressed sequence tag (EST) approach for evolutionary and developmental studies of a new organism. Biol Bull 2008, 214: 284-302. 10.2307/25470670View ArticlePubMed
  30. Chen M, Zou M, Fu B, Li X, Vibranovski MD, Gan X, Wang D, Wang W, Long M, He S: Evolutionary patterns of RNA-based duplication in non-mammalian chordates. PLoS One 2011, 6: e21466. 10.1371/journal.pone.0021466PubMed CentralView ArticlePubMed
  31. Delsuc F, Brinkmann H, Chourrout D, Philippe H: Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature 2006, 439: 965-968. 10.1038/nature04336View ArticlePubMed
  32. Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, Benito-Gutiérrez EL, Dubchak I, Garcia-Fernàndez J, Gibson-Brown JJ, Grigoriev IV, Horton AC, de Jong PJ, Jurka J, Kapitonov VV, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin-I T, et al.: The amphioxus genome and the evolution of the chordate karyotype. Nature 2008, 453: 1064-1071. 10.1038/nature06967View ArticlePubMed
  33. Marin I: GIN transposons: genetic elements linking retrotransposons and genes. Mol Biol Evol 2010, 27: 1903-1911. 10.1093/molbev/msq072View ArticlePubMed
  34. Handler AM, McCombs SD, Fraser MJ, Saul SH: The lepidopteran transposon vector, piggyBac, mediates germ-line transformation in the Mediterranean fruit fly. Proc Natl Acad Sci USA 1998, 95: 7520-7525. 10.1073/pnas.95.13.7520PubMed CentralView ArticlePubMed
  35. Cadinanos J, Bradley A: Generation of an inducible and optimized piggyBac transposon system. Nucleic Acids Res 2007, 35: e87. 10.1093/nar/gkm446PubMed CentralView ArticlePubMed
  36. Rogozin IB, Lyons-Weiler J, Koonin EV: Intron sliding in conserved gene families. Trends Genet 2000, 16: 430-432. 10.1016/S0168-9525(00)02096-5View ArticlePubMed
  37. Wu C, Macleod I, Su AI: BioGPS and organizing online, gene-centric information. Nucleic Acids Res 2013, 41: D561-D565. 10.1093/nar/gks1114PubMed CentralView ArticlePubMed
  38. Visel A, Thaller C, Eichele G: an atlas of gene expression patterns in the mouse embryo. Nucleic Acids Res 2004, 32: D552-D556. 10.1093/nar/gkh029PubMed CentralView ArticlePubMed
  39. Lein ES, Hawrylycz MJ, Ao N, Ayres M, Bensinger A, Bernard A, Boe AF, Boguski MS, Brockway KS, Byrnes EJ, Chen L, Chen L, Chen TM, Chin MC, Chong J, Crook BE, Czaplinska A, Dang CN, Datta S, Dee NR, Desaki AL, Desta T, Diep E, Dolbeare TA, Donelan MJ, Dong HW, Dougherty JG, Duncan BJ, Ebbert AJ, Eichele G, et al.: Genome-wide atlas of gene expression in the adult mouse brain. Nature 2007, 445: 168-176. 10.1038/nature05453View ArticlePubMed
  40. Yokoi M, Mori K, Nakanishi S: Refinement of odor molecule tuning by dendrodendritic synaptic inhibition in the olfactory bulb. Proc Natl Acad Sci USA 1995, 92: 3371-3375. 10.1073/pnas.92.8.3371PubMed CentralView ArticlePubMed
  41. Lledo PM, Saghatelyan A: Integrating new neurons into the adult olfactory bulb: joining the network, life-death decisions, and the effects of sensory experience. Trends Neurosci 2005, 28: 248-254. 10.1016/j.tins.2005.03.005View ArticlePubMed
  42. Nakashiba T, Cushman JD, Pelkey KA, Renaudineau S, Buhl DL, McHugh TJ, Rodriguez Barrera V, Chittajallu R, Iwamoto KS, McBain CJ, Fanselow MS, Tonegawa S: Young dentate granule cells mediate pattern separation, whereas old granule cells facilitate pattern completion. Cell 2012, 149: 188-201. 10.1016/j.cell.2012.01.046PubMed CentralView ArticlePubMed
  43. Ballas N, Mandel G: The many faces of REST oversee epigenetic programming of neuronal genes. Curr Opin Neurobiol 2005, 15: 500-506. 10.1016/j.conb.2005.08.015View ArticlePubMed
  44. Icardi L, Mori R, Gesellchen V, Eyckerman S, De Cauwer L, Verhelst J, Vercauteren K, Saelens X, Meuleman P, Leroux-Roels G, De Bosscher K, Boutros M, Tavernier J: The Sin3a repressor complex is a master regulator of STAT transcriptional activity. Proc Natl Acad Sci USA 2012, 109: 12058-12063. 10.1073/pnas.1206458109PubMed CentralView ArticlePubMed
  45. Abrajano JJ, Qureshi IA, Gokhan S, Molero AE, Zheng D, Bergman A, Mehler MF: Corepressor for element-1-silencing transcription factor preferentially mediates gene networks underlying neural stem cell fate decisions. Proc Natl Acad Sci USA 2010, 107: 16685-16690. 10.1073/pnas.0906917107PubMed CentralView ArticlePubMed
  46. Otto SJ, McCorkle SR, Hover J, Conaco C, Han JJ, Impey S, Yochum GS, Dunn JJ, Goodman RH, Mandel G: A new binding motif for the transcriptional repressor REST uncovers large gene networks devoted to neuronal functions. J Neurosci 2007, 27: 6729-6739. 10.1523/JNEUROSCI.0091-07.2007View ArticlePubMed
  47. Battaglioli E, Andres ME, Rose DW, Chenoweth JG, Rosenfeld MG, Anderson ME, Mandel G: REST repression of neuronal genes requires components of the hSWI.SNF complex. J Biol Chem 2002, 277: 41038-41045. 10.1074/jbc.M205691200View ArticlePubMed
  48. Ernst J, Kellis M: Discovery and characterization of chromatin states for systematic annotation of the human genome. Nat Biotechnol 2010, 28: 817-825. 10.1038/nbt.1662PubMed CentralView ArticlePubMed
  49. Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M, Ku M, Durham T, Kellis M, Bernstein BE: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature 2011, 473: 43-49. 10.1038/nature09906PubMed CentralView ArticlePubMed
  50. Chen ZF, Paquette AJ, Anderson DJ: NRSF/REST is required in vivo for repression of multiple neuronal target genes during embryogenesis. Nat Genet 1998, 20: 136-142. 10.1038/2431View ArticlePubMed
  51. Westbrook TF, Martin ES, Schlabach MR, Leng Y, Liang AC, Feng B, Zhao JJ, Roberts TM, Mandel G, Hannon GJ, Depinho RA, Chin L, Elledge SJ: A genetic screen for candidate tumor suppressors identifies REST. Cell 2005, 121: 837-848. 10.1016/j.cell.2005.03.033View ArticlePubMed
  52. Gould SJ, Vrba E: Exaptation – a missing term in the science of form. Paleobiology 1982, 8: 4-15.
  53. Keith JH, Fraser TS, Fraser MJ Jr: Analysis of the piggyBac transposase reveals a functional nuclear targeting signal in the 94 c-terminal residues. BMC Mol Biol 2008, 9: 72. 10.1186/1471-2199-9-72PubMed CentralView ArticlePubMed
  54. Bobola MS, Silber JR, Ellenbogen RG, Geyer JR, Blank A, Goff RD: O6-methylguanine-DNA methyltransferase, O6-benzylguanine, and resistance to clinical alkylators in pediatric primary brain tumor cell lines. Clin Cancer Res 2005, 11: 2747-2755. 10.1158/1078-0432.CCR-04-2045View ArticlePubMed
  55. Fnu S, Williamson EA, De Haro LP, Brenneman M, Wray J, Shaheen M, Radhakrishnan K, Lee SH, Nickoloff JA, Hromas R: Methylation of histone H3 lysine 36 enhances DNA repair by nonhomologous end-joining. Proc Natl Acad Sci USA 2011, 108: 540-545. 10.1073/pnas.1013571108PubMed CentralView ArticlePubMed
  56. Beck BD, Lee SS, Williamson E, Hromas RA, Lee SH: Biochemical characterization of metnase’s endonuclease activity and its role in NHEJ repair. Biochemistry 2011, 50: 4360-4370. 10.1021/bi200333kPubMed CentralView ArticlePubMed
  57. Markljung E, Jiang L, Jaffe JD, Mikkelsen TS, Wallerman O, Larhammar M, Zhang X, Wang L, Saenz-Vash V, Gnirke A, Lindroth AM, Barrés R, Yan J, Strömberg S, De S, Pontén F, Lander ES, Carr SA, Zierath JR, Kullander K, Wadelius C, Lindblad-Toh K, Göran Andersson G, Hjälm G, Andersson L: ZBED6, a novel transcription factor derived from a domesticated DNA transposon regulates IGF2 expression and muscle growth. PLoS Biol 2009, 7: e1000256. 10.1371/journal.pbio.1000256PubMed CentralView ArticlePubMed
  58. Weiner AM, Gray LT: What role (if any) does the highly conserved CSB-PGBD3 fusion protein play in Cockayne syndrome? Mech Ageing Dev 2013, 134: 225-233. 10.1016/j.mad.2013.01.001PubMed CentralView ArticlePubMed
  59. Nieuwenhuys R: The brain of the lamprey in a comparative perspective. Ann N Y Acad Sci 1977, 299: 97-145. 10.1111/j.1749-6632.1977.tb41902.xView ArticlePubMed
  60. Holowachuk EW, Ruhoff MS: Efficient gene synthesis by Klenow assembly/extension-Pfu polymerase amplification (KAPPA) of overlapping oligonucleotides. PCR Methods Appl 1995, 4: 299-302. 10.1101/gr.4.5.299View ArticlePubMed
  61. Smith JJ, Antonacci F, Eichler EE, Amemiya CT: Programmed loss of millions of base pairs from a vertebrate genome. Proc Natl Acad Sci USA 2009, 106: 11212-11217. 10.1073/pnas.0902358106PubMed CentralView ArticlePubMed
  62. Smith JJ, Baker C, Eichler EE, Amemiya CT: Genetic consequences of programmed genome rearrangement. Curr Biol 2012, 22: 1524-1529. 10.1016/j.cub.2012.06.028PubMed CentralView ArticlePubMed
  63. Smith JJ, Kuraku S, Holt C, Sauka-Spengler T, Jiang N, Campbell MS, Yandell MD, Manousaki T, Meyer A, Bloom OE, Morgan JR, Buxbaum JD, Sachidanandam R, Sims C, Garruss AS, Cook M, Krumlauf R, Wiedemann LM, Sower SA, Decatur WA, Hall JA, Amemiya CT, Saha NR, Buckley KM, Rast JP, Das S, Hirano M, McCurley N, Guo P, Rohner N, et al.: Sequencing of the sea lamprey (Petromyzon marinus) genome provides insights into vertebrate evolution. Nat Genet 2013, 45: 415-421. 421e411-412 10.1038/ng.2568PubMed CentralView ArticlePubMed
  64. Amemiya CT, Alfoldi J, Lee AP, Fan S, Philippe H, Maccallum I, Braasch I, Manousaki T, Schneider I, Rohner N, Organ C, Chalopin D, Smith JJ, Robinson M, Dorrington RA, Gerdol M, Aken B, Biscotti MA, Barucca M, Baurain D, Berlin AM, Blatch GL, Buonocore F, Burmester T, Campbell MS, Canapa A, Cannon JP, Christoffels A, De Moro G, Edkins AL, et al.: The African coelacanth genome provides insights into tetrapod evolution. Nature 2013, 496: 311-316. 10.1038/nature12027PubMed CentralView ArticlePubMed
  65. Nikaido MNH, Nishihara H, Toyoda A, Suzuki Y, Kajitani R, Suzuki H, Okuno M, Aibara M, Ngatunga BP, Mzighani SI, Kalombo HWJ, Masengi KWA, Tudag J, Nogami S, Maedai R, Iwata M, Abe Y, Fujimura K, Okabe M, Amano T, Maeno A, Shiroishi T, Itoh T, Sugano S, Kohara Y, Fujiyama A, Okada N: Coelacanth genomes reveal signatures for evolutionary transition from water to land. Genome Res 2013. In press
  66. Yu JK, Wang MC, Shin IT, Kohara Y, Holland LZ, Satoh N, Satou Y: A cDNA resource for the cephalochordate amphioxus Branchiostoma floridae. Dev Genes Evol 2008, 218: 723-727. 10.1007/s00427-008-0228-xView ArticlePubMed
  67. Ebner B, Panopoulou G, Vinogradov SN, Kiger L, Marden MC, Burmester T, Hankeln T: The globin gene family of the cephalochordate amphioxus: implications for chordate globin evolution. BMC Evol Biol 2010, 10: 370. 10.1186/1471-2148-10-370PubMed CentralView ArticlePubMed
  68. Canestro C, Gonzalez-Duarte R, Albalat R: Minisatellite instability at the Adh locus reveals somatic polymorphism in amphioxus. Nucleic Acids Res 2002, 30: 2871-2876. 10.1093/nar/gkf386PubMed CentralView ArticlePubMed
  69. Yuasa HJ, Cox JA, Takagi T: Diversity of the troponin C genes during chordate evolution. J Biochem 1998, 123: 1180-1190. 10.1093/oxfordjournals.jbchem.a022059View ArticlePubMed


© Pavelitz et al.; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.