Comparative genomic analysis of bacteriocin-producing Weissella cibaria 110.

Weissella cibaria 110 was isolated from plaa-som, a Thai fermented fish product, and known to produce the weissellicin 110 bacteriocin. We carried out comprehensive comparative genomic analysis of W. cibaria 110 with four other non-bacteriocin-producing W. cibaria strains and identified potential antibiotic-resistant genes. We further identified a type III restriction-modification system, a TA system, and a bacteriocin gene cluster that are unique in W. cibaria 110. Genes related to bacteriocin biosynthesis are organized in clusters and are encoded with minimum genetic machinery consisting of structural cognate immunity genes, including ABC transporter and immunity protein. Finally, we predicted W. cibaria 110 to produce a class IId bacteriocin, weissellicin 110, which is 31 amino acids in length and contains a 21-amino-acid N-terminal leader peptide. This is the first bacteriocin-producing sequencing genome in W. cibaria, and we describe the difference between the bacteriocin-producing and non bacteriocin-producing strains from genome point of view.

Unique clustering genes in the bacterial chromosome affecting the type-III secretion of enterohaemorrhagic Escherichia coli.

Bioinformatics analysis was used to search for unknown genes that might influence the phenotypic presentations of enterohaemorrhagic Escherichia coli (EHEC). By so doing and using the known genomic data from EHEC O157 : H7 and K-12, it has been deduced that genes Z4863 to Z4866 of EHEC do not exist in K-12 strains. These four gene sequences have low degrees of homology (18-40 % amino acid identities) to a set of genes in K-12, which have been known to encode fatty acid biosynthesis enzymes. We referred these four consecutive genes as a fasyn cluster and found that deletion of fasyn from EHEC resulted in a defective type-III secretion (T3S). This deletion apparently did not decrease the amounts of the T3S proteins ectopically expressed from plasmids. Examination of the corresponding mRNAs by real-time PCR revealed that the mRNAs readily decreased in the fasyn-deleted mutant and this suppressive effect on the mRNA levels appeared to spread across all lee operons. Complementation with fasyn reverted the T3S-deficient phenotype. Furthermore, this reversion was also seen when the mutant was supplemented with locus of enterocyte effacement activators (Ler or GrlA). Thus, these unique clustering genes located apart from locus of enterocyte effacement on the bacterial chromosome also play a role in affecting T3S of EHEC.

Protein domain repetition is enriched in Streptococcal cell-surface proteins.

Tandem repetition of domain in protein sequence occurs in all three domains of life. It creates protein diversity and adds functional complexity in organisms. In this work, we analyzed 52 streptococcal genomes and found 3748 proteins contained domain repeats. Proteins not harboring domain repeats are significantly enriched in cytoplasm, whereas proteins with domain repeats are significantly enriched in cytoplasmic membrane, cell wall and extracellular locations. Domain repetition occurs most frequently in S. pneumoniae and least in S. thermophilus and S. pyogenes. DUF1542 is the highest repeated domain in a single protein, followed by Rib, CW_binding_1, G5 and HemolysinCabind. 3D structures of 24 repeat-containing proteins were predicted to investigate the structural and functional effect of domain repetition. Several repeat-containing streptococcal cell surface proteins are known to be virulence-associated. Surface-associated tandem domain-containing proteins without experimental functional characterization may be potentially involved in the pathogenesis of streptococci and deserve further investigation.

Highly conserved influenza A virus epitope sequences as candidates of H3N2 flu vaccine targets.

This study focused on identifying the conserved epitopes in a single subtype A (H3N2)-as candidates for vaccine targets. We identified a total of 32 conserved epitopes in four viral proteins [22 HA, 4PB1, 3 NA, 3 NP]. Evaluation of conserved epitopes in coverage during 1968-2010 revealed that (1) 12 HA conserved epitopes were highly present in the circulating viruses; (2) the remaining 10 HA conserved epitopes appeared with lower percentage but a significantly increasing trend after 1989 [p<0.001]; and (3) the conserved epitopes in NA, NP and PB1 are also highly frequent in wild-type viruses. These conserved epitopes also covered an extremely high percentage of the 16 vaccine strains during the 42 year period. The identification of highly conserved epitopes using our approach can also be applied to develop broad-spectrum vaccines.

T4-Like genome organization of the Escherichia coli O157:H7 lytic phage AR1.

