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.

Diversity of free-living marine nematodes (Enoplida) from Baja California assessed by integrative taxonomy

We used morphological and molecular approaches to evaluate the diversity of free-living marine nematodes (order Enoplida) at four coastal sites in the Gulf of California and three on the Pacific coast of Baja California, Mexico. We identified 22 morphological species belonging to six families, of which Thoracostomopsidae and Oncholaimidae were the most diverse. The genus Mesacanthion (Thoracostomopsidae) was the most widespread and diverse. Five allopatric species, genetically and morphologically differentiated, were found in two localities in the Gulf of California (M. sp1 and M. sp2) and three in the Pacific coast (M. sp3, M. sp4 and M. sp5). Overall, we produced 19 and 20 sequences for the 18S and 28S genes, respectively. Neither gene displayed intraspecific polymorphisms, which allowed us to establish that some morphological variation was likely either ontogenetic or due to phenotypic plasticity. Although 18S and 28S phylogenies were topologically congruent (incongruence length difference test, P > 0.05), divergences between species were much higher in the 28S gene. Moreover, this gene possessed a stronger phylogenetic signal to resolve relationships involving Rhabdodemania and Bathylaimus. On the other hand, the close relationship of Pareurystomina (Enchilidiidae) with oncholaimids warrants further study. The 28S sequences (D2D3 domain) may be better suited for DNA barcoding of marine nematodes than those from the 18S rDNA, particularly for differentiating closely related or cryptic species. Finally, our results underline the relevance of adopting an integrative approach encompassing morphological and molecular analyses to improve the assessment of marine nematode diversity and advance their taxonomy.

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.

Allele-Specific Expression of the MAOA Gene and X Chromosome Inactivation in In Vitro Produced Bovine Embryos

During embryogenesis, one of the two X chromosomes is inactivated in embryos. The production of embryos in vitro may affect epigenetic mechanisms that could alter the expression of genes related to embryo development and X chromosome inactivation (XCI). The aim of this study was to understand XCI during in vitro, pre-implantation bovine embryo development by characterizing the allele-specific expression pattern of the X chromosome-linked gene, monoamine oxidase A (MAOA). Two pools of ten embryos, comprised of the 4-, 8- to 16-cell, morula, blastocyst, and expanded blastocyst stages, were collected. Total RNA from embryos was isolated, and the RT-PCR-RFLP technique was used to observe expression of the MAOA gene. The DNA amplicons were also sequenced using the dideoxy sequencing method. MAOA mRNA was detected, and allele-specific expression was identified in each pool of embryos. We showed the presence of both the maternal and paternal alleles in the 4-, 8-to 16-cell, blastocyst and expanded blastocyst embryos, but only the maternal allele was present in the morula stage. Therefore, we can affirm that the paternal X chromosome is totally inactivated at the morula stage and reactivated at the blastocyst stage. To our knowledge, this is the first report of allele-specific expression of an X-linked gene that is subject to XCI in in vitro bovine embryos from the 4-cell to expanded blastocyst stages. We have established a pattern of XCI in our in vitro embryo production system that can be useful as a marker to assist the development of new protocols for in vitro embryo production. Mol. Reprod. Dev. MoL Reprod. Dev. 77: 615-621, 2010. (C) 2010 Wiley-Liss, Inc.