Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Combination of fluorescent reporters for simultaneous monitoring of root colonization and antifungal gene expression by a biocontrol pseudomonad on cereals with flow cytometry.

Some root-associated pseudomonads sustain plant growth by suppressing root diseases caused by pathogenic fungi. We investigated to which extent select cereal cultivars influence expression of relevant biocontrol traits (i.e., root colonization efficacy and antifungal activity) in Pseudomonas fluorescens CHA0. In this representative plant-beneficial bacterium, the antifungal metabolites 2,4-diacetylphloroglucinol (DAPG), pyrrolnitrin (PRN), pyoluteorin (PLT), and hydrogen cyanide (HCN) are required for biocontrol. To monitor host plant effects on the expression of biosynthetic genes for these compounds on roots, we developed fluorescent dual-color reporters suited for flow cytometric analysis using fluorescence-activated cell sorting (FACS). In the dual-label strains, the constitutively expressed red fluorescent protein mCherry served as a cell tag and marker for root colonization, whereas reporter fusions based on the green fluorescent protein allowed simultaneous recording of antifungal gene expression within the same cell. FACS analysis revealed that expression of DAPG and PRN biosynthetic genes was promoted in a cereal rhizosphere, whereas expression of PLT and HCN biosynthetic genes was markedly less sustained. When analyzing the response of the bacterial reporters on roots of a selection of wheat, spelt, and triticale cultivars, we were able to detect subtle species- and cultivar-dependent differences in colonization and DAPG and HCN gene expression levels. The expression of these biocontrol traits was particularly favored on roots of one spelt cultivar, suggesting that a careful choice of pseudomonad-cereal combinations might be beneficial to biocontrol. Our approach may be useful for selective single-cell level analysis of plant effects in other bacteria-root interactions.

Test of four colon cancer risk-scores in formalin fixed paraffin embedded microarray gene expression data.

BACKGROUND: Prognosis prediction for resected primary colon cancer is based on the T-stage Node Metastasis (TNM) staging system. We investigated if four well-documented gene expression risk scores can improve patient stratification. METHODS: Microarray-based versions of risk-scores were applied to a large independent cohort of 688 stage II/III tumors from the PETACC-3 trial. Prognostic value for relapse-free survival (RFS), survival after relapse (SAR), and overall survival (OS) was assessed by regression analysis. To assess improvement over a reference, prognostic model was assessed with the area under curve (AUC) of receiver operating characteristic (ROC) curves. All statistical tests were two-sided, except the AUC increase. RESULTS: All four risk scores (RSs) showed a statistically significant association (single-test, P < .0167) with OS or RFS in univariate models, but with HRs below 1.38 per interquartile range. Three scores were predictors of shorter RFS, one of shorter SAR. Each RS could only marginally improve an RFS or OS model with the known factors T-stage, N-stage, and microsatellite instability (MSI) status (AUC gains < 0.025 units). The pairwise interscore discordance was never high (maximal Spearman correlation = 0.563) A combined score showed a trend to higher prognostic value and higher AUC increase for OS (HR = 1.74, 95% confidence interval [CI] = 1.44 to 2.10, P < .001, AUC from 0.6918 to 0.7321) and RFS (HR = 1.56, 95% CI = 1.33 to 1.84, P < .001, AUC from 0.6723 to 0.6945) than any single score. CONCLUSIONS: The four tested gene expression-based risk scores provide prognostic information but contribute only marginally to improving models based on established risk factors. A combination of the risk scores might provide more robust information. Predictors of RFS and SAR might need to be different.

SOST expression is restricted to the great arteries during embryonic and neonatal cardiovascular development.

