Using gene-network landscape to dissect genotype effects of TCF7L2 genetic variant on diabetes and cardiovascular risk

Vaquero AR, Ferreira NE, Omae SV, Rodrigues MV, Teixeira SK, Krieger JE, Pereira AC. Using gene-network landscape to dissect genotype effects of TCF7L2 genetic variant on diabetes and cardiovascular risk. Physiol Genomics 44: 903-914, 2012. First published August 7, 2012; doi:10.1152/physiolgenomics.00030.2012.-The single nucleotide polymorphism (SNP) within the TCF7L2 gene, rs7903146, is, to date, the most significant genetic marker associated with Type 2 diabetes mellitus (T2DM) risk. Nonetheless, its functional role in disease pathology is poorly understood. The aim of the present study was to investigate, in vascular smooth muscle cells from 92 patients undergoing aortocoronary bypass surgery, the contribution of this SNP in T2DM using expression levels and expression correlation comparison approaches, which were visually represented as gene interaction networks. Initially, the expression levels of 41 genes (seven TCF7L2 splice forms and 40 other T2DM relevant genes) were compared between rs7903146 wild-type (CC) and T2DM-risk (CT + TT) genotype groups. Next, we compared the expression correlation patterns of these 41 genes between groups to observe if the relationships between genes were different. Five TCF7L2 splice forms and nine genes showed significant expression differences between groups. RXR alpha gene was pinpointed as showing the most different expression correlation pattern with other genes. Therefore, T2DM risk alleles appear to be influencing TCF7L2 splice form's expression in vascular smooth muscle cells, and RXR alpha gene is pointed out as a treatment target candidate for risk reduction in individuals with high risk of developing T2DM, especially individuals harboring TCF7L2 risk genotypes.

Association of AXIN2 with Non-syndromic Oral Clefts in Multiple Populations

We have previously shown the association of AXIN2 with oral clefts in a US population. Here, we expanded our study to explore the association of 11 AXIN2 markers in 682 cleft families from multiple populations. Alleles for each AXIN2 marker were tested for transmission distortion with clefts by means of the Family-based Association Test. We observed an association with SNP rs7224837 and all clefts in the combined populations (p = 0.001), and with SNP rs3923086 and cleft lip and palate in Asian populations (p = 0.004). We confirmed our association findings in an additional 528 cleft families from the United States (p < 0.009). We tested for gene-gene interaction between AXIN2 and additional cleft susceptibility loci. We assessed and detected Axin2 mRNA and protein expression during murine palatogenesis. In addition, we also observed co-localization of Axin2 with Irf6 proteins, particularly in the epithelium. Our results continue to support a role for AXIN2 in the etiology of human clefting. Additional studies should be performed to improve our understanding of the biological mechanisms linking AXIN2 to oral clefts.

Preimplantation development regulatory pathway construction through a text-mining approach

Background: The integration of sequencing and gene interaction data and subsequent generation of pathways and networks contained in databases such as KEGG Pathway is essential for the comprehension of complex biological processes. We noticed the absence of a chart or pathway describing the well-studied preimplantation development stages; furthermore, not all genes involved in the process have entries in KEGG Orthology, important information for knowledge application with relation to other organisms. Results: In this work we sought to develop the regulatory pathway for the preimplantation development stage using text-mining tools such as Medline Ranker and PESCADOR to reveal biointeractions among the genes involved in this process. The genes present in the resulting pathway were also used as seeds for software developed by our group called SeedServer to create clusters of homologous genes. These homologues allowed the determination of the last common ancestor for each gene and revealed that the preimplantation development pathway consists of a conserved ancient core of genes with the addition of modern elements. Conclusions: The generation of regulatory pathways through text-mining tools allows the integration of data generated by several studies for a more complete visualization of complex biological processes. Using the genes in this pathway as “seeds” for the generation of clusters of homologues, the pathway can be visualized for other organisms. The clustering of homologous genes together with determination of the ancestry leads to a better understanding of the evolution of such process.

Changes in spatial and temporal gene expression during incompatible interaction between common bean and anthracnose pathogen

