MAGGY, a retrotransposon in the genome of the rice blast fungusMagnaporthe grisea

Full-length copies of a previously described repetitive DNA sequence (CH2-8) were isolated from the genome of theMagnaporthe grisea strain 2539. One copy of the complete element was sequenced and found to resemble agypsy-like LTR retrotransposon. We named this element MAGGY (MAGnaporthe GYpsy-like element). MAGGY contains two internal ORFs putatively encoding Gag, Pol and Env-like proteins which are similar to peptides encoded by retroelements identified in other filamentous fungi. MAGGY was found to be widely distributed amongM. grisea isolates from geographically dispersed locations and different hosts. It was present in high copy number in the genomes of all nine rice-pathogenic isolates examined. By contrast,M. grisea strains isolated from other Gramineae were found to possess varying copy numbers of MAGGY and in some cases the element was completely absent. The wide distribution of MAGGY suggests that this element invaded the genome ofM. grisea prior to the evolution of rice-specific form(s). It may since have been horizontally transmitted to other sub-specific groups. One copy of MAGGY, corresponding to the element we sequenced, was located at identical locations in the genomes of geographically dispersed strains, suggesting that this copy of the element is a relatively ancient insertion.     

Degenerate MAGGY elements in a subgroup of Pyricularia grisea: a possible example of successful capture of a genetic invader by a fungal genome

The LTR-retrotransposon MAGGY is found sporadically in isolates of Pyricularia grisea (Magnaporthe grisea). Based on a dendrogram constructed by RFLP analysis of rDNA, isolates that carry MAGGY elements were classified into a single cluster that comprised four rDNA types. However, in a few members of this cluster, exemplified by isolates from common millet (Panicum miliaceum), the MAGGY element has distinct features. Southern analysis suggested that these isolates possessed a single copy of a MAGGY-related sequence whose restriction map differed from that of MAGGY itself. Sequence analysis revealed that the MAGGY-related sequence was a degenerate form of MAGGY, characterized by numerous C:G to T:A transitions, which have often been reported to result from RIP (Repeat-induced point mutation) or RIP-like processes. However, the favored target site for C:G to T:A transitions in this fungus, determined by examining a total of 501 sites, was (A/T)pCp(A/T), which differs from that for the RIP process originally reported in Neurospora (CpA), and from that reported in Aspergillus (CpG). The fact that certain members of the cluster of MAGGY carriers retain a single copy of a degenerate MAGGY element implies that the ancestor of these isolates successfully “captured” the invading MAGGY element.     

Degenerate MAGGY elements in a subgroup of Pyricularia grisea : a possible example of successful capture of a genetic invader by a fungal genome

The LTR-retrotransposon MAGGY is found sporadically in isolates of Pyricularia grisea (Magnaporthe grisea). Based on a dendrogram constructed by RFLP analysis of rDNA, isolates that carry MAGGY elements were classified into a single cluster that comprised four rDNA types. However, in a few members of this cluster, exemplified by isolates from common millet (Panicum miliaceum), the MAGGY element has distinct features. Southern analysis suggested that these isolates possessed a single copy of a MAGGY-related sequence whose restriction map differed from that of MAGGY itself. Sequence analysis revealed that the MAGGY-related sequence was a degenerate form of MAGGY, characterized by numerous C:G to T:A transitions, which have often been reported to result from RIP (Repeat-induced point mutation) or RIP-like processes. However, the favored target site for C:G to T:A transitions in this fungus, determined by examining a total of 501 sites, was (A/T)pCp(A/T), which differs from that for the RIP process originally reported in Neurospora (CpA), and from that reported in Aspergillus (CpG). The fact that certain members of the cluster of MAGGY carriers retain a single copy of a degenerate MAGGY element implies that the ancestor of these isolates successfully “captured” the invading MAGGY element. Received: 1 February 1999 / Accepted: 8 April 1999     

