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

BackgroundMicrosporidian 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).ResultsOur 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.ConclusionsGenome 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.

[1]  Fabienne Thomarat,et al.  Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi , 2001, Nature.

[3]  T. Embley,et al.  Horizontal gene transfer and the evolution of parasitic protozoa. , 2003, Protist.

[4]  L. Farinelli,et al.  Gain and loss of multiple functionally related, horizontally transferred genes in the reduced genomes of two microsporidian parasites , 2012, Proceedings of the National Academy of Sciences.

[5]  Zeyang Zhou,et al.  Identification of NbME MITE families: potential molecular markers in the microsporidia Nosema bombycis. , 2010, Journal of invertebrate pathology.

[6]  T. Richards,et al.  Gene transfer: anything goes in plant mitochondria , 2010, BMC Biology.

[7]  David Haussler,et al.  Using native and syntenically mapped cDNA alignments to improve de novo gene finding , 2008, Bioinform..

[8]  M. Yamakawa,et al.  Regulation of the innate immune responses in the silkworm, Bombyx mori , 2011 .

[9]  S. Salzberg,et al.  Improved microbial gene identification with GLIMMER. , 1999, Nucleic acids research.

[10]  L. Fang,et al.  The varying microsporidian genome: existence of long-terminal repeat retrotransposon in domesticated silkworm parasite Nosema bombycis. , 2006, International journal for parasitology.

[11]  Zeyang Zhou,et al.  Characterization of a transcriptionally active Tc1-like transposon in the microsporidian Nosema bombycis , 2010, Acta Parasitologica.

[12]  Robert C. Edgar,et al.  MUSCLE: a multiple sequence alignment method with reduced time and space complexity , 2004, BMC Bioinformatics.

[13]  M. Sogin,et al.  Genes coding for reverse transcriptase, DNA-directed RNA polymerase, and chitin synthase from the microsporidian Spraguea lophii. , 1997, The Biological bulletin.

[14]  Steven J. M. Jones,et al.  Circos: an information aesthetic for comparative genomics. , 2009, Genome research.

[15]  D. Ebert,et al.  Draft genome sequence of the Daphnia pathogen Octosporea bayeri: insights into the gene content of a large microsporidian genome and a model for host-parasite interactions , 2009, Genome Biology.

[16]  Anastasios D. Tsaousis,et al.  A novel route for ATP acquisition by the remnant mitochondria of Encephalitozoon cuniculi , 2008, Nature.

[17]  L. Weiss,et al.  The Microsporidia and Microsporidiosis , 1999 .

[18]  Alexandros Stamatakis,et al.  RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models , 2006, Bioinform..

[19]  R. Moyer,et al.  Poxvirus immune modulators: functional insights from animal models. , 2002, Virus research.

[20]  P. Keeling,et al.  Bacterial Catalase in the Microsporidian Nosema locustae: Implications for Microsporidian Metabolism and Genome Evolution , 2003, Eukaryotic Cell.

[21]  L. Farinelli,et al.  Acquisition of an animal gene by microsporidian intracellular parasites , 2011, Current Biology.

[22]  Cédric Feschotte,et al.  Plant transposable elements: where genetics meets genomics , 2002, Nature Reviews Genetics.

[23]  P. Keeling,et al.  Microsporidia: biology and evolution of highly reduced intracellular parasites. , 2002, Annual review of microbiology.

[24]  L. Farinelli,et al.  The complete sequence of the smallest known nuclear genome from the microsporidian Encephalitozoon intestinalis , 2010, Nature communications.

[25]  J. Palmer,et al.  Horizontal gene transfer in eukaryotic evolution , 2008, Nature Reviews Genetics.

[26]  Shi Lei,et al.  Genomic Survey of the Non-Cultivatable Opportunistic Human Pathogen, Enterocytozoon bieneusi , 2009, PLoS pathogens.

[27]  Michael C. Schatz,et al.  Genomic Analyses of the Microsporidian Nosema ceranae, an Emergent Pathogen of Honey Bees , 2009, PLoS pathogens.

[28]  C. Slamovits,et al.  Genome Compaction and Stability in Microsporidian Intracellular Parasites , 2004, Current Biology.

[29]  R. Modigliani,et al.  Occurrence of a new microsporidan: Enterocytozoon bieneusi n.g., n. sp., in the enterocytes of a human patient with AIDS. , 1985, The Journal of protozoology.

[30]  K. Snowden.  Zoonotic microsporidia from animals and arthropods with a discussion of human infections , 2004 .

[31]  Erin E. Gill,et al.  Stripped-down DNA repair in a highly reduced parasite , 2007, BMC Molecular Biology.

[32]  Y. Kawakami,et al.  Identification of a chromosome harboring the small subunit ribosomal RNA gene of Nosema bombycis. , 1994, Journal of invertebrate pathology.

[33]  M. Borodovsky,et al.  Gene identification in novel eukaryotic genomes by self-training algorithm , 2005, Nucleic acids research.

[34]  Terry Gaasterland,et al.  DarkHorse: a method for genome-wide prediction of horizontal gene transfer , 2007, Genome Biology.

[35]  L. F. Kashkarova,et al.  [Range of the hosts of the causative agent of pébrine (Nosema bombycis) in the mulberry silkworm]. , 1980, Parazitologiia.

[36]  M. Wittner.  Historic Perspective on the Microsporidia: Expanding Horizons , 1999 .

[37]  S. Tzipori,et al.  The Reduced Genome of the Parasitic Microsporidian Enterocytozoon bieneusi Lacks Genes for Core Carbon Metabolism , 2010, Genome biology and evolution.

[38]  Daojun Cheng,et al.  A Genome-Wide Survey for Host Response of Silkworm, Bombyx mori during Pathogen Bacillus bombyseptieus Infection , 2009, PloS one.

[39]  C. Slamovits,et al.  The intriguing nature of microsporidian genomes. , 2011, Briefings in functional genomics.

[40]  R. Nielsen,et al.  Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. , 1998, Genetics.

[41]  Nick Goldman,et al.  Accuracy and Power of Statistical Methods for Detecting Adaptive Evolution in Protein Coding Sequences and for Identifying Positively Selected Sites , 2004, Genetics.

[42]  Ziheng Yang.  PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[43]  P. Keeling,et al.  Genome sequence surveys of Brachiola algerae and Edhazardia aedis reveal microsporidia with low gene densities , 2008, BMC Genomics.

[44]  L. Pasteur.  Études sur la maladie des vers a soie, moyen pratique assuré de la combattre et d'en prévenir le retour , 1870 .

[45]  Wei Qian,et al.  Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. , 2000, Molecular biology and evolution.

[46]  Hilary G. Morrison,et al.  Patterns of Genome Evolution among the Microsporidian Parasites Encephalitozoon cuniculi, Antonospora locustae and Enterocytozoon bieneusi , 2007, PloS one.

[47]  A. Lucas,et al.  Serpins, the vasculature, and viral therapeutics. , 2006, Frontiers in bioscience : a journal and virtual library.

[48]  Richard D. Emes,et al.  Duplicated Paralogous Genes Subject to Positive Selection in the Genome of Trypanosoma brucei , 2008, PloS one.

[49]  Huanming Yang,et al.  De novo assembly of human genomes with massively parallel short read sequencing. , 2010, Genome research.