Evidence of neutral transcriptome evolution in plants

The transcriptome of an organism is its set of gene transcripts (mRNAs) at a defined spatial and temporal locus. Because gene expression is affected markedly by environmental and developmental perturbations, it is widely assumed that transcriptome divergence among taxa represents adaptive phenotypic selection. This assumption has been challenged by neutral theories which propose that stochastic processes drive transcriptome evolution. To test for evidence of neutral transcriptome evolution in plants, we quantified 18 494 gene transcripts in nonsenescent leaves of 14 taxa of Brassicaceae using robust cross-species transcriptomics which includes a two-step physical and in silico-based normalization procedure based on DNA similarity among taxa. Transcriptome divergence correlates positively with evolutionary distance between taxa and with variation in gene expression among samples. Results are similar for pseudogenes and chloroplast genes evolving at different rates. Remarkably, variation in transcript abundance among root-cell samples correlates positively with transcriptome divergence among root tissues and among taxa. Because neutral processes affect transcriptome evolution in plants, many differences in gene expression among or within taxa may be nonfunctional, reflecting ancestral plasticity and founder effects. Appropriate null models are required when comparing transcriptomes in space and time.

A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes

Whole-genome transcriptome profiling is revealing how biological systems are regulated at the transcriptional level. This study reports the development of a robust method to profile and compare the transcriptomes of two nonmodel plant species, Thlaspi caerulescens, a zinc (Zn) hyperaccumulator, and Thlaspi arvense, a nonhyperaccumulator, using Affymetrix Arabidopsis thaliana ATH1-121501 GeneChip (R) arrays (Affymetrix, Santa Clara, CA, USA). Transcript abundance was quantified in the shoots of agar- and compost-grown plants of both species. Analyses were optimized using a genomic DNA (gDNA)-based probe-selection strategy based on the hybridization efficiency of Thlaspi gDNA with corresponding A. thaliana probes. In silico alignments of GeneChip (R) probes with Thlaspi gene sequences, and quantitative real-time PCR, confirmed the validity of this approach. Approximately 5000 genes were differentially expressed in the shoots of T. caerulescens compared with T. arvense, including genes involved in Zn transport and compartmentalization. Future functional analyses of genes identified as differentially expressed in the shoots of these closely related species will improve our understanding of the molecular mechanisms of Zn hyperaccumulation.

Genetic responses to phosphorus deficiency

Background: Phosphorus (P) is an essential macronutrient for plants. Plants take up P as phosphate (Pi) from the soil solution. Since little Pi is available in most soils, P fertilizers are applied to crops. However, the use of P fertilizers is unsustainable and may cause pollution. Consequently, there is a need to develop more P-use-efficient (PUE) crops and precise methods to monitor crop P-status. Scope: Manipulating the expression of genes to improve the PUE of crops could reduce their P fertilizer requirement. This has stimulated research towards the identification of genes and signalling cascades involved in plant responses to P deficiency. Genes that respond to P deficiency can be grouped into 'early' genes that respond rapidly and often non-specifically to P deficiency, or 'late' genes that impact on the morphology, physiology or metabolism of plants upon Prolonged P deficiency. Summary: The use of micro-array technology has allowed researchers to catalogue the genetic responses of plants to P deficiency. Genes whose expression is altered by P deficiency include various transcription factors, which are thought to coordinate plant responses to P deficiency, and other genes involved in P acquisition and tissue P economy. Several common cis-regulatory elements have been identified in the promoters of these genes, suggesting that their expression might be coordinated. It is suggested that knowledge of the genes whose expression changes in response to P deficiency might allow the development of crops with improved PUE, and could be used in diagnostic techniques to monitor P deficiency in crops either directly using 'smart' indicator plants or indirectly through transcript profiling. The development of crops with improved PUE and the adoption of diagnostic technology could reduce production costs, minimize the use of a non-renewable resource, reduce pollution and enhance biodiversity.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Tackling drought stress: RECEPTOR-LIKE KINASES present new approaches

