A Genome-Wide Profiling Strategy as an Aid for Searching Unique Identification Biomarkers for <Emphasis Type="Italic">Streptococcus</Emphasis>

The use of rrs (16S rRNA) gene is widely regarded as the “gold standard” for identifying bacteria and determining their phylogenetic relationships. Nevertheless, multiple copies of this gene in a genome is likely to give an overestimation of the bacterial diversity. In each of the 50 Streptococcus genomes (16 species, 50 strains), 4–7 copies of rrs are present. The nucleotide sequences of these rrs genes show high similarity within and among genomes, which did not allow unambiguous identification. A genome-wide search revealed the presence of 27 gene sequences common to all the Streptococcus species. Digestion of these 27 gene sequences with 10 type II restriction endonucleases (REs) showed that unique RE digestion in purH gene is sufficient for clear cut identification of 30 genomes belonging to 16 species. Additional gene-RE combinations allowed identification of another 15 strains belonging to S. pneumoniae, S. pyogenes, and S. suis. For the rest 5 strains, a combination of 2 genes was required for identifying them. The proposed strategy is likely to prove helpful in proper detection of pathogens like Streptococcus.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0561-5) contains supplementary material, which is available to authorized users.     

Insights into the Origin of Clostridiumbotulinum Strains: Evolution of Distinct Restriction Endonuclease Sites in rrs (16S rRNA gene)

Diversity analysis of Clostridium botulinum strains is complicated by high microheterogeneity caused by the presence of 9–22 copies of rrs (16S rRNA gene). The need is to mine genetic markers to identify very closely related strains. Multiple alignments of the nucleotide sequences of the 212 rrs of 13 C. botulinum strains revealed intra- and inter-genomic heterogeneity. Low intragenomic heterogeneity in rrs was evident in strains 230613, Alaska E43, Okra, Eklund 17B, Langeland, 657, Kyoto, BKT015925, and Loch Maree. The most heterogenous rrs sequences were those of C. botulinum strains ATCC 19397, Hall, H04402065, and ATCC 3502. In silico restriction mapping of these rrs sequences was observable with 137 type II Restriction endonucleases (REs). Nucleotide changes (NC) at these RE sites resulted in appearance of distinct and additional sites, and loss in certain others. De novo appearances of RE sites due to NC were recorded at different positions in rrs gene. A nucleotide transition A>G in rrs of C. botulinum Loch Maree and 657 resulted in the generation of 4 and 10 distinct RE sites, respectively. Transitions A>G, G>A, and T>C led to the loss of RE sites. A perusal of the entire NC and in silico RE mapping of rrs of C. botulinum strains provided insights into their evolution. Segregation of strains on the basis of RE digestion patterns of rrs was validated by the cladistic analysis involving six house keeping genes: dnaN, gyrB, metG, prfA, pyrG, and Rho.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0514-z) contains supplementary material, which is available to authorized users.     

Searching Biomarkers in the Sequenced Genomes of <Emphasis Type="Italic">Staphylococcus</Emphasis> for their Rapid Identification

Bacterial identification using rrs (16S rRNA) gene is widely reported. Bacteria possessing multiple copies of rrs lead to overestimation of its diversity. Staphylococcus genomes carries 5–6 copies of rrs showing high similarity in their nucleotide sequences, which lead to ambiguous results. The genomes of 31 strains of Staphylococcus representing 7 species were searched for the presence of common genes. In silico digestion of 34 common genes using 10 restriction endonucleases (REs) lead to select gene-RE combinations, which could be used as biomarkers. RE digestion of recA allowed unambiguous identification of 13 genomes representing all the 7 species. In addition, a few more genes (argH, argR, cysS, gyrB, purH, and pyrE) and RE combinations permitted further identification of 12 strains. By employing additional RE and genes unique to a particular strain, it was possible to identify the rest 6 Staphylococcus aureus strains. This approach has the potential to be utilized for rapid detection of Staphylococcus strains.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-016-0565-9) contains supplementary material, which is available to authorized users.     

Genome Wide Search for Biomarkers to Diagnose Yersinia Infections

