Metabolite analysis reveals flavonoids accumulation during flower development in Rhododendron pulchrum sweet (Ericaceae).

The azalea (Rhododendron simsii Planch.) is an important ornamental woody plant with various medicinal properties due to its phytochemical compositions and components. However little information on the metabolite variation during flower development in Rhododendron has been provided. In our study, a comparative analysis of the flavonoid profile was performed in Rhododendron pulchrum sweet at three stages of flower development, bud (stage 1), partially open flower (stage 2), and full bloom (stage 3). A total of 199 flavonoids, including flavone, flavonol, flavone C-glycosides, flavanone, anthocyanin, and isoflavone were identified. In hierarchical clustering analysis (HCA) and principal component analysis (PCA), the accumulation of flavonoids displayed a clear development stage variation. During flower development, 78 differential accumulated metabolites (DAMs) were identified, and most were enriched to higher levels at the full bloom stage. A total of 11 DAMs including flavone (chrysin, chrysoeriol O-glucuronic acid, and chrysoeriol O-hexosyl-O-pentoside), isoflavone (biochanin A), and flavonol (3,7-di-O-methyl quercetin and isorhamnetin) were significantly altered at three stages. In particular, 3,7-di-O-methyl quercetin was the top increased metabolite during flower development. Furthermore, integrative analyses of metabolomic and transcriptomic were conducted, revealing that the contents of isoflavone, biochanin A, glycitin, and prunetin were correlated with the expression of 2-hydroxyisoflavanone dehydratase (HIDH), which provide insight into the regulatory mechanism that controls isoflavone biosynthesis in R. pulchrum. This study will provide a new reference for increasing desired metabolites effectively by more accurate or appropriate genetic engineering strategies.

OsIPK2, a Rice Inositol Polyphosphate Kinase Gene, Is Involved in Phosphate Homeostasis and Root Development.

Phosphorus (P) is a growth-limiting nutrient for plants, which is taken up by root tissue from the environment as inorganic phosphate (Pi). To maintain an appropriate status of cellular Pi, plants have developed sophisticated strategies to sense the Pi level and modulate their root system architecture (RSA) under the ever-changing growth conditions. However, the molecular basis underlying the mechanism remains elusive. Inositol polyphosphate kinase (IPK2) is a key enzyme in the inositol phosphate metabolism pathway, which catalyzes the phosphorylation of IP3 into IP5 by consuming ATP. In this study, the functions of a rice inositol polyphosphate kinase gene (OsIPK2) in plant Pi homeostasis and thus physiological response to Pi signal were characterized. As a biosynthetic gene for phytic acid in rice, overexpression of OsIPK2 led to distinct changes in inositol polyphosphate profiles and an excessive accumulation of Pi levels in transgenic rice under Pi-sufficient conditions. The inhibitory effects of OsIPK2 on root growth were alleviated by Pi-deficient treatment compared with wild-type plants, suggesting the involvement of OsIPK2 in the Pi-regulated reconstruction of RSA. In OsIPK2-overexpressing plants, the altered acid phosphatase (APase) activities and misregulation of Pi-starvation-induced (PSI) genes were observed in roots under different Pi supply conditions. Notably, the expression of OsIPK2 also altered the Pi homeostasis and RSA in transgenic Arabidopsis. Taken together, our findings demonstrate that OsIPK2 plays an important role in Pi homeostasis and RSA adjustment in response to different environmental Pi levels in plants.

Comparative proteomics analysis of the responses to selenium in selenium-enriched alfalfa (Medicago sativa L.) leaves.

The mass of leaves and the chlorophyll and selenium content of alfalfa can be increased by the foliar spraying of selenite. To better understand the relationship between changes in the expression of specific proteins and the various metabolic and regulatory pathways affected by selenium treatment, labeling with Tandem Mass Tags (TMT) was used as a proteomics technique to compare control leaves with those enriched with Se. A total of 8,411 proteins were identified, the expression levels of 195 of which were significantly modified, 67 significantly up-regulated and 128 significantly down-regulated. Using gene functional classification and metabolic pathway annotation, selenium treatment was found to have a significant impact on metabolic processes. The energy and substances produced by the metabolic processes of a variety of carbohydrates, lipids, and amino acids, and the metabolism of carbon may be responsible for increasing the yield of alfalfa leaves. Administration of selenium substantially influenced Se-responsive proteins, including ABC transporter G family member 36, Probable glutathione S-transferase and cysteine tRNA ligase. Selenium treatment may also enhance photosynthesis and the defense response of cells. Furthermore, protein ubiquitination also played an important role in the selenium response of alfalfa leaves. In summary, a basic analysis of the selenium response pathway in alfalfa leaves at the proteomics level was conducted, which may assist in a more detailed elucidation of selenium enrichment in alfalfa in the future.

