Arrow Plot: a new graphical tool for selecting up and down regulated genes and genes differentially expressed on samples subgroups

Background: A common task in analyzing microarray data is to determine which genes are differentially expressed across two (or more) kind of tissue samples or samples submitted under experimental conditions. Several statistical methods have been proposed to accomplish this goal, generally based on measures of distance between classes. It is well known that biological samples are heterogeneous because of factors such as molecular subtypes or genetic background that are often unknown to the experimenter. For instance, in experiments which involve molecular classification of tumors it is important to identify significant subtypes of cancer. Bimodal or multimodal distributions often reflect the presence of subsamples mixtures. Consequently, there can be genes differentially expressed on sample subgroups which are missed if usual statistical approaches are used. In this paper we propose a new graphical tool which not only identifies genes with up and down regulations, but also genes with differential expression in different subclasses, that are usually missed if current statistical methods are used. This tool is based on two measures of distance between samples, namely the overlapping coefficient (OVL) between two densities and the area under the receiver operating characteristic (ROC) curve. The methodology proposed here was implemented in the open-source R software. Results: This method was applied to a publicly available dataset, as well as to a simulated dataset. We compared our results with the ones obtained using some of the standard methods for detecting differentially expressed genes, namely Welch t-statistic, fold change (FC), rank products (RP), average difference (AD), weighted average difference (WAD), moderated t-statistic (modT), intensity-based moderated t-statistic (ibmT), significance analysis of microarrays (samT) and area under the ROC curve (AUC). On both datasets all differentially expressed genes with bimodal or multimodal distributions were not selected by all standard selection procedures. We also compared our results with (i) area between ROC curve and rising area (ABCR) and (ii) the test for not proper ROC curves (TNRC). We found our methodology more comprehensive, because it detects both bimodal and multimodal distributions and different variances can be considered on both samples. Another advantage of our method is that we can analyze graphically the behavior of different kinds of differentially expressed genes. Conclusion: Our results indicate that the arrow plot represents a new flexible and useful tool for the analysis of gene expression profiles from microarrays.

Differentially expressed gene profile in the 6-hydroxy-dopamine-induced cell culture model of Parkinson's disease.

Parkinson's disease (PD) is a chronic neurodegenerative disorder characterized by progressive loss of dopaminergic (DA) neurons of the substantia nigra pars compacta with unknown aetiology. 6-Hydroxydopamine (6-OHDA) treatment of neuronal cells is an established in vivo model for mimicking the effect of oxidative stress found in PD brains. We examined the effects of 6-OHDA treatment on human neuroblastoma cells (SH-SY5Y) and primary mesencephalic cultures. Using a reverse arbitrarily primed polymerase chain reaction (RAP-PCR) approach we generated reproducible genetic fingerprints of differential expression levels in cell cultures treated with 6-OHDA. Of the resulting sequences, 23 showed considerable homology to known human coding sequences. The results of the RAP-PCR were validated by reverse transcription PCR, real-time PCR and, for selected genes, by Western blot analysis and immunofluorescence. In four cases, [tomoregulin-1 (TMEFF-1), collapsin response mediator protein 1 (CRMP-1), neurexin-1, and phosphoribosylaminoimidazole synthetase (GART)], a down-regulation of mRNA and protein levels was detected. Further studies will be necessary on the physiological role of the identified proteins and their impact on pathways leading to neurodegeneration in PD.

Programmed Cell Death in Procyclic Form Trypanosoma brucei rhodesiense - Identification of Differentially Expressed Genes during Con A Induced Death

