Modeling stochasticity and robustness in gene regulatory networks.

MOTIVATION: Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations. RESULTS: In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs. AVAILABILITY: Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/ approximately garg/genysis.html.

Changes in gene expression profile in human subcutaneous adipose tissue during significant weight loss.

Objective: To analyze the expression of peroxisome proliferator-activated receptor-γ1 and 2 (PPARγ1 and 2), 11β-hydroxysteroid dehydrogenase type 1 (11βHSD1), and leptin in adipose tissue (AT) of obese women during weight loss following Roux-en-Y gastric bypass (RYGB) and to compare these levels with those obtained in AT of nonobese subjects. Methods: Gene expression was determined by real-time RT-PCR prior to surgery and at 3, 6, and 12 months after RYGB. Results: All obese patients lost weight, reaching a mean BMI of 29.3 ± 1.0 kg/m(2) at 1 year after surgery (-33.9 ± 1.5% of their initial body weight). In obese subjects leptin and 11βHSD1 were over-expressed, whereas PPARγ1 was expressed at lower levels compared to controls. After surgery, leptin and 11βHSD1 gene expression decreased, whereas PPARγ1 expression increased. At 12 months after RYGB, these 3 genes had reached levels similar to the controls. In contrast, PPARγ2 gene expression was not different between groups and types of tissue and remained unchanged during weight loss. We found a positive correlation between BMI and levels of gene expression of leptin and 11βHSD1. Conclusion: Gene expression of leptin, PPARγ1, and 11βHSD1 in AT is modified in human obesity. This default is completely corrected by RYGB. Copyright © 2012 S. Karger GmbH, Freiburg.

Regulatory elements involved in the post-transcriptional control of stage-specific gene expression in Trypanosoma cruzi: a review

Trypanosoma cruzi, a protozoan parasite that causes Chagas disease, exhibits unique mechanisms for gene expression such as constitutive polycistronic transcription of protein-coding genes, RNA editing and trans-splicing. In the absence of mechanism controlling transcription initiation, organized subsets of T. cruzi genes must be post-transcriptionally co-regulated in response to extracellular signals. The mechanisms that regulate stage-specific gene expression in this parasite have become much clearer through sequencing its whole genome as well as performing various proteomic and microarray analyses, which have demonstrated that at least half of the T. cruzi genes are differentially regulated during its life cycle. In this review, we attempt to highlight the recent advances in characterising cis and trans-acting elements in the T. cruzi genome that are involved in its post-transcriptional regulatory machinery.

VEGF gene expression in adult human thymus fat: a correlative study with hypoxic induced factor and cyclooxygenase-2.

It is well known that the adult human thymus degenerates into fat tissue; however, it has never been considered as a potential source of angiogenic factors. Recently, we have described that this fat (TAT) produces angiogenic factors and induces human endothelial cell proliferation and migration, indicating its potential angiogenic properties. DESIGN Adult thymus fat and subcutaneous adipose tissue specimens were obtained from 28 patients undergoing cardiac surgery, making this tissue readily available as a prime source of adipose tissue. We focused our investigation on determining VEGF gene expression and characterizing the different genes, mediators of inflammation and adipogenesis, and which are known to play a relevant role in angiogenesis regulation. RESULTS We found that VEGF-A was the isoform most expressed in TAT. This expression was accompanied by an upregulation of HIF-1alpha, COX-2 and HO-1 proteins, and by increased HIF-1 DNA binding activity, compared to SAT. Furthermore, we observed that TAT contains a high percentage of mature adipocytes, 0.25% of macrophage cells, 15% of endothelial cells and a very low percentage of thymocyte cells, suggesting the cellular variability of TAT, which could explain the differences in gene expression observed in TAT. Subsequently, we showed that the expression of genes known as adipogenic mediators, including PPARgamma1/gamma2, FABP-4 and adiponectin was similar in both TAT and SAT. Moreover the expression of these latter genes presented a significantly positive correlation with VEGF, suggesting the potential association between VEGF and the generation of adipose tissue in adult thymus. CONCLUSION Here we suggest that this fat has a potential angiogenic function related to ongoing adipogenesis, which substitutes immune functions within the adult thymus. The expression of VEGF seems to be associated with COX-2, HO-1 and adipogenesis related genes, suggesting the importance that this new fat has acquired in research in relation to adipogenesis and angiogenesis.

