Helitron and CACTA DNA transposons actively reshape the common wheat - AK58 genome

Transposable elements (TEs) play a pivotal role in reshaping the plant genome. Helitrons represent a new class of transposable elements recently discovered in animals and plants. Helitrons, DNA transposons that replicate via a rolling-circle replication mechanism, are a major driving force behind genome evolution. Since the recent divergence of the modern cultivars (e.g., AK58) and landraces (e.g., Chinese Spring), Helitrons appear to have contributed greatly to genome variability. We first identified 214,665 Helitrons in AK58 by HelitronScanner software and further detected 18,668 tandem duplicated Helitron regions (TDHRs) from all the Helitrons identified. There are 39% of TDHRs (7289) translocated since the divergence of the AK58 and Chinese Spring genomes. What interested us even more are the 462 TDHRs exclusive to the AK58 genome. We also found 235 TDHRs in the 21 centromeric regions and these TDHRs contributed to centromere plasticity. Another very interesting DNA transposon, CACTA, accounting for 15% of AK58 genome, was also the focus of this study because they often inserted into gene rich regions. We found that CACTAs have inserted into many agronomically important genes, such as seed dormancy gene TaMFT and vernalization gene TaVrn1, indicating the important role of CACTAs in modern wheat adaptation.     

Rapid proliferation and nucleolar organizer targeting centromeric retrotransposons in cotton

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D‐genome) of the allopolyploid (AD‐genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin‐immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon‐related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A‐genome and related diploid species (B‐, F‐ and G‐genomes), indicating that they colonized the centromeres of D‐genome lineage after the divergence of the A‐ and D‐ ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A‐ and D‐subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D‐ and AD‐genome species, yet localized to just the NORs in A‐, B‐, F‐, and G‐genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.     

Centromere and telomere sequence alterations reflect the rapid genome evolution within the carnivorous plant genus Genlisea

Linear chromosomes of eukaryotic organisms invariably possess centromeres and telomeres to ensure proper chromosome segregation during nuclear divisions and to protect the chromosome ends from deterioration and fusion, respectively. While centromeric sequences may differ between species, with arrays of tandemly repeated sequences and retrotransposons being the most abundant sequence types in plant centromeres, telomeric sequences are usually highly conserved among plants and other organisms. The genome size of the carnivorous genus Genlisea (Lentibulariaceae) is highly variable. Here we study evolutionary sequence plasticity of these chromosomal domains at an intrageneric level. We show that Genlisea nigrocaulis (1C = 86 Mbp; 2n = 40) and G. hispidula (1C = 1550 Mbp; 2n = 40) differ as to their DNA composition at centromeres and telomeres. G. nigrocaulis and its close relative G. pygmaea revealed mainly 161 bp tandem repeats, while G. hispidula and its close relative G. subglabra displayed a combination of four retroelements at centromeric positions. G. nigrocaulis and G. pygmaea chromosome ends are characterized by the Arabidopsis‐type telomeric repeats (TTTAGGG); G. hispidula and G. subglabra instead revealed two intermingled sequence variants (TTCAGG and TTTCAGG). These differences in centromeric and, surprisingly, also in telomeric DNA sequences, uncovered between groups with on average a > 9‐fold genome size difference, emphasize the fast genome evolution within this genus. Such intrageneric evolutionary alteration of telomeric repeats with cytosine in the guanine‐rich strand, not yet known for plants, might impact the epigenetic telomere chromatin modification.     

A novel class of <Emphasis Type="Italic">Helitron</Emphasis>- related transposable elements in maize contain portions of multiple pseudogenes

We recently described a maize mutant caused by an insertion of a Helitron type transposable element (Lal, S.K., Giroux, M.J., Brendel, V., Vallejos, E. and Hannah, L.C., 2003, Plant Cell, 15: 381–391). Here we describe another Helitron insertion in the barren stalk1 gene of maize. The termini of a 6525 bp insertion in the proximal promoter region of the mutant reference allele of maize barren stalk1 gene (ba1-ref) shares striking similarity to the Helitron insertion we reported in the Shrunken-2 gene. This insertion is embedded with pseudogenes that differ from the pseudogenes discovered in the mutant Shrunken-2 insertion. Using the common terminal ends of the mutant insertions as a query, we discovered other Helitron insertions in maize BAC clones. Based on the comparison of the insertion site and PCR amplified genomic sequences, these elements inserted between AT dinucleotides. These putative non-autonomous Helitroninsertions completely lacked sequences similar to RPA (replication protein A) and DNA Helicases reported in other species. A blastn analysis indicated that both the 5 and 3 termini of Helitrons are repeated in the maize genome. These data provide strong evidence that Helitron type transposable elements are active and may have played an essential role in the evolution and expansion of the maize genome.     

Horizontal transfer of a non-autonomous Helitron among insect and viral genomes

Background

The movement of mobile elements among species by horizontal transposon transfer (HTT) influences the evolution of genomes through the modification of structure and function. Helitrons are a relatively new lineage of DNA-based (class II) transposable elements (TEs) that propagate by rolling-circle replication, and are capable of acquiring host DNA. The rapid spread of Helitrons among animal lineages by HTT is facilitated by shuttling in viral particles or by unknown mechanisms mediated by close organism associations (e.g. between hosts and parasites).

Results

A non-autonomous Helitron independently annotated as BmHel-2 from Bombyx mori and the MITE01 element from Ostrinia nubilalis was predicted in the genomes of 24 species in the insect Order Lepidoptera. Integrated Helitrons retained ≥ 65% sequence identity over a 250 bp consensus, and were predicted to retain secondary structures inclusive of a 3′-hairpin and a 5′-subterminal inverted repeat. Highly similar Hel-2 copies were predicted in the genomes of insects and associated viruses, which along with a previous documented case of real-time virus-insect cell line transposition suggests that this Helitron has likely propagated by HTT.

Conclusions

These findings provide evidence that insect virus may mediate the HTT of Helitron-like TEs. This movement may facilitate the shuttling of DNA elements among insect genomes. Further sampling is required to determine the putative role of HTT in insect genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1318-6) contains supplementary material, which is available to authorized users.     