Genome organization of Magnaporthe grisea: integration of genetic maps, clustering of transposable elements and identification of genome duplications and rearrangements

 A high-density genetic map of the rice blast fungus Magnaporthe grisea (Guy11×2539) was constructed by adding 87 cosmid-derived RFLP markers to previously generated maps. The new map consists of 203 markers representing 132 independently segregating loci and spans approximately 900 cM with an average resolution of 4.5 cM. Mapping of 33 cosmid probes from the genetic map generated by Sweigard et al. has allowed the integration of two M. grisea maps. The integrated map showed that the linear order of markers along all seven chromosomes in both maps is in good agreement. Thirty of eighty seven markers were derived from cosmid clones that contained the retrotransposon MAGGY (M. grisea gypsy element). Mapping of single-copy DNA sequences associated with the MAGGY cosmids indicated that MAGGY elements are scattered throughout the fungal genome. In eight cases, the probes associated with MAGGY elements showed abnormal segregation patterns. This suggests that MAGGY may be involved in genomic rearrangements. Two RFLP probes linked to MAGGY elements, and another flanking other repetitive DNA elements, identified sequences that were duplicated in the Guy11 genome. Most of the MAGGY cosmids also contained other classes of repetitive DNA suggesting that repetitive DNA sequences tend to cluster in the M. grisea genome. Received: 17 February 1997 / Accepted: 21 February 1997     

TheMagnaporthe grisea DNA fingerprinting probe MGR586 contains the 3′ end of an inverted repeat transposon

TheMagnaporthe grisea repeat (MGR) sequence MGR586 has been widely used for population studies of the rice blast fungus, and has enabled classification of the fungal population into hundreds of genetic lineages. While studying the distribution of MGR586 sequences in strains ofM. grisea, we discovered that the plasmid probe pCB586 contains a significant amount of single-copy DNA. To define precisely the boundary of the repetitive DNA in pCB586, this plasmid and four cosmid clones containing MGR586 were sequenced. Only 740 bp of one end of the 2.6-bp insert in the pCB586 plasmid was common to all clones. DNA sequence analysis of cosmid DNA revealed that all the cosmids contained common sequences beyond the cloning site in pCB586, indicating that the repetitive DNA in the fingerprinting clone is part of a larger element. The entire repetitive element was sequenced and found to resemble an inverted repeat transposon. This putative transposon is 1.86 kb in length and has perfect terminal repeats of 42 bp, which themselves contain direct repeats of 16 bp. The internal region of the transposon possesses one open reading frame which shows similarity at the peptide level to the Pot2 transposon fromM. grisea and Fot1 fromFusarium oxysporum. Hybridization studies using the entire element as a probe revealed that some strains ofM. grisea, whose DNA hybridized to the pCB586 probe, entirely lacked MGR586 transposon sequences.     

Lineage-Specific Distribution of Insertion Sequence Excision Enhancer in Enterotoxigenic Escherichia coli Isolated from Swine

Insertion sequences (ISs) are the simplest transposable elements and are widely distributed in bacteria; however, they also play important roles in genome evolution. We recently identified a protein called IS excision enhancer (IEE) in enterohemorrhagic Escherichia coli (EHEC) O157. IEE promotes the excision of IS elements belonging to the IS3 family, such as IS629, as well as several other families. IEE-mediated IS excision generates various genomic deletions that lead to the diversification of the bacterial genome. IEE has been found in a broad range of bacterial species; however, among sequenced E. coli strains, IEE is primarily found in EHEC isolates. In this study, we investigated non-EHEC pathogenic E. coli strains isolated from domestic animals and found that IEE is distributed in specific lineages of enterotoxigenic E. coli (ETEC) strains of serotypes O139 or O149 isolated from swine. The iee gene is located within integrative elements that are similar to SpLE1 of EHEC O157. All iee-positive ETEC lineages also contained multiple copies of IS629, a preferred substrate of IEE, and their genomic locations varied significantly between strains, as observed in O157. These data suggest that IEE may have been transferred among EHEC and ETEC in swine via SpLE1 or SpLE1-like integrative elements. In addition, IS629 is actively moving in the ETEC O139 and O149 genomes and, as in EHEC O157, is promoting the diversification of these genomes in combination with IEE.     