Adipose tissue gene expression of factors related to lipid processing in obesity.

BACKGROUND Adipose tissue lipid storage and processing capacity can be a key factor for obesity-related metabolic disorders such as insulin resistance and diabetes. Lipid uptake is the first step to adipose tissue lipid storage. The aim of this study was to analyze the gene expression of factors involved in lipid uptake and processing in subcutaneous (SAT) and visceral (VAT) adipose tissue according to body mass index (BMI) and the degree of insulin resistance (IR). METHODS AND PRINCIPAL FINDINGS VLDL receptor (VLDLR), lipoprotein lipase (LPL), acylation stimulating protein (ASP), LDL receptor-related protein 1 (LRP1) and fatty acid binding protein 4 (FABP4) gene expression was measured in VAT and SAT from 28 morbidly obese patients with Type 2 Diabetes Mellitus (T2DM) or high IR, 10 morbidly obese patients with low IR, 10 obese patients with low IR and 12 lean healthy controls. LPL, FABP4, LRP1 and ASP expression in VAT was higher in lean controls. In SAT, LPL and FABP4 expression were also higher in lean controls. BMI, plasma insulin levels and HOMA-IR correlated negatively with LPL expression in both VAT and SAT as well as with FABP4 expression in VAT. FABP4 gene expression in SAT correlated inversely with BMI and HOMA-IR. However, multiple regression analysis showed that BMI was the main variable contributing to LPL and FABP4 gene expression in both VAT and SAT. CONCLUSIONS Morbidly obese patients have a lower gene expression of factors related with lipid uptake and processing in comparison with healthy lean persons.

Methods for the ex vivo characterization of human CD8+ T subsets based on gene expression and replicative history analysis.

The generation of an antigen-specific T-lymphocyte response is a complex multi-step process. Upon T-cell receptor-mediated recognition of antigen presented by activated dendritic cells, naive T-lymphocytes enter a program of proliferation and differentiation, during the course of which they acquire effector functions and may ultimately become memory T-cells. A major goal of modern immunology is to precisely identify and characterize effector and memory T-cell subpopulations that may be most efficient in disease protection. Sensitive methods are required to address these questions in exceedingly low numbers of antigen-specific lymphocytes recovered from clinical samples, and not manipulated in vitro. We have developed new techniques to dissect immune responses against viral or tumor antigens. These allow the isolation of various subsets of antigen-specific T-cells (with major histocompatibility complex [MHC]-peptide multimers and five-color FACS sorting) and the monitoring of gene expression in individual cells (by five-cell reverse transcription-polymerase chain reaction [RT-PCR]). We can also follow their proliferative life history by flow-fluorescence in situ hybridization (FISH) analysis of average telomere length. Recently, using these tools, we have identified subpopulations of CD8+ T-lymphocytes with distinct proliferative history and partial effector-like properties. Our data suggest that these subsets descend from recently activated T-cells and are committed to become differentiated effector T-lymphocytes.

Modulation of MDR1 and MRP3 Gene Expression in Lung Cancer Cells after Paclitaxel and Carboplatin Exposure

Carboplatin-paclitaxel is a reference regimen in the treatment of locally advanced or disseminated non-small cell lung cancer (NSCLC). This paper discusses the multidrug resistance developed with this drug combination, which is one of the major obstacles to successful treatment. In order to understand and overcome the drug resistance pattern of NSCLC after carboplatin plus paclitaxel exposure, levels of mRNA expression of multidrug resistance 1 (MDR1) and multidrug resistance-associated protein 3 (MRP3) were investigated in primary NSCLC cell lines (A-549 and A-427) and a metastasis-derived NSCLC cell line (NODO). Our results showed that exposure of the three NSCLC lines to plasma concentrations of paclitaxel (5 μM) produced an increase in MDR1 expression, while MRP3 showed no alteration in expression. By contrast, the same cells exposed to carboplatin plasma concentrations (30 μM) showed overexpression of MRP3. In these cells, MDR1 showed no expression changes. Interestingly, the combination of both paclitaxel and carboplatin caused increased expression of the MDR1 drug resistance gene rather than the individual treatments. These results suggest that carboplatin and paclitaxel may induce drug resistance mediated by MDR1 and MRP3, which may be enhanced by the simultaneous use of both drugs.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Zinc-alpha 2-glycoprotein gene expression in adipose tissue is related with insulin resistance and lipolytic genes in morbidly obese patients.

