Large distribution and high sequence identity of a Copia-type retrotransposon in angiosperm families

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A novel linkage map of sugarcane with evidence for clustering of retrotransposon-based markers

Background: The development of sugarcane as a sustainable crop has unlimited applications. The crop is one of the most economically viable for renewable energy production, and CO2 balance. Linkage maps are valuable tools for understanding genetic and genomic organization, particularly in sugarcane due to its complex polyploid genome of multispecific origins. The overall objective of our study was to construct a novel sugarcane linkage map, compiling AFLP and EST-SSR markers, and to generate data on the distribution of markers anchored to sequences of scIvana_1, a complete sugarcane transposable element, and member of the Copia superfamily. Results: The mapping population parents ('IAC66-6' and 'TUC71-7') contributed equally to polymorphisms, independent of marker type, and generated markers that were distributed into nearly the same number of co-segregation groups (or CGs). Bi-parentally inherited alleles provided the integration of 19 CGs. The marker number per CG ranged from two to 39. The total map length was 4,843.19 cM, with a marker density of 8.87 cM. Markers were assembled into 92 CGs that ranged in length from 1.14 to 404.72 cM, with an estimated average length of 52.64 cM. The greatest distance between two adjacent markers was 48.25 cM. The scIvana_1-based markers (56) were positioned on 21 CGs, but were not regularly distributed. Interestingly, the distance between adjacent scIvana_1-based markers was less than 5 cM, and was observed on five CGs, suggesting a clustered organization. Conclusions: Results indicated the use of a NBS-profiling technique was efficient to develop retrotransposon-based markers in sugarcane. The simultaneous maximum-likelihood estimates of linkage and linkage phase based strategies confirmed the suitability of its approach to estimate linkage, and construct the linkage map. Interestingly, using our genetic data it was possible to calculate the number of retrotransposonscIvana_1 (similar to 60) copies in the sugarcane genome, confirming previously reported molecular results. In addition, this research possibly will have indirect implications in crop economics e. g., productivity enhancement via QTL studies, as the mapping population parents differ in response to an important fungal disease.

Radiation of the Tnt1 retrotransposon superfamily in three Solanaceae genera

Abstract Background Tnt1 was the first active plant retrotransposon identified in tobacco after nitrate reductase gene disruption. The Tnt1 superfamily comprises elements from Nicotiana (Tnt1 and Tto1) and Lycopersicon (Retrolyc1 and Tlc1) species. The study presented here was conducted to characterise Tnt1-related sequences in 20 wild species of Solanum and five cultivars of Solanum tuberosum. Results Tnt1-related sequences were amplified from total genomic DNA using a PCR-based approach. Purified fragments were cloned and sequenced, and clustering analysis revealed three groups that differ in their U3 region. Using a network approach with a total of 453 non-redundant sequences isolated from Solanum (197), Nicotiana (140) and Lycopersicon (116) species, it is demonstrated that the Tnt1 superfamily can be treated as a population to resolve previous phylogenetic multifurcations. The resulting RNAseH network revealed that sequences group according to the Solanaceae genus, supporting a strong association with the host genome, whereas tracing the U3 region sequence association characterises the modular evolutionary pattern within the Tnt1 superfamily. Within each genus, and irrespective of species, nearly 20% of Tnt1 sequences analysed are identical, indicative of being part of an active copy. The network approach enabled the identification of putative "master" sequences and provided evidence that within a genus these master sequences are associated with distinct U3 regions. Conclusion The results presented here support the hypothesis that the Tnt1 superfamily was present early in the evolution of Solanaceae. The evidence also suggests that the RNAseH region of Tnt1 became fixed at the host genus level whereas, within each genus, propagation was ensured by the diversification of the U3 region. Different selection pressures seemed to have acted on the U3 and RNAseH modules of ancestral Tnt1 elements, probably due to the distinct functions of these regions in the retrotransposon life cycle, resulting in both co evolution and adaptation of the element population with its host.

Transposable elements: structure and dynamic of the autonomous retrotransposon R2 in Arthropoda with non-canonical reproduction

