Full-length Long Terminal Repeat Research Articles

Characterisation of LTR-Retrotransposons of Stevia rebaudiana and Their Use for the Analysis of Genetic Variability

Stevia rebaudiana is one of the most important crops belonging to the Asteraceae family. Stevia is cultivated all over the world as it represents a valid natural alternative to artificial sweeteners thanks to its leaves, which produce steviol glycosides that have high sweetening power and reduced caloric value. In this work, the stevia genome sequence was used to isolate and characterise full-length long-terminal repeat retrotransposons (LTR-REs), which account for more than half of the genome. The Gypsy retrotransposons were twice as abundant as the Copia ones. A disproportionate abundance of elements belonging to the Chromovirus/Tekay lineage was observed among the Gypsy elements. Only the SIRE and Angela lineages represented significant portions of the genome among the Copia elements. The dynamics with which LTR-REs colonised the stevia genome were also estimated; all isolated full-length elements turned out to be relatively young, with a proliferation peak around 1–2 million years ago. However, a different analysis conducted by comparing sequences encoding retrotranscriptase showed the occurrence of an older period in which there was a lot of LTR-RE proliferation. Finally, a group of isolated full-length elements belonging to the lineage Angela was used to analyse the genetic variability in 25 accessions of S. rebaudiana using the Inter-Retrotransposon Amplified Polymorphism (IRAP) protocol. The obtained fingerprints highlighted a high degree of genetic variability and were used to study the genomic structures of the different accessions. It was hypothesised that there are four ancestral subpopulations at the root of the analysed accessions, which all turned out to be admixed. Overall, these data may be useful for genome sequence annotations and for evaluating genetic variability in this species, which may be useful in stevia breeding.

Reactivation of transposable elements following hybridization in fission yeast.

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this “genomic shock hypothesis” in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human-induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo LTRs before the split of the Sp and Sk lineages. Most full-length elements were inserted more recently, after hybridization. With the exception of a single full-length element with signs of positive selection, both solo LTRs and, in particular, full-length elements carry signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increases with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements, and provides evidence for the “genomic shock hypothesis.”

Correlation in Expression between LTR Retrotransposons and Potential Host Cis-Targets during Infection of Antherea pernyi with ApNPV Baculovirus.

The published genome sequence of Antheraea yamamai (Saturnnidae) was used to construct a library of long terminal repeat (LTR)-retrotransposons that is representative of the wild silkmoth (Antherea) genus, and that includes 22,666 solo LTRs and 541 full-length LTRs. The LTR retrotransposons of Antheraea yamamai (AyLTRs) could be classified into the three canonical groups of Gypsy, Copia and Belpao. Eleven AyLTRs contained the env gene element, but the relationship with the env element of baculovirus, particularly A. yamamai and pernyi nucleopolyhedrovirus (AyNPV and ApNPV), was distant. A total of 251 “independent” full-length AyLTRs were identified that were located within 100 kb distance (downstream or upstream) of 406 neighboring genes in A. yamamai. Regulation of these genes might occur in cis by the AyLTRs, and the neighboring genes were found to be enriched in GO terms such as “response to stimulus”, and KEGG terms such as “mTOR signaling pathway” among others. Furthermore, the library of LTR-retrotransposons and the A. yamamai genome were used to identify and analyze the expression of LTR-retrotransposons and genes in ApNPV-infected and non-infected A. pernyi larval midguts, using raw data of a published transcriptome study. Our analysis demonstrates that 93 full-length LTR-retrotransposons are transcribed in the midgut of A. pernyi of which 12 significantly change their expression after ApNPV infection (differentially expressed LTR-retrotransposons or DELs). In addition, the expression of differentially expressed genes (DEGs) and neighboring DELs on the chromosome following ApNPV infection suggests the possibility of regulation of expression of DEGs by DELs through a cis mechanism, which will require experimental verification. When examined in more detail, it was found that genes involved in Notch signaling and stress granule (SG) formation were significantly up-regulated in ApNPV-infected A. pernyi larval midgut. Moreover, several DEGs in the Notch and SG pathways were found to be located in the neighborhood of particular DELs, indicating the possibility of DEG-DEL cross-regulation in cis for these two pathways.