We report the genome organization and analysis of the first completely sequenced T4-like phage, AR1, of Escherichia coli O157:H7. Unlike most of the other sequenced phages of O157:H7, which belong to the temperate Podoviridae and Siphoviridae families, AR1 is a T4-like phage known to efficiently infect this pathogenic bacterial strain. The 167,435-bp AR1 genome is currently the largest among all the sequenced E. coli O157:H7 phages. It carries a total of 281 potential open reading frames (ORFs) and 10 putative tRNA genes. Of these, 126 predicted proteins could be classified into six viral orthologous group categories, with at least 18 proteins of the structural protein category having been detected by tandem mass spectrometry. Comparative genomic analysis of AR1 and four other completely sequenced T4-like genomes (RB32, RB69, T4, and JS98) indicated that they share a well-organized and highly conserved core genome, particularly in the regions encoding DNA replication and virion structural proteins. The major diverse features between these phages include the modules of distal tail fibers and the types and numbers of internal proteins, tRNA genes, and mobile elements. Codon usage analysis suggested that the presence of AR1-encoded tRNAs may be relevant to the codon usage of structural proteins. Furthermore, protein sequence analysis of AR1 gp37, a potential receptor binding protein, indicated that eight residues in the C terminus are unique to O157:H7 T4-like phages AR1 and PP01. These residues are known to be located in the T4 receptor recognition domain, and they may contribute to specificity for adsorption to the O157:H7 strain.

A Deep Learning Approach for Missing Data Imputation of Rating Scales Assessing Attention-Deficit Hyperactivity Disorder.

A variety of tools and methods have been used to measure behavioral symptoms of attention-deficit/hyperactivity disorder (ADHD). Missing data is a major concern in ADHD behavioral studies. This study used a deep learning method to impute missing data in ADHD rating scales and evaluated the ability of the imputed dataset (i.e., the imputed data replacing the original missing values) to distinguish youths with ADHD from youths without ADHD. The data were collected from 1220 youths, 799 of whom had an ADHD diagnosis, and 421 were typically developing (TD) youths without ADHD, recruited in Northern Taiwan. Participants were assessed using the Conners' Continuous Performance Test, the Chinese versions of the Conners' rating scale-revised: short form for parent and teacher reports, and the Swanson, Nolan, and Pelham, version IV scale for parent and teacher reports. We used deep learning, with information from the original complete dataset (referred to as the reference dataset), to perform missing data imputation and generate an imputation order according to the imputed accuracy of each question. We evaluated the effectiveness of imputation using support vector machine to classify the ADHD and TD groups in the imputed dataset. The imputed dataset can classify ADHD vs. TD up to 89% accuracy, which did not differ from the classification accuracy (89%) using the reference dataset. Most of the behaviors related to oppositional behaviors rated by teachers and hyperactivity/impulsivity rated by both parents and teachers showed high discriminatory accuracy to distinguish ADHD from non-ADHD. Our findings support a deep learning solution for missing data imputation without introducing bias to the data.

Characterization of a new plasmid-like prophage in a pandemic Vibrio parahaemolyticus O3:K6 strain.

Vibrio parahaemolyticus is a common food-borne pathogen that is normally associated with seafood. In 1996, a pandemic O3:K6 strain abruptly appeared and caused the first pandemic of this pathogen to spread throughout many Asian countries, America, Europe, and Africa. The role of temperate bacteriophages in the evolution of this pathogen is of great interest. In this work, a new temperate phage, VP882, from a pandemic O3:K6 strain of V. parahaemolyticus was purified and characterized after mitomycin C induction. VP882 was a Myoviridae bacteriophage with a polyhedral head and a long rigid tail with a sheath-like structure. It infected and lysed high proportions of V. parahaemolyticus, Vibrio vulnificus, and Vibrio cholerae strains. The genome of phage VP882 was sequenced and was 38,197 bp long, and 71 putative open reading frames were identified, of which 27 were putative functional phage or bacterial genes. VP882 had a linear plasmid-like genome with a putative protelomerase gene and cohesive ends. The genome does not integrate into the host chromosome but was maintained as a plasmid in the lysogen. Analysis of the reaction sites of the protelomerases in different plasmid-like phages revealed that VP882 and PhiHAP-1 were highly similar, while N15, PhiKO2, and PY54 made up another closely related group. The presence of DNA adenine methylase and quorum-sensing transcriptional regulators in VP882 may play a specific role in this phage or regulate physiological or virulence-associated traits of the hosts. These genes may also be remnants from the bacterial chromosome following transduction.

A common virulence plasmid in biotype 2 Vibrio vulnificus and its dissemination aided by a conjugal plasmid.