Spatial-temporal regulation of bone morphogenetic protein (BMP) and Wnt activity is essential for normal cardiovascular development, and altered activity of these growth factors causes maldevelopment of the cardiac outflow tract and great arteries. In the present study, we show that SOST, a Dan family member reported to antagonize BMP and Wnt activity, is expressed within the medial vessel wall of the great arteries containing smooth muscle cells. The ascending aorta, aortic arch, brachiocephalic artery, common carotids, and pulmonary trunk were all associated with SOST expressing smooth muscle cells, while the heart itself, including the valves, and more distal arteries, that is, pulmonary arteries, subclavian arteries, and descending aorta, were negative. SOST was expressed from embryonic day 15.5 up to the neonatal period. SOST expression, however, did not correspond with inhibition of Smad-dependent BMP activity or beta-catenin-dependent Wnt activity in the great arteries. Activity of both signaling pathways was already down-regulated before induction of SOST expression.

MAR elements regulate the probability of epigenetic switching between active and inactive gene expression.

Gene expression often cycles between active and inactive states in eukaryotes, yielding variable or noisy gene expression in the short-term, while slow epigenetic changes may lead to silencing or variegated expression. Understanding how cells control these effects will be of paramount importance to construct biological systems with predictable behaviours. Here we find that a human matrix attachment region (MAR) genetic element controls the stability and heritability of gene expression in cell populations. Mathematical modeling indicated that the MAR controls the probability of long-term transitions between active and inactive expression, thus reducing silencing effects and increasing the reactivation of silent genes. Single-cell short-terms assays revealed persistent expression and reduced expression noise in MAR-driven genes, while stochastic burst of expression occurred without this genetic element. The MAR thus confers a more deterministic behavior to an otherwise stochastic process, providing a means towards more reliable expression of engineered genetic systems.

Reactive electrophile species activate defense gene expression in Arabidopsis.

Compounds containing alpha,beta-unsaturated carbonyl groups are increasingly implicated as potent regulators of gene expression; some are powerful cytotoxins known to accumulate at the site of lesion formation in host-pathogen interactions. We used a robust measurement of photosynthetic efficiency to quantify the toxicity of a variety of lipid derivatives in Arabidopsis leaves. Small alpha,beta-unsaturated carbonyl compounds (e.g. acrolein and methyl vinyl ketone) were highly active and proved to be potent stimulators of expression of the pathogenesis-related gene HEL (PR4). These small volatile electrophiles were far more active than larger alkenal homologs like 2(E)-hexenal, and activated HEL expression in a manner independent of salicylate, ethylene, and jasmonate production/perception. Electrophile treatment massively increased the levels of unesterified cyclopentenone jasmonates, which themselves are electrophiles. Patterns of gene expression in response to electrophile treatment and in response to avirulent bacteria were compared, which revealed strikingly similar transcript profiles. The results broaden the range of known biologic effects of reactive electrophile species to include the activation of a pathogenesis-related gene (HEL) and genes involved in metabolism. Electrophiles can act as mediators of both genetic and biochemical effects on core defense signal transduction.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.

aIr-1, a newly found locus on mouse chromosome 16 encoding a trans-acting activator factor for MHC class II gene expression.

RJ 2.2.5 is a human B cell line that has lost the capacity to express MHC class II genes. The human class II-positive phenotype is restored in somatic cell hybrids between RJ 2.2.5 and mouse spleen cells. By karyotype and molecular studies of an informative family of hybrids we have now shown that the reexpression of human class II gene products, as well as the maintenance of the mouse class II-positive phenotype, correlates with the presence of mouse chromosome 16. Thus, the existence on this mouse chromosome of a newly found locus, designated by us aIr-1, that determines a trans-acting activator function for class II gene expression, is established. Possible implications of this finding are discussed.

Critical role for Ets, AP-1 and GATA-like transcription factors in regulating mouse Toll-like receptor 4 (Tlr4) gene expression.