Common bean, one of the most important legumes for human consumption, may have drastic reduction in yield due to anthracnose, a disease caused by the fungus Colletotrichum lindemuthianum. Rapid induction of the plant defense mechanisms is essential to establish an incompatible interaction with this pathogenic fungus. In this study, we evaluated spatial (leaves, epicotyls and hypocotyls) and temporal (24, 48, 72 and 96 hours after inoculation [HAI]) relative expression (RE) of 12 defense-related transcripts selected from previously developed ESTs libraries, during incompatible interaction between the resistant common bean genotype SEL 1308 and the avirulent anthracnose pathogen race 73, using real time quantitative RT-PCR (RT-qPCR) analysis. All selected transcripts, including the ones coding for pathogenesis-related (PR) proteins (PR1a, PR1b, PR2, and PR16a and PR16b) were differentially regulated upon pathogen inoculation. The expression levels of these transcripts were dependent on the tissue and time post inoculation. This study contributes to a better understanding of the kinetics of induced defenses against a fungal pathogen of common bean and may be used as a base line to study defenses against a broad range of pathogens including bacteria as well as non-host resistance. (C) 2012 Elsevier GmbH. All rights reserved.

Assessing the gain of biological data integration in gene networks inference

Background: A current challenge in gene annotation is to define the gene function in the context of the network of relationships instead of using single genes. The inference of gene networks (GNs) has emerged as an approach to better understand the biology of the system and to study how several components of this network interact with each other and keep their functions stable. However, in general there is no sufficient data to accurately recover the GNs from their expression levels leading to the curse of dimensionality, in which the number of variables is higher than samples. One way to mitigate this problem is to integrate biological data instead of using only the expression profiles in the inference process. Nowadays, the use of several biological information in inference methods had a significant increase in order to better recover the connections between genes and reduce the false positives. What makes this strategy so interesting is the possibility of confirming the known connections through the included biological data, and the possibility of discovering new relationships between genes when observed the expression data. Although several works in data integration have increased the performance of the network inference methods, the real contribution of adding each type of biological information in the obtained improvement is not clear. Methods: We propose a methodology to include biological information into an inference algorithm in order to assess its prediction gain by using biological information and expression profile together. We also evaluated and compared the gain of adding four types of biological information: (a) protein-protein interaction, (b) Rosetta stone fusion proteins, (c) KEGG and (d) KEGG+GO. Results and conclusions: This work presents a first comparison of the gain in the use of prior biological information in the inference of GNs by considering the eukaryote (P. falciparum) organism. Our results indicates that information based on direct interaction can produce a higher improvement in the gain than data about a less specific relationship as GO or KEGG. Also, as expected, the results show that the use of biological information is a very important approach for the improvement of the inference. We also compared the gain in the inference of the global network and only the hubs. The results indicates that the use of biological information can improve the identification of the most connected proteins.

Expressão do gene da leptina e seu receptor Ob-Rb no parênquima mamário de novilhas leiteiras

Objetivou-se com este trabalho avaliar os efeitos de uma dieta de alto nível de energia e proteína combinada com a aplicação de bST no perfil de expressão dos genes da leptina e de seu receptor Ob-Rb no parênquima mamário de novilhas leiteiras. Foram utilizadas amostras de parênquima mamário de 32 novilhas holandesas distribuídas aleatoriamente em quatro tratamentos (n=8): dieta com alto ou baixo teor de energia e proteína combinada ou não com a aplicação de bST. O delineamento utilizado foi em blocos casualizados com arranjo de tratamentos em esquema fatorial 2 × 2. A extração do RNA total das amostras de tecido foi feita e o nível de expressão gênica foi analisado por qRT-PCR utilizando-se o gene da glicuronidase β como controle, pelo método 2-ΔΔCt. Animais que receberam a dieta com alto conteúdo de energia e proteína apresentaram maior expressão de mRNA de leptina, com aumento de 56%, e menor expressão de mRNA do receptor Ob-Rb, com redução de 18%. Por outro lado, a aplicação de bST resultou em diminuição da expressão do mRNA de leptina e do receptor Ob-Rb em 74% e 23%, respectivamente. Não houve interação entre dieta e aplicação de bST. O aumento na expressão de leptina pode explicar, ao menos em parte, os efeitos negativos da dieta de alta energia e proteína, oferecida no período pré-púbere, sobre a produção de leite de novilhas leiteiras.

