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 protein-coding and non-coding RNA expression profiles in CD34(+) and in stromal cells in refractory anemia with ringed sideroblasts

Background: Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods: In this study, gene expression profiles of CD34(+) cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results: In CD34(+) cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value <= 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value <= 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions: These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34(+) cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

The Essential Nucleolar Yeast Protein Nop8p Controls the Exosome Function during 60S Ribosomal Subunit Maturation

The yeast nucleolar protein Nop8p has previously been shown to interact with Nip7p and to be required for 60S ribosomal subunit formation. Although depletion of Nop8p in yeast cells leads to premature degradation of rRNAs, the biochemical mechanism responsible for this phenotype is still not known. In this work, we show that the Nop8p amino-terminal region mediates interaction with the 5.8S rRNA, while its carboxyl-terminal portion interacts with Nip7p and can partially complement the growth defect of the conditional mutant strain Dnop8/GAL:NOP8. Interestingly, Nop8p mediates association of Nip7p to pre-ribosomal particles. Nop8p also interacts with the exosome subunit Rrp6p and inhibits the complex activity in vitro, suggesting that the decrease in 60S ribosomal subunit levels detected upon depletion of Nop8p may result from degradation of pre-rRNAs by the exosome. These results strongly indicate that Nop8p may control the exosome function during pre-rRNA processing.

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.

Conservation of dual-targeted proteins in Arabidopsis and rice points to a similar pattern of gene-family evolution

Gene duplication followed by acquisition of specific targeting information and dual targeting were evolutionary strategies enabling organelles to cope with overlapping functions. We examined the evolutionary trend of dual-targeted single-gene products in Arabidopsis and rice genomes. The number of paralogous proteins encoded by gene families and the dual-targeted orthologous proteins were analysed. The number of dual-targeted proteins and the corresponding gene-family sizes were similar in Arabidopsis and rice irrespective of genome sizes. We show that dual targeting of methionine aminopeptidase, monodehydroascorbate reductase, glutamyl-tRNA synthetase, and tyrosyl-tRNA synthetase was maintained despite occurrence of whole-genome duplications in Arabidopsis and rice as well as a polyploidization followed by a diploidization event (gene loss) in the latter.

Atorvastatin effects on SREBF1a and SCAP gene expression in mononuclear cells and its relation with lowering-lipids response

Background: The transcription factors SREBP1 and SCAP are involved in intracellular cholesterol homeostasis. Polymorphisms of these genes have been associated with variations on serum lipid levels and response to statins that are potent cholesterol-lowering drugs. We evaluated the effects of atorvastatin on SREBF1a and SCAP mRNA expression in peripheral blood mononuclear cells (PBMC) and a possible association with gene polymorphisms and lowering-cholesterol response. Methods: Fifty-nine hypercholesterolemic patients were treated with atorvastatin (10 mg/day for 4 weeks). Serum lipid profile and mRNA expression in PBMC were assessed before and after the treatment. Gene expression was quantified by real-time PCR using GAPD as endogenous reference and mRNA expression in HepG2 cells as calibrator. SREBF1 -36delG and SCAP A2386G polymorphisms were detected by PCR-RFLP. Results: Our results showed that transcription of SREBF1a and SCAP was coordinately regulated by atorvastatin (r=0.595, p<0.001), and that reduction in SCAP transcription was associated with the 2386AA genotype (p=0.019). Individuals who responded to atorvastatin with a downregulation of SCAP had also a lower triglyceride compared to those who responded to atorvastatin with an upregulation of SCAP. Conclusion: Atorvastatin has differential effects on SREBF1a and SCAP mRNA expression in PBMC that are associated with baseline transcription levels, triglycerides response to atorvastatin and SCAP A2386G polymorphism. (c) 2008 Elsevier B.V. All rights reserved.

ABCB1 and ABCC1 expression in peripheral mononuclear cells is influenced by gene polymorphisms and atorvastatin treatment

This study investigated the effects of atorvastatin on ABCB1 and ABCC1 mRNA expression on peripheral blood mononuclear cells (PBMC) and their relationship with gene polymorphisms and lowering-cholesterol response. one hundred and thirty-six individuals with hypercholesterolemia were selected and treated with atorvastatin (10 mg/day/4 weeks). Blood samples were collected for serum lipids and apolipoproteins measurements and DNA and RNA extraction. ABCB1 (C3435T and G2677T/A) and ABCC1 (G2012T) gene polymorphisms were identified by polymerase chain reaction-restriction (PCR)-RFLP and mRNA expression was measured in peripheral blood mononuclear cells by singleplex real-time PCR. ABCB1 polymorphisms were associated with risk for coronary artery disease (CAD) (p < 0.05). After atorvastatin treatment, both ABCB1 and ABCC1 genes showed 50% reduction of the mRNA expression (p < 0.05). Reduction of ABCB1 expression was associated with ABCB1 G2677T/A polymorphism (p = 0.039). Basal ABCB1 mRNA in the lower quartile (<0.024) was associated with lower reduction rate of serum low-density lipoprotein (LDL) cholesterol (33.4 +/- 12.4%) and apolipoprotein B (apoB) (17.0 +/- 31.3%) when compared with the higher quartile (>0.085: LDL-c = 40.3 +/- 14.3%; apoB = 32.5 +/- 10.7%; p < 0.05). ABCB1 substrates or inhibitors did not affect the baseline expression, while ABCB1 inhibitors reversed the effects of atorvastatin on both ABCB1 and ABCC1 transporters. In conclusion, ABCB1 and ABCC1 mRNA levels in PBMC are modulated by atorvastatin and ABCB1 G2677T/A polymorphism. and ABCB1 baseline expression is related to differences in serum LDL cholesterol and apoB in response to atorvastatin. (C) 2008 Elsevier Inc. All rights reserved.

