Genetic relationships of Bacillus anthracis and closely related species based on variable-number tandem repeat analysis and BOX-PCR genomic fingerprinting

Variable-number tandem repeats (VNTR) analysis and BOX-repeat-based PCR (BOX-PCR) genomic fingerprinting were performed on 25 Bacillus strains to investigate the genetic relatedness of Bacillus anthracis to the closely related species. Based on VNTR analysis, all B. anthracis strains could be assigned to (VNTR)(4), which is the most commonly found type in the world. Interestingly, a (VNTR)(2) was also observed in Bacillus cereus KCTC 1661 and with an exact match to the tandem repeats found in B. anthracis. This finding has never been reported before in the closely related species. According to the BOX-PCR, B. anthracis strains clustered together and separated reliably from the closely related species. However, B. cereus KCTC 1661 was linked to the B. anthracis cluster and showed close relationships with B. anthracis strains. These results indicated that there was a strong correlation between VNTR analysis and BOX-PCR genomic fingerprinting.     

Multiple-locus variable-number tandem repeat analysis reveals genetic relationships within Bacillus anthracis

Bacillus anthracis is one of the most genetically homogeneous pathogens described, making strain discrimination particularly difficult. In this paper, we present a novel molecular typing system based on rapidly evolving variable-number tandem repeat (VNTR) loci. Multiple-locus VNTR analysis (MLVA) uses the combined power of multiple alleles at several marker loci. In our system, fluorescently labeled PCR primers are used to produce PCR amplification products from eight VNTR regions in the B. anthracis genome. These are detected and their sizes are determined using an ABI377 automated DNA sequencer. Five of these eight loci were discovered by sequence characterization of molecular markers (vrrC(1), vrrC(2), vrrB(1), vrrB(2), and CG3), two were discovered by searching complete plasmid nucleotide sequences (pXO1-aat and pXO2-at), and one was known previously (vrrA). MLVA characterization of 426 B. anthracis isolates identified 89 distinct genotypes. VNTR markers frequently identified multiple alleles (from two to nine), with Nei's diversity values between 0.3 and 0.8. Unweighted pair-group method arithmetic average cluster analysis identified six genetically distinct groups that appear to be derived from clones. Some of these clones show worldwide distribution, while others are restricted to particular geographic regions. Human commerce doubtlessly has contributed to the dispersal of particular clones in ancient and modern times.     

Mass spectrometry provides accurate characterization of two genetic marker types in Bacillus anthracis

Epidemiological and forensic analyses of bioterrorism events involving Bacillus anthracis could be improved if both variable number tandem repeats (VNTRs) and single nucleotide polymorphisms (SNPs) could be combined on a single analysis platform. Here we present the use of electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) to characterize 24 alleles from 6 VNTR loci and 11 alleles from 7 SNP loci in B. anthracis. The results obtained with ESI-FTICR-MS were consistent with independent results obtained from traditional approaches using electrophoretic detection of fluorescent products. However, ESI-FTICR-MS improves on the traditional approaches because it does not require fluorescent labeling of PCR products, minimizes post-PCR processing, obviates electrophoresis, and provides unambiguous base composition of both SNP and VNTR PCR products. In addition, ESI-FTICR-MS allows both marker types to be examined simultaneously and at a rate of approximately 1 sample per min. This technology represents a significant advance in our ability to rapidly characterize B. anthracis isolates using VNTR and SNP loci.     

Molecular diversity in Bacillus anthracis

Molecular typing of Bacillus anthracis has been extremely difficult due to the lack of polymorphic DNA markers. We have identified nine novel variable number tandemly repeated loci from previously known amplified fragment length polymorphism markers or from the DNA sequence. In combination with the previously known vrrA locus, these markers provide discrimination power to genetically characterize B. anthracis isolates. The variable number tandem repeat (VNTR) loci are found in both gene coding (genic) and non-coding (non-genic) regions. The genic differences are 'in frame' and result in additions or deletion of amino acids to the predicted proteins. Due the rarity of molecular differences, the VNTR changes represent a significant portion of the genetic variation found within B. anthracis. This variation could represent an important adaptive mechanism. Marker similarity and differences among diverse isolates have identified seven major diversity groups that may represent the only world-wide B. anthracis clones. The lineages reconstructed using these data may reflect the dispersal and evolution of this pathogen.     