Eukaryotic ribosomal DNA constitutes a multi gene family organized in a cluster called nucleolar organizer region (NOR); this region is composed usually by hundreds to thousands of tandemly repeated units. Ribosomal genes, being repeated sequences, evolve following the typical pattern of concerted evolution. The autonomous retroelement R2 inserts in the ribosomal gene 28S, leading to defective 28S rDNA genes. R2 element, being a retrotransposon, performs its activity in the genome multiplying its copy number through a “copy and paste” mechanism called target primed reverse transcription. It consists in the retrotranscription of the element’s mRNA into DNA, then the DNA is integrated in the target site. Since the retrotranscription can be interrupted, but the integration will be carried out anyway, truncated copies of the element will also be present in the genome. The study of these truncated variants is a tool to examine the activity of the element. R2 phylogeny appears, in general, not consistent with that of its hosts, except some cases (e.g. Drosophila spp. and Reticulitermes spp.); moreover R2 is absent in some species (Fugu rubripes, human, mouse, etc.), while other species have more R2 lineages in their genome (the turtle Mauremys reevesii, the Japanese beetle Popilia japonica, etc). R2 elements here presented are isolated in 4 species of notostracan branchiopods and in two species of stick insects, whose reproductive strategies range from strict gonochorism to unisexuality. From sequencing data emerges that in Triops cancriformis (Spanish gonochoric population), in Lepidurus arcticus (two putatively unisexual populations from Iceland) and in Bacillus rossius (gonochoric population from Capalbio) the R2 elements are complete and encode functional proteins, reflecting the general features of this family of transposable elements. On the other hand, R2 from Italian and Austrian populations of T. cancriformis (respectively unisexual and hermaphroditic), Lepidurus lubbocki (two elements within the same Italian population, gonochoric but with unfunctional males) and Bacillus grandii grandii (gonochoric population from Ponte Manghisi) have sequences that encode incomplete or non-functional proteins in which it is possible to recognize only part of the characteristic domains. In Lepidurus couesii (Italian gonochoric populations) different elements were found as in L. lubbocki, and the sequencing is still in progress. Two hypothesis are given to explain the inconsistency of R2/host phylogeny: vertical inheritance of the element followed by extinction/diversification or horizontal transmission. My data support previous study that state the vertical transmission as the most likely explanation; nevertheless horizontal transfer events can’t be excluded. I also studied the element’s activity in Spanish populations of T. cancriformis, in L. lubbocki, in L. arcticus and in gonochoric and parthenogenetic populations of B. rossius. In gonochoric populations of T. cancriformis and B. rossius I found that each individual has its own private set of truncated variants. The situation is the opposite for the remaining hermaphroditic/parthenogenetic species and populations, all individuals sharing – in the so far analyzed samples - the majority of variants. This situation is very interesting, because it isn’t concordant with the Muller’s ratchet theory that hypothesizes the parthenogenetic populations being either devoided of transposable elements or TEs overloaded. My data suggest a possible epigenetic mechanism that can block the retrotransposon activity, and in this way deleterious mutations don’t accumulate.

Effectiveness of AFLPs and retrotransposon-based markers for the identification of portuguese grapevine cultivars and clones

Grapevine germplasm, including 38 of the main Portuguese cultivars and three foreign cultivars, Pinot Noir, Pinot Blanc and Chasselas, used as a reference, and 37 true-to-type clones from the Alvarinho, Arinto, Loureiro, Moscatel Galego Branco, Trajadura and Vinhão cultivars were studied using AFLP and three retrotransposon-based molecular techniques, IRAP, REMAP and SSAP. To study the retrotransposon-based polymorphisms, 18 primers based on the LTR sequences of Tvv1, Gret1 and Vine-1 were used. In the analysis of 41 cultivars, 517 IRAP, REMAP, AFLP and SSAP fragments were obtained, 83% of which were polymorphic. For IRAP, only the Tvv1Fa primer amplified DNA fragments. In the REMAP analysis, the Tvv1Fa-Ms14 primer combination only produced polymorphic bands, and the Vine-1 primers produced mainly ISSR fragments. The highest number of polymorphic fragments was found for AFLP. Both AFLP and SSAP showed a greater capacity for identifying clones, resulting in 15 and 9 clones identified, respectively. Together, all of the techniques allowed for the identification of 54% of the studied clones, which is an important step in solving one of the challenges that viticulture currently faces.

Evidence for the recent horizontal transfer of long terminal repeat retrotransposon

The evolutionary dynamics existing between transposable elements (TEs) and their host genomes have been likened to an “arms race.” The selfish drive of TEs to replicate, in turn, elicits the evolution of host-mediated regulatory mechanisms aimed at repressing transpositional activity. It has been postulated that horizontal (cross-species) transfer may be one effective strategy by which TEs and other selfish genes can escape host-mediated silencing mechanisms over evolutionary time; however, to date, the most definitive evidence that TEs horizontally transfer between species has been limited to class II or DNA-type elements. Evidence that the more numerous and widely distributed retroelements may also be horizontally transferred between species has been more ambiguous. In this paper, we report definitive evidence for a recent horizontal transfer of the copia long terminal repeat retrotransposon between Drosophila melanogaster and Drosophila willistoni.