Bacterial identification on the basis of the highly conserved 16S rRNA (rrs) gene is limited by its presence in multiple copies and a very high level of similarity among them. The need is to look for other genes with unique characteristics to be used as biomarkers. Fifty-one sequenced genomes belonging to 10 different Yersinia species were used for searching genes common to all the genomes. Out of 304 common genes, 34 genes of sizes varying from 0.11 to 4.42 kb, were selected and subjected to in silico digestion with 10 different Restriction endonucleases (RE) (4–6 base cutters). Yersinia species have 6–7 copies of rrs per genome, which are difficult to distinguish by multiple sequence alignments or their RE digestion patterns. However, certain unique combinations of other common gene sequences—carB, fadJ, gluM, gltX, ileS, malE, nusA, ribD, and rlmL and their RE digestion patterns can be used as markers for identifying 21 strains belonging to 10 Yersinia species: Y. aldovae, Y. enterocolitica, Y. frederiksenii, Y. intermedia, Y. kristensenii, Y. pestis, Y. pseudotuberculosis, Y. rohdei, Y. ruckeri, and Y. similis. This approach can be applied for rapid diagnostic applications.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0552-6) contains supplementary material, which is available to authorized users.     

Genome Wide Analysis for Rapid Identification of Vibrio Species

The highly conserved 16S rRNA (rrs) gene is generally used for bacterial identification. In organisms possessing multiple copies of rrs, high intra-genomic heterogeneity does not allow easy distinction among different species. In order to identify Vibrio species, a wide range of genes have been employed. There is an urgent requirement of a consensus gene, which can be used as biomarker for rapid identification. Eight sequenced genomes of Vibrio species were screened for selecting genes which were common among all the genomes. Out of 108 common genes, 24 genes of sizes varying from 0.11 to 3.94 kb were subjected to in silico digestion with 10 type II restriction endonucleases (RE). A few unique genes—dapF, fadA, hisD, ilvH, lpxC, recF, recR, rph and ruvB in combination with certain REs provided unique digestion patterns, which can be used as biomarkers. This protocol can be exploited for rapid diagnosis of Vibrio species.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0553-5) contains supplementary material, which is available to authorized users.     

Genome Wide Analysis for Searching Novel Markers to Rapidly Identify Clostridium Strains

Microbial classification is based largely on the 16S rRNA (rrs) gene sequence, which is conserved throughout the prokaryotic domain. The Ribosomal Database Project (RDP) has become a reference point for almost all practical purposes. The use of this gene is limited by the fact that it can be used to identify only to the extent to what has been known and is available in the RDP. In order to identify an organism whose rrs is not present in the RDP database, we need to generate novel markers to place the unknown on the evolutionary map. Here, sequenced genomes of 27 Clostridium strains belonging to 9 species have been used to identify two sets of genes: (1) common to most of the species, and (2) unique to a species. Combinations of genes (recN, dnaJ, secA, mutS, and/or grpE) and their unique restriction endonuclease digestion (AluI, BfaI and/or Tru9I) patterns have been established to rapidly identify Clostridium species. This strategy for identifying novel markers can be extended to all other organisms and diagnostic applications.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0535-7) contains supplementary material, which is available to authorized users.     

Analysis of the unexplored features of rrs (16S rDNA) of the Genus Clostridium

Background

Bacterial taxonomy and phylogeny based on rrs (16S rDNA) sequencing is being vigorously pursued. In fact, it has been stated that novel biological findings are driven by comparison and integration of massive data sets. In spite of a large reservoir of rrs sequencing data of 1,237,963 entries, this analysis invariably needs supplementation with other genes. The need is to divide the genetic variability within a taxa or genus at their rrs phylogenetic boundaries and to discover those fundamental features, which will enable the bacteria to naturally fall within them. Within the large bacterial community, Clostridium represents a large genus of around 110 species of significant biotechnological and medical importance. Certain Clostridium strains produce some of the deadliest toxins, which cause heavy economic losses. We have targeted this genus because of its high genetic diversity, which does not allow accurate typing with the available molecular methods.

Results

Seven hundred sixty five rrs sequences (> 1200 nucleotides, nts) belonging to 110 Clostridium species were analyzed. On the basis of 404 rrs sequences belonging to 15 Clostridium species, we have developed species specific: (i) phylogenetic framework, (ii) signatures (30 nts) and (iii) in silico restriction enzyme (14 Type II REs) digestion patterns. These tools allowed: (i) species level identification of 95 Clostridium sp. which are presently classified up to genus level, (ii) identification of 84 novel Clostridium spp. and (iii) potential reduction in the number of Clostridium species represented by small populations.

Conclusions

