Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

Comparative genomic analyses of Mycoplasma hyopneumoniae pathogenic 168 strain and its high-passaged attenuated strain

Wei Liu1, Shaobo Xiao1, Mao Li1, Shaohua Guo1, Sha Li1, Rui Luo1, Zhixin Feng2, Bin Li2, Zhemin Zhou3, Guoqing Shao2, Huanchun Chen1 and Liurong Fang1*

Author affiliations

1 Division of Animal Infectious Diseases, State Key Laboratory of Agricultural Microbiology, College of Veterinary Medicine, Huazhong Agricultural University, Wuhan, 430070, People’s Republic of China

2 Institute of Veterinary Medicine, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, People’s Republic of China

3 Environmental Research Institute, University College Cork, Cork, Ireland

For all author emails, please log on.

Citation and License

BMC Genomics 2013, 14:80  doi:10.1186/1471-2164-14-80


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/14/80


Received:28 October 2012
Accepted:31 January 2013
Published:5 February 2013

© 2013 Liu et al.; licensee 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.

Abstract

Background

Mycoplasma hyopneumoniae is the causative agent of porcine enzootic pneumonia (EP), a mild, chronic pneumonia of swine. Despite presenting with low direct mortality, EP is responsible for major economic losses in the pig industry. To identify the virulence-associated determinants of M. hyopneumoniae, we determined the whole genome sequence of M. hyopneumoniae strain 168 and its attenuated high-passage strain 168-L and carried out comparative genomic analyses.

Results

We performed the first comprehensive analysis of M. hyopneumoniae strain 168 and its attenuated strain and made a preliminary survey of coding sequences (CDSs) that may be related to virulence. The 168-L genome has a highly similar gene content and order to that of 168, but is 4,483 bp smaller because there are 60 insertions and 43 deletions in 168-L. Besides these indels, 227 single nucleotide variations (SNVs) were identified. We further investigated the variants that affected CDSs, and compared them to reported virulence determinants. Notably, almost all of the reported virulence determinants are included in these variants affected CDSs. In addition to variations previously described in mycoplasma adhesins (P97, P102, P146, P159, P216, and LppT), cell envelope proteins (P95), cell surface antigens (P36), secreted proteins and chaperone protein (DnaK), mutations in genes related to metabolism and growth may also contribute to the attenuated virulence in 168-L. Furthermore, many mutations were located in the previously described repeat motif, which may be of primary importance for virulence.

Conclusions

We studied the virulence attenuation mechanism of M. hyopneumoniae by comparative genomic analysis of virulent strain 168 and its attenuated high-passage strain 168-L. Our findings provide a preliminary survey of CDSs that may be related to virulence. While these include reported virulence-related genes, other novel virulence determinants were also detected. This new information will form the foundation of future investigations into the pathogenesis of M. hyopneumoniae and facilitate the design of new vaccines.

Keywords:
Mycoplasma hyopneumoniae; Genetic variation; Virulence attenuation; Sequence analysis; Repetitive sequences; Virulence factors

Background

Mycoplasma hyopneumoniae causes porcine enzootic pneumonia, which is a mild, chronic pneumonia of swine [1]. This highly infectious organism has a worldwide distribution. The primary mycoplasmal infection often becomes complicated by secondary bacterial and viral infections [2], resulting in more severe lung lesions and production losses. Relative control has been achieved through active vaccination programs, but porcine enzootic pneumonia continues to be a major economic problem in the swine industry. While progress has been made in understanding the molecular basis of some Mycoplasma diseases [3], advances in M. hyopneumoniae research have been hampered by its fastidious growth condition and the lack of genetic tools and transformation protocols. To date, few virulence determinants or virulence-associated determinants have been identified. Attachment to the respiratory epithelium is a prerequisite for host colonization and is mediated by the membrane protein P97 [4]. This protein is located on the outer membrane surface, and its role in adherence has been firmly established. The general region of P97 that mediates adherence to swine cilia is thought to be the R1 region, near the C-terminus of the protein [5]. To bind cilia, a minimum of eight tandem copies of the pentapeptide sequence (AAKPV/E) in R1 are required [5]. Although the function of R2 in vivo is unknown, both it and R1 are required to bind heparin [6]. The P97 genes of M. hyopneumoniae strains 7448, 232, and J code for proteins with 10, 15, and 9 of the previously described R1 repeating units (AAKPV/E), respectively; all three strains had more than the minimum number of tandem copies (8 tandem copies) required for cilium binding [7]. Moreover, monoclonal antibodies F1B6 and F2G5, which both react predominantly with P97 [4,5], only partially block adherence of M. hyopneumoniae to receptors on epithelial cell cilia [8]. These observations indicate that molecules other than P97 play a role in facilitating adherence of M. hyopneumoniae to swine cilia. Comparative transcriptomic and proteomic studies are also performed to study transcriptional changes that occur during disease and investigate differentially expressed proteins in pathogenic and non-pathogenic strains [9-11]. Several M. hyopneumoniae proteins, including immunodominant proteins (P36 [12], P46 [13], and P65 [14]), adhesin-related proteins (P102 [15], P146 [16], P159 [17], P216 [18], and LppT [16]), and a 54-kDa cytotoxic factor [19], have been characterized; however, the biological functions of these proteins in pathogenesis are not well understood.

Comparative genomic analysis has previously revealed mechanisms of M. hyopneumoniae pathogenicity [7] and predicted unidentified virulence factors, including genes involved in secretion and/or traffic between host and pathogen cells, or with evasion and/or modulation of the host immune system [20,21]. In 2005, Vasconcelos et al. sequenced a pathogenic and a non-pathogenic strain of M. hyopneumoniae and performed a comparative genomics approach to identify putative virulence genes [7]. They identified various CDSs that could be considered candidate virulence genes, including cilium adhesin homologs, lipoproteins, and other components which might contribute to virulence [7]. However, comparative genomic analysis of a virulent M. hyopneumoniae strain versus its attenuated strain is lacking.

