Integrating Constitutive Gene Expression and Chemoactivity: Mining the NCI60 Anticancer Screen

Studies into the genetic origins of tumor cell chemoactivity pose significant challenges to bioinformatic mining efforts. Connections between measures of gene expression and chemoactivity have the potential to identify clinical biomarkers of compound response, cellular pathways important to efficacy and potential toxicities; all vital to anticancer drug development. An investigation has been conducted that jointly explores tumor-cell constitutive NCI60 gene expression profiles and small-molecule NCI60 growth inhibition chemoactivity profiles, viewed from novel applications of self-organizing maps (SOMs) and pathway-centric analyses of gene expressions, to identify subsets of over- and under-expressed pathway genes that discriminate chemo-sensitive and chemo-insensitive tumor cell types. Linear Discriminant Analysis (LDA) is used to quantify the accuracy of discriminating genes to predict tumor cell chemoactivity. LDA results find 15% higher prediction accuracies, using ∼30% fewer genes, for pathway-derived discriminating genes when compared to genes derived using conventional gene expression-chemoactivity correlations. The proposed pathway-centric data mining procedure was used to derive discriminating genes for ten well-known compounds. Discriminating genes were further evaluated using gene set enrichment analysis (GSEA) to reveal a cellular genetic landscape, comprised of small numbers of key over and under expressed on- and off-target pathway genes, as important for a compound’s tumor cell chemoactivity. Literature-based validations are provided as support for chemo-important pathways derived from this procedure. Qualitatively similar results are found when using gene expression measurements derived from different microarray platforms. The data used in this analysis is available at http://pubchem.ncbi.nlm.nih.gov/and http://www.ncbi.nlm.nih.gov/projects/geo ({"type":"entrez-geo","attrs":{"text":"GPL96","term_id":"96"}}GPL96, {"type":"entrez-geo","attrs":{"text":"GSE32474","term_id":"32474"}}GSE32474).     

Genome re-sequencing suggested a weedy rice origin from domesticated indica-japonica hybridization: a case study from southern China

Main conclusion

Whole-genome re-sequencing of weedy rice from southern China reveals that weedy rice can originate from hybridization of domesticated indica and japonica rice.

Abstract

Weedy rice (Oryza sativa f. spontanea Rosh.), which harbors phenotypes of both wild and domesticated rice, has become one of the most notorious weeds in rice fields worldwide. While its formation is poorly understood, massive amounts of rice genomic data may provide new insights into this issue. In this study, we determined genomes of three weedy rice samples from the lower Yangtze region, China, and investigated their phylogenetics, population structure and chromosomal admixture patterns. The phylogenetic tree and principle component analysis based on 46,005 SNPs with 126 other Oryza accessions suggested that the three weedy rice accessions were intermediate between japonica and indica rice. An ancestry inference study further demonstrated that weedy rice had two dominant genomic components (temperate japonica and indica). This strongly suggests that weedy rice originated from indica-japonica hybridization. Furthermore, 22,443 novel fixed single nucleotide polymorphisms were detected in the weedy genomes and could have been generated after indica-japonica hybridization for environmental adaptation.     

Fitting hidden Markov models of protein domains to a target species: application to Plasmodium falciparum

Background

The identification of gene sets that are significantly impacted in a given condition based on microarray data is a crucial step in current life science research. Most gene set analysis methods treat genes equally, regardless how specific they are to a given gene set.

Results

In this work we propose a new gene set analysis method that computes a gene set score as the mean of absolute values of weighted moderated gene t-scores. The gene weights are designed to emphasize the genes appearing in few gene sets, versus genes that appear in many gene sets. We demonstrate the usefulness of the method when analyzing gene sets that correspond to the KEGG pathways, and hence we called our method P athway A nalysis with D own-weighting of O verlapping G enes (PADOG). Unlike most gene set analysis methods which are validated through the analysis of 2-3 data sets followed by a human interpretation of the results, the validation employed here uses 24 different data sets and a completely objective assessment scheme that makes minimal assumptions and eliminates the need for possibly biased human assessments of the analysis results.

