Sequence Analysis of a Novel Insertion Site of Transposon IS10

A sacB mutant was obtained by transposon IS10 inactivation of a plasmid pXT3sacB carrying the sacB gene. Sequencing of this mutant plasmid DNA (GenBank accession No. AY580883.1) showed that the IS10 flanking the 22 bp inverted repeats were 5′-CTGAGAGATCCCCTCATAATTT-3′ and 5′-AAATCATTAGGGGATTCATCAG-3′, which were the similar to those published in reports previously. However, the target sequence adjacent to IS10 was 5′-TGCTTGGTT-3′ instead of the previously reported 5′-NGCTNAGCN-3′. To our knowledge, this is the first report on the novel insertion site of IS10. In addition, Southern blot hybridization confirmed that the mobile IS10 originated from the chromosomal DNA of the host strain Escherichia coli DH5α and that there were two copies in the DH5α genome.     

Caenorhabditis elegans: A Genetic Guide to Parasitic Nematode Biology

The advent of parasite genome sequencing projects, as well as an increase in biology-directed gene discovery, promises to reveal genes encoding many of the key molecules required for nematode-host interactions. However, distinguishing parasitism genes from those merely required for nematode viability remains a substantial challenge. Although this will ultimately require a functional test in the host or parasite, the free-living nematode Caenorhabditis elegans can be exploited as a heterologous system to determine function of candidate parasitism genes. Studies of C. elegans also have revealed genetic networks, such as the dauer pathway, that may also be important adaptations for parasitism. As a more directed means of identifying parasitism traits, we developed classical genetics for Heterodera glycines and have used this approach to map genes conferring host resistance-breaking phenotypes. It is likely that the C. elegans and H. glycines genomes will be at least partially syntenic, thus permitting predictive physical mapping of H. glycines genes of interest.     

Genomic sequence of the aflatoxigenic filamentous fungus Aspergillus nomius

Background

Aspergillus nomius is an opportunistic pathogen and one of the three most important producers of aflatoxins in section Flavi. This fungus has been reported to contaminate agricultural commodities, but it has also been sampled in non-agricultural areas so the host range is not well known. Having a similar mycotoxin profile as A. parasiticus, isolates of A. nomius are capable of secreting B- and G- aflatoxins.

Results

In this study we discovered that the A. nomius type strain (NRRL 13137) has a genome size of approximately 36 Mb which is comparable to other Aspergilli whose genomes have been sequenced. Its genome encompasses 11,918 predicted genes, 72 % of which were assigned GO terms using BLAST2GO. More than 1,200 of those predicted genes were identified as unique to A. nomius, and the most significantly enriched GO category among the unique genes was oxidoreducatase activity. Phylogenomic inference shows NRRL 13137 as ancestral to the other aflatoxigenic species examined from section Flavi. This strain contains a single mating-type idiomorph designated as MAT1-1.

Conclusions

This study provides a preliminary analysis of the A. nomius genome. Given the recently discovered potential for A. nomius to undergo sexual recombination, and based on our findings, this genome sequence provides an additional evolutionary reference point for studying the genetics and biology of aflatoxin production.     

Genomics and genetics of Sulfolobus islandicus LAL14/1, a model hyperthermophilic archaeon

The 2 465 177 bp genome of Sulfolobus islandicus LAL14/1, host of the model rudivirus SIRV2, was sequenced. Exhaustive comparative genomic analysis of S. islandicus LAL14/1 and the nine other completely sequenced S. islandicus strains isolated from Iceland, Russia and USA revealed a highly syntenic common core genome of approximately 2 Mb and a long hyperplastic region containing most of the strain-specific genes. In LAL14/1, the latter region is enriched in insertion sequences, CRISPR (clustered regularly interspaced short palindromic repeats), glycosyl transferase genes, toxin–antitoxin genes and MITE (miniature inverted-repeat transposable elements). The tRNA genes of LAL14/1 are preferential targets for the integration of mobile elements but clusters of atypical genes (CAG) are also integrated elsewhere in the genome. LAL14/1 carries five CRISPR loci with 10 per cent of spacers matching perfectly or imperfectly the genomes of archaeal viruses and plasmids found in the Icelandic hot springs. Strikingly, the CRISPR_2 region of LAL14/1 carries an unusually long 1.9 kb spacer interspersed between two repeat regions and displays a high similarity to pING1-like conjugative plasmids. Finally, we have developed a genetic system for S. islandicus LAL14/1 and created ΔpyrEF and ΔCRISPR_1 mutants using double cross-over and pop-in/pop-out approaches, respectively. Thus, LAL14/1 is a promising model to study virus–host interactions and the CRISPR/Cas defence mechanism in Archaea.     

When parasitic wasps hijacked viruses: genomic and functional evolution of polydnaviruses

