ProbFAST: Probabilistic Functional Analysis System Tool

Background: The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis. Results: We proposed a web platform named ProbFAST for MDE analysis which uses Bayesian inference to identify key genes that are intuitively prioritized by means of probabilities. A simulated study revealed that our method gives a better performance when compared to other approaches and when applied to public expression data, we demonstrated its flexibility to obtain relevant genes biologically associated with normal and abnormal biological processes. Conclusions: ProbFAST is a free accessible web-based application that enables MDE analysis on a global scale. It offers an efficient methodological approach for MDE analysis of a set of genes that are turned on and off related to functional information during the evolution of a tumor or tissue differentiation. ProbFAST server can be accessed at http://gdm.fmrp.usp.br/probfast.

Functional characterization of the sciarid BhC4-1 core promoter in transgenic Drosophila

Background: Core promoters are cis-regulatory modules to which bind the basal transcriptional machinery and which participate in the regulation of transcription initiation. Although core promoters have not been extensively investigated through functional assays in a chromosomal context, the available data suggested that the response of a given core promoter might vary depending on the promoter context. Previous studies suggest that a (-57/+40) fragment constitutes the core promoter of the BhC4-1 gene which is located in DNA puff C4 of the sciarid fly Bradysia hygida. Here we tested this (-57/+40) fragment in distinct regulatory contexts in order to verify if promoter context affects its core promoter activity. Results: Consistent with the activity of a core promoter, we showed that in the absence of upstream regulatory sequences the (-57/+40) fragment drives low levels of reporter gene mRNA expression throughout development in transgenic Drosophila. By assaying the (-57/+40) fragment in two distinct regulatory contexts, either downstream of the previously characterized Fbp1 enhancer or downstream of the UAS element, we showed that the BhC4-1 core promoter drives regulated transcription in both the germline and in various tissues throughout development. Furthermore, the use of the BhC4-1 core promoter in a UAS construct significantly reduced salivary gland ectopic expression in third instar larvae, which was previously described to occur in the context of the GAL4/UAS system. Conclusions: Our results from functional analysis in transgenic Drosophila show that the BhC4-1 core promoter drives gene expression regardless of the promoter context that was assayed. New insights into the functioning of the GAL4/UAS system in Drosophila were obtained, indicating that the presence of the SV40 sequence in the 3' UTR of a UAS construct does not preclude expression in the germline. Furthermore, our analysis indicated that ectopic salivary gland expression in the GAL4/UAS system does not depend only on sequences present in the GAL4 construct, but can also be affected by the core promoter sequences in the UAS construct. In this context, we propose that the sciarid BhC4-1 core promoter constitutes a valuable core promoter which can be employed in functional assays in insects.

The gene transformer-2 of Anastrepha fruit flies (Diptera, Tephritidae) and its evolution in insects

Background: In the tephritids Ceratitis, Bactrocera and Anastrepha, the gene transformer provides the memory device for sex determination via its auto-regulation; only in females is functional Tra protein produced. To date, the isolation and characterisation of the gene transformer-2 in the tephritids has only been undertaken in Ceratitis, and it has been shown that its function is required for the female-specific splicing of doublesex and transformer pre-mRNA. It therefore participates in transformer auto-regulatory function. In this work, the characterisation of this gene in eleven tephritid species belonging to the less extensively analysed genus Anastrepha was undertaken in order to throw light on the evolution of transformer-2. Results: The gene transformer-2 produces a protein of 249 amino acids in both sexes, which shows the features of the SR protein family. No significant partially spliced mRNA isoform specific to the male germ line was detected, unlike in Drosophila. It is transcribed in both sexes during development and in adult life, in both the soma and germ line. The injection of Anastrepha transformer-2 dsRNA into Anastrepha embryos caused a change in the splicing pattern of the endogenous transformer and doublesex pre-mRNA of XX females from the female to the male mode. Consequently, these XX females were transformed into pseudomales. The comparison of the eleven Anastrepha Transformer-2 proteins among themselves, and with the Transformer-2 proteins of other insects, suggests the existence of negative selection acting at the protein level to maintain Transformer-2 structural features. Conclusions: These results indicate that transformer-2 is required for sex determination in Anastrepha through its participation in the female-specific splicing of transformer and doublesex pre-mRNAs. It is therefore needed for the auto-regulation of the gene transformer. Thus, the transformer/transfomer-2 > doublesex elements at the bottom of the cascade, and their relationships, probably represent the ancestral state ( which still exists in the Tephritidae, Calliphoridae and Muscidae lineages) of the extant cascade found in the Drosophilidae lineage ( in which tra is just another component of the sex determination gene cascade regulated by Sex-lethal). In the phylogenetic lineage that gave rise to the drosophilids, evolution co-opted for Sex-lethal, modified it, and converted it into the key gene controlling sex determination.