Silencing of mitochondrial alternative oxidase gene of Aspergillus fumigatus enhances reactive oxygen species production and killing of the fungus by macrophages

We previously demonstrated that conidia from Aspergillus fumigatus incubated with menadione and paraquat increases activity and expression of cyanide-insensitive alternative oxidase (AOX). Here, we employed the RNA silencing technique in A. fumigatus using the vector pALB1/aoxAf in order to down-regulate the aox gene. Positive transformants for aox gene silencing of A. fumigatus were more susceptible both to an imposed in vitro oxidative stress condition and to macrophages killing, suggesting that AOX is required for the A. fumigatus pathogenicity, mainly for the survival of the fungus conidia during host infection and resistance to reactive oxygen species generated by macrophages.

Leishmania (Viannia) braziliensis transfectants overexpressing the miniexon gene lose virulence in vivo

The miniexon gene has a central role in the processing of polycistronic pre-mRNA of kinetoplastids. It is added to the 5` extremity of each mRNA, supplying the 5`-capped structure to the molecule. Previous studies in Leishmania (Leishmania) major showed that the overexpression of the miniexon array attenuates the Virulence of the parasite in in vivo assays. The results presented here extend those findings to Vionnia subgenus. Leishmania (Vionnia) braziliensis was transfected with a cosmid harboring a tandem array of one hundred miniexon gene copies and then characterized by Northern blot analysis. The overexpression of the exogenous gene was confirmed and its effect on the virulence of L (V.) braziliensis was investigated in hamsters. In BALB/c mice we could not detect parasites during the course of 15 weeks of infection. In addition, hamsters infected with transfectants overexpressing the miniexon gene exhibited only a minor footpad swelling of late onset and failed to develop progressive lesion, these attenuated parasites could be recovered from the inoculation site 1 year after infection. The persistence of parasites in the host indicates that a stable line overexpressing the miniexon may be tested as live vaccine against leishmaniasis. (C) 2008 Elsevier Ireland Ltd. All rights reserved.

The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: Fivefold tandem repetition of a coding region

We sequenced across all of the gene boundaries in the mitochondrial genome of the cattle tick, Boophilus microplus, to determine the arrangement of its genes. The mtDNA of B. microplus has a coding region, composed of tRNA(Glu) and 60 bp of the 3' end of ND1, that is repeated five times. Boophilus microplus is the first coelomate animal known to have more than two copies of a coding sequence. The mitochondrial genome of B, microplus has other unusual features, including (1) reduced T arms in tRNAs, (2) an AT bias in codon use, (3) two control regions that have evolved in concert, (4) three gene rearrangements, and (5) a stem-loop between tRNA(Gln) and tRNA(Phe). The short T arms and small control regions (CRs) of B. microplus and other ticks suggest strong selection for small genomes. Imprecise termination of replication beyond its origin, which can account for the evolution of tandem repeats of coding regions in other mitochondrial genomes, cannot explain the evolution of the fivefold repeated sequence in the mitochondrial genome of B. microplus. Instead, slipped-strand mispairing or recombination are the most plausible explanations for the evolution of these tandem repeats.

Lentiviral vector-mediated tyro sinase-related protein 2 gene transfer to dendritic cells for the therapy of melanoma