Analysis of transposable elements and organellar DNA in male and female genomes of a species with a huge Y chromosome reveals distinct Y centromeres

Few angiosperms have distinct Y chromosomes. Among those that do are Silene latifolia (Caryophyllaceae), Rumex acetosa (Polygonaceae) and Coccinia grandis (Cucurbitaceae), the latter having a male/female difference of 10% of the total genome (female individuals have a 0.85 pg genome, male individuals 0.94 pg), due to a Y chromosome that arose about 3 million years ago. We compared the sequence composition of male and female C. grandis plants and determined the chromosomal distribution of repetitive and organellar DNA with probes developed from 21 types of repetitive DNA, including 16 mobile elements. The size of the Y chromosome is largely due to the accumulation of certain repeats, such as members of the Ty1/copia and Ty3/gypsy superfamilies, an unclassified element and a satellite, but also plastome‐ and chondriome‐derived sequences. An abundant tandem repeat with a unit size of 144 bp stains the centromeres of the X chromosome and the autosomes, but is absent from the Y centromere. Immunostaining with pericentromere‐specific markers for anti‐histone H3Ser10ph and H2AThr120ph revealed a Y‐specific extension of these histone marks. That the Y centromere has a different make‐up from all the remaining centromeres raises questions about its spindle attachment, and suggests that centromeric or pericentromeric chromatin might be involved in the suppression of recombination.     

Characterization of repeat arrays in ultra‐long nanopore reads reveals frequent origin of satellite DNA from retrotransposon‐derived tandem repeats

Amplification of monomer sequences into long contiguous arrays is the main feature distinguishing satellite DNA from other tandem repeats, yet it is also the main obstacle in its investigation because these arrays are in principle difficult to assemble. Here we explore an alternative, assembly‐free approach that utilizes ultra‐long Oxford Nanopore reads to infer the length distribution of satellite repeat arrays, their association with other repeats and the prevailing sequence periodicities. Using the satellite DNA‐rich legume plant Lathyrus sativus as a model, we demonstrated this approach by analyzing 11 major satellite repeats using a set of nanopore reads ranging from 30 to over 200 kb in length and representing 0.73× genome coverage. We found surprising differences between the analyzed repeats because only two of them were predominantly organized in long arrays typical for satellite DNA. The remaining nine satellites were found to be derived from short tandem arrays located within LTR‐retrotransposons that occasionally expanded in length. While the corresponding LTR‐retrotransposons were dispersed across the genome, this array expansion occurred mainly in the primary constrictions of the L. sativus chromosomes, which suggests that these genome regions are favourable for satellite DNA accumulation.     

Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution

Background

Centromeres are essential for chromosome segregation, yet their DNA sequences evolve rapidly. In most animals and plants that have been studied, centromeres contain megabase-scale arrays of tandem repeats. Despite their importance, very little is known about the degree to which centromere tandem repeats share common properties between different species across different phyla. We used bioinformatic methods to identify high-copy tandem repeats from 282 species using publicly available genomic sequence and our own data.

Results

Our methods are compatible with all current sequencing technologies. Long Pacific Biosciences sequence reads allowed us to find tandem repeat monomers up to 1,419 bp. We assumed that the most abundant tandem repeat is the centromere DNA, which was true for most species whose centromeres have been previously characterized, suggesting this is a general property of genomes. High-copy centromere tandem repeats were found in almost all animal and plant genomes, but repeat monomers were highly variable in sequence composition and length. Furthermore, phylogenetic analysis of sequence homology showed little evidence of sequence conservation beyond approximately 50 million years of divergence. We find that despite an overall lack of sequence conservation, centromere tandem repeats from diverse species showed similar modes of evolution.

Conclusions

While centromere position in most eukaryotes is epigenetically determined, our results indicate that tandem repeats are highly prevalent at centromeres of both animal and plant genomes. This suggests a functional role for such repeats, perhaps in promoting concerted evolution of centromere DNA across chromosomes.     

Genome size‐dependent pcna gene copy number in dinoflagellates and molecular evidence of retroposition as a major evolutionary mechanism

Proliferating cell nuclear antigen (PCNA) plays critical roles in eukaryotic DNA replication and replication‐associated processes. It is typically encoded by one or two gene copies (pcna) in eukaryotic genomes. Recently reported higher copy numbers of pcna in some dinoflagellates raised a question of how this gene has uniquely evolved in this phylum. Through real‐time PCR quantification, we found a wide range of pcna copy number (2–287 copies) in 11 dinoflagellate species (n = 38), and a strong positive correlation between pcna copy number and genome size (log₁₀–log₁₀ transformed). Intraspecific pcna diverged up to 21% and are dominated by nonsynonymous substitutions, indicating strong purifying selection pressure on and hence functional necessity of this gene. By surveying pcna copy numbers in eukaryotes, we observed a genome size threshold at 4 pg DNA, above which more than two pcna copies are found. To examine whether retrotransposition is a mechanism of pcna duplication, we measured the copy number of retroposed pcna, taking advantage of the 22‐nt dinoflagellate‐specific spliced leader (DinoSL) capping the 5′ end of dinoflagellate nuclear‐encoded mRNAs, which would exist in the upstream region of a retroposed gene copy. We found that retroposed pcna copy number increased with total pcna copy number and genome size. These results indicate co‐evolution of dinoflagellate pcna copy number with genome size, and retroposition as a major mechanism of pcna duplication in dinoflagellates. Furthermore, we posit that the demand of faithful replication and maintenance of the large dinoflagellate genomes might have favored the preservation of the retroposed pcna as functional genes.     

Identification of chloroplast genome loci suitable for high‐resolution phylogeographic studies of Colocasia esculenta (L.) Schott (Araceae) and closely related taxa

Recently, we reported the chloroplast genome‐wide association of oligonucleotide repeats, indels and nucleotide substitutions in aroid chloroplast genomes. We hypothesized that the distribution of oligonucleotide repeat sequences in a single representative genome can be used to identify mutational hotspots and loci suitable for population genetic, phylogenetic and phylogeographic studies. Using information on the location of oligonucleotide repeats in the chloroplast genome of taro (Colocasia esculenta), we designed 30 primer pairs to amplify and sequence polymorphic loci. The primers have been tested in a range of intra‐specific to intergeneric comparisons, including ten taro samples (Colocasia esculenta) from diverse geographical locations, four other Colocasia species (C. affinis, C. fallax, C. formosana, C. gigantea) and three other aroid genera (represented by Remusatia vivipara, Alocasia brisbanensis and Amorphophallus konjac). Multiple sequence alignments for the intra‐specific comparison revealed nucleotide substitutions (point mutations) at all 30 loci and microsatellite polymorphisms at 14 loci. The primer pairs reported here reveal levels of genetic variation suitable for high‐resolution phylogeographic and evolutionary studies of taro and other closely related aroids. Our results confirm that information on repeat distribution can be used to identify loci suitable for such studies, and we expect that this approach can be used in other plant groups.     

