Characterization of serogroup C meningococci isolated from 14 provinces of China during 1966-2005 using comparative genomic hybridization

Neisseria meningitidis is a major cause of bacterial meningitis and septicemia worldwide. In China, serogroup A strains were responsible for over 95% of the cases, while serogroup B strains were mainly the cause of localized outbreaks and sporadic cases. Before 2003, serogroup C strains were only recovered from a few sporadic cases. However, a sudden increase in the number of cases due to serogroup C strains occurred during 2003-2005 in Anhui Province, China. Many cases were found in other provinces at the same time. Multilocus sequence typing (MLST) results indicated that the unique sequence type 4821 clone meningococci, a new hyper-virulent lineage, was responsible for the serogroup C meningitis outbreaks. We have completed the project of sequencing the whole genome of the Chinese N. meningitidis serogroup C representative isolate 053442. We fabricated a whole-genome microarray of N. meningitidis isolate 053442 and analyzed the genome composition differences among 81 serogroup C isolates which were isolated from 14 provinces of China during 1966-2005. The comparative genomic hybridization (CGH) result shows that the genome compositions of nearly all serogroup C isolates are similar to that of 053442. The products of many absent open reading frames (ORFs) are conserved hypothetical proteins. The results will provide a valuable resource from which one can analyze the genome composition and genetic background of serogroup C meningococci in China.     

Wavelet change-point prediction of transmembrane proteins

MOTIVATION: A non-parametric method, based on a wavelet data-dependent threshold technique for change-point analysis, is applied to predict location and topology of helices in transmembrane proteins. A new propensity scale generated from a transmembrane helix database is proposed. RESULTS: We show that wavelet change-point performs well for smoothing hydropathy and transmembrane profiles generated using different scales. We investigate which wavelet bases and threshold functions are overall most appropriate to detect transmembrane segments. Prediction accuracy is based on the analysis of two data sets used as standard benchmarks for transmembrane prediction algorithms. The analysis of a test set of 83 proteins results in accuracy per segment equal to 98.2%; the analysis of a 48 proteins blind-test set, i.e. containing proteins not used to generate the propensity scales, results in accuracy per segment equal to 97.4%. We believe that this method can also be applied to the detection of boundaries of other patterns such as G + Cisochores and dot-plots. AVAILABILITY: The transmembrane database, TMALN and source code are available upon request from the authors.     

Rapid laboratory detection of meningococcal disease outbreaks caused by serogroup C Neisseria meningitidis

Molecular subtyping is of significant importance to the recognition of outbreaks of meningococcal disease caused by serogroup C Neisseria meningitidis. We describe the application of multilocus variable number tandem repeat analysis (MLVA) for the molecular subtyping of N. meningitidis and compare its performance to that of pulsed-field gel electrophoresis (PFGE). For MLVA, a multiplex PCR assay targeting five variable number tandem repeat regions was developed and evaluated using a panel of sporadic and outbreak-associated serogroup C N. meningitidis isolates. MLVA was highly reproducible and provided results within 6 h. Overall, the discriminatory power of MLVA was equivalent to that of PFGE. The utilization of MLVA for subtyping N. meningitidis isolates provides a rapid and safer alternative to PFGE for identifying outbreaks of meningococcal disease. As such, it may provide public health officials with timely information that may minimize the spread of outbreak-related cases through prophylaxis.     

Detecting anomalous gene clusters and pathogenicity islands in diverse bacterial genomes.

A gene in a genome is defined as putative alien (pA) if its codon usage difference from the average gene exceeds a high threshold and codon usage differences from ribosomal protein genes, chaperone genes and protein-synthesis-processing factors are also high. pA gene clusters in bacterial genomes are relevant for detecting genomic islands (GIs), including pathogenicity islands (PAIs). Four other analyses appropriate to this task are G+C genome variation (the standard method); genomic signature divergences (dinucleotide bias); extremes of codon bias; and anomalies of amino acid usage. For example, the cagA domain of Helicobacter pylori is highly deviant in its genome signature and codon bias from the rest of the genome. Using these methods we can detect two potential PAIs in the Neisseria meningitidis genome, which contain hemagglutinin and/or hemolysin-related genes. Additionally, G+C variation and genome signature differences of the Mycobacterium tuberculosis genome indicate two pA gene clusters.     