Activation of endogenous p53 by combined p19Arf gene transfer and nutlin-3 drug treatment modalities in the murine cell lines B16 and C6

Background: Reactivation of p53 by either gene transfer or pharmacologic approaches may compensate for loss of p19Arf or excess mdm2 expression, common events in melanoma and glioma. In our previous work, we constructed the pCLPG retroviral vector where transgene expression is controlled by p53 through a p53-responsive promoter. The use of this vector to introduce p19Arf into tumor cells that harbor p53wt should yield viral expression of p19Arf which, in turn, would activate the endogenous p53 and result in enhanced vector expression and tumor suppression. Since nutlin-3 can activate p53 by blocking its interaction with mdm2, we explored the possibility that the combination of p19Arf gene transfer and nutlin-3 drug treatment may provide an additive benefit in stimulating p53 function. Methods: B16 (mouse melanoma) and C6 (rat glioma) cell lines, which harbor p53wt, were transduced with pCLPGp19 and these were additionally treated with nutlin-3 or the DNA damaging agent, doxorubicin. Viral expression was confirmed by Western, Northern and immunofluorescence assays. p53 function was assessed by reporter gene activity provided by a p53-responsive construct. Alterations in proliferation and viability were measured by colony formation, growth curve, cell cycle and MTT assays. In an animal model, B16 cells were treated with the pCLPGp19 virus and/or drugs before subcutaneous injection in C57BL/6 mice, observation of tumor progression and histopathologic analyses. Results: Here we show that the functional activation of endogenous p53wt in B16 was particularly challenging, but accomplished when combined gene transfer and drug treatments were applied, resulting in increased transactivation by p53, marked cell cycle alteration and reduced viability in culture. In an animal model, B16 cells treated with both p19Arf and nutlin-3 yielded increased necrosis and decreased BrdU marking. In comparison, C6 cells were quite susceptible to either treatment, yet p53 was further activated by the combination of p19Arf and nutlin-3. Conclusions: To the best of our knowledge, this is the first study to apply both p19Arf and nutlin-3 for the stimulation of p53 activity. These results support the notion that a p53 responsive vector may prove to be an interesting gene transfer tool, especially when combined with p53- activating agents, for the treatment of tumors that retain wild-type p53.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Functional microarray analysis suggests repressed cell-cell signaling and cell survival-related modules inhibit progression of head and neck squamous cell carcinoma

