Arboretum: reconstruction and analysis of the evolutionary history of condition-specific transcriptional modules.

Comparative functional genomics studies the evolution of biological processes by analyzing functional data, such as gene expression profiles, across species. A major challenge is to compare profiles collected in a complex phylogeny. Here, we present Arboretum, a novel scalable computational algorithm that integrates expression data from multiple species with species and gene phylogenies to infer modules of coexpressed genes in extant species and their evolutionary histories. We also develop new, generally applicable measures of conservation and divergence in gene regulatory modules to assess the impact of changes in gene content and expression on module evolution. We used Arboretum to study the evolution of the transcriptional response to heat shock in eight species of Ascomycota fungi and to reconstruct modules of the ancestral environmental stress response (ESR). We found substantial conservation in the stress response across species and in the reconstructed components of the ancestral ESR modules. The greatest divergence was in the most induced stress, primarily through module expansion. The divergence of the heat stress response exceeds that observed in the response to glucose depletion in the same species. Arboretum and its associated analyses provide a comprehensive framework to systematically study regulatory evolution of condition-specific responses.

Bypassing Iron Storage in Endodermal Vacuoles Rescues the Iron Mobilization Defect in the natural resistance associated-macrophage protein3natural resistance associated-macrophage protein4 Double Mutant.

To improve seed iron (Fe) content and bioavailability, it is crucial to decipher the mechanisms that control Fe storage during seed development. In Arabidopsis (Arabidopsis thaliana) seeds, most Fe is concentrated in insoluble precipitates, with phytate in the vacuoles of cells surrounding the vasculature of the embryo. NATURAL RESISTANCE ASSOCIATED-MACROPHAGE PROTEIN3 (AtNRAMP3) and AtNRAMP4 function redundantly in Fe retrieval from vacuoles during germination. When germinated under Fe-deficient conditions, development of the nramp3nramp4 double mutant is arrested as a consequence of impaired Fe mobilization. To identify novel genes involved in seed Fe homeostasis, we screened an ethyl methanesulfonate-mutagenized population of nramp3nramp4 seedlings for mutations suppressing their phenotypes on low Fe. Here, we report that, among the suppressors, two independent mutations in the VACUOLAR IRON TRANSPORTER1 (AtVIT1) gene caused the suppressor phenotype. The AtVIT1 transporter is involved in Fe influx into vacuoles of endodermal and bundle sheath cells. This result establishes a functional link between Fe loading in vacuoles by AtVIT1 and its remobilization by AtNRAMP3 and AtNRAMP4. Moreover, analysis of subcellular Fe localization indicates that simultaneous disruption of AtVIT1, AtNRAMP3, and AtNRAMP4 limits Fe accumulation in vacuolar globoids.

The role of nucleosome positioning in the evolution of gene regulation.

Chromatin organization plays a major role in gene regulation and can affect the function and evolution of new transcriptional programs. However, it can be difficult to decipher the basis of changes in chromatin organization and their functional effect on gene expression. Here, we present a large-scale comparative genomic analysis of the relationship between chromatin organization and gene expression, by measuring mRNA abundance and nucleosome positions genome-wide in 12 Hemiascomycota yeast species. We found substantial conservation of global and functional chromatin organization in all species, including prominent nucleosome-free regions (NFRs) at gene promoters, and distinct chromatin architecture in growth and stress genes. Chromatin organization has also substantially diverged in both global quantitative features, such as spacing between adjacent nucleosomes, and in functional groups of genes. Expression levels, intrinsic anti-nucleosomal sequences, and trans-acting chromatin modifiers all play important, complementary, and evolvable roles in determining NFRs. We identify five mechanisms that couple chromatin organization to evolution of gene regulation and have contributed to the evolution of respiro-fermentation and other key systems, including (1) compensatory evolution of alternative modifiers associated with conserved chromatin organization, (2) a gradual transition from constitutive to trans-regulated NFRs, (3) a loss of intrinsic anti-nucleosomal sequences accompanying changes in chromatin organization and gene expression, (4) re-positioning of motifs from NFRs to nucleosome-occluded regions, and (5) the expanded use of NFRs by paralogous activator-repressor pairs. Our study sheds light on the molecular basis of chromatin organization, and on the role of chromatin organization in the evolution of gene regulation.

Gene duplication and the evolution of ribosomal protein gene regulation in yeast.