A novel method to calculate the G+C content of genomic DNA sequences.

The base composition of a DNA fragment or genome is usually measured by the proportion of A+T or G+C in the sequence. The G+C content along genomic sequences is usually calculated using an overlapping or non-overlapping sliding window method. The result and accuracy of such an approach depends on the size of the window and the moving distance adopted. In this paper, a novel windowless technique to calculate the G+C content of genomic sequences is proposed. By this method, the G+C content can be calculated at different "resolution". In an extreme case, the G+C content may be computed at a specific point, rather than in a window of finite size. This is particularly useful to analyze the fine variation of base composition along genomic sequences. As the first example, the variation of G+C content along each of 16 yeast chromosomes is analyzed. The G+C-rich regions with length larger than 5 kb sequences are detected and listed in details. It is found that each chromosome consists of several G+C-rich and G+C-poor regions alternatively, i.e., a mosaic structure. Another example is to analyze the G+C content for each of the two chromosomes of the Vibrio cholerae genome. Based on the variations of the G+C content in each chromosome, it is shown that some fragments in the Vibrio cholerae genome may have been transferred from other species. Especially, the position and size of the large integron island on the smaller chromosome was precisely predicted. This method would be a useful tool for analyzing genomic sequences.     

Significance tests and weighted values for AFLP similarities, based on Arabidopsis in silico AFLP fragment length distributions

Many AFLP studies include relatively unrelated genotypes that contribute noise to data sets instead of signal. We developed: (1) estimates of expected AFLP similarities between unrelated genotypes, (2) significance tests for AFLP similarities, enabling the detection of unrelated genotypes, and (3) weighted similarity coefficients, including band position information. Detection of unrelated genotypes and use of weighted similarity coefficients will make the analysis of AFLP data sets more informative and more reliable. Test statistics and weighted coefficients were developed for total numbers of shared bands and for Dice, Jaccard, Nei and Li, and simple matching (dis)similarity coefficients. Theoretical and in silico AFLP fragment length distributions (FLDs) were examined as a basis for the tests. The in silico AFLP FLD based on the Arabidopsis thaliana genome sequence was the most appropriate for angiosperms. The G + C content of the selective nucleotides in the in silico AFLP procedure significantly influenced the FLD. Therefore, separate test statistics were calculated for AFLP procedures with high, average, and low G + C contents in the selective nucleotides. The test statistics are generally applicable for angiosperms with a G + C content of approximately 35-40%, but represent conservative estimates for genotypes with higher G + C contents. For the latter, test statistics based on a rice genome sequence are more appropriate.     

HmbR, a hemoglobin-binding outer membrane protein of Neisseria meningitidis, undergoes phase variation

Neisseria meningitidis uses hemoglobin (Hb) as an iron source via two TonB-dependent outer membrane receptors, HmbR and HpuB. Analysis of 25 epidemiologically unrelated clinical isolates from serogroups A, B, C, and Y revealed that 64% strains possessed both Hb receptor genes. Examination of the hmbR expression pattern in strains in which the hpuB gene was genetically inactivated revealed two distinct Hb utilization phenotypes. Five strains retained the ability to grow as a confluent lawn, while seven grew only as single colonies around Hb discs. The single-colony phenotype observed for some hpuB mutants is suggestive of phase variation of hmbR. The length of the poly(G) tract starting at position +1164 of hmbR absolutely correlated with the two Hb utilization phenotypes. All five strains that grew as confluent lawns around Hb discs possessed either 9 or 12 consecutive G residues. All seven strains that grew as single colonies around Hb discs had poly(G) tracts of a length other than 9 or 12. These single-colony variants that arose around the Hb discs had poly(G) tracts with either 9 or 12 consecutive G residues restoring the hmbR reading frame. Inactivation of hmbR in these strains resulted in a loss of Hb utilization, demonstrating that the change in the hmbR gene was responsible for the phenotypic switch. The switching rates from hmbR phase off to phase on were approximately 5 x 10(-4) in four serogroup C strains, 2 x 10(-2) in the serogroup A isolate, and 7 x 10(-6) in the serogroup B isolate.     