The need to control the spread of M. hyopneumoniae prompted the development of live attenuated vaccine strains. M. hyopneumoniae strain 168-L has been extensively used as vaccine against M. hyopneumoniae in China [22,23]. This attenuated vaccine strain is derived from the virulent parent strain 168. Strain 168 was originally isolated in 1974, from an Er-hua-nian pig (a Chinese local breed very sensitive to M. hyopneumoniae) with typical clinical and pathogenic characteristics of mycoplasmal pneumonia of swine (MPS) [24]. This field strain was gradually attenuated by more than 300 continuous passages through KM2 cell-free medium (a modified Friis medium) and the 380th passage was named strain 168-L. Currently, the genetic basis for the attenuation of virulence in 168-L is poorly understood.

To gain new insight into the components that contribute to virulence and the mechanisms by which M. hyopneumoniae causes disease, we sequenced the genomes of strains 168 and 168-L. This allowed us to perform the first comprehensive analysis of virulent and attenuated strains, and identify CDSs that may be related to virulence. We further investigated these putative virulence related CDSs and compared them with reported virulence determinants. Notably, almost all reported virulence determinants were found in putative virulence related CDSs. Besides the reported virulence determinants, other candidate virulence genes were also identified. The study of these candidate virulence genes and their corresponding products will be important to better comprehend the pathogenesis of M. hyopneumoniae.

Results and discussion

Genomic features of M. hyopneumoniae 168-L and its global comparison with pathogenic strain 168

The complete genome of M. hyopneumoniae 168-L consists of a 921,093 bp (GC content 28.46%) single circular chromosome (GenBank accession number CP003131). A total of 689 protein-encoding genes were predicted. The average protein size is 378 amino acids and the mean coding percentage is 84.8%. Approximately 51% of genes were assigned to specific functional clusters of orthologous groups (COGs), and 28% were assigned an enzyme classification (EC) number (Figure 1). Comparison with the M. hyopneumoniae 168 genome (GenBank accession CP002274) revealed a highly conserved gene content and order between the two strains. The 168-L genome is 4,483 bp smaller than that of 168 (925,576 bp), because there are 60 insertions and 43 deletions (indels; insertions and deletions of any size) in 168-L relative to 168 (see Additional file 1: Table S1; Additional file 2: Table S2; Additional file 3: Table S3). Among these, 33 indels are located in predicted CDSs, and 70 are in noncoding regions. Besides these indels, 227 single nucleotide variations (SNVs) were identified between 168 and 168-L (Additional file 4: Table S4). While 31 SNVs were mapped to intergenic regions, 196 were in coding regions, inducing amino acid substitutions, frame shifts, and translational stops.

Additional file 1: Table S1. Nucleotide differences identified between genomes of 168 and 168-L.

Format: XLS Size: 2.1MB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 2:Table S2. Insertions detected in 168-L compared to 168.

Format: XLS Size: 56KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 3: Table S3. Deletions detected in 168-L compared to 168.

Format: XLS Size: 52KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 4: Table S4. SNVs detected in 168-L compared to 168.

Format: XLS Size: 86KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

thumbnailFigure 1. Genome architecture. The dnaA gene is at position zero. Moving inside, the first circle shows the genome length (units in M. hyopneumoniae); the second and the third circles show the locations of the predicted CDSs on the plus and minus strands, respectively, which were color-coded by COG categories (the color codes for the functional assignments are shown in the key); the fourth circle shows tRNAs (purple) and rRNAs (red); the fifth circle shows the centered GC (G+C) content of each CDS (blue: above mean and cyan: below mean); and the sixth circle shows the GC (G+C) skew plot (red: above zero and pink: below zero). Circles 7–10 show comparative amino acids analysis of 168 with amino acids identities color-coded according to the similarity shown in the key to strains 168-L (seventh circle), 232 (eight circle), J (ninth circle), 7448 (tenth circle).

ISMHp1-Related genetic variations between 168 and 168-L

The difference between the genome sizes of strains 168 and 168-L is mainly due to differences in the duplication of Insertion Sequence (IS) elements. IS elements are distributed stochastically across the entire genome of both strains. The 168-L genome contains nine complete and one disrupted IS elements, which is almost identical to that of 168 except for slight differences in ISMHp1. There are nine complete copies of ISMHp1 in 168-L, but 12 copies in 168. The difference in ISMHp1 copy number between 168 and 168-L is due to three complete ISMHp1 deletions (located 690 kb, 870 kb, and 900 kb from oriC) and one complete ISMHp1 inversion, which was originally located at 378 kb from oriC, but was inverted in 168-L (1656 bp, located 372 kb from oriC) (Figure 2a).

thumbnailFigure 2. Alignment of the whole genomes and inverted regions. (a) Alignment of the M. hyopneumoniae 168 and 168-L genomes. (b) Alignment of the inverted regions of 168 relative to 168-L. The gray bars represent the forward and reverse strands. Green triangles represent ISMHP1 elements. DNA BLASTN alignments (BLASTN matches) between the two sequences are indicated by a red (same strand) or blue (opposite strand) line.

Other than the IS elements, notable large-scale genomic differences were also indicated. Compared to strain 168, a genomic deletion of approximately 1.36 kb (locus 1: between MHP168L_311 and MHP168L_729) was identified, which had been substituted with an approximately 2.32 kb novel insertion sequence (locus 2) that was joined to a complete ISMHp1 element in 168-L (Figure 2b). This 168-L-related insertion fragment was also observed in strains 7448 and 232.

Molecular analysis of integrative conjugative element (ICE)

The integrative conjugative element (ICE) is a mobile DNA that is probably involved in genomic recombination events and in pathogenicity. The ICEH elements are more divergent than the typical similarity of other chromosomal locus in M. hyopneunomiae[25], suggesting an accelerated evolution of these constins [26]. During a survey of specific sequences, a specific 26.9-kb region with similarity to the integrative conjugal element of M. fermentans (ICEF) [27] was found in strain 168, which was designated ICEH (for integrative conjugal element of M. hyopneumoniae). Unlike ICEH in strains 7448 and 232, which consist of nineteen and twenty two CDSs, respectively, the ICEH168 consist of 20 CDSs (Additional file 5: Table S5). The organization of these elements is very similar. Some CDSs present similarity to tra genes, which are usually associated with the bacteria conjugative plasmids such as traK, traI, traE[26]. The ICEH168 has three tra genes, with one traG and two copies of the traE gene. Besides, a CDS encoding for a single strand binding protein (SSB) that is essential for the transfer process is also observed.