Metabolite profiling of violet, white and pink flowers revealing flavonoids composition patterns in Rhododendron pulchrum Sweet.

Flower color is the major characteristics and critical breeding program for most Rhododendron species. However, little is known about their coloration mechanism and color inheritance. In this study, petal pigment constituents of three Rhododendron pulchrum Sweet cultivars with different colors were clarified based on LCESI- MS/MS method. Using a broad-targeted metabolomic approach, a total of 149 flavonoids and their glycosylated or methylated derivatives were identified, including 18 anthocyanins (Pg, Cy, Dp, Pn, Pt, and Mv) and 32 flavonols (mainly kaempferol 3-O-glycosides and quercetin 3-O-glycosides). Moreover, anthocyanins were mainly represented by anthocyanidin-3-O-glycosides (glucoside, rutinoside, galactoside, and di-glycosides). Flavone and C-glycosylated flavone were major second metabolites responsible for the difference among three different R. pulchrum cultivars. The accumulation of total flavonoids displayed a clear phenotypic variation: cultivars 'zihe' and 'fenhe' were clustered together, while 'baihe' was clustered alone in the HCA analysis. The composition and content of anthocyanins were more complex in colored flowers ('zehe' and 'fenhe') than in white flower ('baihe'). This study further enhanced our understanding on the flavonoids profile of flower coloration and will provide biochemical basis for further genetic breeding in Rhododendron species.

Comprehensive Profiling and Inheritance Patterns of Metabolites in Foxtail Millet.

Metabolomics aims at determining a sample's metabolites profile and hence provides a straight functional statement of an organism's physiological condition. Here, we investigated comprehensive profiling, natural variation and species-specific accumulation of both primary and secondary metabolites in foxtail millet using LC-MS, and inheritance patterns of metabolome in millet hybrids. The application of a broad target metabolomics method facilitated the simultaneous identification and quantification of more than 300 metabolites. The metabolic analysis of these compounds, such as flavonoids, phenolamides, hydrocinnamoyl derivatives, vitamins and LPCs, revealed their developmentally controlled accumulation, and natural variation in different tissues/varieties. Species-specific accumulation of secondary metabolites was observed based on a comparative metabolic analysis between millet and rice, such as flavonoid O-rutinosides/neohesperidosides and malonylated flavonoid O-glycosides. In analyzing the metabolic variations between hybrid progenies and their parental lines, including a photothermo-sensitive genic male sterility line and five Zhangzagu varieties, metabolic overdominant, and dominant patterns of inheritance could be observed. For example, hydrocinnamoyl derivatives and feruloylated flavonoids were identified as over-parent heterosis (overdominant) metabolites in milet hybrids. Our work paves the way for developing predictors of hybrid performance and the future analysis of the biosynthesis and regulation of relevant metabolic pathways in millet.

Gene coexpression network analysis combined with metabonomics reveals the resistance responses to powdery mildew in Tibetan hulless barley.

Powdery mildew is a fungal disease that represents a ubiquitous threat to crop plants. Transcriptomic and metabolomic analyses were used to identify molecular and physiological changes in Tibetan hulless barley in response to powdery mildew. There were 3418 genes and 405 metabolites differentially expressed between the complete resistance cultivar G7 and the sensitive cultivar Z13. Weighted gene coexpression network analysis was carried out, and the differentially expressed genes were enriched in five and four major network modules in G7 and Z13, respectively. Further analyses showed that phytohormones, photosynthesis, phenylpropanoid biosynthesis, and flavonoid biosynthesis pathways were altered during Qingke-Blumeria graminis (DC.) f.sp. hordei (Bgh) interaction. Comparative analyses showed a correspondence between gene expression and metabolite profiles, and the activated defenses resulted in changes of metabolites involved in plant defense response, such as phytohormones, lipids, flavone and flavonoids, phenolamides, and phenylpropanoids. This study enabled the identification of Bgh responsive genes and provided new insights into the dynamic physiological changes that occur in Qingke during response to powdery mildew. These findings greatly improve our understanding of the mechanisms of induced defense response in Qingke and will provide new clues for the development of resistant Tibetan hulless barley varieties.