Bovine leukemia virus long terminal repeat variability: identification of single nucleotide polymorphisms in regulatory sequences

BackgroundLimited data are available on the incidence of variations in nucleotide sequences of long terminal repeat (LTR) regions of Bovine Leukemia Virus (BLV). Consequently, the possible impact of SNPs on BLV LTR function are poorly elucidated. Thus, a detailed and representative study of full-length LTR sequences obtained from sixty-four BLV isolates from different geographical regions of Poland, Moldova, Croatia, Ukraine and Russia were analyzed for their genetic variability.MethodsOverlap extension PCR, sequencing and Bayesian phylogenetic reconstruction of LTR sequences were performed. These analyses were followed by detailed sequence comparison, estimation of genetic heterogeneity and identification of transcription factor binding site (TFBS) modifications.ResultsPhylogenetic analysis of curated LTR sequences and those available in the GenBank database reflected the acknowledged env gene classification of BLV into 10 genotypes, and further clustered analysed sequences into three genotypes - G4, G7 and G8. Additional molecular studies revealed the presence of 97 point mutations distributed at 89 positions throughout all 64 LTR sequences. The highest rate of variability was noted in U3 and U5 subregions. However, the variability in regulatory sequences (VR) was assessed as lower than the variability within non-regulatory sequences (VNR) for both, U3 and U5 subregions. In contrast, VR value for R subregion, as well as for the total LTR, was higher than the VNR suggesting the existence of positive selection. Twelve unique SNPs for these LTR sequences localized in regulatory and non-regulatory elements were identified. The presence of different types of substitutions lead to the abrogation of present or to the creation of additional TFBS.ConclusionThis study represents the largest study of LTR genetic variability of BLV field isolates from Eastern part of Europe. Phylogenetic analysis of LTRs supports the clustering BLV variants based on their geographic origin. The SNP screening showed variations modifying LTR regulatory sequences, as well as altering TFBS. These features warrant further exploration as they could be related to proviral load and distinctive regulation of BLV transcription and replication.

The long terminal repeat negative control region is a critical element for insertional oncogenesis after gene transfer into hematopoietic progenitors with Moloney murine leukemia viral vectors

Integrating vectors based on γ-retroviruses and containing full-length long terminal repeats (LTRs) have been associated with activation of oncogene expression and leukemogenesis in human gene therapy trials. Identification of the specific molecular elements of the LTRs that have a role in insertional oncogenesis events is important as it can lead to the development of safer gene transfer vectors. The negative control region (NCR) of the LTR is a particularly well-conserved sequence among mammalian γ-retroviruses with demonstrated regulatory activity of gene transcription in hematopoietic cells, which led us to hypothesize that this region may have a role in insertional oncogenesis after γ-retroviral vector (GV)-mediated gene transfer into hematopoietic progenitors. We used an in vitro assay of murine bone marrow cell immortalization to compare the immortalization capabilities of a series of GVs carrying murine leukemia virus (MLV) LTR deletion mutants. Compared with GV carrying the full-length MLV LTR, deletion of the complete LTR enhancer sequence showed significant reduction of immortalization rates. However, the use of a mutant LTR deleted of the enhancer sequence, with exception of the NCR, did not affect immortalization. Importantly, the inclusion of an LTR mutant devoid only of the NCR did show significant reduction of immortalization rates compared with the full LTR sequence. Therefore, our data point to the NCR as a key element for immortalization and justify additional studies to evaluate its specific role in MLV-mediated insertional oncogenesis.

Genomic Landscape of Long Terminal Repeat Retrotransposons (LTR-RTs) and Solo LTRs as Shaped by Ectopic Recombination in Chicken and Zebra Finch.

