LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over-or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

Differential gene expression between the biotrophic-like and saprotrophic mycelia of the witches` broom pathogen Moniliophthora perniciosa

Moniliophthora perniciosa is a hemibiotrophic fungus that causes witches` broom disease (WBD) in cacao. Marked dimorphism characterizes this fungus, showing a monokaryotic or biotrophic phase that causes disease symptoms and a later dikaryotic or saprotrophic phase. A combined strategy of DNA microarray, expressed sequence tag, and real-time reverse-transcriptase polymerase chain reaction analyses was employed to analyze differences between these two fungal stages in vitro. In all, 1,131 putative genes were hybridized with cDNA from different phases, resulting in 189 differentially expressed genes, and 4,595 reads were clusterized, producing 1,534 unigenes. The analysis of these genes, which represent approximately 21% of the total genes, indicates that the biotrophic-like phase undergoes carbon and nitrogen catabollite repression that correlates to the expression of phytopathogenicity genes. Moreover, downregulation of mitochondrial oxidative phosphorylation and the presence of a putative ngr1 of Saccharomyces cerevisiae could help explain its lower growth rate. In contrast, the saprotrophic mycelium expresses genes related to the metabolism of hexoses, ammonia, and oxidative phosphorylation, which could explain its faster growth. Antifungal toxins were upregulated and could prevent the colonization by competing fungi. This work significantly contributes to our understanding of the molecular mechanisms of WBD and, to our knowledge, is the first to analyze differential gene expression of the different phases of a hemibiotrophic fungus.

Cloning, characterization and functional analysis of a 1-FEH cDNA from Vernonia herbacea (Vell.) Rusby

Variations in the inulin contents have been detected in rhizophores of Vernonia herbacea during the phenological cycle. These variations indicate the occurrence of active inulin synthesis and depolymerization throughout the cycle and a role for this carbohydrate as a reserve compound. 1-Fructan exohydrolase (1-FEH) is the enzyme responsible for inulin depolymerization, and its activity has been detected in rhizophores of sprouting plants. Defoliation and low temperature are enhancer conditions of this 1-FEH activity. The aim of the present work was the cloning of this enzyme. Rhizophores were collected from plants induced to sprout, followed by storage at 5C. A full length 1-FEH cDNA sequence was obtained by PCR and inverse PCR techniques, and expressed in Pichia pastoris. Cold storage enhances FEH gene expression. Vh1-FEH was shown to be a functional 1-FEH, hydrolyzing predominantly -2,1 linkages, sharing high identity with chicory FEH sequences, and its activity was inhibited by 81 in the presence of 10 mM sucrose. In V. herbacea, low temperature and sucrose play a role in the control of fructan degradation. This is the first study concerning the cloning and functional analysis of a 1-FEH cDNA of a native species from the Brazilian Cerrado. Results will contribute to understanding the role of fructans in the establishment of a very successful fructan flora of the Brazilian Cerrado, subjected to water limitation and low temperature during winter.

Antiviral and antiparasite properties of an L-amino acid oxidase from the Snake Bothrops jararaca: Cloning and identification of a complete cDNA sequence (Retracted article. See vol. 80, pg. 288, 2010)

L-Amino acid oxidases (LAAOs, EC 1.4.3.2) are flavoenzymes that catalyze the stereospecific oxidative deamination of an L-amino acid substrate to the corresponding a-ketoacid with hydrogen peroxide and ammonia production. The present work describes the first report on the antiviral (Dengue virus) and antiprotozoal (trypanocidal and leishmanicide) activities of a Bothrops jararaca L-amino acid oxidase (BjarLAAO-I) and identify its cDNA sequence. Antiparasite effects were inhibited by catalase, suggesting that they are mediated by H(2)O(2) production. Cells infected with DENV-3 virus previously treated with BjarLAAO-I, showed a decrease in viral titer (13-83-fold) when compared with cells infected with untreated viruses. Untreated and treated promastigotes (T. cruzi and L. amazonensis) were observed by transmission electron microscopy with different degrees of damage. Its complete cDNA sequence, with 1452 bp, encoded an open reading frame of 484 amino acid residues with a theoretical molecular weight and pl of 54,771.8 and 5.7, respectively. The cDNA-deduced amino acid sequence of BjarLAAO shows high identity to LAAOs from other snake venoms. Further investigations will be focused on the related molecular and functional correlation of these enzymes. Such a study should provide valuable information for the therapeutic development of new generations of microbicidal drugs. (C) 2008 Elsevier Inc. All rights reserved.

Effect of the immunosuppressants on hepatocyte proliferation and apoptosis in a young animal model of liver regeneration: An immunohistochemical study using tissue microarrays