TLR4 (Toll-like receptor 4) is essential for sensing the endotoxin of Gram-negative bacteria. Mutations or deletion of the TLR4 gene in humans or mice have been associated with altered predisposition to or outcome of Gram-negative sepsis. In the present work, we studied the expression and regulation of the Tlr4 gene of mouse. In vivo, TLR4 levels were higher in macrophages compared with B, T or natural killer cells. High basal TLR4 promoter activity was observed in RAW 264.7, J774 and P388D1 macrophages transfected with a TLR4 promoter reporter vector. Analysis of truncated and mutated promoter constructs identified several positive [two Ets (E twenty-six) and one AP-1 (activator protein-1) sites] and negative (a GATA-like site and an octamer site) regulatory elements within 350 bp upstream of the transcriptional start site. The myeloid and B-cell-specific transcription factor PU.1 bound to the proximal Ets site. In contrast, none among PU.1, Ets-1, Ets-2 and Elk-1, but possibly one member of the ESE (epithelium-specific Ets) subfamily of Ets transcription factors, bound to the distal Ets site, which was indispensable for Tlr4 gene transcription. Endotoxin did not affect macrophage TLR4 promoter activity, but it decreased TLR4 steady-state mRNA levels by increasing the turnover of TLR4 transcripts. TLR4 expression was modestly altered by other pro- and anti-inflammatory stimuli, except for PMA plus ionomycin which strongly increased promoter activity and TLR4 mRNA levels. The mouse and human TLR4 genes were highly conserved. Yet, notable differences exist with respect to the elements implicated in gene regulation, which may account for species differences in terms of tissue expression and modulation by microbial and inflammatory stimuli.

Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures.

INTRODUCTION: Breast cancer subtyping and prognosis have been studied extensively by gene expression profiling, resulting in disparate signatures with little overlap in their constituent genes. Although a previous study demonstrated a prognostic concordance among gene expression signatures, it was limited to only one dataset and did not fully elucidate how the different genes were related to one another nor did it examine the contribution of well-known biological processes of breast cancer tumorigenesis to their prognostic performance. METHOD: To address the above issues and to further validate these initial findings, we performed the largest meta-analysis of publicly available breast cancer gene expression and clinical data, which are comprised of 2,833 breast tumors. Gene coexpression modules of three key biological processes in breast cancer (namely, proliferation, estrogen receptor [ER], and HER2 signaling) were used to dissect the role of constituent genes of nine prognostic signatures. RESULTS: Using a meta-analytical approach, we consolidated the signatures associated with ER signaling, ERBB2 amplification, and proliferation. Previously published expression-based nomenclature of breast cancer 'intrinsic' subtypes can be mapped to the three modules, namely, the ER-/HER2- (basal-like), the HER2+ (HER2-like), and the low- and high-proliferation ER+/HER2- subtypes (luminal A and B). We showed that all nine prognostic signatures exhibited a similar prognostic performance in the entire dataset. Their prognostic abilities are due mostly to the detection of proliferation activity. Although ER- status (basal-like) and ERBB2+ expression status correspond to bad outcome, they seem to act through elevated expression of proliferation genes and thus contain only indirect information about prognosis. Clinical variables measuring the extent of tumor progression, such as tumor size and nodal status, still add independent prognostic information to proliferation genes. CONCLUSION: This meta-analysis unifies various results of previous gene expression studies in breast cancer. It reveals connections between traditional prognostic factors, expression-based subtyping, and prognostic signatures, highlighting the important role of proliferation in breast cancer prognosis.

Na(+)-K(+)-ATPase gene expression during in vitro development of rat fetal forebrain.

