A Seriation Approach for Visualization-Driven Discovery of Co-Expression Patterns in Serial Analysis of Gene Expression (SAGE) Data

Background: Serial Analysis of Gene Expression (SAGE) is a DNA sequencing-based method for large-scale gene expression profiling that provides an alternative to microarray analysis. Most analyses of SAGE data aimed at identifying co-expressed genes have been accomplished using various versions of clustering approaches that often result in a number of false positives. Principal Findings: Here we explore the use of seriation, a statistical approach for ordering sets of objects based on their similarity, for large-scale expression pattern discovery in SAGE data. For this specific task we implement a seriation heuristic we term ‘progressive construction of contigs ’ that constructs local chains of related elements by sequentially rearranging margins of the correlation matrix. We apply the heuristic to the analysis of simulated and experimental SAGE data and compare our results to those obtained with a clustering algorithm developed specifically for SAGE data. We show using simulations that the performance of seriation compares favorably to that of the clustering algorithm on noisy SAGE data. Conclusions: We explore the use of a seriation approach for visualization-based pattern discovery in SAGE data. Using both simulations and experimental data, we demonstrate that seriation is able to identify groups of co-expressed genes more accurately than a clustering algorithm developed specifically for SAGE data. Our results suggest that seriation is a usefu

Effect of the down-regulation of the high Grain Protein Content (GPC) genes on the wheat transcriptome during monocarpic senescence

<p>Abstract</p> <p>Background</p> <p>Increasing the nutrient concentration of wheat grains is important to ameliorate nutritional deficiencies in many parts of the world. Proteins and nutrients in the wheat grain are largely derived from the remobilization of degraded leaf molecules during monocarpic senescence. The down-regulation of the NAC transcription factor <it>Grain Protein Content </it>(<it>GPC</it>) in transgenic wheat plants delays senescence (>3 weeks) and reduces the concentration of protein, Zn and Fe in the grain (>30%), linking senescence and nutrient remobilization.</p> <p>Based on the early and rapid up-regulation of <it>GPC </it>in wheat flag leaves after anthesis, we hypothesized that this transcription factor is an early regulator of monocarpic senescence. To test this hypothesis, we used high-throughput mRNA-seq technologies to characterize the effect of the <it>GPC </it>down-regulation on the wheat flag-leaf transcriptome 12 days after anthesis. At this early stage of senescence <it>GPC </it>transcript levels are significantly lower in transgenic GPC-RNAi plants than in the wild type, but there are still no visible phenotypic differences between genotypes.</p> <p>Results</p> <p>We generated 1.4 million 454 reads from early senescing flag leaves (average ~350 nt) and assembled 1.2 million into 30,497 contigs that were used as a reference to map 145 million Illumina reads from three wild type and four GPC-RNAi plants. Following normalization and statistical testing, we identified a set of 691 genes differentially regulated by <it>GPC </it>(431 ≥ 2-fold change). Transcript level ratios between transgenic and wild type plants showed a high correlation (<it>R </it>= 0.83) between qRT-PCR and Illumina results, providing independent validation of the mRNA-seq approach. A set of differentially expressed genes were analyzed across an early senescence time-course.</p> <p>Conclusions</p> <p>Monocarpic senescence is an active process characterized by large-scale changes in gene expression which begins considerably before the appearance of visual symptoms of senescence. The mRNA-seq approach used here was able to detect small differences in transcript levels during the early stages of senescence. This resulted in an extensive list of <it>GPC</it>-regulated genes, which includes transporters, hormone regulated genes, and transcription factors. These <it>GPC</it>-regulated genes, particularly those up-regulated during senescence, provide valuable entry points to dissect the early stages of monocarpic senescence and nutrient remobilization in wheat.</p

Transcriptomic Analysis of the Salivary Glands of an Invasive Whitefly

<div><h3>Background</h3><p>Some species of the whitefly <em>Bemisia tabaci</em> complex cause tremendous losses to crops worldwide through feeding directly and virus transmission indirectly. The primary salivary glands of whiteflies are critical for their feeding and virus transmission. However, partly due to their tiny size, research on whitefly salivary glands is limited and our knowledge on these glands is scarce.</p> <h3>Methodology/Principal Findings</h3><p>We sequenced the transcriptome of the primary salivary glands of the Mediterranean species of <em>B. tabaci</em> complex using an effective cDNA amplification method in combination with short read sequencing (Illumina). In a single run, we obtained 13,615 unigenes. The quantity of the unigenes obtained from the salivary glands of the whitefly is at least four folds of the salivary gland genes from other plant-sucking insects. To reveal the functions of the primary glands, sequence similarity search and comparisons with the whole transcriptome of the whitefly were performed. The results demonstrated that the genes related to metabolism and transport were significantly enriched in the primary salivary glands. Furthermore, we found that a number of highly expressed genes in the salivary glands might be involved in secretory protein processing, secretion and virus transmission. To identify potential proteins of whitefly saliva, the translated unigenes were put into secretory protein prediction. Finally, 295 genes were predicted to encode secretory proteins and some of them might play important roles in whitefly feeding.</p> <h3>Conclusions/Significance:</h3><p>The combined method of cDNA amplification, Illumina sequencing and <em>de novo</em> assembly is suitable for transcriptomic analysis of tiny organs in insects. Through analysis of the transcriptome, genomic features of the primary salivary glands were dissected and biologically important proteins, especially secreted proteins, were predicted. Our findings provide substantial sequence information for the primary salivary glands of whiteflies and will be the basis for future studies on whitefly-plant interactions and virus transmission.</p> </div