Strains of Vibrio vulnificus, a marine bacterial species pathogenic for humans and eels, are divided into three biotypes, and those virulent for eels are classified as biotype 2. All biotype 2 strains possess one or more plasmids, which have been shown to harbor the biotype 2-specific DNA sequences. In this study we determined the DNA sequences of three biotype 2 plasmids: pR99 (68.4 kbp) in strain CECT4999 and pC4602-1 (56.6 kb) and pC4602-2 (66.9 kb) in strain CECT4602. Plasmid pC4602-2 showed 92% sequence identity with pR99. Curing of pR99 from strain CECT4999 resulted in loss of resistance to eel serum and virulence for eels but had no effect on the virulence for mice, an animal model, and resistance to human serum. Plasmids pC4602-2 and pR99 could be transferred to the plasmid-cured strain by conjugation in the presence of pC4602-1, which was self-transmissible, and acquisition of pC4602-2 restored the virulence of the cured strain for eels. Therefore, both pR99 and pC4602-2 were virulence plasmids for eels but not mice. A gene in pR99, which encoded a novel protein and had an equivalent in pC4602-2, was further shown to be essential, but not sufficient, for the resistance to eel serum and virulence for eels. There was evidence showing that pC4602-2 may form a cointegrate with pC4602-1. An investigation of six other biotype 2 strains for the presence of various plasmid markers revealed that they all harbored the virulence plasmid and four of them possessed the conjugal plasmid in addition.

Bacterial Composition and Diversity in Breast Milk Samples from Mothers Living in Taiwan and Mainland China.

Human breast milk is widely recognized as the best source of nutrients for healthy growth and development of infants; it contains a diverse microbiota. Here, we characterized the diversity of the microbiota in the breast milk of East Asian women and assessed whether delivery mode influenced the microbiota in the milk of healthy breast-feeding mothers. We profiled the microbiota in breast milk samples collected from 133 healthy mothers in Taiwan and in six regions of mainland China (Central, East, North, Northeast, South, and Southwest China) by using 16S rRNA pyrosequencing. Lactation stage (months postpartum when the milk sample was collected) and maternal body mass index did not influence the breast milk microbiota. Bacterial composition at the family level differed significantly among samples from the seven geographical regions. The five most predominant bacterial families were Streptococcaceae (mean relative abundance: 24.4%), Pseudomonadaceae (14.0%), Staphylococcaceae (12.2%), Lactobacillaceae (6.2%), and Oxalobacteraceae (4.8%). The microbial profiles were classified into three clusters, driven by Staphylococcaceae (abundance in Cluster 1: 42.1%), Streptococcaceae (Cluster 2: 48.5%), or Pseudomonadaceae (Cluster 3: 26.5%). Microbial network analysis at the genus level revealed that the abundances of the Gram-positive Staphylococcus, Streptococcus, and Rothia were negatively correlated with those of the Gram-negative Acinetobacter, Bacteroides, Halomonas, Herbaspirillum, and Pseudomonas. Milk from mothers who had undergone Caesarian section (C-section group) had a significantly higher abundance of Lactobacillus (P < 0.05) and a higher number of unique unclassified operational taxonomic units (OTUs) (P < 0.001) than that from mothers who had undergone vaginal delivery (vaginal group). These findings revealed that (i) geographic differences in the microbial profiles were found in breast milk from mothers living in Taiwan and mainland China, (ii) the predominant bacterial families Streptococcaceae, Staphylococcaceae, and Pseudomonadaceae were key components for forming three respective clusters, and (iii) a significantly greater number of unique OTUs was found in the breast milk from mothers who had undergone C-section than from those who had delivered vaginally.

Si-Wu-Tang and its constituents promote mammary duct cell proliferation by up-regulation of HER-2 signaling.