The Polydnaviridae (PDV), including the Bracovirus (BV) and Ichnovirus genera, originated from the integration of unrelated viruses in the genomes of two parasitoid wasp lineages, in a remarkable example of convergent evolution. Functionally active PDVs represent the most compelling evolutionary success among endogenous viral elements (EVEs). BV evolved from the domestication by braconid wasps of a nudivirus 100 Ma. The nudivirus genome has become an EVE involved in BV particle production but is not encapsidated. Instead, BV genomes have co-opted virulence genes, used by the wasps to control the immunity and development of their hosts. Gene transfers and duplications have shaped BV genomes, now encoding hundreds of genes. Phylogenomic studies suggest that BVs contribute largely to wasp diversification and adaptation to their hosts. A genome evolution model explains how multidirectional wasp adaptation to different host species could have fostered PDV genome extension. Integrative studies linking ecological data on the wasp to genomic analyses should provide new insights into the adaptive role of particular BV genes. Forthcoming genomic advances should also indicate if the associations between endoparasitoid wasps and symbiotic viruses evolved because of their particularly intimate interactions with their hosts, or if similar domesticated EVEs could be uncovered in other parasites.     

Molecular Transfer of Nematode Resistance Genes

Recombinant DNA techniques have been used to introduce agronomically valuable traits, including resistance to viruses, herbicides, and insects, into crop plants. Introduction of these genes into plants frequently involves Agrobacterium-mediated gene transfer. The potential exists for applying this technology to nematode control by introducing genes conferring resistance to nematodes. Transferred genes could include those encoding products detrimental to nematode development or reproduction as well as cloned host resistance genes. Host genes that confer resistance to cyst or root-knot nematode species have been identified in many plants. The best characterized is Mi, a gene that confers resistance to root-knot nematodes in tomato. A map-based cloning approach is being used to isolate the gene. For development of a detailed map of the region of the genome surrounding Mi, DNA markers genetically linked to Mi have been identified and analyzed in tomato lines that have undergone a recombination event near Mi. The molecular map will be used to identify DNA corresponding to Mi. We estimate that a clone of Mi will be obtained in 2-5 years. An exciting prospect is that introduction of this gene will confer resistance in plant species without currently available sources of resistance.     

A parasitic selfish gene that affects host promiscuity

Selfish genes demonstrate transmission bias and invade sexual populations despite conferring no benefit to their hosts. While the molecular genetics and evolutionary dynamics of selfish genes are reasonably well characterized, their effects on hosts are not. Homing endonuclease genes (HEGs) are one well-studied family of selfish genes that are assumed to be benign. However, we show that carrying HEGs is costly for Saccharomyces cerevisiae, demonstrating that these genetic elements are not necessarily benign but maybe parasitic. We estimate a selective load of approximately 1–2% in ‘natural’ niches. The second aspect we examine is the ability of HEGs to affect hosts'' sexual behaviour. As all selfish genes critically rely on sex for spread, then any selfish gene correlated with increased host sexuality will enjoy a transmission advantage. While classic parasites are known to manipulate host behaviour, we are not aware of any evidence showing a selfish gene is capable of affecting host promiscuity. The data presented here show a selfish element may increase the propensity of its eukaryote host to undergo sex and along with increased rates of non-Mendelian inheritance, this may counterbalance mitotic selective load and promote spread. Demonstration that selfish genes are correlated with increased promiscuity in eukaryotes connects with ideas suggesting that selfish genes promoted the evolution of sex initially.     

Single-Nucleotide-Specific Targeting of the Tf1 Retrotransposon Promoted by the DNA-Binding Protein Sap1 of Schizosaccharomyces pombe

Transposable elements (TEs) constitute a substantial fraction of the eukaryotic genome and, as a result, have a complex relationship with their host that is both adversarial and dependent. To minimize damage to cellular genes, TEs possess mechanisms that target integration to sequences of low importance. However, the retrotransposon Tf1 of Schizosaccharomyces pombe integrates with a surprising bias for promoter sequences of stress-response genes. The clustering of integration in specific promoters suggests that Tf1 possesses a targeting mechanism that is important for evolutionary adaptation to changes in environment. We report here that Sap1, an essential DNA-binding protein, plays an important role in Tf1 integration. A mutation in Sap1 resulted in a 10-fold drop in Tf1 transposition, and measures of transposon intermediates support the argument that the defect occurred in the process of integration. Published ChIP-Seq data on Sap1 binding combined with high-density maps of Tf1 integration that measure independent insertions at single-nucleotide positions show that 73.4% of all integration occurs at genomic sequences bound by Sap1. This represents high selectivity because Sap1 binds just 6.8% of the genome. A genome-wide analysis of promoter sequences revealed that Sap1 binding and amounts of integration correlate strongly. More important, an alignment of the DNA-binding motif of Sap1 revealed integration clustered on both sides of the motif and showed high levels specifically at positions +19 and −9. These data indicate that Sap1 contributes to the efficiency and position of Tf1 integration.     

Dissecting the fungal biology of Bipolaris papendorfii: from phylogenetic to comparative genomic analysis