Silent chromatin determines target preference of the Saccharomyces retrotransposon Ty5

The HML and HMR mating loci of Saccharomyces cerevisiae are bound in silent chromatin, which is assembled at the flanking E and I transcriptional silencers. The retrotransposon Ty5 preferentially integrates into regions of silent chromatin, and Ty5 insertions near the HMR-E silencer account for ≈2% of total transposition events. Most Ty5 insertions occur within 800 bp on either side of the autonomously replicating consensus sequence within HMR-E. Ty5 target preference is determined by silent chromatin, because integration near HMR-E is abolished in strains with silencer mutations that alleviate transcriptional repression. The recognition of specific DNA sequences per se does not direct integration, rather, it is the protein complex assembled at the silencers. As demonstrated here for Ty5, recognition of specific chromatin domains may be a general mechanism by which retrotransposons and retroviruses determine integration sites.

Methylation is not the main force repressing the retrotransposon MAGGY in Magnaporthe grisea

We have introduced the LTR-retrotransposon MAGGY into a naive genome of Magnaporthe grisea and estimated the copy number of MAGGY in a cell by serial isolation of fungal protoplasts at certain time intervals. The number of MAGGY elements rapidly increased for a short period following introduction. However, it did not increase geometrically and reached equilibrium at 20–30 copies per genome, indicating that MAGGY was repressed or silenced during proliferation. De novo methylation of MAGGY occurred immediately following invasion into the genome but the degree of methylation was constant and did not correlate with the repression of MAGGY. 5-Azacytidine treatment demethylated and transcriptionally activated the MAGGY element in regenerants but did not affect transpositional frequency, suggesting that post-transcriptional suppression, not methylation, is the main force that represses MAGGY proliferation in M.grisea. Support for this conclusion was also obtained by examining the methylation status of MAGGY sequences in field isolates of M.grisea with active or inactive MAGGY elements. Methylation of the MAGGY sequences was detected in some isolates but not in others. However, the methylation status did not correlate with the copy numbers and activity of the elements.

Retrotransposon insertion induces an isozyme of sn-glycerol-3-phosphate dehydrogenase in Drosophila melanogaster.

The insertion of the blood retrotransposon into the untranslated region of exon 7 of the sn-glycerol-3-phosphate dehydrogenase-encoding gene (Gpdh) in Drosophila melanogaster induces a GPDH isozyme-GPDH-4-and alters the pattern of expression of the three normal isozymes-GPDH-1 to GPDH-3. The process of transcript terminus formation inside the retrotransposon insertion reduces the level of the Gpdh transcript that contains exon 8 and increases the level of the transcript that contains exons 1-7. The induced GPDH-4 isozyme is a translation product of the three transcripts that contain fragments of the blood retrotransposon. The mechanism of mutagenesis by the blood insertion is postulated to involve the pause or termination of transcription within the blood sequence, which in turn is caused by the interference of a DNA-binding protein with the RNA polymerase. Thus, we show the formation of a new functional GPDH protein by the insertion of a transposable element and discuss the evolutionary significance of this phenomenon.

Retrotransposon expression profiling: Unveiling the hidden SINE transcriptome through Next-Generation Sequencing data analysis

Of the ~1.7 million SINE elements in the human genome, only a tiny number are estimated to be active in transcription by RNA polymerase (Pol) III. Tracing the individual loci from which SINE transcripts originate is complicated by their highly repetitive nature. By exploiting RNA-Seq datasets and unique SINE DNA sequences, we devised a bioinformatic pipeline allowing us to identify Pol III-dependent transcripts of individual SINE elements. When applied to ENCODE transcriptomes of seven human cell lines, this search strategy identified ~1300 Alu loci and ~1100 MIR loci corresponding to detectable transcripts, with ~120 and ~60 respectively Alu and MIR loci expressed in at least three cell lines. In vitro transcription of selected SINEs did not reflect their in vivo expression properties, and required the native 5’-flanking region in addition to internal promoter. We also identified a cluster of expressed AluYa5-derived transcription units, juxtaposed to snaR genes on chromosome 19, formed by a promoter-containing left monomer fused to an Alu-unrelated downstream moiety. Autonomous Pol III transcription was also revealed for SINEs nested within Pol II-transcribed genes raising the possibility of an underlying mechanism for Pol II gene regulation by SINE transcriptional units. Moreover the application of our bioinformatic pipeline to both RNA-seq data of cells subjected to an in vitro pro-oncogenic stimulus and of in vivo matched tumor and non-tumor samples allowed us to detect increased Alu RNA expression as well as the source loci of such deregulation. The ability to investigate SINE transcriptomes at single-locus resolution will facilitate both the identification of novel biologically relevant SINE RNAs and the assessment of SINE expression alteration under pathological conditions.