OBJECTIVE: The consumption of over-the-counter natural products by perimenopausal women remains a challenging problem. It is our aim to investigate the proliferative effect of Si-Wu-Tang (SWT) and its constituents on MCF7 breast cancer cells as well as the SWT-modulated cell signaling and HER-2 gene expression. DESIGN: By using the MCF7 (ER+, HER-2 low), BT474 (ER+, HER-2 high), MDAMB231 (ER-, HER-2 low), and SKBR3 (ER-, HER-2 high) mammary duct cell lines as our in vitro model, the mitogenic effects of SWT and its constituents were assessed by trypan blue dye exclusion assay and DNA flow cytometry. SWT-modulated cell signaling and HER-2 gene expression were evaluated in the MCF7 line by Western blot and reverse transcriptase-polymerase chain reaction analyses. RESULTS: The results showed that SWT and some of its constituents dose-dependently stimulated cell proliferation of MCF7 cells. The activation of HER-2, its downstream signaling molecules AKT and ERK1/2, as well as HER-2 gene up-regulation were involved in SWT-stimulated cell proliferation. The addition of neutralizing antibody against HER-2 abrogated the SWT-up-regulated HER-2 expression, indicating a positive feedback control for the action of HER-2 in this setting. Ferulic acid, one of the major compounds in SWT, not only promoted cell proliferation of MCF7, BT474, MDAMB231, and SKBR3 cells, but also increased the phosphorylation of HER-2, AKT, and ERK1/2, as well as overexpression of HER-2, on MCF7 cells. CONCLUSIONS: We conclude that SWT and its active constituents stimulate mammary duct cell proliferation by modulating HER-2, PI3K/AKT, and MAPK signaling and the positive feedback of HER-2 gene expression. This provides important information for perimenopausal women who are at risk of or have breast cancer or other growths in breast tissue.

Lessons from the Largest Epidemic of Avian Influenza Viruses in Taiwan, 2015.

The largest epidemic of avian influenza (AI) in history attacked poultry and wild birds throughout Taiwan starting January 6, 2015. This study analyzed surveillance results, epidemiologic characteristics, and viral sequences by using government-released information, with the intention to provide recommendations to minimize future pandemic influenza. The H5 clade 2.3.4.4 highly pathogenic AI viruses (HPAIVs) had not been detected in Taiwan before 2015. During this epidemic, four types of etiologic agents were identified: the three novel subtypes H5N2, H5N8, and H5N3 clade 2.3.4.4 HPAIVs and one endemic chicken H5N2 subtype (Mexican-like lineage) of low pathogenic AI viruses. Cocirculation of mixed subtypes also occurred, with H5N2 clade 2.3.4.4 HPAIVs accompanied by the H5N8 and H5N3 subtypes or old H5N2 viruses in the same farm. More than 90% of domestic geese died from this AI epidemic; geese were affected the most at the early outbreaks. The epidemic peaked in mid-January for all three novel H5 subtypes. Spatial epidemiology found that most affected areas were located in southwestern coastal areas. In terrestrial poultry (mostly chickens), different geographic distributions of AI virus subtypes were detected, with hot spots of H5N2 clade 2.3.4.4 vs. past-endemic old H5N2 viruses in Changhwa (P = 0.03) and Yunlin (P = 0.007) counties, respectively, of central Taiwan. Phylogenetic and sequence analyses of all the early 10 Taiwan H5 clade 2.3.4.4 isolates covering the three subtypes showed that they were very different from the HA of the past local H5 viruses from domestic ducks (75%-80%) and chickens (70%-75%). However, they had the highest sequence identity percentages (99.53%-100%), with the HA of A/crane/Kagoshima/KU13/2014(H5N8) isolated on December 7, 2014, in Japan being higher than those of recent American and Korean H5 HPAIVs [A/Northern pintail/Washington/40964/2014 (H5N2) and A/gyrfalcon/Washington/41088-6/2014 (H5N8): 99.02%-99.54% and A/Baikal teal/Korea/Donglim3/2014 (H5N8): 98.61%-99.08%], implying a likely common ancestor of these H5 clade 2.3.4.4 viruses. The multiple subtypes of H5 clade 2.3.4.4 HPAIVs imply high viral reassortment. We recommend establishing an integrated surveillance system, involving clinical, virologic, and serologic surveillance in poultry and wild birds, swine and other mammals prevalent on multiple-animal mixed-type traditional farms, and high-risk human populations, as a crucially important step to minimize future pandemic influenza.

ErbB2 activation of ESX gene expression.

Various members of the Ets multigene family exhibit diverse roles in development, cell differentiation, tissue-specific gene expression and human malignancy. In the search for Ets factors involved in mammary gland development and malignancy, ESX was found to be upregulated in a subset of breast tumours and cell lines. We report the transcriptional regulation of ESX in epithelial breast cancer cells. Transient reporter assays using the ESX promoter show that ESX transcription is regulated by ErbB receptor signalling. In cell lines and in 45 primary ductal breast cancers we show that ESX transcript expression significantly correlates with ErbB2 transcript levels. Moreover, expression of ErbB2 in cells upregulates ESX promoter activity while inhibition of ErbB2 or its downstream signaling pathways decrease both ESX promoter activity and endogenous ESX protein levels. These results indicate that the ESX promoter represents a transcriptional target of ErbB2, and ESX expression may represent a downstream mediator of ErbB2 signaling and ErbB2-induced gene expression.

