Identification and molecular mapping of Rps11, a novel gene conferring resistance to Phytophthora sojae in soybean.

KEY MESSAGE: Rps11 confers excellent resistance to predominant Phytophthora sojae isolates capable of defeating major Rps genes deployed into soybean production, representing a novel source of resistance for soybean cultivar enhancement. ABSTRACT: Phytophthora root and stem rot (PRSR), caused by the soil-borne pathogen Phytophthora sojae, is a devastating disease of soybean [Glycine max (L.) Merr.] throughout the world. Deploying resistant soybean cultivars is the most effective and environmentally friendly approach to managing this disease. The soybean landrace PI 594527 was found to carry excellent resistance to all P. sojae isolates examined, some of which were capable of overcoming the major Rps genesp, such as Rps1-k, Rps1-c, and Rps3-a, predominantly used for soybean protection in the past decades. A mapping population consisting of 58 F2 individuals and 209 F2:3 families derived from a cross between PI 594527 and the susceptible cultivar 'Williams' was used to characterize the inheritance pattern of the resistance to P. soja (Rps) in PI 594527. It was found that the resistance was conferred by a single Rps gene, designated Rps11, which was initially defined as an ~5 Mb genomic region at the beginning of chromosome 7 by bulked segregant analysis (BSA) with a nucleotide polymorphism (SNP) chip comprising 7039 SNP markers. Subsequently, simple sequence repeat (SSR) markers in the defined region were used to genotype the F2:3 mapping population to map Rps11 to a 225.3 kb genomic region flanked by SSR markers BARCSOYSSR_07_0286 and BARCSOYSSR_07_0300, according to the soybean reference genome sequence. Particularly, an SSR marker (i.e., BARCSOYSSR_07_0295) was found to tightly co-segregate with Rps11 in the mapping population and can be effectively used for marker-assisted selection of this gene for development of resistant soybean cultivars.

Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus.

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

A Sorghum bicolor expression atlas reveals dynamic genotype-specific expression profiles for vegetative tissues of grain, sweet and bioenergy sorghums.

BACKGROUND: Effective improvement in sorghum crop development necessitates a genomics-based approach to identify functional genes and QTLs. Sequenced in 2009, a comprehensive annotation of the sorghum genome and the development of functional genomics resources is key to enable the discovery and deployment of regulatory and metabolic genes and gene networks for crop improvement. RESULTS: This study utilizes the first commercially available whole-transcriptome sorghum microarray (Sorgh-WTa520972F) to identify tissue and genotype-specific expression patterns for all identified Sorghum bicolor exons and UTRs. The genechip contains 1,026,373 probes covering 149,182 exons (27,577 genes) across the Sorghum bicolor nuclear, chloroplast, and mitochondrial genomes. Specific probesets were also included for putative non-coding RNAs that may play a role in gene regulation (e.g., microRNAs), and confirmed functional small RNAs in related species (maize and sugarcane) were also included in our array design. We generated expression data for 78 samples with a combination of four different tissue types (shoot, root, leaf and stem), two dissected stem tissues (pith and rind) and six diverse genotypes, which included 6 public sorghum lines (R159, Atlas, Fremont, PI152611, AR2400 and PI455230) representing grain, sweet, forage, and high biomass ideotypes. CONCLUSIONS: Here we present a summary of the microarray dataset, including analysis of tissue-specific gene expression profiles and associated expression profiles of relevant metabolic pathways. With an aim to enable identification and functional characterization of genes in sorghum, this expression atlas presents a new and valuable resource to the research community.

An attenuated Lassa vaccine in SIV-infected rhesus macaques does not persist or cause arenavirus disease but does elicit Lassa virus-specific immunity.