Transcriptome analysis and identification of genes associated with flower development in Rhododendron pulchrum Sweet (Ericaceae).

Flowering process is essential for plant development. However, the molecular mechanisms driving flower development of ornamental woody Rhododendron pulchrum Sweet are difficult to elucidate due to the lack of genomic data. In this research, high-throughput sequencing and comparative transcriptome analyses of R. pulchrum flowers collected at three key stages were performed: floral bud stage, early flowering stage, and full-flowering stage. Furthermore, expression of genes involved in flower development was also validated with quantitative real-time PCR (qRT-PCR). RNA-seq yielded 96,350,697 bp of clean reads, which were assembled into 98,610 unigenes with an average length of 717 bp. 58,279 (59.10%) unigenes could be annotated, including 324 major unigenes associated with floral development. In addition, ten modules (20,443 mRNAs) were dissected in the co-expression network. Especially, Flowering Locus (FLC) and Flowering Locus T (FT) were co-expressed. 9493 differentially expressed genes (DEGs) were scanned among three stages, and most DEGs existed between flower bud stage and early flowering stage. In particular, 79 DGEs associated with flowering process were enriched in 28 GO terms. Moreover, the expression levels of MYC2, EIN3, and ARR-B were all lowest at early flowering stage, while transcripts of MYC2, TIR1, CYCD3, COL-1, and EIN3 were all peaked at flower bud stage. Transcriptome profile presented here will benefit deep insights into molecular mechanism underlying R. pulchrum flowering process.

The Role of Arabidopsis Inositol Polyphosphate Kinase AtIPK2ß in Glucose Suppression of Seed Germination and Seedling Development.

Seed germination and subsequent seedling development are critical phases in plants. These processes are regulated by a complex molecular network in which sugar has been reported to play an essential role. However, factors affecting sugar responses remain to be fully elucidated. In this study, we demonstrate that AtIPK2ß, known to participate in the synthesis of myo-inositol 1,2,3,4,5,6-hexakisphosphate (IP6, phytate), affects Arabidopsis responses to glucose during seed germination. The loss-of-function mutant atipk2ß showed increased sensitivity to 6% glucose and paclobutrazol (PAC). Yeast two-hybrid assay showed that AtIPK2ß interacts with sucrose non-fermenting-1-related protein kinase (SnRK1.1), and bimolecular fluorescence complementation (BiFC) and pull-down assay further confirmed this interaction. Moreover, AtIPK2ß was phosphorylated by SnRK1.1 in vitro, and the effect of restoring AtIPK2ß to yeast cells lacking IPK2 (Δipk2) was abolished by catalytically active SnRK1.1. Further analysis indicated that IP6 reduces the suppression of seed germination caused by glucose, accompanied by altered expression levels of glucose-/hormone-responsive genes. Collectively, these findings indicate that AtIPK2ß and IP6 are involved in glucose suppression of seed germination and that AtIPK2ß enzyme activity is likely to be regulated by SnRK1.1.

Arabidopsis inositol polyphosphate multikinase delays flowering time through mediating transcriptional activation of FLOWERING LOCUS C.

