Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis.

Wounding in multicellular eukaryotes results in marked changes in gene expression that contribute to tissue defense and repair. Using a cDNA microarray technique, we analyzed the timing, dynamics, and regulation of the expression of 150 genes in mechanically wounded leaves of Arabidopsis. Temporal accumulation of a group of transcripts was correlated with the appearance of oxylipin signals of the jasmonate family. Analysis of the coronatine-insensitive coi1-1 Arabidopsis mutant that is also insensitive to jasmonate allowed us to identify a large number of COI1-dependent and COI1-independent wound-inducible genes. Water stress was found to contribute to the regulation of an unexpectedly large fraction of these genes. Comparing the results of mechanical wounding with damage by feeding larvae of the cabbage butterfly (Pieris rapae) resulted in very different transcript profiles. One gene was specifically induced by insect feeding but not by wounding; moreover, there was a relative lack of water stress-induced gene expression during insect feeding. These results help reveal a feeding strategy of P. rapae that may minimize the activation of a subset of water stress-inducible, defense-related genes.

Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes

Mammalian sex chromosomes stem from ancestral autosomes and have substantially differentiated. It was shown that X-linked genes have generated duplicate intronless gene copies (retrogenes) on autosomes due to this differentiation. However, the precise driving forces for this out-of-X gene "movement" and its evolutionary onset are not known. Based on expression analyses of male germ-cell populations, we here substantiate and extend the hypothesis that autosomal retrogenes functionally compensate for the silencing of their X-linked housekeeping parental genes during, but also after, male meiotic sex chromosome inactivation (MSCI). Thus, sexually antagonistic forces have not played a major role for the selective fixation of X-derived gene copies in mammals. Our dating analyses reveal that although retrogenes were produced ever since the common mammalian ancestor, selectively driven retrogene export from the X only started later, on the placental mammal (eutherian) and marsupial (metatherian) lineages, respectively. Together, these observations suggest that chromosome-wide MSCI emerged close to the eutherian-marsupial split approximately 180 million years ago. Given that MSCI probably reflects the spread of the recombination barrier between the X and Y, crucial for their differentiation, our data imply that these chromosomes became more widely differentiated only late in the therian ancestor, well after the divergence of the monotreme lineage. Thus, our study also provides strong independent support for the recent notion that our sex chromosomes emerged, not in the common ancestor of all mammals, but rather in the therian ancestor, and therefore are much younger than previously thought

Inducible gene expression in fetal thymic epithelium: A new BAC transgenic model.

The thymus is the site of T cell development. Several stromal and hematopoietic cell types are necessary for the proper function of thymic selection and eventually peripheral immunity. Thymic epithelial cells (TECs) are essential for T cell lineage commitment, expansion, and maturation in the thymus. We were interested in developing an in vivo model in which exogenous gene expression could be transiently induced in embryonic TEC (Tet-On system). To this end, we have generated a bacterial artificial chromosome (BAC) transgenic mouse line in which the reverse tetracycline-dependent transactivator (rtTA) is expressed under the control of the Foxn1 promoter, a transcriptional factor indispensable for TEC development. To analyze the expression pattern and efficiency of this novel mouse model, we crossed the Foxn1-rtTA founder with a Tet-Responsive Element (TRE)-LacZ GFP mouse reporter to obtain a double transgenic mouse. In the presence of doxycycline, rtTA can interact with TRE and induce the expression of GFP and LacZ. In this double transgenic mouse, we observed that GFP expression was high, inducible and limited to TEC in fetal thymus. In contrast, in adult thymus, when TEC development and maturation is completed, GFP was barely detectable. Therefore, Foxn1-rtTA represents a new and efficient transgenic mouse model to induce genes of interest specifically in fetal thymic epithelium. genesis 51:717-724. © 2013 Wiley Periodicals, Inc.

Steady-state expression of self-regulated genes.

