The emergence and transmission dynamics of HIV-1 CRF07_BC in Mainland China.

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset ('full1155'), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences ('backbone_295') and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989-August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

Fast Evaluation of Viral Emerging Risks (FEVER): A computational tool for biosurveillance, diagnostics, and mutation typing of emerging viral pathogens.

Viral pathogens can rapidly evolve, adapt to novel hosts, and evade human immunity. The early detection of emerging viral pathogens through biosurveillance coupled with rapid and accurate diagnostics are required to mitigate global pandemics. However, RNA viruses can mutate rapidly, hampering biosurveillance and diagnostic efforts. Here, we present a novel computational approach called FEVER (Fast Evaluation of Viral Emerging Risks) to design assays that simultaneously accomplish: 1) broad-coverage biosurveillance of an entire group of viruses, 2) accurate diagnosis of an outbreak strain, and 3) mutation typing to detect variants of public health importance. We demonstrate the application of FEVER to generate assays to simultaneously 1) detect sarbecoviruses for biosurveillance; 2) diagnose infections specifically caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2); and 3) perform rapid mutation typing of the D614G SARS-CoV-2 spike variant associated with increased pathogen transmissibility. These FEVER assays had a high in silico recall (predicted positive) up to 99.7% of 525,708 SARS-CoV-2 sequences analyzed and displayed sensitivities and specificities as high as 92.4% and 100% respectively when validated in 100 clinical samples. The D614G SARS-CoV-2 spike mutation PCR test was able to identify the single nucleotide identity at position 23,403 in the viral genome of 96.6% SARS-CoV-2 positive samples without the need for sequencing. This study demonstrates the utility of FEVER to design assays for biosurveillance, diagnostics, and mutation typing to rapidly detect, track, and mitigate future outbreaks and pandemics caused by emerging viruses.

Optical Biosensor Platforms Display Varying Sensitivity for the Direct Detection of Influenza RNA.

Detection methods that do not require nucleic acid amplification are advantageous for viral diagnostics due to their rapid results. These platforms could provide information for both accurate diagnoses and pandemic surveillance. Influenza virus is prone to pandemic-inducing genetic mutations, so there is a need to apply these detection platforms to influenza diagnostics. Here, we analyzed the Fast Evaluation of Viral Emerging Risks (FEVER) pipeline on ultrasensitive detection platforms, including a waveguide-based optical biosensor and a flow cytometry bead-based assay. The pipeline was also evaluated in silico for sequence coverage in comparison to the U.S. Centers for Disease Control and Prevention's (CDC) influenza A and B diagnostic assays. The influenza FEVER probe design had a higher tolerance for mismatched bases than the CDC's probes, and the FEVER probes altogether had a higher detection rate for influenza isolate sequences from GenBank. When formatted for use as molecular beacons, the FEVER probes detected influenza RNA as low as 50 nM on the waveguide-based optical biosensor and 1 nM on the flow cytometer. In addition to molecular beacons, which have an inherently high background signal we also developed an exonuclease selection method that could detect 500 pM of RNA. The combination of high-coverage probes developed using the FEVER pipeline coupled with ultrasensitive optical biosensors is a promising approach for future influenza diagnostic and biosurveillance applications.

Prevalence of gp160 polymorphisms known to be related to decreased susceptibility to temsavir in different subtypes of HIV-1 in the Los Alamos National Laboratory HIV Sequence Database.