Hepatocyte proliferation and apoptosis (programmed cell death) occur during the liver parenchyma regeneration and the liver size modeling is mainly controlled by hepatocyte apoptosis. The purpose of the present study was to verify the influence of immunosuppressant drugs on these phenomena by utilizing tissue microarray techniques. Thirty-six weaning rats (age 21-23 days, weight 30-50 g) were divided into six groups: control, sham, hepatectomy, hepatectomy plus solumedrol, hepatectomy plus CsA, and hepatectomy plus Tac. The animals were killed one day after hepatectomy, and the remnant livers were weighed and harvested for tissue microarray sections. Liver cell proliferation was evaluated by staining for PCNA and apoptosis was detected by the TUNEL method. It was verified that CsA promoted a decrease in the liver weight, Tac and CsA decreased the proliferation index of hepatocytes, and glucocorticoid had no significant effects. The apoptosis index was not altered by hepatectomy or immunosuppressants. Our data indicate that, in the growing rat, CsA and Tac have negative effects on hepatocyte proliferation and have no effect on the hepatocyte apoptosis.

cDNA microarray analysis of differentially expressed genes in penile tissue after treatment with tadalafil

To evaluate differential gene expression in penile tissue after treatment with the phosphodiesterase 5 (PDE5) inhibitor tadalafil, as of the three clinically available PDE5 inhibitors (sildenafil, tadalafil, and vardenafil) used for the treatment of erectile dysfunction (ED), tadalafil has a long half-life and low incidence of side-effects. In all, 32 adult rats were divided into two groups. The control group received 0.5 mL of drinking water alone, while the tadalafil group was treated with tadalafil at a dose of 0.27 mg/kg. At 4 h after treatment with water or tadalafil the rats were killed and the penile tissue was removed. The total RNA was isolated from the penile tissue from both groups and differentially expressed genes were identified by cDNA microarray analysis. To validate the expression data from the microarray analysis, quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry were used. In all, 153 genes were differentially expressed between the control group and the tadalafil group. We validated the microarray results by quantitative PCR for the insulin-like growth factor binding protein 6 (IGFBP-6) gene and the neuronal calcium sensor 1 (NCS-1) gene, both of which were up-regulated in the tadalafil group, and for the natriuretic peptide receptor 1 (NPR-1) gene that was down-regulated in this group. Immunohistochemistry showed localization of the NCS-1 protein in sinusoid trabeculae of the corpus cavernosum in control and tadalafil-treated rats. There was differential expression in 153 genes after tadalafil treatment. Some of these genes such as IGFBP-6, NPR-1 and NCS-1, might result in new targets in the treatment of ED.

Molecular cloning, sequencing, and expression analysis of presenilin cDNA from Schistosoma mansoni

Presenilins (PS) are integral membrane proteins involved, among other functions, in regulated intramembrane proteolysis. In this study, we report the identification and characterization of a complementary DNA from Schistosoma mansoni exhibiting a significant homology to human and nonvertebrate presinilins. S. mansoni contained a 1,485 bp open reading frame encoding a predicted protein of 494 amino acids. Alignment of predicted amino acid sequence of S. mansoni with PS (SmPS) from other species revealed up to 40% similarity shared among the investigated organisms. In addition, phylogenetic analyses demonstrated SmPS being closely related to its orthologues found in Schistosoma japonicum and Caenorhabditis elegans. Expression analysis of SmPS using quantitative real-time PCR revealed that the transcript is up-regulated in the egg stage. We hypothesize that the high level of SmPS in the S. mansoni embryo correlates to an important role during cellular signaling associated to larval development. To our knowledge, this study represents the first attempt to investigate the existence and abundance of PS from a helminth parasite.

Magnetic purification of biotinylated cDNA removes false priming and ensures strand-specificity of RT-PCR for enteroviral RNAs

The detection of replicative intermediate RNAs as markers of active replication of RNA viruses is an essential tool to investigate pathogenesis in acute viral infections, as well as in their long-term sequelae. In this regard, strand-specific PCR has been used widely to distinguish (-) and (+) enteroviral RNAs in pathogenesis studies of diseases such as dilated cardiomyopathy. It has been generally assumed that oligonucleotide-primed reverse transcription of a given RNA generates only the corresponding specific cDNA, thus assuring the specificity of a PCR product amplified from it. Nevertheless, such assumed strand-specificity is a fallacy, because falsely primed cDNAs can be produced by RNA reverse transcription in the absence of exogenously added primers, (cDNA(primer)(-)), and such falsely primed cDNAs are amplifiable by PCR in the same way as the correctly primed cDNAs. Using as a prototype the coxsackievirus B5 (CVB5), a (+) strand RNA virus, it was shown that cDNA(primer)(-) renders the differential detection of viral (-) and (+) RNAs by conventional PCR virtually impossible, due to gross non-specificity. Using in vitro transcribed CVB5 RNAs (+) and (-), it was shown that cDNA(primer)(-) could be removed effectively by magnetic physical separation of correctly primed biotinylated cDNA. Such strategy enabled truly strand-specific detection of RNA (-) and (+), not only for CVB5, but also for other non-polio enteroviruses. These findings indicate that previous conclusions supporting a role for the persistence of actively replicating enterovirus in the pathogenesis of chronic myocarditis should be regarded with strong skepticism and purification of correctly primed cDNA should be used for strand-specific PCR of viral RNA in order to obtain reliable information on this important subject. (C) 2009 Elsevier B.V. All rights reserved.