This integrated approach is quite sensitive and can be easily extended as a molecular tool for diagnostic and taxonomic identification of any microbe of importance to food industries and health services. Since rapid and correct identification allows quicker diagnosis and consequently treatment as well, it is likely to lead to reduction in economic losses and mortality rates.     

New Pseudomonas spp. Are Pathogenic to Citrus

Five putative novel Pseudomonas species shown to be pathogenic to citrus have been characterized in a screening of 126 Pseudomonas strains isolated from diseased citrus leaves and stems in northern Iran. The 126 strains were studied using a polyphasic approach that included phenotypic characterizations and phylogenetic multilocus sequence analysis. The pathogenicity of these strains against 3 cultivars of citrus is demonstrated in greenhouse and field studies. The strains were initially grouped phenotypically and by their partial rpoD gene sequences into 11 coherent groups in the Pseudomonas fluorescens phylogenetic lineage. Fifty-three strains that are representatives of the 11 groups were selected and analyzed by partial sequencing of their 16S rRNA and gyrB genes. The individual and concatenated partial sequences of the three genes were used to construct the corresponding phylogenetic trees. The majority of the strains were identified at the species level: P. lurida (5 strains), P. monteilii (2 strains), P. moraviensis (1 strain), P. orientalis (16 strains), P. simiae (7 strains), P. syringae (46 strains, distributed phylogenetically in at least 5 pathovars), and P. viridiflava (2 strains). This is the first report of pathogenicity on citrus of P. orientalis, P. simiae, P. lurida, P. moraviensis and P. monteilii strains. The remaining 47 strains that could not be identified at the species level are considered representatives of at least 5 putative novel Pseudomonas species that are not yet described.     

Analyses of phylogeny, evolution, conserved sequences and genome-wide expression of the ICK/KRP family of plant CDK inhibitors

Background and Aims

The cell cycle is controlled by cyclin-dependent kinases (CDKs), and CDK inhibitors are major regulators of their activities. The ICK/KRP family of CDK inhibitors has been reported in several plants, with seven members in arabidopsis; however, the phylogenetic relationship among members in different species is unknown. Also, there is a need to understand how these genes and proteins are regulated. Furthermore, little information is available on the functional differences among ICK/KRP family members.

Methods

We searched publicly available databases and identified over 120 unique ICK/KRP protein sequences from more than 60 plant species. Phylogenetic analysis was performed using 101 full-length sequences from 40 species and intron–exon organization of ICK/KRP genes in model species. Conserved sequences and motifs were analysed using ICK/KRP protein sequences from arabidopsis (Arabidopsis thaliana), rice (Orysa sativa) and poplar (Populus trichocarpa). In addition, gene expression was examined using microarray data from arabidopsis, rice and poplar, and further analysed by RT-PCR for arabidopsis.

Key Results and Conclusions

Phylogenetic analysis showed that plant ICK/KRP proteins can be grouped into three major classes. Whereas the C-class contains sequences from dicotyledons, monocotyledons and gymnosperms, the A- and B-classes contain only sequences from dicotyledons or monocotyledons, respectively, suggesting that the A- and B-classes might have evolved from the C-class. This classification is also supported by exon–intron organization. Genes in the A- and B- classes have four exons, whereas genes in the C-class have only three exons. Analysis of sequences from arabidopsis, rice and poplar identified conserved sequence motifs, some of which had not been described previously, and putative functional sites. The presence of conserved motifs in different family members is consistent with the classification. In addition, gene expression analysis showed preferential expression of ICK/KRP genes in certain tissues. A model has been proposed for the evolution of this gene family in plants.     

Identification of signature and primers specific to genus <Emphasis Type="Italic">Pseudomonas</Emphasis> using mismatched patterns of 16S rDNA sequences

Background

Pseudomonas, a soil bacterium, has been observed as a dominant genus that survives in different habitats with wide hostile conditions. We had a basic assumption that the species level variation in 16S rDNA sequences of a bacterial genus is mainly due to substitutions rather than insertion or deletion of bases. Keeping this in view, the aim was to identify a region of 16S rDNA sequence and within that focus on substitution prone stretches indicating species level variation and to derive patterns from these stretches that are specific to the genus.

Results

