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.     

Rolling‐circle amplification of centromeric Helitrons in plant genomes

The unusual eukaryotic Helitron transposons can readily capture host sequences and are, thus, evolutionarily important. They are presumed to amplify by rolling‐circle replication (RCR) because some elements encode predicted proteins homologous to RCR prokaryotic transposases. In support of this replication mechanism, it was recently shown that transposition of a bat Helitron generates covalently closed circular intermediates. Another strong prediction is that RCR should generate tandem Helitron concatemers, yet almost all Helitrons identified to date occur as solo elements in the genome. To investigate alternative modes of Helitron organization in present‐day genomes, we have applied the novel computational tool HelitronScanner to 27 plant genomes and have uncovered numerous tandem arrays of partially decayed, truncated Helitrons in all of them. Strikingly, most of these Helitron tandem arrays are interspersed with other repeats in centromeres. Many of these arrays have multiple Helitron 5′ ends, but a single 3′ end. The number of repeats in any one array can range from a handful to several hundreds. We propose here an RCR model that conforms to the present Helitron landscape of plant genomes. Our study provides strong evidence that plant Helitrons amplify by RCR and that the tandemly arrayed replication products accumulate mostly in centromeres.     

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.     

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.     

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.     

Detection and characterization of Wolbachia infection in silkworm

Wolbachia naturally infects a wide variety of arthropods, where it plays important roles in host reproduction. It was previously reported that Wolbachia did not infect silkworm. By means of PCR and sequencing we found in this study that Wolbachia is indeed present in silkworm. Phylogenetic analysis indicates that Wolbachia infection in silkworm may have occurred via transfer from parasitic wasps. Furthermore, Southern blotting results suggest a lateral transfer of the wsp gene into the genomes of some wild silkworms. By antibiotic treatments, we found that tetracycline and ciprofloxacin can eliminate Wolbachia in the silkworm and Wolbachia is important to ovary development of silkworm. These results provide clues towards a more comprehensive understanding of the interaction between Wolbachia and silkworm and possibly other lepidopteran insects.     

家蚕LTR逆转录转座子的鉴定、分类及系统发育分析

转座子是真核生物基因组的重要组成成分。为了研究家蚕Bombyx mori长末端重复序列 (long terminal repeat, LTR)逆转录转座子的分类及进化, 本研究采用de novo预测和同源性搜索相结合的方法, 在家蚕基因组中共鉴定出了38个LTR逆转录转座子家族, 序列长度占整个基因组的0.64%, 远小于先前预测的11.8%, 其中有6个家族为本研究的新发现。38个家族中, 26个家族有表达序列标签 (expression sequence tag, EST)证据, 表明这些家族具有潜在的活性。对有EST证据的6个家族和没有EST证据的5个家族用RT-PCR进行了组织表达谱实验, 结果表明这11个家族在一些组织中有表达, 这进一步证实了这些家族具有转录活性, 基于此我们推测家蚕中大部分的LTR逆转录转座子家族很可能具有潜在活性。对转座子的插入时间进行估计, 结果表明绝大部分元件都是最近1百万年内插入到家蚕基因组中的。我们还比较了黑腹果蝇Drosophila melanogaster、 冈比亚按蚊Anopheles gambiae和家蚕B. mori中Ty3/Gypsy超家族分支的差异, 结果表明不同枝在不同昆虫中有着不同的扩张。家蚕中LTR逆转录转座子的鉴定和系统分析有助于我们理解逆转录转座子在昆虫进化中的作用。     

Identification of Recombination and Positively Selected Genes in Brucella

Brucella is a facultative intracellular bacterium belongs to the class alpha proteobacteria. It causes zoonotic disease brucellosis to wide range of animals. Brucella species are highly conserved in nucleotide level. Here, we employed a comparative genomics approach to examine the role of homologous recombination and positive selection in the evolution of Brucella. For the analysis, we have selected 19 complete genomes from 8 species of Brucella. Among the 1599 core genome predicted, 24 genes were showing signals of recombination but no significant breakpoint was found. The analysis revealed that recombination events are less frequent and the impact of recombination occurred is negligible on the evolution of Brucella. This leads to the view that Brucella is clonally evolved. On other hand, 56 genes (3.5 % of core genome) were showing signals of positive selection. Results suggest that natural selection plays an important role in the evolution of Brucella. Some of the genes that are responsible for the pathogenesis of Brucella were found positively selected, presumably due to their role in avoidance of the host immune system.     