Conclusions

PADOG significantly improves gene set ranking and boosts sensitivity of analysis using information already available in the gene expression profiles and the collection of gene sets to be analyzed. The advantages of PADOG over other existing approaches are shown to be stable to changes in the database of gene sets to be analyzed. PADOG was implemented as an R package available at: http://bioinformaticsprb.med.wayne.edu/PADOG/or http://www.bioconductor.org.     

Different non-host resistance responses of two rice subspecies,japonica and indica,to Puccinia striiformis f. sp. tritici

Key message

Japonica and indica have different non-host resistance (NHR) abilities to Puccinia striiformis f. sp. tritici ( Pst ), and hydrogen peroxide (H ₂ O ₂ ) has a positive function in NHR to japonica against Pst.

Abstract

Non-host interactions between Puccinia striiformis f. sp. tritici (Pst) and two rice subspecies were characterized using 23 rice varieties, including 11 japonica and 12 indica. Results showed that the infected fungal structures were easily produced in the leaves of indica, whereas only several substomatal vesicles and primary infection hyphae were observed in the leaves of japonica. This result indicated that indica is less resistant or more susceptible to Pst than japonica. Hydrogen peroxide accumulated in the initial phase of japonica–Pst interaction but not in indica–Pst interaction. A set of reactive oxygen species (ROS)-related genes was also induced in response to Pst infection, suggesting that ROS activation is one of the major mechanisms of non-host resistance of rice to Pst.     

Gene Expression Browser: large-scale and cross-experiment microarray data integration, management, search & visualization

Background

In the last decade, a large amount of microarray gene expression data has been accumulated in public repositories. Integrating and analyzing high-throughput gene expression data have become key activities for exploring gene functions, gene networks and biological pathways. Effectively utilizing these invaluable microarray data remains challenging due to a lack of powerful tools to integrate large-scale gene-expression information across diverse experiments and to search and visualize a large number of gene-expression data points.

Results

Gene Expression Browser is a microarray data integration, management and processing system with web-based search and visualization functions. An innovative method has been developed to define a treatment over a control for every microarray experiment to standardize and make microarray data from different experiments homogeneous. In the browser, data are pre-processed offline and the resulting data points are visualized online with a 2-layer dynamic web display. Users can view all treatments over control that affect the expression of a selected gene via Gene View, and view all genes that change in a selected treatment over control via treatment over control View. Users can also check the changes of expression profiles of a set of either the treatments over control or genes via Slide View. In addition, the relationships between genes and treatments over control are computed according to gene expression ratio and are shown as co-responsive genes and co-regulation treatments over control.

Conclusion