Characterization of the variable-number tandem repeats in vrrA from different Bacillus anthracis isolates.

PCR analysis of 198 Bacillus anthracis isolates revealed a variable region of DNA sequence differing in length among the isolates. Five polymorphisms differed by the presence of two to six copies of the 12-bp tandem repeat 5'-CAATATCAACAA-3'. This variable-number tandem repeat (VNTR) region is located within a larger sequence containing one complete open reading frame that encodes a putative 30-kDa protein. Length variation did not change the reading frame of the encoded protein and only changed the copy number of a 4-amino-acid sequence (QYQQ) from 2 to 6. The structure of the VNTR region suggests that these multiple repeats are generated by recombination or polymerase slippage. Protein structures predicted from the reverse-translated DNA sequence suggest that any structural changes in the encoded protein are confined to the region encoded by the VNTR sequence. Copy number differences in the VNTR region were used to define five different B. anthracis alleles. Characterization of 198 isolates revealed allele frequencies of 6.1, 17.7, 59.6, 5.6, and 11.1% sequentially from shorter to longer alleles. The high degree of polymorphism in the VNTR region provides a criterion for assigning isolates to five allelic categories. There is a correlation between categories and geographic distribution. Such molecular markers can be used to monitor the epidemiology of anthrax outbreaks in domestic and native herbivore populations.     

Phenotypic and genetic features of cultural-morphologic variants of Bacillus anthracis

Comparative analysis of MVLA-genotypes of 6 Bacillus anthracis strains and 40 their variants differing on capsule- and toxin synthesis, hemolytic, proteolytic and lecitinase activity, nutritional requirements, susceptibility to anthrax bacteriophages, virulence, immunogenicity, and presence of genes for capsule and toxin synthesis was performed. Results of phylogenetic analysis of 5 chromosome locuses and plasmid locus pXO1aat which are variable for this sample of B. anthracis cultures showed that all strains divided on 2 main clusters - A and B. Cluster A consisted of 5 genotypes whereas cluster B - of 1 genotype. All highly virulent original strains and variants with characteristic phenotype Cap(CO2)(+)(O2)(-)Tox(+)ProtA(+)Hly(+) Lec(-)Trp(+) had identical genotype in 4 groups and in 5th group differences were present only in vrrA locus. All original strains and variants with the most atypical complex of phenotypic characteristics Cap (CO2)(+)(O2)(+)Tox(-)ProtA(-)Hly(-)Lec(-)Trp(-) also had the same genotype belonging to cluster B and diverged on characteristic of 5 chromosomal VNTR locuses and pXO1aat locus from typical strains. Absence of toxin production in vitro was not related to loss of genetic determinants of toxin components. Cultures with typical characteristics, one of which was ability to produce toxin in vitro, had larger sizes of amplicons of pXO1aat locus (135 and 132 nbp), whereas atoxigenic original strains and variants with complex of atypical characteristics and identical chromosome genotype had the smallest sizes (123 bnp). All original cultures were isolated in Russia, their genotypes are described for the first time.     

A subset of skin-expressed microRNAs with possible roles in goat and sheep hair growth based on expression profiling of mammalian microRNAs