Trypanosoma brucei rhodesiense can be induced to undergo apoptosis after stimulation with Con A. As cell death in these parasites is associated with de novo gene expression we have applied a differential display technique, Randomly Amplified Differential Expressed Sequence-Polymerase Chain Reaction (RADES-PCR) to the study of gene expression during Con A induced cell death in these organisms. Twenty-two differentially displayed products have been cloned and sequenced. These represent the first endogenous genes to be identified as implicated in cellular death in trypanosomatids (the most primitive eukaryote in which apoptosis has been described). Evidence for an ancestral death machinery, `proto-apoptosis' in single celled organisms is discussed.

Characterization of constitutive and putative differentially expressed mRNAs by means of expressed sequence tags, differential display reverse transcriptase-PCR and randomly amplified polymorphic DNA-PCR from the sand fly vector Lutzomyia longipalpis

Molecular studies of insect disease vectors are of paramount importance for understanding parasite-vector relationship. Advances in this area have led to important findings regarding changes in vectors' physiology upon blood feeding and parasite infection. Mechanisms for interfering with the vectorial capacity of insects responsible for the transmission of diseases such as malaria, Chagas disease and dengue fever are being devised with the ultimate goal of developing transgenic insects. A primary necessity for this goal is information on gene expression and control in the target insect. Our group is investigating molecular aspects of the interaction between Leishmania parasites and Lutzomyia sand flies. As an initial step in our studies we have used random sequencing of cDNA clones from two expression libraries made from head/thorax and abdomen of sugar fed L. longipalpis for the identification of expressed sequence tags (EST). We applied differential display reverse transcriptase-PCR and randomly amplified polymorphic DNA-PCR to characterize differentially expressed mRNA from sugar and blood fed insects, and, in one case, from a L. (V.) braziliensis-infected L. longipalpis. We identified 37 cDNAs that have shown homology to known sequences from GeneBank. Of these, 32 cDNAs code for constitutive proteins such as zinc finger protein, glutamine synthetase, G binding protein, ubiquitin conjugating enzyme. Three are putative differentially expressed cDNAs from blood fed and Leishmania-infected midgut, a chitinase, a V-ATPase and a MAP kinase. Finally, two sequences are homologous to Drosophila melanogaster gene products recently discovered through the Drosophila genome initiative.

Identification of a differentially expressed mRNA in axenic Leishmania panamensis amastigotes

Differential display technique was applied in order to identify transcripts which are present in axenic amastigotes but not in promastigotes of the Leishmania panamensis parasites. One of them was cloned and the sequence reveals an open reading frame of 364 amino acids (aprox. 40 kDa). The deduced protein is homologous to the serine/threonine protein kinases and specially to the mitogen activates protein kinases from eukaryotic species. Southern blot analysis suggest that this transcript, named lpmkh, is present in the genome of the parasite as a single copy gene. These results could imply that lpmkh could be involved in the differentiation process or the preservation of amastigotes in axenic conditions.

MicroRNA Profile of Circulating CD4-positive Regulatory T Cells in Human Adults and Impact of Differentially Expressed MicroRNAs on Expression of Two Genes Essential to Their Function.

Regulatory T cells (Tregs) are characterized by a high expression of IL-2 receptor α chain (CD25) and of forkhead box P3 (FOXP3), the latter being essential for their development and function. Another major player in the regulatory function is the cytotoxic T-lymphocyte associated molecule-4 (CTLA-4) that inhibits cytotoxic responses. However, the regulation of CTLA-4 expression remains less well explored. We therefore studied the microRNA signature of circulating CD4(+) Tregs isolated from adult healthy donors and identified a signature composed of 15 differentially expressed microRNAs. Among those, miR-24, miR-145, and miR-210 were down-regulated in Tregs compared with controls and were found to have potential target sites in the 3'-UTR of FOXP3 and CTLA-4; miR-24 and miR-210 negatively regulated FOXP3 expression by directly binding to their two target sites in its 3'-UTR. On the other hand, miR-95, which is highly expressed in adult peripheral blood Tregs, positively regulated FOXP3 expression via an indirect mechanism yet to be identified. Finally, we showed that miR-145 negatively regulated CTLA-4 expression in human CD4(+) adult peripheral blood Tregs by binding to its target site in CTLA-4 transcript 3'-UTR. To our knowledge, this is the first identification of a human adult peripheral blood CD4(+) Treg microRNA signature. Moreover, unveiling one mechanism regulating CTLA-4 expression is novel and may lead to a better understanding of the regulation of this crucial gene.

FTZ-F1-related orphan receptors in Xenopus laevis: transcriptional regulators differentially expressed during early embryogenesis.

Orphan receptors of the FTZ-F1-related group of nuclear receptors (xFF1r) were identified in Xenopus laevis by isolation of cDNAs from a neurula stage library. Two cDNAs were found, which encode full length, highly related receptor proteins, xFF1rA and B, whose closet relative known so far is the murine LRH-1 orphan receptor. xFF1rA protein expressed by a recombinant vaccinia virus system specifically binds to FTZ-F1 response elements (FRE; PyCAAGGPyCPu). In cotransfection studies, xFF1rA constitutively activates transcription, in a manner dependent on the number of FREs. The amounts of at least four mRNAs encoding full-length receptors greatly increase between gastrula and early tailbud stages and decrease at later stages. At early tailbud stages, xFTZ-F1-related antigens are found in all nuclei of the embryo.

Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues

MicroRNAs (miRNAs) are short non-coding RNA molecules playing regulatory roles by repressing translation or cleaving RNA transcripts. Although the number of verified human miRNA is still expanding, only few have been functionally described. However, emerging evidences suggest the potential involvement of altered regulation of miRNA in pathogenesis of cancers and these genes are thought to function as both tumours suppressor and oncogenes. In our study, we examined by Real-Time PCR the expression of 156 mature miRNA in colorectal cancer. The analysis by several bioinformatics algorithms of colorectal tumours and adjacent non-neoplastic tissues from patients and colorectal cancer cell lines allowed identifying a group of 13 miRNA whose expression is significantly altered in this tumor. The most significantly deregulated miRNA being miR-31, miR-96, miR-133b, miR-135b, miR-145, and miR-183. In addition, the expression level of miR-31 was correlated with the stage of CRC tumor. Our results suggest that miRNA expression profile could have relevance to the biological and clinical behavior of colorectal neoplasia.

Networking of differentially expressed genes in human cancer cells resistant to methotrexate

BACKGROUND: The need for an integrated view of data obtained from high-throughput technologies gave rise to network analyses. These are especially useful to rationalize how external perturbations propagate through the expression of genes. To address this issue in the case of drug resistance, we constructed biological association networks of genes differentially expressed in cell lines resistant to methotrexate (MTX). METHODS: Seven cell lines representative of different types of cancer, including colon cancer (HT29 and Caco2), breast cancer (MCF-7 and MDA-MB-468), pancreatic cancer (MIA PaCa-2), erythroblastic leukemia (K562) and osteosarcoma (Saos-2), were used. The differential expression pattern between sensitive and MTX-resistant cells was determined by whole human genome microarrays and analyzed with the GeneSpring GX software package. Genes deregulated in common between the different cancer cell lines served to generate biological association networks using the Pathway Architect software. RESULTS: Dikkopf homolog-1 (DKK1) is a highly interconnected node in the network generated with genes in common between the two colon cancer cell lines, and functional validations of this target using small interfering RNAs (siRNAs) showed a chemosensitization toward MTX. Members of the UDP-glucuronosyltransferase 1A (UGT1A) family formed a network of genes differentially expressed in the two breast cancer cell lines. siRNA treatment against UGT1A also showed an increase in MTX sensitivity. Eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) was overexpressed among the pancreatic cancer, leukemia and osteosarcoma cell lines, and siRNA treatment against EEF1A1 produced a chemosensitization toward MTX. CONCLUSIONS: Biological association networks identified DKK1, UGT1As and EEF1A1 as important gene nodes in MTX-resistance. Treatments using siRNA technology against these three genes showed chemosensitization toward MTX.

Networking of differentially expressed genes in human cancer cells resistant to methotrexate

BACKGROUND: The need for an integrated view of data obtained from high-throughput technologies gave rise to network analyses. These are especially useful to rationalize how external perturbations propagate through the expression of genes. To address this issue in the case of drug resistance, we constructed biological association networks of genes differentially expressed in cell lines resistant to methotrexate (MTX). METHODS: Seven cell lines representative of different types of cancer, including colon cancer (HT29 and Caco2), breast cancer (MCF-7 and MDA-MB-468), pancreatic cancer (MIA PaCa-2), erythroblastic leukemia (K562) and osteosarcoma (Saos-2), were used. The differential expression pattern between sensitive and MTX-resistant cells was determined by whole human genome microarrays and analyzed with the GeneSpring GX software package. Genes deregulated in common between the different cancer cell lines served to generate biological association networks using the Pathway Architect software. RESULTS: Dikkopf homolog-1 (DKK1) is a highly interconnected node in the network generated with genes in common between the two colon cancer cell lines, and functional validations of this target using small interfering RNAs (siRNAs) showed a chemosensitization toward MTX. Members of the UDP-glucuronosyltransferase 1A (UGT1A) family formed a network of genes differentially expressed in the two breast cancer cell lines. siRNA treatment against UGT1A also showed an increase in MTX sensitivity. Eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) was overexpressed among the pancreatic cancer, leukemia and osteosarcoma cell lines, and siRNA treatment against EEF1A1 produced a chemosensitization toward MTX. CONCLUSIONS: Biological association networks identified DKK1, UGT1As and EEF1A1 as important gene nodes in MTX-resistance. Treatments using siRNA technology against these three genes showed chemosensitization toward MTX.

Hunting for differentially expressed genes

Differentially expressed genes are usually identified by comparing steady-state mRNA concentrations. Several methods have been used for this purpose, including differential hybridization, cDNA subtraction, differential display and, more recently, DNA chips. Subtractive hybridization has significantly improved after the polymerase chain reaction was incorporated into the original method and many new protocols have been established. Recently, the availability of the well-known coding sequences for some organisms has greatly facilitated gene expression analysis using high-density microarrays. Here, we describe some of these modifications and discuss the benefits and drawbacks of the various methods corresponding to the main advances in this field.

The DNA puff BhB10-1 gene is differentially expressed in various tissues of Bradysia hygida late larvae and constitutively transcribed in transgenic Drosophila

We extended the characterization of the DNA puff BhB10-1 gene of Bradysia hygida by showing that, although its mRNA is detected only at the end of the fourth larval instar, BhB10-1 expression is not restricted to the salivary gland, the tissue in which this gene is amplified. Different amounts of BhB10-1 mRNA were detected in other larval tissues such as gut, Malpighian tubules, fat body, brain and cuticle, suggesting that this gene is expressed differentially in the various tissues analyzed. Analysis of transgenic Drosophila carrying the BhB10-1 transcription unit and flanking sequences revealed that the tested fragment promotes transcription in a constitutive manner. We suggest that either cis-regulatory elements are missing in the transgene or factors that temporally regulate the BhB10-1 gene in B. hygida are not conserved in Drosophila.

Global transcriptome analysis and identification of differentially expressed genes after infection of Fragaria vesca with powdery mildew (Podosphaera aphanis)

Background: Podosphaera aphanis, the causal agent of strawberry powdery mildew causes significant economic loss worldwide. Methods: We used the diploid strawberry species Fragaria vesca as a model to study plant pathogen interactions. RNA-seq was employed to generate a transcriptome dataset from two accessions, F. vesca ssp. vesca Hawaii 4 (HW) and F. vesca f. semperflorens Yellow Wonder 5AF7 (YW) at 1 d (1 DAI) and 8 d (8 DAI) after infection. Results: Of the total reads identified about 999 million (92%) mapped to the F. vesca genome. These transcripts were derived from a total of 23,470 and 23,464 genes in HW and YW, respectively from the three time points (control, 1 and 8 DAI). Analysis identified 1,567, 1,846 and 1,145 up-regulated genes between control and 1 DAI, control and 8 DAI, and 1 and 8 DAI, respectively in HW. Similarly, 1,336, 1,619 and 968 genes were up-regulated in YW. Also 646, 1,098 and 624 down-regulated genes were identified in HW, while 571, 754 and 627 genes were down-regulated in YW between all three time points, respectively. Conclusion: Investigation of differentially expressed genes (log2 fold changes �5) between control and 1 DAI in both HW and YW identified a large number of genes related to secondary metabolism, signal transduction; transcriptional regulation and disease resistance were highly expressed. These included flavonoid 3´-monooxygenase, peroxidase 15, glucan endo-1,3-β-glucosidase 2, receptor-like kinases, transcription factors, germin-like proteins, F-box proteins, NB-ARC and NBS-LRR proteins. This is the first application of RNA-seq to any pathogen interaction in strawberry

Screening the Schistosoma mansoni transcriptome for genes differentially expressed in the schistosomulum stage in search for vaccine candidates

Schistosomiasis affects more than 200 million people worldwide; another 600 million are at risk of infection. The schistosomulum stage is believed to be the target of protective immunity in the attenuated cercaria vaccine model. In an attempt to identify genes up-regulated in the schistosomulum stage in relation to cercaria, we explored the Schistosoma mansoni transcriptome by looking at the relative frequency of reads in EST libraries from both stages. The 400 genes potentially up-regulated in schistosomula were analyzed as to their Gene Ontology categorization, and we have focused on those encoding-predicted proteins with no similarity to proteins of other organisms, assuming they could be parasite-specific proteins important for survival in the host. Up-regulation in schistosomulum relative to cercaria was validated with real-time reverse transcription polymerase chain reaction (RT-PCR) for five out of nine selected genes (56%). We tested their protective potential in mice through immunization with DNA vaccines followed by a parasite challenge. Worm burden reductions of 16-17% were observed for one of them, indicating its protective potential. Our results demonstrate the value and caveats of using stage-associated frequency of ESTs as an indication of differential expression coupled to DNA vaccine screening in the identification of novel proteins to be further investigated as potential vaccine candidates.

Sugarcane genes differentially expressed during water deficit

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)