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.

Identification of COL6A1 as a differentially expressed gene in human astrocytomas

Diffuse infiltrating gliomas are the most common tumors of the central nervous system. Gliomas are classified by the WHO according to their histopathological and clinical characteristics into four classes: grade I (pilocytic astrocytoma), grade II (diffuse astrocytoma), grade III (anaplastic astrocytoma), and grade IV (glioblastoma multiforme). Several genes have already been correlated with astrocytomas, but many others are yet to be uncovered. By analyzing the public SAGE data from 21 patients, comprising low malignant grade astrocytomas and glioblastomas, we found COL6A1 to be differentially expressed, confirming this finding by real time RT-PCR in 66 surgical samples. To the best of our knowledge, COL6A1 has never been described in gliomas. The expression of this gene has significantly different means when normal glia is compared with low-grade astrocytomas (grades I and II) and high-grade astrocytomas (grades III and IV), with a tendency to be greater in higher grade samples, thus rendering it a powerful tumor marker.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Construct and Compare Gene Coexpression Networks with DAPfinder and DAPview

Background: DAPfinder and DAPview are novel BRB-ArrayTools plug-ins to construct gene coexpression networks and identify significant differences in pairwise gene-gene coexpression between two phenotypes. Results: Each significant difference in gene-gene association represents a Differentially Associated Pair (DAP). Our tools include several choices of filtering methods, gene-gene association metrics, statistical testing methods and multiple comparison adjustments. Network results are easily displayed in Cytoscape. Analyses of glioma experiments and microarray simulations demonstrate the utility of these tools. Conclusions: DAPfinder is a new friendly-user tool for reconstruction and comparison of biological networks.

Beta-2 adrenergic receptor gene polymorphisms Gln27Glu, Arg16Gly in patients with heart failure

Background -: Beta-2 adrenergic receptor gene polymorphisms Gln27Glu, Arg16Gly and Thr164Ile were suggested to have an effect in heart failure. We evaluated these polymorphisms relative to clinical characteristics and prognosis of alarge cohort of patients with heart failure of different etiologies. Methods -: We studied 501 patients with heart failure of different etiologies. Mean age was 58 years (standard deviation 14.4 years), 298 (60%) were men. Polymorphisms were identified by polymerase chain reaction-restriction fragment length polymorphism. Results -: During the mean follow-up of 12.6 months (standard deviation 10.3 months), 188 (38%) patients died. Distribution of genotypes of polymorphism Arg16Gly was different relative to body mass index (chi(2) = 9.797; p = 0.04). Overall the probability of survival was not significantly predicted by genotypes of Gln27Glu, Arg16Gly, or Thr164Ile. Allele and haplotype analysis also did not disclose any significant difference regarding mortality. Exploratory analysis through classification trees pointed towards a potential association between the Gln27Glu polymorphism and mortality in older individuals. Conclusion -: In this study sample, we were not able to demonstrate an overall influence of polymorphisms Gln27Glu and Arg16Gly of beta-2 receptor gene on prognosis. Nevertheless, Gln27Glu polymorphism may have a potential predictive value in older individuals.

Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice

Lacchini S, Heimann AS, Evangelista FS, Cardoso L, Silva GJ, Krieger JE. Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice. Physiol Genomics 37: 225-230, 2009. First published March 3, 2009; doi:10.1152/physiolgenomics.90288.2008.-We tested the hypothesis that small changes in angiotensin I-converting enzyme (ACE) expression can alter the vascular response to injury. Male mice containing one, two, three, and four copies of the Ace gene with no detectable vascular abnormality or changes in blood pressure were submitted to cuff-induced femoral artery injury. Femoral thickening was higher in 3- and 4-copy mice (42.4 +/- 4.3% and 45.7 +/- 6.5%, respectively) compared with 1- and 2-copy mice (8.3 +/- 1.3% and 8.5 +/- 0.9%, respectively). Femoral ACE levels from control and injured vessels were assessed in 1- and 3-copy Ace mice, which represent the extremes of the observed response. ACE vascular activity was higher in 3- vs. 1-copy Ace mice (2.4-fold, P < 0.05) in the control uninjured vessel. Upon injury, ACE activity significantly increased in both groups [2.41-fold and 2.14-fold (P < 0.05) for 1- and 3-copy groups, respectively] but reached higher levels in 3- vs. 1-copy Ace mice (P < 0.05). Pharmacological interventions were then used as a counterproof and to indirectly assess the role of angiotensin II (ANG II) on this response. Interestingly, ACE inhibition (enalapril) and ANG II AT(1) receptor blocker (losartan) reduced intima thickening in 3-copy mice to 1-copy mouse values (P < 0.05) while ANG II treatment significantly increased intima thickening in 1-copy mice to 3-copy mouse levels (P < 0.05). Together, these data indicate that small physiologically relevant changes in ACE, not associated with basal vascular abnormalities or blood pressure levels, do influence the magnitude of cuff-induced neointima thickening in mice.