A Systematic Analysis on DNA Methylation and the Expression of Both mRNA and microRNA in Bladder Cancer

Background: DNA methylation aberration and microRNA (miRNA) deregulation have been observed in many types of cancers. A systematic study of methylome and transcriptome in bladder urothelial carcinoma has never been reported. Methodology/Principal Findings: The DNA methylation was profiled by modified methylation-specific digital karyotyping (MMSDK) and the expression of mRNAs and miRNAs was analyzed by digital gene expression (DGE) sequencing in tumors and matched normal adjacent tissues obtained from 9 bladder urothelial carcinoma patients. We found that a set of significantly enriched pathways disrupted in bladder urothelial carcinoma primarily related to "neurogenesis" and "cell differentiation" by integrated analysis of -omics data. Furthermore, we identified an intriguing collection of cancer-related genes that were deregulated at the levels of DNA methylation and mRNA expression, and we validated several of these genes (HIC1, SLIT2, RASAL1, and KRT17) by Bisulfite Sequencing PCR and Reverse Transcription qPCR in a panel of 33 bladder cancer samples. Conclusions/Significance: We characterized the profiles between methylome and transcriptome in bladder urothelial carcinoma, identified a set of significantly enriched key pathways, and screened four aberrantly methylated and expressed genes. Conclusively, our findings shed light on a new avenue for basic bladder cancer research

Transcriptome Analysis of H2O2-Treated Wheat Seedlings Reveals a H2O2-Responsive Fatty Acid Desaturase Gene Participating in Powdery Mildew Resistance

Hydrogen peroxide (H2O2) plays important roles in plant biotic and abiotic stress responses. However, the effect of H2O2 stress on the bread wheat transcriptome is still lacking. To investigate the cellular and metabolic responses triggered by H2O2, we performed an mRNA tag analysis of wheat seedlings under 10 mM H2O2 treatment for 6 hour in one powdery mildew (PM) resistant (PmA) and two susceptible (Cha and Han) lines. In total, 6,156, 6,875 and 3,276 transcripts were found to be differentially expressed in PmA, Han and Cha respectively. Among them, 260 genes exhibited consistent expression patterns in all three wheat lines and may represent a subset of basal H2O2 responsive genes that were associated with cell defense, signal transduction, photosynthesis, carbohydrate metabolism, lipid metabolism, redox homeostasis, and transport. Among genes specific to PmA, ‘transport’ activity was significantly enriched in Gene Ontology analysis. MapMan classification showed that, while both up- and down- regulations were observed for auxin, abscisic acid, and brassinolides signaling genes, the jasmonic acid and ethylene signaling pathway genes were all up-regulated, suggesting H2O2-enhanced JA/Et functions in PmA. To further study whether any of these genes were involved in wheat PM response, 19 H2O2-responsive putative defense related genes were assayed in wheat seedlings infected with Blumeria graminis f. sp. tritici (Bgt). Eight of these genes were found to be co-regulated by H2O2 and Bgt, among which a fatty acid desaturase gene TaFAD was then confirmed by virus induced gene silencing (VIGS) to be required for the PM resistance. Together, our data presents the first global picture of the wheat transcriptome under H2O2 stress and uncovers potential links between H2O2 and Bgt responses, hence providing important candidate genes for the PM resistance in wheat

Genome-Wide Analysis of Heteroduplex DNA in Mismatch Repair–Deficient Yeast Cells Reveals Novel Properties of Meiotic Recombination Pathways