Gene Expression Browser is composed of a set of software tools, including a data extraction tool, a microarray data-management system, a data-annotation tool, a microarray data-processing pipeline, and a data search & visualization tool. The browser is deployed as a free public web service (http://www.ExpressionBrowser.com) that integrates 301 ATH1 gene microarray experiments from public data repositories (viz. the Gene Expression Omnibus repository at the National Center for Biotechnology Information and Nottingham Arabidopsis Stock Center). The set of Gene Expression Browser software tools can be easily applied to the large-scale expression data generated by other platforms and in other species.     

Meta-expression analysis of unannotated genes in rice and approaches for network construction to suggest the probable roles

Key message

This work suggests 2020 potential candidates in rice for the functional annotation of unannotated genes using meta-analysis of anatomical samples derived from microarray and RNA-seq technologies and this information will be useful to identify novel morphological agronomic traits.

Abstract

Although the genome of rice (Oryza sativa) has been sequenced, 14,365 genes are considered unannotated because they lack putative annotation information. According to the Rice Genome Annotation Project Database (http://rice.plantbiology.msu.edu/), the proportion of functionally characterized unannotated genes (0.35%) is quite limited when compared with the approximately 3.9% of annotated genes with assigned putative functions. Researchers require additional information to help them investigate the molecular mechanisms associated with those unannotated genes. To determine which of them might regulate morphological or physiological traits in the rice genome, we conducted a meta-analysis of expression data that covered a wide range of tissue/organ samples. Overall, 2020 genes showed cultivar-, tissue-, or organ-preferential patterns of expression. Representative candidates from featured groups were validated by RT-PCR, and the GUS reporter system was used to validate the expression of genes that were clustered according to their leaf or root preference. Taking a molecular and genetics approach, we examined meta-expression data and found that 127 genes were differentially expressed between japonica and indica rice cultivars. This is potentially significant for future agronomic applications. We also used a T-DNA insertional mutant and performed a co-expression network analysis of Sword shape dwarf1 (SSD1), a gene that regulates cell division. This network was refined via RT-PCR analysis. Our results suggested that SSD1 represses the expression of four genes related to the processes of DNA replication or cell division and provides insight into possible molecular mechanisms. Together, these strategies present a valuable tool for in-depth characterization of currently unannotated genes.

     

GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists

Background

Although expression microarrays have become a standard tool used by biologists, analysis of data produced by microarray experiments may still present challenges. Comparison of data from different platforms, organisms, and labs may involve complicated data processing, and inferring relationships between genes remains difficult.

Results

S TAR N ET 2 is a new web-based tool that allows post hoc visual analysis of correlations that are derived from expression microarray data. S TAR N ET 2 facilitates user discovery of putative gene regulatory networks in a variety of species (human, rat, mouse, chicken, zebrafish, Drosophila, C. elegans, S. cerevisiae, Arabidopsis and rice) by graphing networks of genes that are closely co-expressed across a large heterogeneous set of preselected microarray experiments. For each of the represented organisms, raw microarray data were retrieved from NCBI's Gene Expression Omnibus for a selected Affymetrix platform. All pairwise Pearson correlation coefficients were computed for expression profiles measured on each platform, respectively. These precompiled results were stored in a MySQL database, and supplemented by additional data retrieved from NCBI. A web-based tool allows user-specified queries of the database, centered at a gene of interest. The result of a query includes graphs of correlation networks, graphs of known interactions involving genes and gene products that are present in the correlation networks, and initial statistical analyses. Two analyses may be performed in parallel to compare networks, which is facilitated by the new H EAT S EEKER module.

Conclusion

S TAR N ET 2 is a useful tool for developing new hypotheses about regulatory relationships between genes and gene products, and has coverage for 10 species. Interpretation of the correlation networks is supported with a database of previously documented interactions, a test for enrichment of Gene Ontology terms, and heat maps of correlation distances that may be used to compare two networks. The list of genes in a S TAR N ET network may be useful in developing a list of candidate genes to use for the inference of causal networks. The tool is freely available at http://vanburenlab.medicine.tamhsc.edu/starnet2.html, and does not require user registration.     

Fine mapping of qHd1, a minor heading date QTL with pleiotropism for yield traits in rice (Oryza sativa L.)

Key message

A minor QTL for heading date located on the long arm of rice chromosome 1 was delimitated to a 95.0-kb region using near isogenic lines with sequential segregating regions.

Abstract

Heading date and grain yield are two key factors determining the commercial potential of a rice variety. In this study, rice populations with sequential segregating regions were developed and used for mapping a minor QTL for heading date, qHd1. A total of 18 populations in six advanced generations through BC₂F₆ to BC₂F₁₁ were derived from a single BC₂F₃ plant of the indica rice cross Zhenshan 97 (ZS97)///ZS97//ZS97/Milyang 46. The QTL was delimitated to a 95.0-kb region flanked by RM12102 and RM12108 in the terminal region of the long arm of chromosome 1. Results also showed that qHd1 was not involved in the photoperiodic response, having an additive effect ranging from 2.4 d to 2.9 d observed in near isogenic lines grown in the paddy field and under the controlled conditions of either short day or long day. The QTL had pleiotropic effects on yield traits, with the ZS97 allele delaying heading and increasing the number of spikelets per panicle, the number of grains per panicle and grain yield per plant. The candidate region contains ten annotated genes including two genes with functional information related to the control of heading date. These results lay a foundation for the cloning of qHd1. In addition, this kind of minor QTLs could be of great significance in rice breeding for allowing minor adjustment of heading date and yield traits.     

Mitochondrial DNA genetic polymorphism in thirteen rice cytoplasmic male sterile lines

Key message

Thirteen rice CMS lines derived from different cytoplasms were classified into eight groups by PCR amplification on mtDNA. The orf79 gene, which causes Boro II CMS, possibly results in Dian1-CMS.

Abstract

Thirteen rice cytoplasmic male sterile (CMS) lines derived from different cytoplasms are widely used for hybrid rice breeding. Based on 27 loci on mitochondrial DNA, including single nucleotide polymorphisms and segmental sequence variations between typical indica and japonica as well as high-polymorphism segmental sequence variations and single nucleotide polymorphisms among rice CMS lines, the 13 rice CMS lines were classified into eight groups: (I) wild-abortive CMS, Indonesian Shuitiangu CMS, K-CMS, Gang CMS, D-CMS and dwarf abortive CMS; (II) Maxie-CMS; (III) Honglian CMS; (IV) Boro II CMS; (V) Dian1-CMS; (VI) Liao-CMS; (VII) Lead CMS; and (VIII) Chinese wild rice CMS. According to their pollen abortion phenotypes, groups I and II (including 7 CMS lines) were classified as sporophytic CMS lines, the cytoplasmic genetic relationships among which were very close. They could have originated from similar, or even the same, cytoplasm donors. Groups III–VIII (including 6 CMS lines) were categorized as gametophytic CMS lines, the cytoplasms of which differed from one another, with some having relatively far genetic relationships. Dian1-CMS was found to harbor the orf79 gene, which causes Boro II CMS, whereas Liao-CMS had an orf79 structure that does not result in Lead CMS. Therefore, we speculated that orf79 is associated with Dian1-CMS but not with Liao-CMS. The atp6–orf79 structure related to sterility was also found to experience multiple evolutionary turnovers. All sporophytic CMS lines were indica-like. Except the Honglian CMS line, which was indica-like, all gametophytic CMS lines were japonica-like.     

Development of universal genetic markers based on single-copy orthologous (COSII) genes in Poaceae

Key message

We develop a set of universal genetic markers based on single-copy orthologous (COSII) genes in Poaceae.

Abstract

Being evolutionary conserved, single-copy orthologous (COSII) genes are particularly useful in comparative mapping and phylogenetic investigation among species. In this study, we identified 2,684 COSII genes based on five sequenced Poaceae genomes including rice, maize, sorghum, foxtail millet, and brachypodium, and then developed 1,072 COSII markers whose transferability and polymorphism among five bamboo species were further evaluated with 46 pairs of randomly selected primers. 91.3 % of the 46 primers obtained clear amplification in at least one bamboo species, and 65.2 % of them produced polymorphism in more than one species. We also used 42 of them to construct the phylogeny for the five bamboo species, and it might reflect more precise evolutionary relationship than the one based on the vegetative morphology. The results indicated a promising prospect of applying these markers to the investigation of genetic diversity and the classification of Poaceae. To ease and facilitate access of the information of common interest to readers, a web-based database of the COSII markers is provided (http://www.sicau.edu.cn/web/yms/PCOSWeb/PCOS.html).     

Multiple sclerosis and pregnancy: what does the patient think? a questionnaire study

Background

The first objective of a DNA microarray experiment is typically to generate a list of genes or probes that are found to be differentially expressed or represented (in the case of comparative genomic hybridizations and/or copy number variation) between two conditions or strains. Rank Products analysis comprises a robust algorithm for deriving such lists from microarray experiments that comprise small numbers of replicates, for example, less than the number required for the commonly used t-test. Currently, users wishing to apply Rank Products analysis to their own microarray data sets have been restricted to the use of command line-based software which can limit its usage within the biological community.

Findings

Here we have developed a web interface to existing Rank Products analysis tools allowing users to quickly process their data in an intuitive and step-wise manner to obtain the respective Rank Product or Rank Sum, probability of false prediction and p-values in a downloadable file.

Conclusions

The online interactive Rank Products analysis tool RankProdIt, for analysis of any data set containing measurements for multiple replicated conditions, is available at: http://strep-microarray.sbs.surrey.ac.uk/RankProducts     

RNAstrand: reading direction of structured RNAs in multiple sequence alignments

Motivation

Genome-wide screens for structured ncRNA genes in mammals, urochordates, and nematodes have predicted thousands of putative ncRNA genes and other structured RNA motifs. A prerequisite for their functional annotation is to determine the reading direction with high precision.

Results

While folding energies of an RNA and its reverse complement are similar, the differences are sufficient at least in conjunction with substitution patterns to discriminate between structured RNAs and their complements. We present here a support vector machine that reliably classifies the reading direction of a structured RNA from a multiple sequence alignment and provides a considerable improvement in classification accuracy over previous approaches.

Software

RNAstrand is freely available as a stand-alone tool from http://www.bioinf.uni-leipzig.de/Software/RNAstrand and is also included in the latest release of RNAz, a part of the Vienna RNA Package.     

PlantGM: a database for genetic markers in rice (Oryza sativa) and Chinese cabbage (Brassica rapa)

The Plant Genetic Map Database (PlantGM) has been developed as a web-based system which provides information about genetic markers in rice (Oryza sativa) and Chinese cabbage (Brassica rapa). The database has three major parts and functions; (1) Map Search, (2) Marker Search, and (3) QTL Search. At present, the database provides characterization information for about 3258 genetic markers. It has 2800 RFLP and 112 QTL markers related to rice in addition to 321 RFLP and 25 PCR-based markers for Chinese cabbage. In addition, a genetic linkage map was also constructed by using 1,054 markers from 2,912 markers in rice.

Availability

The database is available for free at http://www.niab.go.kr/nabic/PlantGM     

Fine mapping and candidate gene analysis of the novel thermo-sensitive genic male sterility tms9-1 gene in rice

Key message

Fine mapping of the novel thermo-sensitive genic male sterility locus tms9 - 1 in the traditional TGMS line HengnongS-1 revealed that the MALE STERILITY1 homolog OsMS1 is the candidate gene.

Abstract

Photoperiod-thermo-sensitive genic male sterility (P/TGMS) has been widely used in the two-line hybrid rice breeding system. HengnongS-1 is one of the oldest TGMS lines and is often used in indica two-line breeding programs in China. In this study, our genetic analysis showed that the TGMS gene in HengnongS-1 was controlled by a single recessive gene that was non-allelic with the other TGMS loci identified, including C815S, Zhu1S and Y58S. Using SSR markers and bulked segregant analysis, we located the TGMS locus on chromosome 9 and named the gene tms9-1. Fine mapping further narrowed the tms9-1 loci to a 162 kb interval between two dCAPS markers. Sequence analysis revealed that a T to C substitution results in an amino acid change in the tms9-1 candidate gene (Os09g27620) in HengnongS-1 as compared to Minghui63. Sequencing of other rice accessions, including six P/TGMS lines, seven indica varieties and nine japonica varieties, showed that this SNP was exclusive to HengnongS-1. With multiple sequence alignment and expression pattern analyses, the rice MALE STERILITY1 homolog OsMS1 gene was identified as the candidate gene for tms9-1. Therefore, our study identified a novel TGMS locus and will facilitate the functional identification of the tms9-1 gene. Moreover, the markers linked to the tms9-1 gene will provide useful tools for the development of new TGMS lines by marker-assisted selection in two-line hybrid rice breeding programs.