The microRNAs (miRNAs) are an extensive class of small noncoding RNAs (18-25 nucleotides) with probable roles in the regulation of gene expression. Due to the fact that miRNAs are conserved in closely related eukaryotes and some are also conserved across a larger evolutionary distance, their potential functions in mammalian development are under active study. In order to identify mammalian miRNAs that might function in hair growth, we characterized the expression of 159 miRNAs in adult body side skin and ear skin from goat and sheep using microarray analysis. Of these, 19 miRNAs were specifically expressed or greatly enriched in body side skin in goats and sheep. This suggests hair growth-specific functions for miRNAs. Of the coexpressed 105 miRNAs, the degree of correlation within species is higher than interspecies. Nine of the expressed miRNAs were detected exclusively in the goat body side skin area where more cashmere was grown than coat hair; mmu-miR-720 and mmu-miR-199b were expressed primarily in goat skin. The identification of 105 of skin-expressed miRNAs whose expression behavior is conserved in goats and sheep differentiating hair follicles implicates these miRNAs have functions in mammalian hair follicles growth and development. We demonstrate that miRNAs previously associated with hair follicles in the mouse are also expressed in the adult skin of goats and sheep. In addition, 69 more conserved miRNAs cross-species were discerned in the study. Of them, the let-7, mir-17, mir-30, mir-15, and mir-8 gene families were expressed in high frequency. These results reveal critical roles of them in skin and hair follicle development and function.     

High-throughput sequencing of hair follicle development-related micrornas in cashmere goat at various fetal periods

Inner Mongolia cashmere goat marks a precious gerplasm genetic resource due to its excellent cashmere traits. Therefore, it is of crucial importance to investigate the cashmere development mechanism of cashmere goat and to search for the important cashmere growth-related candidate genes. Fetal skin samples at 10 different periods of cashmere goat were collected in this research. Moreover, high-throughput sequencing was conducted on RNA samples from side skin of cashmere goat fetuses collected at three critical periods of skin hair follicle initiation, growth and development (namely, 45, 55 and 65?days) after balanced mix in line with the previous research results. Meanwhile, 3 samples at corresponding periods were used as the biological duplications. Data regarding microRNA and mRNA expression in skin and hair follicles of cashmere goats at various fetal periods were obtained using the high-throughput sequencing method. The results indicated that microRNAs in the oar-let-7 and oar-miR-200 families in 55?days and 66?days of pregnancy samples had been notably up-regulated relative to those in 45?days of pregnancy samples. This revealed that they might be the critical microRNAs in hair follicle development.     

Establishment and characterization of outer root sheath (ORS) cell line from Jining grey goat

A new line of outer root sheath (ORS) cells was established from hair follicles of Jining grey goat by using a mechanical separation combined with enzyme digestion. Cell morphology is described at different phases. The chromosome analysis of ORS cells, identification of the ORS cells and morphological reversion test were detected at the 4th and 40th passages. The ORS cells were healthy and the growth characteristics were stable with a population doubling time of 52 h. Chromosome analysis showed that >58% of cells were diploid. Test for ORS cell line CK19 expression was positive. This newly established ORS cell line not only lays the foundation for further studying on the growth, regeneration, development law of goat hair follicle but also provides a mirror for the research of human hair in medical field.     

Comparative study on seasonal hair follicle cycling by analysis of the transcriptomes from cashmere and milk goats

Guard hair and cashmere undercoat are developed from primary and secondary hair follicle, respectively. Little is known about the gene expression differences between primary and secondary hair follicle cycling. In this study, we obtained RNA-seq data from cashmere and milk goats grown at four different seasons. We studied the differentially expressed genes (DEGs) during the yearly hair follicle cycling, and between cashmere and milk goats. WNT, NOTCH, MAPK, BMP, TGFβ and Hedgehog signaling pathways were involved in hair follicle cycling in both cashmere and milk goat. However, Milk goat DEGs between different months were significantly more than cashmere goat DEGs, with the largest difference being identified in December. Some expression dynamics were confirmed by quantitative PCR and western blot, and immunohistochemistry. This study offers new information sources related to hair follicle cycling in milk and cashmere goats, which could be applicable to improve the wool production and quality.     

Exposure to ram wool stimulates gonadotropin-releasing hormone pulse generator activity in the female goat