Global climate change and a growing population require tackling the reduction in arable land and improving biomass production and seed yield per area under varying conditions. One of these conditions is suboptimal water availability. Here, we review some of the classical approaches to dealing with plant response to drought stress and we evaluate how research on RECEPTOR-LIKE KINASES (RLKs) can contribute to improving plant performance under drought stress. RLKs are considered as key regulators of plant architecture and growth behavior, but they also function in defense and stress responses. The available literature and analyses of available transcript profiling data indeed suggest that RLKs can play an important role in optimizing plant responses to drought stress. In addition, RLK pathways are ideal targets for nontransgenic approaches, such as synthetic molecules, providing a novel strategy to manipulate their activity and supporting translational studies from model species, such as Arabidopsis thaliana, to economically useful crops.

Response to zinc deficiency of two rice lines with contrasting tolerance is determined by root growth maintenance and organic acid exudation rates, and not by zinc-transporter activity

Zinc (Zn)-deficient soils constrain rice (Oryza sativa) production and cause Zn malnutrition. The identification of Zn-deficiency-tolerant rice lines indicates that breeding might overcome these constraints. Here, we seek to identify processes underlying Zn-deficiency tolerance in rice at the physiological and transcriptional levels. A Zn-deficiency-tolerant line RIL46 acquires Zn more efficiently and produces more biomass than its nontolerant maternal line (IR74) at low Zn(ext) under field conditions. We tested if this was the result of increased expression of Zn(2+) transporters; increased root exudation of deoxymugineic acid (DMA) or low-molecular-weight organic acids (LMWOAs); and/or increased root production. Experiments were performed in field and controlled environment conditions. There was little genotypic variation in transcript abundance of Zn-responsive root Zn(2+)-transporters between the RIL46 and IR74. However, root exudation of DMA and LMWOA was greater in RIL46, coinciding with increased root expression of putative ligand-efflux genes. Adventitious root production was maintained in RIL46 at low Zn(ext), correlating with altered expression of root-specific auxin-responsive genes. Zinc-deficiency tolerance in RIL46 is most likely the result of maintenance of root growth, increased efflux of Zn ligands, and increased uptake of Zn-ligand complexes at low Zn(ext); these traits are potential breeding targets.

A novel plant-fungus symbiosis benefits the host without forming mycorrhizal structures

Most terrestrial plants form mutually beneficial symbioses with specific soil-borne fungi known as mycorrhiza. In a typical mycorrhizal association, fungal hyphae colonize plant roots, explore the soil beyond the rhizosphere and provide host plants with nutrients that might be chemically or physically inaccessible to root systems. Here, we combined nutritional, radioisotopic (33P) and genetic approaches to describe a plant growth promoting symbiosis between the basidiomycete fungus Austroboletus occidentalis and jarrah (Eucalyptus marginata), which has quite different characteristics. We show that the fungal partner does not colonize plant roots; hyphae are localized to the rhizosphere soil and vicinity and consequently do not transfer nutrients located beyond the rhizosphere. Transcript profiling of two high-affinity phosphate (Pi) transporter genes (EmPHT1;1 and EmPHT1;2) and hyphal-mediated 33Pi uptake suggest that the Pi uptake shifts from an epidermal to a hyphal pathway in ectomycorrhizal plants (Scleroderma sp.), similar to arbuscular mycorrhizal symbioses, whereas A. occidentalis benefits its host indirectly. The enhanced rhizosphere carboxylates are linked to growth and nutritional benefits in the novel symbiosis. This work is a starting point for detailed mechanistic studies on other basidiomycete–woody plant relationships, where a continuum between heterotrophic rhizosphere fungi and plant beneficial symbioses is likely to exist.

Effects of oxidative stress on the cardiac myocyte proteome: modifications to peroxiredoxins and small heat shock proteins

