Identification of NbME MITE families: potential molecular markers in the microsporidia Nosema bombycis

Genome analysis and phylogenetic characterization of two deformed wing virus strains from Apis cerana in Vietnam

Background

Deformed wing virus (DWV) is a virulent virus that causes honeybee disease. DWV can exist as a latent infection in honeybees, outbreak into epidemics, and cause serious damage to beekeeping cross the world, including Vietnam.

Methods

The two DWV strains circulating in Vietnamese honeybee, Apis cerana, were first isolated from adult honeybees in North Vietnam (DWV-NVN) and South Vietnam (DWV-SVN). Their complete nucleotide sequences were determined, aligned, and compared with other DWV strains.

Results

The two Vietnamese DWV strains comprised 10,113 bp and contained a large single open reading frame (ORF) of 2,893 amino acids, initiating at nucleotide 1,130 and terminating at nucleotide 9,812. Multiple nucleotide sequence alignment between these two DWV-VN strains and DWV strains in A. mellifera was performed. The DWV-VN strains showed a low genetic identity (from 91.4% to 92.0%) with almost of these strains, but lower identities (89.2% and 89.4%) with UK2 and (89.6%) with the China2 strain. Low identities (91.7% and 91.9%) were also observed between the China3 strain (in A. cerana) and the DWV-VN strains, respectively. The deduced amino acid sequence alignment showed high genetic similarities (97.0%–97.9%) when the USA1, Chile, Italy1, France, UK1, UK2, Japan, Korea2, China1, China2 and China3 strains were compared to the DWV-VN strains. This ratio was 96.7% and 96.8% when the Korea1 strain was compared to the DWV-SVN and DWV-NVN strains, respectively. Numerous amino acid substitutions were identified in the L, VP3, and RdRp sequences. Notably, we observed six substitutions positioned at amino acids 27 (E > I), 98 (S > T), 120 (A > V), 153 (M > T), 170 (D > F), and 174 (Y > F) in the L protein, two amino acid changes at positions 980 (S > A) and 1032 (E > T) in VP3, and one amino acid change at position 2627 (R > C) unique to the DWV-VN strains. Phylogenetic analysis based on complete genome sequences, RdRp sequences and Simplot analysis indicated that there was a significant difference between DWV-VN strains in A. cerana and DWV strains in A. mellifera. The results suggested that the genetic variations of the DWV-VN strains in A. cerana help them to adapt geographical conditions and may lead to change the viral pathogenicity of DWV-VN strains.

The genomic survey of Tc1-like elements in the silkworm microsporidia Nosema bombycis

BACKGROUND

Microsporidia Nosema bombycis is the destructive pathogen in the production of sericulture. The Tc1/mariner elements belong to important component of DNA transposon.

METHODS

The genomic data of N. bombycis and related Nosema species were screened to identify the Tc1-like elements and analyzed the phylogenetic relationship, based on bioinformational analysis. High-throughput data of transcriptomes and small RNAs were used to evaluate the expressed level and potential rasiRNAs for the Tc1-like elements of N. bombycis.

RESULTS

Twelve complete Tc1-like elements belonging to DD34,E clade is confirmed in the whole genome of N. bombycis, and divided into two branches. Six of them are sole in N. bombycis and thereby would be the molecular marker to differentiate this species from others Nosema spp. Most of the elements have the transcriptional active and are the source of sRNAs.

CONCLUSION

Abundant Tc1-like elements in N. bombycis reflect the expansion of transposons for this genomic characters, comparing with others Nosema spp. The finding of distribution, phylogeny and potential functional activity for Tc1Nbs in N. bombycis will help understanding the role of the DNA transposon in genomic evolution of microsporidia.

Evolutionary dynamics of hAT DNA transposon families in Saccharomycetaceae

Transposable elements (TEs) are widespread in eukaryotes but uncommon in yeasts of the Saccharomycotina subphylum, in terms of both host species and genome fraction. The class II elements are especially scarce, but the hAT element Rover is a noteworthy exception that deserves further investigation. Here, we conducted a genome-wide analysis of hAT elements in 40 ascomycota. A novel family, Roamer, was found in three species, whereas Rover was detected in 15 preduplicated species from Kluyveromyces, Eremothecium, and Lachancea genera, with up to 41 copies per genome. Rover acquisition seems to have occurred by horizontal transfer in a common ancestor of these genera. The detection of remote Rover copies in Naumovozyma dairenensis and in the sole Saccharomyces cerevisiae strain AWRI1631, without synteny, suggests that two additional independent horizontal transfers took place toward these genomes. Such patchy distribution of elements prevents any anticipation of TE presence in incoming sequenced genomes, even closely related ones. The presence of both putative autonomous and defective Rover copies, as well as their diversification into five families, indicate particular dynamics of Rover elements in the Lachancea genus. Especially, we discovered the first miniature inverted-repeat transposable elements (MITEs) to be described in yeasts, together with their parental autonomous copies. Evidence of MITE insertion polymorphism among Lachancea waltii strains suggests their recent activity. Moreover, 40% of Rover copies appeared to be involved in chromosome rearrangements, showing the large structural impact of TEs on yeast genome and opening the door to further investigations to understand their functional and evolutionary consequences.