Bipolaris papendorfii has been reported as a fungal plant pathogen that rarely causes opportunistic infection in humans. Secondary metabolites isolated from this fungus possess medicinal and anticancer properties. However, its genetic fundamental and basic biology are largely unknown. In this study, we report the first draft genome sequence of B. papendorfii UM 226 isolated from the skin scraping of a patient. The assembled 33.4 Mb genome encodes 11,015 putative coding DNA sequences, of which, 2.49% are predicted transposable elements. Multilocus phylogenetic and phylogenomic analyses showed B. papendorfii UM 226 clustering with Curvularia species, apart from other plant pathogenic Bipolaris species. Its genomic features suggest that it is a heterothallic fungus with a putative unique gene encoding the LysM-containing protein which might be involved in fungal virulence on host plants, as well as a wide array of enzymes involved in carbohydrate metabolism, degradation of polysaccharides and lignin in the plant cell wall, secondary metabolite biosynthesis (including dimethylallyl tryptophan synthase, non-ribosomal peptide synthetase, polyketide synthase), the terpenoid pathway and the caffeine metabolism. This first genomic characterization of B. papendorfii provides the basis for further studies on its biology, pathogenicity and medicinal potential.     

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.     

Virulence related sequences; insights provided by comparative genomics of Streptococcus uberis of differing virulence

Background

Streptococcus uberis, a Gram-positive, catalase-negative member of the family Streptococcaceae is an important environmental pathogen responsible for a significant proportion of subclinical and clinical bovine intramammary infections. Currently, the genome of only a single reference strain (0140J) has been described. Here we present a comparative analysis of complete draft genome sequences of an additional twelve S. uberis strains.

Results

Pan and core genome analysis revealed the core genome common to all strains to be 1,550 genes in 1,509 orthologous clusters, complemented by 115-246 accessory genes present in one or more S. uberis strains but absent in the reference strain 0140J. Most of the previously predicted virulent genes were present in the core genome of all 13 strains but gene gain/loss was observed between the isolates in CDS associated with clustered regularly interspaced short palindromic repeats (CRISPRs), prophage and bacteriocin production. Experimental challenge experiments confirmed strain EF20 as non-virulent; only able to infect in a transient manner that did not result in clinical mastitis. Comparison of the genome sequence of EF20 with the validated virulent strain 0140J identified genes associated with virulence, however these did not relate clearly with clinical/non-clinical status of infection.

Conclusion

The gain/loss of mobile genetic elements such as CRISPRs and prophage are a potential driving force for evolutionary change. This first “whole-genome” comparison of strains isolated from clinical vs non-clinical intramammary infections including the type virulent vs non-virulent strains did not identify simple gene gain/loss rules that readily explain, or be confidently associated with, differences in virulence. This suggests that a more complex dynamic determines infection potential and clinical outcome not simply gene content.

Electronic supplementary material

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

Draft Genome Sequencing and Comparative Analysis of Aspergillus sojae NBRC4239

We conducted genome sequencing of the filamentous fungus Aspergillus sojae NBRC4239 isolated from the koji used to prepare Japanese soy sauce. We used the 454 pyrosequencing technology and investigated the genome with respect to enzymes and secondary metabolites in comparison with other Aspergilli sequenced. Assembly of 454 reads generated a non-redundant sequence of 39.5-Mb possessing 13 033 putative genes and 65 scaffolds composed of 557 contigs. Of the 2847 open reading frames with Pfam domain scores of >150 found in A. sojae NBRC4239, 81.7% had a high degree of similarity with the genes of A. oryzae. Comparative analysis identified serine carboxypeptidase and aspartic protease genes unique to A. sojae NBRC4239. While A. oryzae possessed three copies of α-amyalse gene, A. sojae NBRC4239 possessed only a single copy. Comparison of 56 gene clusters for secondary metabolites between A. sojae NBRC4239 and A. oryzae revealed that 24 clusters were conserved, whereas 32 clusters differed between them that included a deletion of 18 508 bp containing mfs1, mao1, dmaT, and pks-nrps for the cyclopiazonic acid (CPA) biosynthesis, explaining the no productivity of CPA in A. sojae. The A. sojae NBRC4239 genome data will be useful to characterize functional features of the koji moulds used in Japanese industries.     

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus,Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.     

Complete Mitochondrial Genome of Anoplocephala magna Solidifying the Species

The 2 species of the genus Anoplocephala (Anoplocephalidae), A. perfoliata and A. magna, are among the most important equine cestode parasites. However, there is little information about their differences at the molecular level. The present study revealed that the mitochondrial (mt) genome of A. magna was 13,759 bp in size and 700 bp shorter than that of A. perfoliata. The 2 species includes 2 rRNA, 22 tRNA, and 12 protein-coding genes each. The size of each of the 36 genes was the same as that of A. perfoliata, except for cox1, rrnL, trnC, trnS2(UCN), trnG, trnH, trnQ, and trnP. In the full mitochondrial genome, the sequence similarity was 87.1%. The divergence in the nucleotide and amino acid sequences of individual protein-coding genes ranged from 11.1% to 16% and 6.8% to 16.4%, respectively. The 2 noncoding regions of the mt genome of A. magna were 199 bp and 271 bp in length, while the equivalent regions in A. perfoliata were 875 bp and 276 bp, respectively. The results of this study support the proposal that A. magna and A. perfoliata are separate species, consistent with previous morphological analyses.