Endogenous oxidative stress is a likely cause of cardiac myocyte death in vivo. We examined the early (0-2 h) changes in the proteome of isolated cardiac myocytes from neonatal rats exposed to H2O2 (0.1 mM), focussing on proteins with apparent molecular masses of between 20 and 30 kDa. Proteins were separated by two-dimensional gel electrophoresis (2DGE), located by silver-staining and identified by mass spectrometry. Incorporation of [35S]methionine or 32Pi was also studied. For selected proteins, transcript abundance was examined by reverse transcriptase-polymerase chain reaction. Of the 38 protein spots in the region, 23 were identified. Two families showed changes in 2DGE migration or abundance with H2O2 treatment: the peroxiredoxins and two small heat shock protein (Hsp) family members: heat shock 27 kDa protein 1 (Hsp25) and alphaB-crystallin. Peroxiredoxins shifted to lower pI values and this was probably attributable to 'over-oxidation' of active site Cys-residues. Hsp25 also shifted to lower pI values but this was attributable to phosphorylation. alphaB-crystallin migration was unchanged but its abundance decreased. Transcripts encoding peroxiredoxins 2 and 5 increased significantly. In addition, 10 further proteins were identified. For two (glutathione S-transferase pi, translationally-controlled tumour protein), we could not find any previous references indicating their occurrence in cardiac myocytes. We conclude that exposure of cardiac myocytes to oxidative stress causes post-translational modification in two protein families involved in cytoprotection. These changes may be potentially useful diagnostically. In the short term, oxidative stress causes few detectable changes in global protein abundance as assessed by silver-staining.

Cardiac myocyte gene expression profiling during H2O2-induced apoptosis

High levels of oxidative stress promote cardiac myocyte death, though lower levels are potentially cytoprotective/anabolic. We examined the changes in gene expression in rat neonatal cardiac myocytes exposed to apoptotic (0.2 mM) or nontoxic (0.04 mM) concentrations of H2O2 (2, 4, or 24 h) using Affymetrix microarrays. Using U34B arrays, we identified a ubiquitously expressed, novel H2O2-responsive gene [putative peroxide-inducible transcript 1 (Perit1)], which generates two alternatively spliced transcripts. Using 230 2.0 arrays, H2O2 (0.04 mM) promoted significant changes in expression of only 32 genes, all of which were seen with 0.2 mM H2O2. We failed to detect any increase in the rate of protein synthesis in cardiac myocytes exposed to <0.1 mM H2O2, further suggesting that global, low concentrations of H2O2 are not anabolic in this system. H2O2 (0.2 mM) promoted significant (P < 0.05, >1.75-fold) changes in expression of 649 mRNAs and 187 RNAs corresponding to no established gene. Of the mRNAs, 114 encoded transcriptional regulators including Krüppel-like factors (Klfs). Quantitative PCR independently verified the changes in Klf expression. Thus, H2O2-induced cardiac myocyte apoptosis is associated with dynamic changes in gene expression. The expression of these genes and their protein products potentially influences the progression of the apoptotic response.

Soil bacterial phoD gene abundance and expression in response to applied phosphorus and long-term management

Bacterial transformation of phosphorus (P) compounds in soil is largely dependent on soil microbial community function, and is therefore sensitive to anthropogenic disturbances such as fertilization or cropping systems. However, the effect of soil management on the transcription of bacterial genes that encode phosphatases, such as phoD, is largely unknown. This greenhouse study examined the effect of long-term management and P amendment on potential alkaline phosphatase (ALP) activity and phoD gene (DNA) and transcript (RNA) abundance. Soil samples (0–15 cm) were collected from the Glenlea Long-term Rotation near Winnipeg, Manitoba, to compare organic, conventional and prairie management systems. In the greenhouse, pots of soil from each management system were amended with P as either soluble mineral fertilizer or cattle manure and then planted with Italian ryegrass (Lolium multiforum). Soils from each pot were sampled for analysis immediately and after 30 and 106 days. Significant differences among the soil/P treatments were detected for inorganic P, but not the organic P in NaHCO3-extracts. At day 0, ALP activity was similar among the soil/P treatments, but was higher after 30 days for all P amendments in soil from organically managed plots. In contrast, ALP activity in soils under conventional and prairie management responded to increasing rates of manure only, with significant effects from medium and high manure application rates at 30 and 106 days. Differences in ALP activity at 30 days corresponded to the abundance of bacterial phoD genes, which were also significantly higher in soils under organic management. However, this correlation was not significant for transcript abundance. Next-generation sequencing allowed the identification of 199 unique phoD operational taxonomic units (OTUs) from the metagenome (soil DNA) and 35 unique OTUs from the metatranscriptome (soil RNA), indicating that a subset of phoD genes was being transcribed in all soils.

