Average GC Content Research Articles

Isolation and complete genome sequencing of the virulent phage vB_EcoS_XY3 infecting multidrug-resistant Escherichia coli.

A virulent phage, named vB_EcoS_XY3, was isolated from hospital wastewater in Xiangyang, China. Its morphological characteristics, growth parameters, adsorption rate, and pH and temperature stability were determined. Phage vB_EcoS_XY3 was found to be able to infect Escherichia coli laboratory strains and also some multidrug-resistant E. coli strains. Its complete genome consists of 51,345 base pairs of double-stranded DNA with an average GC content of 55.24% and 85 putative protein-coding genes. Forty-four genes were annotated with known functions. These results will not only provide further insights into E. coli phages but also have implications for the development of potential biocontrol agents.

Complete Genome Sequence of the Acetic Acid Bacterium Acetobacter aceti NBRC 14818.

We report here the complete genome sequence of the acetic acid bacterium Acetobacter aceti type strain NBRC 14818. The genome comprises a chromosome of 3,596,270 bp with 57.1% GC content and four plasmids/phages of 63,279 bp, 25,755 bp, 4,858 bp, and 2,964 bp, with an average GC content of 57.0%.

Complete Chloroplast Genome of Argania spinosa: Structural Organization and Phylogenetic Relationships in Sapotaceae.

Argania spinosa (Sapotaceae), an important endemic Moroccan oil tree, is a primary source of argan oil, which has numerous dietary and medicinal proprieties. The plant species occupies the mid-western part of Morocco and provides great environmental and socioeconomic benefits. The complete chloroplast (cp) genome of A. spinosa was sequenced, assembled, and analyzed in comparison with those of two Sapotaceae members. The A. spinosa cp genome is 158,848 bp long, with an average GC content of 36.8%. The cp genome exhibits a typical quadripartite and circular structure consisting of a pair of inverted regions (IR) of 25,945 bp in length separating small single-copy (SSC) and large single-copy (LSC) regions of 18,591 and 88,367 bp, respectively. The annotation of A. spinosa cp genome predicted 130 genes, including 85 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. A total of 44 long repeats and 88 simple sequence repeats (SSR) divided into mononucleotides (76), dinucleotides (7), trinucleotides (3), tetranucleotides (1), and hexanucleotides (1) were identified in the A. spinosa cp genome. Phylogenetic analyses using the maximum likelihood (ML) method were performed based on 69 protein-coding genes from 11 species of Ericales. The results confirmed the close position of A. spinosa to the Sideroxylon genus, supporting the revisiting of its taxonomic status. The complete chloroplast genome sequence will be valuable for further studies on the conservation and breeding of this medicinally and culinary important species and also contribute to clarifying the phylogenetic position of the species within Sapotaceae.

Demographic history and adaptive synonymous and nonsynonymous variants of nuclear genes in Rhododendron oldhamii (Ericaceae)

Demographic events are important in shaping the population genetic structure and exon variation can play roles in adaptive divergence. Twelve nuclear genes were used to investigate the species-level phylogeography of Rhododendron oldhamii, test the difference in the average GC content of coding sites and of third codon positions with that of surrounding non-coding regions, and test exon variants associated with environmental variables. Spatial expansion was suggested by R2 index of the aligned intron sequences of all genes of the regional samples and sum of squared deviations statistic of the aligned intron sequences of all genes individually and of all genes of the regional and pooled samples. The level of genetic differentiation was significantly different between regional samples. Significantly lower and higher average GC contents across 94 sequences of the 12 genes at third codon positions of coding sequences than that of surrounding non-coding regions were found. We found seven exon variants associated strongly with environmental variables. Our results demonstrated spatial expansion of R. oldhamii in the late Pleistocene and the optimal third codon position could end in A or T rather than G or C as frequent alleles and could have been important for adaptive divergence in R. oldhamii.

Full-length transcriptome sequencing analysis and development of EST-SSR markers for the endangered species Populus wulianensis