Arrangement of gene pairs, retrotransposon insertions, and regulation of gene expression in plants

Plant genomes are extremely complex. Myriad factors contribute to their evolution and organization, as well as to the expression and regulation of individual genes. Here we present investigations into several such factors and their influence on genome structure and gene expression: the arrangement of pairs of physically adjacent genes, retrotransposons closely associated with genes, and the effect of retrotransposons on gene pair evolution. All sequenced plant genomes contain a significant fraction of retrotransposons, including that of rice. We investigated the effects of retrotransposons within rice genes and within a 1 kb putative promoter region upstream of each gene. We found that approximately one-sixth of all rice genes are closely associated with retrotransposons. Insertions within a gene’s promoter region tend to block gene expression, while retrotransposons within genes promote the existence of alternative splicing forms. We also identified several other trends in retrotransposon insertion and its effects on gene expression. Several studies have previously noted a connection among genes between physical proximity and correlated expression profiles. To determine the degree to which this correlation depends on an exact physical arrangement, we studied the expression and interspecies conservation of convergent and divergent gene pairs in rice, Arabidopsis, and Populus trichocarpa. Correlated expression among gene pairs was quite common in all three species, yet conserved arrangement was rare. However, conservation of gene pair arrangement was significantly more common among pairs with strongly correlated expression levels. In order to uncover additional properties of gene pair conservation and rearrangement, we performed a comparative analysis of convergent, divergent, and tandem gene pairs in rice, sorghum, maize, and Brachypodium. We noted considerable differences between gene pair types and species. We also constructed a putative evolutionary history for each pair, which led to several interesting discoveries. To further elucidate the causes of gene pair conservation and rearrangement, we identified retrotransposon insertions in and near rice gene pairs. Retrotransposon-associated pairs are less likely to be conserved, although there are significant differences in the possible effect of different types and locations of retrotransposon insertions. The three types of gene pair also varied in their susceptibility to retrotransposon-associated evolutionary changes.

Apple skin patterning is associated with differential expression of MYb10

Background Some apple (Malus × domestica Borkh.) varieties have attractive striping patterns, a quality attribute that is important for determining apple fruit market acceptance. Most apple cultivars (e.g. 'Royal Gala') produce fruit with a defined fruit pigment pattern, but in the case of 'Honeycrisp' apple, trees can produce fruits of two different kinds: striped and blushed. The causes of this phenomenon are unknown. Results Here we show that striped areas of 'Honeycrisp' and 'Royal Gala' are due to sectorial increases in anthocyanin concentration. Transcript levels of the major biosynthetic genes and MYB10, a transcription factor that upregulates apple anthocyanin production, correlated with increased anthocyanin concentration in stripes. However, nucleotide changes in the promoter and coding sequence of MYB10 do not correlate with skin pattern in 'Honeycrisp' and other cultivars differing in peel pigmentation patterns. A survey of methylation levels throughout the coding region of MYB10 and a 2.5 Kb region 5' of the ATG translation start site indicated that an area 900 bp long, starting 1400 bp upstream of the translation start site, is highly methylated. Cytosine methylation was present in all three contexts, with higher methylation levels observed for CHH and CHG (where H is A, C or T) than for CG. Comparisons of methylation levels of the MYB10 promoter in 'Honeycrisp' red and green stripes indicated that they correlate with peel phenotypes, with an enrichment of methylation observed in green stripes. Conclusions Differences in anthocyanin levels between red and green stripes can be explained by differential transcript accumulation of MYB10. Different levels of MYB10 transcript in red versus green stripes are inversely associated with methylation levels in the promoter region. Although observed methylation differences are modest, trends are consistent across years and differences are statistically significant. Methylation may be associated with the presence of a TRIM retrotransposon within the promoter region, but the presence of the TRIM element alone cannot explain the phenotypic variability observed in 'Honeycrisp'. We suggest that methylation in the MYB10 promoter is more variable in 'Honeycrisp' than in 'Royal Gala', leading to more variable color patterns in the peel of this cultivar.

A copia-like element in Pisum demonstrates the uses of dispersed repeated sequences in genetic analysis

A DNA sequence between two legumin genes in Pisum is a member of the copia-like class of retrotransposons and represents one member of a polymorphic and heterogeneous dispersed repeated sequence family in Pisum. This sequence can be exploited in genetic studies either by RFLP analysis where several markers can be scored together, or the segregation of individual elements can be followed after PCR amplification of specific members.