Transposable elements (TEs) are nearly ubiquitous among eukaryotic genomes, but TE contents vary dramatically among phylogenetic lineages. Several mechanisms have been proposed as drivers of TE dynamics in genomes, including the fixation/loss of a particular TE insertion by selection or drift as well as structural changes in the genome due to mutation (e.g., recombination). In particular, recombination can have a significant and directional effect on the genomic TE landscape. For example, ectopic recombination removes internal regions of long terminal repeat retrotransposons (LTR-RTs) as well as one long terminal repeat (LTR), resulting in a solo LTR. In this study, we focus on the intra-species dynamics of LTR-RTs and solo LTRs in bird genomes. The distribution of LTR-RTs and solo LTRs in birds is intriguing because avian recombination rates vary widely within a given genome. We used published linkage maps and whole genome assemblies to study the relationship between recombination rates and LTR-removal events in the chicken and zebra finch. We hypothesized that regions with low recombination rates would harbor more full-length LTR-RTs (and fewer solo LTRs) than regions with high recombination rates. We tested this hypothesis by comparing the ratio of full-length LTR-RTs and solo LTRs across chromosomes, across non-overlapping megabase windows, and across physical features (i.e., centromeres and telomeres). The chicken data statistically supported the hypothesis that recombination rates are inversely correlated with the ratio of full-length to solo LTRs at both the chromosome level and in 1-Mb non-overlapping windows. We also found that the ratio of full-length to solo LTRs near chicken telomeres was significantly lower than those ratios near centromeres. Our results suggest a potential role of ectopic recombination in shaping the chicken LTR-RT genomic landscape.

Prediction of retrotransposons and assessment of genetic variability based on developed retrotransposon-based insertion polymorphism (RBIP) markers in Pyrus L.

Interspecific hybridization has been considered the major mode of evolution in Pyrus (pear), and thus, the genetic relationships within this genus have not been well documented. Retrotransposons are ubiquitous components of plant genomes and 42.4% of the pear genome was reported to be long terminal repeat (LTR) retrotransposons, implying that retrotransposons might be significant in the evolution of Pyrus. In this study, 1,836 putative full-length LTR retrotransposons were isolated and 196 retrotransposon-based insertion polymorphism (RBIP) primers were developed, of which 24 pairs to the Ppcr1 subfamily of copia retrotransposons were used to analyze genetic diversity among 110 Pyrus accessions from Eurasia. Our results showed that Ppcr1 replicated many times in the development of cultivated Asian pears. The genetic structure analysis and the unweighted pair group method with arithmetic mean (UPGMA) dendrogram indicated that all accessions could be divided into Oriental and Occidental groups. In Oriental pears, wild pea pears clustered separately into independent groups in accordance with their morphological classifications. Cultivars of P. ussuriensis Maxim, P. pyrifolia Nakai, and P. pyrifolia Chinese white pear were mingled together, which inferred that hybridization events occurred during the development of the cultivated Asian pears. In Occidental pears, two clades were obtained in the UPGMA dendrogram in accordance with their geographical distribution; one contained the European species and the other included species from North Africa and West Asia. New findings in this study will be important to further understand the phylogeny of Pyrus and origins of cultivated pears.

Quadruplex-forming sequences occupy discrete regions inside plant LTR retrotransposons

Retrotransposons with long terminal repeats (LTR) form a significant proportion of eukaryotic genomes, especially in plants. They have gag and pol genes and several regulatory regions necessary for transcription and reverse transcription. We searched for potential quadruplex-forming sequences (PQSs) and potential triplex-forming sequences (PTSs) in 18 377 full-length LTR retrotransposons collected from 21 plant species. We found that PQSs were often located in LTRs, both upstream and downstream of promoters from which the whole retrotransposon is transcribed. Upstream-located guanine PQSs were dominant in the minus DNA strand, whereas downstream-located guanine PQSs prevailed in the plus strand, indicating their role both at transcriptional and post-transcriptional levels. Our circular dichroism spectroscopy measurements confirmed that these PQSs readily adopted guanine quadruplex structures—some of them were paralell-stranded, while others were anti-parallel-stranded. The PQS often formed doublets at a mutual distance of up to 400 bp. PTSs were most abundant in 3′UTR (but were also present in 5′UTR). We discuss the potential role of quadruplexes and triplexes as the regulators of various processes participating in LTR retrotransposon life cycle and as potential recombination sites during post-insertional retrotransposon-based genome rearrangements.

Genome-Wide Survey and Comparative Analysis of LTR Retrotransposons and Their Captured Genes in Rice and Sorghum