Background: Cancer shows a great diversity in its clinical behavior which cannot be easily predicted using the currently available clinical or pathological markers. The identification of pathways associated with lymph node metastasis (N+) and recurrent head and neck squamous cell carcinoma (HNSCC) may increase our understanding of the complex biology of this disease. Methods: Tumor samples were obtained from untreated HNSCC patients undergoing surgery. Patients were classified according to pathologic lymph node status (positive or negative) or tumor recurrence (recurrent or non-recurrent tumor) after treatment (surgery with neck dissection followed by radiotherapy). Using microarray gene expression, we screened tumor samples according to modules comprised by genes in the same pathway or functional category. Results: The most frequent alterations were the repression of modules in negative lymph node (N0) and in non-recurrent tumors rather than induction of modules in N+ or in recurrent tumors. N0 tumors showed repression of modules that contain cell survival genes and in non-recurrent tumors cell-cell signaling and extracellular region modules were repressed. Conclusions: The repression of modules that contain cell survival genes in N0 tumors reinforces the important role that apoptosis plays in the regulation of metastasis. In addition, because tumor samples used here were not microdissected, tumor gene expression data are represented together with the stroma, which may reveal signaling between the microenvironment and tumor cells. For instance, in non-recurrent tumors, extracellular region module was repressed, indicating that the stroma and tumor cells may have fewer interactions, which disable metastasis development. Finally, the genes highlighted in our analysis can be implicated in more than one pathway or characteristic, suggesting that therapeutic approaches to prevent tumor progression should target more than one gene or pathway, specially apoptosis and interactions between tumor cells and the stroma.

The Penicillium chrysogenum aclA gene encodes a broad-substrate-specificity acyl-coenzyme A ligase involved in activation of adipic acid, a side-chain precursor for cephem antibiotics

Activation of the cephalosporin side-chain precursor to the corresponding CoA-thioester is an essential step for its incorporation into the P-lactam backbone. To identify an acyl-CoA ligase involved in activation of adipate, we searched in the genome database of Penicillium chrysogenum for putative structural genes encoding acyl-CoA ligases. Chemostat-based transcriptome analysis was used to identify the one presenting the highest expression level when cells were grown in the presence of adipate. Deletion of the gene renamed aclA, led to a 32% decreased specific rate of adipate consumption and a threefold reduction of adipoyl-6-aminopenicillanic acid levels, but did not affect penicillin V production. After overexpression in Escherichia coli, the purified protein was shown to have a broad substrate range including adipate. Finally, protein-fusion with cyan-fluorescent protein showed co-localization with microbody-borne acyl-transferase. Identification and functional characterization of aclA may aid in developing future metabolic engineering strategies for improving the production of different cephalosporins. (C) 2009 Elsevier Inc. All rights reserved.

Extension of the core map of common bean with EST-SSR, RGA, AFLP, and putative functional markers

Microsatellites and gene-derived markers are still underrepresented in the core molecular linkage map of common bean compared to other types of markers. In order to increase the density of the core map, a set of new markers were developed and mapped onto the RIL population derived from the `BAT93` x `Jalo EEP558` cross. The EST-SSR markers were first characterized using a set of 24 bean inbred lines. On average, the polymorphism information content was 0.40 and the mean number of alleles per locus was 2.7. In addition, AFLP and RGA markers based on the NBS-profiling method were developed and a subset of the mapped RGA was sequenced. With the integration of 282 new markers into the common bean core map, we were able to place markers with putative known function in some existing gaps including regions with QTL for resistance to anthracnose and rust. The distribution of the markers over 11 linkage groups is discussed and a newer version of the common bean core linkage map is proposed.

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.

Functional characterization of the Aspergillus fumigatus PHO80 homologue

Phosphate is an ion that is essential for fungal growth. The systems for inorganic phosphate (Pi) acquisition in eukaryotic cells (PHO) have been characterized as a low-affinity (that assures a supply of Pi at normal or high external Pi concentrations) and a high-affinity (activated in response to Pi starvation). Here, as an initial step to understand the PHO pathway in Aspergillus fumigatus, we characterized the PH080 homologue, PhoB(PHO80). We show that the Delta phoB(PHO80) mutant has a polar growth defect (i.e., a delayed germ tube emergence) and, by phenotypic and phosphate uptake analyses, establish a link between PhoB(PHO80), calcineurin and calcium metabolism. Microarray hybridizations carried out with RNA obtained from wild-type and Delta phoB(PHO80) mutant cells identify Afu4g03610 (phoD(PHO84)), Afu7g06350 (phoE(PHO89)), Afu4g06020 (phoC(PHO81)), and Afu2g09040 (vacuolar transporter Vtc4) as more expressed both in the Delta phoB(PHO80) mutant background and under phosphate-limiting conditions of 0.1 mM P-i. Epifluorescence microscopy revealed accumulation of poly-phosphate in Delta phoB(PHO80) vacuoles, which was independent of extracellular phosphate concentration. Surprisingly, a phoD(PHO84) deletion mutant is indistinguishable phenotypically from the corresponding wild-type strain. mRNA analyses suggest that protein kinase A absence supports the expression of PHO genes in A. fumigatus. Furthermore, Delta phoB(PHO80) and Delta phoD(PHO84) mutant are fully virulent in a murine low dose model for invasive aspergillosis. (C) 2008 Elsevier Inc. All rights reserved.