Additional file 5: Table S5. Annotated CDS of the elements ICEH168 and ICEH168-L.

Format: XLS Size: 51KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The ICE analysis of three M. hyopneumoniae genomes (7448, J and 232) carried out previously, revealed that the ICEH is present in the two pathogenic strains (7448 and 232) but is absent from the non-pathogenic one (J strain) [26]. Interest has therefore shifted to questions of whether the ICEH is present in the attenuated vaccine strain 168-L. Interestingly, the ICEH was also observed in strain 168-L. Moreover, the ICEH168 and ICEH168-L are almost the same, except for a missense mutation (G192E) identified in ICEH-ORF3 (MHP168_235). Our analyses indicate that the ICEH may not only present in pathogenic strains of M. hyopneumoniae.

Mutations affecting epithelium adhesion

In our previous study, the ability of adherence and damage to the cilia between strains 168 and 168-L were compared by using scanning electron microscopy. The results showed that the pathogenic strain 168 adheres to cilia inducing tangling, clumping, and longitudinal splitting of cilia, while the strain 168-L does not cause ciliary damage comparing to control group [28]. The adherence of M. hyopneumoniae to porcine ciliated respiratory cells is essential for the organism to colonize the respiratory epithelium and cause pneumonia [4]. The adherence process is mainly mediated by receptor-ligand interactions, and the M. hyopneumoniae proteins possibly involved in these interactions are obvious candidates as virulence factors [8,29-31]. We investigated the genetic variation between strains 168 and 168-L (Table 1; Additional file 6: Table S6), and compared mutations affecting CDSs corresponding to previously described mycoplasma adhesins (P97, P102, P146, P159, P216, MgPa, LppS, and LppT) [7]. Notably, almost all the reported mycoplasma adhesins are included in the CDSs affected by mutations (Table 1).

Additional file 6: Table S6. Complete list of "168-L-specific" genetic variations in intergenic.

Format: XLS Size: 63KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 1. Complete list of "168-L-specific" genetic variations in CDSs

In 168-L, three transversions were identified in the R1 region, near the C-terminus, of P97 (MHP168_110/MHP168L_110), which encodes cilium adhesin. In M. hyopneumoniae, attachment to the respiratory epithelium is mainly mediated by the membrane protein P97 [4]. This protein is located on the outer membrane surface, and its role in adherence has been firmly established. To bind cilia, a minimum of eight tandem copies of the pentapeptide sequence (AAKPV/E) in R1 are required [5]. Notably, all three transversion mutations were located in the tandem repeat unit (AAKPV/E), causing an E863V substitution. Significant alteration in this critical repeat unit might partly affect the adhesion reaction in 168-L.

Previous studies have demonstrated that P102 binds fibronectin and contributes to the recruitment of plasmin(ogen) to the M. hyopneumoniae cell surface [15]. P102 is commonly linked to P97 cilium adhesin, forming a two-gene operon [32]. Both P97 and P102 have several paralogs within the M. hyopneumoniae genome. However, the paralogs have only part of the complete sequence. Interestingly, P102, the companion gene in this operon, was truncated at 564 bp by a single base insertion in strain 168. Another intact copy of P102 (99% identity) was found 85 kb from this operon. Conversely, in strain 168-L, the original truncated P102 was reverted to the intact one, while another intact copy of P102, 85 kb away, was truncated.

The P146 adhesin-like protein of M. hyopneumoniae shows strong similarity to the LppS lipoprotein of Mycoplasma conjunctivae, which is involved in in vitro adhesion [16]. In addition, the N-terminus region of P146 also shows strong similarity to the P97 adhesin, and has a strongly hydrophobic region (amino acids 7–29), indicating a transmembrane region, and suggesting that the protein is expressed on the surface of M. hyopneumoniae cells [33]. Compared with its counterpart MHP168_676 in 168, MHP168L_676 from 168-L has an in-frame insertion of one amino acid (Q) at the N-terminus of P146. The enormous intra-specific diversity shown for the P146 encoding gene is at least partly because of differences between several repeat regions present in the gene, most notably a polyserine chain of variable length, and a [Q]n[(P/S)Q]m repeat region [34]. Interestingly, this one in-frame insertion (Q) was located in the [Q]n[(P/S)Q]m repeat region. Polyserine chains often function as a spacer region in proteins involved in complex carbohydrate degradation [35], while sequences rich in both proline and glutamine are not uncommon and can form a conformation known as a polyproline II helix [36,37]. Such proline-rich sequences are often involved in binding processes and are highly immunogenic [37]. However, because the function of the P146 protein remains unknown, correlations with virulence or adhesion are speculative and need further investigation.

P159 is a proteolytically processed surface adhesin of M. hyopneumoniae[17]. Three proteins with apparent molecular masses of 27 (P27), 52 (P52), and 110 (P110) kDa were identified through proteomic analysis of M. hyopneumoniae lysates [17], with each representing a different region spanning P159. These cleavage fragments are located on the cell surface and present at all growth stages. In 168-L, MHP168L_504 (P159) has a missense mutation resulting in a G240A replacement in the (S)(S)G(G)S repeat region of P159. Although this (S)(S)G(G)S repeat region has been reported, its biological function is unknown.

P216 (MHP168L_503/MHP168_503) is a proteolytically processed cilium and heparin binding protein of M. hyopneumoniae[18]. This surface protein is post-translationally processed to generate N-terminal P120 and C-terminal P85 fragments, both of which can bind cilia [18]. The 168-L P216 gene has an in-frame four amino acid deletion in a poly Q motif near the C-terminus. Previous studies have suggested that poly Q and KEKE motifs may play a role in maintaining P85 on the cell surface [18,34]. Collectively, the deletion mutations affecting P216 may affect its cilium adhesion and may be associated with virulence attenuation in 168-L.