OBJECTIVE Zinc-α(2) glycoprotein (ZAG) stimulates lipid loss by adipocytes and may be involved in the regulation of adipose tissue metabolism. However, to date no studies have been made in the most extreme of obesity. The aims of this study are to analyze ZAG expression levels in adipose tissue from morbidly obese patients, and their relationship with lipogenic and lipolytic genes and with insulin resistance (IR). METHODS mRNA expression levels of PPARγ, IRS-1, IRS-2, lipogenic and lipolytic genes and ZAG were quantified in visceral (VAT) and subcutaneous adipose tissue (SAT) of 25 nondiabetic morbidly obese patients, 11 with low IR and 14 with high IR. Plasma ZAG was also analyzed. RESULTS The morbidly obese patients with low IR had a higher VAT ZAG expression as compared with the patients with high IR (p = 0.023). In the patients with low IR, the VAT ZAG expression was greater than that in SAT (p = 0.009). ZAG expression correlated between SAT and VAT (r = 0.709, p<0.001). VAT ZAG expression was mainly predicted by insulin, HOMA-IR, plasma adiponectin and expression of adiponectin and ACSS2. SAT ZAG expression was only predicted by expression of ATGL. CONCLUSIONS ZAG could be involved in modulating lipid metabolism in adipose tissue and is associated with insulin resistance. These findings suggest that ZAG may be a useful target in obesity and related disorders, such as diabetes.

Real-time recording of circadian liver gene expression in freely moving mice reveals the phase-setting behavior of hepatocyte clocks.

The mammalian circadian timing system consists of a master pacemaker in the suprachiasmatic nucleus (SCN) in the hypothalamus, which is thought to set the phase of slave oscillators in virtually all body cells. However, due to the lack of appropriate in vivo recording technologies, it has been difficult to study how the SCN synchronizes oscillators in peripheral tissues. Here we describe the real-time recording of bioluminescence emitted by hepatocytes expressing circadian luciferase reporter genes in freely moving mice. The technology employs a device dubbed RT-Biolumicorder, which consists of a cylindrical cage with reflecting conical walls that channel photons toward a photomultiplier tube. The monitoring of circadian liver gene expression revealed that hepatocyte oscillators of SCN-lesioned mice synchronized more rapidly to feeding cycles than hepatocyte clocks of intact mice. Hence, the SCN uses signaling pathways that counteract those of feeding rhythms when their phase is in conflict with its own phase.

The MRC1/CD68 Ratio Is Positively Associated with Adipose Tissue Lipogenesis and with Muscle Mitochondrial Gene Expression in Humans.

BACKGROUND Alternative macrophages (M2) express the cluster differentiation (CD) 206 (MCR1) at high levels. Decreased M2 in adipose tissue is known to be associated with obesity and inflammation-related metabolic disturbances. Here we aimed to investigate MCR1 relative to CD68 (total macrophages) gene expression in association with adipogenic and mitochondrial genes, which were measured in human visceral [VWAT, n = 147] and subcutaneous adipose tissue [SWAT, n = 76] and in rectus abdominis muscle (n = 23). The effects of surgery-induced weight loss were also longitudinally evaluated (n = 6). RESULTS MCR1 and CD68 gene expression levels were similar in VWAT and SWAT. A higher proportion of CD206 relative to total CD68 was present in subjects with less body fat and lower fasting glucose concentrations. The ratio MCR1/CD68was positively associated with IRS1gene expression and with the expression of lipogenic genes such as ACACA, FASN and THRSP, even after adjusting for BMI. The ratio MCR1/CD68 in SWAT increased significantly after the surgery-induced weight loss (+44.7%; p = 0.005) in parallel to the expression of adipogenic genes. In addition, SWAT MCR1/CD68ratio was significantly associated with muscle mitochondrial gene expression (PPARGC1A, TFAM and MT-CO3). AT CD206 was confirmed by immunohistochemistry to be specific of macrophages, especially abundant in crown-like structures. CONCLUSION A decreased ratio MCR1/CD68 is linked to adipose tissue and muscle mitochondrial dysfunction at least at the level of expression of adipogenic and mitochondrial genes.