Exercise training improves muscle vasodilatation in individuals with T786C polymorphism of endothelial nitric oxide synthase gene

Negrão M.V, Alves CR, Alves G.B, Pereira A.C, Dias R.G, Laterza M.C, Mota G.F, Oliveira E.M, Bassaneze V, Krieger J.E, Negrão C.E, Rondon M.U.P. Exercise training improves muscle vasodilatation in individuals with T786C polymorphism of endothelial nitric oxide synthase gene. Physiol Genomics 42A: 71-77, 2010. First published July 6, 2010; doi:10.1152/physiolgenomics.00145.2009.-Allele T at promoter region of the eNOS gene has been associated with an increase in coronary disease mortality, suggesting that this allele increases susceptibility for endothelial dysfunction. In contrast, exercise training improves endothelial function. Thus, we hypothesized that: 1) Muscle vasodilatation during exercise is attenuated in individuals homozygous for allele T, and 2) Exercise training improves muscle vasodilatation in response to exercise for TT genotype individuals. From 133 preselected healthy individuals genotyped for the T786C polymorphism, 72 participated in the study: TT (n = 37; age 27 +/- 1 yr) and CT + CC (n = 35; age 26 +/- 1 yr). Forearm blood flow (venous occlusion plethysmography) and blood pressure (oscillometric automatic cuff) were evaluated at rest and during 30% handgrip exercise. Exercise training consisted of three sessions per week for 18 wk, with intensity between anaerobic threshold and respiratory compensation point. Resting forearm vascular conductance (FVC, P = 0.17) and mean blood pressure (P = 0.70) were similar between groups. However, FVC responses during handgrip exercise were significantly lower in TT individuals compared with CT + CC individuals (0.39 +/- 0.12 vs. 1.08 +/- 0.27 units, P = 0.01). Exercise training significantly increased peak VO(2) in both groups, but resting FVC remained unchanged. This intervention significantly increased FVC response to handgrip exercise in TT individuals (P = 0.03), but not in CT + CC individuals (P = 0.49), leading to an equivalent FVC response between TT and CT + CC individuals (1.05 +/- 0.18 vs. 1.59 +/- 0.27 units, P = 0.27). In conclusion, exercise training improves muscle vasodilatation in response to exercise in TT genotype individuals, demonstrating that genetic variants influence the effects of interventions such as exercise training.

Molecular characterization of a miraculin-like gene differentially expressed during coffee development and coffee leaf miner infestation

The characterization of a coffee gene encoding a protein similar to miraculin-like proteins, which are members of the plant Kunitz serine trypsin inhibitor (STI) family of proteinase inhibitors (PIs), is described. PIs are important proteins in plant defence against insects and in the regulation of proteolysis during plant development. This gene has high identity with the Richadella dulcifica taste-modifying protein miraculin and with the tomato protein LeMir; and was named as CoMir (Coffea miraculin). Structural protein modelling indicated that CoMir had structural similarities with the Kunitz STI proteins, but suggested specific folding structures. CoMir was up-regulated after coffee leaf miner (Leucoptera coffella) oviposition in resistant plants of a progeny derived from crosses between C. racemosa (resistant) and C. arabica (susceptible). Interestingly, this gene was down-regulated during coffee leaf miner herbivory in susceptible plants. CoMir expression was up-regulated after abscisic acid application and wounding stress and was prominent during the early stages of flower and fruit development. In situ hybridization revealed that CoMir transcripts accumulated in the anther tissues that display programmed cell death (tapetum, endothecium and stomium) and in the metaxylem vessels of the petals, stigma and leaves. In addition, the recombinant protein CoMir shows inhibitory activity against trypsin. According to the present results CoMir may act in proteolytic regulation during coffee development and in the defence against L. coffeella. The similarity of CoMir with other Kunitz STI proteins and the role of CoMir in plant development and plant stress are discussed.

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.