Whole genome sequencing to investigate the emergence of clonal complex 23 Neisseria meningitidis serogroup Y disease in the United States

In the United States, serogroup Y, ST-23 clonal complex Neisseria meningitidis was responsible for an increase in meningococcal disease incidence during the 1990s. This increase was accompanied by antigenic shift of three outer membrane proteins, with a decrease in the population that predominated in the early 1990s as a different population emerged later in that decade. To understand factors that may have been responsible for the emergence of serogroup Y disease, we used whole genome pyrosequencing to investigate genetic differences between isolates from early and late N. meningitidis populations, obtained from meningococcal disease cases in Maryland in the 1990s. The genomes of isolates from the early and late populations were highly similar, with 1231 of 1776 shared genes exhibiting 100% amino acid identity and an average π(N) = 0.0033 and average π(S) = 0.0216. However, differences were found in predicted proteins that affect pilin structure and antigen profile and in predicted proteins involved in iron acquisition and uptake. The observed changes are consistent with acquisition of new alleles through horizontal gene transfer. Changes in antigen profile due to the genetic differences found in this study likely allowed the late population to emerge due to escape from population immunity. These findings may predict which antigenic factors are important in the cyclic epidemiology of meningococcal disease.     

The longest (A+T) and (G+C) blocks in the human and other genomes

We analysed complete or almost complete nucleotide sequences of the human, chimp, mouse, rat, chicken, dog, and other genomes to find that they contain extremely long (A+T) a (G+C) blocks that do not occur at all in the corresponding randomized sequences. The longest is an (A+T) block containing 1040 consecutive AT pairs that occurs in the 16th human chromosome. The longest human (G+C) block has 261 bp in length. About a half of the longest blocks occur in introns. The (A+T) blocks are discrete units whereas the (G+C) blocks are diffuse. They are imbedded in the genome through connectors longer than 1 kilobase where the (G+C) content gradually decreases to the value of 50%. Remarkably, the (A+T) as well as (G+C) blocks are substantially shorter in the chimp genome. Chicken is characteristic by very long (G+C) blocks that are even longer than in the human genome. Though much shorter, long (G+C) and especially (A+T) blocks occur in lower organisms as well, which means that AT and GC pair clustering is an ancient property that has evolved into large scales in higher eukaryote genomes and the human genome in particular. Very long (A+T) and (G+C) blocks confer specific biophysical properties on DNA that are likely to influence genome folding in cell nuclei and its functional properties.     

Bordetella pertussis fimbriae are effective carrier proteins in Neisseria meningitidis serogroup C conjugate vaccines

Serogroup C meningococcal conjugate vaccines generally use diphtheria or tetanus toxoids as the protein carriers. The use of alternative carrier proteins may allow multivalent conjugate vaccines to be formulated into a single injection and circumvent potential problems of immune suppression in primed individuals. Bordetella pertussis fimbriae were assessed as carrier proteins for Neisseria meningitidis serogroup C polysaccharide. Fimbriae were conjugated to the polysaccharide using modifications of published methods and characterised by size exclusion chromatography; co-elution of protein and polysaccharide moieties confirmed conjugation. The conjugates elicited boostable IgG responses to fimbriae and serogroup C polysaccharide in mice, and IgG:IgM ratios indicated that the responses were thymus-dependent. High bactericidal antibody titres against a serogroup C strain of N. meningitidis were also observed. In a mouse infection model, the conjugate vaccine protected against lethal infection with N. meningitidis. Therefore, B. pertussis fimbriae are effective carrier proteins for meningococcal serogroup C polysaccharide and could produce a vaccine to protect against meningococcal disease and to augment protection against pertussis.     

Whole-genome sequence of the transformable Neisseria meningitidis serogroup A strain WUE2594

Serogroup A meningococci are a leading cause of bacterial meningitis in children and young adults worldwide. However, the genetic basis of serogroup A strains' virulence and their epidemiological properties remain poorly understood. Therefore, we sequenced the complete genome of the transformable Neisseria meningitidis serogroup A strain WUE2594.     

