Re-sequencing and optical mapping reveals misassemblies and real inversions on Corynebacterium pseudotuberculosis genomes

Resequencing and characterization of the first Corynebacterium pseudotuberculosis genome isolated from camel

Background

Corynebacterium pseudotuberculosis is a zoonotic Gram-positive bacterial pathogen known to cause different diseases in many mammals, including lymph node abscesses in camels. Strains from biovars equi and ovis of C. pseudotuberculosis can infect camels. Comparative genomics could help to identify features related to host adaptation, and currently strain Cp162 from biovar equi is the only one from camel with a sequenced genome.

Methods

In this work, we compared the quality of three genome assemblies of strain Cp162 that used data from the DNA sequencing platforms SOLiD v3 Plus, IonTorrent PGM, and Illumina HiSeq 2500 with an optical map and investigate the unique features of this strain. For this purpose, we applied comparative genomic analysis on the different Cp162 genome assembly versions and included other 129 genomes from the same species.

Results

Since the first version of the genome, there was an increase of 88 Kbp and 121 protein-coding sequences, a decrease of pseudogenes from 139 to 53, and two inversions and one rearrangement corrected. We identified 30 virulence genes, none associated to the camel host, and the genes rpob2 and rbpA predicted to confer resistance to rifampin. In comparison to 129 genomes of the same species, strain Cp162 has four genes exclusively present, two of them code transposases and two truncated proteins, and the three exclusively absent genes lysG, NUDIX domain protein, and Hypothetical protein. All 130 genomes had the rifampin resistance genes rpob2 and rbpA. Our results found no unique gene that could be associated with tropism to camel host, and further studies should include more genomes and genome-wide association studies testing for genes and SNPs.

Molecular Characterization of Corynebacterium pseudotuberculosis Isolated over a 15-Year Period in Switzerland

Corynebacterium pseudotuberculosis biovar Ovis is the etiological agent of the contagious and chronic disease caseous lymphadenitis (CLA) in sheep and goats. The economic impact of CLA in Switzerland remains largely unknown, and the transmission modalities, as well as the genetic diversity of circulating strains, are poorly understood. This work presents further characterization data for 215 C. pseudotuberculosis isolates from sheep, goats and a dromedary originating from Switzerland and the Principality of Liechtenstein, collected over a 15-year period. The isolates were classified into the two biovars Ovis and Equi, analyzed for the presence of the diphtheria-like toxin gene and characterized using MLSA. All sheep and goat isolates were classified as C. pseudotuberculosis biovar Ovis. The isolate from a dromedary was classified as biovar Equi. No isolates harboring the diphtheria-like toxin gene were detected. Phylogenetic analysis of the concatenated sequences of four genes revealed the existence of 24 clusters. There was no correlation between MLSA sequence types, year of isolation and the geographical origin of the isolates. These findings confirm the presence of several MLSA sequence types in the study area and over a 15-year period. Moreover, no sheep- and goat-specific MLSA sequence types were found.

Cell wall glycolipids from Corynebacterium pseudotuberculosis strains with different virulences differ in terms of composition and immune recognition

Caseous lymphadenitis (CLA) is an infectious disease caused by Corynebacterium pseudotuberculosis in small ruminants and is characterized by the development of granulomas in the lymph nodes, spleen, liver, and lungs. Although little is known about the host-pathogen relationship of this bacterium, it was previously reported that the pathogen's lipids are important for its taxonomic classification and survival inside macrophages. However, there are no studies regarding the composition of these molecules. In this study, cell wall glycolipids from two C. pseudotuberculosis strains presenting different virulence profiles were purified and its composition was characterized. A difference was observed between the electrophoretic and chromatogram profiles for cell wall components from the two strains, mainly among molecules with low molecular weights. IgM from sheep with acute CLA recognized antigens with an estimated molecular weight of 11 kDa of the low-pathogenicity strain, while low-molecular weight antigens from the high-pathogenicity strain presented a lower recognition by these antibodies. Mass spectrometry analysis showed that the cell wall of the high-pathogenicity strain contained glycolipids with high amounts of unsaturated fatty acids and glycerophosphoinositols, which may contribute to the capacity of this strain to cause severe disease. In conclusion, it is indicated that cell wall non-protein antigens can play a key role in C. pseudotuberculosis virulence.

Improved reconstruction and comparative analysis of chromosome 12 to rectify Mis-assemblies in Gossypium arboreum

BACKGROUND

Genome sequencing technologies have been improved at an exponential pace but precise chromosome-scale genome assembly still remains a great challenge. The draft genome of cultivated G. arboreum was sequenced and assembled with shotgun sequencing approach, however, it contains several misassemblies. To address this issue, we generated an improved reassembly of G. arboreum chromosome 12 using genetic mapping and reference-assisted approaches and evaluated this reconstruction by comparing with homologous chromosomes of G. raimondii and G. hirsutum.

RESULTS

In this study, we generated a high quality assembly of the 94.64 Mb length of G. arboreum chromosome 12 (A_A12) which comprised of 144 scaffolds and contained 3361 protein coding genes. Evaluation of results using syntenic and collinear analysis of reconstructed G. arboreum chromosome A_A12 with its homologous chromosomes of G. raimondii (D_D08) and G. hirsutum (AD_A12 and AD_D12) confirmed the significant improved quality of current reassembly as compared to previous one. We found major misassemblies in previously assembled chromosome 12 (A_Ca9) of G. arboreum particularly in anchoring and orienting of scaffolds into a pseudo-chromosome. Further, homologous chromosomes 12 of G. raimondii (D_D08) and G. arboreum (A_A12) contained almost equal number of transcription factor (TF) related genes, and showed good collinear relationship with each other. As well, a higher rate of gene loss was found in corresponding homologous chromosomes of tetraploid (AD_A12 and AD_D12) than diploid (A_A12 and D_D08) cotton, signifying that gene loss is likely a continuing process in chromosomal evolution of tetraploid cotton.

CONCLUSION

This study offers a more accurate strategy to correct misassemblies in sequenced draft genomes of cotton which will provide further insights towards its genome organization.