Meiotic DNA double-strand breaks (DSBs) initiate crossover (CO) recombination, which is necessary for accurate chromosome segregation, but DSBs may also repair as non-crossovers (NCOs). Multiple recombination pathways with specific intermediates are expected to lead to COs and NCOs. We revisited the mechanisms of meiotic DSB repair and the regulation of CO formation, by conducting a genome-wide analysis of strand-transfer intermediates associated with recombination events. We performed this analysis in a SK1 × S288C Saccharomyces cerevisiae hybrid lacking the mismatch repair (MMR) protein Msh2, to allow efficient detection of heteroduplex DNAs (hDNAs). First, we observed that the anti-recombinogenic activity of MMR is responsible for a 20% drop in CO number, suggesting that in MMR–proficient cells some DSBs are repaired using the sister chromatid as a template when polymorphisms are present. Second, we observed that a large fraction of NCOs were associated with trans–hDNA tracts constrained to a single chromatid. This unexpected finding is compatible with dissolution of double Holliday junctions (dHJs) during repair, and it suggests the existence of a novel control point for CO formation at the level of the dHJ intermediate, in addition to the previously described control point before the dHJ formation step. Finally, we observed that COs are associated with complex hDNA patterns, confirming that the canonical double-strand break repair model is not sufficient to explain the formation of most COs. We propose that multiple factors contribute to the complexity of recombination intermediates. These factors include repair of nicks and double-stranded gaps, template switches between non-sister and sister chromatids, and HJ branch migration. Finally, the good correlation between the strand transfer properties observed in the absence of and in the presence of Msh2 suggests that the intermediates detected in the absence of Msh2 reflect normal intermediates

Comparative mRNA and microRNA Expression Profiling of Three Genitourinary Cancers Reveals Common Hallmarks and Cancer-Specific Molecular Events

Genome-wide gene expression profile using deep sequencing technologies can drive the discovery of cancer biomarkers and therapeutic targets. Such efforts are often limited to profiling the expression signature of either mRNA or microRNA (miRNA) in a single type of cancer.Here we provided an integrated analysis of the genome-wide mRNA and miRNA expression profiles of three different genitourinary cancers: carcinomas of the bladder, kidney and testis.Our results highlight the general or cancer-specific roles of several genes and miRNAs that may serve as candidate oncogenes or suppressors of tumor development. Further comparative analyses at the systems level revealed that significant aberrations of the cell adhesion process, p53 signaling, calcium signaling, the ECM-receptor and cell cycle pathways, the DNA repair and replication processes and the immune and inflammatory response processes were the common hallmarks of human cancers. Gene sets showing testicular cancer-specific deregulation patterns were mainly implicated in processes related to male reproductive function, and general disruptions of multiple metabolic pathways and processes related to cell migration were the characteristic molecular events for renal and bladder cancer, respectively. Furthermore, we also demonstrated that tumors with the same histological origins and genes with similar functions tended to group together in a clustering analysis. By assessing the correlation between the expression of each miRNA and its targets, we determined that deregulation of 'key' miRNAs may result in the global aberration of one or more pathways or processes as a whole.This systematic analysis deciphered the molecular phenotypes of three genitourinary cancers and investigated their variations at the miRNA level simultaneously. Our results provided a valuable source for future studies and highlighted some promising genes, miRNAs, pathways and processes that may be useful for diagnostic or therapeutic applications

LNCaP Atlas: Gene expression associated with in vivo progression to castration-recurrent prostate cancer

<p>Abstract</p> <p>Background</p> <p>There is no cure for castration-recurrent prostate cancer (CRPC) and the mechanisms underlying this stage of the disease are unknown.</p> <p>Methods</p> <p>We analyzed the transcriptome of human LNCaP prostate cancer cells as they progress to CRPC <it>in vivo </it>using replicate LongSAGE libraries. We refer to these libraries as the LNCaP atlas and compared these gene expression profiles with current suggested models of CRPC.</p> <p>Results</p> <p>Three million tags were sequenced using <it>in vivo </it>samples at various stages of hormonal progression to reveal 96 novel genes differentially expressed in CRPC. Thirty-one genes encode proteins that are either secreted or are located at the plasma membrane, 21 genes changed levels of expression in response to androgen, and 8 genes have enriched expression in the prostate. Expression of 26, 6, 12, and 15 genes have previously been linked to prostate cancer, Gleason grade, progression, and metastasis, respectively. Expression profiles of genes in CRPC support a role for the transcriptional activity of the androgen receptor (<it>CCNH, CUEDC2, FLNA, PSMA7</it>), steroid synthesis and metabolism (<it>DHCR24, DHRS7</it>, <it>ELOVL5, HSD17B4</it>, <it>OPRK1</it>), neuroendocrine (<it>ENO2, MAOA, OPRK1, S100A10, TRPM8</it>), and proliferation (<it>GAS5</it>, <it>GNB2L1</it>, <it>MT-ND3</it>, <it>NKX3-1</it>, <it>PCGEM1</it>, <it>PTGFR</it>, <it>STEAP1</it>, <it>TMEM30A</it>), but neither supported nor discounted a role for cell survival genes.</p> <p>Conclusions</p> <p>The <it>in vivo </it>gene expression atlas for LNCaP was sequenced and support a role for the androgen receptor in CRPC.</p