Interaction of plant growth regulators and reactive oxygen species to regulate petal senescence in wallflowers (Erysimum linifolium)

Background In many species floral senescence is coordinated by ethylene. Endogenous levels rise, and exogenous application accelerates senescence. Furthermore, floral senescence is often associated with increased reactive oxygen species, and is delayed by exogenously applied cytokinin. However, how these processes are linked remains largely unresolved. Erysimum linifolium (wallflower) provides an excellent model for understanding these interactions due to its easily staged flowers and close taxonomic relationship to Arabidopsis. This has facilitated microarray analysis of gene expression during petal senescence and provided gene markers for following the effects of treatments on different regulatory pathways. Results In detached Erysimum linifolium (wallflower) flowers ethylene production peaks in open flowers. Furthermore senescence is delayed by treatments with the ethylene signalling inhibitor silver thiosulphate, and accelerated with ethylene released by 2-chloroethylphosphonic acid. Both treatments with exogenous cytokinin, or 6-methyl purine (which is an inhibitor of cytokinin oxidase), delay petal senescence. However, treatment with cytokinin also increases ethylene biosynthesis. Despite the similar effects on senescence, transcript abundance of gene markers is affected differentially by the treatments. A significant rise in transcript abundance of WLS73 (a putative aminocyclopropanecarboxylate oxidase) was abolished by cytokinin or 6-methyl purine treatments. In contrast, WFSAG12 transcript (a senescence marker) continued to accumulate significantly, albeit at a reduced rate. Silver thiosulphate suppressed the increase in transcript abundance both of WFSAG12 and WLS73. Activity of reactive oxygen species scavenging enzymes changed during senescence. Treatments that increased cytokinin levels, or inhibited ethylene action, reduced accumulation of hydrogen peroxide. Furthermore, although auxin levels rose with senescence, treatments that delayed early senescence did not affect transcript abundance of WPS46, an auxin-induced gene. Conclusions A model for the interaction between cytokinins, ethylene, reactive oxygen species and auxin in the regulation of floral senescence in wallflowers is proposed. The combined increase in ethylene and reduction in cytokinin triggers the initiation of senescence and these two plant growth regulators directly or indirectly result in increased reactive oxygen species levels. A fall in conjugated auxin and/or the total auxin pool eventually triggers abscission.

Conceptus-mediated endometrial vascular changes during early pregnancy in mares: an anatomic, histomorphometric, and vascular endothelial growth factor receptor system immunolocalization and gene expression study