BACKGROUND: Lassa hemorrhagic fever (LHF) is a rodent-borne viral disease that can be fatal for human beings. In this study, an attenuated Lassa vaccine candidate, ML29, was tested in SIV-infected rhesus macaques for its ability to elicit immune responses without instigating signs pathognomonic for arenavirus disease. ML29 is a reassortant between Lassa and Mopeia viruses that causes a transient infection in non-human primates and confers sterilizing protection from lethal Lassa viral challenge. However, since the LHF endemic area of West Africa also has high HIV seroprevalence, it is important to determine whether vaccination could be safe in the context of HIV infection. RESULTS: SIV-infected and uninfected rhesus macaques were vaccinated with the ML29 virus and monitored for specific humoral and cellular immune responses, as well as for classical and non-classical signs of arenavirus disease. Classical disease signs included viremia, rash, respiratory distress, malaise, high liver enzyme levels, and virus invasion of the central nervous system. Non-classical signs, derived from profiling the blood transcriptome of virulent and non-virulent arenavirus infections, included increased expression of interferon-stimulated genes (ISG) and decreased expression of COX2, IL-1ß, coagulation intermediates and nuclear receptors needed for stress signaling. All vaccinated monkeys showed ML29-specific antibody responses and ML29-specific cell-mediated immunity. CONCLUSION: SIV-infected and uninfected rhesus macaques responded similarly to ML29 vaccination, and none developed chronic arenavirus infection. Importantly, none of the macaques developed signs, classical or non-classical, of arenavirus disease.

Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis.

A synergistic combination of two next-generation sequencing platforms with a detailed comparative BAC physical contig map provided a cost-effective assembly of the genome sequence of the domestic turkey (Meleagris gallopavo). Heterozygosity of the sequenced source genome allowed discovery of more than 600,000 high quality single nucleotide variants. Despite this heterozygosity, the current genome assembly (∼1.1 Gb) includes 917 Mb of sequence assigned to specific turkey chromosomes. Annotation identified nearly 16,000 genes, with 15,093 recognized as protein coding and 611 as non-coding RNA genes. Comparative analysis of the turkey, chicken, and zebra finch genomes, and comparing avian to mammalian species, supports the characteristic stability of avian genomes and identifies genes unique to the avian lineage. Clear differences are seen in number and variety of genes of the avian immune system where expansions and novel genes are less frequent than examples of gene loss. The turkey genome sequence provides resources to further understand the evolution of vertebrate genomes and genetic variation underlying economically important quantitative traits in poultry. This integrated approach may be a model for providing both gene and chromosome level assemblies of other species with agricultural, ecological, and evolutionary interest.

An essential role for the antiviral endoribonuclease, RNase-L, in antibacterial immunity.

Type I IFNs were discovered as the primary antiviral cytokines and are now known to serve critical functions in host defense against bacterial pathogens. Accordingly, established mediators of IFN antiviral activity may mediate previously unrecognized antibacterial functions. RNase-L is the terminal component of an RNA decay pathway that is an important mediator of IFN-induced antiviral activity. Here, we identify a role for RNase-L in the host antibacterial response. RNase-L(-/-) mice exhibited a dramatic increase in mortality after challenge with Bacillus anthracis and Escherichia coli; this increased susceptibility was due to a compromised immune response resulting in increased bacterial load. Investigation of the mechanisms of RNase-L antibacterial activity indicated that RNase-L is required for the optimal induction of proinflammatory cytokines that play essential roles in host defense from bacterial pathogens. RNase-L also regulated the expression of the endolysosomal protease, cathepsin-E, and endosome-associated activities, that function to eliminate internalized bacteria and may contribute to RNase-L antimicrobial action. Our results reveal a unique role for RNase-L in the antibacterial response that is mediated through multiple mechanisms. As a regulator of fundamental components of the innate immune response, RNase-L represents a viable therapeutic target to augment host defense against diverse microbial pathogens.

A giant NLR gene confers broad-spectrum resistance to Phytophthora sojae in soybean.

Phytophthora root and stem rot caused by P. sojae is a destructive soybean soil-borne disease found worldwide. Discovery of genes conferring broad-spectrum resistance to the pathogen is a need to prevent the outbreak of the disease. Here, we show that soybean Rps11 is a 27.7-kb nucleotide-binding site-leucine-rich repeat (NBS-LRR or NLR) gene conferring broad-spectrum resistance to the pathogen. Rps11 is located in a genomic region harboring a cluster of large NLR genes of a single origin in soybean, and is derived from rounds of unequal recombination. Such events result in promoter fusion and LRR expansion that may contribute to the broad resistance spectrum. The NLR gene cluster exhibits drastic structural diversification among phylogenetically representative varieties, including gene copy number variation ranging from five to 23 copies, and absence of allelic copies of Rps11 in any of the non-Rps11-donor varieties examined, exemplifying innovative evolution of NLR genes and NLR gene clusters.