Timely flowering is critical for successful reproduction and seed yield in plants. A diverse range of regulators have been found to control flowering time in response to environmental and endogenous signals. Among these regulators, FLOWERING LOCUS C (FLC) acts as a central repressor of floral transition by blocking the expression of flowering integrator genes. Here, we report that Arabidopsis inositol polyphosphate multikinase (AtIPK2ß) functions in flowering time control by mediating transcriptional regulation of FLC at the chromatin level. The atipk2ß mutant flowers earlier, and AtIPK2ß overexpressing plants exhibit late-flowering phenotypes. Quantitative reverse transcription-PCR (qRT-PCR) revealed that AtIPK2ß promotes FLC expression. We performed chromatin immunoprecipitation-qPCR (ChIP-qPCR) assays and found that AtIPK2ß binds to FLC chromatin. Further analysis showed that AtIPK2ß interacts with FVE, a key repressor required for epigenetic silencing of FLC. qRT-PCR, ChIP-qPCR, and genetic analysis demonstrated that AtIPK2ß is involved in FVE-mediated transcriptional regulation of FLC by repressing the accumulation of FVE on FLC. Moreover, we found that AtIPK2ß associates with HDA6, an interaction partner of FVE mediating FLC chromatin silencing, and attenuates HDA6 accumulation at the FLC locus. Taken together, these findings suggest that AtIPK2ß negatively regulates flowering time by blocking chromatin silencing of FLC.

Rice Inositol Polyphosphate Kinase (OsIPK2) Directly Interacts with OsIAA11 to Regulate Lateral Root Formation.

The plant hormone auxin controls many aspects of plant growth and development by promoting the degradation of Auxin/Indole-3-acetic acid (Aux/IAA) proteins. The domain II (DII) of Aux/IAA proteins is sufficient for eliciting the degradation by directly interacting with the auxin receptor F-box protein TIR1 to form a TIR1/AFBs-Aux/IAA complex in an auxin-dependent manner. However, the underlying mechanisms of fine-tuning Aux/IAA degradation by auxin stimuli remain to be elucidated. Here, we show that OsIPK2, a rice (Oryza sativa) inositol polyphosphate kinase, directly interacts with an Aux/IAA protein OsIAA11 to repress its degradation. In a rice protoplast transient expression system, the auxin-induced degradation of Myc-OsIAA11 fusion was delayed by co-expressed GFP-OsIPK2 proteins. Furthermore, expressing additional OsIPK2 or its N-terminal amino acid sequence enhanced the accumulation of OsIAA11 proteins in transgenic plants, which in turn caused defects in lateral root formation and auxin response. Taken together, we identify a novel co-factor of Aux/IAA in auxin signaling and demonstrate its role in regulating lateral root development.

Arabidopsis inositol polyphosphate kinase AtIpk2ß is phosphorylated by CPK4 and positively modulates ABA signaling.

Arabidopsis inositol polyphosphate kinase 2ß (AtIpk2ß) has multiple functions in plant development and in responding to abiotic stress. Although some related clues suggested a potential role of AtIpk2ß in ABA signaling, the defined evidence was still lack. Here we discovered that a key ABA signaling component calcium-dependent protein kinase 4 (CPK4) can interact with AtIpk2ß under ABA treated conditions through affinity purification and mass spectrometry detection. The interaction between CPK4 and AtIpk2ß were further confirmed by yeast two hybrid and bimolecular fluorescence complementation assays. Expression of AtIpk2ß also can be rapidly induced by ABA. In addition, we found that CPK4 can phosphorylate AtIpk2ß in vitro and identified five novel phosphorylation sites of AtIpk2ß by CPK4 kinase, including Tyr46, Ser48, Ser51, Thr128, Ser147. Overexpression of AtIpk2ß in Arabidopsis was more sensitive to ABA in seed germination, primary root inhibition, ABA-responsive gene expression than wild type plants, whereas knockout mutant atipk2ß exhibited no significant difference. The AtIpk2ß variants containing Tyr46, Thr128, Ser147 mutated to Ala cannot complement the yeast mutant ipk2 growth in high temperature, suggesting that those three amino acid residues are critical for AtIpk2ß. These findings provide insight into the modulation of ABA signaling by AtIpk2ß.

Light-Switchable Self-Healing Hydrogel Based on Host-Guest Macro-Crosslinking.

A self-healing hydrogel is prepared by crosslinking acrylamide with a host-guest macro-crosslinker assembled from poly(ß-cyclodextrin) nanogel and azobenzeneacrylamide. The photoisomerizable azobenzene moiety can change its binding affinity with ß-cyclodextrin, therefore the crosslinking density and rheology property of the hydrogel can be tuned with light stimulus. The hydrogel can repair its wound autonomously through the dynamic host-guest interaction. In addition, the wounded hydrogel will lose its ability of self-healing when exposed to ultraviolet light, and the self-healing behavior can be recovered upon the irradiation of visible light. The utilizing of host-guest macro-crosslinking approach manifests the as-prepared hydrogel reversible and light-switchable self-healing property, which would broaden the potential applications of self-healing polymers.