Comparative genomic analysis of two Burkholderia glumae strains from different geographic origins reveals a high degree of plasticity in genome structure associated with genomic islands

Burkholderia glumae is the major causal agent of bacterial panicle blight of rice, a growing disease problem in global rice production. To better understand its genome-scale characteristics, the genome of the highly virulent B. glumae strain 336gr-1 isolated from Louisiana, USA was sequenced using the Illumina Genome Analyser II system. De novo assembled 336gr-1 contigs were aligned and compared with the previously sequenced genome of B. glumae strain BGR1, which was isolated from an infected rice plant in South Korea. Comparative analysis of the whole genomes of B. glumae 336gr-1 and B. glumae BGR1 revealed numerous unique genomic regions present only in one of the two strains. These unique regions contained accessory genes including mobile elements and phage-related genes, and some of the unique regions in B. glumae BGR1 corresponded to predicted genomic islands. In contrast, little variation was observed in known and potential virulence genes between the two genomes. The considerable amount of plasticity largely based on accessory genes and genome islands observed from the comparison of the genomes of these two strains of B. glumae may explain the versatility of this bacterial species in various environmental conditions and geographic locations.     

Genome organization of Magnaporthe grisea: genetic map,electrophoretic karyotype,and occurrence of repeated DNAs

A genetic map of Magnaporthe grisea (anamorph=Pyricularia oryzae and P. grisea), the causal agent of rice blast disease, was generated from segregation data utilizing 97 RFLP markers, two isoenzyme loci and the mating type locus among progeny of a cross between parental strains Guy 11 and 2539. Of the seven chromosomes of M. Grisea, three were resolved by contour-clamped homogeneous electric field (CHEF) electrophoresis, while the remaining four migrated as two doublet bands. By utilizing differences between CHEF mobilities of unresolved chromosomes from the parental strains, Southern analysis with selected markers allowed the chromosomal assignment of all linkage groups. A small translocation involving 1 marker was found in the parental strains used to produce the segregating population from which the map was constructed. Nine classes of repetitive DNA elements were found in the genome of a fungal isolate pathogenic to rice. These occurred only a few times or not at all in the genomes of isolates showing reduced virulence on rice. One repetitive DNA was shown to have structural similarity to the Alu sequences found in primates, a sequence similarity to the copia-like elements of Drosophila, and peptide similarity to transposable elements found in Drosophila, other fungi, and higher plants.The mention of a trademark, proprietary product, or vendor anywhere in this paper does not constitute a guarantee or warranty of the product by the USDA-ARS and does not imply its approval to the exclusion of other products or vendors that also may be suitable     

Adaptation in Toxic Environments: Arsenic Genomic Islands in the Bacterial Genus Thiomonas

Acid mine drainage (AMD) is a highly toxic environment for most living organisms due to the presence of many lethal elements including arsenic (As). Thiomonas (Tm.) bacteria are found ubiquitously in AMD and can withstand these extreme conditions, in part because they are able to oxidize arsenite. In order to further improve our knowledge concerning the adaptive capacities of these bacteria, we sequenced and assembled the genome of six isolates derived from the Carnoulès AMD, and compared them to the genomes of Tm. arsenitoxydans 3As (isolated from the same site) and Tm. intermedia K12 (isolated from a sewage pipe). A detailed analysis of the Tm. sp. CB2 genome revealed various rearrangements had occurred in comparison to what was observed in 3As and K12 and over 20 genomic islands (GEIs) were found in each of these three genomes. We performed a detailed comparison of the two arsenic-related islands found in CB2, carrying the genes required for arsenite oxidation and As resistance, with those found in K12, 3As, and five other Thiomonas strains also isolated from Carnoulès (CB1, CB3, CB6, ACO3 and ACO7). Our results suggest that these arsenic-related islands have evolved differentially in these closely related Thiomonas strains, leading to divergent capacities to survive in As rich environments.     