Reciprocal gene flow between conventional and genetically modified soybean cultivars

The objective of this work was to determine the reciprocal gene flow between two soybean cultivars, one tolerant and the other sensitive to glyphosate, as well as to use estimators to determine the outcrossing rate in the population and the number of hybrid seeds in the progeny. The experiment was composed of four blocks of 40 soybean rows, of which 20 rows of each cultivar (CD217 and CD219RR). At the R8 stage, five rows, distant 0.5, 1, 2, 4 and 8 m from the adjacent cultivar, were harvested, threshed and analyzed as for the occurrence of gene flow. As phenotypical markers, the trait color of flowers, hypocotil and pubescence, as well as the tolerance to glyphosate were used. The cultivars contrast for all the analyzed traits, each one conditioned by a single gene with two alleles, in a complete dominance interaction. In the tolerant cultivar progeny, the largest outcross rate was 0.27%, and in the sensitive cultivar progeny, 0.83% was identified; by the dilution effect hypothesis, the outcross rates in natural populations would be 0.104 and 0.388%, respectively. The reciprocal gene flow between CD217 and CD219RR cultivars is not the same in both directions. The proposed estimators are useful for determining the hybrid rates in seed samples.

Association analysis of clinical aspects and vitamin D receptor gene polymorphism with external apical root resorption in orthodontic patients

Introduction: Vitamin D is responsible for the regulation of certain genes at the transcription level, via interaction with the vitamin D receptor, and influences host immune responses and aspects of bone development, growth, and homeostasis. Our aim was to investigate the association of TaqI vitamin D receptor gene polymorphism with external apical root resorption during orthodontic treatment. Methods: Our subjects were 377 patients with Class II Division 1 malocclusion, divided into 3 groups: (1) 160 with external apical root resorption <= 1.43 mm, (2) 179 with external apical root resorption >1.43 mm), and (3) 38 untreated subjects. External apical root resorption of the maxillary incisors was evaluated on periapical radiographs taken before and after 6 months of treatment. After DNA collection and purification, vitamin D receptor TaqI polymorphism analysis was performed by polymerase chain reaction-restriction fragment length polymorphism. Univariate and multivariate analyses were performed to verify the association of clinical and genetic variables with external apical root resorption (P <0.05). Results: There was a higher proportion of external apical root resorption in orthodontically treated patients compared with the untreated subjects. In patients orthodontically treated, age higher than 14 years old, initial size of the maxillary incisor root superior to 30 mm, and premolar extraction were associated with increased external apical root resorption. Genotypes containing the C allele were weakly associated with protection against external apical root resorption (CC + CT x TT [odds ratio, 0.29; 95% confidence interval, 0.07-1.23; P = 0.091]) when treated orthodontic patients were compared to untreated individuals. Conclusions: Clinical factors and vitamin D receptor TaqI polymorphism were associated with external apical root resorption in orthodontic patients. (Am J Orthod Dentofacial Orthop 2012; 142: 339-47)

Interaction among nitric oxide (NO)-related genes in migraine susceptibility