Populus wulianensis is an endangered species endemic to Shandong Province, China. Despite the economic and ornamental value of this species, few genomics and genetic studies have been performed. In this study, we performed a relevant analysis of the full-length transcriptome sequencing data of P. wulianensis and obtained expressed sequence tag (EST)-simple sequence repeat (SSR) markers with polymorphisms that can be used for further genetic research. In total, 8.18 Gb (3,521,665) clean reads with an average GC content of 42.12% were obtained. From the corrected 64,737 high-quality isoforms, 42,323 transcript sequences were obtained after redundancy analysis with CD-HIT. Among these transcript sequences, 41,876 sequences were annotated successfully. A total of 23,539 potential EST-SSRs were identified from 16,057 sequences. Excluding mononucleotides, the most abundant motifs were trinucleotide SSRs (47.80%), followed by di- (46.80%), tetra- (2.98%), hexa- (1.58%) and pentanucleotide SSRs (0.84%). Among the 100 designed EST-SSRs, 18 were polymorphic with high PIC values (0.721 and 0.683) and could be used for analyses of the genetic diversity and population structure of P. wulianensis. These full-length transcriptome sequencing data will facilitate gene discovery and functional genomics research in P. wulianensis, and the novel EST-SSRs developed in our study will promote molecular-assisted breeding, genetic diversity and conservation biology research in this species.

Analysis of six chloroplast genomes provides insight into the evolution of Chrysosplenium (Saxifragaceae)

BackgroundChrysosplenium L. (Saxifragaceae) is a genus of plants widely distributed in Northern Hemisphere and usually found in moist, shaded valleys and mountain slopes. This genus is ideal for studying plant adaptation to low light conditions. Although some progress has been made in the systematics and biogeography of Chrysosplenium, its chloroplast genome evolution remains to be investigated.ResultsTo fill this gap, we sequenced the chloroplast genomes of six Chrysosplenium species and analyzed their genome structure, GC content, and nucleotide diversity. Moreover, we performed a phylogenetic analysis and calculated non-synonymous (Ka) /synonymous (Ks) substitution ratios using the combined protein-coding genes of 29 species within Saxifragales and two additional species as outgroups, as well as a pair-wise estimation for each gene within Chrysosplenium. Compared with the outgroups in Saxifragaceae, the six Chrysosplenium chloroplast genomes had lower GC contents; they also had conserved boundary regions and gene contents, as only the rpl32 gene was lost in four of the Chrysosplenium chloroplast genomes. Phylogenetic analyses suggested that the Chrysosplenium separated to two major clades (the opposite group and the alternate group). The selection pressure estimation (Ka/Ks ratios) of genes in the Chrysosplenium species showed that matK and ycf2 were subjected to positive selection.ConclusionThis study provides genetic resources for exploring the phylogeny of Chrysosplenium and sheds light on plant adaptation to low light conditions. The lower average GC content and the lacking gene of rpl32 indicated selective pressure in their unique habitats. Different from results previously reported, our selective pressure estimation suggested that the genes related to photosynthesis (such as ycf2) were under positive selection at sites in the coding region.

Genome Mining Reveals the Biosynthetic Pathways of Polyhydroxyalkanoate and Ectoines of the Halophilic Strain Salinivibrio proteolyticus M318 Isolated from Fermented Shrimp Paste.

Salinivibrio proteolyticus M318, a halophilic bacterium isolated from fermented shrimp paste, is able to produce polyhydroxyalkanoate (PHA) from different carbon sources. In this study, we report the whole-genome sequence of strain M138, which comprises 2 separated chromosomes and 2 plasmids, and the complete genome contains 3,605,935bp with an average GC content of 49.9%. The genome of strain M318 contains 3341 genes, 98 tRNA genes, and 28 rRNA genes. The 16S rRNA gene sequence and average nucleotide identity analysis associated with morphological and biochemical tests showed that this strain has high homology to the reference strain Salinivibrio proteolyticus DSM 8285. The genes encoding key enzymes for PHA and ectoine synthesis were identified from the bacterial genome. In addition, the TeaABC transporter responsible for ectoine uptake from the environment and the operon doeABXCD responsible for the degradation of ectoine were also detected. Strain M318 was able to produce poly(3-hydroxybutyrate) [P(3HB)] from different carbon sources such as glycerol, maltose, glucose, fructose, and starch. The ability to produce ectoines at different NaCl concentrations was investigated. High ectoine content of 26.2% of cell dry weight was obtained by this strain at 18% NaCl. This report provides genetic information regarding adaptive mechanisms of strain M318 to stress conditions, as well as new knowledge to facilitate the application of this strain as a bacterial cell factory for the production of PHA and ectoine.