Swine and poultry pathogens: the complete genome sequences of two strains of Mycoplasma hyopneumoniae and a strain of Mycoplasma synoviae

This work reports the results of analyses of three complete mycoplasma genomes, a pathogenic (7448) and a nonpathogenic (J) strain of the swine pathogen Mycoplasma hyopneumoniae and a strain of the avian pathogen Mycoplasma synoviae; the genome sizes of the three strains were 920,079 bp, 897,405 bp, and 799,476 bp, respectively. These genomes were compared with other sequenced mycoplasma genomes reported in the literature to examine several aspects of mycoplasma evolution. Strain-specific regions, including integrative and conjugal elements, and genome rearrangements and alterations in adhesin sequences were observed in the M. hyopneumoniae strains, and all of these were potentially related to pathogenicity. Genomic comparisons revealed that reduction in genome size implied loss of redundant metabolic pathways, with maintenance of alternative routes in different species. Horizontal gene transfer was consistently observed between M. synoviae and Mycoplasma gallisepticum. Our analyses indicated a likely transfer event of hemagglutinin-coding DNA sequences from M. gallisepticum to M. synoviae.     

Comparison of type A enterotoxigenic Staphylococcus aureus strains isolated from geographically far distant locations by pulsed field gel electrophoresis

Staphylococcal enterotoxin A (SEA) is one of the major staphylococcal enterotoxins which may cause food-borne outbreaks. In order to investigate the difference in genomic types and to elucidate the most disseminated strains for enterotoxin A-producing strains of Staphylococcus aureus , a total of 60 SEA Staph. aureus strains isolated from food and clinical samples in Taiwan and 30 strains of the same enterotoxigenic type of strains obtained from geographically far distant locations were compared for their pulsed field gel electrophoresis (PFGE) patterns. The rare cutting endonuclease Sma I generated 10 distinct genome patterns for the 60 local SEA isolates and 15 and eight genome patterns, respectively, for the 20 and 10 SEA strains originally isolated from the USA and other countries. The local isolates are less diverse in genome patterns as compared to the US isolates. Of all these PFGE patterns, a certain pattern, such as pattern 3, is shared by the food and clinical isolates and the local and foreign isolates. Thus, although SEA Staph. aureus strains from geographically far distant locations showed considerable genetic diversity, PFGE pattern 3 strain might be one of the most disseminated strains.     

Genomic Characterization of Campylobacter jejuni Strain M1

Campylobacter jejuni strain M1 (laboratory designation 99/308) is a rarely documented case of direct transmission of C. jejuni from chicken to a person, resulting in enteritis. We have sequenced the genome of C. jejuni strain M1, and compared this to 12 other C. jejuni sequenced genomes currently publicly available. Compared to these, M1 is closest to strain 81116. Based on the 13 genome sequences, we have identified the C. jejuni pan-genome, as well as the core genome, the auxiliary genes, and genes unique between strains M1 and 81116. The pan-genome contains 2,427 gene families, whilst the core genome comprised 1,295 gene families, or about two-thirds of the gene content of the average of the sequenced C. jejuni genomes. Various comparison and visualization tools were applied to the 13 C. jejuni genome sequences, including a species pan- and core genome plot, a BLAST Matrix and a BLAST Atlas. Trees based on 16S rRNA sequences and on the total gene families in each genome are presented. The findings are discussed in the background of the proven virulence potential of M1.     

Comparative DNA sequence analyses of Pyramimonas parkeae (Prasinophyceae) chloroplast genomes