BACKGROUND: Fostemsavir, a prodrug of the gp120-directed attachment inhibitor temsavir, is indicated for use in heavily treatment-experienced individuals with MDR HIV-1. Reduced susceptibility to temsavir in the clinic maps to discrete changes at amino acid positions in gp160: S375, M426, M434 and M475. OBJECTIVES: To query the Los Alamos National Laboratory (LANL) HIV Sequence Database for the prevalence of polymorphisms at gp160 positions of interest. METHODS: Full-length gp160 sequences (N = 7560) were queried for amino acid polymorphisms relative to the subtype B consensus at positions of interest; frequencies were reported for all sequences and among subtypes/circulating recombinant forms (CRFs) with ≥10 isolates in the database. RESULTS: Among 239 subtypes in the database, the 5 most prevalent were B (n = 2651, 35.1%), C (n = 1626, 21.5%), CRF01_AE (n = 674, 8.9%), A1 (n = 273, 3.6%) and CRF02_AG (n = 199, 2.6%). Among all 7560 sequences, the most prevalent amino acids at positions of interest (S375, 73.5%; M426, 82.1%; M434, 88.2%; M475, 89.9%) were the same as the subtype B consensus. Specific polymorphisms with the potential to decrease temsavir susceptibility (S375H/I/M/N/T/Y, M426L/P, M434I/K and M475I) were found in <10% of isolates of subtypes D, G, A6, BC, F1, CRF07_BC, CRF08_BC, 02A, CRF06_cpx, F2, 02G and 02B. S375H and M475I were predominant among CRF01_AE (S375H, 99.3%; M475I, 76.3%; consistent with previously reported low temsavir susceptibility of this CRF) and 01B (S375H, 71.7%; M475I, 49.5%). CONCLUSIONS: Analysis of the LANL HIV Sequence Database found a low prevalence of gp160 amino acid polymorphisms with the potential to reduce temsavir susceptibility overall and among most of the common subtypes.

HIV-1 and SARS-CoV-2: Patterns in the evolution of two pandemic pathogens.

Humanity is currently facing the challenge of two devastating pandemics caused by two very different RNA viruses: HIV-1, which has been with us for decades, and SARS-CoV-2, which has swept the world in the course of a single year. The same evolutionary strategies that drive HIV-1 evolution are at play in SARS-CoV-2. Single nucleotide mutations, multi-base insertions and deletions, recombination, and variation in surface glycans all generate the variability that, guided by natural selection, enables both HIV-1's extraordinary diversity and SARS-CoV-2's slower pace of mutation accumulation. Even though SARS-CoV-2 diversity is more limited, recently emergent SARS-CoV-2 variants carry Spike mutations that have important phenotypic consequences in terms of both antibody resistance and enhanced infectivity. We review and compare how these mutational patterns manifest in these two distinct viruses to provide the variability that fuels their evolution by natural selection.

Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years.

We aimed to investigate whether the sequence length of HIV-1 increases over time. We performed a longitudinal analysis of full-length coding region sequences (FLs) during an HIV-1 outbreak among patients with hemophilia and local controls infected with the Korean subclade B of HIV-1 (KSB). Genes were amplified by overlapping RT-PCR or nested PCR and subjected to direct sequencing. Overall, 141 FLs were sequentially determined over 30 years in 62 KSB-infected patients. Phylogenetic analysis indicated that within KSB, two FLs from plasma donors O and P comprised two clusters, together with 8 and 12 patients with hemophilia, respectively. Signature pattern analysis of the KSB of HIV-1 revealed 91 signature nucleotide residues (1.1%). In total, 48 and 43 signature nucleotides originated from clusters O and P, respectively. Six positions contained 100% specific nucleotide(s) in clusters O and P. In-depth FL analysis for over 30 years indicated that the KSB FL significantly increased over time before combination antiretroviral therapy (cART) and decreased with cART. This increase occurred due to the significant increase in env and nef genes, originating in the variable regions of both genes. The increase in sequence length of HIV-1 over time suggests an evolutionary direction.

Evolutionary history, potential intermediate animal host, and cross-species analyses of SARS-CoV-2.

To investigate the evolutionary history of the recent outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in China, a total of 70 genomes of virus strains from China and elsewhere with sampling dates between 24 December 2019 and 3 February 2020 were analyzed. To explore the potential intermediate animal host of the SARS-CoV-2 virus, we reanalyzed virome data sets from pangolins and representative SARS-related coronaviruses isolates from bats, with particular attention paid to the spike glycoprotein gene. We performed phylogenetic, split network, transmission network, likelihood-mapping, and comparative analyses of the genomes. Based on Bayesian time-scaled phylogenetic analysis using the tip-dating method, we estimated the time to the most recent common ancestor and evolutionary rate of SARS-CoV-2, which ranged from 22 to 24 November 2019 and 1.19 to 1.31 × 10-3 substitutions per site per year, respectively. Our results also revealed that the BetaCoV/bat/Yunnan/RaTG13/2013 virus was more similar to the SARS-CoV-2 virus than the coronavirus obtained from the two pangolin samples (SRR10168377 and SRR10168378). We also identified a unique peptide (PRRA) insertion in the human SARS-CoV-2 virus, which may be involved in the proteolytic cleavage of the spike protein by cellular proteases, and thus could impact host range and transmissibility. Interestingly, the coronavirus carried by pangolins did not have the RRAR motif. Therefore, we concluded that the human SARS-CoV-2 virus, which is responsible for the recent outbreak of COVID-19, did not come directly from pangolins.