Coexpression of genes within a functional module can be conserved at great evolutionary distances, whereas the associated regulatory mechanisms can substantially diverge. For example, ribosomal protein (RP) genes are tightly coexpressed in Saccharomyces cerevisiae, but the cis and trans factors associated with them are surprisingly diverged across Ascomycota fungi. Little is known, however, about the functional impact of such changes on actual expression levels or about the selective pressures that affect them. Here, we address this question in the context of the evolution of the regulation of RP gene expression by using a comparative genomics approach together with cross-species functional assays. We show that an activator (Ifh1) and a repressor (Crf1) that control RP gene regulation in normal and stress conditions in S. cerevisiae are derived from the duplication and subsequent specialization of a single ancestral protein. We provide evidence that this regulatory innovation coincides with the duplication of RP genes in a whole-genome duplication (WGD) event and may have been important for tighter control of higher levels of RP transcripts. We find that subsequent loss of the derived repressor led to the loss of a stress-dependent repression of RPs in the fungal pathogen Candida glabrata. Our comparative computational and experimental approach shows how gene duplication can constrain and drive regulatory evolution and provides a general strategy for reconstructing the evolutionary trajectory of gene regulation across species.

The Arabidopsis MTP8 transporter determines the localization of manganese and iron in seeds.

Understanding how seeds obtain and store nutrients is key to developing crops with higher agronomic and nutritional value. We have uncovered unique patterns of micronutrient localization in seeds using synchrotron X-ray fluorescence (SXRF). Although all four members of the Arabidopsis thaliana Mn-CDF family can transport Mn, here we show that only mtp8-2 has an altered Mn distribution pattern in seeds. In an mtp8-2 mutant, Mn no longer accumulates in hypocotyl cortex cells and sub-epidermal cells of the embryonic cotyledons, but rather accumulates with Fe in the cells surrounding the vasculature, a pattern previously shown to be determined by the vacuolar transporter VIT1. We also show that MTP8, unlike the other three Mn-CDF family members, can transport Fe and is responsible for localization of Fe to the same cells that store Mn. When both the VIT1 and MTP8 transporters are non-functional, there is no accumulation of Fe or Mn in specific cell types; rather these elements are distributed amongst all cell types in the seed. Disruption of the putative Fe binding sites in MTP8 resulted in loss of ability to transport Fe but did not affect the ability to transport Mn.

Inference and Evolutionary Analysis of Genome-Scale Regulatory Networks in Large Phylogenies.

Changes in transcriptional regulatory networks can significantly contribute to species evolution and adaptation. However, identification of genome-scale regulatory networks is an open challenge, especially in non-model organisms. Here, we introduce multi-species regulatory network learning (MRTLE), a computational approach that uses phylogenetic structure, sequence-specific motifs, and transcriptomic data, to infer the regulatory networks in different species. Using simulated data from known networks and transcriptomic data from six divergent yeasts, we demonstrate that MRTLE predicts networks with greater accuracy than existing methods because it incorporates phylogenetic information. We used MRTLE to infer the structure of the transcriptional networks that control the osmotic stress responses of divergent, non-model yeast species and then validated our predictions experimentally. Interrogating these networks reveals that gene duplication promotes network divergence across evolution. Taken together, our approach facilitates study of regulatory network evolutionary dynamics across multiple poorly studied species.

Mn-euvering manganese: the role of transporter gene family members in manganese uptake and mobilization in plants.

Manganese (Mn), an essential trace element, is important for plant health. In plants, Mn serves as a cofactor in essential processes such as photosynthesis, lipid biosynthesis and oxidative stress. Mn deficient plants exhibit decreased growth and yield and are more susceptible to pathogens and damage at freezing temperatures. Mn deficiency is most prominent on alkaline soils with approximately one third of the world's soils being too alkaline for optimal crop production. Despite the importance of Mn in plant development, relatively little is known about how it traffics between plant tissues and into and out of organelles. Several gene transporter families have been implicated in Mn transport in plants. These transporter families include NRAMP (natural resistance associated macrophage protein), YSL (yellow stripe-like), ZIP (zinc regulated transporter/iron-regulated transporter [ZRT/IRT1]-related protein), CAX (cation exchanger), CCX (calcium cation exchangers), CDF/MTP (cation diffusion facilitator/metal tolerance protein), P-type ATPases and VIT (vacuolar iron transporter). A combination of techniques including mutant analysis and Synchrotron X-ray Fluorescence Spectroscopy can assist in identifying essential transporters of Mn. Such knowledge would vastly improve our understanding of plant Mn homeostasis.

Methodological approaches for using synchrotron X-ray fluorescence (SXRF) imaging as a tool in ionomics: examples from Arabidopsis thaliana.