Repeating elements that are highly conserved across different species of Pseudomonas were considered as guiding markers to locate a region within the 16S gene. Four repeating patterns showing more than 80% consistency across fifty different species of Pseudomonas were identified. The sub-sequences between the repeating patterns yielded a continuous region of 495 bases. The sub-sequences after alignment and using Shanon's entropy measure yielded a consensus pattern. A stretch of 24 base positions in this region, showing maximum variations across the sampled sequences was focused for possible genus specific patterns. Nine patterns in this stretch showed nearly 70% specificity to the target genus. These patterns were further used to obtain a signature that is highly specific to Pseudomonas. The signature region was used to design PCR primers, which yielded a PCR product of 150 bp whose specificity was validated through a sample experiment.

Conclusions

The developed approach was successfully applied to genus Pseudomonas. It could be tried in other bacterial genera to obtain respective signature patterns and thereby PCR primers, for their rapid tracking in the environmental samples.

     

MALDI-TOF mass spectrometry as a tool for differentiation of Bradyrhizobium species: Application to the identification of Lupinus nodulating strains

Genus Bradyrhizobium includes slow growing bacteria able to nodulate different legumes as well as species isolated from plant tumours. The slow growth presented by the members of this genus and the phylogenetic closeness of most of its species difficults their identification. In the present work we applied for the first time Matrix-Assisted Laser Desorption Ionization-Time-of-Flight Mass Spectrometry (MALDI-TOF MS) to the analysis of Bradyrhizobium species after the extension of MALDI Biotyper 2.0 database with the currently valid species of this genus. With this methodology it was possible to identify strains belonging to phylogenetically closely related species of genus Bradyrhizobium allowing the discrimination among species with rrs gene identities higher than 99%. The application of MALDI-TOF MS to strains isolated from nodules of different Lupinus species in diverse geographical locations allowed their correct identification when comparing with the results of rrs gene and ITS analyses. The nodulation of Lupinus gredensis, an endemic species of the west of Spain, by B. canariense supports the European origin of this species.     

The endemic Genista versicolor from Sierra Nevada National Park in Spain is nodulated by putative new Bradyrhizobium species and a novel symbiovar (sierranevadense)

Genista versicolor is an endemic legume from Sierra Nevada National Park which constitutes one of the UNESCO-recognized Biosphere Reserves. In the present study, a collection of strains nodulating this legume was analysed in characteristic soils of this ecosystem. Most strains nodulating G. versicolor belonged to rrs group I within the genus Bradyrhizobium and only one strain, named GV137, belonged to rrs group II from which only a single species, B. retamae, has been described in Europe to date. Strain GV137, and some strains from rrs group I, belonged to putative new species of Bradyrhizobium, although most strains from group I belonged to B. canariense, according to the ITS fragment and atpD gene analysis. This result contrasted with those obtained in Genista tinctoria in Northeast Europe whose endosymbionts were identified as B. japonicum. The analysis of the symbiotic nodC and nifH genes carried by G. versicolor-nodulating strains showed that most of them belonged to symbiovar genistearum, as did those isolated from G. tinctoria. Nevertheless, strain GV137, belonging to rrs group II, formed a divergent lineage that constituted a novel symbiovar within the genus Bradyrhizobium for which the name sierranevadense is proposed. This finding showed that the Genisteae are not restrictive legumes only nodulated by symbiovar genistearum, since Genista is a promiscuous legume nodulated by at least two symbiovars of Bradyrhizobium, as occurs in Retama species.     

Detection of Streptomycin Resistance in Mycobacterium tuberculosis Clinical Isolates Using Four Molecular Methods in China

To evaluate the relationship between mutations in rpsL or rrs genes and streptomycin (SM) resistance, we compared four molecular methods for their clinical value in the detection of SM resistance. Genotypic analysis of SM resistance in 167 M. tuberculosis clinical strains isolated from Chinese patients was performed by direct DNA sequencing, SSCP, RFLP, and reverse dot-blot hybridization (RDBH) assays. Of the 98 SM-resistant isolates, 78 (79.6%) had missense mutations in codon 43 or 88 of rpsL resulting in a Lys to Arg substitution, 6 (6.1%) had mutations of the rrs gene at positions 513 A to C or T or 516 C to T, and 14 (14.3%) had the wild-type sequence. None of the 69 SM-susceptible isolates examined had alterations in rpsL or rrs. The results of the SSCP, RFLP, and RDBH analyses for these mutations and wild-type sequences were completely consistent with DNA sequencing data. Five distinct single-nucleotide substitutions in codon 43 or 88 of rpsL gene or in position 513 or 516 of rrs gene were correctly identified in 84 of 98 (85.7%) phenotypically SM-resistant isolates by RDBH assay. Molecular analyses of the rpsL and rrs genes are useful for rapid prediction of SM resistance in most clinical strains of M. tuberculosis. Reverse dot-blot hybridization assay is a rapid, simple, and reliable method for the detection of drug resistance.     