Long terminal repeat (LTR) retrotransposons are the major class I mobile elements in plants. They play crucial roles in gene expansion, diversification and evolution. However, their captured genes are yet to be genome-widely identified and characterized in most of plants although many genomes have been completely sequenced. In this study, we have identified 7,043 and 23,915 full-length LTR retrotransposons in the rice and sorghum genomes, respectively. High percentages of rice full-length LTR retrotransposons were distributed near centromeric region in each of the chromosomes. In contrast, sorghum full-length LTR retrotransposons were not enriched in centromere regions. This dissimilarity could be due to the discrepant retrotransposition during and after divergence from their common ancestor thus might be contributing to species divergence. A total of 672 and 1,343 genes have been captured by these elements in rice and sorghum, respectively. Gene Ontology (GO) and gene set enrichment analysis (GSEA) showed that no over-represented GO term was identified in LTR captured rice genes. For LTR captured sorghum genes, GO terms with functions in DNA/RNA metabolism and chromatin organization were over-represented. Only 36% of LTR captured rice genes were expressed and expression divergence was estimated as 11.9%. Higher percentage of LTR captured rice genes have evolved into pseudogenes under neutral selection. On the contrary, higher percentage of LTR captured sorghum genes were under purifying selection and 72.4% of them were expressed. Thus, higher percentage of LTR captured sorghum genes was functional. Small RNA analysis suggested that some of LTR captured genes in rice and sorghum might have been involved in negative regulation. On the other hand, positive selection has been observed in both rice and sorghum LTR captured genes and some of them were still expressed and functional. The data suggest that some of these LTR captured genes might have evolved into new gene functions.

Comparison of Immortalization Potential of Gamma-Retroviral, Lentiviral and Foamy Virus Gene Transfer Vectors

Abstract 3116Therapeutic gene transfer has been used successfully to treat a variety of human genetic diseases. Although protocols have shown positive clinical outcomes, the success of clinical trials were tempered by adverse events, in which integration of the viral vectors increased transcription of cancer-related genes and thereby contributed to development of leukemia. In all documented cases of insertional mutagenesis, the viral vectors used contained full-length, gamma-retrovirus long terminal repeats (LTRs) that are able to exert strong promoter and enhancer activity driving not only the expression of the transgene carried by the vector, but also that of genes neighboring the insertion site. Assessing safety of integrating viral vectors for future clinical use is therefore of paramount importance.In preparation for gene therapy approaches for the Wiskott-Aldrich syndrome (WAS), we used an in vitro assay of murine bone marrow (BM) cell immortalization to compare the consequences of transduction by four different kinds of viral vectors, including a full-length LTR Moloney leukemia virus (MLV), a self-inactivating MLV carrying a 3' LTR deleted of the viral enhancer region (SIN), a lentivirus (LV), and a foamy virus (FV) construct. All vectors carried EGFP under the control of a ubiquitously acting chromatin opening element (UCOE) or the WAS endogenous promoter (WASp). In this assay, BM cells are harvested from C57BL6 mice, exposed to retroviral supernatants and cultured long-term. Derived lines are considered immortalized based on their ability to continue to grow in vitro for more than 6 weeks in the presence of interleukin-3 and stem cell factor.To date, MLV and SIN transduction of 123 and 132 cultures, respectively, gave rise to 48 and 43 immortalized lines (39.0% and 32.6% immortalization rate). The difference in immortalization rate between MLV and SIN vectors was not statistically significant (Chi square: p=0.30). As expected, immortalized cells were negative/low for IgER, cKit and Sca1 expression, and positive to different degrees for expression of the myeloid markers CD11b/Mac1 and Ly6g/Gr1. Transduction of 114 and 62 cultures with LV and FV vectors, respectively, resulted in no immortalized lines.Real-time PCR was performed to evaluate transduction efficiency of bone marrow cells and immortalized lines. Integrated vector sequences in bone marrow cells transduced by LV and FV were detected in significantly higher quantities than in cells transduced with MLV and SIN vectors. However, expression of the EGFP transgene was markedly reduced in LV- and FV-transduced cells compared to cells exposed to MLV vectors (MFI: 14.0, 1.88, 93.2, respectively).These preliminary results confirm that gamma-retroviral gene transfer vectors are prone to causing immortalization of hematopoietic cells and suggest the vectors based on LV and FV backbones may be safer alternatives for WAS and other genetic disorders, provided that effective gene expression levels can be achieved in the biological model system under study.TableSummary of immortalization results using MLV, SIN, LV and FV vectorsVirusMOITransduction efficiencyTransduction experiments% ImmortalizationMLV/UCOE/EGFP2067%5643MLV/WASp/EGFP2065%6736SIN/UCOE/EGFP532%7236SIN/WASp/EGFP537%6028LV/UCOE/EGFP1040%570LV/WASp/EGFP1049%570FV/UCOE/EGFP1027%280FV/WASp/EGFP1022%340negativeN.A.N.A.720 Disclosures:No relevant conflicts of interest to declare.