The MHP168L_424 gene and its gene product LppT were analyzed in detail because they showed approximately 22% identity to the LppT protein from M. conjunctivae. LppT is the second gene in a two gene operon with LppS, which was reported to be an adhesin in M. conjunctivae[16]. The LppT gene lacked a promoter and is likely to be co-transcribed with LppS, thus suggesting a functional relationship between LppS and LppT [16]. In M. hyopneumoniae, LppT encoded a protein of 954 aa with a calculated molecular mass 108 kDa. The gene product encoded by LppT is also a membrane protein, with a signal sequence of 34 aa at the amino-terminal end and a transmembrane structure. Notably, one of the amino acid substitutions in 168-L occurs near the C-terminus of LppT, resulting in a L814F replacement.

Mutations altering the cell envelope and genes encoding secreted proteins

Cell envelope proteins and secreted proteins are involved in virulence, host cell interaction, and immune responses [14,38,39]. Outer membrane protein-P95 is a cell envelope protein in M. hyopneumoniae. In 168-L, MHP168L_103 (P95) is truncated by a nonsense mutation compared to 168, resulting in an E965* termination near the C-terminus of P95. Significant alteration in this outer membrane protein could conceivably cause a truncation in coding region, and in turn alter the function of P95 outer membrane protein.

M. hyopneumoniae contains an abbreviated membrane protein secretory system [1]. The pathway consists of secA (MHP168L_088), secY (MHP168L_128), secD (MHP168L_259), prsA (MHP168L_664), dnaK (MHP168L_069), trigger factor (MHP168L_154), and lepA (MHP168L_076). It has recently been demonstrated that some pathogenic bacteria use a type IV secretion system, composed of subunits related to the conjugation machinery, to deliver effector molecules to host cells [40], and that this system may be involved in pathogenesis [41]. We found no pathogenic mutations in the protein secretory system, except a synonymous substitution (P407P) in MHP168L_069, which encodes chaperone protein DnaK.

Mutations affecting antigens

Studies on the antigenic properties of M. hyopneumoniae revealed several immunodominant proteins, including the P36 cytosolic protein [12,42], the P46, P65, and P74 membranous proteins [43-46], the elongation factor Tu [47], the chaperone protein DnaK[47], the pyruvate dehydrogenase E1-beta subunit [47], and the P97 adhesin [4]. The functions of these proteins have not been well elucidated, but specific reactants may eventually be useful tools to diagnose M. hyopneumoniae[48].

The cytosolic P36 protein is a lactate dehydrogenase [49] that induced an early immune response in pigs that are experimentally and naturally infected by M. hyopneumoniae[50]. Comparative studies with other Mycoplasmas commonly found in pigs demonstrated that the P36 proteins carry highly conserved species-specific antigenic determinants for M. hyopneumoniae[42]. Hyperimmune sera produced against recombinant P36 protein showed no reactivity against other porcine Mycoplasmas, including M. flocculare, M. hyorhinis, and Acholeplasma laidlawii[12]. Notably, one of the observed amino acid substitutions in 168-L occurs near the C-terminus of P36 (MHP168L_167), resulting in a N204D replacement.

P65 is an immunodominant surface lipoprotein of M. hyopneumoniae that is specifically recognized during disease [14]. Analysis of the translated amino acid sequence of the gene encoding p65 revealed similarity to the GDSL family of lipolytic enzymes [14]. The monospecific antibodies against heat shock protein-like P42 antigen, part of P65, can block the growth of M. hyopneumoniae[51]. In 168-L, MHP168L_668 (P65) has a missense mutation resulting in a T138A replacement.

Mutations affecting transport proteins

As Mycoplasmas are dependent on the exogenous supply of many nutrients, it has been predicted that they may need many transport systems [3]. Motif analysis revealed a family of proteins with a phosphotransferase (PTS) motif. The open reading frames (ORFs) included sgaA (MHP168L_422), sgaB (MHP168L_421), sgaT (MHP168L_563), mtlF (MHP168L_739), mtlA (MHP168L_561), nagE (MHP168L_582), and licA (MHP168L_041) [7]. However, no mutations were identified in this PTS transporter family. There are approximately 30 genes with ABC transporter family signatures in the genome of M. hyopneumoniae 168-L (Additional file 7: Table S7), and five missense mutations and one synonymous substitution were identified in this group. These included a cobalt import ATP-binding protein (MHP168L_284, K64E), an ABC transporter permease protein (MHP168L_394, A492G), a xylose ABC transporter ATP-binding protein (MHP168L_523, S8S), and three ABC transporter ATP-binding proteins (MHP168L_413, D118Y; MHP168L_462, V219L; MHP168L_631, R60W). Interestingly, the expression of MHP168L_394 and MHP168L_413 was reported to be up-regulated in vivo during disease relative to in vitro-grown [11]. The variability between strains 168 and 168-L in multi-transport proteins indicates that they may affect growth and survival in different hosts or host tissues.

Additional file 7: Table S7. Mutations affecting genes with ABC transporter family signatures.

Format: XLS Size: 47KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Mutations affecting genes directly related to metabolism and in vivo growth

M. hyopneumoniae strain 168 encodes 695 genes, approximately one quarter of which are involved in metabolism and in vivo growth. Of particular interest were the mutations observed in the genes involved in various metabolic pathways (Figure 3), including glycolysis/gluconeogenesis (MHP168_167, MHP168_186), purine metabolism (MHP168_289, MHP168_639), pyrimidine metabolism (MHP168_086), glycerophospholipid metabolism (MHP168_596), oxidative phosphorylation (MHP168_085), aminoacyl-tRNA biosynthesis (MHP168_058), and the pentose phosphate pathway (MHP168_142, MHP168_152). An in-frame insertion of two amino acids (TG), a missense mutation (E393G) and a nonsense mutation (N456*) were identified in 168-L near the C-terminus of MHP168L_085, which encodes a NADH oxidase involved in oxidative phosphorylation. Mycoplasma genomes are deficient in genes coding for components of intermediary and energy metabolism [3]. Thus, Mycoplasmas depend mostly on glycolysis to synthesize ATP [3]. In the glycolysis pathway, missense mutations in both L-lactate and pyruvate dehydrogenase were observed, resulting in N204D and S194G replacements, respectively. Iron deprivation, is a prominent feature of the host innate immune response, and most certainly impacts growth of Mycoplasmas in vivo[52]. Through transcriptome analysis, MHP168_639 was identified to be down-regulated during iron limiting conditions [52]. This suggests that MHP168_639 may play a role in M. hyopneumoniae’s response to iron stress. In 168-L, MHP168L_639 has a missense mutation resulting in an E411K replacement. Mutations in these metabolism-related genes accumulated over 300 in vitro passages likely affect growth and survival within host cells.