Exploring comprehensive within-motif dependence of transcription factor binding in Escherichia coli.

Modeling the binding of transcription factors helps to decipher the control logic behind transcriptional regulatory networks. Position weight matrix is commonly used to describe a binding motif but assumes statistical independence between positions. Although current approaches take within-motif dependence into account for better predictive performance, these models usually rely on prior knowledge and incorporate simple positional dependence to describe binding motifs. The inability to take complex within-motif dependence into account may result in an incomplete representation of binding motifs. In this work, we applied association rule mining techniques and constructed models to explore within-motif dependence for transcription factors in Escherichia coli. Our models can reflect transcription factor-DNA recognition where the explored dependence correlates with the binding specificity. We also propose a graphical representation of the explored within-motif dependence to illustrate the final binding configurations. Understanding the binding configurations also enables us to fine-tune or design transcription factor binding sites, and we attempt to present the configurations through exploring within-motif dependence.

A17, the First Sequenced Strain of Lactococcus lactis subsp. cremoris with Potential Immunomodulatory Functions.

Lactococcus lactis subsp. cremoris A17, isolated from Taiwan fermented cabbage, is the first sequenced strain of L. lactis subsp. cremoris with immunomodulatory activity and antiallergic functions. The resulting A17 draft genome contains 2,679,936 bp and indicates that A17 is a potential exopolysaccharide-producing strain without any known virulence gene.

Behavior training reverses asymmetry in hippocampal transcriptome of the cav3.2 knockout mice.

Homozygous Cav3.2 knockout mice, which are defective in the pore-forming subunit of a low voltage activated T-type calcium channel, have been documented to show impaired maintenance of late-phase long-term potentiation (L-LTP) and defective retrieval of context-associated fear memory. To investigate the role of Cav3.2 in global gene expression, we performed a microarray transcriptome study on the hippocampi of the Cav3.2-/- mice and their wild-type littermates, either naïve (untrained) or trace fear conditioned. We found a significant left-right asymmetric effect on the hippocampal transcriptome caused by the Cav3.2 knockout. Between the naive Cav3.2-/- and the naive wild-type mice, 3522 differentially expressed genes (DEGs) were found in the left hippocampus, but only 4 DEGs were found in the right hippocampus. Remarkably, the effect of Cav3.2 knockout was partially reversed by trace fear conditioning. The number of DEGs in the left hippocampus was reduced to 6 in the Cav3.2 knockout mice after trace fear conditioning, compared with the wild-type naïve mice. To our knowledge, these results demonstrate for the first time the asymmetric effects of the Cav3.2 and its partial reversal by behavior training on the hippocampal transcriptome.

Genome architecture of Lactobacillus plantarum PS128, a probiotic strain with potential immunomodulatory activity.

BACKGROUND: Clinical and preclinical observations indicate that Lactobacillus plantarum has anti-inflammatory activity and may regulate the immune responses of its hosts when ingested. Recently, modification of teichoic acids (TAs) produced by L. plantarum was reported as a key to regulating the systemic immune response in mice. However, data linking TA-related genetic determinants and the immunomodulatory effect are limited. To provide genomic information for elucidating the underlying mechanism of immunomodulation by L. plantarum, we sequenced the genome of L. plantarum strain PS128. RESULTS: The PS128 genome contains 11 contigs (3,325,806 bp; 44.42% GC content) after hybrid assembly of sequences derived with Illumina MiSeq and PacBio RSII systems. The most abundant functions of the protein-coding genes are carbohydrate, amino acid, and protein metabolism. The 16S rDNA sequences of PS128 are closest to the sequences of L. plantarum WCFS1 and B21; these three strains form a distinct clade based on 16S rDNA sequences. PS128 shares core genes encoding the metabolism, transport, and modification of TAs with other sequenced L. plantarum strains. Compared with the TA-related genes of other completely sequenced L. plantarum strains, the PS128 contains more lipoteichoic acid exporter genes. CONCLUSIONS: We determined the draft genome sequence of PS128 and compared its TA-related genes with those of other L. plantarum strains. Shared genomic features with respect to TA-related subsystems may be important clues to the mechanism by which L. plantarum regulates its host immune responses, but unique TA-related genetic determinants should be further investigated to elucidate strain-specific immunomodulatory effects.