BuT2 is a member of the third major group of hAT transposons and is involved in horizontal transfer events in the genus Drosophila

The hAT superfamily comprises a large and diverse array of DNA transposons found in all supergroups of eukaryotes. Here we characterized the Drosophila buzzatii BuT2 element and found that it harbors a five-exon gene encoding a 643-aa putatively functional transposase. A phylogeny built with 85 hAT transposases yielded, in addition to the two major groups already described, Ac and Buster, a third one comprising 20 sequences that includes BuT2, Tip100, hAT-4_BM, and RP-hAT1. This third group is here named Tip. In addition, we studied the phylogenetic distribution and evolution of BuT2 by in silico searches and molecular approaches. Our data revealed BuT2 was, most often, vertically transmitted during the evolution of genus Drosophila being lost independently in several species. Nevertheless, we propose the occurrence of three horizontal transfer events to explain its distribution and conservation among species. Another aspect of BuT2 evolution and life cycle is the presence of short related sequences, which contain similar 5' and 3' regions, including the terminal inverted repeats. These sequences that can be considered as miniature inverted repeat transposable elements probably originated by internal deletion of complete copies and show evidences of recent mobilization.

Comparative genomics of parasitic silkworm microsporidia reveal an association between genome expansion and host adaptation

Background

Microsporidian Nosema bombycis has received much attention because the pébrine disease of domesticated silkworms results in great economic losses in the silkworm industry. So far, no effective treatment could be found for pébrine. Compared to other known Nosema parasites, N. bombycis can unusually parasitize a broad range of hosts. To gain some insights into the underlying genetic mechanism of pathological ability and host range expansion in this parasite, a comparative genomic approach is conducted. The genome of two Nosema parasites, N. bombycis and N. antheraeae (an obligatory parasite to undomesticated silkworms Antheraea pernyi), were sequenced and compared with their distantly related species, N. ceranae (an obligatory parasite to honey bees).

Results

Our comparative genomics analysis show that the N. bombycis genome has greatly expanded due to the following three molecular mechanisms: 1) the proliferation of host-derived transposable elements, 2) the acquisition of many horizontally transferred genes from bacteria, and 3) the production of abundnant gene duplications. To our knowledge, duplicated genes derived not only from small-scale events (e.g., tandem duplications) but also from large-scale events (e.g., segmental duplications) have never been seen so abundant in any reported microsporidia genomes. Our relative dating analysis further indicated that these duplication events have arisen recently over very short evolutionary time. Furthermore, several duplicated genes involving in the cytotoxic metabolic pathway were found to undergo positive selection, suggestive of the role of duplicated genes on the adaptive evolution of pathogenic ability.

Conclusions

Genome expansion is rarely considered as the evolutionary outcome acting on those highly reduced and compact parasitic microsporidian genomes. This study, for the first time, demonstrates that the parasitic genomes can expand, instead of shrink, through several common molecular mechanisms such as gene duplication, horizontal gene transfer, and transposable element expansion. We also showed that the duplicated genes can serve as raw materials for evolutionary innovations possibly contributing to the increase of pathologenic ability. Based on our research, we propose that duplicated genes of N. bombycis should be treated as primary targets for treatment designs against pébrine.

Miniature inverted-repeat transposable elements: discovery, distribution, and activity

Eukaryotic organisms have dynamic genomes, with transposable elements (TEs) as a major contributing factor. Although the large autonomous TEs can significantly shape genomic structures during evolution, genomes often harbor more miniature nonautonomous TEs that can infest genomic niches where large TEs are rare. In spite of their cut-and-paste transposition mechanisms that do not inherently favor copy number increase, miniature inverted-repeat transposable elements (MITEs) are abundant in eukaryotic genomes and exist in high copy numbers. Based on the large number of MITE families revealed in previous studies, accurate annotation of MITEs, particularly in newly sequenced genomes, will identify more genomes highly rich in these elements. Novel families identified from these analyses, together with the currently known families, will further deepen our understanding of the origins, transposase sources, and dramatic amplification of these elements.

Characterization of active ribosomal RNA harboring MITEs insertion in microsporidian Nosema bombycis genome