The low temperature-responsive, Solanum CBF1 genes maintain high identity in their upstream regions in a genomic environment undergoing gene duplications, deletions, and rearrangements.

Some plants like Arabidopsis thaliana increase in freezing tolerance when exposed to low nonfreezing temperatures, a process known as cold acclimation. Other plants including tomato, Solanum lycopersicum, are chilling sensitive and incur injury during prolonged low temperature exposure. A key initial event that occurs upon low temperature exposure is the induction of genes encoding the CBF transcription factors. In Arabidopsis three CBF genes, present in a tandemly-linked cluster, are induced by low temperatures. Tomato also harbors three tandemly-linked CBF genes, Sl-CBF3-CBF1-CBF2, but only one of these, Sl-CBF1, is low-temperature responsive. Here we report that Solanum species that are closely-allied to cultivated tomato essentially share this structural organization, but the locus is in a dynamic state of flux. Additional paralogs and in-frame deletions between adjacent genes occur, and the genomic regions flanking the CBF genes are dissimilar across Solanum species. Nevertheless, the CBF1 upstream region remains intact and highly conserved. This feature differed for CBF2 and CBF3, whose upstream regions were far less conserved. CBF1 was also the only low-temperature responsive gene in the cluster and its expression was greatly affected by a circadian clock. The tuber-bearing S. tuberosum and S. commersonii also harbored a fourth gene, CBF4, which was also low temperature responsive. CBF4 was physically linked to CBF5 in S. tuberosum, but CBF5 was absent from S. commersonii. Phylogenic analyses suggest that CBF5-CBF4 resulted from the duplication of the CBF3-CBF1-CBF2 cluster. DNA sequence motifs shared between the Solanum CBF1 and CBF4 upstream regions were identified, portions of which were also present in the Arabidopsis CBF1-3 upstream regions. These results suggest that much greater functional constraints are placed upon the Solanum CBF1 upstream regions over the other CBF upstream regions and that CBF4 has retained the capacity for low temperature responsiveness following the duplication event that gave rise to CBF4.

Consistent dissection of the protein interaction network by combining global and local metrics.

We propose a new network decomposition method to systematically identify protein interaction modules in the protein interaction network. Our method incorporates both a global metric and a local metric for balance and consistency. We have compared the performance of our method with several earlier approaches on both simulated and real datasets using different criteria, and show that our method is more robust to network alterations and more effective at discovering functional protein modules.

A parallel edge-betweenness clustering tool for Protein-Protein Interaction networks.

The increasing availability of protein-protein interaction graphs (PPI) requires new efficient tools capable of extracting valuable biological knowledge from these networks. Among the wide range of clustering algorithms, Girvan and Newman's edge betweenness algorithm showed remarkable performances in discovering clustering structures in several real-world networks. Unfortunately, their algorithm suffers from high computational cost and it is impractical for inputs of the size of large PPI networks. Here we report on a novel parallel implementation of Girvan and Newman's clustering algorithm that achieves almost linear speed-up for up to 32 processors. The tool is available in the public domain from the authors' website.

Pervasive adaptive evolution in mammalian fertilization proteins.

Mammalian fertilization exhibits species specificity, and the proteins mediating sperm-egg interactions evolve rapidly between species. In this study, we demonstrate that the evolution of seven genes involved in mammalian fertilization is promoted by positive Darwinian selection by using likelihood ratio tests (LRTs). Several of these proteins are sperm proteins that have been implicated in binding the mammalian egg coat zona pellucida glycoproteins, which were shown previously to be subjected to positive selection. Taken together, these represent the major candidates involved in mammalian fertilization, indicating positive selection is pervasive amongst mammalian reproductive proteins. A new LRT is implemented to determine if the d(N)/d(S) ratio is significantly greater than one. This is a more refined test of positive selection than the previous LRTs which only identified if there was a class of sites with a d(N)/d(S) ratio >1 but did not test if that ratio was significantly greater than one.