Corky, a gypsy-like retrotransposon is differentially transcribed in Quercus suber tissues
© Rocheta et al.; licensee BioMed Central Ltd. 2012
Received: 18 March 2012
Accepted: 2 August 2012
Published: 13 August 2012
Transposable elements (TEs) make up a large part of eukaryotic genomes. Due to their repetitive nature and to the fact that they harbour regulatory signals, TEs can be responsible for chromosomal rearrangements, movement of gene sequences and evolution of gene regulation and function. Retrotransposon ubiquity raises the question about their function in genomes and most are transcriptionally inactive due to rearrangements that compromise their activity. However, the activity of TEs is currently considered to have been one of the major processes in genome evolution.
We report on the characterization of a transcriptionally active gypsy-like retrotransposon (named Corky) from Quercus suber, in a comparative and quantitative study of expression levels in different tissues and distinct developmental stages through RT-qPCR. We observed Corky’s differential transcription levels in all the tissues analysed.
These results document that Corky’s transcription levels are not constant. Nevertheless, they depend upon the developmental stage, the tissue analysed and the potential occurring events during an individuals’ life span. This modulation brought upon by different developmental and environmental influences suggests an involvement of Corky in stress response and during development.
KeywordsQuercus suber Corky LTR retroelement RT-qPCR Expression
Retrotransposons are generally the most abundant class of Transposable Elements (TEs), concerning their proportion in the genomes and, are widely distributed among eukaryotic genomes, especially in plants . Due to their wide distribution and the diverse types of induced mutations, TEs are thought to have contributed significantly to eukaryotic genes and genomes evolution . The increasing number of data obtained from genome-wide sequencing projects indicate that TEs take part in major events and are a potential pool of promoter regions for host regulatory sequences . TE regulatory regions are known to be sequences of extremely rapid evolution, a characteristic of eukaryotic regulatory regions attributed to having to cope with changing genomic environments . LTR-retrotransposons are 'copy-and-paste' (class I) TEs that replicate via an RNA intermediate. Like animal retroviruses, these retrotransposons have two LTRs, with signals for transcription initiation and termination, flanking an internal region (gag-pol) that typically contains genes and other features necessary for autonomous retrotransposition. Retrotransposon ubiquity raises the question about their function in genomes. Retrotransposon insertions in, or next to coding regions, generate mutations that can lead to changes in gene expression. For instance, Tnt 1A transposition preferentially targets genic regions, suggesting that the activity of transposable elements can modulate genic functions and represent a natural source of phenotypic diversity . Furthermore, run-off transcription from retrotransposons can lead to overexpression or suppression of nearby genes . Transcription activity detected in several retrotransposons during certain stages of development seems to point to a potential role of these elements during plant growth [7, 8]. Additionally, some biotic and abiotic stresses can increase transcript levels of retroelements, such as tobacco Tnt1, Tto1, Tto2, rice´s Tos17 and Rtsp-1 from sweet potato . An overall picture of retrotransposon expression is however difficult to establish due to the absence of exhaustive comparative studies in different tissues. Several Gypsy and Copia- like retroelements are known to be well represented in the Mediterranean Quercus suber[13, 14].
IFG 7  is one of the most representative Gypsy-like elements in coniferous genomes such as in several Pines [15–18]and Taxodium distichum, and sometimes is considered as a conifer-specific LTR retroelement . However, elements like IFG 7 were not yet identified in Angiosperms. In order to study the possible occurrence of a conifer derived LTR retroelement in a distant related Angiosperm tree species, as well as its potential active transcriptional activity in this species, we used IFG 7 as a Gypsy representative element.
The key aims of this work were the molecular characterization of a new retrotransposon in the Quercus suber genome which is homologue to the previously identified IFG 7 from Pinus radiata and PpRT 1 from Pinus pinaster and, the quantification of its transcriptional activity in different tissues and distinct developmental stages and conditions. Together, the data presented here clearly show that this retrotransposon, named Corky, makes up a dynamic component of the cork oak genome.
Organization and structure of Corky
Corky is a gypsy retroelement that was isolated throughout genome walking in Q. suber genome. All generated DNA fragments were sequenced and further analysed.
Each LTR is 333 bp long and is flanked by a short 7 bp direct repeat 5’- CTCGATG-3’ (Figure 1b), probably representing a duplication of the genomic target site produced by the insertion of a Corky copy, such as it has been reported for other retroelements . Both LTRs begin and end with a 5 bp inverted repeat 5’TGTTA…TAACA-3’ including the retroviral consensus 5’-TG…CA-3’. LTRs inverted repeats are present in all retroviruses and are thought to be important for their integration  (Figure 1b). 5’LTR´s Corky sequence analysis (Additional file 2) revealed two characteristic patterns of repeating motifs: one is a simple pattern of short tandem sequence motifs (a..a..) TA(G)TGATTACCCC(A)T(T)(A) and TA(T)TG(T)ATTA(TA)CCCC(T)T(A)(T), while the other one, more complex, has two adjacent heterologous motifs (TATTGTTA, TTATATT), repeated twice as a group (ab..ab), as present within the HIV-1 and gypsy enhancers . Both patterns are dispersed between the two TATA sequences (TATATATA) (Additional file 2). Enhancers typically consist of a series short repeated sequence motifs that are often associated with regulatory protein binding domains .
Quantification of Corky expression
Our data show good evidence that a retrotransposon (Corky) has escaped from host silencing mechanisms. The differential expression in several plant tissues in different developmental stages suggests, at least, an involvement of this retrotransposon in stress response and in developmental processes. It is likely that retroelements do not increase plasticity in an evolutionarily active way but they might play a crucial role in response to developmental/environmental challenges. Together, these results set the need to further investigate both regulation and control mechanisms that implicate retrotransposons and development.
Materials and methods
Acorns of Quercus suber L. produced by open pollination and pollen used in this study were collected in a natural population at Alcácer do Sal (Portugal). The plants used in this study were obtained from those acorns and grown in the greenhouse until they were used (at 2.5 years old). Plant tissues were frozen in liquid nitrogen. Genomic DNA was extracted from samples using DNeasy® Plant Mini Kit (Qiagen®), according to the manufacturer’s instructions.
Initial DNA amplification strategy
The first set of primers [Forward- 5’ttcaactgagtcaaatttc3’ and Reverse- 5’ctgtcaacccaagaaatcctcgcag 3’] (Additional file 3) used, were constructed by the assumption that the RVT sequence in Q. suber has sufficient similarities with the previous retrotransposon amplified in P. pinaster (named PpRT 1) . For this part of the work only DNA from young leaves was handled. A set of primers was designed to guarantee that we are in the presence of the same copy of Corky (Additional file 3). The PCR protocol consisted of the subsequent steps: an initial denaturation period at 94°C for 4 min., 30 cycles of amplification, each of which consisted of 45 s of denaturation at 94°C, 45 s of annealing at 57°C, and 90 s of elongation at 72°C with a final elongation step of 4 min at 72°C. After purification with the QIAquick® PCR purification kit, the amplified fragment was cloned using pCR 2.1-TOPO vector (Invitrogen®) and sequenced.
Genome walking was performed using the Genome Walker® kit (Clontech®) components according to the manufacturer’s instructions. The amplification of upstream and downstream regions of RVT sequences from the libraries was performed also according to the Genome Walker® Kit protocol and the primers melting temperature (Additional file 3). All the PCR amplifications were performed with the proofreading enzyme Phusion (New England Biolabs®). The major PCR products obtained were gel extracted by the Gel Extraction® Kit (Qiagen®) additionally inserted in pCR 2.1-TOPO vector® (Invitrogen®), sequenced and aligned using the online service of National Center for Biotechnology Information (NCBI) . To guarantee that all sequences belong to the same retroelement we performed numerous amplifications for the same region with different sets of primers. Additionally, primers were designed assuring that all fragments amplified overlap. Thus, all the fragments obtained were used to assemble the entire retroelement. Conversely, without other resources such as Bacterial Artificial Chromosomes (BACs), we cannot say that we have isolated the same genomic element. Although, the high overlap of the individual sequences ensures that we have got the same element, we cannot discard the hypothesis that we have reconstructed a chimeric sequence. The assembled sequence was used to search all the retrotransposon regions between both LTRs, according to the conserved motives.
RNA isolation and cDNA preparation for RT-qPCR
Total RNA was extracted from secondary roots, old leaves (one year old) and young leaves (from the year) from ten 2,5 year old plants, from ten dormant embryos, from the primordia of leaves and roots of ten germinated embryos and from ten different pollen samples, each replicate corresponding to tissue originating from one single plant and also from ten wounded leaves (leaves were pierced with a needle 240 min prior to freezing), using the RNAqueous® kit (Ambion®), according to the manufacturer’s instructions. Nucleic acid concentration of each sample was quantified by spectrophotometry using the software Gen5 1.09 (Synergy HT, Bio-Tek Instruments, Winooski, USA). Total RNA quality was assessed by the A260/A280 and A260/A230. Only RNA samples with A260/A280 between 1.8 and 2.1 and A260/A230 between 2.0 and 2.2 were accepted for the experience. Total RNA integrity was tested through 1% agarose gel electrophoresis under denaturing conditions.
RNA samples were treated with RQ1 RNase-Free DNase (Promega, Madison, WI). cDNA was synthesized from 2 μg of total RNA using random hexamers and Superscript II RNase H- reverse transcriptase (Invitrogen®, Carlsbad, CA), according to the manufacturer’s recommendations followed by PCR amplification using specific primers for the RVT and a region between Integrase and the chromodomain of Corky (Figure 2). As expected, amplification products were not obtained in RNA samples not yielded to reverse transcription prior to PCR. cDNA was stored at −20°C.
Transcriptional activity of Corky
RT-qPCR was performed in a 96 well white reaction plates (Bio-Rad®, Hercules, CA), using an IQ5 Real Time PCR (Bio-Rad®, Hercules, CA) with ten biological replicates and two technical replicates. For amplification specific primers corresponding to the RVT domain of Corky and a region between the Integrase and the chromodomain were used (Figure 2). Each 20 μL reaction mixture well contained 10.0 μL of 2x master mix iQ SYBR Green Supermix®, 2.0 μL of HPLC-purified primers (10 μM), 7.0 μL of PCR-grade H2O and 1.0 μL target DNA solution. PCR amplification products were monitored via intercalation of SYBR-Green (included in the master mix). The PCR protocol consisted of an initial denaturation step at 95°C for 3 min, 40 cycles of amplification, each of which consisted of 15 s of denaturation at 95°C, 20 s of annealing at 57°C and 50 s of elongation at 72°C. As expected, amplification products were not obtained in RNA samples not subjected to the reverse transcription step prior to PCR.
To assess the primers amplification efficiency, identical volumes of cDNA samples were diluted and used to generate five-point standard curves based on a five-fold dilution series (1; 1:5; 1:25; 1:125; 1:625), in triplicate. Amplification efficiency (E) is calculated as E = 10(−1/a)-1, “a” being the slope of the linear regression curve (y = a log (x) + b) fitted over the log-transformed data of the input cDNA dilution (y) plotted against the respective quantification cycle (Cq) values (x). E-values of the target genes were considered comparable when they did not exceed 100 ± 10%, corresponding to a standard curve slope of 3.3 ± 0.33. All cDNA samples were diluted 50 fold and were amplified in duplicate in two independent PCR runs.
To generate a baseline-subtracted plot of the logarithmic increase in fluorescence signal (ΔRn) versus cycle number, baseline data were collected between the cycles 5 and 17. All amplification plots were analysed with an R n threshold of 0.2, at the beginning of the region of exponential amplification, to obtain Cq and the data obtained were exported into a MS Excel workbook (Microsoft® Inc.) for further analysis. In order to compare data from different PCR runs or cDNA samples, Cq values were normalized to the Cq value of actin, a housekeeping gene expressed at a relatively high and constant level . Gene expression was calculated using the ΔΔCq method . Results are expressed as fold variation of each tissue relative to each of the other.
reverse transcription real time PCR
We thank “Fundação para a Ciência e Tecnologia” for the post-doc grant SFRH/BPD/5707/2001 to L.C., for post-doc grant SFRH/BPD/64905/2009 to M. R. and for the “Plurianual” funding to CBAA and PTDC/AGR-GFL/104197/2008. We are extremely grateful to Sara Amâncio for the cheer encouragement and helpful suggestions. We also give our thanks’ Quirina Santos-Costa for the English revision.
- Hua-Van A, Rouzic AL, Maisonhaute C, Capy C: Abundance, distribution and dynamics of retrotransposable elements and transposons: similarities and differences. Cytogenet Genome Res. 2005, 110: 426-440. 10.1159/000084975.PubMedView ArticleGoogle Scholar
- Kazazian HH: Mobile elements: drivers of genome evolution. Science. 2004, 303: 1626-1632. 10.1126/science.1089670.PubMedView ArticleGoogle Scholar
- Fablet M, Souames S, Biémont C, Vieira C: Evolutionary pathways of the tirant LTR retrotransposon in the Drosophila melanogaster subgroup of species. J Mol Evol. 2007, 64: 438-447. 10.1007/s00239-006-0108-9.PubMedView ArticleGoogle Scholar
- Ludwig MZ, Bergman C, Patel NH, Kreitman M: Evidence for stabilizing selection in a eukaryotic enhancer element. Nature. 2000, 403: 564-567. 10.1038/35000615.PubMedView ArticleGoogle Scholar
- Grandbastien MA, Audeon CE, Bonnivard JM, Casacuberta B, Chalhoub APP, Costa QH, Lea D, Melayah M, Petit C, Poncet SM: Stress activation and genomic impact of Tnt 1 retrotransposons in Solanaceae. Cytogenet Genome Res. 2005, 110: 229-241. 10.1159/000084957.PubMedView ArticleGoogle Scholar
- Kashkush K, Feldman M, Levy AA: Transcriptional activation of retrotransposons alters the expression of adjacent genes in wheat. Nat Genet. 2003, 33: 102-106. 10.1038/ng1063.PubMedView ArticleGoogle Scholar
- Tahara M, Aoki T, Suzuka S, Yamashita H, Tanaka M, Matsunaga S, Kokumai S: Isolation of an active element from a high-copy-number family of retrotransposons in the sweetpotato genome Mol Genet Genomics. 2004, 272: 116-127.Google Scholar
- Rico-Cabanas L, Martínez-Izquierdo JA: Cire 1, a novel transcriptionally active Ty 1-copia retrotransposon from Citrus sinensis. Mol Genet Genomics. 2007, 277: 365-377. 10.1007/s00438-006-0200-2.PubMedView ArticleGoogle Scholar
- Beguiristain T: Grandbastien MA, Puigdomenech P, Casacuberta JM: Three Tnt 1 subfamilies show different stress-associated patterns of expression in tobacco. Consequences for retrotransposon control and evolution in plants. Plant Physiol. 2001, 127: 212-221.Google Scholar
- Takeda S, Sugimoto K, Otsuki H, Hirochika H: A 13-bp cis-regulatory element in the LTR promotor of the tobacco retrotransposob Tto 1 is involved in responsiveness to tissue culture, wounding, methyl jasmonate and fungal elicitors. Plant J. 1999, 18: 383-393. 10.1046/j.1365-313X.1999.00460.x.PubMedView ArticleGoogle Scholar
- Hirochika H: Activation of tobacco retrotransposons during tissue culture. Embo J. 1993, 12: 2521-2528.PubMedPubMed CentralGoogle Scholar
- Hirochika H, Sugimoto K, Otsuki Y, Tsugawa H, Kanda M: Retrotransposons of rice involved in mutations induced by tissue culture. Proc Natl Acad Sci USA. 1996, 93: 7783-7788. 10.1073/pnas.93.15.7783.PubMedPubMed CentralView ArticleGoogle Scholar
- Alves S, Ribeiro T, Inácio V, Rocheta M, Morais-Cecílio L: Genomic organization and dynamics of repetitive DNA sequences in representatives of three Fagaceae genera. Genome. 2012, 55: 348-359. 10.1139/g2012-020.PubMedView ArticleGoogle Scholar
- Carvalho M, Ribeiro T, Viegas W, Morais-Cecílio L, Rocheta M: Presence of env-like sequences in Quercus suber retrotransposons. J Appl Genet. 2010, 51: 461-467. 10.1007/BF03208875.PubMedView ArticleGoogle Scholar
- Kossack DS, Kinlaw CS: IFG, a gypsy-like retrotransposon in Pinus (Pinaceae), has an extensive history in pines. Plant Mol Biol. 1999, 39: 417-426. 10.1023/A:1006115732620.PubMedView ArticleGoogle Scholar
- Kovach A, Wegrzyn JL, Parra G, Holt C, Bruening GE, Loopstra CA, Hartigan J, Yandell M, Langley CH, Korf I: The Pinus taeda genome is characterized by diverse and highly diverged repetitive sequences. BMC Genomics. 2010, 11: 420-10.1186/1471-2164-11-420.PubMedPubMed CentralView ArticleGoogle Scholar
- Magbanua ZV, Ozkan S, Bartlett BD, Chouvarine P, Saski CA, Liston A, Cronn RC, Nelson CD: Peterson DG. Adventures in the enormous: A 1.8 million clone BAC library for the 21.7 Gb genome of loblolly pine. PLoS One. 2011, 6: e16214-Google Scholar
- Rocheta M, Cordeiro J, Oliveira M, Miguel C: PpRT 1: the first complete gypsy-like retrotransposon isolated in Pinus pinaster. Planta. 2006, 225: 551-562.PubMedView ArticleGoogle Scholar
- Liu WTS, Sehgal S, Chouvarine P, Peterson D: Characterization of the genome of bald cypress. BMC Genomics. 2011, 12: 553-10.1186/1471-2164-12-553.PubMedPubMed CentralView ArticleGoogle Scholar
- Kovach AWJ, Parra G, Holt C, Bruening GE, Loopstra CA, Hartigan J, Yandell M, Langley CH, Korf I, Neale DB: The Pinus taeda genome is characterized by diverse and highly diverged repetitive sequences. BMC Genomics. 2010, 11: 420-10.1186/1471-2164-11-420.PubMedPubMed CentralView ArticleGoogle Scholar
- Singleton TL: Levin. HL: A long terminal repeat retrotransposon of fission yeast has strong preferences for specific sites of insertion. Eukaryot Cell. 2002, 1: 44-55. 10.1128/EC.01.1.44-55.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Lynch C, Tristem M: A co-opted gypsy-type LTR-retrotransposon is conserved in the genomes of humans, sheep, mice, and rats. Curr Biol. 2003, 13: 1518-1523. 10.1016/S0960-9822(03)00618-3.PubMedView ArticleGoogle Scholar
- Eissenberg JC: Molecular biology of the chromo domain: an ancient chromatin module comes of age. Gene. 2001, 275: 19-29. 10.1016/S0378-1119(01)00628-X.PubMedView ArticleGoogle Scholar
- Konieczny A, Voytas DF, Cummingst MP, Ausubel FM: A superfamily of Arabidopsis thaliana retrotransposons. Genetics. 1991, 127: 801-809.PubMedPubMed CentralGoogle Scholar
- Hindmarsh P, Leis J: Retroviral DNA integration. Microbiol Mol Biol Rev. 1999, 63: 836-843.PubMedPubMed CentralGoogle Scholar
- McDonald JF, Matyunina LV, Wilson S, Jordan IK, Bowen NJ, Miller WJ: LTR retrotransposons and the evolution of eukaryotic enhancers. Genetica. 1997, 100: 3-13. 10.1023/A:1018392117410.PubMedView ArticleGoogle Scholar
- Atchison ML: Enhancers: mechanism of action and cell specificity. Annu Rev Cell Biol. 1998, 4: 127-153.View ArticleGoogle Scholar
- Vicient CM: Transcriptional activity of transposable elements in maize. BMC Genomics. 2010, 11: 601-10.1186/1471-2164-11-601.PubMedPubMed CentralView ArticleGoogle Scholar
- Tapia G, Verdugo I, Yañez M, Ahumada I, Theoduloz C, Cordero C, Poblete F, González E: Ruiz-Lara S. Involvement of ethylene in stress-induced expression of the TLC1.1 retrotransposon from Lycopersicon chilense Dun. Plant Physiol. 2005, 138: 2075-2086.Google Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Kim BR, Nam HY, Kim SU, Kim SI, Chang YJ: Normalization of reverse transcription quantitative-PCR with housekeeping genes in rice. Biotechnol Lett. 2003, 25: 1869-1872.PubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
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.