This work examined how the conceptus modulates endometrial tissue remodeling and vascular development prior to implantation in mares. A macroscopic uterine examination was completed at day 21 of pregnancy. In situ morphology revealed that the endometrium involved in encroachment is restricted to the dorsal endometrium immediately overlying the yolk sac. The amount of stromal area occupied by blood vessels and the number of endometrial glands were increased during early pregnancy. Endometrial histomorphometry as well as the endometrial mRNA abundance and immunolocalization of VEGF, VEGFR1, VEGFR2, and Ki-67 was completed at days 14 and 21 of pregnancy, at day 10 of the estrous cycle, and during estrus. No obvious differences in VEGF and VEGFR1 protein localization were detected between pregnant and cycling mares but differential staining pattern for VEGFR2 and Ki-67 was observed. VEGFR2 localized to luminal and glandular epithelium of pregnant mares, while luminal epithelium was negative in cycling mares. Ki-67 staining was weak during the luteal phase but exhibited prominent luminal epithelium staining during estrus. In pregnant mares, all endometrial layers were Ki-67 positive. Quantitative RT-PCR revealed a greater abundance of VEGF mRNA during pregnancy. VEGFR2 transcript abundance was greatest in pregnant mares on day 21. This study supports the concept that the conceptus plays an active role in directing vasculogenesis within the uterus and thereby establishing hemotrophic nutrition that supports pregnancy after implantation. Reproduction (2011) 142 593-603

Axillary bud development in pineapple nodal segments correlates with changes on cell cycle gene expression, hormone level, and sucrose and glutamate contents

During the process of lateral organ development after plant decapitation, cell division and differentiation occur in a balanced manner initiated by specific signaling, which triggers the reentrance into the cell cycle. Here, we investigated short-term variations in the content of some endogenous signals, such as auxin, cytokinins (Cks), and other mitogenic stimuli (sucrose and glutamate), which are likely correlated with the cell cycle reactivation in the axillary bud primordium of pineapple nodal segments. Transcript levels of cell cycle-associated genes, CycD2;1, and histone H2A were analyzed. Nodal segments containing the quiescent axillary meristem cells were cultivated in vitro during 24 h after the apex removal and de-rooting. From the moment of stem apex and root removal, decapitated nodal segment (DNS) explants showed a lower indol-3-acetic acid (IAA) concentration than control explants, and soon after, an increase of endogenous sucrose and iP-type Cks were detected. The decrease of IAA may be the primary signal for cell cycle control early in G1 phase, leading to the upregulation of CycD2;1 gene in the first h. Later, the iP-type Cks and sucrose could have triggered the progression to S-phase since there was an increase in H2A expression at the eighth h. DNS explants revealed substantial increase in Z-type Cks and glutamate from the 12th h, suggesting that these mitogens could also operate in promoting pineapple cell cycle progression. We emphasize that the use of non-synchronized tissue rather than synchronous cell suspension culture makes it more difficult to interpret the results of a dynamic cell division process. However, pineapple nodal segments cultivated in vitro may serve as an interesting model to shed light on apical dominance release and the reentrance of quiescent axillary meristem cells into the cell cycle.

Gene Expression Profile of Mesenchymal Stem Cells from Paired Umbilical Cord Units: Cord is Different from Blood

Mesenchymal stem cells (MSC) are multipotent cells which can be obtained from several adult and fetal tissues including human umbilical cord units. We have recently shown that umbilical cord tissue (UC) is richer in MSC than umbilical cord blood (UCB) but their origin and characteristics in blood as compared to the cord remains unknown. Here we compared, for the first time, the exonic protein-coding and intronic noncoding RNA (ncRNA) expression profiles of MSC from match-paired UC and UCB samples, harvested from the same donors, processed simultaneously and under the same culture conditions. The patterns of intronic ncRNA expression in MSC from UC and UCB paired units were highly similar, indicative of their common donor origin. The respective exonic protein-coding transcript expression profiles, however, were significantly different. Hierarchical clustering based on protein-coding expression similarities grouped MSC according to their tissue location rather than original donor. Genes related to systems development, osteogenesis and immune system were expressed at higher levels in UCB, whereas genes related to cell adhesion, morphogenesis, secretion, angiogenesis and neurogenesis were more expressed in UC cells. These molecular differences verified in tissue-specific MSC gene expression may reflect functional activities influenced by distinct niches and should be considered when developing clinical protocols involving MSC from different sources. In addition, these findings reinforce our previous suggestion on the importance of banking the whole umbilical cord unit for research or future therapeutic use.