Development and Application of Microsatellites in Candidate Genes Related to Wood Properties in the Chinese White Poplar (Populus tomentosa Carr.)

Gene-derived simple sequence repeats (genic SSRs), also known as functional markers, are often preferred over random genomic markers because they represent variation in gene coding and/or regulatory regions. We characterized 544 genic SSR loci derived from 138 candidate genes involved in wood formation, distributed throughout the genome of Populus tomentosa, a key ecological and cultivated wood production species. Of these SSRs, three-quarters were located in the promoter or intron regions, and dinucleotide (59.7%) and trinucleotide repeat motifs (26.5%) predominated. By screening 15 wild P. tomentosa ecotypes, we identified 188 polymorphic genic SSRs with 861 alleles, 2–7 alleles for each marker. Transferability analysis of 30 random genic SSRs, testing whether these SSRs work in 26 genotypes of five genus Populus sections (outgroup, Salix matsudana), showed that 72% of the SSRs could be amplified in Turanga and 100% could be amplified in Leuce. Based on genotyping of these 26 genotypes, a neighbour-joining analysis showed the expected six phylogenetic groupings. In silico analysis of SSR variation in 220 sequences that are homologous between P. tomentosa and Populus trichocarpa suggested that genic SSR variations between relatives were predominantly affected by repeat motif variations or flanking sequence mutations. Inheritance tests and single-marker associations demonstrated the power of genic SSRs in family-based linkage mapping and candidate gene-based association studies, as well as marker-assisted selection and comparative genomic studies of P. tomentosa and related species.     

Pseudomonas aestusnigri sp. nov., isolated from crude oil-contaminated intertidal sand samples after the Prestige oil spill

Strains VGXO14^T and Vi1 were isolated from the Atlantic intertidal shore from Galicia, Spain, after the Prestige oil spill. Both strains were Gram-negative rod-shaped bacteria with one polar inserted flagellum, strictly aerobic, and able to grow at 18–37 °C, pH 6–10 and 2–10% NaCl. A preliminary analysis of the 16S rRNA and the partial rpoD gene sequences indicated that these strains belonged to the Pseudomonas genus but were distinct from any known Pseudomonas species. A polyphasic taxonomic approach including phylogenetic, chemotaxonomic, phenotypic and genotypic data confirmed that the strains belonged to the Pseudomonas pertucinogena group. In a multilocus sequence analysis, the similarity of VGXO14^T and Vi1 to the closest type strain of the group, Pseudomonas pachastrellae, was 90.4%, which was lower than the threshold of 97% established to discriminate species in the Pseudomonas genus. The DNA–DNA hybridisation similarity between strains VGXO14^T and Vi1 was 79.6%, but below 70% with the type strains in the P. pertucinogena group. Therefore, the strains should be classified within the genus Pseudomonas as a novel species, for which the name Pseudomonas aestusnigri is proposed. The type strain is VGXO14^T (=CCUG 64165^T = CECT 8317^T).     

New species of pathogenic Pseudomonas isolated from citrus in Tunisia: Proposal of Pseudomonas kairouanensis sp. nov. and Pseudomonas nabeulensis sp. nov.

A collection of Pseudomonas strains was isolated in different regions of Tunisia in the period 2016–2017 from the fruits and leaves of Citrus sinensis cv. ‘Valencia Late’ and Citrus limon cv. ‘Eureka’ plants with symptoms of blast and black pit disease. A phylogenetic analysis of the housekeeping gene rpoD was used for strain identification at the species level. The results demonstrated the affiliation of these strains with the genus Pseudomonas and revealed the presence of 11 strains representing two putative new species in two monophyletic branches. These strains were analyzed morphologically and genotypically by multilocus sequence analyses of the rpoD, gyrB and 16S rRNA (rrs) gene sequences, and their phenotypic characteristics by API 20NE and Biolog GEN III. Plant pathogenic properties were confirmed on fruits and detached leaves of C. limon cv. ‘Eureka’. Fatty acids and WC MALDI-TOF MS major protein profiles were determined. The genomes of both representatives were sequenced. The average nucleotide index and genome-to-genome distance from KC12^T and E10B^T are below the cut-off established for a described species. These results support the conclusion that the strains KC12^T, KC17, KC20, KC22, KC24A, KC25 and KC26 represent a novel species of Pseudomonas, for which the name of Pseudomonas kairouanensis is proposed. The type strain is KC12^T (=CECT9766 and CFBP 8662). The strains E10B^T, E10AB, E10CB1 and Iy3BA represent another novel species of Pseudomonas for which the name of Pseudomonas nabeulensis is proposed; the type strain is E10B^T (=CECT9765 and CFBP 8661).     