thumbnailFigure 3. Metabolic potential. The metabolic pathways of M. hyopneumoniae strains 168 and 168-L were mapped and analyzed using KEGG Pathway Database. Those pathways, containing mutations affected metabolic-related genes, are shown in green.

In addition, approximately 41% of mutations affected genes coding for hypothetical proteins. Despite the lack of functional annotations for these genes, their disruption in 168-L makes them obvious targets for investigation as potential virulence factors. Further molecular genetics and in vivo studies are required to confirm and assess the relative importance of these genes in the attenuation of virulence in 168-L.

Conclusions

We successfully used a combination of sequencing genomics and comparative genomics strategies to provide a comprehensive analysis of virulent and attenuated M. hyopneumoniae strains to identify determinants involved in pathogenesis. The genome of the attenuated high-passage derivative strain 168-L was sequenced and compared to virulent strain 168, revealing mutations in numerous CDSs. These mutations affected CDSs are likely to be associated with virulence. We then compared these putative virulence factor CDSs to reported virulence determinants. Notably, almost all of the reported M. hyopneumoniae virulence determinants were included in the list of putative virulence factor CDSs. Variations in the previously described mycoplasma adhesins (P97, P102, P146, P159, P216 and LppT), cell envelope proteins (P95), cell surface antigens (P36), secreted proteins, chaperone protein (DnaK), and genes directly related to metabolism and in vivo growth may contribute to loss of virulence in 168-L. We then proceeded to characterize the alterations in gene functions caused by mutations at the protein level, and compared those mutations with previously described repeat motifs that may be of primary importance for virulence [34]. Interestingly, we found that many mutations were located in the virulence associated motifs of the various proteins. To bind cilia, a minimum of eight tandem copies of the pentapeptide sequence (AAKPV/E) in the R1 region of P97 are required [5]. We identified three mutations in the tandem repeat unit (AAKPV/E), causing an E863V substitution. A similar situation was also observed in several other virulence associated genes (P146, P159, P216, and LppT). We hypothesize that the cumulative effect of mutations in virulence associated genes may account for the attenuation of virulence in 168-L. In this study, a total of 330 genetic variations were identified. While these included reported virulence-related genes, other novel virulence determinants were also identified. However, further molecular genetics and in vivo studies are required to confirm and assess the relative importance of these suspected novel virulence determinants in the attenuation of virulence. The comparative genomic analysis presented here will not only provide insights into the basis of attenuation of virulence in 168-L, but may also provide targets for mutagenesis in the pursuit of development of a more efficacious vaccine.

Methods

Bacterial strains, growth conditions, and DNA extraction

Clonal isolates of M. hyopneumoniae strain 168 and 168-L were selected for sequencing. Both of the strains were grown in KM2 cell-free medium at 37°C. The culture was harvested from 100 mL KM2 cell-free medium by centrifugation at 1,200×g for 30 min, and then total genomic DNA was extracted from mycoplasma cultures using a TIANamp Bacteria DNA Kit (Tiangen, Beijing, China) according to the manufacturer’s instructions.

Genome sequencing and assembly

Genomic libraries containing 8 kb inserts were constructed according to the manufacturer’s protocols. Whole-genome sequencing of strain 168-L was performed by combining GS FLX and Solexa paired-end sequencing technologies. A total of 242,507 reads (67.4% paired ends) were produced with the GS FLX system, giving 44.5-fold coverage of the genome. Eighty-eight percent (215,346) of reads were assembled into one large scaffold using Newbler (454 Life Sciences, Branford, CT, USA). A total of 1,971,358 reads were generated with an Illumina Solexa genome analyzer IIx (Illumina, San Diego, CA, USA) and were mapped to the scaffold with the Burrows-Wheeler Alignment (BWA) tool [53]. Gaps were filled by local assembly of the Solexa/Roche 454 reads or sequencing PCR products using a Prism 3730 capillary sequencer (Applied Biosystems, Foster City, CA, USA). All repeated DNA regions and low-quality regions were verified by PCR and sequencing of the product amplified from genomic DNA.

Annotation and sequence analyses