Prasinophytes form a paraphyletic assemblage of early diverging green algae, which have the potential to reveal the traits of the last common ancestor of the main two green lineages: (i) chlorophyte algae and (ii) streptophyte algae. Understanding the genetic composition of prasinophyte algae is fundamental to understanding the diversification and evolutionary processes that may have occurred in both green lineages. In this study, we sequenced the chloroplast genome of Pyramimonas parkeae NIES254 and compared it with that of P. parkeae CCMP726, the only other fully sequenced P. parkeae chloroplast genome. The results revealed that P. parkeae chloroplast genomes are surprisingly variable. The chloroplast genome of NIES254 was larger than that of CCMP726 by 3,204 bp, the NIES254 large single copy was 288 bp longer, the small single copy was 5,088 bp longer, and the IR was 1,086 bp shorter than that of CCMP726. Similarity values of the two strains were almost zero in four large hot spot regions. Finally, the strains differed in copy number for three protein‐coding genes: ycf20, psaC, and ndhE. Phylogenetic analyses using 16S and 18S rDNA and rbcL sequences resolved a clade consisting of these two P. parkeae strains and a clade consisting of these plus other Pyramimonas isolates. These results are consistent with past studies indicating that prasinophyte chloroplast genomes display a higher level of variation than is commonly found among land plants. Consequently, prasinophyte chloroplast genomes may be less useful for inferring the early history of Viridiplantae than has been the case for land plant diversification.     

Illumina TruSeq Synthetic Long-Reads Empower De Novo Assembly and Resolve Complex,Highly-Repetitive Transposable Elements

High-throughput DNA sequencing technologies have revolutionized genomic analysis, including the de novo assembly of whole genomes. Nevertheless, assembly of complex genomes remains challenging, in part due to the presence of dispersed repeats which introduce ambiguity during genome reconstruction. Transposable elements (TEs) can be particularly problematic, especially for TE families exhibiting high sequence identity, high copy number, or complex genomic arrangements. While TEs strongly affect genome function and evolution, most current de novo assembly approaches cannot resolve long, identical, and abundant families of TEs. Here, we applied a novel Illumina technology called TruSeq synthetic long-reads, which are generated through highly-parallel library preparation and local assembly of short read data and which achieve lengths of 1.5–18.5 Kbp with an extremely low error rate (0.03% per base). To test the utility of this technology, we sequenced and assembled the genome of the model organism Drosophila melanogaster (reference genome strain y; cn, bw, sp) achieving an N50 contig size of 69.7 Kbp and covering 96.9% of the euchromatic chromosome arms of the current reference genome. TruSeq synthetic long-read technology enables placement of individual TE copies in their proper genomic locations as well as accurate reconstruction of TE sequences. We entirely recovered and accurately placed 4,229 (77.8%) of the 5,434 annotated transposable elements with perfect identity to the current reference genome. As TEs are ubiquitous features of genomes of many species, TruSeq synthetic long-reads, and likely other methods that generate long-reads, offer a powerful approach to improve de novo assemblies of whole genomes.     

Pot2, an inverted repeat transposon from the rice blast fungus Magnaporthe grisea

We report the cloning and characterisation of Pot2, a putative transposable element from Magnaporthe grisea. The element is 1857 by in size, has 43-bp perfect terminal inverted repeats (TIRs) and 16-bp direct repeats within the TIRs. A large open reading frame, potentially coding for a transposase-like protein, was identified. This putative protein coding region showed extensive identity to that of Fott, a transposable element from another phytopathogenic fungus, Fusarium oxysporum. Pot2, like the transposable elements Tc1 and Mariner of Caenorhabditis elegans and Drosophila, respectively, duplicates the dinucleotide TA at the target insertion site. Sequence analysis of DNA flanking 12 Pot2 elements revealed similarity to the consensus insertion sequence of Tct. Pot2 is present at a copy number of approximately 100 per haploid genome and represents one of the major repetitive DNAs shared by both rice and non-rice pathogens of M. grisea.