Genome-wide mRNA expression correlates of viral control in CD4+ T-cells from HIV-1-infected individuals.

There is great interindividual variability in HIV-1 viral setpoint after seroconversion, some of which is known to be due to genetic differences among infected individuals. Here, our focus is on determining, genome-wide, the contribution of variable gene expression to viral control, and to relate it to genomic DNA polymorphism. RNA was extracted from purified CD4+ T-cells from 137 HIV-1 seroconverters, 16 elite controllers, and 3 healthy blood donors. Expression levels of more than 48,000 mRNA transcripts were assessed by the Human-6 v3 Expression BeadChips (Illumina). Genome-wide SNP data was generated from genomic DNA using the HumanHap550 Genotyping BeadChip (Illumina). We observed two distinct profiles with 260 genes differentially expressed depending on HIV-1 viral load. There was significant upregulation of expression of interferon stimulated genes with increasing viral load, including genes of the intrinsic antiretroviral defense. Upon successful antiretroviral treatment, the transcriptome profile of previously viremic individuals reverted to a pattern comparable to that of elite controllers and of uninfected individuals. Genome-wide evaluation of cis-acting SNPs identified genetic variants modulating expression of 190 genes. Those were compared to the genes whose expression was found associated with viral load: expression of one interferon stimulated gene, OAS1, was found to be regulated by a SNP (rs3177979, p = 4.9E-12); however, we could not detect an independent association of the SNP with viral setpoint. Thus, this study represents an attempt to integrate genome-wide SNP signals with genome-wide expression profiles in the search for biological correlates of HIV-1 control. It underscores the paradox of the association between increasing levels of viral load and greater expression of antiviral defense pathways. It also shows that elite controllers do not have a fully distinctive mRNA expression pattern in CD4+ T cells. Overall, changes in global RNA expression reflect responses to viral replication rather than a mechanism that might explain viral control.

Phospholipase gene expression during Paracoccidioides brasiliensis morphological transition and infection

Phospholipase is an important virulence factor for pathogenic fungi. In this study, we demonstrate the following: (i) the Paracoccidioides brasiliensis pld gene is preferentially expressed in mycelium cells, (ii) the plb1 gene is mostly up-regulated by infection after 6 h of co-infection of MH-S cells or during BALB/c mice lung infection, (iii) during lung infection, plb1, plc and pld gene expression are significantly increased 6-48 h post-infection compared to 56 days after infection, strongly suggesting that phospholipases play a role in the early events of infection, but not during the chronic stages of pulmonary infection by P. brasiliensis.

Effects of ploidy and sex-locus genotype on gene expression patterns in the fire ant Solenopsis invicta.

Males in many animal species differ greatly from females in morphology, physiology and behaviour. Ants, bees and wasps have a haplodiploid mechanism of sex determination whereby unfertilized eggs become males while fertilized eggs become females. However, many species also have a low frequency of diploid males, which are thought to develop from diploid eggs when individuals are homozygous at one or more sex determination loci. Diploid males are morphologically similar to haploids, though often larger and typically sterile. To determine how ploidy level and sex-locus genotype affect gene expression during development, we compared expression patterns between diploid males, haploid males and females (queens) at three developmental timepoints in Solenopsis invicta. In pupae, gene expression profiles of diploid males were very different from those of haploid males but nearly identical to those of queens. An unexpected shift in expression patterns emerged soon after adult eclosion, with diploid male patterns diverging from those of queens to resemble those of haploid males, a pattern retained in older adults. The finding that ploidy level effects on early gene expression override sex effects (including genes implicated in sperm production and pheromone production/perception) may explain diploid male sterility and lack of worker discrimination against them during development.