In the sheep and goat, exposure of anestrous females to a conspecific male odor enhances reproductive activity. Interestingly, a previous report indicated that male goat hair stimulated pulsatile luteinizing hormone (LH) secretion in the ewe. In the present study, we addressed whether ram wool affects the gonadotropin-releasing hormone (GnRH) pulse generator activity in the female goat. Five ovariectomized (OVX) goats were chronically implanted with recording electrodes in the mediobasal hypothalamus, and manifestations of the GnRH pulse generator were monitored as characteristic increases in multiple-unit activity (MUA volleys). Wool or hair samples were collected from a mature ram, ewe and male goat, and their effects on the MUA volley were examined. The exposure to ram wool induced an MUA volley within 1 min in all five OVX goats, as did the exposure to male goat hair. The ewe wool had no effect on the timing of an MUA volley occurrence. An invariable association of MUA volleys with LH pulses in the peripheral circulation was also confirmed in two OVX goats exposed to ram wool. The present results clearly indicate that exposure to ram wool stimulates pulsatile GnRH/LH release in the female goat. Since exposure to male goat hair enhances pulsatile LH secretion in the ewe, it is likely that very similar, if not identical, molecules are contained in the male-effect pheromone in the sheep and goat.     

Comparison of minisatellite polymorphisms in the Bacillus cereus complex: a simple assay for large-scale screening and identification of strains most closely related to Bacillus anthracis

Polymorphism of five tandem repeats that are monomorphic in Bacillus anthracis was investigated in 230 isolates of the B. cereus group and in 5 sequenced B. cereus genomes in search for markers allowing identification of B. cereus and B. thuringiensis strains most closely related to B. anthracis. Using this multiple-locus variable number of tandem repeat analysis (MLVA), a cluster of 30 strains was selected for further characterization. Eventually, six of these were characterized by multilocus sequence type analysis. One of the strains is only six point mutations (of almost 3,000 bp) away from B. anthracis and was also proposed to be closest to B. anthracis by MLVA analysis. However, this strain remains separated from B. anthracis by a number of significant genetic events observed in B. anthracis, including the loss of the hemolysin activity, the presence of four prophages, and the presence of the two virulence plasmids, pXO1 and pXO2. One particular minisatellite marker provides an efficient assay to identify the subset of B. cereus and B. thuringiensis strains closely related to B. anthracis. Based on these results, a very simple assay is proposed that allows the screening of hundreds of strains from the B. cereus complex, with modest equipment and at a low cost, to eventually fill the gap with B. anthracis and better understand the origin and making of this dangerous pathogen.     

Application of Variable-Number Tandem-Repeat Typing To Discriminate Ralstonia solanacearum Strains Associated with English Watercourses and Disease Outbreaks

Variable-number tandem-repeat (VNTR) analysis was used for high-resolution discrimination among Ralstonia solanacearum phylotype IIB sequevar 1 (PIIB-1) isolates and further evaluated for use in source tracing. Five tandem-repeat-containing loci (comprising six tandem repeats) discriminated 17 different VNTR profiles among 75 isolates from potato, geranium, bittersweet (Solanum dulcamara), tomato, and the environment. R. solanacearum isolates from crops at three unrelated outbreak sites where river water had been used for irrigation had distinct VNTR profiles that were shared with PIIB-1 isolates from infected bittersweet growing upriver of each site. The VNTR profiling results supported the implication that the source of R. solanacearum at each outbreak was contaminated river water. Analysis of 51 isolates from bittersweet growing in river water at different locations provided a means to evaluate the technique for studying the epidemiology of the pathogen in the environment. Ten different VNTR profiles were identified among bittersweet PIIB-1 isolates from the River Thames. Repeated findings of contiguous river stretches that produced isolates that shared single VNTR profiles supported the hypothesis that the pathogen had disseminated from infected bittersweet plants located upriver. VNTR profiles shared between bittersweet isolates from two widely separated Thames tributaries (River Ray and River Colne) suggested they were independently contaminated with the same clonal type. Some bittersweet isolates had VNTR profiles that were shared with potato isolates collected outside the United Kingdom. It was concluded that VNTR profiling could contribute to further understanding of R. solanacearum epidemiology and assist in control of future disease outbreaks.     