Detection of Human Immunodeficiency Virus Type 1 (HIV-1) Antisense Protein (ASP) RNA Transcripts in Patients by Strand-Specific RT-PCR.

In retroviruses, antisense transcription has been described in both human immunodeficiency virus type 1 (HIV-1) and human T-lymphotropic virus 1 (HTLV-1). In HIV-1, the antisense protein ASP gene is located on the negative strand of env, in the reading frame -2, spanning the junction gp120/gp41. In the sense orientation, the 3' end of the ASP open reading frame overlaps with gp120 hypervariable regions V4 and V5. The study of ASP RNA has been thwarted by a phenomenon known as RT-self-priming, whereby RNA secondary structures have the ability to prime RT in absence of the specific primer, generating non-specific cDNAs. The combined use of high RNA denaturation with biotinylated reverse primers in the RT reaction, together with affinity purification of the cDNA onto streptavidin-coated magnetic beads, has allowed us to selectively amplify ASP RNA in CD4+ T cells derived from individuals infected with HIV-1. Our method is relatively low-cost, simple to perform, highly reliable, and easily reproducible. In this respect, it represents a powerful tool for the study of antisense transcription not only in HIV-1 but also in other biological systems.

Genetic Analysis of the Full-Length gag Gene from the Earliest Korean Subclade B of HIV-1: An Outbreak among Korean Hemophiliacs.

We determined the earliest full-length HIV-1 gag gene sequences in 110 patients with HIV-1, including 20 hemophiliacs (HPs) and 90 local controls (LCs). The gag gene from stored sera was amplified using RT-PCR, and was subjected to direct sequencing. Phylogenetic analysis indicated that 94 and 16 sequences belonged to the Korean subclade of HIV-1 subtype B (KSB) and subtype B, respectively. A total of 12 signature pattern amino acids were found within the KSB, distinct from the worldwide consensus of subtype B. Within the KSB, the gag gene sequences from donors O and P and those from the 20 HPs comprised two subclusters. In particular, sequences from donor O strongly clustered with those of eight HPs. Moreover, signature pattern analysis indicated that 14 signature nucleotides were shared between the HPs and LCs within KSB (p < 0.01). Among the 14 nucleotides, positions 9 and 5 belonged to clusters O and P, respectively. In conclusion, signature pattern analysis for the gag gene revealed 12 signature pattern residues within the KSB and also confirmed the previous conclusion that the 20 HPs were infected with viruses due to incompletely inactivated clotting factor IX. This study is the first genetic analysis of the HIV-1 gag gene in Korea.

Impact of HIV-1 subtype and Korean Red Ginseng on AIDS progression: comparison of subtype B and subtype D.

BACKGROUND: To date, no study has described disease progression in Asian patients infected with HIV-1 subtype D. METHODS: To determine whether the disease progression differs in patients infected with subtypes D and B prior to starting combination antiretroviral therapy, the annual decline (AD) in CD4+ T cell counts over 96 ± 59 months was retrospectively analyzed in 163 patients and compared in subtypes D and B based on the nef gene. RESULTS: CD4+ T cell AD was significantly higher in the six subtype D-infected patients than in the 157 subtype B-infected patients irrespective of Korean Red Ginseng (KRG) treatment (p < 0.001). Of these, two subtype D-infected patients and 116 subtype B-infected patients had taken KRG. AD was significantly lower in patient in the KRG-treated group than in those in the KRG-naïve group irrespective of subtype (p < 0.05). To control for the effect of KRG, patients not treated with KRG were analyzed, with AD found to be significantly greater in subtype D-infected patients than in subtype B-infected patients (p < 0.01). KRG treatment had a greater effect on AD in subtype D-infected patients than in subtype B-infected patients (4.5-fold vs. 1.6-fold). Mortality rates were significantly higher in both the 45 KRG-naïve (p < 0.001) and all 163 (p < 0.01) patients infected with subtype D than subtype B. CONCLUSION: Subtype D infection is associated with a >2-fold higher risk of death and a 2.9-fold greater rate of progression than subtype B, regardless of KRG treatment.