The Longest (A+T) and (G+C) Blocks in the Human and Other Genomes

Abstract

We analysed complete or almost complete nucleotide sequences of the human, chimp, mouse, rat, chicken, dog, and other genomes to find that they contain extremely long (A+T) a (G+C) blocks that do not occur at all in the corresponding randomized sequences. The longest is an (A+T) block containing 1040 consecutive AT pairs that occurs in the 16^th human chromosome. The longest human (G+C) block has 261 bp in length. About a half of the longest blocks occur in introns. The (A+T) blocks are discrete units whereas the (G+C) blocks are diffuse. They are embeeded in the genome through connectors longer than 1 kilobase where the (G+C) content gradually decreases to the value of 50%. Remarkably, the (A+T) as well as (G+C) blocks are substantially shorter in the chimp genome. Chicken is characteristic by very long (G+C) blocks that are even longer than in the human genome. Though much shorter, long (G+C) and especially (A+T) blocks occur in lower organisms as well, which means that AT and GC pair clustering is an ancient property that has evolved into large scales in higher eukaryote genomes and the human genome in particular. Very long (A+T) and (G+C) blocks confer specific biophysical properties on DNA that are likely to influence genome folding in cell nuclei and its functional properties.     

Conserved virulence of C to B capsule switched Neisseria meningitidis clinical isolates belonging to ET-37/ST-11 clonal complex

Capsule switching in Neisseria meningitidis is thought to occur by horizontal DNA exchange between meningococcal strains. Antigenic variants may be generated by allelic replacement of the siaD gene; the variants may then be selected by specific immunity against the capsular antigen. There were several vaccination campaigns against serogroup C in France in 2002, following an increase in the prevalence of invasive isolates of serogroup C of the phenotype C:2a:P1.5 and C:2a:P1.5,2 belonging to the ET-37/ST-11 clonal complex. We evaluated the emergence of capsule variants by the detection of B:2a:P1.5 and B:2a:P1.5,2 meningococcal isolates of the ET-37/ST-11 clonal complex. These isolates were significantly more frequent after the year 2002. Pulsed field gel electrophoresis profiles of the serogroup B (ET-37/ST-11) isolates differed from that of serogroup C (ET-37/ST-11) isolates by the bands that harbor the siaD genes responsible for the serogroup specificity. However, serogroup B and C, ET37/ST-11 isolates both express similar virulence as assessed from colonization and invasiveness in a mouse model. Our results indicate that capsule switching events within the same clonal complex can arise frequently with no alteration in virulence. This justifies an enhanced system of surveillance by molecular typing of such isolates, particularly after serogroup-specific vaccination.     

Application of Multiple-Locus Variable Number Tandem Repeat Analysis to Identify Outbreak-Associated Neisseria meningitides Serogroup C Sequence Type 4821 in China

Neisseria meningitidis (N. meningitidis) serogroup C sequence type (ST)-4821 caused an outbreak in 2010 in Shandong province of China. Twenty-one non-outbreak-associated strains were isolated, along with twenty-eight N. meningitides serogroup C ST-4821 isolates. Therefore, it’s essential to identify and clarify characterization of the real outbreak-associated strains with a rapid method during an outbreak investigation. In this study, multiple-locus variable number tandem repeat analysis (MLVA) was applied to analyze 84 N. meningitidis strains, among which 58 were recovered from two outbreaks and 26 were sporadic isolates. Three MLVA schemes with different combination of VNTR loci were tested, and two of them were suitable for isolates from China: scheme 2 with six loci was found to separate ST into finer resolution, and scheme 3 with five loci can be used to identify outbreak-associated isolates from the same outbreak that caused by N. meningitidis serogroup C ST-4821.     

HmbR outer membrane receptors of pathogenic Neisseria spp.: iron-regulated, hemoglobin-binding proteins with a high level of primary structure conservation.

