Horizontal gene transfer in silkworm, Bombyx mori
© Zhu et al; licensee BioMed Central Ltd. 2011
Received: 21 October 2010
Accepted: 19 May 2011
Published: 19 May 2011
The domesticated silkworm, Bombyx mori, is the model insect for the order Lepidoptera, has economically important values, and has gained some representative behavioral characteristics compared to its wild ancestor. The genome of B. mori has been fully sequenced while function analysis of BmChi-h and BmSuc1 genes revealed that horizontal gene transfer (HGT) maybe bestow a clear selective advantage to B. mori. However, the role of HGT in the evolutionary history of B. mori is largely unexplored. In this study, we compare the whole genome of B. mori with those of 382 prokaryotic and eukaryotic species to investigate the potential HGTs.
Ten candidate HGT events were defined in B. mori by comprehensive sequence analysis using Maximum Likelihood and Bayesian method combining with EST checking. Phylogenetic analysis of the candidate HGT genes suggested that one HGT was plant-to- B. mori transfer while nine were bacteria-to- B. mori transfer. Furthermore, functional analysis based on expression, coexpression and related literature searching revealed that several HGT candidate genes have added important characters, such as resistance to pathogen, to B. mori.
Results from this study clearly demonstrated that HGTs play an important role in the evolution of B. mori although the number of HGT events in B. mori is in general smaller than those of microbes and other insects. In particular, interdomain HGTs in B. mori may give rise to functional, persistent, and possibly evolutionarily significant new genes.
The silkworm, Bombyx mori, was domesticated over 5,000 years ago and is well-known for its industrial importance in sericulture . B. mori has become a model organism for studying other Lepidoptera insects that cause serious agricultural damage and is also an important model for scientific discovery in the areas of microbiology, physiology, and genetics. However, compared to its wild ancestor, B. mori has gained some representative behavioral characteristics such as tolerance to human proximity and handling, as well as extensive crowding and lost other traits such as flight, predators, and diseases avoidance .
Horizontal gene transfer (HGT) has been not only regarded as a driving force in the innovation and evolution of genomes in prokaryotes, but also plays an important role in eukaryotes [3, 4]. In some specific eukaryotes (such as rotifers), HGT also serves as a important evolutionary impetus [5, 6]. However, in most cases, the transferred genes in insect genomes are absence of function . In contrast, Daimon et al. (2003, 2008) found that HGT maybe bestow a clear selective advantage to B. mori based on function analysis of BmChi-h and BmSuc1 genes, which have been cloned from B. mori[8, 9] The first gene contributed the potential fungi resistance to B. mori while the latter gene serves as a sugar-digesting enzyme which can degrade the alkaloidal sugar mimic glycosidase inhibitors that are toxic to B. mori. These results indicated that the HGTs in B. mori are different from those of other insects, which may play important function in the evolution of B. mori.
Some methods such as GC content and codon usage have been applied in the detection of gene transfers, but these methods have been demonstrated to be unreliable without phylogenetic analysis [10, 11]. Indeed, the detection of gene transfers is best achieved by generation of a strongly supported phylogenetic tree which contradicts the known species phylogeny [4, 10]. Luckily, the genome of B. mori has been recently fully sequenced . This provides a strong basis for the overall understanding of the existence and functions of HGTs in B. mori based on phylogenetic analysis.
In this study, the HGT events between B. mori and other unrelated species were investigated by sequence comparison of the entire predicted proteomes from B. mori with the genome sequences of 87 eukaryotic and 295 prokaryotic (Additional file 1), which represent a wide taxonomic diversity, and generated individual phylogenetic trees by different models combining with other evidences for every gene that showed a higher similarity to Arthropoda than those of other genomes. This process identified 10 putative gene transfers. Analysis of expression and coexpression revealed that these HGTs could enhance the disease resistance ability, nutrient and energy metabolism and toxin degradation. These studies give us a first glance to understand the HGTs in the evolution of B. mori by whole genome analysis.
The genome of Bombyx mori has a size of 428.7 Mb and is composed of 28 chromosomes . Gene candidates for HGT in B. mori were identified by screening the 14,623 coding genes downloaded from SilkDB .
Local database generating
In order to construct the phylogenetic trees for each sequence identified with the target taxon/similarity profile automatically and preciously, a local database (300 cores Processor BladeSystem in Zhejiang University, IBM-Biocomputing Laboratory) containing predicted protein sequences (2,385,947 coding sequences in total) from 382 species representing a wide diversity of eukaryotes and prokaryote taxa including seven insects, seven plants, 43 fungi, 30 other eukaryotes and 295 prokaryotes (Additional file 1).
The remaining candidate HGT genes were selected to search against GenBank non-redundant protein database (nr). Search strategy was the same as described above. These sequences were aligned using ClustalW , the conserved region of each alignment was trimmed using Gblocks  with a stringent settings described previously . Maximum Likelihood (ML) phylogenies were constructed by Phyml  using a JTT + Γ + I substitution model (Γ + I parameters were estimated by Phyml). The proportion of invariant sites was estimated from the data. For Bayesian phylogenies generated by MrBayes , two independent Metropolis-coupled Monte Carlo Markov Chain (MCMC) runs, each with one cold and three heated chains (heat parameter = default), were analyzed for one million generations after a burn-in of 25,000 samples and allowed mixed models of amino acid substitution. For Maximum Likelihood phylogenies, 1,000 bootstraps were performed to gain the branch support values. To be considered as an indicative of a potential HGT event, the ML branch support about B. mori should ≥ 80% while Bayesian posterior probability should ≥ 85%.
To test the support for contentious topology, we performed nonparametric branch support tests based on a Shimodaira-Hasegawa-like procedure using Tree-Puzzle . In addition, we used this software to test the statistical significance (5%) of specific topology over a collapsed version of the same branching relationship.
To test whether the transferred genes were associated with the domestication of B. mori, single nucleotide polymorphism (SNP) genelist  produced by resequencing 40 B. mori genomes was used to search if the candidate genes are on the list. Meanwhile, expression and coexpression analysis about the transferred genes were also conducted. A Pearson correlation coefficient was calculated based on microarray gene expression profiling, which were downloaded from SilkDB  to investigate the function of transferred genes analyzed in this paper. Background correction and data normalization was done by RMA algorithm in Bioconductor  and the poor annotated probe-sets were removed. Measurements for unique gene was calculated from means of the probes belong to same gene. The gene coexpressed with one of either gene in ten putative transferred genes with the absolute correlation |r| > 0.5 was further selected to analyze the function. For each transferred gene, we obtained the top 300 genes in each list sorted by absolute correlation value r referred to as coexpressed genes in descending order. The distribution of coexpressed genes with each transferred gene was mapped onto pathways. A score reflecting the extent of coexpression in each gene was assigned to each pathway by using following formula: , where G is the number of genes in a pathway present in the top 300 of all same transferred genes. R i is the rank of the ith gene of the specified pathway in a list, r i is the absolute correlation value r.
Results and Discussion
The process of HGT candidate search
HGT candidates in the Bombyx mori genome.
Top Hit organism *
SH test $
Glycerophosphoryl diester phosphodiesterase
aromatic ring-opening dioxygenase LigB subunit
In this study, stringent filters were used to identify the real HGT events. Although this procedure could identify the really HGT event undoubtedly, it may neglect several potential HGT events. First, the HGT in the common ancestor of Insecta may be considered as vertical inheritance in this study if none of the other sampled taxons has homolog. As no homologs were detected outside insects, it might be assumed that this HGT was from an extinct linkages . However, this assumption can only been proved by fossil evidence. Second, our filter could also remove genes transferred from other animal to B. mori as we only retained the genes showed a higher BLAST similarity score to a gene from any other species rather than Arthropoda (Material and Methods). However, fewer instances of potential HGTs between eukaryotes are known [4, 26], and most of these HGTs were plastids-eukaryote transfer. Finally, we only retained the transferred genes that have putative functions as functional genes may play more important role in physical and biochemical reaction in B. mori.
Comparative phylogenetic analysis, which included models that account for site rate heterogeneity and, where appropriate, comparative topology tests were used in our study. We finally identified 10 genes that was the most consistent explanation to have been transferred from the genome of bacteria to the genome of B. mori. The number of transferred genes is nearly the same as recently Acyrthosiphon pisum sequencing paper . The direction of transfer from prokaryote to B. mori is also highly suggested in our analysis. As in the phylogenetic trees, B. mori nested within a bacterial clade, and many bacterial species formed basel branches. Meanwhile, all of these 10 genes are intronless in the B. mori. Based on these results, the direction of transfer from prokaryote to B. mori is mainly confirmed. Recent findings suggested that substantial HGT between prokaryotic endosymbionts and their insect hosts were always happened [7, 27–32]. However, in this study, none of the putative transferred gene was acquired from a previous commonly reported endosymbiont (Wolbachia or Buchnera).
The HGT events in B. mori are rare but recent
It has long been described that many insects such as pea aphid, mosquitoes, and fruit flies have been influenced by HGT from Wolbachia. Nikoh et al. have suggested that nearly 30% of a Wolbachia genome was found on the X-chromosome of the insect, probably as the result of a single HGT event . In this study, we only get 10 HGTs in silkworm. Compared to previous reports, we may think that the HGT events are rare compared to other insects. The CDS of 10 putative transferred genes were used to search against wildsilkbase , which is a BLAST searchable catalogue of ESTs generated from three major wild silkworms, Antheraea assama, Samia cynthia ricini and Antheraea mylitta. We found that the transferred genes in B. mori were actually absent in three wild silkmoths. This phenomenon can be explained by two hypotheses. The first suggests that the aliened genome fragment was transferred to the B. mori chromosome after speciation of Bombycoidea, while the latter suggests an ancient HGT in the common ancestor with several losses. However, the "losses" hypothesis is highly unlikely for two reasons. It requires massive independent losses to explain the presence of these genes only in B. mori. Furthermore, this hypothesis has difficulty to explain the remarkable sequence identity between B. mori and its candidate donors (Table 1), while the general gene similarity between wild and domesticated silkworm is 58.6% based on the protein sequences of Antheraea assama deposited in Genbank. Taking into account, it appears most likely that the aliened genome fragment was transferred to the B. mori chromosome after speciation of Bombycoidea. The result suggested that these HGT events may be happened much recently. The recovery of only 10 potential HGT in B. mori suggests that HGT events happened in this species have played only a very minor role in its evolution (0.065% in the B. mori protein coding genes). Interestingly, none of the HGT events involved donor linkages from the common ancestor of Insecta. Instead, all of the HGTs involved transfer from bacteria or plant directly to the Bombycidae except chitinase, which has been demonstrated to be an ancient HGT event . It is worth noting that half of these transfers originate from the Enterobacteriaceae species. Many Enterobacteriaceae species have been demonstrated to be as symbionts of insects . It may facilitate to exchange the genetic material between symbionts and hosts.
Putative functional assignment of 10 HGTs
The potential identities and biological functions of the 10 HGT candidates were investigated based on the similarity of their putatively encoded products to known proteins. We found that the 10 putative transferred genes were all enzymes, which suggested a strong functional trend among these HGTs in B. mori. Furthermore, hypergeometric distribution analysis of the candidate HGT genes indicated catalytic activity and hydrolase activity genes are significantly overrepresented (with p-values 4.6e-05 and 6.4e-06, respectively). Compared with our study, none of the putative functional gene was identified in Aphid genome sequence analysis . The acquisition of such genes by HGT may therefore have been selected by a competitive advantage in the recipient to degrade or take up substrates with greater efficiency.
A putative glycerophosphoryl diester phosphodiesterase (BGIBMGA007767, pfam03009, e-value 1e-52) was encoded by an HGT candidate from the Pseudomonas species to B. mori (Additional file 3). This enzyme display broad specificity for glycerophosphodiesters, glycerophosphocholine, glycerophosphoethanolamine, glycerophosphoglycerol and bis (glycerophosphoglycerol), all of which are hydrolysed by this enzyme. This gene also has been reported to be involved in HGT between spider and bacteria . Expression analysis suggested that this gene was not only highly expressed in integument and fat body but also in midgut and hemocyte. Coexpression analysis indicated that the function of this gene may be consistent with BGIBMGA010285.
Two genes, encoded chitinase and sucrose-6-phosphate hydrolase, are the only two genes which have been reported to be transferred from bacteria to B. mori previously [8, 9] were also identified in this study. B. mori acquired the first gene from a bacterium or an ancestral baculovirus and the latter gene from a bacterium [8, 9]. The results suggested that both genes bestow a clear selective advantage to B. mori. The discovery of both genes indicated that our strategy in searching candidate HGTs in B. mori is reliable.
Another candidate HGT was a gene transferred from Edwardsiella tarda (Additional file 4). Homology analysis suggested that this gene was a glucose-1-phosphatase/inositol phosphatase (BGIBMGA011204, pfam00328, e-value 1.4e-17). This enzyme has been reported to enhance the phytase activity in bacteria . Expression pattern revealed that this gene is extremely highly expressed in head (Figure 2). Coexpression analysis suggested that the function of this gene may be most associated with Butirosin and neomycin biosynthesis (Additional file 2). As the gene in this pathway is involved in aminoglycoside antibiotics biosynthesis, which is effective against several types of bacteria . The result indicated that this gene may be associated with the resistance of B. mori to bacterial pathogen.
The sixth HGT candidate (Additional file 5) showed sequence similarity to alginate lyase (BGIBMGA005615, pfam08240, e-value 2.1e-08), which can catalyze the degradation of alginate by a β-elimination mechanism . In bacteria, this enzyme has been reported to destroy the cell detachment, which is the key process in biofilm formation . Unfortunately, there is no corresponding probe to detect this gene in microarray data.
The seventh, eighth and ninth HGT are putative kynureninase (BGIBMGA007146, pfam00266, e-value 2.6e-18), N-methyltryptophan oxidase (BGIBMGA008215, pfam01266, e-value 4.2e-54) and gamma- glutamyltranspeptidase (pfam01019, e-value 4.3e-125), respectively. The seventh gene encoded by an HGT candidate from Listeria grayi to B. mori (Additional file 6). The mutation of this gene has been reported to be associated with abnormal body coloration in B. mori. The eighth and ninth genes were both encoded by an HGT candidate from Enterobacteriaceae to B. mori (Additional file 7 and 8). In both genes, the eighth one has been reported to be involved in substrate specificity determination  while the ninth gene plays a key role in the gamma-glutamyl cycle, a pathway for the synthesis and degradation of glutathione and drug and xenobiotic detoxification . All of the three genes are highly expressed in malpighian tubule (Figure 2), and the coexpression patterns are similar (Additional file 2). The results indicated that these genes may be associated with valine, leucine and isoleucine degradation pathway (Additional file 2). It has been reported that after B. mori was infected by its bacterial pathogen Bacillus bombyseptieus, the expression level of the genes which are associated with valine, leucine and isoleucine degradation pathway were more than three times the baseline levels in the malpighian tubule. Huang et al. (2009) speculated that these genes might be involved in the detoxification of malpighian tubules . The results suggested that these genes may be associated with the resistance of B. mori to bacterial pathogen.
The tenth gene which encoded as aromatic ring-opening dioxygenase LigB subunit (BGIBMGA003842, pfam02900, e-value 1.2e-51) was transferred from a plant donor (Additional file 9). No tissue specific expression pattern was observed in this gene (Figure 2). However, coexpression analysis suggested that this gene may be most associated with Butirosin and neomycin biosynthesis (Additional file 2). It has been reported that Butirosin and neomycin belong to a family of clinically valuable 2-deoxystreptamine (2DOS)-containing aminoglycoside antibiotics . The result indicated that this gene may be associated with the resistance of B. mori to bacterial pathogen.
In this study, the coexpressed genes involved in HGT are mapped onto biological pathway to investigate the potential function in B. mori. This scheme provides a comprehensive way to identify previously unknown functional patterns in sets of genes with known functions. It has been well known that the domesticated B. mori has gained or lost some representative behavioral characteristics compared to its wild ancestor. A lot of researches have focused on the search for 'domestication genes' in mammals by the polygenic nature of most behavioral traits . Indeed, the linking of coat-color genes in mammals to brain biochemistry and behaviors has been found to be able to facilitate domestication and ease of handling . In our study, SNP data searching indicated that all of the 10 genes have no direct relationship with human domestication although all of these genes are suggested to be transferred recently after the divergence from the wild silkworm. Therefore, we could suggest that HGT in B. mori maybe mainly depend on its nature requirement, but not on human activity. However, our study strongly suggested that these transferred genes increase the ability of B. mori to natural survival and adaption.
In this study, we conclude that HGT is both rare and recent between B. mori and other unrelated species based on large-scale genome sequence data analysis using a strict and highly conservative set of ML and Bayesian methods combining with other analysis such as EST checking for inferring HGTs. Results from this study clearly indicated that HGT in B. mori has occurred and may have provided advantageous gene functions that could enhance the disease resistance ability, nutrient and energy metabolism and toxin degradation.
This project was supported by Zhejiang Provincial Natural Science Foundation of China (Y3090150), Zhejiang Provincial Project (2010R10091), the Fundamental Research Funds for the Central Universities, Specialized Research Fund for the Doctoral Program of Higher Education (20090101120083) and the Agricultural Ministry of China (nyhyzx 201003029; 201003066),.
- Goldsmith MR, Shimada T, Abe H: The genetics and genomics of the silkworm, Bombyx mori. Annu Rev Entomol. 2005, 50: 71-100. 10.1146/annurev.ento.50.071803.130456.View ArticlePubMedGoogle Scholar
- Xia QY, Guo YR, Zhang Z, Li D, Xuan ZL, Li Z, Dai FY, Li YR, Cheng DJ, Li RQ, et al: Complete Resequencing of 40 Genomes Reveals Domestication Events and Genes in Silkworm (Bombyx). Science. 2009, 326: 433-436. 10.1126/science.1176620.View ArticlePubMedPubMed CentralGoogle Scholar
- Boucher Y, Douady CJ, Papke RT, Walsh DA, Boudreau MER, Nesbo CL, Case RJ, Doolittle WF: Lateral gene transfer and the origins of prokaryotic groups. Annu Rev Genet. 2003, 37: 283-328. 10.1146/annurev.genet.37.050503.084247.View ArticlePubMedGoogle Scholar
- Keeling PJ, Palmer JD: Horizontal gene transfer in eukaryotic evolution. Nat Rev Genet. 2008, 9: 605-618. 10.1038/nrg2386.View ArticlePubMedGoogle Scholar
- Kaessmann H: Origins, evolution, and phenotypic impact of new genes. Genome research. 2010, 20: 1313-1326. 10.1101/gr.101386.109.View ArticlePubMedPubMed CentralGoogle Scholar
- Gladyshev EA, Meselson M, Arkhipova IR: Massive horizontal gene transfer in bdelloid rotifers. Science. 2008, 320: 1210-1213. 10.1126/science.1156407.View ArticlePubMedGoogle Scholar
- Nikoh N, McCutcheon JP, Kudo T, Miyagishima S, Moran NA, Nakabachi A: Bacterial Genes in the Aphid Genome: Absence of Functional Gene Transfer from Buchnera to Its Host. PLoS Genet. 2010, 6:Google Scholar
- Daimon T, Taguchi T, Meng Y, Katsuma S, Mita K, Shimada T: Beta-fructofuranosidase genes of the silkworm, Bombyx mori - Insights into enzymatic adaptation of B. mori to toxic alkaloids in mulberry latex. J Biol Chem. 2008, 283: 15271-15279. 10.1074/jbc.M709350200.View ArticlePubMedPubMed CentralGoogle Scholar
- Daimon T, Hamada K, Mita K, Okano K, Suzuki MG, Kobayashi M, Shimada T: A Bombyx mori gene, BmChi-h, encodes a protein homologous to bacterial and baculovirus chitinases. Insect Biochem Mol Biol. 2003, 33: 749-759. 10.1016/S0965-1748(03)00084-5.View ArticlePubMedGoogle Scholar
- Richards TA, Soanes DM, Foster PG, Leonard G, Thomton CR, Talbot NJ: Phylogenomic Analysis Demonstrates a Pattern of Rare and Ancient Horizontal Gene Transfer between Plants and Fungi. Plant Cell. 2009, 21: 1897-1911. 10.1105/tpc.109.065805.View ArticlePubMedPubMed CentralGoogle Scholar
- Ragan MA: On surrogate methods for detecting lateral gene transfer. FEMS Microbiol Lett. 2001, 201: 187-191. 10.1111/j.1574-6968.2001.tb10755.x.View ArticlePubMedGoogle Scholar
- Xia QY, Zhou ZY, Lu C, Cheng DJ, Dai FY, Li B, Zhao P, Zha XF, Cheng TC, Chai CL, et al: A draft sequence for the genome of the domesticated silkworm (Bombyx mori). Science. 2004, 306: 1937-1940.View ArticlePubMedGoogle Scholar
- Wang J, Xia QY, He XM, Dai MT, Ruan J, Chen J, Yu G, Yuan HF, Hu YF, Li RQ, et al: SilkDB: a knowledgebase for silkworm biology and genomics. Nucleic Acids Res. 2005, 33: D399-D402.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.View ArticlePubMedPubMed CentralGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: ClustalW - improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.View ArticlePubMedPubMed CentralGoogle Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17: 540-552.View ArticlePubMedGoogle Scholar
- Ciccarelli FD, Doerks T, von Mering C, Creevey CJ, Snel B, Bork P: Toward automatic reconstruction of a highly resolved tree of life. Science. 2006, 311: 1283-1287. 10.1126/science.1123061.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.View ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge YC, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5:Google Scholar
- Li L, Stoeckert CJ, Roos DS: OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13: 2178-2189. 10.1101/gr.1224503.View ArticlePubMedPubMed CentralGoogle Scholar
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17: 1246-1247. 10.1093/bioinformatics/17.12.1246.View ArticlePubMedGoogle Scholar
- Mita K, Morimyo M, Okano K, Koike Y, Nohata J, Kawasaki H, Kadono-Okuda K, Yamamoto K, Suzuki MG, Shimada T, et al: The construction of an EST database for Bombyx mori and its application. Proc Natl Acad Sci USA. 2003, 100: 14121-14126. 10.1073/pnas.2234984100.View ArticlePubMedPubMed CentralGoogle Scholar
- Fournier GP, Huang JL, Gogarten JP: Horizontal gene transfer from extinct and extant lineages: biological innovation and the coral of life. Philos Trans R Soc Lond B Biol Sci. 2009, 364: 2229-2239. 10.1098/rstb.2009.0033.View ArticlePubMedPubMed CentralGoogle Scholar
- Andersson JO: Lateral gene transfer in eukaryotes. Cell Mol Life Sci. 2005, 62: 1182-1197. 10.1007/s00018-005-4539-z.View ArticlePubMedGoogle Scholar
- Woolfit M, Iturbe-Ormaetxe I, McGraw EA, O'Neill SL: An Ancient Horizontal Gene Transfer between Mosquito and the Endosymbiotic Bacterium Wolbachia pipientis. Mol Biol Evol. 2009, 26: 367-374. 10.1093/molbev/msn253.View ArticlePubMedGoogle Scholar
- Kondo N, Nikoh N, Ijichi N, Shimada M, Fukatsu T: Genome fragment of Wolbachia endosymbiont transferred to X chromosome of host insect. Proc Natl Acad Sci USA. 2002, 99: 14280-14285. 10.1073/pnas.222228199.View ArticlePubMedPubMed CentralGoogle Scholar
- Hotopp JCD, Clark ME, Oliveira D, Foster JM, Fischer P, Torres MC, Giebel JD, Kumar N, Ishmael N, Wang SL, et al: Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science. 2007, 317: 1753-1756. 10.1126/science.1142490.View ArticleGoogle Scholar
- Nikoh N, Tanaka K, Shibata F, Kondo N, Hizume M, Shimada M, Fukatsu T: Wolbachia genome integrated in an insect chromosome: Evolution and fate of laterally transferred endosymbiont genes. Genome Res. 2008, 18: 272-280. 10.1101/gr.7144908.View ArticlePubMedPubMed CentralGoogle Scholar
- Klasson L, Kambris Z, Cook PE, Walker T, Sinkins SP: Horizontal gene transfer between Wolbachia and the mosquito Aedes aegypti. BMC Genomics. 2009, 10:Google Scholar
- Nikoh N, Nakabachi A: Aphids acquired symbiotic genes via lateral gene transfer. BMC Biol. 2009, 7:Google Scholar
- Hotopp D, Julie C: Horizontal gene transfer between bacteria and animals. Trends Genet. 2011, 27: 157-163. 10.1016/j.tig.2011.01.005.View ArticleGoogle Scholar
- Arunkumar KP, Tomar A, Daimon T, Shimada T, Nagaraju J: WildSilkbase: An EST database of wild silkmoths. BMC Genomics. 2008, 9:Google Scholar
- Moran NA, Russell JA, Koga R, Fukatsu T: Evolutionary relationships of three new species of Enterobacteriaceae living as symbionts of aphids and other insects. Appl Environ Microbiol. 2005, 71: 3302-3310. 10.1128/AEM.71.6.3302-3310.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Anderson I, Sorokin A, Kapatral V, Reznik G, Bhattacharya A, Mikhailova N, Burd H, Joukov V, Kaznadzey D, Walunas T, et al: Comparative genome analysis of Bacillus cereus group genomes with Bacillus subtilis. FEMS Microbiol Lett. 2005, 250: 175-184. 10.1016/j.femsle.2005.07.008.View ArticlePubMedGoogle Scholar
- Nosenko T, Bhattacharya D: Horizontal gene transfer in chromalveolates. BMC Evol Biol. 2007, 7:Google Scholar
- Arrese EL, Soulages JL: Insect Fat Body: Energy, Metabolism, and Regulation. Annu Rev Entomol. 2010, 55: 207-225. 10.1146/annurev-ento-112408-085356.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou ZH, Yang HJ, Chen M, Lou CF, Zhang YZ, Chen KP, Wang Y, Yu ML, Yu F, Li JY, Zhong BX: Comparative Proteomic Analysis between the Domesticated Silkworm (Bombyx mori) Reared on Fresh Mulberry Leaves and on Artificial Diet. J Proteome Res. 2008, 7: 5103-5111. 10.1021/pr800383r.View ArticlePubMedGoogle Scholar
- Cordes MHJ, Binford GJ: Lateral gene transfer of a dermonecrotic toxin between spiders and bacteria. Bioinformatics. 2006, 22: 264-268. 10.1093/bioinformatics/bti811.View ArticlePubMedGoogle Scholar
- Herter T, Berezina OV, Zinin NV, Velikodvorskaya GA, Greiner R, Borriss R: Glucose-1-phosphatase (AgpE) from Enterobacter cloacae displays enhanced phytase activity. Appl Microbiol Biotechnol. 2006, 70: 60-64. 10.1007/s00253-005-0024-8.View ArticlePubMedGoogle Scholar
- Kudo F, Eguchi T: Biosynthetic genes for aminoglycoside antibiotics. J Antibiot. 2009, 62: 471-481. 10.1038/ja.2009.76.View ArticlePubMedGoogle Scholar
- Wong TY, Preston LA, Schiller NL: Alginate lyase: Review of major sources and enzyme characteristics, structure-function analysis, biological roles, and applications. Annu Rev Microbiol. 2000, 54: 289-340. 10.1146/annurev.micro.54.1.289.View ArticlePubMedGoogle Scholar
- Boyd A, Chakrabarty AM: RRole of alginate lyase in cell detachment of Pseudomonas aeruginosa. Appl Environ Microbiol. 1994, 60: 2355-2359.PubMedPubMed CentralGoogle Scholar
- Meng Y, Katsuma S, Mita K, Shimada T: Abnormal red body coloration of the silkworm, Bombyx mori, is caused by a mutation in a novel kynureninase. Genes Cells. 2009, 14: 129-140. 10.1111/j.1365-2443.2008.01257.x.View ArticlePubMedGoogle Scholar
- Ilari A, Bonamore A, Franceschini S, Fiorillo A, Boffi A, Colotti G: The X-ray structure of N-methyltryptophan oxidase reveals the structural determinants of substrate specificity. Proteins Struct Funct Bioinformat. 2008, 71: 2065-2075. 10.1002/prot.21898.View ArticleGoogle Scholar
- Courtay C, Oster T, Michelet F, Visvikis A, Diederich M, Wellman M, Siest G: Gamma-glutamyl-transferase - nucleotide-sequence of the human pancreatic cdna - evidence for a ubiquitous gamma-glutamyl-transferase polypeptide in human tissues. Biochem Pharmacol. 1992, 43: 2527-2533. 10.1016/0006-2952(92)90140-E.View ArticlePubMedGoogle Scholar
- Huang LL, Cheng TC, Xu PZ, Cheng DJ, Fang T, Xia QY: A Genome-Wide Survey for Host Response of Silkworm, Bombyx mori during Pathogen Bacillus bombyseptieus Infection. PLoS One. 2009, 4:Google Scholar
- Kudo F, Eguchi T: Biosynthetic enzymes for the aminoglycosides butirosin and neomycin. Methods Enzymol. 2009, 459: 493-591.View ArticlePubMedGoogle Scholar
- Price EO: Animal domestication and behavior. 2002, Wallingford Oxon: Cab InternationalView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.