Identification of MicroRNAs in Helicoverpa armigera and Spodoptera litura Based on Deep Sequencing and Homology Analysis

The current identification of microRNAs (miRNAs) in insects is largely dependent on genome sequences. However, the lack of available genome sequences inhibits the identification of miRNAs in various insect species. In this study, we used a miRNA database of the silkworm Bombyx mori as a reference to identify miRNAs in Helicoverpa armigera and Spodoptera litura using deep sequencing and homology analysis. Because all three species belong to the Lepidoptera, the experiment produced reliable results. Our study identified 97 and 91 conserved miRNAs in H. armigera and S. litura, respectively. Using the genome of B. mori and BAC sequences of H. armigera as references, 1 novel miRNA and 8 novel miRNA candidates were identified in H. armigera, and 4 novel miRNA candidates were identified in S. litura. An evolutionary analysis revealed that most of the identified miRNAs were insect-specific, and more than 20 miRNAs were Lepidoptera-specific. The investigation of the expression patterns of miR-2a, miR-34, miR-2796-3p and miR-11 revealed their potential roles in insect development. miRNA target prediction revealed that conserved miRNA target sites exist in various genes in the 3 species. Conserved miRNA target sites for the Hsp90 gene among the 3 species were validated in the mammalian 293T cell line using a dual-luciferase reporter assay. Our study provides a new approach with which to identify miRNAs in insects lacking genome information and contributes to the functional analysis of insect miRNAs.     

Characterization and potential evolutionary impact of transposable elements in the genome of Cochliobolus heterostrophus

Background

Cochliobolus heterostrophus is a dothideomycete that causes Southern Corn Leaf Blight disease. There are two races, race O and race T that differ by the absence (race O) and presence (race T) of ~ 1.2-Mb of DNA encoding genes responsible for the production of T-toxin, which makes race T much more virulent than race O. The presence of repetitive elements in fungal genomes is considered to be an important source of genetic variability between different species.

Results

A detailed analysis of class I and II TEs identified in the near complete genome sequence of race O was performed. In total in race O, 12 new families of transposons were identified. In silico evidence of recent activity was found for many of the transposons and analyses of expressed sequence tags (ESTs) demonstrated that these elements were actively transcribed. Various potentially active TEs were found near coding regions and may modify the expression and structure of these genes by acting as ectopic recombination sites. Transposons were found on scaffolds carrying polyketide synthase encoding genes, responsible for production of T-toxin in race T. Strong evidence of ectopic recombination was found, demonstrating that TEs can play an important role in the modulation of genome architecture of this species. The Repeat Induced Point mutation (RIP) silencing mechanism was shown to have high specificity in C. heterostrophus, acting only on transposons near coding regions.

Conclusions

New families of transposons were identified. In C. heterostrophus, the RIP silencing mechanism is efficient and selective. The co-localization of effector genes and TEs, therefore, exposes those genes to high rates of point mutations. This may accelerate the rate of evolution of these genes, providing a potential advantage for the host. Additionally, it was shown that ectopic recombination promoted by TEs appears to be the major event in the genome reorganization of this species and that a large number of elements are still potentially active. So, this study provides information about the potential impact of TEs on the evolution of C. heterostrophus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-536) contains supplementary material, which is available to authorized users.     

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.     

Genome-wide identification and evolution of <Emphasis Type="Italic">TC1</Emphasis>/<Emphasis Type="Italic">Mariner</Emphasis> in the silkworm (<Emphasis Type="Italic">Bombyx mori</Emphasis>) genome