Differentiation of Bacillus anthracis, B. cereus, and B. thuringiensis by using pulsed-field gel electrophoresis

A pulsed-field gel electrophoresis (PFGE) method was developed for discriminating Bacillus anthracis from B. cereus and B. thuringiensis. A worldwide collection of 25 B. anthracis isolates showed high-profile homology, and these isolates were unambiguously distinguished from B. cereus and B. thuringiensis isolates by cluster analysis of the whole-genome macrorestriction enzyme digestion patterns generated by NotI.     

Global genetic population structure of Bacillus anthracis

Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064-6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.     

Distribution of S-layers on the surface of Bacillus cereus strains: phylogenetic origin and ecological pressure

Bacillus anthracis , Bacillus cereus and Bacillus thuringiensis have been described as members of the Bacillus cereus group but are, in fact, one species. B. anthracis is a mammal pathogen, B. thuringiensis an entomopathogen and B. cereus a ubiquitous soil bacterium and an occasional human pathogen. In two clinical isolates of B. cereus , in some B. thuringiensis strains and in B. anthracis , an S-layer has been described. We investigated how the S-layer is distributed in B. cereus , and whether phylogeny or ecology could explain its presence on the surface of some but not all strains. We first developed a simple biochemical assay to test for the presence of the S-layer. We then used the assay with 51 strains of known genetic relationship: 26 genetically diverse B. cereus and 25 non- B. anthracis of the B. anthracis cluster. When present, the genetic organization of the S-layer locus was analysed further. It was identical in B. cereus and B. anthracis . Nineteen strains harboured an S-layer, 16 of which belonged to the B. anthracis cluster. All 19 were B. cereus clinical isolates or B. thuringiensis , except for one soil and one dairy strain. These findings suggest a common phylogenetic origin for the S-layer at the surface of B. cereus strains and, presumably, ecological pressure on its maintenance.     

Identification of genes preferentially expressed in goat hair follicle anagen-catagen transition using suppression subtractive hybridization

Suppressive subtraction hybridization (SSH) was used to identify differentially expressed genes in goat (Capra hircus) hair follicle anagen-catagen transition. The cDNA fragments, derived from SSH positive subtractive library (tester: anagen-catagen transition, driver: later anagen), were cloned into pEGM-T vector. Two hundred cDNA fragments screened from this library were subjected to identify forty-five unregulated isolates. Sequence analysis revealed that these fragments represented twenty-three genes. Blasting analysis with database in GenBank showed that twenty genes were previously clearly annotated, two were homologous to un-annotated expressed sequence tag (ESTs), and one might be novel. To identify characters of gene expression, seven genes in later anagen and anagen-catagen transition skin tissues were chosen for quantitative real-time PCR. Results indicated that expression of these seven genes varied much, reaching threefold among them, furthering indicating that expression of those genes was up-regulation in the anagen-catagen transition. We characterized expression levels of this potential novel gene and the goat ectodysplasin A during differential stages of hair cycle. These profiles suggested that these two genes might play a role in the goat secondary hair follicle cycle.     

Meso-scale ecology of anthrax in southern Africa: a pilot study of diversity and clustering

It has only recently been possible to detect sufficient genetic diversity among anthrax isolates to allow genotype grouping (Keim et al. 1997). Early results of such grouping suggest that the southern African subcontinent may be the geographical origin of Bacillus anthracis. This report describes a pilot investigation of the genetic diversity of a study group of isolates from the Kruger National Park, South Africa, and efforts to detect spatio-temporal clustering within the study group. This study has also served as further validation for the newly developed Multi-Locus VNTR Analysis (MLVA), designed to simplify genotyping of B. anthracis isolates. The results reveal a diverse range of genotypes within the park allied with three genotype reference groups, and show that the MLVA procedure is robust for rapid analysis of B. anthracis genotypes. We also observed multiple genotype groups within epidemics and between geographically and temporally close epidemic episodes. This is in contrast to earlier characterizations of anthrax epidemics. The result of a Mantel test for time-space clustering indicates clustering of the anthrax isolates selected for the study.