A computational study of the dynamics of LTR retrotransposons in the Populus trichocarpa genome

Retrotransposons are an ubiquitous component of plant genomes, especially abundant in species with large genomes. Populus trichocarpa has a relatively small genome, which was entirely sequenced; however, studies focused on poplar retrotransposons dynamics are rare. With the aim to study the retrotransposon component of the poplar genome, we have scanned the complete genome sequence searching full-length long-terminal repeat (LTR) retrotransposons, i.e., characterised by two long terminal repeats at the 5′ and 3′ ends. A computational approach based on detection of conserved structural features, on building multiple alignments, and on similarity searches was used to identify 1,479 putative full-length LTR retrotransposons. Ty1-copia elements were more numerous than Ty3-gypsy. However, many LTR retroelements were not assigned to any superfamily because lacking of diagnostic features and non-autonomous. LTR retrotransposon remnants were by far more numerous than full-length elements, indicating that during the evolution of poplar, large amplification of these elements was followed by DNA loss. Within superfamilies, Ty3-gypsy families are made of more members than Ty1-copia ones. Retrotransposition occurred with increasing frequency following the separation of Populus sections, with different waves of retrotransposition activity between Ty3-gypsy and Ty1-copia elements. Recently inserted elements appear more frequently expressed than older ones. Finally, different levels of activity of retrotransposons were observed according to their position and their density in the linkage groups. On the whole, the results support the view of retrotransposons as a community of different organisms in the genome, whose activity (both retrotransposition and DNA loss) has heavily impacted and probably continues to impact poplar genome structure and size.

Structure and evolution of full-length LTR retrotransposons in rice genome

The long terminal repeat (LTR) retrotransposons are the most abundant class of transposable elements in plant genomes and play important roles in genome divergence and evolution. Their accumulation is the main factor influencing genome size increase in plants. Rice (Oryza sativa L.) is a model monocot and is the focus of much biological research due to its economic importance. We conducted a comprehensive survey of full-length LTR retrotransposons based on the completed genome of japonica rice variety Nipponbare (TIGR Release 5), with the newly published tool LTR-FINDER. The elements could be categorized into 29 structural domain categories (SDCs), and their total copy number identified was estimated at >6,000. Most of them were relatively young: more than 90% were less than 10 My. There existed a high level of activity among them as a whole at 0–1 Mya, but different categories possessed distinct amplification patterns. Most recently inserted elements were specific to the rice genome, while a few were conserved across species. This study provides new insights into the structure and evolutionary history of the full-length retroelements in the rice genome.

Large-scale transcriptome data reveals transcriptional activity of fission yeast LTR retrotransposons