We have recently cloned and characterized the hemoglobin receptor gene from Neisseria meningitidis serogroup C. N. meningitidis cells expressing HmbR protein were able to bind biotinylated hemoglobin, and the binding was specifically inhibited by unlabeled hemoglobin and not heme. The HmbR-mediated hemoglobin binding activity of N. meningitidis cells was shown to be iron regulated. The presence of hemoglobin but not heme in the growth medium stimulated HmbR-mediated hemoglobin binding activity. The efficiency of utilization of different hemoglobins by the HmbR-expressing N. meningitidis cells was shown to be species specific; human hemoglobin was the best source of iron, followed by horse, rat, turkey, dog, mouse, and sheep hemoglobins, The phenotypic characterization of HmbR mutants of some clinical strains of N. meningitidis suggested the existence of two unrelated hemoglobin receptors. The HmbR-unrelated hemoglobin receptor was shown to be identical to Hpu, the hemoglobin-haptoglobin receptor of N. meningitidis. The Hpu-dependent hemoglobin utilization system was not able to distinguish between different sources of hemoglobin; all animal hemoglobins were utilized equally well. HmbR-like genes are also present in N. meningitidis serogroups A and B, Neisseria gonorrhoeae MS11 and FA19, Neisseria perflava, and Neisseria polysaccharea. The hemoglobin receptor genes from N. meningitidis serogroups A and B and N. gonorrhoeae MS11 were cloned, and their nucleotide sequences were determined. The nucleotide sequence identity ranged between 86.5% (for N. meningitidis serogroup B hmbR and MS11 hmbR) and 93.4% (for N. meningitidis serogroup B hmbR and N. meningitidis serogroup C hmbR). The deduced amino acid sequences of these neisserial hemoglobin receptors were also highly related, with overall 84.7% conserved amino acid residues. A stop codon was found in the hmbR gene of N. gonorrhoeae MS11. This strain was still able to use hemoglobin and hemoglobin-haptoglobin complexes as iron sources, indicating that some gonococci may express only the HmbR-independent hemoglobin utilization system.     

Isochore structures in the chicken genome

The availability of the complete chicken genome sequence provides an unprecedented opportunity to study the global genome organization at the sequence level. Delineating compositionally homogeneous G + C domains in DNA sequences can provide much insight into the understanding of the organization and biological functions of the chicken genome. A new segmentation algorithm, which is simple and fast, has been proposed to partition a given genome or DNA sequence into compositionally distinct domains. By applying the new segmentation algorithm to the draft chicken genome sequence, the mosaic organization of the chicken genome can be confirmed at the sequence level. It is shown herein that the chicken genome is also characterized by a mosaic structure of isochores, long DNA segments that are fairly homogeneous in the G + C content. Consequently, 25 isochores longer than 2 Mb (megabases) have been identified in the chicken genome. These isochores have a fairly homogeneous G + C content and often correspond to meaningful biological units. With the aid of the technique of cumulative GC profile, we proposed an intuitive picture to display the distribution of segmentation points. The relationships between G + C content and the distributions of genes (CpG islands, and other genomic elements) were analyzed in a perceivable manner. The cumulative GC profile, equipped with the new segmentation algorithm, would be an appropriate starting point for analyzing the isochore structures of higher eukaryotic genomes.     

In vivo expression of Neisseria meningitidis proteins homologous to the Haemophilus influenzae Hap and Hia autotransporters

The genome sequences of Neisseria meningitidis serogroup B strain MC58 and serogroup A strain Z2491 were systematically searched for open reading frames (ORFs) encoding autotransporters. Eight ORFs were identified, six of which were present in both genomes, whereas two were specific for MC58. Among the identified ORFs was the gene encoding the known autotransporter IgA1 protease. The deduced amino acid sequences of the other identified ORFs were homologous to known autotransporters and found to contain an N-terminal signal sequence and a C-terminal domain that could constitute a beta-barrel in the outer membrane. The ORFs NMB1985 and NMB0992, encoding homologs of the Hap (for Haemophilus adhesion and penetration protein) and Hia (for Haemophilus influenzae adherence protein) autotransporters of H. influenzae, were cloned from serogroup B strain H44/76 and expressed in Escherichia coli. Western blots revealed that all sera of patients (n=14) and healthy carriers (n=3) tested contained antibodies against at least one of the recombinant proteins. These results indicate that both genes are widely distributed among N. meningitidis isolates and expressed during colonization and infection.     