Here we present approaches for using multi-elemental imaging (specifically synchrotron X-ray fluorescence microscopy, SXRF) in ionomics, with examples using the model plant Arabidopsis thaliana. The complexity of each approach depends on the amount of a priori information available for the gene and/or phenotype being studied. Three approaches are outlined, which apply to experimental situations where a gene of interest has been identified but has an unknown phenotype (phenotyping), an unidentified gene is associated with a known phenotype (gene cloning) and finally, a screening approach, where both gene and phenotype are unknown. These approaches make use of open-access, online databases with which plant molecular genetics researchers working in the model plant Arabidopsis will be familiar, in particular the Ionomics Hub and online transcriptomic databases such as the Arabidopsis eFP browser. The approaches and examples we describe are based on the assumption that altering the expression of ion transporters can result in changes in elemental distribution. We provide methodological details on using elemental imaging to aid or accelerate gene functional characterization by narrowing down the search for candidate genes to the tissues in which elemental distributions are altered. We use synchrotron X-ray microprobes as a technique of choice, which can now be used to image all parts of an Arabidopsis plant in a hydrated state. We present elemental images of leaves, stem, root, siliques and germinating hypocotyls.

Evolutionary principles of modular gene regulation in yeasts.

Divergence in gene regulation can play a major role in evolution. Here, we used a phylogenetic framework to measure mRNA profiles in 15 yeast species from the phylum Ascomycota and reconstruct the evolution of their modular regulatory programs along a time course of growth on glucose over 300 million years [corrected]. We found that modules have diverged proportionally to phylogenetic distance, with prominent changes in gene regulation accompanying changes in lifestyle and ploidy, especially in carbon metabolism. Paralogs have significantly contributed to regulatory divergence, typically within a very short window from their duplication. Paralogs from a whole genome duplication (WGD) event have a uniquely substantial contribution that extends over a longer span. Similar patterns occur when considering the evolution of the heat shock regulatory program measured in eight of the species, suggesting that these are general evolutionary principles. DOI:http://dx.doi.org/10.7554/eLife.00603.001.

Correction: Evolutionary principles of modular gene regulation in yeasts.

[This corrects the article DOI: 10.7554/eLife.00603.].

Determinants of divergent adaptation and Dobzhansky-Muller interaction in experimental yeast populations.

Divergent adaptation can be associated with reproductive isolation in speciation [1]. We recently demonstrated the link between divergent adaptation and the onset of reproductive isolation in experimental populations of the yeast Saccharomyces cerevisiae evolved from a single progenitor in either a high-salt or a low-glucose environment [2]. Here, whole-genome resequencing and comparative genome hybridization of representatives of three populations revealed 17 mutations, six of which explained the adaptive increases in mitotic fitness. In two populations evolved in high salt, two different mutations occurred in the proton efflux pump gene PMA1 and the global transcriptional repressor gene CYC8; the ENA genes encoding sodium efflux pumps were overexpressed once through expansion of this gene cluster and once because of mutation in the regulator CYC8. In the population from low glucose, one mutation occurred in MDS3, which modulates growth at high pH, and one in MKT1, a global regulator of mRNAs encoding mitochondrial proteins, the latter recapitulating a naturally occurring variant. A Dobzhansky-Muller (DM) incompatibility between the evolved alleles of PMA1 and MKT1 strongly depressed fitness in the low-glucose environment. This DM interaction is the first reported between experimentally evolved alleles of known genes and shows how reproductive isolation can arise rapidly when divergent selection is strong.

An erythropoietin autocrine/paracrine axis modulates the growth and survival of human prostate cancer cells.

Erythropoietin receptors have been identified on a variety of cancer-derived cell lines and primary cancer cells, including those of prostate cancer. The functional status of these extrahematopoietic erythropoietin receptors remains a matter of some dispute. The publication of several important clinical trials suggesting a direct effect of erythropoietin on the growth and survival of primary tumors adds further importance to the question of whether erythropoietin receptors on cancer cells are functional. We have reported previously that human prostate cancer cell lines and primary prostate cancer cells express functional erythropoietin receptors that respond to exogenous erythropoietin by increased cell proliferation and STAT5 phosphorylation. We now show that prostate cancer cell lines express both the EPO gene and the biologically active erythropoietin. The coexpression of functional receptor and biologically active ligand in the cells has led us to hypothesize an autocrine/paracrine mechanism, driven by endogenous erythropoietin, which may modulate the growth and progression of prostate cancer. To test our hypothesis, we have knocked down, independently, erythropoietin receptor and erythropoietin on prostate cancer cells by transfection with short hairpin RNAs. Erythropoietin receptor knockdown cells grow significantly more slowly than their erythropoietin receptor-bearing counterparts in monolayer culture, produce fewer, smaller colonies in soft agar, and do not exhibit erythropoietin-induced signaling. Erythropoietin knockdown cells exhibit dramatically slower rates of growth, which could be restored by transfecting the cells with a murine erythropoietin gene. Taken together, our data suggest that the coordinated regulation of a functional erythropoietin/erythropoietin receptor axis in prostate cancer cells may be integral to the growth and progression of prostate cancer.