BackgroundRetrotransposons are transposable elements that proliferate within eukaryotic genomes through a process involving reverse transcription. The numbers of retrotransposons within genomes and differences between closely related species may yield insight into the evolutionary history of the elements. Less is known about the ongoing dynamics of retrotransposons, as analysis of genome sequences will only reveal insertions of retrotransposons that are fixed - or near fixation - in the population or strain from which genetic material has been extracted for sequencing. One pre-requisite for retrotransposition is transcription of the elements. Given their intrinsic sequence redundancy, transcriptome-level analyses of transposable elements are scarce. We have used recently published transcriptome data from the fission yeast Schizosaccharomyces pombe to assess the ability to detect and describe transcriptional activity from Long Terminal Repeat (LTR) retrotransposons. LTR retrotransposons are normally flanked by two LTR sequences. However, the majority of LTR sequences in S. pombe exist as solitary LTRs, i.e. as single terminal repeat sequences not flanking a retrotransposon. Transcriptional activity was analysed for both full-length LTR retrotransposons and solitary LTRs.ResultsTwo independent sets of transcriptome data reveal the presence of full-length, polyadenylated transcripts from LTR retrotransposons in S. pombe during growth phase in rich medium. The redundancy of retrotransposon sequences makes it difficult to assess which elements are transcriptionally active, but data strongly indicates that only a subset of the LTR retrotransposons contribute significantly to the detected transcription. A considerable level of reverse strand transcription is also detected. Equal levels of transcriptional activity are observed from both strands of solitary LTR sequences. Transcriptome data collected during meiosis suggests that transcription of solitary LTRs is correlated with the transcription of nearby protein-coding genes.ConclusionsPresumably, the host organism negatively regulates proliferation of LTR retrotransposons. The finding of considerable transcriptional activity of retrotransposons suggests that part of this regulation is likely to take place at a post-transcriptional level. Alternatively, the transcriptional activity may signify a hitherto unrecognized activity level of retrotransposon proliferation. Our findings underline the usefulness of transcriptome data in elucidating dynamics in retrotransposon transcription.

LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons

BackgroundTransposable elements are abundant in eukaryotic genomes and it is believed that they have a significant impact on the evolution of gene and chromosome structure. While there are several completed eukaryotic genome projects, there are only few high quality genome wide annotations of transposable elements. Therefore, there is a considerable demand for computational identification of transposable elements. LTR retrotransposons, an important subclass of transposable elements, are well suited for computational identification, as they contain long terminal repeats (LTRs).ResultsWe have developed a software tool LTRharvest for the de novo detection of full length LTR retrotransposons in large sequence sets. LTRharvest efficiently delivers high quality annotations based on known LTR transposon features like length, distance, and sequence motifs. A quality validation of LTRharvest against a gold standard annotation for Saccharomyces cerevisae and Drosophila melanogaster shows a sensitivity of up to 90% and 97% and specificity of 100% and 72%, respectively. This is comparable or slightly better than annotations for previous software tools. The main advantage of LTRharvest over previous tools is (a) its ability to efficiently handle large datasets from finished or unfinished genome projects, (b) its flexibility in incorporating known sequence features into the prediction, and (c) its availability as an open source software.ConclusionLTRharvest is an efficient software tool delivering high quality annotation of LTR retrotransposons. It can, for example, process the largest human chromosome in approx. 8 minutes on a Linux PC with 4 GB of memory. Its flexibility and small space and run-time requirements makes LTRharvest a very competitive candidate for future LTR retrotransposon annotation projects. Moreover, the structured design and implementation and the availability as open source provides an excellent base for incorporating novel concepts to further improve prediction of LTR retrotransposons.

De novo identification of LTR retrotransposons in eukaryotic genomes

BackgroundLTR retrotransposons are a class of mobile genetic elements containing two similar long terminal repeats (LTRs). Currently, LTR retrotransposons are annotated in eukaryotic genomes mainly through the conventional homology searching approach. Hence, it is limited to annotating known elements.ResultsIn this paper, we report a de novo computational method that can identify new LTR retrotransposons without relying on a library of known elements. Specifically, our method identifies intact LTR retrotransposons by using an approximate string matching technique and protein domain analysis. In addition, it identifies partially deleted or solo LTRs using profile Hidden Markov Models (pHMMs). As a result, this method can de novo identify all types of LTR retrotransposons. We tested this method on the two pairs of eukaryotic genomes, C. elegans vs. C. briggsae and D. melanogaster vs. D. pseudoobscura. LTR retrotransposons in C. elegans and D. melanogaster have been intensively studied using conventional annotation methods. Comparing with previous work, we identified new intact LTR retroelements and new putative families, which may imply that there may still be new retroelements that are left to be discovered even in well-studied organisms. To assess the sensitivity and accuracy of our method, we compared our results with a previously published method, LTR_STRUC, which predominantly identifies full-length LTR retrotransposons. In summary, both methods identified comparable number of intact LTR retroelements. But our method can identify nearly all known elements in C. elegans, while LTR_STRUCT missed about 1/3 of them. Our method also identified more known LTR retroelements than LTR_STRUCT in the D. melanogaster genome. We also identified some LTR retroelements in the other two genomes, C. briggsae and D. pseudoobscura, which have not been completely finished. In contrast, the conventional method failed to identify those elements. Finally, the phylogenetic and chromosomal distributions of the identified elements are discussed.ConclusionWe report a novel method for de novo identification of LTR retrotransposons in eukaryotic genomes with favorable performance over the existing methods.