Draft Genome Sequence of Enterobacter roggenkampii Strain OS53, Isolated from Corroded Pipework at an Offshore Oil Production Facility.

Here, we report the genome sequence of Enterobacter roggenkampii strain OS53, isolated from corroded pipework at an offshore oil production facility. The draft genome sequence comprises 6 contigs and contains 5,194,507 bp with an average GC content of 55.90%.

Complete genome of a unicellular parasite (Antonospora locustae) and transcriptional interactions with its host locust.

Microsporidia are a large group of unicellular parasites that infect insects and mammals. The simpler life cycle of microsporidia in insects provides a model system for understanding their evolution and molecular interactions with their hosts. However, no complete genome is available for insect-parasitic microsporidian species. The complete genome of Antonospora locustae, a microsporidian parasite that obligately infects insects, is reported here. The genome size of A. locustae is 3 170 203 nucleotides, composed of 17 chromosomes onto which a total of 1857 annotated genes have been mapped and detailed. A unique feature of the A. locustae genome is the presence of an ultra-low GC region of approximately 25 kb on 16 of the 17 chromosomes, in which the average GC content is only 20 %. Transcription profiling indicated that the ultra-low GC region of the parasite could be associated with differential regulation of host defences in the fat body to promote the parasite’s survival and propagation. Phylogenetic gene analysis showed that A. locustae, and the microsporidian family in general, is likely at an evolutionarily transitional position between prokaryotes and eukaryotes, and that it evolved independently. Transcriptomic analysis showed that A. locustae can systematically inhibit the locust phenoloxidase PPO, TCA and glyoxylate cycles, and PPAR pathways to escape melanization, and can activate host energy transfer pathways to support its reproduction in the fat body, which is an insect energy-producing organ. Our study provides a platform and model for studies of the molecular mechanisms of microsporidium–host interactions in an energy-producing organ and for understanding the evolution of microsporidia.

Analysis on the virulomes and resistomes of multi-drug resistance clinical Escherichia coli isolates, as well as the interactome with gut microbiome

Escherichia coli is one of the most diverse microbial species. Pathogenic E. coli is capable of causing various diseases in humans, including several types of diarrhea, urinary tract infections, sepsis, and meningitis. This study focused on the antibiotic susceptibility profile and genomic analysis of a clinical E. coli Guangzhou-Eco330 isolated from a hospitalized 8-year-old female patient suffered from pulmonary infection in 2017. Susceptibility to 15 antibiotics were determined using Vitek2™ Automated Susceptibility System and Etest strips and interpreted based on CLSI guidelines. The genome was sequenced using Illumina Hiseq 2500 platform and assembled de novo using Velvet, followed by bioinformatics analysis. The genome has a length of 5,132,642 bp and contains 4989 predicted genes with an average GC content of 50.51%. The carriage of rfbE gene suggested the strain belonging to O157. In the genome, 70 non-coding RNAs, 50 repeat sequences, 18 transposons, 78 GIs, 9 CRISPRs, and 3 large prophages were identified. 37 PHI related genes and 108 virulence genes were determined to contribute to its pathogenicity. Specifically, the acquisition of multiple antibiotic resistance genes including blaCTX-M-55, blaOXA-10, blaCMY-48, tetB, and qnrS1 contributed to its resistance to penicillins, telracyclines, cephalosporin, and quinolones. The understanding of the genome may aid in further study on the clinical control of multi-drug resistance E. coli.

Identification of species in the genus Nitraria L. (Nitrariaceae) based on nucleotide variability of nuclear ribosomal DNA.