Transcriptome sequencing and comparative analysis of cucumber flowers with different sex types.

BACKGROUND: Cucumber, Cucumis sativus L., is an economically and nutritionally important crop of the Cucurbitaceae family and has long served as a primary model system for sex determination studies. Recently, the sequencing of its whole genome has been completed. However, transcriptome information of this species is still scarce, with a total of around 8,000 Expressed Sequence Tag (EST) and mRNA sequences currently available in GenBank. In order to gain more insights into molecular mechanisms of plant sex determination and provide the community a functional genomics resource that will facilitate cucurbit research and breeding, we performed transcriptome sequencing of cucumber flower buds of two near-isogenic lines, WI1983G, a gynoecious plant which bears only pistillate flowers, and WI1983H, a hermaphroditic plant which bears only bisexual flowers. RESULT: Using Roche-454 massive parallel pyrosequencing technology, we generated a total of 353,941 high quality EST sequences with an average length of 175bp, among which 188,255 were from gynoecious flowers and 165,686 from hermaphroditic flowers. These EST sequences, together with approximately 5,600 high quality cucumber EST and mRNA sequences available in GenBank, were clustered and assembled into 81,401 unigenes, of which 28,452 were contigs and 52,949 were singletons. The unigenes and ESTs were further mapped to the cucumber genome and more than 500 alternative splicing events were identified in 443 cucumber genes. The unigenes were further functionally annotated by comparing their sequences to different protein and functional domain databases and assigned with Gene Ontology (GO) terms. A biochemical pathway database containing 343 predicted pathways was also created based on the annotations of the unigenes. Digital expression analysis identified approximately 200 differentially expressed genes between flowers of WI1983G and WI1983H and provided novel insights into molecular mechanisms of plant sex determination process. Furthermore, a set of SSR motifs and high confidence SNPs between WI1983G and WI1983H were identified from the ESTs, which provided the material basis for future genetic linkage and QTL analysis. CONCLUSION: A large set of EST sequences were generated from cucumber flower buds of two different sex types. Differentially expressed genes between these two different sex-type flowers, as well as putative SSR and SNP markers, were identified. These EST sequences provide valuable information to further understand molecular mechanisms of plant sex determination process and forms a rich resource for future functional genomics analysis, marker development and cucumber breeding.

Rapid and reliable single nucleotide polymorphism-based differentiation of Brucella live vaccine strains from field strains.

The reliable differentiation of live Brucella vaccine strains from field isolates is an important element in brucellosis control programs. We describe the design, validation, and implementation of a novel single nucleotide polymorphism (SNP)-based typing platform that offers a rapid, reliable, and robust tool to achieve this with improved diagnostic accuracy compared to existing molecular tests. Furthermore, the assays described are designed such that they supplement, and can be run as an intrinsic part of, a previously described assay identifying Brucella isolates to the species level (K. K. Gopaul, C. J. Smith, M. S. Koylass, and A. M. Whatmore, BMC Microbiol. 8:86), giving a comprehensive molecular typing platform.

Suppression of LPS-induced Interferon-gamma and nitric oxide in splenic lymphocytes by select estrogen-regulated microRNAs: a novel mechanism of immune modulation.

MicroRNAs (miRNAs), recently identified noncoding small RNAs, are emerging as key regulators in homeostasis of the immune system. Therefore, aberrant expression of miRNAs may be linked to immune dysfunction, such as in chronic inflammation and autoimmunity. In this study, we investigated the potential role of miRNAs in estrogen-mediated regulation of innate immune responses, as indicated by up-regulation of lipopolysaccharide (LPS)-induced interferon-gamma (IFNgamma), inducible nitric oxide synthase (iNOS), and nitric oxide in splenic lymphocytes from estrogen-treated mice. We found that miR-146a, a negative regulator of Toll-like receptor (TLR) signaling, was decreased in freshly isolated splenic lymphocytes from estrogen-treated mice compared with placebo controls. Increasing the activity of miR-146a significantly inhibited LPS-induced IFNgamma and iNOS expression in mouse splenic lymphocytes. Further, miRNA microarray and real-time reverse transcriptase-polymerase chain reaction (RT-PCR) analysis revealed that estrogen selectively up-regulates/down-regulates the expression of miRNAs in mouse splenic lymphocytes. miR-223, which is markedly enhanced by estrogen, regulates LPS-induced IFNgamma, but not iNOS or nitric oxide in splenic lymphocytes. Inhibition of miR-223 activity decreased LPS-induced IFNgamma in splenic lymphocytes from estrogen-treated mice. Our data are the first to demonstrate the selective regulation of miRNA expression in immune cells by estrogen and are indicative of an important role of miRNAs in estrogen-mediated immune regulation.