Detection of antisense protein (ASP) RNA transcripts in individuals infected with human immunodeficiency virus type 1 (HIV-1).

The detection of antisense RNA is hampered by reverse transcription (RT) non-specific priming, due to the ability of RNA secondary structures to prime RT in the absence of specific primers. The detection of antisense RNA by conventional RT-PCR does not allow assessment of the polarity of the initial RNA template, causing the amplification of non-specific cDNAs. In this study we have developed a modified protocol for the detection of human immunodeficiency virus type 1 (HIV-1) antisense protein (ASP) RNA. Using this approach, we have identified ASP transcripts in CD4+ T cells isolated from five HIV-infected individuals, either untreated or under suppressive therapy. We show that ASP RNA can be detected in stimulated CD4+ T cells from both groups of patients, but not in unstimulated cells. We also show that in untreated patients, the patterns of expression of ASP and env are very similar, with the levels of ASP RNA being markedly lower than those of env. Treatment of cells from one viraemic patient with α-amanitin greatly reduces the rate of ASP RNA synthesis, suggesting that it is associated with RNA polymerase II, the central enzyme in the transcription of protein-coding genes. Our data represent the first nucleotide sequences obtained in patients for ASP, demonstrating that its transcription indeed occurs in those HIV-1 lineages in which the ASP open reading frame is present.

Persistent Viral Reservoirs in Lymphoid Tissues in SIV-Infected Rhesus Macaques of Chinese-Origin on Suppressive Antiretroviral Therapy.

Understanding HIV latent reservoirs in tissues is essential for the development of new strategies targeting these sites for eradication. Here, we assessed the size of latent reservoirs and the source of residual viruses in multiple lymphoid tissues of SIV-infected and fully suppressed rhesus macaques of Chinese-origin (cRMs). Eight cRMs were infected with SIVmac251 and treated with tenofovir and emtricitabine daily for 24 weeks initiated 4 weeks post-infection. Four of the eight animals reached sustained full viral suppression with undetectable viremia. The levels of cell-associated SIV DNA varied in peripheral blood mononuclear cells (PBMCs) and multiple lymphoid tissues, but with higher levels in the mesenteric lymph nodes (MesLNs). The levels of cell-associated SIV RNA also varied in different tissues. The higher frequency of viral RNA detection in the MesLNs was also observed by in situ hybridization. Consistently, the infection unit per million cells (IUPM) in the MesLNs was higher than in PBMCs and other tested lymphoid tissues by quantitative viral outgrowth assay (QVOA). Furthermore, env gp120 from tissue SIV RNA was amplified by single genome amplification. Phylogenetic analysis revealed diverse variants from tissues parallel to the viral inoculum in all viral suppressed animals. These results demonstrate that the latency and viral reservoirs in the lymphoid tissues still exist in aviremic macaques under full suppressive therapy. Moreover, the size of viral latent reservoirs differs in various lymphoid tissues with a relatively larger size in the MesLNs.

HIV-2 Infection in a Migrant from Gambia: The History of the Disease Combined with Phylogenetic Analysis Revealed the Real Source of Infection.

Human immunodeficiency virus type 2 (HIV-2) infection prevalence is increasing in some European countries. The increasing migratory flow from countries where HIV-2 is endemic has facilitated the spread of the virus into Europe and other regions. We describe a case of HIV-2 infection in a migrant individual in the Asylum Seekers Centre (ASC) in Italy. The patient's virus was sequenced and found to be a typical HIV-2 genotype A virus. Bayesian evolutionary analysis revealed that the HIV-2 sequence from migrant dated back to 1986 in a subcluster, including sequences from Guinea Bissau. This was coherent with the history of the migrant who lived in Guinea Bissau from his birth until 1998 when he was 13 years old. Monitoring for HIV-2 infection in migrants from western Africa is necessary using adequate molecular tools to improve the diagnosis and understand the real origin of infection.