Open reading frames containing more than 30 amino acid residues were predicted using Glimmer 3.0 [54] with modified genetic code 4 and verified manually using the strain 168 annotation. Loci discrepancies between the 168 and 168-L consensus sequences were manually examined for support at the trace data level. Transfer RNA (tRNA) and ribosomal RNA (rRNA) genes were predicted using the tRNAscan-SE program [55] or by observing similarities with the M. hyopneumoniae strain 232 and strain J rRNA genes. Artemis (release 12) [56] was used to collate and annotate data. Functional predictions were based on BLASTP similarity searches against the UniProtKB [57], GenBank [58], Swiss-Prot protein [59], and COG [60] databases. EC numbers were assigned using the Kyoto Encyclopedia of Genes and Genomes (KEGG) [67] and metabolic pathways were mapped and analyzed using KEGG Pathway Database (http://www.genome.jp/kegg/pathway.html webcite). Pseudogenes were detected by BLASTN analysis, comparing the genome sequences of 168-L with those of 232 and J, and then the annotation was revised manually.

Single nucleotide polymorphism (SNP) analysis

Nucleotide comparisons and single nucleotide polymorphism (SNP) analysis for strains 168 and 168-L were performed using the Artemis Comparison Tool (ACT) [61] and Mauve 2.3.1 genome alignment software [62]. ORF graphical visualization and manual annotation were carried out using Artemis, release 12 [56]. Screening for unusual coding differences between the 168 and 168-L genomes (stops and frame shifts) was conducted using FASTA program packages [63,64] and BLAST [65]. The coding differences between the 168 and 168-L genomes were checked manually.

Accession numbers

M. hyopneumoniae strains 168 and 168-L genome sequences have been deposited in GenBank under accession numbers CP002274.1 and CP003131, respectively.

Abbreviations

EP: Porcine Enzootic Pneumonia; SNVs: Single Nucleotide Variations; MPS: Mycoplasmal Pneumonia Of Swine; COG: Clusters Of Orthologous Groups; EC: Enzyme Classification; IS: Insertion Sequence; ICE: Integrative Conjugative Element; SSB: Single Strand Binding Protein; PTS: Phosphotransferase; BWA: Burrows-Wheeler Alignment; tRNA: Transfer RNA; rRNA: ribosomal RNA; KEGG: Kyoto Encyclopedia of Genes and Genomes; ACT: Artemis Comparison Tool.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: WL SX LF HC. Performed the experiments: WL ML SL SG ZF. Analyzed the data: WL ZZ GS RL BL. Contributed reagents/materials/analysis tools: WL GS RL BL. Wrote the paper: WL LF SX. All authors read and approved the final manuscript.

Acknowledgements

Dr. Lei Wang and members of his lab at Nankai University are thanked for their assistance with sequencing and analysis. We would also like to thank Dr. Maojun Liu at Jiangsu Academy of Agricultural Sciences for aid in Preparation of Materials. This study was supported by the National Natural Sciences Foundation of China (31170160, 31001080, and 31001066).

References

  1. Minion FC, Lefkowitz EJ, Madsen ML, Cleary BJ, Swartzell SM, Mahairas GG: The genome sequence of Mycoplasma hyopneumoniae strain 232, the agent of swine mycoplasmosis.

    J Bacteriol 2004, 186:7123-7133. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Maes D, Segales J, Meyns T, Sibila M, Pieters M, Haesebrouck F: Control of Mycoplasma hyopneumoniae infections in pigs.

    Vet Microbiol 2008, 126:297-309. PubMed Abstract | Publisher Full Text OpenURL

  3. Razin S, Yogev D, Naot Y: Molecular biology and pathogenicity of mycoplasmas.

    Microbiol Mol Biol Rev: MMBR 1998, 62:1094-1156. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Zhang Q, Young TF, Ross RF: Identification and characterization of a Mycoplasma hyopneumoniae adhesin.

    Infect Immun 1995, 63:1013-1019. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Minion FC, Adams C, Hsu T: R1 region of P97 mediates adherence of Mycoplasma hyopneumoniae to swine cilia.

    Infect Immun 2000, 68:3056-3060. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Jenkins C, Wilton JL, Minion FC, Falconer L, Walker MJ, Djordjevic SP: Two domains within the Mycoplasma hyopneumoniae cilium adhesin bind heparin.

    Infect Immun 2006, 74:481-487. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Vasconcelos AT, Ferreira HB, Bizarro CV, Bonatto SL, Carvalho MO, Pinto PM, Almeida DF, Almeida LG, Almeida R, Alves-Filho L, et al.: Swine and poultry pathogens: the complete genome sequences of two strains of Mycoplasma hyopneumoniae and a strain of Mycoplasma synoviae.

    J Bacteriol 2005, 187:5568-5577. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Zhang Q, Young TF, Ross RF: Microtiter plate adherence assay and receptor analogs for Mycoplasma hyopneumoniae.

    Infect Immun 1994, 62:1616-1622. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Li YZ, Ho YP, Chen ST, Chiou TW, Li ZS, Shiuan D: Proteomic comparative analysis of pathogenic strain 232 and avirulent strain J of Mycoplasma hyopneumoniae.

    Biochemistry Biokhimiia 2009, 74:215-220. PubMed Abstract | Publisher Full Text OpenURL

  10. Pinto PM, Klein CS, Zaha A, Ferreira HB: Comparative proteomic analysis of pathogenic and non-pathogenic strains from the swine pathogen Mycoplasma hyopneumoniae.

    Proteome sci 2009, 7:45. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  11. Madsen ML, Puttamreddy S, Thacker EL, Carruthers MD, Minion FC: Transcriptome changes in Mycoplasma hyopneumoniae during infection.

    Infect Immun 2008, 76:658-663. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Strasser M, Frey J, Bestetti G, Kobisch M, Nicolet J: Cloning and expression of a species-specific early immunogenic 36-kilodalton protein of Mycoplasma hyopneumoniae in Escherichia coli.

    Infect Immun 1991, 59:1217-1222. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Futo S, Seto Y, Okada M, Sato S, Suzuki T, Kawai K, Imada Y, Mori Y: Recombinant 46-kilodalton surface antigen (P46) of Mycoplasma hyopneumoniae expressed in Escherichia coli can be used for early specific diagnosis of mycoplasmal pneumonia of swine by enzyme-linked immunosorbent assay.

    J Clin Microbiol 1995, 33:680-683. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Schmidt JA, Browning GF, Markham PF: Mycoplasma hyopneumoniae p65 surface lipoprotein is a lipolytic enzyme with a preference for shorter-chain fatty acids.

    J Bacteriol 2004, 186:5790-5798. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Seymour LM, Jenkins C, Deutscher AT, Raymond BB, Padula MP, Tacchi JL, Bogema DR, Eamens GJ, Woolley LK, Dixon NE, et al.: Mhp182 (P102) binds fibronectin and contributes to the recruitment of plasmin(ogen) to the Mycoplasma hyopneumoniae cell surface.

    Cell Microbiol 2012, 14:81-94. PubMed Abstract | Publisher Full Text OpenURL

  16. Belloy L, Vilei EM, Giacometti M, Frey J: Characterization of LppS, an adhesin of Mycoplasma conjunctivae.

    Microbiology 2003, 149:185-193. PubMed Abstract | Publisher Full Text OpenURL

  17. Burnett TA, Dinkla K, Rohde M, Chhatwal GS, Uphoff C, Srivastava M, Cordwell SJ, Geary S, Liao X, Minion FC, et al.: P159 is a proteolytically processed, surface adhesin of Mycoplasma hyopneumoniae: defined domains of P159 bind heparin and promote adherence to eukaryote cells.

    Mol Microbiol 2006, 60:669-686. PubMed Abstract | Publisher Full Text OpenURL

  18. Wilton J, Jenkins C, Cordwell SJ, Falconer L, Minion FC, Oneal DC, Djordjevic MA, Connolly A, Barchia I, Walker MJ, et al.: Mhp493 (P216) is a proteolytically processed, cilium and heparin binding protein of Mycoplasma hyopneumoniae.

    Mol Microbiol 2009, 71:566-582. PubMed Abstract | Publisher Full Text OpenURL

  19. Geary SJ, Walczak EM: Isolation of a cytopathic factor from Mycoplasma hyopneumoniae.

    Infect Immun 1985, 48:576-578. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Szczepanek SM, Tulman ER, Gorton TS, Liao X, Lu Z, Zinski J, Aziz F, Frasca S Jr, Kutish GF, Geary SJ: Comparative genomic analyses of attenuated strains of Mycoplasma gallisepticum.

    Infect Immun 2010, 78:1760-1771. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Zheng H, Lu L, Wang B, Pu S, Zhang X, Zhu G, Shi W, Zhang L, Wang H, Wang S, et al.: Genetic basis of virulence attenuation revealed by comparative genomic analysis of Mycobacterium tuberculosis strain H37Ra versus H37Rv.

    PLoS One 2008, 3:e2375. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Feng Z, Shao G, Liu M, Wu X, Zhou Y, Gan Y: Immune responses to the attenuated Mycoplasma hyopneumoniae 168 strain vaccine by intrapulmonic immunization in piglets.

    Agric Sci China 2010, 9:423-431. Publisher Full Text OpenURL

  23. Li Y, Li P, Wang X, Yu Q, Yang Q: Co-administration of attenuated Mycoplasma hyopneumoniae 168 strain with bacterial DNA enhances the local and systemic immune response after intranasal vaccination in pigs.

    Vaccine 2012, 30:2153-2158. PubMed Abstract | Publisher Full Text OpenURL

  24. Liu W, Feng Z, Fang L, Zhou Z, Li Q, Li S, Luo R, Wang L, Chen H, Shao G, et al.: Complete genome sequence of Mycoplasma hyopneumoniae strain 168.

    J Bacteriol 2011, 193:1016-1017. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Madsen ML, Oneal MJ, Gardner SW, Strait EL, Nettleton D, Thacker EL, Minion FC: Array-based genomic comparative hybridization analysis of field strains of Mycoplasma hyopneumoniae.

    J Bacteriol 2007, 189:7977-7982. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Pinto PM, de Carvalho MO, Alves-Junior L, Brocchi M, Schrank IS: Molecular analysis of an integrative conjugative element, ICEH, present in the chromosome of different strains of Mycoplasma hyopneumoniae.

    Genet Mol Biol 2007, 30:256-263. Publisher Full Text OpenURL

  27. Calcutt MJ, Lewis MS, Wise KS: Molecular genetic analysis of ICEF, an integrative conjugal element that is present as a repetitive sequence in the chromosome of Mycoplasma fermentans PG18.

    J Bacteriol 2002, 184:6929-6941. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Shao G, Liu M, Sun P, Wang J, Du G, Zhou Y, Liu D: The establishment of an experimental swine model of swine mycoplasma pneumonia.

    J Microbes Infec 2007, 2:215-218. OpenURL

  29. Zhang Q, Young TF, Ross RF: Glycolipid receptors for attachment of Mycoplasma hyopneumoniae to porcine respiratory ciliated cells.

    Infect Immun 1994, 62:4367-4373. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Zielinski GC, Ross RF: Adherence of Mycoplasma hyopneumoniae to porcine ciliated respiratory tract cells.

    Amj Vet Res 1993, 54:1262-1269. OpenURL

  31. Zielinski GC, Young T, Ross RF, Rosenbusch RF: Adherence of Mycoplasma hyopneumoniae to cell monolayers.

    Amj Vet Res 1990, 51:339-343. OpenURL

  32. Adams C, Pitzer J, Minion FC: In vivo expression analysis of the P97 and P102 paralog families of Mycoplasma hyopneumoniae.

    Infect Immun 2005, 73:7784-7787. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Stakenborg T, Vicca J, Maes D, Peeters J, de Kruif A, Haesebrouck F, Butaye P: Comparison of molecular techniques for the typing of Mycoplasma hyopneumoniae isolates.

    J Microbiol Meth 2006, 66:263-275. Publisher Full Text OpenURL

  34. de Castro LA, Rodrigues Pedroso T, Kuchiishi SS, Ramenzoni M, Kich JD, Zaha A, Henning Vainstein M, Bunselmeyer Ferreira H: Variable number of tandem aminoacid repeats in adhesion-related CDS products in Mycoplasma hyopneumoniae strains.

    Vet Microbiol 2006, 116:258-269. PubMed Abstract | Publisher Full Text OpenURL

  35. Howard MB, Ekborg NA, Taylor LE, Hutcheson SW, Weiner RM: Identification and analysis of polyserine linker domains in prokaryotic proteins with emphasis on the marine bacterium Microbulbifer degradans.

    Protein Sci 2004, 13:1422-1425. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Rucker AL, Creamer TP: Polyproline II helical structure in protein unfolded states: lysine peptides revisited.

    Protein Sci 2002, 11:980-985. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Kay BK, Williamson MP, Sudol M: The importance of being proline: the interaction of proline-rich motifs in signaling proteins with their cognate domains.

    FASEB J 2000, 14:231-241. PubMed Abstract | Publisher Full Text OpenURL

  38. Deutscher AT, Tacchi JL, Minion FC, Padula MP, Crossett B, Bogema DR, Jenkins C, Kuit TA, Walker MJ, Djordjevic SP: Mycoplasma hyopneumoniae surface proteins Mhp385 and Mhp384 bind host cilia and glycosaminoglycans and are endoproteolytically processed by proteases that recognize different cleavage motifs.

    J Proteome Res 2012, 11:1924-1936. PubMed Abstract | Publisher Full Text OpenURL

  39. Cheikh Saad Bouh K, Shareck F, Dea S: Monoclonal antibodies to Escherichia coli-expressed P46 and P65 membranous proteins for specific immunodetection of Mycoplasma hyopneumoniae in lungs of infected pigs.

    Clin Diagn Lab Immunol 2003, 10:459-468. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Ding Z, Atmakuri K, Christie PJ: The outs and ins of bacterial type IV secretion substrates.

    Trends Microbiol 2003, 11:527-535. PubMed Abstract | Publisher Full Text OpenURL

  41. Seubert A, Hiestand R, de la Cruz F, Dehio C: A bacterial conjugation machinery recruited for pathogenesis.

    Mol Microbiol 2003, 49:1253-1266. PubMed Abstract | Publisher Full Text OpenURL

  42. Stipkovits L, Nicolet J, Haldimann A, Frey J: Use of antibodies against the P36 protein of Mycoplasma hyopneumoniae for the identification of M. Hyopneumoniae strains.

    Mol Cell Probes 1991, 5:451-457. PubMed Abstract | Publisher Full Text OpenURL

  43. Kim MF, Heidari MB, Stull SJ, McIntosh MA, Wise KS: Identification and mapping of an immunogenic region of Mycoplasma hyopneumoniae p65 surface lipoprotein expressed in Escherichia coli from a cloned genomic fragment.

    Infect Immun 1990, 58:2637-2643. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Klinkert MQ, Herrmann R, Schaller H: Surface proteins of Mycoplasma hyopneumoniae identified from an Escherichia coli expression plasmid library.

    Infect Immun 1985, 49:329-335. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Mori Y, Hamaoka T, Sato S: Use of monoclonal antibody in an enzyme-linked immunosorbent assay (ELISA) for the detection of antibodies against Mycoplasma hyopneumoniae.

    Isr J Med Sci 1987, 23:657-662. PubMed Abstract OpenURL

  46. Brooks E, Faulds D: The Mycoplasma hyopneumoniae 74.5kD antigen elicits neutralizing antibodies and shares sequence similarity with heat-shock proteins.

    Vaccine 1989, 89:265-269. OpenURL

  47. Pinto PM, Chemale G, de Castro LA, Costa AP, Kich JD, Vainstein MH, Zaha A, Ferreira HB: Proteomic survey of the pathogenic Mycoplasma hyopneumoniae strain 7448 and identification of novel post-translationally modified and antigenic proteins.

    Vet Microbiol 2007, 121:83-93. PubMed Abstract | Publisher Full Text OpenURL

  48. Caron J, Sawyer N, Ben Abdel Moumen B, Cheikh Saad Bouh K, Dea S: Species-specific monoclonal antibodies to Escherichia coli-expressed p36 cytosolic protein of Mycoplasma hyopneumoniae.

    Clin Diagn Lab Immunol 2000, 7:528-535. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Haldimann A, Nicolet J, Frey J: DNA sequence determination and biochemical analysis of the immunogenic protein P36, the lactate dehydrogenase (LDH) of Mycoplasma hyopneumoniae.

    J Gen Microbiol 1993, 139:317-323. PubMed Abstract OpenURL

  50. Frey J, Haldimann A, Kobisch M, Nicolet J: Immune response against the L-lactate dehydrogenase of Mycoplasma hyopneumoniae in enzootic pneumonia of swine.

    Microb Pathog 1994, 17:313-322. PubMed Abstract | Publisher Full Text OpenURL

  51. Chou SY, Chung TL, Chen RJ, Ro LH, Tsui PI, Shiuan D: Molecular cloning and analysis of a HSP (heat shock protein)-like 42 kDa antigen gene of Mycoplasma hyopneumoniae.

    Biochem Mol Biol Int 1997, 41:821-831. PubMed Abstract OpenURL

  52. Madsen ML, Nettleton D, Thacker EL, Minion FC: Transcriptional profiling of Mycoplasma hyopneumoniae during iron depletion using microarrays.

    Microbiology 2006, 152:937-944. PubMed Abstract | Publisher Full Text OpenURL

  53. Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform.

    Bioinformatics 2009, 25:1754-1760. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Delcher AL, Bratke KA, Powers EC, Salzberg SL: Identifying bacterial genes and endosymbiont DNA with Glimmer.

    Bioinformatics 2007, 23:673-679. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence.

    Nucleic Acids Res 1997, 25:955-964. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation.

    Bioinformatics 2000, 16:944-945. PubMed Abstract | Publisher Full Text OpenURL

  57. Emmert DB, Stoehr PJ, Stoesser G, Cameron GN: The European Bioinformatics Institute (EBI) databases.

    Nucleic Acids Res 1994, 22:3445-3449. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Rapp BA, Wheeler DL: GenBank.

    Nucleic Acids Res 2000, 28:15-18. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Bairoch A, Apweiler R: The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1999.

    Nucleic Acids Res 1999, 27:49-54. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families.

    Science 1997, 278:631-637. PubMed Abstract | Publisher Full Text OpenURL

  61. Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis Comparison Tool.

    Bioinformatics 2005, 21:3422-3423. PubMed Abstract | Publisher Full Text OpenURL

  62. Darling AC, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements.

    Genome Res 2004, 14:1394-1403. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Pearson WR, Lipman DJ: Improved tools for biological sequence comparison.

    Proc Natl Acad Sci U S A 1988, 85:2444-2448. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Pearson WR, Wood T, Zhang Z, Miller W: Comparison of DNA sequences with protein sequences.

    Genomics 1997, 46:24-36. PubMed Abstract | Publisher Full Text OpenURL

  65. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215:403-410. PubMed Abstract OpenURL