Intragenomic polymorphism of ITS1 and ITS2 of nuclear ribosomal DNA sequences was analysed in 33 samples belonging to the Nitraria species N. schoberi, N. sibirica, and N. komarovii. The nucleotide variability of the ITS region was detected in the Nitraria species as single-nucleotide substitutions (mainly transitions) and single-nucleotide deletion. Information about the nucleotide variability of fragments is given for the first time by us. The ITS1-5.8S-ITS2 region contained 17 phylogenetically informative single-nucleotide polymorphisms. Eleven single-nucleotide substitutions (transitions, C/T) were detected in ITS1. The ITS2 spacer contained 273–274 bp and was more conservative. A total of 5 phylogenetically informative single-nucleotide polymorphisms (4 transitions: C/T, G/A, one transversion: G/C), one single-nucleotide deletion (T/–) were detected in ITS2. The average GC content was 61.5 %. The GC content was lower in N. sibirica (59.2 %) than in N. schoberi and N. komarovii (62.7 %). It has been shown that the shorter ITS2 is a suitable molecular marker separating these species, due to the low interspecific variability and simultaneous available intraspecific variability. Phylogenetic ML and BI trees constructed separately for the ITS1 and ITS2 spacers, as well as separately for the full-size ITS region and the ITS2 spacer, were congruent. The results obtained on the intraspecific differentiation of N. sibirica revealed two main ribotypes among the samples of this species: the main Siberian sibirica-ribotype and the main Kazakh sibiricaribotype. Geographical features of the distribution of N. sibirica ribotypes, as well as the presence of significant differences between the main Siberian and Kazakh sibirica-ribotypes (3 single-nucleotide substitutions) indicated significant inter-population differences and taxonomic heterogeneity of N. sibirica. Most likely, the processes of homogenization of nuclear ribosomal DNA of N. sibirica samples, the origin of which is associated with hybridization and speciation, are currently continuing.

Complete genome assembly data of paenibacillus sp. RUD330, a hypothetical symbiont of euglena gracilis.

An unknown bacterial strain was detected in the cytostome of Euglena gracilis and on the cell surface of Euglena gracilis using transmission electron microscopy. To identify the unknown bacterium and its function, we performed isolation experiments. Here we present the genome sequence of the isolate that was determined to be Paenibacillus sp. The genome of the bacterium was sequenced four times using Illumina technology with pair-end reads, Illumina technology with mate pair reads (inserts 3–4 and 6–8 Kb), and Nanopore technology with long reads (tens of thousands of nucleotides). Assemblies based on Illumina reads including mate-pair reads could not resolve issues caused by long tandem copies of rRNA, other tandem repeats, and extremely GC-rich regions (90–100%). Only long Nanopore reads resolved those gaps and made it possible to complete the entire genome; moreover, we found one plasmid. The length of the genome is 5.56 Mbp, and the average GC content is 59%. The genome of Paenibacillus sp. RUD330 included 8 copies of all the rRNA genes (23S; 16S; 5S), the length of the plasmid was 8.3 Kb.We hope that our genome assembly and the methods used can help other investigators in the assembly of complex genomes. Our reliable assembly could be a good basis for further physiological and genetic engineering studies of similar strains.

Lentzea indica sp. nov., a novel actinobacteria isolated from Indian Himalayan-soil.

A novel actinobacterium, strain PSKA42T was isolated from a soil sample of Kashmir Himalaya, India (latitude: 34°-01' N; longitude: 74°-47' E; altitude: 5328 ft). The morphological, biochemical and molecular characteristics predicted that PSKA42T belongs to the genera Lentzea. Using the Ezbiocloud server it has also been manifested that the 16S rRNA gene sequence of PSKA42 was 98.4% similar with Lentzea guizhouensis DHSC013T, Lentzea albidocapillata subsp. violacea IMSNU50388T; 98.26% similar with Lentzea aerocolonigenes NRRLB-3298T, Lentzea albidocapillata subsp. albidocapillata NRRLB-24057T and 98.19% similar with Lentzea californiensis DSM43393T. Whole-genome average nucleotide identity values between PSKA42T and L. guizhouensis DHSC013T, L. albidocapillata subsp. violaera IMSNU50388T, L. aerocolonigenes NRRLB-3298T, L. albidocapillata subsp. albidocapillata NRRLB-24057T, L. californiensis DSM43393T, L. flaviverucosa As40578T were 84.94%, 85.8%, 87.15%, 87.71%, 79.29, and 87.22%, respectively, suggesting that PSKA42T represented a new species. In-silco genome-to-genome distances analysis and MLSA with selected housekeeping genes also indicated that the isolate should be assigned to a new species under the genus Lentzea. Analysis of the whole-cell hydrolysate found mannose and ribose as major sugar while cell wall amino acid was identified as meso-diaminopimelic acid. Major fatty acids (> 5%) were observed as 15:0 iso (23.81%), and 16:0 iso (26.11%). PSKA42 contains polyketide synthase type-I genes but unable to show antimicrobial activity in laboratory conditions. The average GC content in its genome is 68.3%. Further, based on physicochemical, chemotaxonomic markers and molecular analysis, PSKA42T can be assigned as novel species of Lentzea, we proposed Lentzea indica sp nov. The type strain is PSKA42T (JCM 33729T, MTCC 12936T, MCC 4127T).

