Concatemer-Associated Transgene Expression Patterns in Transgenic Marine Medaka Oryzias dancena Strains
Concatemer-Associated Transgene Expression Patterns in Transgenic Marine Medaka Oryzias dancena Strains
Fisheries and aquatic sciences. 2015. Mar, 18(1): 73-80
Copyright © 2015, The Korean Society of Fisheries and Aquatic Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial Licens ( permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : September 26, 2014
  • Accepted : January 26, 2015
  • Published : March 30, 2015
Export by style
Cited by
About the Authors
Young Sun, Cho
Institute of Marine Living Modified Organisms, Pukyong National University, Busan 608-737, Korea
Dong Soo, Kim
Institute of Marine Living Modified Organisms, Pukyong National University, Busan 608-737, Korea
Yoon Kwon, Nam
Center of Marine-Integrated Biomedical Technology (BK21 Plus Team), Pukyong National University, Busan 608-737, Korea

To examine the interrelationship between transgenic insertion patterns and transgene expression profiles in established transgenic fish lines, four stable transgenic marine medaka Oryzias dancena germlines harboring β-actin regulator-driven RFP reporter constructs were selected. The established transgenic strains were characterized with regard to their transgenic genotypes (insertion pattern, concatemer formation, and transgene copy number based on genomic Southern blot hybridization and qPCR assay) and expression characteristics at the mRNA (qRT-PCR), protein (western blot), and phenotypic (fluorescent appearance) levels. From comparative examinations, it was found that transgenic expression at both the transcription and translation levels could be significantly downregulated in transgenic strains, potentially through methylation-mediated transgene silencing that was particularly associated with the formation of a long tail -to- head tandem concatemer in the chromosomal integration site(s). When this occurred, an inverse relationship between the transgene copy number and fluorescence intensity was observed in the resultant transgenic fish. However, with the other transgenic genotype, transgenic individuals with an identical Southern blot hybridization pattern, containing a tandem concatemer(s), had very different expression levels (highly robust vs . low expression strengths), which was possibly related to the differential epigenetic modifications and/or degrees of methylation. The concatemer-dependent downregulation of transgene activity could be induced in transgenic fish, but the overall pattern was strain-specific. Our data suggest that neither a low (or single) transgene copy number nor tandem transgene concatemerization is indicative of strong or silenced transgene expression in transgenic fish carrying a ubiquitous transgene. Hence, a sufficient number of transgenic lineages, with different genotypes, should be considered to ensure the establishment of the best-performance transgenic line(s) for practical applications.
Experimental transgenesis in fish is useful in genomics research to utilize the recombinant animal model for acquiring the unique genotype and phenotype ( Hackett and Alvarez, 2000 ; Dunham and Winn, 2014 ). In addition, the transgenic manipulation of aquaculture fish species using gene(s) related to quantitative traits has received attention due to its potential as a powerful tool in the molecular genetic breeding of farmed fishes for productivity enhancement ( Nam et al., 2001 , 2007 ). The microinjection of DNA into fertilized eggs has long been utilized to produce transgenic fish strains, although novel methodologies to induce homologous and/or targeted transgenesis have been challenged recently in some fish models ( Rembold et al., 2006 ; Grabher and Wittbrodt, 2007 ; Casco-Robles et al., 2010 ). In microinjection-mediated gene delivery, it is widely accepted that the exogenously introduced DNA integrates into vertebrate genomes in an approximately random fashion. Multiple copies of the injected construct usually form a concatemer(s) in tandem arrays, in which the copy sequences are often incomplete at the ends, and DNA nibbling is sometimes observed between repeat units ( Bishop, 1996 ; Baup et al., 2010 ). The repeat number of transgene copies is reported to be unregulated and varies widely between different microinjection experiments within a given host species injected with identical constructs ( Würtele et al., 2003 ).
The context of transgene integration in the host genome provides the basis for the expression and functionality of transgenic individuals. Transgene expression has been found to be correlated with the integrant copy number in some cases, but studies in mammals have suggested it is related to heterochromatin-associated gene silencing ( Koetsier et al., 1996 ; Kong et al., 2009 ). Mammalian transgenics have shown that high-copy number transgenes arranged in a large tandem concatemer(s) are prone to structural modification, with a repressive influence on transgene expression, and the resulting transgenic individuals have often displayed inconsistent expression ( Dorer and Henikoff, 1997 ; Geurts et al., 2006 ). However, this inverse relationship between repetitive copy numbers and transgene expression has also been disputed ( Pena et al., 2004 ).
In contrast to the relatively rich information available in mammals, the interrelationship between transgene insertion patterns and expression profiles has been less well studied in established transgenic fish germlines, with only a few previous studies reported ( Wu et al., 2005 ; Uh et al., 2006 ). Unlike mammalian transgenics, high-copy number-related depressive regulation has not always been observed in transgenic fish strains. High-copy number transgenic fishes, carrying a long concatemer, have been reported to drive robust transgene expression, sometimes with markedly higher levels than those of single-copy or low-copy number transgenic fishes, particularly when measured absolutely ( Nam et al., 2001 ; Cho et al., 2013a ). For this reason, the copy number- and concatemer-dependent expression patterns that are common in mammalian transgenics might not apply to fish transgenics, and hence empirical data on the consequences of non-homologous, illegitimate recombination of transgenes are required in transgenic fish. The paucity of such information may be preventing the development of an effective method to establish stable and functional transgenic germlines.
Herein, we report the transgenic strain-specific patterns of either faithful or repressive expression, with the involvement of transgene concatemerization, as exemplified by marine medaka Oryzias dancena transgenics. Because the transgenically induced fluorescent phenotype could be easily profiled in a qualitative and quantitative fashion, we selected the best fit models from fluorescent transgenic marine medaka lines carrying a β-actin-driven red fluorescent protein (RFP) reporter construct, to scrutinize the potential concatemer-dependency of transgene expression patterns.
Materials and Methods
- Transgenic strains
Various transgenic marine medaka founders were generated by the microinjection of a fluorescent reporter construct consisting of the RFP gene (DsRed2; Clontech Laboratories Inc., Mountain View, CA, USA) under the control of endogenous β-actin regulator (podβ-actRFP; Cho et al., 2011 ). Microinjection was performed on one-cell embryos using the Cla I-linearized podβ-actRFP construct suspended at a concentration of 50 μg/ml. The functional ability of the β-actin regulator to drive ubiquitous expression of the red fluorescent reporter in this fish species has been described previously ( Cho et al., 2011 ). Details of the expression patterns of the endogenous β-actin gene, as well as the β-actin promoter-driven RFP expression in marine medaka, were provided in our previous work ( Cho et al., 2011 ). The β-actin-driven fluorescent signal in transgenic marine medaka usually becomes clearly visible by the embryonic stage (onset often occurs at late gastrulation), with the RFP phenotype lasting persistently to adulthood. Using various adult founder fishes, a series of preliminary evaluations were performed to select the appropriate transgenic line models for further characterization. Transgenic strains with extraordinarily high transgene copy numbers (>100 copies per cell) were not considered for further evaluation, because a comprehensive explanation of the transgene integration pattern at the Southern blot hybridization level would not have been feasible. Transgenic lines including those suspected to possess extrachromosomally persistent transgene copies were disregarded. Transgenic lines showing any physiological abnormalities, such as depressed fecundity or viability, were also excluded. Based on preliminary examinations, four transgenic genotype groups were selected for further characterization regarding the interrelationship between transgene integration pattern and transgene expression level.
- Characterization of the transgenic genotypes of selected strains
To examine the genomic integration and inheritance patterns of the transgene in the selected strains, a series of Southern blot hybridization analyses was performed. First, the absence of an extrachromosomally persistent transgene copy in the selected transgenic lines was confirmed by Southern blot hybridization using Bam HI, a non-cutter of the podβ-actRFP transgene. Second, a hybridization blot using Hin dIII (two-cutter of the transgene) was conducted to visualize transgene concatemerization and the potential junction site(s) between the transgene and host chromosome. Third, an internal transgene fragment of known size expected to be spliced from the podβ-actRFP irrespective of concatemerization was visualized by blot hybridization using Dra II (four-cutter). For all the Southern blot hybridization assays performed in this study, restriction enzyme-digested DNA was separated on 1% agarose gels, transferred to a positively charged nylon membrane (Roche Applied Science, Manheim, Germany), and hybridized using a 483-bp rfp fragment probe. Probe labeling with digoxygenin-11-UTP, hybridization, stringent washing, and signal detection were performed using a DIG DNA labeling and detection kit (Roche Applied Science) according to the manufacturer’s recommendations.
To assess the transgene copy numbers, qPCR was performed using 10 ng genomic DNA as a template. The thermal cycling reaction included 2 × iQ™ SYBR ® Green Supermix (Bio-Rad, Hercules, CA, USA) and a PCR primer pair (qDsRed-1F: CACAACACCGTGAAGCTGAA and qDsRed-1R: GAAGGACAGCTTCTTGTAG). The amplicon was a 141-bp internal rfp fragment ( Cho et al., 2013b ). The amplification signal was detected using an iCycler ® Real-Time Optical Module (Bio-Rad). Amplification was normalized to the endogenous β-actin gene actb (GenBank accession number HM347346 ) to ensure similar DNA input across samples ( Cho et al., 2013b ). Standard curves were prepared using 4-log dilutions (1 ‒ 1,000 copies) of positive podβ-actRFP plasmid, and a PCR efficiency greater than 90% was confirmed for each amplification reaction. Four F1 individuals from each transgenic group were examined using triplicate assays per DNA sample.
- Analysis of transgenic expression
Transgene expression in the selected lines was characterized at the mRNA (qRT-PCR assay of rfp transcripts), protein (Western blot of RFP protein) and fluorescent phenotype (external appearance and fluorescence microscopy) levels. Total RNA extracted from whole-body juveniles was reverse transcribed into cDNA using the Omniscript Reverse Transcription kit (Qiagen, Hilden, Germany), including oligo (dT) 20 primers and an O. dancena 18S rRNA reverse primer ( Cho et al., 2013a ). The prepared cDNA template was subjected to real-time thermal cycling to amplify either the 142-bp rfp fragment using a pair of oligonucleotides (qRT-RFP 1F: 5′-CACAACACCGTGAAGCTGAAG-3′ and qRT-RFP 1R: 5′-GGAAGGACAGCTTCTTGTAGT-3′) or the 188-bp 18S rRNA normalization control using qOD18S 1F (5′-AAGCTCGTAGTTGGATCTCG-3′) and qOD18S 1R (5′-CCTAGCTGCGGTATTCAGGC-3′) primers. The relative expression levels of rfp transcripts were based on normalization against 18S rRNA, as described previously ( Cho et al., 2011 ). The average PCR efficiency ranged from 92 to 99%. Differences between the transgenic genotypes were assessed using an ANOVA, followed by Duncan’s multiple range test at P = 0.05, based on triplicate examinations.
For Western blot analysis, the amount of soluble protein in the homogenate supernatant prepared with whole body fry was measured using a BCA protein assay kit (Pierce Biotechnology, Woburn, MA, USA). Resolved proteins from sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDSPAGE) were electrotransferred onto a nitrocellulose membrane (Amersham Pharmacia Biotech, Piscataway, NJ, USA) using a Trans-Blot ® SD semi-dry transfer cell (Bio-Rad) and probed with Living Colors ® DsRed monoclonal antibody (Clontech) diluted 1:2,000, and then with alkaline-phosphatase conjugated anti-mouse IgG secondary antibody (Santa Cruz Biotechnology Inc., Santa Cruz, CA, USA) diluted 1:4,000. Positive signals were visualized using a colorimetric method with nitroblue tetrazolium salt (Sigma-Aldrich, St Louis, MO, USA) and 5-bromo-4-chloro-3-idolyl phosphate (Sigma-Aldrich) and then analyzed using Quantity-One image analysis software (Bio-Rad).
At the fluorescent phenotypic level, the RFP signals in transgenic models were analyzed using an AZ100 fluorescence microscope equipped with NIS-Elements BR image analysis software (ver. 3.1; Nikon Corporation Instruments Company, Tokyo, Japan) as described previously ( Cho et al., 2011 ). The external appearance of RFP phenotypes in transgenic individuals were photographed using an EOD 5D Mark II digital camera (Canon, Tokyo, Japan) under either normal daylight or ultraviolet (UV) light (black light lamp; 352 nm; 15W; Sankyo Denki Co., Ltd., Tokyo, Japan). RFP signals in the internal organs (brain, eye lens, fin, gill, heart, intestine, kidney, liver, spleen, and gonads) were examined by fluorescence microscopy as described above.
Results and Discussion
Two germline positive founders (ACT201 and ACT203) passed on their fluorescent transgene to F1 progeny, with respective frequencies of 25 and 37% for RFP (+) F1 embryos, suggesting that the mosaic status in their germ cells was similar to those of most transgenic founders developed by microinjection ( Iyengar et al., 1996 ; Kim et al., 2004 ). However, interestingly, the expression pattern achieved in F1 individuals from each of the two founders was not uniform. Strikingly different RFP levels were apparent among F1 individuals from a given founder, in which some F1 offspring displayed very strong RFP signals whereas others exhibited only faint RFP expression by external appearance. From the ACT201 founder, the proportion of F1 progeny with strong RFP expression (labeled ACT201-S) of the total RFP (+) F1 fish was approximately 20%, while the remaining 80% RFP (+) F1 displayed only weak expression (labeled ACT201-W). Conversely, the ACT203 founder resulted in similar incidences (1:1 ratio) of strong (ACT203-S) and weak (ACT203-W) expression. The transgenic F1 fish selected from ACT201-S, ACT201-W, ACT203-S, or ACT203-W lines were able to transmit their phenotypic characteristics to F2 fish, exhibiting Mendelian frequency. Differences in the external fluorescent appearance between strong and weak RFP types within a given lineage were clearly apparent ( Fig. 1 ). The red fluorescent signals visualized over the whole body in transgenic individuals belonging to either the ACT201-S or ACT203-S lines were stunningly bright, and red fluorescence could be observed easily with the unaided eye under normal daylight conditions. Conversely, the ACT201-W and ACT203-W lines had only moderate or weak red fluorescence in their external appearance, although they could easily be distinguished from non-transgenic, wild-type individuals. Notably, several individuals belonging to the ACT203-W line lost the clear RFP signal in their eyes (data not shown). From microscopic examinations of the internal organs, transgenic individuals usually retained the ubiquitous distribution pattern over a wide array of tissues examined ( Fig. 1 ), indicating that the transgenically introduced β-actin promoter could induce expression resembling the pattern of the endogenous β-actin gene ( Cho et al., 2011 ; Lee et al., 2012 ). However, in several transgenic individuals belonging to the ACT203-W line, the presence of RFP (-) organs was determined under fluorescence microscopy conditions, although the organs lacking RFP signals (eye, heart, spleen, and/or testis depending on individuals; photograph not shown) were not uniform among these transgenic individuals.
PPT Slide
Lager Image
Representative image to show the phenotypic fluorescence achieved in each transgenic genotype group (ACT201-S, ACT201-W, ACT203-S and ACT203-W). Mosaic RFP expression patterns of the two founder fishes (ACT201-F0 and ACT203-F0) are also shown on the top along with their F2 progeny. Photograph for external appearance of each fish is taken under normal daylight conditions. Fluorescence microscopic images for selected organs are also provided at the bottom, as exemplified by the fish belonging to the line ACT201-S.
Transgenic fish belonging to one of the four phenotypic groups ( i.e ., ACT201-S, ACT201-W, ACT203-S, and ACT203-W) were subjected to genomic Southern blot hybridization assays using different restriction strategies ( Fig. 2 ). Southern blot hybridization using a non-cutter enzyme Bam HI confirmed that none of the transgenic individuals (F1 fish) showed a notable sign of the presence of extrachromosomally persistent copies, as evidenced by the presence of a single hybridization signal in the high molecular weight range. Southern blot hybridization using Hin dIII suggested that the ACT201 founder experienced two different, independent genomic integration events during its embryonic development, consequently leading to the chimeric status in germ cell populations at adulthood. In the ACT201 lineage, the ACT201-S line showed two hybridization-positive bands at positions 3.0 and 4.7 kb (total copy number of the transgene = 2.41 based on a qPCR assay), while the ACT201-W line showed six bands in the range of 2.3 to 9.0 kb (total copy number of the transgene = 26.25 based on a qPCR assay). Considering the restriction sites recognized by Hin dIII, the relatively strong bands at 4.5 kb seen in ACT201-W F1 individuals could be a sign of tail -to- head concatemerization ( i.e ., joining of the 2.31-kb tail to the 2.18-kb head ), suggesting that the transgenic insertion site contained a large number of tandem copies ( Fig. 2 ). Because the F1 fishes transmitted the fluorescent transgene to the next generation following the Mendelian ratio of single gene inheritance, these multiples in the concatemer(s) might insert into a single chromosomal site, or into very close neighboring sites within a chromosome ( Nam et al., 1999 ; Cho et al., 2013a ). In contrast, the ACT201-S F1 individuals displayed two hybridization signals. One was visualized slightly above the 4.5-kb tail -to- head concatemer-associated band seen in the ACT201-W group. The other hybridization signal was visualized at 3.0 kb in the ACT201-S group. Taking into account the similar hybridization intensities between the two bands, the approximate copy number of integrants was close to 2, which was in agreement with the copy number estimated by qPCR (average of 2.41 copies per cell). These bands might reflect the junction regions between transgene and host chromosomes, although direct cloning is required to verify the detailed structure of the junction sites. Each of the ACT201 F1 types ( i.e ., -S and -W genotypes) was successfully transmitted from the F1 to F2 generations, without structural modifications seen at the Southern blot levels, suggesting that the potential hemizygous genotype for the transgene could be established in each phenotypic lineage.
PPT Slide
Lager Image
Genomic Southern blot hybridization analyses of marine medaka Oryzias dancena transgenic lines ACT201-S, ACT201-W, ACT203-S and ACT203-W harboring the podβ-actRFP transgene construct. Molecular size markers (in kb) are based on lambda-HindIII digests and/or l kb ladder plus. BamHI is a non-cutter of the podβ-actRFP transgene. Relative positions of HindIII and DraII sites on the podβ-actRFP transgene (ClaI-linearized form) are indicated in a schematic drawing on the top to show an example of tail-to-head concatemerization. In that drawing, the fragments targeted by the probe (labeled P) either in DraII (1.24-kb) or HindIII (4.49-kb and > 2.31-kb) digestion are also indicated. In the HindIII-blot, the 4.5-kb hybridization bands resulting from the tail-to-head concatemerization appeared in ACT201-W and ACT203-S/W lines are indicated by arrows. The genomic DNA run in the first lane is the same DNA from F1 transgenic fish from ACT203-S (marked with open circles) to confirm that the same size of the concatemer-band from ACT201-W and ACT203-S/W lines.
Unlike the ACT201 lineages, only a single hybridization pattern was found in the Hin dIII blot for F1 transgenic individuals from the ACT203 founder, irrespective of the RFP phenotype. The Hin dIII digestion of the ACT203 F1 transgenic line produced at least eight recognized bands ranging from 2.2 to 8.0 kb without notable variability among F1 individuals (average copy number = 42.10 based on qPCR). Similarly to the pattern from the ACT201-W group, all the F1 transgenic individuals belonging to either the ACT203-S or ACT203-W groups showed a common, hybridization signal at 4.5 kb, which corresponded to the band resulting from tandem tail -to- head concatemerization of multiple copies. The intensity of the concatemer-originating band in the ACT203 lineage was stronger than that in the ACT201-W line, suggesting that the ACT203 transgenic individuals might possess a larger concatemer(s) compared with the ACT201-W group. The F1 fish belonging to the ACT203-S group were successful in transmitting their hybridization pattern to F2 fish without any changes, while some F1 individuals from the ACT203-W group showed a small modification of the hybridization pattern at 5.0 kb during transmission to the F2 generation
Another Southern blot assay using the Dra II enzyme showed that the 1.24-kb internal fragments were spliced out successfully from most transgene copies for each transgenic strain. However, unlike ACT201-S, in the other three lines (ACT201-W, ACT203-S, and ACT203-W containing a long transgene concatemer), some transgene copies were resistant to digestion by Dra II. In addition to the expected 1.24-kb band, the ACT201-W group showed a weak hybridization signal at 6.6 kb and two bands at 2.3 kb. Four unexpected hybridization bands at approximately 9.0, 6.0, 5.5, and 2.5 kb were visualized in both the ACT203-S and ACT203-W transgenic lines. Dra II (but not Hin dIII) is a methylation-sensitive enzyme, inhibited at purine-G/GNC * C-pyrimidine by 5´-methylcytosine as indicated (*) ( Schnetz and Rak, 1988 ). Hence, the resistance of such copies to Dra II digestion, as detected in the concatemer-containing transgenic lines, might reflect methylation of the transgene, although the possibility of other undefined modification(s), including the destruction of Dra II sites in certain transgene copies, cannot be discounted. Concatemer-mediated methylation of other transgenes has also been reported, and methylation-associated epigenetic modification is considered to be an important factor affecting the downregulation and silencing of transgene activities ( Matzke et al., 2000 ). After examination up to the F2 generation, the ACT201-S line was selected and further propagated to verify consistent transmission of its transgenic genotype through numerous subsequent generations to confirm the stable establishment of the functional transgenic insertion site.
The transgenic lines ACT201-S/W and ATC203-S/W were subjected to expression analyses at the mRNA (qRT-PCR), protein (Western blot), and phenotype (fluorescent appearance) levels. qRT-PCR of rfp transcripts indicated the highest mRNA expression in the ACT201-S line, followed by ACT203-S ( P <0.05). The ACT201-W and ACT203-W lines showed similarly weak expression of rfp transcripts ( P >0.05) ( Fig. 3A ). The mRNA expression pattern was also well correlated with the protein expression pattern observed by Western blot; in terms of the significantly abundant expression of RFP in ACT201-S and ACT203-S compared with the low expression in ACT201-W and ACT203-W ( Fig. 3B ). Collectively, the expression data assessed at the mRNA and protein levels were congruent with the phenotypic fluorescent characteristics acquired in the transgenic individuals.
PPT Slide
Lager Image
Quantitative analyses of rfp transcripts using qRT-PCR (A) and RFP protein using western blot (B). In (A), the rfp mRNA levels in F2 generations (n = 6 fish per genotype group) were normalized against 18S rRNA control. Means ± SDs with different letters are significantly different based on ANOVA followed by Duncan’s multiple range tests at P = 0.05. In (B), the representative blot image to show the differential amounts of RFP (25.7 kDa) across the transgene genotype groups (ACT201-S, ACT201-W, ACT203-S and ACT203-W). To validate the similar amount of protein loading, the actin protein was visualized at 48 kDa as an internal control by using universal anti-ACTIN antibody (Santa Cruz Biotechnology Inc.).
Taken together, the data presented here suggest strain-dependent patterns for the interrelationship between insertion structure and expression characteristics of the transgene expressed in marine medaka. As evidenced in the ACT201-S/W families, transgenic expression could be downregulated in transgenic fish strains by the formation of a long concatemer. Low-copy number transgenic fish without a long concatemer (ACT201-S) had a significantly higher level of transgene expression than that of their siblings with high transgene copy numbers in the tail -to- head concatemer (ACT201-W). Because these two unique transgenic genotypes originated from a single founder, the genetic background of the hosts does not need to be seriously contemplated as a potential factor affecting transgene expression. The main reason for the inverse relationship between transgene copy number and expression strength in the ACT201 lineages might be the long concatemer-associated epigenetic methylation at the transgene insertion sites (as supported by Dra II digestion data), consequently leading to significantly depressed transgene activity in the transgenic fish carrying a long tandem array of transgene copies. However, this concatemer-mediated regulation pattern is not always displayed in all transgenic lines ( Cho et al., 2011 , 2013b ). As evidenced by the ACT203 families in this study, the identical transgene genotype seen by Southern blot might not be indicative of similar expression characteristics among the transgenic siblings. Unfortunately, the specific mechanism behind the different expression patterns between the ACT203-S and ACT203-W lineages is not yet understood. It is possible that the differential degree of methylation, especially in the promoter region, between the two lineages, could not be visualized simply by the Dra II-blot. The exclusive occurrence of RFP-negative organs in certain individuals of the ACT203-W line might be, at least in part, related to differences in methylation patterns among individuals of this line. A more quantitative approach to measure the degree of methylation, including bisulfite sequencing potentially, is needed to obtain a deeper insight into the transgene methylation status ( Feng et al., 2011 ). In vivo treatment of methyltransferase inhibitor(s) could also be a rational approach to determine whether the expression profile of the transgene can be modified in certain transgenic lines by methylation-associated silencing of the transgene ( Gibbs et al., 1994 ).
This study was supported by Korea Institute of Marine Science & Technology Promotion (Project #20130369).
Baup D , Fraga L , Pernot E , Acker AV , Vanherck AS , Breckpot K , Thielemans K , Schurmans S , Moser M , Leo O 2010 Variegation and silencing in a lentiviral-based murine transgenic model Transgenic Res 19 399 - 414    DOI : 10.1007/s11248-009-9318-4
Bishop JO 1996 Chromosomal insertion of foreign DNA Repd Nutr Dev 36 607 - 618
Casco-Robles MM , Yamada S , Miura T , Chiba C 2010 Simple and efficient transgenesis with I-SceI meganuclease in the newt, Cynops pyrrhogaster Dev Dyn 239 3275 - 3284    DOI : 10.1002/dvdy.22463
Cho YS , Kim DS , Nam YK 2013 Characterization of estrogen-responsive transgenic marine medaka Oryzias dancena germlines harboring red fluorescent protein gene under the control by endogenous choriogenin H promoter Transgenic Res 22 501 - 517    DOI : 10.1007/s11248-012-9650-y
Cho YS , Lee SY , Kim DS , Nam YK 2013 Characterization of stable fluorescent transgenic marine medaka (Oryzias dancena) lines carrying red fluorescent protein gene driven by myosin light chain 2 promoter Transgenic Res 22 849 - 859    DOI : 10.1007/s11248-012-9675-2
Cho YS , Lee SY , Kim YK , Kim DS , Nam YK 2011 Functional ability of cytoskeletal β-actin regulator to drive constitutive and ubiquitous expression of a fluorescent reporter throughout the life cycle of transgenic marine medaka Oryzias dancena Transgenic Res 20 1333 - 1335    DOI : 10.1007/s11248-011-9501-2
Dorer DR , Henikoff S 1997 Transgene repeat arrays interact with distant heterochromatin and cause silencing in cis and trans Genetics 147 1181 - 1190
Dunham RA , Winn RN , Pinkert CA 2014 Transgenic animal technology - A laboratory handbook 3rd ed Elsevier Inc. Production of transgenic fish 305 - 334
Feng S , Rubbi L , Sacobsen SE , Pellegrini M 2011 Determining DNA methylation profiles using sequencing Methods Mol Biol 733 223 - 238
Geurts AM , Collier LS , Geurts JL , Oseth LL , Bell ML , Mu D , Lucito R , Godbout SA , Green LE , Lowe SW , Hirsch BA , Leinwand LA , Largaespada DA 2006 Gene mutations and genomic rearrangements in the mouse as a result of transposon mobilization from chromosomal concatemers PLoS Genet 2 e156 -    DOI : 10.1371/journal.pgen.0020156
Gibbs PD , Peek A , Thorgaard G 1994 An in vivo screen for the luciferase transgene in zebrafish Mol Mar Biol Biotechnol 3 307 - 316
Grabher C , Wittbrodt J 2007 Meganuclease and transposon mediated transgenesis in medaka Genome Biol 8 S10 -    DOI : 10.1186/gb-2007-8-s1-s10
Hackett PB , Alvarez MC , Fingerman M , Nagabhushanam R 2000 Recent advances in marine biotechnology Science Publishers Enfield The molecular genetics of transgenic fish 77 - 145
Iyengar A , Muller F , Maclean N 1996 Regulation and expression of transgenes in fish - a review Transgenic Res 5 147 - 166    DOI : 10.1007/BF01969704
Kim DS , Kim BS , Lee SJ , Park IS , Nam YK 2004 Comparative analysis of inherited patterns of the transgene in transgenic mud loach Misgurnus mizolepis lines carrying the CAT reporter gene Fish Sci 70 201 - 210    DOI : 10.1111/j.1444-2906.2003.00792.x
Koetsier PA , Mangel L , Schmitz B , Doerfler W 1996 Stability of transgene methylation patterns in mice: position effects, strain specificity and cellular mosaicism Transgenic Res 54 235 - 244
Kong Q , Wu M , Huan Y , Zhang L , Liu H , Bou G , Luo Y , Mu Y , Liu Z 2009 Transgene expression is associated with copy number and cytomegalovirus promoter methylation in transgenic pigs PLoS One 4 e6679 -    DOI : 10.1371/journal.pone.0006679
Lee SY , Kim DS , Nam YK 2012 Molecular characterization of cytoskeletal beta-actin and its promoter in the Javanese ricefish Oryzias javanicus Fish Aquat Sci 15 317 - 324
Matzke MA , Mette MF , Matzke AJM 2000 Transgene silencing by the host genome defense: implications for the evolution of epigenetic control mechanisms in plants and vertebrates Plant Mol Biol 43 401 - 415    DOI : 10.1023/A:1006484806925
Nam YK , Maclean N , Fu C , Pandian TJ , Eguia MRR , Kapuscinski AR , Hayes KR , Li S , Dana G 2007 Environmental risk assessment for genetically modified organisms CABI Press Cambridge Development of transgenic fish: scientific background 61 - 94
Nam YK , Noh JK , Cho YS , Cho HJ , Cho KN , Kim CG , Kim DS 2001 Dramatically accelerated growth and extraordinary gigantism of transgenic mud loach Misgurnus mizolepis Transgenic Res 10 353 - 362    DOI : 10.1023/A:1016696104185
Nam YK , Noh CH , Kim DS 1999 Transmission and expression of an integrated reporter construct in three generations of transgenic mud loach (Misgurnus mizolepis) Aquaculture 172 229 - 245    DOI : 10.1016/S0044-8486(98)00433-5
Pena RN , Webster J , Kwan S , Korbel J , Whitelaw BA 2004 Transgene methylation in mice reflects copy number but not expression level Mol Biotechnol 26 215 - 219    DOI : 10.1385/MB:26:3:215
Rembold M , Lahiri K , Foulkes NS , Wittbrodt J 2006 Transgenesis in fish: efficient selection of transgenic fish by co-injection with a fluorescent reporter construct Nat Protoc 1 1133 - 1139    DOI : 10.1038/nprot.2006.165
Schnetz K , Rak B 1988 Cleavage by EcoO109 and DraII is inhibited by overlapping dcm methylation Nucleic Acids Res 16 1623 -    DOI : 10.1093/nar/16.4.1623
Uh M , Khattra J , Devlin RH 2006 Transgene constructs in coho salmon (Oncorhynchus kisutch) are repeated in a head-to-tail fashion and can be integrated adjacent to horizontally-transmitted parasite DNA Transgenic Res 15 711 - 727    DOI : 10.1007/s11248-006-9016-4
Wu B , Sun YH , Wang YW , Wang YP , Zhu ZY 2005 Characterization of transgene integration pattern in F4 hGH-transgenic common carp (Cyprinus carpio L.) Cell Res 15 447 - 454    DOI : 10.1038/
Würtele H , Little KCE , Chartrand P 2003 Illegitimate DNA integration in mammalian cells Gene Therapy 10 1791 - 1799    DOI : 10.1038/