MOTIVATION: Regulatory gene networks contain generic modules such as feedback loops that are essential for the regulation of many biological functions. The study of the stochastic mechanisms of gene regulation is instrumental for the understanding of how cells maintain their expression at levels commensurate with their biological role, as well as to engineer gene expression switches of appropriate behavior. The lack of precise knowledge on the steady-state distribution of gene expression requires the use of Gillespie algorithms and Monte-Carlo approximations. METHODOLOGY: In this study, we provide new exact formulas and efficient numerical algorithms for computing/modeling the steady-state of a class of self-regulated genes, and we use it to model/compute the stochastic expression of a gene of interest in an engineered network introduced in mammalian cells. The behavior of the genetic network is then analyzed experimentally in living cells. RESULTS: Stochastic models often reveal counter-intuitive experimental behaviors, and we find that this genetic architecture displays a unimodal behavior in mammalian cells, which was unexpected given its known bimodal response in unicellular organisms. We provide a molecular rationale for this behavior, and we implement it in the mathematical picture to explain the experimental results obtained from this network.

Effects of mitotane on gene expression in the adrenocortical cell line NCI-H295R: a microarray study.

Aim: The adrenolytic agent mitotane is widely used in the treatment of adrenocortical cancer; however, its mechanism of action is poorly elucidated. We have studied mitotane-induced mRNA expression changes in the NCI-H295R adrenocortical cancer cell line. Materials & methods: Cell viability and hormone assays were used to select the optimal mitotane concentration effectively inhibiting hormone secretion without affecting cell viability. RNA isolated from cultures treated for 48 and 72 h was subjected to Agilent 4×44K microarray platforms. Microarray results were validated by quantitative reverse-transcription PCR. Results: Altogether, 117 significantly differentially expressed genes were detected at 48 h and 72 h (p < 0.05) in mitotane-treated samples relative to controls. Three significantly underexpressed genes involved in steroid hormone biosynthesis (HSD3B1, HSD3B2 and CYP21A2) and four significantly overexpressed genes (GDF15, ALDH1L2, TRIB3 and SERPINE2) have been validated. Conclusion: Gene-expression changes might be involved in the adrenal action of mitotane and in the inhibition of hormone secretion. Original submitted 20 January 2012; Revision submitted 17 May 2012.

Myelin gene expression during demyelination and remyelination in aggregating brain cell cultures.

Remyelination can be studied in aggregating rat brain cell cultures after limited demyelination. Demyelination was induced using a monoclonal antibody against myelin/oligodendrocyte glycoprotein (MOG mAb), in the presence of complement. De- and remyelination were assessed by measuring myelin basic protein (MBP). Two days after removing the MOG mAb, MBP levels reached 50% of controls and after 7 days 93%. During this period, cell proliferation determined by [14C]thymidine incorporation was similar in remyelinating and control cultures. Hormones and growth factors were tested for possible stimulatory effect on remyelinating cultures. Bovine growth hormone (bGH), triiodothyronine (T3), basic fibroblast growth factor (bFGF) and platelet-derived growth factor (PDGF) did not improve remyelination. Only epidermal growth factor (EGF) increased the level of remyelination. PDGF increased the rate of cell proliferation in both control and remyelinating cultures. A significant proportion of oligodendrocytes entered the cell division cycle and were not available for remyelination. The results obtained with PDGF and FGF (inhibition) support the idea that a pool of progenitor cells was still present and able to proliferate and differentiate into myelinating oligodendrocytes. The levels of myelin protein mRNAs were investigated during de- and remyelination. During demyelination, myelin protein mRNA levels decreased to approximately 50% of control cultures and returned to normal during remyelination. These preliminary results indicate that normal levels of gene transcription are sufficient to meet the increased need for newly synthesized myelin proteins during remyelination.

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Analysis of HIV-1 expression level and sense of transcription by high-throughput sequencing of the infected cell.

Next-generation sequencing offers an unprecedented opportunity to jointly analyze cellular and viral transcriptional activity without prerequisite knowledge of the nature of the transcripts. SupT1 cells were infected with a vesicular stomatitis virus G envelope protein (VSV-G)-pseudotyped HIV vector. At 24 h postinfection, both cellular and viral transcriptomes were analyzed by serial analysis of gene expression followed by high-throughput sequencing (SAGE-Seq). Read mapping resulted in 33 to 44 million tags aligning with the human transcriptome and 0.23 to 0.25 million tags aligning with the genome of the HIV-1 vector. Thus, at peak infection, 1 transcript in 143 is of viral origin (0.7%), including a small component of antisense viral transcription. Of the detected cellular transcripts, 826 (2.3%) were differentially expressed between mock- and HIV-infected samples. The approach also assessed whether HIV-1 infection modulates the expression of repetitive elements or endogenous retroviruses. We observed very active transcription of these elements, with 1 transcript in 237 being of such origin, corresponding on average to 123,123 reads in mock-infected samples (0.40%) and 129,149 reads in HIV-1-infected samples (0.45%) mapping to the genomic Repbase repository. This analysis highlights key details in the generation and interpretation of high-throughput data in the setting of HIV-1 cellular infection.

