Dataset Information

Comparative genomic analysis of Mycoplasma anatis strains.

ABSTRACT:

Background

The Gram-negative intracellular bacterium Mycoplasma anatis is a pathogen of respiratory infectious diseases in ducks and has caused significant economic losses in the poultry industry.

Objective

This study, as the first report of the structure and function of the pan-genome of Mycoplasma anatis, may provide a valuable genetic basis for many aspects of future research on the pathogens of waterfowl.

Methods

We sequenced the whole genomes of 15 Mycoplasma anatis isolated from ducks in China. Draft genome sequencing was carried out and whole-genome sequencing was performed by the sequencers of the PacBio Sequel and an IonTorrent Personal Genome Machine (PGM). Then the common genic elements of protein-coding genes, tRNAs, and rRNAs of Mycoplasma anatis genomes were predicted by using the pipeline Prokka v1.13.7. To investigate homologous protein clusters across Mycoplasma anatis genomes, we adopted Roary v3.13.0 to cluster orthologous genes (OGs) based on the following criteria.

Results

We obtained one complete genome and 14 genome sketches. Microbial mobile genetic element analysis revealed the distribution of insertion sequences (IS30, IS3, and IS1634), prophage regions, and CRISPR arrays in the genome of Mycoplasma anatis. Comparative genomic analysis decoded the genetic components and functional classification of the pan-genome of Mycoplasma anatis that comprised 646 core genes, 231 dispensable genes and among them 110 was strain-specific. Virulence-related gene profiles of Mycoplasma anatis were systematically identified, and the products of these genes included bacterial ABC transporter systems, iron transport proteins, toxins, and secretion systems.

Conclusion

A complete virulence-related gene profile of Mycoplasma anatis has been identified, most of the genes are highly conserved in all strains. Sequencing results are relevant to the molecular mechanisms of drug resistance, adaptive evolution of pathogens, population structure, and vaccine development.

SUBMITTER: Zhou Q

PROVIDER: S-EPMC8237044 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Mycoplasma hyopneumoniae is the causative agent of porcine enzootic pneumonia and a major factor in the porcine respiratory disease complex. A clear understanding of the mechanisms of pathogenesis does not exist although it is clear that M. hyopneumoniae adheres to porcine ciliated epithelium by action of a protein called P97. Previous studies have shown variation in the gene encoding the P97cilium adhesin within different strains of M. hyopneumoniae, but the extent of genetic variation among field strains across the genome is not known. Since M. hyopneumoniae is a worldwide problem, it is reasonable to expect that a wide range of genetic variability may exist given all of the different breed and housing conditions. This variation may impact the overall virulence of a single strain. Using microarray technology, this study examined potential variation of fourteen field strains in comparison to strain 232 on which the array was based. Genomic DNA was obtained, amplified with TempliPhiâ¢, and labeled indirectly with Alexa dyes. Post genomic hybridization, the arrays were scanned and data analyzed using a linear statistical model. Results indicate that genetic variation could be detected in all fourteen field strains but across different loci, suggesting that variation occurs throughout the genome. Fifty-nine percent of the variable loci were hypothetical genes. Twenty-two percent of the lipoprotein genes showed variation in at least one field strain. A permutation test identified a location in M. hyopneumoniae genome where spatial clustering of variability between the field strains and strain 232 exists. TempliPhiâ¢ amplified DNA samples from field isolates were compared to control strain 232 using a two-color experimental microarray design. Independent samples from one isolate labeled with one dye were paired with control samples labeled with the alternate dye; the samples were mixed and hybridized to the microarray. For nine of the fourteen isolates (95MP1501, 95MP1502, 95MP1503, 95MP1504, 95MP1508, 95MP1509, 97MP0001, 00MP1301, and 05MP2301), four independent field isolate DNA samples were paired with four independent DNA samples from control 232. In two of the four arrays, the control sample was labeled with Alexa 555 dye and compared to the field isolate sample labeled with Alexa 647 dye (Molecular Probes, Inc., Eugene, Ore.). The dye assignment to control and treated samples was reversed for the other two arrays (dye swap). The arrays were hybridized under identical conditions as described below. This procedure was repeated for isolates 95MP1510 for a total of four arrays, where the control sample was labeled with Alexa 647 dye in three of the arrays and Alexa 555 dye for the fourth array. For isolates 95MP1505, 95MP1506 and 95MP1507, a total of five arrays each, including two dye swaps, were done; and for isolate 00MP1502, a total of six arrays were done with the control labeled with Alexa 555 dye for four of the arrays and Alexa 647 for the other two arrays.