Functional characterization of the Aspergillus fumigatus CRZ1 homologue, CrzA

The protein phosphatase calcineurin is an important mediator connecting calcium-dependent signalling to various cellular responses in multiple organisms. In fungi calcineurin acts largely through regulating Crz1p-like transcription factors. Here we characterize an Aspergillus fumigatus CRZ1 homologue, CrzA and demonstrate its mediation of cellular tolerance to increased concentrations of calcium and manganese. In addition to acute sensitivitiy to these ions, and decreased conidiation, the crzA null mutant suffers altered expression of calcium transporter mRNAs under high concentrations of calcium, and loss of virulence when compared with the corresponding complemented and wild-type strains. We use multiple expression analyses to probe the transcriptional basis of A. fumigatus calcium tolerance identifying several genes having calA and/or crzA dependent mRNA accumulation patterns. We also demonstrate that contrary to previous findings, the gene encoding the Aspergillus nidulans calcineurin subunit homologue, cnaA, is not essential and that the cnaA deletion mutant shares the morphological phenotypes observed in the corresponding A. fumigatus mutant, Delta calA. Exploiting the A. nidulans model system, we have linked calcineurin activity with asexual developmental induction, finding that CrzA supports appropriate developmental induction in a calcineurin and brlA-dependent manner in both species.

Functional characterization of the putative Aspergillus nidulans DNA damage binding protein homologue DdbA