Transcriptional Response of Peripheral Lymphocytes to Early Fibrosarcoma: A Model System for Cancer Detection Based on Hybridization Signatures

Since circulating leukocytes, mainly B and T cells, continuously maintain vigilant and comprehensive immune surveillance, these cells could be used as reporters for signs of infection or other pathologies, including cancer. Activated lymphocyte clones trigger a sensitive transcriptional response, which could be identified by gene expression profiling. To assess this hypothesis, we conducted microarray analysis of the gene expression profile of lymphocytes isolated from immunocompetent BALB/c mice subcutaneously injected with different numbers of tumorigenic B61 fibrosarcoma cells. Flow cytometry demonstrated that the number of circulating T (CD3(+)CD4(+) or CD3(+)CD8(+)) or B (CD19(+)) cells did not change. However, the lymphocytes isolated from tumor cell-injected animals expressed a unique transcriptional profile that was identifiable before the development of a palpable tumor mass. This finding demonstrates that the transcriptional response appears before alterations in the main lymphocyte subsets and that the gene expression profile of peripheral lymphocytes can serve as a sensitive and accurate method for the early detection of cancer. Exp Biol Med 234:802-812, 2009

Production of DNA microarray and expression analysis of genes from Xylella fastidiosa in different culture media

DNA Microarray was developed to monitor the expression of many genes from Xylella fastidiosa, allowing the side by-side comparison of two situations in a single experiment. The experiments were performed using X. fastidiosa cells grown in two culture media: BCYE and XDM2. The primers were synthesized, spotted onto glass slides and the array was hybridized against fluorescently labeled cDNAs. The emitted signals were quantified, normalized and the data were statistically analyzed to verify the differentially expressed genes. According to the data, 104 genes were differentially expressed in XDM2 and 30 genes in BCYE media. The present study showed that DNA microarray technique efficiently differentiate the expressed genes under different conditions.

Cloning, purification and comparative characterization of two digestive lysozymes from Musca domestica larvae

cDNA coding for two digestive lysozymes (MdL1 and MdL2) of the Musca domestica housefly was cloned and sequenced. MdL2 is a novel minor lysozyme, whereas MdL1 is the major lysozyme thus far purified from M. domestica midgut. MdL1 and MdL2 were expressed as recombinant proteins in Pichia pastoris, purified and characterized. The lytic activities of MdL1 and MdL2 upon Micrococcus lysodeikticus have an acidic pH optimum (4.8) at low ionic strength (μ = 0.02), which shifts towards an even more acidic value, pH 3.8, at a high ionic strength (μ = 0.2). However, the pH optimum of their activities upon 4-methylumbelliferyl N-acetylchitotrioside (4.9) is not affected by ionic strength. These results suggest that the acidic pH optimum is an intrinsic property of MdL1 and MdL2, whereas pH optimum shifts are an effect of the ionic strength on the negatively charged bacterial wall. MdL2 affinity for bacterial cell wall is lower than that of MdL1. Differences in isoelectric point (pI) indicate that MdL2 (pI = 6.7) is less positively charged than MdL1 (pI = 7.7) at their pH optima, which suggests that electrostatic interactions might be involved in substrate binding. In agreement with that finding, MdL1 and MdL2 affinities for bacterial cell wall decrease as ionic strength increases.

The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Searching for molecular markers in head and neck squamous cell carcinomas (HNSCC) by statistical and bioinformatic analysis of larynx-derived SAGE libraries

