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.

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.

The genome of woodland strawberry (Fragaria vesca).

The woodland strawberry, Fragaria vesca (2n = 2x = 14), is a versatile experimental plant system. This diminutive herbaceous perennial has a small genome (240 Mb), is amenable to genetic transformation and shares substantial sequence identity with the cultivated strawberry (Fragaria × ananassa) and other economically important rosaceous plants. Here we report the draft F. vesca genome, which was sequenced to ×39 coverage using second-generation technology, assembled de novo and then anchored to the genetic linkage map into seven pseudochromosomes. This diploid strawberry sequence lacks the large genome duplications seen in other rosids. Gene prediction modeling identified 34,809 genes, with most being supported by transcriptome mapping. Genes critical to valuable horticultural traits including flavor, nutritional value and flowering time were identified. Macrosyntenic relationships between Fragaria and Prunus predict a hypothetical ancestral Rosaceae genome that had nine chromosomes. New phylogenetic analysis of 154 protein-coding genes suggests that assignment of Populus to Malvidae, rather than Fabidae, is warranted.

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.

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).

Members of the Entamoeba histolytica transmembrane kinase family play non-redundant roles in growth and phagocytosis.

Entamoeba histolytica contains a large and novel family of transmembrane kinases (TMKs). The expression patterns of the E. histolytica TMKs in individual trophozoites and the roles of the TMKs for sensing and responding to extracellular cues were incompletely characterised. Here we provide evidence that single cells express multiple TMKs and that TMK39 and TMK54 likely serve non-redundant cellular functions. Laser-capture microdissection was used in conjunction with microarray analysis to demonstrate that single trophozoites express more than one TMK gene. Anti-peptide antibodies were raised against unique regions in the extracellular domains of TMK39, TMK54 and PaTMK, and TMK expression was analysed at the protein level. Flow cytometric assays revealed that populations of trophozoites homogeneously expressed TMK39, TMK54 and PaTMK, while confocal microscopy identified different patterns of cell surface expression for TMK39 and TMK54. The functions of TMK39 and TMK54 were probed by the inducible expression of dominant-negative mutants. While TMK39 co-localised with ingested beads and expression of truncated TMK39 interfered with trophozoite phagocytosis of apoptotic lymphocytes, expression of a truncated TMK54 inhibited growth of amoebae and altered the surface expression of the heavy subunit of the E. histolytica Gal/GalNAc lectin. Overall, our data indicates that multiple members of the novel E. histolytica TMK family are utilised for non-redundant functions by the parasite.

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.

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.

Genome sequence of Brucella abortus vaccine strain S19 compared to virulent strains yields candidate virulence genes.

The Brucella abortus strain S19, a spontaneously attenuated strain, has been used as a vaccine strain in vaccination of cattle against brucellosis for six decades. Despite many studies, the physiological and molecular mechanisms causing the attenuation are not known. We have applied pyrosequencing technology together with conventional sequencing to rapidly and comprehensively determine the complete genome sequence of the attenuated Brucella abortus vaccine strain S19. The main goal of this study is to identify candidate virulence genes by systematic comparative analysis of the attenuated strain with the published genome sequences of two virulent and closely related strains of B. abortus, 9-941 and 2308. The two S19 chromosomes are 2,122,487 and 1,161,449 bp in length. A total of 3062 genes were identified and annotated. Pairwise and reciprocal genome comparisons resulted in a total of 263 genes that were non-identical between the S19 genome and any of the two virulent strains. Amongst these, 45 genes were consistently different between the attenuated strain and the two virulent strains but were identical amongst the virulent strains, which included only two of the 236 genes that have been implicated as virulence factors in literature. The functional analyses of the differences have revealed a total of 24 genes that may be associated with the loss of virulence in S19. Of particular relevance are four genes with more than 60 bp consistent difference in S19 compared to both the virulent strains, which, in the virulent strains, encode an outer membrane protein and three proteins involved in erythritol uptake or metabolism.

Early blood profiles of virus infection in a monkey model for Lassa fever.

Acute arenavirus disease in primates, like Lassa hemorrhagic fever in humans, begins with flu-like symptoms and leads to death approximately 2 weeks after infection. Our goal was to identify molecular changes in blood that are related to disease progression. Rhesus macaques (Macaca mulatta) infected intravenously with a lethal dose of lymphocytic choriomeningitis virus (LCMV) provide a model for Lassa virus infection of humans. Blood samples taken before and during the course of infection were used to monitor gene expression changes that paralleled disease onset. Changes in blood showed major disruptions in eicosanoid, immune response, and hormone response pathways. Approximately 12% of host genes alter their expression after LCMV infection, and a subset of these genes can discriminate between virulent and non-virulent LCMV infection. Major transcription changes have been given preliminary confirmation by quantitative PCR and protein studies and will be valuable candidates for future validation as biomarkers for arenavirus disease.

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.

Comprehensive transcript profiling of Pto- and Prf-mediated host defense responses to infection by Pseudomonas syringae pv. tomato.

The disease resistance gene Pto encodes a serine/threonine protein kinase that confers resistance in tomato to Pseudomonas syringae pv. tomato strains that express the effector protein AvrPto. Pto-mediated resistance to bacterial speck disease also requires Prf, a protein with leucine-rich repeats and a putative nucleotide-binding site, although the role of Prf in the defense pathway is not known. We used GeneCalling, an open-architecture, mRNA-profiling technology, to identify genes that are either induced or suppressed in leaves 4 h after bacterial infection in the Pto- and Prf-mediated tomato-Pseudomonas(avrPto) interaction. Over 135 000 individual cDNA fragments representing an estimated 90% of the transcripts expressed in tomato leaves were examined and 432 differentially expressed genes were identified. The genes encode over 25 classes of proteins including 11 types of transcription factors and many signal transduction components. Differential expression of 91% of the genes required both Pto and Prf. Interestingly, differential expression of 32 genes did not require Pto but was dependent on Prf. Thus, our data support a role for Prf early in the Pto pathway and indicate that Prf can also function as an independent host recognition determinant of bacterial infection. Comprehensive expression profiling of the Pto-mediated defense response allows the development of many new hypotheses about the molecular basis of resistance to bacterial speck disease.