Hierarchical clustering of breast cancer methylomes revealed differentially methylated and expressed breast cancer genes.

Oncogenic transformation of normal cells often involves epigenetic alterations, including histone modification and DNA methylation. We conducted whole-genome bisulfite sequencing to determine the DNA methylomes of normal breast, fibroadenoma, invasive ductal carcinomas and MCF7. The emergence, disappearance, expansion and contraction of kilobase-sized hypomethylated regions (HMRs) and the hypomethylation of the megabase-sized partially methylated domains (PMDs) are the major forms of methylation changes observed in breast tumor samples. Hierarchical clustering of HMR revealed tumor-specific hypermethylated clusters and differential methylated enhancers specific to normal or breast cancer cell lines. Joint analysis of gene expression and DNA methylation data of normal breast and breast cancer cells identified differentially methylated and expressed genes associated with breast and/or ovarian cancers in cancer-specific HMR clusters. Furthermore, aberrant patterns of X-chromosome inactivation (XCI) was found in breast cancer cell lines as well as breast tumor samples in the TCGA BRCA (breast invasive carcinoma) dataset. They were characterized with differentially hypermethylated XIST promoter, reduced expression of XIST, and over-expression of hypomethylated X-linked genes. High expressions of these genes were significantly associated with lower survival rates in breast cancer patients. Comprehensive analysis of the normal and breast tumor methylomes suggests selective targeting of DNA methylation changes during breast cancer progression. The weak causal relationship between DNA methylation and gene expression observed in this study is evident of more complex role of DNA methylation in the regulation of gene expression in human epigenetics that deserves further investigation.

Diagnostics for stochastic genome-scale modeling via model slicing and debugging.

Modeling of biological behavior has evolved from simple gene expression plots represented by mathematical equations to genome-scale systems biology networks. However, due to obstacles in complexity and scalability of creating genome-scale models, several biological modelers have turned to programming or scripting languages and away from modeling fundamentals. In doing so, they have traded the ability to have exchangeable, standardized model representation formats, while those that remain true to standardized model representation are faced with challenges in model complexity and analysis. We have developed a model diagnostic methodology inspired by program slicing and debugging and demonstrate the effectiveness of the methodology on a genome-scale metabolic network model published in the BioModels database. The computer-aided identification revealed specific points of interest such as reversibility of reactions, initialization of species amounts, and parameter estimation that improved a candidate cell's adenosine triphosphate production. We then compared the advantages of our methodology over other modeling techniques such as model checking and model reduction. A software application that implements the methodology is available at http://gel.ym.edu.tw/gcs/.

Microbial diversity analysis of fermented mung beans (Lu-Doh-Huang) by using pyrosequencing and culture methods.

In Taiwanese alternative medicine Lu-doh-huang (also called Pracparatum mungo), mung beans are mixed with various herbal medicines and undergo a 4-stage process of anaerobic fermentation. Here we used high-throughput sequencing of the 16S rRNA gene to profile the bacterial community structure of Lu-doh-huang samples. Pyrosequencing of samples obtained at 7 points during fermentation revealed 9 phyla, 264 genera, and 586 species of bacteria. While mung beans were inside bamboo sections (stages 1 and 2 of the fermentation process), family Lactobacillaceae and genus Lactobacillus emerged in highest abundance; Lactobacillus plantarum was broadly distributed among these samples. During stage 3, the bacterial distribution shifted to family Porphyromonadaceae, and Butyricimonas virosa became the predominant microbial component. Thereafter, bacterial counts decreased dramatically, and organisms were too few to be detected during stage 4. In addition, the microbial compositions of the liquids used for soaking bamboo sections were dramatically different: Exiguobacterium mexicanum predominated in the fermented soybean solution whereas B. virosa was predominant in running spring water. Furthermore, our results from pyrosequencing paralleled those we obtained by using the traditional culture method, which targets lactic acid bacteria. In conclusion, the microbial communities during Lu-doh-huang fermentation were markedly diverse, and pyrosequencing revealed a complete picture of the microbial consortium.

Draft Genome Sequence of Lactobacillus pobuzihii E100301T.

Lactobacillus pobuzihii E100301(T) is a novel Lactobacillus species previously isolated from pobuzihi (fermented cummingcordia) in Taiwan. Phylogenetically, this strain is closest to Lactobacillus acidipiscis, but its phenotypic characteristics can be clearly distinguished from those of L. acidipiscis. We present the draft genome sequence of strain L. pobuzihii E100301(T).