Characterization of AFLAV, a Tf1/Sushi retrotransposon from Aspergillus flavus

The plasmid, pAF28, a genomic clone from Aspergillus flavus NRRL 6541, has been used as a hybridization probe to fingerprint A. flavus strains isolated in corn and peanut fields. The insert of pAF28 contains a 4.5 kb region which encodes a truncated retrotransposon (AfRTL-1). In search for a full-length and intact copy of retrotransposon, we exploited a novel PCR cloning strategy by amplifying a 3.4 kb region from the genomic DNA of A. flavus NRRL 6541. The fragment was cloned into pCR 4-TOPO. Sequence analysis confirmed that this region encoded putative domains of partial reverse transcriptase, RNase H, and integrase of the predicted retrotransposon. The two flanking long terminal repeats (LTRs) and the sequence between them comprise a putative full-length LTR retrotransposon of 7799 bp in length. This intact retrotransposon sequence is named AFLAV (A. flavus Retrotransposon). The order of the predicted catalytic domains in the polyprotein (Pol) placed AFLAV in the Tf1/sushi subgroup of the Ty3/gypsy retrotransposon family. Primers derived from AFLAV sequence were used to screen this retrotransposon in other strains of A. flavus. More than fifty strains of A. flavus isolated from different geological origins were surveyed and the results show that many strains have extensive deletions in the regions encoding the capsid (Gag) and Pol.

CCAAT/Enhancer-binding Protein-mediated Inhibition of HTLV-1 Viral Gene Expression

Recently, CCAAT/enhancer-binding protein (C/EBP) transcription factors have been shown to form heterodimers with cAMP-responsive element binding protein 2 (CREB2), another bZIP transcriptional regulator known to play a critical role in driving basal and Tax-mediated transactivation of the HTLV-1 long terminal repeat (LTR). Herein, we describe that overexpression of C/EBPa and C/EBPb, including the endogenous isoforms liver-enriched activation protein (LAP) and liver-enriched inhibitory protein (LIP) inhibits Tax-mediated transactivation of the HTLV-1 LTR. C/EBP-mediated inhibition was not the result of competition with the viral oncoprotein Tax for recruitment of co-activators CREB-binding protein (CBP)/p300 to the viral promoter. Electrophoretic mobility shift analyses demonstrated that C/EBP proteins derived from U937 monocytic nuclear extracts directly interacted with the Tax-responsive element 1 repeat III as well as Taxresponsive element 2. However, deletion of these sequences within the context of the full-length LTR did not prevent C/EBP-mediated inhibition. Disruption of C/ EBPb/CREB-2 heterodimerization by deletion of the C/ EBPb leucine zipper prevented C/EBP inhibition of Taxmediated transactivation of the LTR. These results suggest that C/EBP binding to Tax was not a requirement for C/ EBP-mediated inhibition. from 2005 International Meeting of The Institute of Human Virology Baltimore, USA, 29 August – 2 September 2005

Evolutionary history of Oryza sativa LTR retrotransposons: a preliminary survey of the rice genome sequences.