High Prevalence of Non-B HIV-1 Subtypes in Overseas Sailors and Prostitutes in Korea.

There have been no studies related to groups at the highest risk for HIV-1 infection in Korea before 1993. In this study, for the first time, we report the distribution of HIV subtypes in overseas sailors (OSs) and prostitutes who worked in brothels near U.S. military bases in Korea. We retrospectively determined the sequences of nef in 131 patients using reverse transcription polymerase chain reaction (RT-PCR). These patients composed of 102 OSs, 14 OS spouses, and 15 prostitutes. Phylogenetic analysis was performed using 128 Korean OSs, OS spouses, and prostitutes. The distribution of non-B subtypes (n = 105) was as follows: 39, CRF02_AG; 15, CRF01_AE; 7, A1; 7, A2; 6, D; 2, CRF06_cpx; 3, C; 6, G; 11, untypable; and 1 each for CRF09_cpx, CRF12_BF, CRF50_A1D, A3, AFG, H, F1, F2, and A. Of the 116 OSs and OS spouses, 101 (87%), 11 (9%), and 4 (3%) subjects had non-B, Western B, and Korean subclade B (KSB) HIV-1s, respectively. Among the 15 prostitutes, 10 had Western B (67%), 4 non-B (27%), and 1 KSB (7%) HIV-1s. All 14 couples, each comprising of an OS and his spouse, had the same subtype. KSB (5%) was detected in OSs and prostitutes in 1990 and 1994, respectively. Of the 131 patients analyzed in this study, 105 (80%), 21 (16%), and 5 (4%) were infected with the non-B, Western B, and KSB subtypes of HIV, respectively. In future, these data may provide an important foundation for analysis of HIV-1 subtypes in Korea.

Persistence of SIV in the brain of SIV-infected Chinese rhesus macaques with or without antiretroviral therapy.

Persistence of HIV-1 reservoirs in the central nervous system (CNS) is an obstacle to cure strategies. However, little is known about residual viral distribution, viral replication levels, and genetic diversity in different brain regions of HIV-infected individuals on combination antiretroviral therapy (cART). Because myeloid cells particularly microglia are likely major reservoirs in the brain, and more microglia exist in white matter than gray matter in a human brain, we hypothesized the major viral reservoirs in the brain are the white matter reflected by higher levels of viral DNA. To address the issue, we used the Chinese rhesus macaque (ChRM) model of SIV infection, and treated 11 SIVmac251-infected animals including long-term nonprogressors with cART for up to 24 weeks. SIV reservoirs were assessed by SIV DNA levels in 16 specific regions of the brain and 4 regions of spinal cord. We found relatively high frequencies of SIV in basal ganglia and brain stem compared to other regions. cART-receiving animals had significantly lower SIV DNA levels in the gray matter than white matter. Moreover, a shortened envelope gp120 with 21 nucleotide deletions and guanine-to-adenine hypermutations were observed. These results demonstrate that SIV enters the CNS in SIV-infected ChRM with a major reservoir in the white matter after cART; the SIV/ChRM/cART is an appropriate model for studying HIV CNS reservoirs and testing new eradication strategies. Further, examining multiple regions of the CNS may be needed when assessing whether an agent is successful in reducing the size of SIV reservoirs in the CNS.

Signature pattern analysis for the full-length env gene of the earliest Korean subclade B of HIV-1: outbreak among Korean hemophiliacs.

The epidemiological link in the hypervariable env gene between viruses infecting HIV-positive hemophiliacs (HPs) and plasma donors was not studied. We determined full-length env gene sequences in 20 HPs, 3 plasma donors whose plasma was used for domestic clotting factor (DCF) production, and 54 local controls (LCs). Env genes from viruses in frozen stored sera obtained 1-3 years after diagnosis and from samples collected several years after infection were amplified via RT-PCR and subjected to direct sequencing. Phylogenetic analysis indicated that all sequences were subtype B, including 133 sequences from 77 cases (20 HPs, 3 plasma donors, and 54 LCs) belonging to the Korean subclade B (KSB) and 6 sequences from 5 cases that did not belong to the KSB. Env gene sequences from donors O and P and those of the 20 HPs comprised 2 subclusters within the KSB, although phylogenetic analysis did not support significant bootstrap values. In contrast, signature pattern analysis indicated signature nucleotides at 43 positions between the HPs and LCs (P < 0.05). In particular, specific signature nucleotides at 4 positions were fully conserved in the HPs, but not in the LCs (P < 0.0001). Furthermore, there were 26 signature residues within the KSB and were distinct from the worldwide consensus for subtype B. In conclusion, signature pattern analysis for the hypervariable env gene revealed an epidemiological link that the 20 HPs in this study had been infected with viruses from the DCF used for treatment, consistent with our previous finding.