TC1/Mariner transposons belong to class II transposable elements (TEs) that use DNA-mediated “cut and paste” mechanism to transpose, and they have been identified in almost all organisms. Although silkworm (Bombyx mori) has a large amount of TC1/Mariner elements, the genome wide information of this superfamily in the silkworm is unknown. In this study, we have identified 2670 TC1/Mariner (Bmmar) elements in the silkworm genome. All the TEs were classified into 22 families by means of fgclust, a tool of repetitive sequence classification, seven of which was first reported in this study. Phylogenetic and structure analyses based on the catalytic domain (DDxD/E) of transposase sequences indicated that all members of TC1/Mariner were grouped into five subgroups: Mariner, Tc1, maT, DD40D and DD41D/E. Of these five subgroups, maT rather than Mariner possessed most members of TC1/Mariner (51.23%) in the silkworm genome. In particular, phylogenetic analysis and structure analysis revealed that Bmmar15 (DD40D) formed a new basal subgroup of TC1/Mariner element in insects, which was referred to as bmori. Furthermore, we concluded that DD40D appeared to intermediate between mariner and Tc1. Finally, we estimated the insertion time for each copy of TC1/Mariner in the silkworm and found that most of members were dramatically amplified during a period from 0 to 1 mya. Moreover, the detailed functional data analysis showed that Bmmar1, Bmmar6 and Bmmar9 had EST evidence and intact transposases. These implied that TC1/Mariner might have potential transpositional activity. In conclusion, this study provides some new insights into the landscape, origin and evolution of TC1/Mariner in the insect genomes.     

Bursts of transposition from non-long terminal repeat retrotransposon families of the RTE clade in Schistosoma mansoni

The genus Schistosoma is composed of blood flukes that infect vertebrates, from which three species are major causative agents of human schistosomiasis, a tropical disease that affects more than 200 million people. Current models of the recent evolution of Schistosoma indicate multiple events of migration and speciation from an Asian ancestral species. Transposable elements are important drivers of genome evolution and have been hypothesised to have an important role in speciation. In this work, we describe a comprehensive inventory of Schistosoma mansoni and Schistosoma japonicum retrotransposons, based on their recently published genomic data. We find a considerable difference in retrotransposon representation between the two species (22% and 13%, respectively). A large part of this difference can be attributed to higher representation of two previously described families of S. mansoni retrotransposons (SR2 and Perere-3/SR3), compared with the representation of their closest relative families in S. japonicum. A more detailed analysis suggests that these two S. mansoni families were the subject of recent bursts of transposition that were not paralleled by their S. japonicum counterparts. We hypothesise that these bursts could be a consequence of the evolutionary pressure resulting from migration of Schistosoma from Asia to Africa and their establishment in this new environment, helping both speciation and adaptation.     

Effects of egt gene transfer on the development of Bombyx mori

This study investigated the effects of gain of ecdysteroid UDP-glucosyltransferase (EGT) gene function mutation on the development of the silkworm, Bombyx mori. A novel piggyBac-derived plasmid containing the egt gene from B. mori nucleopolyhedrovirus (BmNPV) driven by a heat-shock protein (hsp) 23.7 promoter, with a neomycin-resistance gene (neo) controlled by the BmNPV ie-1 promoter and a green fluorescent protein gene (gfp) under the control of the B. mori actin 3 (A3) promoter was constructed. The vector was transferred into silkworm eggs by sperm-mediated gene transfer. Transgenic silkworms were produced after screening for neo and gfp genes and gene transfer was verified by polymerase chain reaction, dot-blot hybridization and western blotting. The hatching rate of G1 generation silkworm eggs was about 60% lower than that of normal silkworm eggs. The duration of the G1 generation larval period was extended, and the G2 generation pupal stage lasted four days longer than that in non-transgenic silkworms. The ecdysone blood level in G2 silkworms in the third instar molting stage was reduced by up to 90%. These results show that EGT suppressed transgenic silkworm molting, and that egt expression in egt-transgenic silkworms resulted in arrest of metamorphosis from pupae to moths.     

Diverse retrotransposon families and an AT-rich satellite DNA revealed in giant genomes of Fritillaria lilies

Background and Aims

The genus Fritillaria (Liliaceae) comprises species with extremely large genomes (1C = 30 000–127 000 Mb) and a bicontinental distribution. Most North American species (subgenus Liliorhiza) differ from Eurasian Fritillaria species by their distinct phylogenetic position and increased amounts of heterochromatin. This study examined the contribution of major repetitive elements to the genome obesity found in Fritillaria and identified repeats contributing to the heterochromatin arrays in Liliorhiza species.

Methods