BackgroundLTR Retrotransposons transpose through reverse transcription of an RNA intermediate and are ubiquitous components of all eukaryotic genomes thus far examined. Plant genomes, in particular, have been found to be comprised of a remarkably high number of LTR retrotransposons. There is a significant body of direct and indirect evidence that LTR retrotransposons have contributed to gene and genome evolution in plants.ResultsTo explore the evolutionary history of long terminal repeat (LTR) retrotransposons and their impact on the genome of Oryza sativa, we have extended an earlier computer-based survey to include all identifiable full-length, fragmented and solo LTR elements in the rice genome database as of April 2002. A total of 1,219 retroelement sequences were identified, including 217 full-length elements, 822 fragmented elements, and 180 solo LTRs. In order to gain insight into the chromosomal distribution of LTR-retrotransposons in the rice genome, a detailed examination of LTR-retrotransposon sequences on Chromosome 10 was carried out. An average of 22.3 LTR-retrotransposons per Mb were detected in Chromosome 10.ConclusionsGypsy-like elements were found to be >4 × more abundant than copia-like elements. Eleven of the thirty-eight investigated LTR-retrotransposon families displayed significant subfamily structure. We estimate that at least 46.5% of LTR-retrotransposons in the rice genome are older than the age of the species (< 680,000 years). LTR-retrotransposons present in the rice genome range in age from those just recently inserted up to nearly 10 million years old. Approximately 20% of LTR retrotransposon sequences lie within putative genes. The distribution of elements across chromosome 10 is non-random with the highest density (48 elements per Mb) being present in the pericentric region.

CDP binding to multiple sites in the mouse mammary tumor virus long terminal repeat suppresses basal and glucocorticoid-induced transcription.

Mouse mammary tumor virus (MMTV) is transcribed at high levels in the lactating mammary gland to ensure transmission of virus from the milk of infected female mice to susceptible offspring. We previously have shown that the transcription factor CCAAT displacement protein (CDP) is expressed in high amounts in virgin mammary gland, yet DNA-binding activity for the MMTV long terminal repeat (LTR) disappears as mammary tissue differentiates during lactation. CDP is a repressor of MMTV expression and, therefore, MMTV expression is suppressed during early mammary gland development. In this study, we have shown using DNase I footprinting and electrophoretic mobility shift assays that there are at least five CDP-binding sites in the MMTV LTR upstream of those previously described in the promoter-proximal negative regulatory element (NRE). Single mutations in two of these upstream sites (+691 or +692 and +735 relative to the first base of the LTR) reduced CDP binding to the cognate sites and elevated reporter gene expression from the full-length MMTV LTR. Combination of a mutation in the promoter-distal NRE with a mutation in the proximal NRE gave approximately additive increases in LTR-reporter gene activity, suggesting that these binding sites act independently. Mutations in several different CDP-binding sites allowed elevation of reporter gene activity from the MMTV promoter in the absence and presence of glucocorticoids, hormones that contribute to high levels of MMTV transcription during lactation by activation of hormone receptor binding to the LTR. In addition, overexpression of CDP in transient-transfection assays suppressed both basal and glucocorticoid-induced LTR-mediated transcription in a dose-dependent manner. These data suggest that multiple CDP-binding sites contribute independently to regulate binding of positive factors, including glucocorticoid receptor, to the MMTV LTR during mammary gland development.

Characterization of Interactions between Transcription Factors and a Regulatory Region Spanning nt –320 to –281 of the HIV-1 LTR in T-Lymphoid and Non-T-Lymphoid Cells

HIV-1 gene expression is regulated by the interplay of transcription factors with multiple binding motifs present within the U3, R and U5 regions of the long terminal repeat (LTR). Here we report novel DNA binding complexes (termed 9a, 9b and 9c) between nuclear proteins from T-lymphoid and non-T-lymphoid cells and a region of the U3 LTR between nucleotides (nts) –320 to –281 in the HIV strain HXB2. Complex 9b bound a motif predicted to bind E-box or c-Myb proteins and a partially overlapping dyad symmetrical motif, and included basic helix-loop-helix proteins (E12, E47 or ITF-1) but surprisingly not c-Myb. Complex 9c, which bound to a pair of GATA sites, included GATA-3 and GATA-4 in Jurkat and MT-2 cells, respectively. We also demonstrate that the c-Myb/E-box and GATA sites form a bipartite motif required for the formation of complex 9a. Transient transfection experiments with T cells revealed that in the context of a minichromosome assembled full-length LTR, mutation of region –320 to –281 increased basal and PMA-stimulated LTR activity. These findings suggest that this region is an important component of the HIV-1 LTR required for response to different cellular transcription factors.