Dendritic cells (DCs) are the most potent professional antigen-presenting cells (APCs), which play a vital role in primary immune responses. Introducing genes into DCs will allow constitutive expression of the encoded proteins and thus prolong the presentation of the antigens derived therefrom. In addition, multiple and unidentified epitopes encoded by the entire tumor-associated antigen (TAA) gene may enhance T cell activation. This study demonstrated that an HIV-1-based lentiviral vector conferred efficient gene transfer to DCs. The transgene, murine tyrosinase-related protein 2 (mTRP-2), encodes a clinically relevant melanoma-associated antigen (MAA), which has been found to be a tumor rejection antigen for B16 melanoma. The transfer and proper processing of mTRP-2 in DCs, in terms of RNA transcription activity and protein expression, were verified by RT-PCR and specific antibody, respectively. Administration of mTRP-2 gene-modified DCs (DC-HR'CmT2) to C57BL/6 mice evoked strong protection against tumor challenge, for which the presence of CD4(+) and CD8(+) cells during both the priming and challenge phase was essential. In a therapy model, our results showed that four of seven mice with preestablished tumor remained tumor free for 80 days after therapeutic vaccination. Given the results shown in this study, mTRP-2 gene transfer to DCs provides a potential therapeutic strategy for the management of melanoma, especially in the early stage of the disease.

A preliminary phylogeny of the scale insects (Hemiptera : Sternorrhyncha : Coccoidea) based on nuclear small-subunit ribosomal DNA

Scale insects (Hemiptera: Sternorrhyncha: Coccoidea) are a speciose and morphologically specialized group of plant-feeding bugs in which evolutionary relationships and thus higher classification are controversial. Sequences derived from nuclear small-subunit ribosomal DNA were used to generate a preliminary molecular phylogeny for the Coccoidea based on 39 species representing 14 putative families. Monophyly of the archaeococcoids (comprising Ortheziidae, Margarodidae sensu lato, and Phenacoleachia) was equivocal, whereas monophyly of the neococcoids was supported. Putoidae, represented by Puto yuccae, was found to be outside the remainder of the neococcoid clade. These data are consistent with a single origin (in the ancestor of the neococcoid clade) of a chromosome system involving paternal genome elimination in males. Pseudococcidae (mealybugs) appear to be sister to the rest of the neococcoids and there are indications that Coccidae (soft scales) and Kerriidae (lac scales) are sister taxa. The Eriococcidae (felt scales) was not recovered as a monophyletic group and the eriococcid genus Eriococcus sensu lato was polyphyletic. (C) 2002 Elsevier Science (USA). All rights reserved.

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.

RNAi-mediated silencing of vitellogenin gene function turns honeybee (Apis mellifera) workers into extremely precocious foragers

The switch from within-hive activities to foraging behavior is a major transition in the life cycle of a honeybee (Apis mellifera) worker. A prominent regulatory role in this switch has long been attributed to juvenile hormone (JH), but recent evidence also points to the yolk precursor protein vitellogenin as a major player in behavioral development. In the present study, we injected vitellogenin double-stranded RNA (dsVg) into newly emerged worker bees of Africanized genetic origin and introduced them together with controls into observation hives to record flight behavior. RNA interference-mediated silencing of vitellogenin gene function shifted the onset of long-duration flights (> 10 min) to earlier in life (by 3-4 days) when compared with sham and untreated control bees. In fact, dsVg bees were observed conducting such flights extremely precociously, when only 3 days old. Short-duration flights (< 10 min), which bees usually perform for orientation and cleaning, were not affected. Additionally, we found that the JH titer in dsVg bees collected after 7 days was not significantly different from the controls. The finding that depletion of the vitellogenin titer can drive young bees to become extremely precocious foragers could imply that vitellogenin is the primary switch signal. At this young age, downregulation of vitellogenin gene activity apparently had little effect on the JH titer. As this unexpected finding stands in contrast with previous results on the vitellogenin/JH interaction at a later age, when bees normally become foragers, we propose a three-step sequence in the constellation of physiological parameters underlying behavioral development.

Contractile activity per se induces transcriptional activation of SLC2A4 gene in soleus muscle: involvement of MEF2D, HIF-1a, and TR alpha transcriptional factors

Lima GA, Anhe GF, Giannocco G, Nunes MT, Correa-Giannella ML, Machado UF. Contractile activity per se induces transcriptional activation of SLC2A4 gene in soleus muscle: involvement of MEF2D, HIF-1a, and TR alpha transcriptional factors. Am J Physiol Endocrinol Metab 296: E132-E138, 2009. First published October 28, 2008; doi: 10.1152/ajpendo.90548.2008.-Skeletal muscle is a target tissue for approaches that can improve insulin sensitivity in insulin-resistant states. In muscles, glucose uptake is performed by the GLUT-4 protein, which is encoded by the SLC2A4 gene. SLC2A4 gene expression increases in response to conditions that improve insulin sensitivity, including chronic exercise. However, since chronic exercise improves insulin sensitivity, the increased SLC2A4 gene expression could not be clearly attributed to the muscle contractile activity per se and/or to the improved insulin sensitivity. The present study was designed to investigate the role of contractile activity per se in the regulation of SLC2A4 gene expression as well as in the participation of the transcriptional factors myocyte enhancer factor 2D (MEF2D), hypoxia inducible factor 1a (HIF-1a), and thyroid hormone receptor-alpha (TR alpha). The performed in vitro protocol excluded the interference of metabolic, hormonal, and neural effects. The results showed that, in response to 10 min of electrically induced contraction of soleus muscle, an early 40% increase in GLUT-4 mRNA (30 min) occurred, with a subsequent 65% increase (120 min) in GLUT-4 protein content. EMSA and supershift assays revealed that the stimulus rapidly increased the binding activity of MEF2D, HIF-1a, and TR alpha into the SLC2A4 gene promoter. Furthermore, chromatin immunoprecipitation assay confirmed, in native nucleosome, that contraction induced an approximate fourfold (P < 0.01) increase in MEF2D and HIF-1a-binding activity. In conclusion, muscle contraction per se enhances SLC2A4 gene expression and that involves MEF2D, HIF-1a, and TR alpha transcription factor activation. This finding reinforces the importance of physical activity to improve glycemic homeostasis independently of other additional insulin sensitizer approaches.