The establishment of Neisseria meningitidis serogroup W135 of the clonal complex ET-37/ST-11 as an epidemic clone and the persistence of serogroup A isolates in Burkina Faso

We analyzed 48 invasive isolates of Neisseria meningitidis that were isolated from meningitis cases in Burkina Faso (April 2002 to April 2003). Thirty-nine of these isolates had the phenotype (serogroup:serotype:serosubtype) W135:2a:P1.5,2, eight isolates were A:4:P1.9 and one isolate was nongroupable:nonserotypable:nonserosubtypable. Genotyping of meningococcal isolates showed that W135 isolates belonged to the sequence type (ST)-11. The nongroupable isolate was of genogroup W135 and belonged to ST-192. Isolates of serogroup A belonged to ST-2859 (a member of the subgroup III/ST-5 clonal complex). W135 (ST-11) isolates involved in meningitis outbreaks in Burkina Faso differed from those involved in the Hajj-2000 associated outbreak by their pulsed-field gel electrophoresis profile. These data confirm the changing epidemiology of meningococcal infection in Burkina Faso with the establishment and expansion of serogroup W135 N. meningitidis strains of the ET-37/ST-11 clonal complex, as well as the emergence of a new clone within the subgroup III/ST-5 clonal complex.     

The Neisseria meningitidis serogroup A capsular polysaccharide O-3 and O-4 acetyltransferase

Neisseria meningitidis serogroup A capsular polysaccharide (CPS) is composed of a homopolymer of O-acetylated, alpha1-->6-linked ManNAc 1-phosphate that is distinct from the capsule structures of the other meningococcal disease-causing serogroups, B, C, Y, and W-135. The serogroup A capsule biosynthetic genetic cassette consists of four open reading frames, mynA-D (sacA-D), that are specific to serogroup A, but the functions of these genes have not been well characterized. mynC was found to encode an inner membrane-associated acetyltransferase that is responsible for the O-acetylation of the CPS of serogroup A. The wild-type CPS as revealed by 1H NMR had 60-70% O-acetylated ManNAc residues that contained acetyl groups at O-3, with some species acetylated at O-4 and at both O-3 and O-4. A non-polar mynC mutant generated by introducing an aphA-3 kanamycin resistance cassette produced CPS with no O-acetylation. A serogroup A capsule-specific monoclonal antibody was shown to recognize the wild-type O-acetylated CPS, but not the CPS of the mynC mutant, which lacked O-acetylation. MynC was C-terminally His-tagged and overexpressed in Escherichia coli to obtain the predicted approximately 26-kDa protein. The acetyltransferase activity of purified MynC was demonstrated in vitro using [14C]acetyl-CoA. MynC O-acetylated the O-acetylated CPS of the mynC mutant and further acetylated the wild-type CPS of serogroup A meningococci, but not the CPS of serogroup B or C meningococci. Genetic complementation of the mynC mutant confirmed the function of MynC as the serogroup A CPS O-3 and O-4 acetyltransferase. MynC represents a new subclass of O-acetyltransferases that utilize acetyl-CoA to decorate the D-mannosamine capsule of N. meningitidis serogroup A.     

Searching for RNA genes using base-composition statistics

The hypothesis that genomic regions rich in non-protein-coding RNAs (ncRNAs) can be identified using local variations in single-base and dinucleotide statistics has been investigated. (G+C)%, (G-C)% difference, (A-T)% difference and dinucleotide-frequency statistics were compared among seven classes of ncRNAs and three genomes. Significant variations were observed in (G+C)% and, in Methanococcus jannaschii, in the frequency of the dinucleotide 'CG'. Screening programs based on these two base-composition statistics were developed. With (G+C)% screening alone, a 1% fraction of the M.jannaschii genome containing all 44 known transfer RNAs, ribosomal RNAs and signal recognition particle RNAs could be identified. When (G+C)% combined with CG dinucleotide-frequency screening was used, 43 of the 44 known M.jannaschii structural ncRNAs were again identified, while the number of presumably false hits overlapping a known or putative protein-coding gene was reduced from 15 to 6. In addition, 19 candidate ncRNAs were identified including one with significant homology to several known archaeal RNaseP RNAs.