Two Fritillaria species of similar genome size were selected for detailed analysis, one from each phylogeographical clade: F. affinis (1C = 45·6 pg, North America) and F. imperialis (1C = 43·0 pg, Eurasia). Fosmid libraries were constructed from their genomic DNAs and used for identification, sequence characterization, quantification and chromosome localization of clones containing highly repeated sequences.

Key Results and Conclusions

Repeats corresponding to 6·7 and 4·7 % of the F. affinis and F. imperialis genome, respectively, were identified. Chromoviruses and the Tat lineage of Ty3/gypsy group long terminal repeat retrotransposons were identified as the predominant components of the highly repeated fractions in the F. affinis and F. imperialis genomes, respectively. In addition, a heterogeneous, extremely AT-rich satellite repeat was isolated from F. affinis. The FriSAT1 repeat localized in heterochromatic bands makes up approx. 26 % of the F. affinis genome and substantial genomic fractions in several other Liliorhiza species. However, no evidence of a relationship between heterochromatin content and genome size variation was observed. Also, this study was unable to reveal any predominant repeats which tracked the increasing/decreasing trends of genome size evolution in Fritillaria. Instead, the giant Fritillaria genomes seem to be composed of many diversified families of transposable elements. We hypothesize that the genome obesity may be partly determined by the failure of removal mechanisms to counterbalance effectively the retrotransposon amplification.     

Identification,Diversity and Evolution of MITEs in the Genomes of Microsporidian Nosema Parasites

Miniature inverted-repeat transposable elements (MITEs) are short, non-autonomous DNA transposons, which are widespread in most eukaryotic genomes. However, genome-wide identification, origin and evolution of MITEs remain largely obscure in microsporidia. In this study, we investigated structural features for de novo identification of MITEs in genomes of silkworm microsporidia Nosema bombycis and Nosema antheraeae, as well as a honeybee microsporidia Nosema ceranae. A total of 1490, 149 and 83 MITE-related sequences from 89, 17 and five families, respectively, were found in the genomes of the above-mentioned species. Species-specific MITEs are predominant in each genome of microsporidian Nosema, with the exception of three MITE families that were shared by N. bombycis and N. antheraeae. One or multiple rounds of amplification occurred for MITEs in N. bombycis after divergence between N. bombycis and the other two species, suggesting that the more abundant families in N. bombycis could be attributed to the recent amplification of new MITEs. Significantly, some MITEs that inserted into the homologous protein-coding region of N. bombycis were recruited as introns, indicating that gene expansion occurred during the evolution of microsporidia. NbS31 and NbS24 had polymorphisms in different geographical strains of N. bombycis, indicating that they could still be active. In addition, several small RNAs in the MITEs in N. bombycis are mainly produced from both ends of the MITEs sequence.     

Genome-wide analysis of tandem repeats in Tribolium castaneum genome reveals abundant and highly dynamic tandem repeat families with satellite DNA features in euchromatic chromosomal arms

Although satellite DNAs are well-explored components of heterochromatin and centromeres, little is known about emergence, dispersal and possible impact of comparably structured tandem repeats (TRs) on the genome-wide scale. Our bioinformatics analysis of assembled Tribolium castaneum genome disclosed significant contribution of TRs in euchromatic chromosomal arms and clear predominance of satellite DNA-typical 170 bp monomers in arrays of ≥5 repeats. By applying different experimental approaches, we revealed that the nine most prominent TR families Cast1–Cast9 extracted from the assembly comprise ∼4.3% of the entire genome and reside almost exclusively in euchromatic regions. Among them, seven families that build ∼3.9% of the genome are based on ∼170 and ∼340 bp long monomers. Results of phylogenetic analyses of 2500 monomers originating from these families show high-sequence dynamics, evident by extensive exchanges between arrays on non-homologous chromosomes. In addition, our analysis shows that concerted evolution acts more efficiently on longer than on shorter arrays. Efficient genome-wide distribution of nine TR families implies the role of transposition only in expansion of the most dispersed family, and involvement of other mechanisms is anticipated. Despite similarities in sequence features, FISH experiments indicate high-level compartmentalization of centromeric and euchromatic tandem repeats.     

The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms

Background

Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods.

Results

We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology.

Conclusions

Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-720) contains supplementary material, which is available to authorized users.