Background: Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancies in humans. The average 5-year survival rate is one of the lowest among aggressive cancers, showing no significant improvement in recent years. When detected early, HNSCC has a good prognosis, but most patients present metastatic disease at the time of diagnosis, which significantly reduces survival rate. Despite extensive research, no molecular markers are currently available for diagnostic or prognostic purposes. Methods: Aiming to identify differentially-expressed genes involved in laryngeal squamous cell carcinoma (LSCC) development and progression, we generated individual Serial Analysis of Gene Expression (SAGE) libraries from a metastatic and non-metastatic larynx carcinoma, as well as from a normal larynx mucosa sample. Approximately 54,000 unique tags were sequenced in three libraries. Results: Statistical data analysis identified a subset of 1,216 differentially expressed tags between tumor and normal libraries, and 894 differentially expressed tags between metastatic and non-metastatic carcinomas. Three genes displaying differential regulation, one down-regulated (KRT31) and two up-regulated (BST2, MFAP2), as well as one with a non-significant differential expression pattern (GNA15) in our SAGE data were selected for real-time polymerase chain reaction (PCR) in a set of HNSCC samples. Consistent with our statistical analysis, quantitative PCR confirmed the upregulation of BST2 and MFAP2 and the downregulation of KRT31 when samples of HNSCC were compared to tumor-free surgical margins. As expected, GNA15 presented a non-significant differential expression pattern when tumor samples were compared to normal tissues. Conclusion: To the best of our knowledge, this is the first study reporting SAGE data in head and neck squamous cell tumors. Statistical analysis was effective in identifying differentially expressed genes reportedly involved in cancer development. The differential expression of a subset of genes was confirmed in additional larynx carcinoma samples and in carcinomas from a distinct head and neck subsite. This result suggests the existence of potential common biomarkers for prognosis and targeted-therapy development in this heterogeneous type of tumor.

The expression of genes coding for distinct types of glycine-rich proteins varies according to the biology of three metastriate ticks, Rhipicephalus (Boophilus) microplus, Rhipicephalus sanguineus and Amblyomma cajennense

Background: Ticks secrete a cement cone composed of many salivary proteins, some of which are rich in the amino acid glycine in order to attach to their hosts' skin. Glycine-rich proteins (GRPs) are a large family of heterogeneous proteins that have different functions and features; noteworthy are their adhesive and tensile characteristics. These properties may be essential for successful attachment of the metastriate ticks to the host and the prolonged feeding necessary for engorgement. In this work, we analyzed Expressed Sequence Tags (ESTs) similar to GRPs from cDNA libraries constructed from salivary glands of adult female ticks representing three hard, metastriate species in order to verify if their expression correlated with biological differences such as the numbers of hosts ticks feed on during their parasitic life cycle, whether one (monoxenous parasite) or two or more (heteroxenous parasite), and the anatomy of their mouthparts, whether short (Brevirostrata) or long (Longirostrata). These ticks were the monoxenous Brevirostrata tick, Rhipicephalus (Boophilus) microplus, a heteroxenous Brevirostrata tick, Rhipicephalus sanguineus, and a heteroxenous Longirostrata tick, Amblyomma cajennense. To further investigate this relationship, we conducted phylogenetic analyses using sequences of GRPs from these ticks as well as from other species of Brevirostrata and Longirostrata ticks. Results: cDNA libraries from salivary glands of the monoxenous tick, R. microplus, contained more contigs of glycine-rich proteins than the two representatives of heteroxenous ticks, R. sanguineus and A. cajennense (33 versus, respectively, 16 and 11). Transcripts of ESTs encoding GRPs were significantly more numerous in the salivary glands of the two Brevirostrata species when compared to the number of transcripts in the Longirostrata tick. The salivary gland libraries from Brevirostrata ticks contained numerous contigs significantly similar to silks of true spiders (17 and 8 in, respectively, R. microplus and R. sanguineus), whereas the Longirostrata tick contained only 4 contigs. The phylogenetic analyses of GRPs from various species of ticks showed that distinct clades encoding proteins with different biochemical properties are represented among species according to their biology. Conclusions: We found that different species of ticks rely on different types and amounts of GRPs in order to attach and feed on their hosts. Metastriate ticks with short mouthparts express more transcripts of GRPs than a tick with long mouthparts and the tick that feeds on a single host during its life cycle contain a greater variety of these proteins than ticks that feed on several hosts.

Gene expression profile of the plant pathogen Fusarium graminearum under the antagonistic effect of Pantoea agglomerans

The pathogenic fungus Fusarium graminearum is an ongoing threat to agriculture, causing losses in grain yield and quality in diverse crops. Substantial progress has been made in the identification of genes involved in the suppression of phytopathogens by antagonistic microorganisms; however, limited information regarding responses of plant pathogens to these biocontrol agents is available. Gene expression analysis was used to identify differentially expressed transcripts of the fungal plant pathogen F. graminearum under antagonistic effect of the bacterium Pantoea agglomerans. A macroarray was constructed, using 1014 transcripts from an F. graminearum cDNA library. Probes consisted of the cDNA of F. graminearum grown in the presence and in the absence of P. agglomerans. Twenty-nine genes were either up (19) or down (10) regulated during interaction with the antagonist bacterium. Genes encoding proteins associated with fungal defense and/or virulence or with nutritional and oxidative stress responses were induced. The repressed genes coded for a zinc finger protein associated with cell division, proteins containing cellular signaling domains, respiratory chain proteins, and chaperone-type proteins. These data give molecular and biochemical evidence of response of F. graminearum to an antagonist and could help develop effective biocontrol procedures for pathogenic plant fungi.