Nucleotide excision repair (NER) eliminates helix-distorting DNA base lesions. Seven XP-deficient genetic complementation groups (XPA to XPG) have already been identified in mammals, and their corresponding genes have been cloned. Hereditary defects in NER are associated with several diseases, including xeroderma pigmentosum (XP). UV-DDB (XPE) is formed by two associated subunits, DDB1 and DDB2. UV-DDB was identified biochemically as a protein factor that exhibits very strong and specific binding to ultraviolet (UV)-treated DNA. As a preliminary step to characterize the components of the NER in the filamentous fungus Aspergillus nidulans, here we identified a putative DDB1 homologue, DdbA. Deletion and expression analysis indicated that A. nidulans ddbA gene is involved in the DNA damage response, more specifically in the UV light response and 4-nitroquinoline oxide (4-NQO) sensitivity. Furthermore, the Delta ddbA strain cannot self-cross and expression analysis showed that ddbA can be induced by oxidative stress and is developmentally regulated in both asexual and sexual processes. The Delta ddbA mutation can genetically interact with uvsB(ATR), atmA(ATM), nkuA(KU70), H2AX-S129A (a replacement of the conserved serine in the C-terminal of H2AX with alanine), and cshB (a mutation in CSB Cockayne`s syndrome protein involved in the transcription-coupled repair subpathway of NER) mutations. Finally, to determine the DdbA cellular localization, we constructed a GFP:DdbA strain. In the presence and absence of DNA damage, DdbA was mostly detected in the nuclei, indicating that DdbA localizes to nuclei and its cellular localization is not affected by the cellular response to DNA damage induced by 4-NQO and UV light.

A functional matrix metalloproteinase (MMP)-9 polymorphism modifies plasma MMP-9 levels in subjects environmentally exposed to mercury

Mercury (Hg) exposure causes health problems including cardiovascular diseases. Although precise mechanisms have not been precisely defined yet, matrix metalloproteinases (MMPs) may be involved. The gene encoding MMP-9 presents genetic polymorphisms which affect the expression and activity level of this enzyme. Two polymorphisms in the promoter region [C(-1562)T and (CA)(n)] are functionally relevant, and are implicated in several diseases. This study aimed at examining how these polymorphisms affect the circulating MMP-9 levels and its endogenous inhibitor, the tissue inhibitor of metalloproteinase-1 (TIMP-1) in 266 subjects environmentally exposed to Hg. Blood and plasma Hg concentrations were determined by inductively coupled plasma-mass spectrometry (ICP-MS). MMP-9 and TIMP-1 concentrations were measured in plasma samples by gelatin zymography and ELISA, respectively. Genotypes for the C(-1562)T and the microsatellite (CA)(n) polymorphisms were determined. We found a positive association (P<0.05) between plasma Hg concentrations and MMP-9/TIMP-1 ratio (an index of net MMP-9 activity). When the subjects were divided into tertiles with basis on their plasma Hg concentrations, we found that the (CA)(n) polymorphism modified MMP-9 concentrations and MMP-9/TIMP-1 ratio in subjects with the lowest Hg concentrations (first tertile), with the highest MMP-9 levels being found in subjects with genotypes including alleles with 21 or more CA repeats (H alleles) (P<0.05). Conversely, this polymorphism had no effects on subjects with intermediate or high plasma Hg levels (second and third tertiles, respectively). The C(-1562)T polymorphism had no effects on MMP-9 levels. These findings suggest a significant interaction between the (CA)(n) polymorphism and low levels of Hg exposure, possibly increasing the risk of developing diseases in subjects with H alleles. (c) 2010 Elsevier B.V. All rights reserved.

Functional and Evolutionary Insights from the Genomes of Three Parasitoid Nasonia Species

We report here genome sequences and comparative analyses of three closely related parasitoid wasps: Nasonia vitripennis, N. giraulti, and N. longicornis. Parasitoids are important regulators of arthropod populations, including major agricultural pests and disease vectors, and Nasonia is an emerging genetic model, particularly for evolutionary and developmental genetics. Key findings include the identification of a functional DNA methylation tool kit; hymenopteran-specific genes including diverse venoms; lateral gene transfers among Pox viruses, Wolbachia, and Nasonia; and the rapid evolution of genes involved in nuclear-mitochondrial interactions that are implicated in speciation. Newly developed genome resources advance Nasonia for genetic research, accelerate mapping and cloning of quantitative trait loci, and will ultimately provide tools and knowledge for further increasing the utility of parasitoids as pest insect-control agents.

A non-invasive method for silencing gene transcription in honeybees maintained under natural conditions

in the Apis mellifera post-genomic era, RNAi protocols have been used in functional approaches. However, sample manipulation and invasive methods such as injection of double-stranded RNA (dsRNA) can compromise physiology and survival. To circumvent these problems, we developed a non-invasive method for honeybee gene knockdown, using a well-established vitellogenin RNAi system as a model. Second instar larvae received dsRNA for vitellogenin (dsVg-RNA) in their natural diet. For exogenous control, larvae received dsRNA for GFP (dsGFP-RNA). Untreated larvae formed another control group. Around 60% of the treated larvae naturally developed until adult emergence when 0.5 mu g of dsVg-RNA or dsGFP-RNA was offered while no larvae that received 3.0 mu g of dsRNA reached pupal stages. Diet dilution did not affect the removal rates. Viability depends not only on the delivered doses but also on the internal conditions of colonies. The weight of treated and untreated groups showed no statistical differences. This showed that RNAi ingestion did not elicit drastic collateral effects. Approximately 90% of vitellogenin transcripts from 7-day-old workers were silenced compared to controls. A large number of samples are handled in a relatively short time and smaller quantities of RNAi molecules are used compared to invasive methods. These advantages culminate in a versatile and a cost-effective approach. (c) 2008 Elsevier Ltd. All rights reserved.