Modulation of hepatic PPAR expression during Ft LVS LPS-induced protection from Francisella tularensis LVS infection.

BACKGROUND: It has been shown previously that administration of Francisella tularensis (Ft) Live Vaccine Strain (LVS) lipopolysaccharide (LPS) protects mice against subsequent challenge with Ft LVS and blunts the pro-inflammatory cytokine response. METHODS: To further investigate the molecular mechanisms that underlie Ft LVS LPS-mediated protection, we profiled global hepatic gene expression following Ft LVS LPS or saline pre-treatment and subsequent Ft LVS challenge using Affymetrix arrays. RESULTS: A large number of genes (> 3,000) were differentially expressed at 48 hours post-infection. The degree of modulation of inflammatory genes by infection was clearly attenuated by pre-treatment with Ft LVS LPS in the surviving mice. However, Ft LVS LPS alone had a subtle effect on the gene expression profile of the uninfected mice. By employing gene set enrichment analysis, we discovered significant up-regulation of the fatty acid metabolism pathway, which is regulated by peroxisome proliferator activated receptors (PPARs). CONCLUSIONS: We hypothesize that the LPS-induced blunting of pro-inflammatory response in mouse is, in part, mediated by PPARs (alpha and gamma).

Transcriptome sequencing of the Microarray Quality Control (MAQC) RNA reference samples using next generation sequencing.

BACKGROUND: Transcriptome sequencing using next-generation sequencing platforms will soon be competing with DNA microarray technologies for global gene expression analysis. As a preliminary evaluation of these promising technologies, we performed deep sequencing of cDNA synthesized from the Microarray Quality Control (MAQC) reference RNA samples using Roche's 454 Genome Sequencer FLX. RESULTS: We generated more that 3.6 million sequence reads of average length 250 bp for the MAQC A and B samples and introduced a data analysis pipeline for translating cDNA read counts into gene expression levels. Using BLAST, 90% of the reads mapped to the human genome and 64% of the reads mapped to the RefSeq database of well annotated genes with e-values

Gene expression in primate liver during viral hemorrhagic fever.

BACKGROUND: Rhesus macaques infected with lymphocytic choriomeningitis virus (LCMV) provide a model for human Lassa fever. Disease begins with flu-like symptoms and progresses rapidly with fatal consequences. Previously, we profiled the blood transcriptome of LCMV-infected monkeys (M. Djavani et al J. Virol. 2007) showing distinct pre-viremic and viremic stages that discriminated virulent from benign infections. In the present study, changes in liver gene expression from macaques infected with virulent LCMV-WE were compared to gene expression in uninfected monkeys as well as to monkeys that were infected but not diseased. RESULTS: Based on a functional pathway analysis of differentially expressed genes, virulent LCMV-WE had a broader effect on liver cell function than did infection with non-virulent LCMV-Armstrong. During the first few days after infection, LCMV altered expression of genes associated with energy production, including fatty acid and glucose metabolism. The transcriptome profile resembled that of an organism in starvation: mRNA for acetyl-CoA carboxylase, a key enzyme of fatty acid synthesis was reduced while genes for enzymes in gluconeogenesis were up-regulated. Expression was also altered for genes associated with complement and coagulation cascades, and with signaling pathways involving STAT1 and TGF-beta. CONCLUSION: Most of the 4500 differentially expressed transcripts represented a general response to both virulent and mild infections. However, approximately 250 of these transcripts had significantly different expression in virulent infections as compared to mild infections, with approximately 30 of these being differentially regulated during the pre-viremic stage of infection. The genes that are expressed early and differently in mild and virulent disease are potential biomarkers for prognosis and triage of acute viral disease.