Evolution of gene expression changes in newborn rats after mechanical ventilation with reversible intubation.

Mechanical ventilation (MV) is life-saving but potentially harmful for lungs of premature infants. So far, animal models dealt with the acute impact of MV on immature lungs, but less with its delayed effects. We used a newborn rodent model including non-surgical and therefore reversible intubation with moderate ventilation and hypothesized that there might be distinct gene expression patterns after a ventilation-free recovery period compared to acute effects directly after MV. Newborn rat pups were subjected to 8 hr of MV with 60% oxygen (O(2) ), 24 hr after injection of lipopolysaccharide (LPS), intended to create a low inflammatory background as often recognized in preterm infants. Animals were separated in controls (CTRL), LPS injection (LPS), or full intervention with LPS and MV with 60% O(2) (LPS + MV + O(2) ). Lungs were recovered either directly following (T:0 hr) or 48 hr after MV (T:48 hr). Histologically, signs of ventilator-induced lung injury (VILI) were observed in LPS + MV + O(2) lungs at T:0 hr, while changes appeared similar to those known from patients with chronic lung disease (CLD) with fewer albeit larger gas exchange units, at T:48 hr. At T:0 hr, LPS + MV + O(2) increased gene expression of pro-inflammatory MIP-2. In parallel anti-inflammatory IL-1Ra gene expression was increased in LPS and LPS + MV + O(2) groups. At T:48 hr, pro- and anti-inflammatory genes had returned to their basal expression. MMP-2 gene expression was decreased in LPS and LPS + MV + O(2) groups at T:0 hr, but no longer at T:48 hr. MMP-9 gene expression levels were unchanged directly after MV. However, at T:48 hr, gene and protein expression increased in LPS + MV + O(2) group. In conclusion, this study demonstrates the feasibility of delayed outcome measurements after a ventilation-free period in newborn rats and may help to further understand the time-course of molecular changes following MV. The differences obtained from the two time points could be interpreted as an initial transitory increase of inflammation and a delayed impact of the intervention on structure-related genes. Pediatr Pulmonol. 2012; 47:1204-1214. © 2012 Wiley Periodicals, Inc.

Small nucleolar RNA expression profiling identifies potential prognostic markers in peripheral T-cell lymphoma.

Peripheral T-cell lymphoma (PTCL) is a rare, heterogeneous type of non-Hodgkin lymphoma (NHL) that, in general, is associated with a poor clinical outcome. Therefore, a current major challenge is the discovery of new prognostic tools for this disease. In the present study, a cohort of 122 patients with PTCL was collected from a multicentric T-cell lymphoma consortium (TENOMIC). We analyzed the expression of 80 small nucleolar RNAs (snoRNAs) using high-throughput quantitative PCR. We demonstrate that snoRNA expression analysis may be useful in both the diagnosis of some subtypes of PTCL and the prognostication of both PTCL-not otherwise specified (PTCL-NOS; n = 26) and angio-immunoblastic T-cell lymphoma (AITL; n = 46) patients treated with chemotherapy. Like miRNAs, snoRNAs are globally down-regulated in tumor cells compared with their normal counterparts. In the present study, the snoRNA signature was robust enough to differentiate anaplastic large cell lymphoma (n = 32) from other PTCLs. For PTCL-NOS and AITL, we obtained 2 distinct prognostic signatures with a reduced set of 3 genes. Of particular interest was the prognostic value of HBII-239 snoRNA, which was significantly over-expressed in cases of AITL and PTCL-NOS that had favorable outcomes. Our results suggest that snoRNA expression profiles may have a diagnostic and prognostic significance for PTCL, offering new tools for patient care and follow-up.

Steroid hormone pulsing drives cyclic gene expression.

Transcriptional cycling of activated glucocorticoid receptor (GR) and ultradian glucocorticoid secretion are well established processes. Ultradian hormone release is now shown to result in pulsatile gene transcription through dynamic exchange of GR with the target-gene promoter and GR cycling through the chaperone machinery.

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.