Wheat centromeric retrotransposons: the new ones take a major role in centromeric structure

The physical map of the hexaploid wheat chromosome 3B was screened using centromeric DNA probes. A 1.1‐Mb region showing the highest number of positive bacterial artificial chromosome (BAC) clones was fully sequenced and annotated, revealing that 96% of the DNA consisted of transposable elements, mainly long terminal repeat (LTR) retrotransposons (88%). Estimation of the insertion time of the transposable elements revealed that CRW (also called Cereba) and Quinta are the youngest elements at the centromeres of common wheat (Triticum spp.) and its diploid ancestors, with Quinta being younger than CRW in both diploid and hexaploid wheats. Chromatin immunoprecipitation experiments revealed that both CRW and Quinta families are targeted by the centromere‐specific histone H3 variant CENH3. Immuno colocalization of retroelements and CENH3 antibody indicated that a higher proportion of Quinta than CRWs was associated with CENH3, although CRWs were more abundant. Long arrays of satellite repeats were also identified in the wheat centromere regions, but they lost the ability to bind with CENH3. In addition to transposons, two functional genes and one pseudogene were identified. The gene density in the centromere appeared to be between three and four times lower than the average gene density of chromosome 3B. Comparisons with related grasses also indicated a loss of microcollinearity in this region. Finally, comparison of centromeric sequences of Aegilops tauschii (DD), Triticum boeoticum (AA) and hexaploid wheat revealed that the centromeres in both the polyploids and diploids are still undergoing dynamic changes, and that the new CRWs and Quintas may have undertaken the core role in kinetochore formation.     

Pif1 Helicases and the Evidence for a Prokaryotic Origin of Helitrons

Helitrons are the only group of rolling-circle transposons that encode a transposase with a helicase domain (Hel), which belongs to the Pif1 family. Because Pif1 helicases are important components of eukaryotic genomes, it has been suggested that Hel domains probably originated after a host eukaryotic Pif1 gene was captured by a Helitron ancestor. However, the few analyses exploring the evolution of Helitron transposases (RepHel) have focused on its Rep domain, which is also present in other mobile genetic elements. Here, we used phylogenetic and nonmetric multidimensional scaling analyses to investigate the relationship between Hel domains and Pif1-like helicases from a variety of organisms. Our results reveal that Hel domains are only distantly related to genomic helicases from eukaryotes and prokaryotes, and thus are unlikely to have originated from a captured Pif1 gene. Based on this evidence, and on recent studies indicating that Rep domains are more closely related to rolling-circle plasmids and phages, we suggest that Helitrons are descendants of a RepHel-encoding prokaryotic plasmid element that invaded eukaryotic genomes before the radiation of its major groups. We discuss how a Pif1-like helicase domain might have favored the transposition of Helitrons in eukaryotes beyond simply unwinding DNA intermediates. Finally, we demonstrate that some examples in the literature describing genomic helicases from eukaryotes actually consist of Hel domains from Helitrons, a finding that underscores how transposons can hamper the analysis of eukaryotic genes. This investigation also revealed that two groups of land plants appear to have lost genomic Pif1 helicases independently.     

Identification and characterization of functional centromeres of the common bean