Characterization of an Entamoeba histolytica high-mobility-group box protein induced during intestinal infection.

The unicellular eukaryote Entamoeba histolytica is a human parasite that causes amebic dysentery and liver abscess. A genome-wide analysis of gene expression modulated by intestinal colonization and invasion identified an upregulated transcript that encoded a putative high-mobility-group box (HMGB) protein, EhHMGB1. We tested if EhHMGB1 encoded a functional HMGB protein and determined its role in control of parasite gene expression. Recombinant EhHMGB1 was able to bend DNA in vitro, a characteristic of HMGB proteins. Core conserved residues required for DNA bending activity in other HMGB proteins were demonstrated by mutational analysis to be essential for EhHMGB1 activity. EhHMGB1 was also able to enhance the binding of human p53 to its cognate DNA sequence in vitro, which is expected for an HMGB1 protein. Confocal microscopy, using antibodies against the recombinant protein, confirmed its nuclear localization. Overexpression of EhHMGB1 in HM1:IMSS trophozoites led to modulation of 33 transcripts involved in a variety of cellular functions. Of these, 20 were also modulated at either day 1 or day 29 in the mouse model of intestinal amebiasis. Notably, four transcripts with known roles in virulence, including two encoding Gal/GalNAc lectin light chains, were modulated in response to EhHMGB1 overexpression. We concluded that EhHMGB1 was a bona fide HMGB protein with the capacity to recapitulate part of the modulation of parasite gene expression seen during adaptation to the host intestine.

Gene expression profiling of nonhuman primates exposed to aerosolized Venezuelan equine encephalitis virus.

Host responses to Venezuelan equine encephalitis viruses (VEEV) were studied in cynomolgus macaques after aerosol exposure to the epizootic virus. Changes in global gene expression were assessed for the brain, lungs, and spleen. In the brain, major histocompatibility complex (MHC) class I transcripts were induced, while the expression of S100b, a factor associated with brain injury, was inhibited, as was expression of the encephalitogenic gene MOG. Cytokine-mediated signals were affected by infection, including those involving IFN-mediated antiviral activity (IRF-7, OAS, and Mx transcripts), and the increased transcription of caspases. Induction of a few immunologically relevant genes (e.g. IFITM1 and STAT1) was common to all tested tissues. Herein, both tissue-specific and nontissue specific transcriptional changes in response to VEEV are described, including induction of IFN-regulated transcripts and cytokine-induced apoptotic factors, in addition to cellular factors in the brain that may be descriptive of the health status of the brain during the infectious process. Altogether, this work provides novel information on common and tissue-specific host responses against VEEV in a nonhuman primate model of aerosol exposure.

Impact of intestinal colonization and invasion on the Entamoeba histolytica transcriptome.

A genome-wide transcriptional analysis of Entamoeba histolytica was performed on trophozoites isolated from the colon of six infected mice and from in vitro culture. An Affymetrix platform gene expression array was designed for this analysis that included probe sets for 9435 open reading frames (ORFs) and 9066 5' and 3' flanking regions. Transcripts were detected for > 80% of all ORFs. A total of 523 transcripts (5.2% of all E. histolytica genes) were significantly changed in amebae isolated from the intestine on Days 1 and 29 after infection: 326 and 109 solely on Days 1 and 29, and 88 on both days. Quantitative real-time reverse transcriptase PCR confirmed these changes in 11/12 genes tested using mRNA isolated from an additional six mice. Adaptation to the intestinal environment was accompanied by increases in a subset of cell signaling genes including transmembrane kinases, ras and rho family GTPases, and calcium binding proteins. Significant decreases in mRNA abundance for genes involved in glycolysis and concomitant increases in lipases were consistent with a change in energy metabolism. Defense against bacteria present in the intestine (but lacking from in vitro culture) was suggested by alterations in mRNA levels of genes similar to the AIG1 plant antibacterial proteins. Decreases in oxygen detoxification pathways were observed as expected in the anaerobic colonic lumen. Of the known virulence factors the most remarkable changes were a 20-35-fold increase in a cysteine proteinase four-like gene, and a 2-3-fold decrease in two members of the Gal/GalNAc lectin light subunit family. Control of the observed changes in mRNA abundance in the intestine might potentially rest with four related proteins with DNA binding domains that were down-regulated 6-16-fold in the intestinal environment. In conclusion, the first genome-wide analysis of the transcriptome of E. histolytica demonstrated that the vast majority of genes are transcribed in trophozoites, and that in the host intestine trophozoites altered the expression of mRNAs for genes implicated in metabolism, oxygen defense, cell signaling, virulence, antibacterial activity, and DNA binding.