Functional characterization of the oxaloacetase encoding gene and elimination of oxalate formation in the beta-lactam producer Penicillium chrysogenum

Penicillium chrysogenum is widely used as an industrial antibiotic producer, in particular in the synthesis of g-lactam antibiotics such as penicillins and cephalosporins. In industrial processes, oxalic acid formation leads to reduced product yields. Moreover, precipitation of calcium oxalate complicates product recovery. We observed oxalate production in glucose-limited chemostat cultures of P. chrysogenum grown with or without addition of adipic acid, side-chain of the cephalosporin precursor adipoyl-6-aminopenicillinic acid (ad-6-APA). Oxalate accounted for up to 5% of the consumed carbon source. In filamentous fungi, oxaloacetate hydrolase (OAH; EC3.7.1.1) is generally responsible for oxalate production. The P. chrysogenum genome harbours four orthologs of the A. niger oahA gene. Chemostat-based transcriptome analyses revealed a significant correlation between extracellular oxalate titers and expression level of the genes Pc18g05100 and Pc22g24830. To assess their possible involvement in oxalate production, both genes were cloned in Saccharomyces cerevisiae, yeast that does not produce oxalate. Only the expression of Pc22g24830 led to production of oxalic acid in S. cerevisiae. Subsequent deletion of Pc22g28430 in P. chrysogenum led to complete elimination of oxalate production, whilst improving yields of the cephalosporin precursor ad-6-APA. (C) 2011 Elsevier Inc. All rights reserved.

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.

Transgenic Sweet Orange (Citrus sinensis L. Osbeck) Expressing the attacin A Gene for Resistance to Xanthomonas citri subsp citri

Genetic transformation with genes that code for antimicrobial peptides has been an important strategy used to control bacterial diseases in fruit crops, including apples, pears, and citrus. Asian citrus canker (ACC) caused by Xanthomonas citri subsp. citri Schaad et al. (Xcc) is a very destructive disease, which affects the citrus industry in most citrus-producing areas of the world. Here, we report the production of genetically transformed Natal, Pera, and Valencia sweet orange cultivars (Citrus sinensis L. Osbeck) with the insect-derived attacin A (attA) gene and the evaluation of the transgenic plants for resistance to Xcc. Agrobacterium tumefaciens Smith and Towns-mediated genetic transformation experiments involving these cultivars led to the regeneration of 23 different lines. Genetically transformed plants were identified by polymerase chain reaction, and transgene integration was confirmed by Southern blot analyses. Transcription of attA gene was detected by Northern blot analysis in all plants, except for one Natal sweet orange transformation event. Transgenic lines were multiplied by grafting onto Rangpur lime rootstock plants (Citrus limonia Osbeck) and spray-inoculated with an Xcc suspension (10(6) cfu mL(-1)). Experiments were repeated three times in a completely randomized design with seven to ten replicates. Disease severity was determined in all transgenic lines and in the control (non-transgenic) plants 30 days after inoculation. Four transgenic lines of Valencia sweet orange showed a significant reduction in disease severity caused by Xcc. These reductions ranged from 58.3% to 77.8%, corresponding to only 0.16-0.30% of leaf diseased area as opposed to 0.72% on control plants. One transgenic line of Natal sweet orange was significantly more resistant to Xcc, with a reduction of 45.2% comparing to the control plants, with only 0.14% of leaf diseased area. Genetically transformed Pera sweet orange plants expressing attA gene did not show a significant enhanced resistance to Xcc, probably due to its genetic background, which is naturally more resistant to this pathogen. The potential effect of attacin A antimicrobial peptide to control ACC may be related to the genetic background of each sweet orange cultivar regarding their natural resistance to the pathogen.

Genetic transformation of Citrus sinensis cv. Hamlin with hrpN gene from Erwinia amylovora and evaluation of the transgenic lines for resistance to citrus canker