Complete genome sequence of the aromatic-hydrocarbon-degrading bacterium Achromobacter xylosoxidans DN002.

Achromobacter xylosoxidans DN002 is capable of utilizing numerous aromatic hydrocarbons as sole carbon and energy resource. In this study, the whole genome of strain DN002 was sequenced and analyzed, which consisted of one circular chromosome of 5,943,204bp and a 278,917bp plasmid with an average GC content of 65.46mol%, 5694 protein-coding genes, 13 rRNA genes and 57 tRNA genes. Analysis of cluster of orthologous group (COG) demonstrated that strain DN002 had remarkable gene abundance foramino acid transport and metabolism, transcription, inorganic ion transport and metabolism, energy production and conversion, and carbohydrate transport and metabolism. Genes related to biodegradation of aromatic hydrocarbons, chemotaxis and flagella were identified from the genome, which will advance our fundamental understanding the molecular mechanism for degradation and metabolizing of aromatic hydrocarbons.

Complete mitochondrial DNA sequence of the yeast Zygosaccharomyces siamensis (Saccharomycetes: Saccharomycetales) from fermented honey of the Apis cerana japonica in Japan

Japanese Apis cerana honey is believed to ferment due to the presence of some yeast. We analyzed the complete mitochondrial genome of Zygosaccharomyces siamensis, a yeast involved in fermenting honey. Zygosaccharomyces siamensis was obtained from the honey collected from A. cerana hives under traditional beekeeping in the forest of the Minabe-Tanabe Ume system in Wakayama Prefecture, Japan. Its mitochondrial genome was identified as a 23,184 bp circular molecule containing 8 protein-coding genes (PCGs), 24 tRNA genes, and 2 rRNA genes. The PCGs contained a common set of genes encoding ATPase subunits (ATP6, ATP8, and ATP9), three subunits of cytochrome C oxidase (COX1, COX2, and COX3), apocytochrome b (CYTB), and ribosome-associated protein (RPS3). In addition, two open-reading frames coding for LAGLIDADG endonucleases were predicted to be about 1100 bps. The average GC content was found to be 48.6%. The heavy strand was predicted to have 7 PCGs, 22 tRNA genes, and 2 rRNA genes, while the light strand was predicted to contain one PCG and two tRNA genes. Molecular phylogenetic analyses of the mitochondrial DNA genes strongly supported the result obtained from the phylogenetic analysis of partial ITS region sequences, grouping the monophyletic species within the genus Zygosaccharomyces. The complete mitochondrial DNA genome of this honey-fermenting yeast will provide useful information for understanding the basis of the honey fermentation process.

Characterization of the complete chloroplast genome of Phalaenopsis lobbii (Orchidaceae), an important horticultural plant in China

ABSTACT An important horticultural moth orchid in China, Phalaenopsis lobbii, has small and lovely flowers with unique labial characteristics. We sequenced the chloroplast genome (cpDNAs) of P. lobbii and analyzed the structural characteristics of its complete cpDNAs sequence. The results showed that P. lobbii had a typical circular double-stranded tetrad structure with the average GC content of 36.60%, and its total length was 144,607 bp. The complete chloroplast was annotated with 125 genes, eight rRNAs and 38 tRNAs. Phylogenetic tree displayed that P. lobbii was closely related to P. equestris, P. aphrodite subsp. formosana, P. ‘Tiny Star’ and P. japonica. This published cpDNAs of P. lobbii will provide information for future phylogenetic and artificial breeding of Phalaenopsis.

Chloroplast characterizations of a Phalaenopsis native to China, Phalaenopsis mannii (Orchidaceae)