Occurrence and genetic diversity of Arcobacter spp. in a spinach-processing plant and evaluation of two Arcobacter-specific quantitative PCR assays

Some species of the genus Arcobacter are considered to be emerging food pathogens. With respect to recent vegetable-borne outbreaks, the aim of this work was to investigate the occurrence and diversity of Arcobacter within the production chain of a spinach-processing plant by a combination of cultivation and molecular methods. Samples including spinach, water, and surface biofilm were taken over a period of three years from the entire processing line. Ten 16S rRNA (rrs) gene clone libraries were constructed and analysed using amplified rRNA gene restriction analysis (ARDRA). Approximately 1200 clones were studied that resulted in 44 operational taxonomic units (OTUs). Sequences with high similarities to Arcobacter cryaerophilus (13% of clones, 3 OTUs), A. ellisii (4%, 6 OTUs), A. suis (15%, 3 OTUs), and the type strain of A. nitrofigilis (1%, 7 OTUs) were identified. This represents the first report of the detection of the recently described species A. ellisii, A. suis and, in addition, A. venerupis from alternative habitats. A total of 67% of the clones (22 OTUs) could not be assigned to a genus, which indicated the presence of uncharacterised Arcobacter species.     

Molecular Identification of Diphyllobothrium nihonkaiense from 3 Human Cases in Heilongjiang Province with a Brief Literature Review in China

Human diphyllobothriasis is a widespread fish-borne zoonosis caused by the infection with broad tapeworms belonging to the genus Diphyllobothrium. In mainland China, so far 20 human cases of Diphyllobothrium infections have been reported, and the etiologic species were identified as D. latum and D. nihonkaiense based on morphological characteristics or molecular analysis. In the present study, proglottids of diphyllobothriid tapeworms from 3 human cases that occurred in Heilongjiang Province, China were identified as D. nihonkaiense by sequencing mitochondrial cytochrome c oxidase subunit I (cox1) and NADH dehydrogenase subunit 5 (nad5) genes. Two different cox1 gene sequences were obtained. One sequence showed 100% homology with those from humans in Japan. The remaining cox1 gene sequence and 2 different nad5 gene sequences obtained were not described previously, and might reflect endemic genetic characterizations. D. nihonkaiense might also be a major causative species of human diphyllobothriasis in China. Meanwhile, the finding of the first pediatric case of D. nihonkaiense infection in China suggests that infants infected with D. nihonkaiense should not be ignored.     

Seasonal and regional diversity of maple sap microbiota revealed using community PCR fingerprinting and 16S rRNA gene clone libraries

An arbitrary primed community PCR fingerprinting technique based on capillary electrophoresis was developed to study maple sap microbial community characteristics among 19 production sites in Québec over the tapping season. Presumptive fragment identification was made with corresponding fingerprint profiles of bacterial isolate cultures. Maple sap microbial communities were subsequently compared using a representative subset of 13 16S rRNA gene clone libraries followed by gene sequence analysis. Results from both methods indicated that all maple sap production sites and flow periods shared common microbiota members, but distinctive features also existed. Changes over the season in relative abundance of predominant populations showed evidence of a common pattern. Pseudomonas (64%) and Rahnella (8%) were the most abundantly and frequently represented genera of the 2239 sequences analyzed. Janthinobacterium, Leuconostoc, Lactococcus, Weissella, Epilithonimonas and Sphingomonas were revealed as occasional contaminants in maple sap. Maple sap microbiota showed a low level of deep diversity along with a high variation of similar 16S rRNA gene sequences within the Pseudomonas genus. Predominance of Pseudomonas is suggested as a typical feature of maple sap microbiota across geographical regions, production sites, and sap flow periods.