Microsporidia are a group of obligate intracellular parasites of medical and agricultural importance, which can infect almost all animals, including human beings. Using the genome data of Nosema bombycis, four families of miniature inverted-repeat transposable elements (MITEs) in ribosomal DNA (rDNA) were characterized in the microsporidian N. bombycis and were named LSUME1, ITSME1, SSUME1, and SSUME2, respectively. The genome-wide investigation of these MITEs shows that these MITEs families distribute randomly in N. bombycis genome. All insertion sequences have conserved characteristics of MITEs, the direct repeat sequence and terminal inverted-repeat sequence at both ends of each MITEs sequence. Additionally, using the CLC RNA Workbench Software, secondary structures of rRNA containing MITEs sequence have been predicted and were located in variable region or expansion segment. Furthermore, using two different probes, one is prepared by MITE sequence only (short probe) and the other is prepared by MITE sequence flanking partial rDNA sequence (long probe); northern blotting and dot blotting have been performed to detect the transcriptional and functional activity of the rDNA containing MITEs insertion. Fortunately, we found that the rDNA, which harbors the MITE, not only can be transcripted but also can form a complete ribosome. This is an interesting thing that one gene can keep active even when it has been inserted with another sequence. But the biological and structural significance of this observation is not readily apparent.

Mobilization of Stowaway-like MITEs in newly formed allohexaploid wheat species

Transposable elements (TEs) dominate the genetic capacity of most eukaryotes, especially plants, where they can account for up to 90 % of the genome, such as in wheat. The relationship between TEs and their hosts and the role of TEs in organismal biology are poorly understood. In this study, we have applied next generation sequencing, together with a transposon display technique in order to test whether a Stowaway-like MITE, termed Minos, transposes following allopolyploidization events in wheat. We have generated a 454-pyrosequencing database of Minos-specific amplicons (transposon display products) from a newly formed wheat allohexaploid and its parental lines and retrieved hundreds of novel MITE insertions in the allohexaploid. Clear mobilization of Minos was also seen by site-specific PCR analysis and sequence validation. In addition, using real-time qPCR analysis we observed an insignificant change in the relative quantity of Minos from the expected value of merging the two parental genomes, indicating that, despite its activation, no significant burst in Minos copy number can be seen in the newly formed allohexaploid. Interestingly, we found that CCGG sites surrounding Minos underwent massive hypermethylation following the allohexaploidization process. Our data suggest that MITEs have maintained their capacity for activity throughout the evolution of wheat and might be epigenetically deregulated in the first generations following allopolyploidization.

Mar, a MITE family of hAT transposons in Drosophila

Background

Miniature inverted-repeat transposable elements (MITEs) are short, nonautonomous DNA elements flanked by subterminal or terminal inverted repeats (TIRs) with no coding capacity. MITEs were originally recognized as important components of plant genomes, where they can attain extremely high copy numbers, and are also found in several animal genomes, including mosquitoes, fish and humans. So far, few MITEs have been described in Drosophila.

Results

Herein we describe the distribution and evolution of Mar, a MITE family of hAT transposons, in Drosophilidae species. In silico searches and PCR screening showed that Mar distribution is restricted to the willistoni subgroup of the Drosophila species, and a phylogenetic analysis of Mar indicates that this element may have originated prior to the diversification of these species. Most of the Mar copies in D. willistoni present conserved target site duplications and TIRs, indicating recent mobilization of these sequences. We also identified relic copies of potentially full-length Mar transposon in D. tropicalis and D. willistoni. The phylogenetic relationship among transposases from the putative full-length Mar and other hAT superfamily elements revealed that Mar is placed into the recently determined Buster group of hAT transposons.

Conclusion

On the basis of the obtained data, we can suggest that the origin of these Mar MITEs occurred before the subgroup willistoni speciation, which started about 5.7 Mya. The Mar relic transposase existence indicates that these MITEs originated by internal deletions and suggests that the full-length transposon was recently functional in D. willistoni, promoting Mar MITEs mobilization.

msechBari, a new MITE-like element in Drosophila sechellia related to the Bari transposon

A few occurrences of miniature inverted-repeat transposable elements (MITEs) have been reported in species of the genus Drosophila. Here, we describe msechBari, a MITE-like element in Drosophila sechellia. The element is short, approximately 90 bp in length, AT-rich and occurs in association with, or close to, genes, characteristics that are typical for MITEs. The identification was performed in silico using the sequenced genome of D. sechellia and confirmed in a laboratory strain. This short element is related to the Bari_DM transposon of Drosophila melanogaster, having terminal inverted repeats (TIRs) of a similar length and a high identity with the full-length Bari_DM element. The estimated recent origin of the element and the homogeneity observed between copies found in the genome suggests that msechBari could be active in D. sechellia.

Extreme reduction and compaction of microsporidian genomes

Microsporidia are fungi-related obligate intracellular parasites with a highly reduced and compact genome, as for Encephalitozoon species which harbor a genome smaller than 3 Mbp. Genome compaction is reflected by high gene density and, for larger microsporidian genomes, size variation is due to repeat elements that do not drastically affect gene density. Furthermore, these pathogens present strong host dependency illustrated by extensive gene loss. Such adaptations associated with genome compaction induced gene size reduction but also simplification of cellular processes such as transcription. Thus, microsporidia are excellent models for eukaryotic genome evolution and gene expression in the context of host-pathogen relationships.