The pathogenic mechanisms involved in migraine are complex and not completely clarified. Because there is evidence for the involvement of nitric oxide (NO) in migraine pathophysiology, candidate gene approaches focusing on genes affecting the endothelial function have been studied including the genes encoding endothelial NO synthase (eNOS), inducible NO synthase (iNOS), and vascular endothelial growth factor (VEGF). However, investigations on gene-gene interactions are warranted to better elucidate the genetic basis of migraine. This study aimed at characterizing interactions among nine clinically relevant polymorphisms in eNOS (T-786C/rs2070744, the 27 bp VNTR in intron 4, the Glu298Asp/rs1799983, and two additional tagSNPs rs3918226 and rs743506), iNOS (C(-1026)A/rs2779249 and G2087A/rs2297518), and VEGF (C(-2578)A/rs699947 and G(-634)C/rs2010963) in migraine patients and control group. Genotypes were determined by real-time polymerase chain reaction using the Taqman(A (R)) allele discrimination assays or PCR and fragment separation by electrophoresis in 99 healthy women without migraine (control group) and in 150 women with migraine divided into two groups: 107 with migraine without aura and 43 with aura. The multifactor dimensionality reduction method was used to detect and characterize gene-gene interactions. We found a significant interaction between eNOS rs743506 and iNOS 2087G/A polymorphisms in migraine patients compared to control group (P < 0.05), suggesting that this combination affect the susceptibility to migraine. Further studies are needed to determine the molecular mechanisms explaining this interaction.

rpb2 is a reliable reference gene for quantitative gene expression analysis in the dermatophyte Trichophyton rubrum

The selection of reference genes used for data normalization to quantify gene expression by real-time PCR amplifications (qRT-PCR) is crucial for the accuracy of this technique. In spite of this, little information regarding such genes for qRT-PCR is available for gene expression analyses in pathogenic fungi. Thus, we investigated the suitability of eight candidate reference genes in isolates of the human dermatophyte Trichophyton rubrum subjected to several environmental challenges, such as drug exposure, interaction with human nail and skin, and heat stress. The stability of these genes was determined by geNorm, NormFinder and Best-Keeper programs. The gene with the most stable expression in the majority of the conditions tested was rpb2 (DNA-dependent RNA polymerase II), which was validated in three T. rubrum strains. Moreover, the combination of rpb2 and chs1 (chitin synthase) genes provided for the most reliable qRT-PCR data normalization in T. rubrum under a broad range of biological conditions. To the best of our knowledge this is the first report on the selection of reference genes for qRT-PCR data normalization in dermatophytes and the results of these studies should permit further analysis of gene expression under several experimental conditions, with improved accuracy and reliability.

Effect of human TGF-beta on the gene expression profile of Schistosoma mansoni adult worms

Schistosoma mansoni is responsible for schistosomiasis, a parasitic disease that affects 200 million people worldwide. Molecular mechanisms of host-parasite interaction are complex and involve a crosstalk between host signals and parasite receptors. TGF-beta signaling pathway has been shown to play an important role in S. mansoni development and embryogenesis. In particular human (h) TGF-beta has been shown to bind to a S. mansoni receptor, transduce a signal that regulates the expression of a schistosome target gene. Here we describe 381 parasite genes whose expression levels are affected by in vitro treatment with hTGF-beta. Among these differentially expressed genes we highlight genes related to morphology, development and cell cycle that could be players of cytokine effects on the parasite. We confirm by qPCR the expression changes detected with microarrays for 5 out of 7 selected genes. We also highlight a set of non-coding RNAs transcribed from the same loci of protein-coding genes that are differentially expressed upon hTCF-beta treatment. These datasets offer potential targets to be explored in order to understand the molecular mechanisms behind the possible role of hTGF-beta effects on parasite biology. (C) 2012 Elsevier B.V. All rights reserved.

Multimodality and flexibility of stochastic gene expression

We consider a general class of mathematical models for stochastic gene expression where the transcription rate is allowed to depend on a promoter state variable that can take an arbitrary (finite) number of values. We provide the solution of the master equations in the stationary limit, based on a factorization of the stochastic transition matrix that separates timescales and relative interaction strengths, and we express its entries in terms of parameters that have a natural physical and/or biological interpretation. The solution illustrates the capacity of multiple states promoters to generate multimodal distributions of gene products, without the need for feedback. Furthermore, using the example of a three states promoter operating at low, high, and intermediate expression levels, we show that using multiple states operons will typically lead to a significant reduction of noise in the system. The underlying mechanism is that a three-states promoter can change its level of expression from low to high by passing through an intermediate state with a much smaller increase of fluctuations than by means of a direct transition.