The existence of at least three isoforms of Na(+)-K(+)-ATPase in adult brain tissues [alpha 1, kidney type; alpha 2 [or alpha(+)]; alpha 3] suggests that these genes might be regulated in a cell-specific and time-dependent manner during development. We have studied this question in serum-free aggregating cell cultures of mechanically dissociated rat fetal telencephalon. At the protein level, the relative rate of synthesis of the pool of alpha 1-, alpha 2-, and alpha 3-subunits increased approximately twofold over 15 days of culture, leading to a marked increase in the immunochemical pool of alpha-subunits as measured by a panspecific polyclonal antibody. Concomitantly, Na(+)-K(+)-ATPase enzyme-specific activity increased three- (lower forebrain) to sixfold (upper forebrain). The transcripts of all three alpha-isoforms and beta-subunit were detected in vitro in similar proportion to the level observed in vivo. alpha 3-mRNA (3.7 kb) was more abundant than alpha 1 (3.7 kb) or alpha 2 (5.3 and 3.4 kb). Cytosine arabinoside (0.4 microM) and cholera toxin (0.1 microM) were used to selectively eliminate glial cells or neurons, respectively. It was found that alpha 2-mRNA is predominantly transcribed in glial cell cultures, whereas alpha 3- and beta 1-mRNA (2.7, 2.3, and 1.8 kb) are predominant in neuronal cultures.

Detection of plant-modulated alterations in antifungal gene expression in Pseudomonas fluorescens CHA0 on roots by flow cytometry.

The biocontrol activity of the root-colonizing Pseudomonas fluorescens strain CHA0 is largely determined by the production of antifungal metabolites, especially 2,4-diacetylphloroglucinol. The expression of these metabolites depends on abiotic and biotic environmental factors, in particular, elements present in the rhizosphere. In this study, we have developed a new method for the in situ analysis of antifungal gene expression using flow cytometry combined with green fluorescent protein (GFP)-based reporter fusions to the phlA and prnA genes essential for the production of the antifungal compounds 2,4-diacetylphloroglucinol and pyrrolnitrin, respectively, in strain CHA0. Expression of phlA-gfp and prnA-gfp in CHA0 cells harvested from the rhizosphere of a set of plant species as well as from the roots of healthy, leaf pathogen-attacked, and physically stressed plants were analyzed using a FACSCalibur. After subtraction of background fluorescence emitted by plant-derived particles and CHA0 cells not carrying the gfp reporters, the average gene expression per bacterial cell could be calculated. Levels of phlA and prnA expression varied significantly in the rhizospheres of different plant species. Physical stress and leaf pathogen infection lowered phlA expression levels in the rhizosphere of cucumber. Our results demonstrate that the newly developed approach is suitable to monitor differences in levels of antifungal gene expression in response to various plant-derived factors. An advantage of the method is that it allows quantification of bacterial gene expression in rhizosphere populations at a single-cell level. To our best knowledge, this is the first study using flow cytometry for the in situ analysis of biocontrol gene expression in a plant-beneficial bacterium in the rhizosphere.

The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

Fas and Fas ligand in embryos and adult mice: ligand expression in several immune-privileged tissues and coexpression in adult tissues characterized by apoptotic cell turnover.

The cell surface receptor Fas (FasR, Apo-1, CD95) and its ligand (FasL) are mediators of apoptosis that have been shown to be implicated in the peripheral deletion of autoimmune cells, activation-induced T cell death, and one of the two major cytolytic pathways mediated by CD8+ cytolytic T cells. To gain further understanding of the Fas system., we have analyzed Fas and FasL expression during mouse development and in adult tissues. In developing mouse embryos, from 16.5 d onwards, Fas mRNA is detectable in distinct cell types of the developing sinus, thymus, lung, and liver, whereas FasL expression is restricted to submaxillary gland epithelial cells and the developing nervous system. Significant Fas and FasL expression were observed in several nonlymphoid cell types during embryogenesis, and generally Fas and FasL expression were not localized to characteristic sites of programmed cell death. In the adult mouse, RNase protection analysis revealed very wide expression of both Fas and FasL. Several tissues, including the thymus, lung, spleen, small intestine, large intestine, seminal vesicle, prostate, and uterus, clearly coexpress the two genes. Most tissues constitutively coexpressing Fas and FasL in the adult mouse are characterized by apoptotic cell turnover, and many of those expressing FasL are known to be immune privileged. It may be, therefore, that the Fas system is implicated in both the regulation of physiological cell turnover and the protection of particular tissues against potential lymphocyte-mediated damage.