Primate immunodeficiency virus classification and nomenclature: Review.

The International Committee for the Taxonomy and Nomenclature of Viruses does not rule on virus classifications below the species level. The definition of species for viruses cannot be clearly defined for all types of viruses. The complex and interesting epidemiology of Human Immunodeficiency Viruses demands a detailed and informative nomenclature system, while at the same time it presents challenges such that many of the rules need to be flexibly applied or modified over time. This review outlines the nomenclature system for primate lentiviruses and provides an update on new findings since the last review was written in 2000.

Comparative genomic analyses reveal broad diversity in botulinum-toxin-producing Clostridia.

BACKGROUND: Clostridium botulinum is a diverse group of bacteria characterized by the production of botulinum neurotoxin. Botulinum neurotoxins are classified into serotypes (BoNT/A-G), which are produced by six species/Groups of Clostridia, but the genetic background of the bacteria remains poorly understood. The purpose of this study was to use comparative genomics to provide insights into the genetic diversity and evolutionary history of bacteria that produce the potent botulinum neurotoxin. RESULTS: Comparative genomic analyses of over 170 Clostridia genomes, including our draft genome assemblies for 59 newly sequenced Clostridia strains from six continents and publicly available genomic data, provided in-depth insights into the diversity and distribution of BoNT-producing bacteria. These newly sequenced strains included Group I and II strains that express BoNT/A,/B,/E, or/F as well as bivalent strains. BoNT-producing Clostridia and closely related Clostridia species were delineated with a variety of methods including 16S rRNA gene, concatenated marker genes, core genome and concatenated multi-locus sequencing typing (MLST) gene phylogenies that related whole genome sequenced strains to publicly available strains and sequence types. These analyses illustrated the phylogenetic diversity in each Group and the diversity of genomic backgrounds that express the same toxin type or subtype. Comparisons of the botulinum neurotoxin genes did not identify novel toxin types or variants. CONCLUSIONS: This study represents one of the most comprehensive analyses of whole genome sequence data for Group I and II BoNT-producing strains. Read data and draft genome assemblies generated for 59 isolates will be a resource to the research community. Core genome phylogenies proved to be a powerful tool for differentiating BoNT-producing strains and can provide a framework for the study of these bacteria. Comparative genomic analyses of Clostridia species illustrate the diversity of botulinum-neurotoxin-producing strains and the plasticity of the genomic backgrounds in which bont genes are found.

Genomic sequences of six botulinum neurotoxin-producing strains representing three clostridial species illustrate the mobility and diversity of botulinum neurotoxin genes.

The whole genomes for six botulinum neurotoxin-producing clostridial strains were sequenced to provide references for under-represented toxin types, bivalent strains or unusual toxin complexes associated with a bont gene. The strains include three Clostridium botulinum Group I strains (CDC 297, CDC 1436, and Prevot 594), a Group II C. botulinum strain (Eklund 202F), a Group IV Clostridium argentinense strain (CDC 2741), and a Group V Clostridium baratii strain (Sullivan). Comparisons of the Group I genomic sequences revealed close relationships and conservation of toxin gene locations with previously published Group I C. botulinum genomes. The bont/F6 gene of strain Eklund 202F was determined to be a chimeric toxin gene composed of bont/F1 and bont/F2. The serotype G strain CDC 2741 remained unfinished in 20 contigs with the bont/G located within a 1.15Mb contig, indicating a possible chromosomal location for this toxin gene. Within the genome of C. baratii Sullivan strain, direct repeats of IS1182 insertion sequence (IS) elements were identified flanking the bont/F7 toxin complex that may be the mechanism of bont insertion into C. baratii. Highlights of the six strains are described and release of their genomic sequences will allow further study of unusual neurotoxin-producing clostridial strains.