ABSTACT Phalaenopsis mannii, one of the native Phalaenopsis in China, is an important parent for breeding new varieties. However, its position has been unclear in Phalaenopsis. The obtained high-quality P. mannii chloroplast genome will provide useful information for phylogenetic and future breeding of Phalaenopsis. Herein, we reported a complete chloroplast genome of P. mannii from Yunnan, China. The sequencing data obtained from BGISEQ-500 platform were assembled. This sequence had a circular molecular length of 148,596 bp and contained a total of 127 genes with an average GC content of 36.7%. Phylogenetic analysis indicated that Phalaenopsis was monophyletic with strong support, in which the P.mannii was the sister-group of Phalaenopsis aphrodite subsp. formosas, Phalaenopsis ‘TinyStar’ and Phalaenopsis equestris.

Dataset of complete genome assembly and analysis of mycobacterium tuberculosis strain SIT745/EAI1-MYS.

In this dataset, we report the genome assembly and data analysis of Mycobacterium tuberculosis strain SIT745/EAI1-MYS. Previously, this strain was isolated from a Malaysian patient with extra-pulmonary tuberculosis, and identification of this strain is done by spoligotype patterns with fifteen known Shared International Type (SITs). Further analysis showed that this strain has a remarkable phylogeographical specificity for Malaysia. Based on the National Center for Biotechnology Information (NCBI) nucleotide database information, the complete genome consists of 150 contigs with various sequence lengths and was not assembled. In this assembly, the aforementioned contigs along with reference sequence from Mycobacterium tuberculosis strain H37Rv and Mycobacterium bovis strain AF2122/97 was used for gap closures, were assembled into a single circular chromosome length of approximately 4.42 Mega bases (Mb) with an average GC content of 65.6%. The single circular chromosome was shown to contain 4,009 protein-coding sequences, 3 ribosomal RNAs, 45 transfer RNAs, and 12 superclasses distributed with 277 subsystems which constitute nearly 1900 genes, respectively. The genome information will provide fundamental knowledge of this organism as well as insight for understanding genomic and proteomic profiling, phylogenetic relationship.

Genome Sequence Resources of Two 'Candidatus Liberibacter asiaticus' Strains from Pakistan.

'Candidatus Liberibacter asiaticus' (CLas) is an unculturable, phloem-restricted αProteobacteria, associated with citrus Huanglongbing (HLB), which is one of the most destructive diseases in citrus production worldwide. Here, we present the genome sequences of CLas strains PA19 and PA20 from HLB-affected kinnow trees in Multan, Punjab Province, Pakistan. The CLas genomes of PA19 and PA20 comprise 1,224,156 bp and 1,226,225 bp, respectively, with an average GC content of 36.4%. Both harbored the Type 2 prophage. In this study, we report two CLas genomes from Pakistan, which extends the sequence database of CLas and will contribute to CLas biology and HLB management.

Comparative Genome Analysis of a Pathogenic Erysipelothrix rhusiopathiae Isolate WH13013 from Pig Reveals Potential Genes Involve in Bacterial Adaptions and Pathogenesis.

Erysipelothrix rhusiopathiae is a common pathogen responsible for pig erysipelas. However, the molecular basis for the pathogenesis of E. rhusiopathiae remains to be elucidated. In this study, the complete genome sequence of the E. rhusiopathiae strain WH13013, a pathogenic isolate from a diseased pig, was generated using a combined strategy of PacBio RSII and Illumina sequencing technologies. The strategy finally generated a single circular chromosome of approximately 1.78 Mb in size for the complete genome of WH13013, with an average GC content of 36.49%. The genome of WH13013 encoded 1633 predicted proteins, 55 tRNAs, as well as 15 rRNAs. It contained four genomic islands and several resistance-associated genes were identified within these islands. Phylogenetic analysis revealed that WH13013 was close to many other sequenced E. rhusiopathiae virulent strains. The comprehensive comparative analysis of eight E. rhusiopathiae virulent strains, including WH13013, identified a total of 1184 core genes. A large proportion (approximately 75.31%) of these core genes participated in nutrition and energy uptake and metabolism as well as the other bioactivities that are necessary for bacterial survival and adaption. The core genes also contained those encoding proteins participating in the biosynthesis and/or the components of the proposed virulence factors of E. rhusiopathiae, including the capsule (cpsA, cpsB, cpsC), neuraminidase (nanH), hyaluronidase (hylA, hylB, hylC), and surface proteins (spaA, rspA, rspB). The obtaining of the complete genome sequence of this virulent strain, WH13013, and this comprehensive comparative genome analysis will help in further studies of the genetic basis of the pathogenesis of E. rhusiopathiae.