Transgenic Citrus sinensis (L.) Osb. cv. Hamlin plants expressing the hrpN gene were obtained by Agrobacterium tumefaciens (Smith and Towns) Conn-mediated transformation. hrpN encodes a harpin protein, which elicits the hypersensitive response and systemic acquired resistance in plants. The gene construct consisted of gst1, a pathogen-inducible promoter, a signal peptide for protein secretion to the apoplast, the selection genes nptI1 or aacC1 and the Nos terminator. The function of gst1 in citrus was evaluated in transgenic C. sinensis cv. Valencia harboring the reporter gene uidA (gus) driven by this promoter. Histochemical analysis for gus revealed that gst1 is activated in citrus leaves by both wounding and inoculation with Xanthomonas axonopodis Starr and Garces pv. citri (Hasse) Vauterin et al. Genetic transformation was confirmed by Southern blot hybridization in eight cv. Hamlin acclimatized plants. RT-PCR confirmed hrpN gene expression in seven cv. Hamlin transgenic lines before pathogen inoculation. Some hrpN transgenic lines showed severe leaf curling and abnormal growth. Six hrpN transgenic lines were propagated and evaluated for susceptibility to X axonopodis pv. citri. RT-PCR confirmed gene expression in all six hrpN transgenic lines after pathogen inoculation. Several of the hrpN transgenic lines showed reduction in susceptibility to citrus canker as compared with non-transgenic plants. One hrpN transgenic line exhibited normal vegetative development and displayed very high resistance to the pathogen, estimated as up to 79% reduction in disease severity. This is the first report of genetic transformation of citrus using a pathogen-inducible promoter and the hrpN gene. Further evaluations of the transgenic plants under field conditions are planned. Nevertheless, the evidence to date suggests that the hrpN gene reduces the susceptibility of citrus plants to the canker disease. (C) 2009 Elsevier B.V. All rights reserved.

Immunological responses and cytokine gene expression analysis to Cooperia punctata infections in resistant and susceptible Nelore cattle

Cellular and humoral immune response, as well as cytokine gene expression, was assessed in Nelore cattle with different degrees of resistance to Cooperia punctata natural infection. One hundred cattle (male, weaned, 11-12 months old), kept together on pasture, were evaluated. Faecal and blood samples were collected for parasitological and immunological assays. Based on nematode faecal egg counts (FEC) and worm burden, the seven most resistant and the eight most susceptible animals were selected. Tissue samples of the small intestine were collected for histological quantification of inflammatory cells and analysis of cytokine gene expression (IL-2, IL-4, IL-8, IL-1 2p35, IL-13, TNF-alpha, IFN-gamma, MCP-1, MCP-2, and MUC- 1) using real-time RT-PCR. Mucus samples were also collected for IgA levels determination. Serum IgG1 mean levels against C. punctata antigens were higher in the resistant group, but significant differences between groups were only observed 14 days after the beginning of the experiment against infective larvae (1-3) and 14 and 84 days against adult antigens. The resistant group also presented higher IgA levels against C. punctata (L3 and adult) antigens with significant difference 14 days after the beginning of the trial (P < 0.05). In the small-intestine mucosa, levels of IgA anti-L3 and anti-adult C. punctata were higher in the resistant group, compared with the susceptible group (P < 0.05). Gene expression of both T(H)2 cytokines (IL-4 and IL-13) in the resistant group and T(H)1 cytokines (IL-2, IL-1 2p35, IFN-gamma and MCP-1) in the susceptible group was up-regulated. Such results suggested that immune response to C. punctata was probably mediated by TH2 cytokines in the resistant group and by T(H)1 cytokines in the susceptible group. (C) 2008 Elsevier B.V. All rights reserved.

A new polymorphism in the Growth and Differentiation Factor 9 (GDF9) gene is associated with increased ovulation rate and prolificacy in homozygous sheep

P>Brazilian Santa Ines (SI) sheep are very well-adapted to the tropical conditions of Brazil and are an important source of animal protein. A high rate of twin births was reported in some SI flocks. Growth and Differentiation Factor 9 (GDF9) and Bone Morphogenetic Protein 15 (BMP15) are the first two genes expressed by the oocyte to be associated with an increased ovulation rate in sheep. All GDF9 and BMP15 variants characterized, until now, present the same phenotype: the heterozygote ewes have an increased ovulation rate and the mutated homozygotes are sterile. In this study, we have found a new allele of GDF9, named FecGE (Embrapa), which leads to a substitution of a phenylalanine with a cysteine in a conservative position of the mature peptide. Homozygote ewes presenting the FecGE allele have shown an increase in their ovulation rate (82%) and prolificacy (58%). This new phenotype can be very useful in better understanding the genetic control of follicular development; the mechanisms involved in the control of ovulation rate in mammals; and for the improvement of sheep production.