Gene expression in primate liver during viral hemorrhagic fever

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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-β.</p> <p>Conclusion</p> <p>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.</p

A Versatile Computational Pipeline for Bacterial Genome Annotation Improvement and Comparative Analysis, with \u3cem\u3eBrucella\u3c/em\u3e as a Use Case

We present a bacterial genome computational analysis pipeline, called GenVar. The pipeline, based on the program GeneWise, is designed to analyze an annotated genome and automatically identify missed gene calls and sequence variants such as genes with disrupted reading frames (split genes) and those with insertions and deletions (indels). For a given genome to be analyzed, GenVar relies on a database containing closely related genomes (such as other species or strains) as well as a few additional reference genomes. GenVar also helps identify gene disruptions probably caused by sequencing errors. We exemplify GenVar’s capabilities by presenting results from the analysis of four Brucella genomes. Brucella is an important human pathogen and zoonotic agent. The analysis revealed hundreds of missed gene calls, new split genes and indels, several of which are species specific and hence provide valuable clues to the understanding of the genome basis of Brucella pathogenicity and host specificity

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

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

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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 ≤ 10^-20. We measured gene expression levels in the A and B samples by counting the numbers of reads that mapped to individual RefSeq genes in multiple sequencing runs to evaluate the MAQC quality metrics for reproducibility, sensitivity, specificity, and accuracy and compared the results with DNA microarrays and Quantitative RT-PCR (QRTPCR) from the MAQC studies. In addition, 88% of the reads were successfully aligned directly to the human genome using the AceView alignment programs with an average 90% sequence similarity to identify 137,899 unique exon junctions, including 22,193 new exon junctions not yet contained in the RefSeq database.</p> <p>Conclusion</p> <p>Using the MAQC metrics for evaluating the performance of gene expression platforms, the ExpressSeq results for gene expression levels showed excellent reproducibility, sensitivity, and specificity that improved systematically with increasing shotgun sequencing depth, and quantitative accuracy that was comparable to DNA microarrays and QRTPCR. In addition, a careful mapping of the reads to the genome using the AceView alignment programs shed new light on the complexity of the human transcriptome including the discovery of thousands of new splice variants.</p

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.Keywords: Family, Disease resistance genes, NBS, Angiosperms, Preservation, Expression, Balance hypothesis, Trichocarpa, Arabidopsis thaliana, Mechanism

Identification of a Common Lupus Disease-Associated microRNA Expression Pattern in Three Different Murine Models of Lupus

Recent reports have shown that microRNAs (miRNAs) regulate vital immunological processes and have emerged as key regulators of immune system development and function. Therefore, it is important to determine miRNA dysregulation and its pathogenic contribution in autoimmune diseases, an aspect not adequately addressed thus far.In this study, we profiled miRNA expressions in splenic lymphocytes from three murine lupus models (MRL-lpr, B6-lpr and NZB/W(F₁)) with different genetic background by miRNA microarray assays and Real-time RT-PCR. Despite the genetic differences among these three lupus stains, a common set of dysregulated miRNAs (miR-182-96-183 cluster, miR-31, and miR-155) was identified in splenocytes when compared with age-matched control mice. The association of these miRNAs with the disease was highlighted by our observation that this miRNA expression pattern was evident in NZB/W mice only at an age when lupus disease is manifested. Further, we have shown that the miRNA dysregulation in MRL-lpr mice was not simply due to the activation of splenocytes. By Real-time RT-PCR, we confirmed that these miRNAs were upregulated in both purified splenic B and T cells from MRL-lpr mice. miR-127 and miR-379, which were greatly upregulated in splenocytes from lpr mice, were moderately increased in diseased NZB/W mice. In addition, Real-time RT-PCR revealed that miR-146a, miR-101a, and miR-17-92 were also markedly upregulated in splenic T, but not B cells from MRL-lpr mice.The identification of common lupus disease-associated miRNAs now forms the basis for the further investigation of the pathogenic contribution of these miRNAs in autoimmune lupus, which will advance our knowledge of the role of miRNAs in autoimmunity. Given that miRNAs are conserved, with regard to both evolution and function, our observation of a common lupus disease-associated miRNA expression pattern in murine lupus models is likely to have significant pathogenic, diagnostic, and/or therapeutic implications in human lupus

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 60bp 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

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

<p>Abstract</p> <p>Background</p> <p>It has been shown previously that administration of <it>Francisella tularensis </it>(<it>Ft</it>) Live Vaccine Strain (LVS) lipopolysaccharide (LPS) protects mice against subsequent challenge with <it>Ft </it>LVS and blunts the pro-inflammatory cytokine response.</p> <p>Methods</p> <p>To further investigate the molecular mechanisms that underlie <it>Ft </it>LVS LPS-mediated protection, we profiled global hepatic gene expression following <it>Ft </it>LVS LPS or saline pre-treatment and subsequent <it>Ft </it>LVS challenge using Affymetrix arrays.</p> <p>Results</p> <p>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 <it>Ft </it>LVS LPS in the surviving mice. However, <it>Ft </it>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).</p> <p>Conclusions</p> <p>We hypothesize that the LPS-induced blunting of pro-inflammatory response in mouse is, in part, mediated by PPARs (α and γ).</p

Targets of the Entamoeba histolytica Transcription Factor URE3-BP

The Entamoeba histolytica transcription factor Upstream Regulatory Element 3-Binding Protein (URE3-BP) is a calcium-responsive regulator of two E. histolytica virulence genes, hgl5 and fdx1. URE3-BP was previously identified by a yeast one-hybrid screen of E. histolytica proteins capable of binding to the sequence TATTCTATT (Upstream Regulatory Element 3 (URE3)) in the promoter regions of hgl5 and fdx1. In this work, precise definition of the consensus URE3 element was performed by electrophoretic mobility shift assays (EMSA) using base-substituted oligonucleotides, and the consensus motif validated using episomal reporter constructs. Transcriptome profiling of a strain induced to produce a dominant-positive URE3-BP was then used to identify additional genes regulated by URE3-BP. Fifty modulated transcripts were identified, and of these the EMSA defined motif T[atg]T[tc][cg]T[at][tgc][tg] was found in over half of the promoters (54% p<0.0001). Fifteen of the URE3-BP regulated genes were potential membrane proteins, suggesting that one function of URE3-BP is to remodel the surface of E. histolytica in response to a calcium signal. Induction of URE3-BP leads to an increase in tranwell migration, suggesting a possible role in the regulation of cellular motility