In higher eukaryotes, centromeres are typically composed of megabase‐sized arrays of satellite repeats that evolve rapidly and homogenize within a species' genome. Despite the importance of centromeres, our knowledge is limited to a few model species. We conducted a comprehensive analysis of common bean (Phaseolus vulgaris) centromeric satellite DNA using genomic data, fluorescence in situ hybridization (FISH), immunofluorescence and chromatin immunoprecipitation (ChIP). Two unrelated centromere‐specific satellite repeats, CentPv1 and CentPv2, and the common bean centromere‐specific histone H3 (PvCENH3) were identified. FISH showed that CentPv1 and CentPv2 are predominantly located at subsets of eight and three centromeres, respectively. Immunofluorescence‐ and ChIP‐based assays demonstrated the functional significance of CentPv1 and CentPv2 at centromeres. Genomic analysis revealed several interesting features of CentPv1 and CentPv2: (i) CentPv1 is organized into an higher‐order repeat structure, named Nazca, of 528 bp, whereas CentPv2 is composed of tandemly organized monomers; (ii) CentPv1 and CentPv2 have undergone chromosome‐specific homogenization; and (iii) CentPv1 and CentPv2 are not likely to be commingled in the genome. These findings suggest that two distinct sets of centromere sequences have evolved independently within the common bean genome, and provide insight into centromere satellite evolution.     

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).     

A lineage‐specific centromere retrotransposon in Oryza brachyantha

Most eukaryotic centromeres contain large quantities of repetitive DNA, such as satellite repeats and retrotransposons. Unlike most transposons in plant genomes, the centromeric retrotransposon (CR) family is conserved over long evolutionary periods among a majority of the grass species. CR elements are highly concentrated in centromeres, and are likely to play a role in centromere function. In order to study centromere evolution in the Oryza (rice) genus, we sequenced the orthologous region to centromere 8 of Oryza sativa from a related species, Oryza brachyantha. We found that O. brachyantha does not have the canonical CRR (CR of rice) found in the centromeres of all other Oryza species. Instead, a new Ty3‐gypsy (Metaviridae) retroelement (FRetro3) was found to colonize the centromeres of this species. This retroelement is found in high copy numbers in the O. brachyantha genome, but not in other Oryza genomes, and based on the dating of long terminal repeats (LTRs) of FRetro3 it was amplified in the genome in the last few million years. Interestingly, there is a high level of removal of FRetro3 based on solo‐LTRs to full‐length elements, and this rapid turnover may have played a role in the replacement of the canonical CRR with the new element by active deletion. Comparison with previously described ChIP cloning data revealed that FRetro3 is found in CENH3‐associated chromatin sequences. Thus, within a single lineage of the Oryza genus, the canonical component of grass centromeres has been replaced with a new retrotransposon that has all the hallmarks of a centromeric retroelement.     

Excision of Helitron Transposons in Maize

Helitrons are novel transposons discovered by bioinformatic analysis of eukaryotic genome sequences. They are believed to move by rolling circle (RC) replication because their predicted transposases are homologous to those of bacterial RC transposons. We report here evidence of somatic Helitron excision in maize, an unexpected finding suggesting that Helitrons can exhibit an excisive mode of transposition.     

Evolution of Tandemly Repeated Sequences: What Happens at the End of an Array?

Tandemly repeated sequences are a major component of the eukaryotic genome. Although the general characteristics of tandem repeats have been well documented, the processes involved in their origin and maintenance remain unknown. In this study, a region on the paternal sex ratio (PSR) chromosome was analyzed to investigate the mechanisms of tandem repeat evolution. The region contains a junction between a tandem array of PSR2 repeats and a copy of the retrotransposon NATE, with other dispersed repeats (putative mobile elements) on the other side of the element. Little similarity was detected between the sequence of PSR2 and the region of NATE flanking the array, indicating that the PSR2 repeat did not originate from the underlying NATE sequence. However, a short region of sequence similarity (11/15 bp) and an inverted region of sequence identity (8 bp) are present on either side of the junction. These short sequences may have facilitated nonhomologous recombination between NATE and PSR2, resulting in the formation of the junction. Adjacent to the junction, the three most terminal repeats in the PSR2 array exhibited a higher sequence divergence relative to internal repeats, which is consistent with a theoretical prediction of the unequal exchange model for tandem repeat evolution. Other NATE insertion sites were characterized which show proximity to both tandem repeats and complex DNAs containing additional dispersed repeats. An ``accretion model' is proposed to account for this association by the accumulation of mobile elements at the ends of tandem arrays and into ``islands' within arrays. Mobile elements inserting into arrays will tend to migrate into islands and to array ends, due to the turnover in the number of intervening repeats. Received: 18 August 1997 / Accepted: 18 September 1998