Open architecture expression profiling of plant transcriptomes and gene discovery using GeneCalling technology.

The recent rapid developments in genomics tools, technologies, and bioinformatics have revolutionized gene expression analysis. It is now routine to measure gene expression modulation at the genomic level. GeneCalling technology is an open architecture system capable of assaying more than 95% of genes expressed in a tissue. Unlike the closed systems, GeneCalling is not dependent upon an existing sequence or clone database. GeneCalling uses as low as 50 pg of the cDNA from samples and identifies cDNA fragments that are differentially modulated within a set of samples. With the use of 96 pairs of restriction enzymes, more than 30,000 cDNA fragments are routinely assayed to identify those that are differentially modulated. Specific processes, such as SeqCalling, Trace Poisoning, and GeneCall Poisoning, are set up to not only confirm the known genes, but also to clone and analyze unknown and novel genes that have an interesting expression profile. GeneCalling has been successfully applied to expression profiling of several plant and fungal species, and resulted in identification and characterization of genes that are useful in commercial applications towards improving agriculturally important traits in plants.

Transcriptome analysis of human peripheral blood mononuclear cells exposed to Lassa virus and to the attenuated Mopeia/Lassa reassortant 29 (ML29), a vaccine candidate.

Lassa virus (LASV) is the causative agent of Lassa Fever and is responsible for several hundred thousand infections and thousands of deaths annually in West Africa. LASV and the non-pathogenic Mopeia virus (MOPV) are both rodent-borne African arenaviruses. A live attenuated reassortant of MOPV and LASV, designated ML29, protects rodents and primates from LASV challenge and appears to be more attenuated than MOPV. To gain better insight into LASV-induced pathology and mechanism of attenuation we performed gene expression profiling in human peripheral blood mononuclear cells (PBMC) exposed to LASV and the vaccine candidate ML29. PBMC from healthy human subjects were exposed to either LASV or ML29. Although most PBMC are non-permissive for virus replication, they remain susceptible to signal transduction by virus particles. Total RNA was extracted and global gene expression was evaluated during the first 24 hours using high-density microarrays. Results were validated using RT-PCR, flow cytometry and ELISA. LASV and ML29 elicited differential expression of interferon-stimulated genes (ISG), as well as genes involved in apoptosis, NF-kB signaling and the coagulation pathways. These genes could eventually serve as biomarkers to predict disease outcomes. The remarkable differential expression of thrombomodulin, a key regulator of inflammation and coagulation, suggests its involvement with vascular abnormalities and mortality in Lassa fever disease.

Microbial diversity in saliva of oral squamous cell carcinoma.

In the oral cavity, chronic inflammation has been observed at various stages of oral squamous cell carcinomas (OSCC). Such inflammation could result from persistent mucosal or epithelial cell colonization by microorganisms. There is increasing evidence of the involvement of oral bacteria in inflammation, warranting further studies on the association of bacteria with the progression of OSCC. The objective of this study was to evaluate the diversity and relative abundance of bacteria in the saliva of subjects with OSCC. Using 454 parallel DNA sequencing, ∼58,000 PCR amplicons that span the V4-V5 hypervariable region of rRNAs from five subjects were sequenced. Members of eight phyla (divisions) of bacteria were detected. The majority of classified sequences belonged to the phyla Firmicutes (45%) and Bacteroidetes (25%). Further, 52 different genera containing approximately 860 (16.51%) known species were identified and 1077 (67%) sequences belonging to various uncultured bacteria or unclassified groups. The species diversity estimates obtained with abundance-based coverage estimators and Chao1 were greater than published analyses of other microbial profiles from the oral cavity. Fifteen unique phylotypes were present in all three OSCC subjects.