A Predictive Coexpression Network Identifies Novel Genes Controlling the Seed-to-Seedling Phase Transition in Arabidopsis thaliana

Developmental pleiotropy of SDP1 from seedling to mature stages in B. napus

Rapeseed, an important oil crop, relies on robust seedling emergence for optimal yields. Seedling emergence in the field is vulnerable to various factors, among which inadequate self-supply of energy is crucial to limiting seedling growth in early stage. SUGAR-DEPENDENT1 (SDP1) initiates triacylglycerol (TAG) degradation, yet its detailed function has not been determined in B. napus. Here, we focused on the effects of plant growth during whole growth stages and energy mobilization during seedling establishment by mutation in BnSDP1. Protein sequence alignment and haplotypic analysis revealed the conservation of SDP1 among species, with a favorable haplotype enhancing oil content. Investigation of agronomic traits indicated bnsdp1 had a minor impact on vegetative growth and no obvious developmental defects when compared with wild type (WT) across growth stages. The seed oil content was improved by 2.0-2.37% in bnsdp1 lines, with slight reductions in silique length and seed number per silique. Furthermore, bnsdp1 resulted in lower seedling emergence, characterized by a shrunken hypocotyl and poor photosynthetic capacity in the early stages. Additionally, impaired seedling growth, especially in yellow seedlings, was not fully rescued in medium supplemented with exogenous sucrose. The limited lipid turnover in bnsdp1 was accompanied by induced amino acid degradation and PPDK-dependent gluconeogenesis pathway. Analysis of the metabolites in cotyledons revealed active amino acid metabolism and suppressed lipid degradation, consistent with the RNA-seq results. Finally, we proposed strategies for applying BnSDP1 in molecular breeding. Our study provides theoretical guidance for understanding trade-off between oil accumulation and seedling energy mobilization in B. napus.

Expression of the Arabidopsis redox-related LEA protein, SAG21 is regulated by ERF, NAC and WRKY transcription factors

SAG21/LEA5 is an unusual late embryogenesis abundant protein in Arabidopsis thaliana, that is primarily mitochondrially located and may be important in regulating translation in both chloroplasts and mitochondria. SAG21 expression is regulated by a plethora of abiotic and biotic stresses and plant growth regulators indicating a complex regulatory network. To identify key transcription factors regulating SAG21 expression, yeast-1-hybrid screens were used to identify transcription factors that bind the 1685 bp upstream of the SAG21 translational start site. Thirty-three transcription factors from nine different families bound to the SAG21 promoter, including members of the ERF, WRKY and NAC families. Key binding sites for both NAC and WRKY transcription factors were tested through site directed mutagenesis indicating the presence of cryptic binding sites for both these transcription factor families. Co-expression in protoplasts confirmed the activation of SAG21 by WRKY63/ABO3, and SAG21 upregulation elicited by oligogalacturonide elicitors was partially dependent on WRKY63, indicating its role in SAG21 pathogen responses. SAG21 upregulation by ethylene was abolished in the erf1 mutant, while wound-induced SAG21 expression was abolished in anac71 mutants, indicating SAG21 expression can be regulated by several distinct transcription factors depending on the stress condition.

Interplay between coding and non-coding regulation drives the Arabidopsis seed-to-seedling transition

Translation of seed stored mRNAs is essential to trigger germination. However, when RNAPII re-engages RNA synthesis during the seed-to-seedling transition has remained in question. Combining csRNA-seq, ATAC-seq and smFISH in Arabidopsis thaliana we demonstrate that active transcription initiation is detectable during the entire germination process. Features of non-coding regulation such as dynamic changes in chromatin accessible regions, antisense transcription, as well as bidirectional non-coding promoters are widespread throughout the Arabidopsis genome. We show that sensitivity to exogenous ABSCISIC ACID (ABA) during germination depends on proximal promoter accessibility at ABA-responsive genes. Moreover, we provide genetic validation of the existence of divergent transcription in plants. Our results reveal that active enhancer elements are transcribed producing non-coding enhancer RNAs (eRNAs) as widely documented in metazoans. In sum, this study defining the extent and role of coding and non-coding transcription during key stages of germination expands our understanding of transcriptional mechanisms underlying plant developmental transitions.

SeedTransNet: a directional translational network revealing regulatory patterns during seed maturation and germination

Seed maturation is the developmental process that prepares the embryo for the desiccated waiting period before germination. It is associated with a series of physiological changes leading to the establishment of seed dormancy, seed longevity, and desiccation tolerance. We studied translational changes during seed maturation and observed a gradual reduction in global translation during seed maturation. Transcriptome and translatome profiling revealed specific reduction in the translation of thousands of genes. By including previously published data on germination and seedling establishment, a regulatory network based on polysome occupancy data was constructed: SeedTransNet. Network analysis predicted translational regulatory pathways involving hundreds of genes with distinct functions. The network identified specific transcript sequence features suggesting separate translational regulatory circuits. The network revealed several seed maturation-associated genes as central nodes, and this was confirmed by specific seed phenotypes of the respective mutants. One of the regulators identified, an AWPM19 family protein, PM19-Like1 (PM19L1), was shown to regulate seed dormancy and longevity. This putative RNA-binding protein also affects the translational regulation of its target mRNA, as identified by SeedTransNet. Our data show the usefulness of SeedTransNet in identifying regulatory pathways during seed phase transitions.

Brassinosteroid gene regulatory networks at cellular resolution in the Arabidopsis root

Brassinosteroids are plant steroid hormones that regulate diverse processes, such as cell division and cell elongation, through gene regulatory networks that vary in space and time. By using time series single-cell RNA sequencing to profile brassinosteroid-responsive gene expression specific to different cell types and developmental stages of the Arabidopsis root, we identified the elongating cortex as a site where brassinosteroids trigger a shift from proliferation to elongation associated with increased expression of cell wall-related genes. Our analysis revealed HOMEOBOX FROM ARABIDOPSIS THALIANA 7 (HAT7) and GT-2-LIKE 1 (GTL1) as brassinosteroid-responsive transcription factors that regulate cortex cell elongation. These results establish the cortex as a site of brassinosteroid-mediated growth and unveil a brassinosteroid signaling network regulating the transition from proliferation to elongation, which illuminates aspects of spatiotemporal hormone responses.

To curve for survival: Apical hook development

Apical hook is a simple curved structure formed at the upper part of hypocotyls when dicot seeds germinate in darkness. The hook structure is transient but essential for seedlings' survival during soil emergence due to its efficient protection of the delicate shoot apex from mechanical injury. As a superb model system for studying plant differential growth, apical hook has fascinated botanists as early as the Darwin age, and significant advances have been achieved at both the morphological and molecular levels to understand how apical hook development is regulated. Here, we will mainly summarize the research progress at these two levels. We will also briefly compare the growth dynamics between apical hook and hypocotyl gravitropic bending at early seed germination phase, with the aim to deduce a certain consensus on their connections. Finally, we will outline the remaining questions and future research perspectives for apical hook development.

Identification and functional validation of super-enhancers in Arabidopsis thaliana

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

AtHB40 modulates primary root length and gravitropism involving CYCLINB and auxin transporters

Gravitropism is a finely regulated tropistic response based on the plant perception of directional cues. Such perception allows them to direct shoot growth upwards, above ground, and root growth downwards, into the soil, anchoring the plant to acquire water and nutrients. Gravity sensing occurs in specialized cells and depends on auxin distribution, regulated by influx/efflux carriers. Here we report that AtHB40, encoding a transcription factor of the homeodomain-leucine zipper I family, was expressed in the columella and the root tip. Athb40 mutants exhibited longer primary roots. Enhanced primary root elongation was in agreement with a higher number of cells in the transition zone and the induction of CYCLINB transcript levels. Moreover, athb40 mutants and AtHB40 overexpressors displayed enhanced and delayed gravitropistic responses, respectively. These phenotypes were associated with altered auxin distribution and deregulated expression of the auxin transporters LAX2, LAX3, and PIN2. Accordingly, lax2 and lax3 mutants also showed an altered gravitropistic response, and LAX3 was identified as a direct target of AtHB40. Furthermore, AtHB40 is induced by AtHB53 when the latter is upregulated by auxin. Altogether, these results indicate that AtHB40 modulates cell division and auxin distribution in the root tip thus altering primary root length and gravitropism.

Seed-to-Seedling Transition in Pisum sativum L.: A Transcriptomic Approach

The seed-to-seedling transition is a crucial step in the plant life cycle. The transition occurs at the end of seed germination and corresponds to the initiation of embryonic root growth. To improve our understanding of how a seed transforms into a seedling, we germinated the Pisum sativum L. seeds for 72 h and divided them into samples before and after radicle protrusion. Before radicle protrusion, seeds survived after drying and formed normally developed seedlings upon rehydration. Radicle protrusion increased the moisture content level in seed axes, and the accumulation of ROS first generated in the embryonic root and plumule. The water and oxidative status shift correlated with the desiccation tolerance loss. Then, we compared RNA sequencing-based transcriptomics in the embryonic axes isolated from pea seeds before and after radicle protrusion. We identified 24,184 differentially expressed genes during the transition to the post-germination stage. Among them, 2101 genes showed more prominent expression. They were related to primary and secondary metabolism, photosynthesis, biosynthesis of cell wall components, redox status, and responses to biotic stress. On the other hand, 415 genes showed significantly decreased expression, including the groups related to water deprivation (eight genes) and response to the ABA stimulus (fifteen genes). We assume that the water deprivation group, especially three genes also belonging to ABA stimulus (LTI65, LTP4, and HVA22E), may be crucial for the desiccation tolerance loss during a metabolic switch from seed to seedling. The latter is also accompanied by the suppression of ABA-related transcription factors ABI3, ABI4, and ABI5. Among them, HVA22E, ABI4, and ABI5 were highly conservative in functional domains and showed homologous sequences in different drought-tolerant species. These findings elaborate on the critical biochemical pathways and genes regulating seed-to-seedling transition.

Morphological, transcriptomic and metabolomic analyses of Sophora davidii mutants for plant height

Sophora davidii is an important plant resource in the karst region of Southwest China, but S. davidii plant-height mutants are rarely reported. Therefore, we performed phenotypic, anatomic structural, transcriptomic and metabolomic analyses to study the mechanisms responsible for S. davidii plant-height mutants. Phenotypic and anatomical observations showed that compared to the wild type, the dwarf mutant displayed a significant decrease in plant height, while the tall mutant displayed a significant increase in plant height. The dwarf mutant cells were smaller and more densely arranged, while those of the wild type and the tall mutant were larger and loosely arranged. Transcriptomic analysis revealed that differentially expressed genes (DEGs) involved in cell wall biosynthesis, expansion, phytohormone biosynthesis, signal transduction pathways, flavonoid biosynthesis and phenylpropanoid biosynthesis were significantly enriched in the S. davidii plant-height mutants. Metabolomic analysis revealed 57 significantly differential metabolites screened from both the dwarf and tall mutants. A total of 8 significantly different flavonoid compounds were annotated to LIPID MAPS, and three metabolites (chlorogenic acid, kaempferol and scopoletin) were involved in phenylpropanoid biosynthesis and flavonoid biosynthesis. These results shed light on the molecular mechanisms of plant height in S. davidii mutants and provide insight for further molecular breeding programs.

The Arabidopsis transcription factors AtPHL1 and AtHB23 act together promoting carbohydrate transport from pedicel-silique nodes to seeds

Carbohydrates are produced in green tissues through photosynthesis and then transported to sink tissues. Carbon partitioning is a strategic process, fine regulated, involving specific sucrose transporters in each connecting tissue. Here we report that a screening of an Arabidopsis transcription factor (TF) library using the homeodomain-leucine zipper I member AtHB23 as bait, allowed identifying the TF AtPHL1 interacting with the former. An independent Y2H assay, and in planta by BiFC, confirmed such interaction. AtHB23 and AtPHL1 coexpressed in the pedicel-silique nodes and the funiculus. Mutant plants (phl1, and amiR23) showed a marked reduction of lipid content in seeds, although lipid composition did not change compared to the wild type. While protein and carbohydrate contents were not significantly different between mutants and control mature seeds, we observed a reduced carbohydrate content in mutant plants young siliques (7 days after pollination). Moreover, using a CFDA probe, we revealed an impaired transport to the seeds, and the gene encoding the carbohydrate transporters SWEET10 and SWEET11, usually expressed in connecting tissues, was repressed in the amiR23 and phl1 mutant plants. Altogether, the results indicated that AtHB23 and AtPHL1 act together, promoting sucrose transport, and the lack of any of them provoked a reduction in seeds lipid content.

The evolutionarily conserved kinase SnRK1 orchestrates resource mobilization during Arabidopsis seedling establishment

The onset of plant life is characterized by a major phase transition. During early heterotrophic seedling establishment, seed storage reserves fuel metabolic demands, allowing the plant to switch to autotrophic metabolism. Although metabolic pathways leading to storage compound mobilization are well-described, the regulatory circuits remain largely unresolved. Using an inducible knockdown approach of the evolutionarily conserved energy master regulator Snf1-RELATED-PROTEIN-KINASE1 (SnRK1), phenotypic studies reveal its crucial function in Arabidopsis thaliana seedling establishment. Importantly, glucose feeding largely restores growth defects of the kinase mutant, supporting its major impact in resource mobilization. Detailed metabolite studies reveal sucrose as a primary resource early in seedling establishment, in a SnRK1-independent manner. Later, SnRK1 orchestrates catabolism of triacylglycerols and amino acids. Concurrent transcriptomic studies highlight SnRK1 functions in controlling metabolic hubs fuelling gluconeogenesis, as exemplified by cytosolic PYRUVATE ORTHOPHOSPHATE DIKINASE (cyPPDK). Here, SnRK1 establishes its function via phosphorylation of the transcription factor BASIC LEUCINE ZIPPER63 (bZIP63), which directly targets and activates the cyPPDK promoter. Taken together, our results disclose developmental and catabolic functions of SnRK1 in seed storage mobilization and describe a prototypic gene regulatory mechanism. As seedling establishment is important for plant vigor and crop yield, our findings are of agronomical importance.

Multi-omics network-based functional annotation of unknown Arabidopsis genes

Unraveling gene function is pivotal to understanding the signaling cascades that control plant development and stress responses. As experimental profiling is costly and labor intensive, there is a clear need for high-confidence computational annotation. In contrast to detailed gene-specific functional information, transcriptomics data are widely available for both model and crop species. Here, we describe a novel automated function prediction method, which leverages complementary information from multiple expression datasets by analyzing study-specific gene co-expression networks. First, we benchmarked the prediction performance on recently characterized Arabidopsis thaliana genes, and showed that our method outperforms state-of-the-art expression-based approaches. Next, we predicted biological process annotations for known (n = 15 790) and unknown (n = 11 865) genes in A. thaliana and validated our predictions using experimental protein-DNA and protein-protein interaction data (covering >220 000 interactions in total), obtaining a set of high-confidence functional annotations. Our method assigned at least one validated annotation to 5054 (42.6%) unknown genes, and at least one novel validated function to 3408 (53.0%) genes with computational annotations only. These omics-supported functional annotations shed light on a variety of developmental processes and molecular responses, such as flower and root development, defense responses to fungi and bacteria, and phytohormone signaling, and help fill the information gap on biological process annotations in Arabidopsis. An in-depth analysis of two context-specific networks, modeling seed development and response to water deprivation, shows how previously uncharacterized genes function within the respective networks. Moreover, our automated function prediction approach can be applied in future studies to facilitate gene discovery for crop improvement.

Dynamic network biomarker analysis discovers IbNAC083 in the initiation and regulation of sweet potato root tuberization

The initiation and development of storage roots (SRs) are intricately regulated by a transcriptional regulatory network. One key challenge is to accurately pinpoint the tipping point during the transition from pre-swelling to SRs and to identify the core regulators governing such a critical transition. To solve this problem, we performed a dynamic network biomarker (DNB) analysis of transcriptomic dynamics during root development in Ipomoea batatas (sweet potato). First, our analysis identified stage-specific expression patterns for a significant proportion (>9%) of the sweet potato genes and unraveled the chronology of events that happen at the early and later stages of root development. Then, the results showed that different root developmental stages can be depicted by co-expressed modules of sweet potato genes. Moreover, we identified the key components and transcriptional regulatory network that determine root development. Furthermore, through DNB analysis an early stage, with a root diameter of 3.5 mm, was identified as the critical period of SR swelling initiation, which is consistent with morphological and metabolic changes. In particular, we identified a NAM/ATAF/CUC (NAC) domain transcription factor, IbNAC083, as a core regulator of this initiation in the DNB-associated network. Further analyses and experiments showed that IbNAC083, along with its associated differentially expressed genes, induced dysfunction of metabolism processes, including the biosynthesis of lignin, flavonol and starch, thus leading to the transition to swelling roots.

Genomic and transcriptomic inventory of membrane transporters in coffee: Exploring molecular mechanisms of metabolite accumulation

The genus Coffea (Rubiaceae) encompasses a group of perennial plant species, including a commodity crop from which seeds are roasted, ground, and infused to make one of the most appreciated beverages in the world. As an important tropical crop restricted to specific regions of the world, coffee production is highly susceptible to the effects of environmental instabilities (i.e., local year-to-year weather fluctuations and global climate change) and threatening pest pressures, not to mention an increasing quality rigor by consumers in industrialized countries. Specialized metabolites are substances that largely affect plant-environment interactions as well as how consumers experience agricultural products. Membrane transporters are key targets, albeit understudied, for understanding and tailoring the spatiotemporal distribution of specialized metabolites as they mediate and control molecular trafficking and substance accumulation. Therefore, we analyzed the transportome of C. canephora encoded within the 25,574 protein-coding genes annotated in the genome of this species and identified 1847 putative membrane transporters. Following, we mined 152 transcriptional profiles of C. canephora and C. arabica and performed a comprehensive co-expression analysis to identify transporters potentially involved in the accumulation of specialized metabolites associated with beverage quality and bioactivity attributes. In toto, this report points to an avenue of possibilities on Coffea genomic and transcriptomic data mining for genetic breeding strategies, which can lead to the development of new, resilient varieties for more sustainable coffee production systems.

An HD-ZIP transcription factor, MxHB13, integrates auxin-regulated and juvenility-determined control of adventitious rooting in Malus xiaojinensis

Adventitious root (AR) formation is a critical factor in the vegetative propagation of forestry and horticultural plants. Competence for AR formation declines in many species during the miR156/SPL-mediated vegetative phase change. Auxin also plays a regulatory role in AR formation. In apple rootstock, both high miR156 expression and exogenous auxin application are prerequisites for AR formation. However, the mechanism by which the miR156/SPL module interacts with auxin in controlling AR formation is unclear. In this paper, leafy cuttings of juvenile (Mx-J) and adult (Mx-A) phase Malus xiaojinensis were used in an RNA-sequencing experiment. The results revealed that numerous genes involved in phytohormone signaling, carbohydrate metabolism, cell dedifferentiation, and reactivation were downregulated in Mx-A cuttings in response to indole butyric acid treatment. Among the differentially expressed genes, an HD-ZIP transcription factor gene, MxHB13, was found to be under negative regulation of MdSPL26 by directly binding to MxHB13 promoter. MxTIFY9 interacts with MxSPL26 and may play a role in co-repressing the expression of MxHB13. The expression of MxTIFY9 was induced by exogenous indole butyric acid. MxHB13 binds to the promoter of MxABCB19-2 and positively affects the expression. A model is proposed in which MxHB13 links juvenility-limited and auxin-limited AR recalcitrance mechanisms in Mx-A.

Heat Maps: Perfect Maps for Quick Reading? Comparing Usability of Heat Maps with Different Levels of Generalization

Recently, due to Web 2.0 and neocartography, heat maps have become a popular map type for quick reading. Heat maps are graphical representations of geographic data density in the form of raster maps, elaborated by applying kernel density estimation with a given radius on point- or linear-input data. The aim of this study was to compare the usability of heat maps with different levels of generalization (defined by radii of 10, 20, 30, and 40 pixels) for basic map user tasks. A user study with 412 participants (16–20 years old, high school students) was carried out in order to compare heat maps that showed the same input data. The study was conducted in schools during geography or IT lessons. Objective (the correctness of the answer, response times) and subjective (response time self-assessment, task difficulty, preferences) metrics were measured. The results show that the smaller radius resulted in the higher correctness of the answers. A larger radius did not result in faster response times. The participants perceived the more generalized maps as easier to use, although this result did not match the performance metrics. Overall, we believe that heat maps, in given circumstances and appropriate design settings, can be considered an efficient method for spatial data presentation.

Characterization of the Heat-Stable Proteome during Seed Germination in Arabidopsis with Special Focus on LEA Proteins

During seed germination, desiccation tolerance is lost in the radicle with progressing radicle protrusion and seedling establishment. This process is accompanied by comprehensive changes in the metabolome and proteome. Germination of Arabidopsis seeds was investigated over 72 h with special focus on the heat-stable proteome including late embryogenesis abundant (LEA) proteins together with changes in primary metabolites. Six metabolites in dry seeds known to be important for seed longevity decreased during germination and seedling establishment, while all other metabolites increased simultaneously with activation of growth and development. Thermo-stable proteins were associated with a multitude of biological processes. In the heat-stable proteome, a relatively similar proportion of fully ordered and fully intrinsically disordered proteins (IDP) was discovered. Highly disordered proteins were found to be associated with functional categories development, protein, RNA and stress. As expected, the majority of LEA proteins decreased during germination and seedling establishment. However, four germination-specific dehydrins were identified, not present in dry seeds. A network analysis of proteins, metabolites and amino acids generated during the course of germination revealed a highly connected LEA protein network.

The underground life of homeodomain-leucine zipper transcription factors

Roots are the anchorage organs of plants, responsible for water and nutrient uptake, exhibiting high plasticity. Root architecture is driven by the interactions of biomolecules, including transcription factors and hormones that are crucial players regulating root plasticity. Multiple transcription factor families are involved in root development; some, such as ARFs and LBDs, have been well characterized, whereas others remain less well investigated. In this review, we synthesize the current knowledge about the involvement of the large family of homeodomain-leucine zipper (HD-Zip) transcription factors in root development. This family is divided into four subfamilies (I-IV), mainly according to structural features, such as additional motifs aside from HD-Zip, as well as their size, gene structure, and expression patterns. We explored and analyzed public databases and the scientific literature regarding HD-Zip transcription factors in Arabidopsis and other species. Most members of the four HD-Zip subfamilies are expressed in specific cell types and several individuals from each group have assigned functions in root development. Notably, a high proportion of the studied proteins are part of intricate regulation pathways involved in primary and lateral root growth and development.

To dry perchance to live: Insights from the genome of the desiccation-tolerant biocrust moss Syntrichia caninervis

With global climate change, water scarcity threatens whole agro/ecosystems. The desert moss Syntrichia caninervis, an extremophile, offers novel insights into surviving desiccation and heat. The sequenced S. caninervis genome consists of 13 chromosomes containing 16 545 protein-coding genes and 2666 unplaced scaffolds. Syntenic relationships within the S. caninervis and Physcomitrella patens genomes indicate the S. caninervis genome has undergone a single whole genome duplication event (compared to two for P. patens) and evidence suggests chromosomal or segmental losses in the evolutionary history of S. caninervis. The genome contains a large sex chromosome composed primarily of repetitive sequences with a large number of Copia and Gypsy elements. Orthogroup analyses revealed an expansion of ELIP genes encoding proteins important in photoprotection. The transcriptomic response to desiccation identified four structural clusters of novel genes. The genomic resources established for this extremophile offer new perspectives for understanding the evolution of desiccation tolerance in plants.

Temporal dynamics of QTL effects on vegetative growth in Arabidopsis thaliana

We assessed early vegetative growth in a population of 382 accessions of Arabidopsis thaliana using automated non-invasive high-throughput phenotyping. All accessions were imaged daily from 7 d to 18 d after sowing in three independent experiments and genotyped using the Affymetrix 250k SNP array. Projected leaf area (PLA) was derived from image analysis and used to calculate relative growth rates (RGRs). In addition, initial seed size was determined. The generated datasets were used jointly for a genome-wide association study that identified 238 marker-trait associations (MTAs) individually explaining up to 8% of the total phenotypic variation. Co-localization of MTAs occurred at 33 genomic positions. At 21 of these positions, sequential co-localization of MTAs for 2-9 consecutive days was observed. The detected MTAs for PLA and RGR could be grouped according to their temporal expression patterns, emphasizing that temporal variation of MTA action can be observed even during the vegetative growth phase, a period of continuous formation and enlargement of seemingly similar rosette leaves. This indicates that causal genes may be differentially expressed in successive periods. Analyses of the temporal dynamics of biological processes are needed to gain important insight into the molecular mechanisms of growth-controlling processes in plants.

Co-expression Networks From Gene Expression Variability Between Genetically Identical Seedlings Can Reveal Novel Regulatory Relationships

Co-expression networks are a powerful tool to understand gene regulation. They have been used to identify new regulation and function of genes involved in plant development and their response to the environment. Up to now, co-expression networks have been inferred using transcriptomes generated on plants experiencing genetic or environmental perturbation, or from expression time series. We propose a new approach by showing that co-expression networks can be constructed in the absence of genetic and environmental perturbation, for plants at the same developmental stage. For this, we used transcriptomes that were generated from genetically identical individual plants that were grown under the same conditions and for the same amount of time. Twelve time points were used to cover the 24-h light/dark cycle. We used variability in gene expression between individual plants of the same time point to infer a co-expression network. We show that this network is biologically relevant and use it to suggest new gene functions and to identify new targets for the transcriptional regulators GI, PIF4, and PRR5. Moreover, we find different co-regulation in this network based on changes in expression between individual plants, compared to the usual approach requiring environmental perturbation. Our work shows that gene co-expression networks can be identified using variability in gene expression between individual plants, without the need for genetic or environmental perturbations. It will allow further exploration of gene regulation in contexts with subtle differences between plants, which could be closer to what individual plants in a population might face in the wild.

Network Analysis Prioritizes DEWAX and ICE1 as the Candidate Genes for Major eQTL Hotspots in Seed Germination of Arabidopsis thaliana

Seed germination is characterized by a constant change of gene expression across different time points. These changes are related to specific processes, which eventually determine the onset of seed germination. To get a better understanding on the regulation of gene expression during seed germination, we performed a quantitative trait locus mapping of gene expression (eQTL) at four important seed germination stages (primary dormant, after-ripened, six-hour after imbibition, and radicle protrusion stage) using Arabidopsis thaliana Bay x Sha recombinant inbred lines (RILs). The mapping displayed the distinctness of the eQTL landscape for each stage. We found several eQTL hotspots across stages associated with the regulation of expression of a large number of genes. Interestingly, an eQTL hotspot on chromosome five collocates with hotspots for phenotypic and metabolic QTL in the same population. Finally, we constructed a gene co-expression network to prioritize the regulatory genes for two major eQTL hotspots. The network analysis prioritizes transcription factors DEWAX and ICE1 as the most likely regulatory genes for the hotspot. Together, we have revealed that the genetic regulation of gene expression is dynamic along the course of seed germination.

Prediction of condition-specific regulatory genes using machine learning

Recent advances in genomic technologies have generated data on large-scale protein–DNA interactions and open chromatin regions for many eukaryotic species. How to identify condition-specific functions of transcription factors using these data has become a major challenge in genomic research. To solve this problem, we have developed a method called ConSReg, which provides a novel approach to integrate regulatory genomic data into predictive machine learning models of key regulatory genes. Using Arabidopsis as a model system, we tested our approach to identify regulatory genes in data sets from single cell gene expression and from abiotic stress treatments. Our results showed that ConSReg accurately predicted transcription factors that regulate differentially expressed genes with an average auROC of 0.84, which is 23.5–25% better than enrichment-based approaches. To further validate the performance of ConSReg, we analyzed an independent data set related to plant nitrogen responses. ConSReg provided better rankings of the correct transcription factors in 61.7% of cases, which is three times better than other plant tools. We applied ConSReg to Arabidopsis single cell RNA-seq data, successfully identifying candidate regulatory genes that control cell wall formation. Our methods provide a new approach to define candidate regulatory genes using integrated genomic data in plants.

Desiccation Tolerance: Avoiding Cellular Damage During Drying and Rehydration

Desiccation of plants is often lethal but is tolerated by the majority of seeds and by vegetative tissues of only a small number of land plants. Desiccation tolerance is an ancient trait, lost from vegetative tissues following the appearance of tracheids but reappearing in several lineages when selection pressures favored its evolution. Cells of all desiccation-tolerant plants and seeds must possess a core set of mechanisms to protect them from desiccation- and rehydration-induced damage. This review explores how desiccation generates cell damage and how tolerant cells assuage the complex array of mechanical, structural, metabolic, and chemical stresses and survive.Likewise, the stress of rehydration requires appropriate mitigating cellular responses. We also explore what comparative genomics, both structural and responsive, have added to our understanding of cellular protection mechanisms induced by desiccation, and how vegetative desiccation tolerance circumvents destructive, stress-induced cell senescence.

Salicylic Acid Suppresses Apical Hook Formation via NPR1-Mediated Repression of EIN3 and EIL1 in Arabidopsis

Salicylic acid (SA) and ethylene (ET) are important phytohormones that regulate numerous plant growth, development, and stress response processes. Previous studies have suggested functional interplay of SA and ET in defense responses, but precisely how these two hormones coregulate plant growth and development processes remains unclear. Our present work reveals antagonism between SA and ET in apical hook formation, which ensures successful soil emergence of etiolated dicotyledonous seedlings. Exogenous SA inhibited ET-induced expression of HOOKLESS1 (HLS1) in Arabidopsis (Arabidopsis thaliana) in a manner dependent on ETHYLENE INSENSITIVE3 (EIN3) and EIN3-LIKE1 (EIL1), the core transcription factors in the ET signaling pathway. SA-activated NONEXPRESSER OF PR GENES1 (NPR1) physically interacted with EIN3 and interfered with the binding of EIN3 to target gene promoters, including the HLS1 promoter. Transcriptomic analysis revealed that NPR1 and EIN3/EIL1 coordinately regulated subsets of genes that mediate plant growth and stress responses, suggesting that the interaction between NPR1 and EIN3/EIL1 is an important mechanism for integrating the SA and ET signaling pathways in multiple physiological processes. Taken together, our findings illuminate the molecular mechanism underlying SA regulation of apical hook formation as well as the antagonism between SA and ET in early seedling establishment and possibly other physiological processes.

HD-ZIP I Transcription Factor (PtHB13) Negatively Regulates Citrus Flowering through Binding to FLOWERING LOCUS C Promoter

For floral induction in adult citrus, low temperature is one of the most important environmental factors. FLOWERING LOCUS C (FLC) plays a very important role in low-temperature-induced Arabidopsis flowering by repressed FLC expression under exposure to prolonged low-temperature conditions. However, little is known about the FLC regulation mechanism in perennial woody plants such as citrus. In this study, the functions of citrus FLC homolog (PtFLC) were investigated by ectopic expression in Arabidopsis. Transcription factor of homeodomain leucine zipper I (HD-ZIP I) as an upstream regulator of PtFLC was identified by yeast one-hybrid screen to regulate its transcription. The HD-ZIP I transcription factor was highly homologous to Arabidopsis ATHB13 and thus was named PtHB13. Ectopically expressed PtHB13 inhibited flowering in transgenic Arabidopsis. Furthermore, the expression of PtFLC and PtHB13 showed a seasonal change during the floral induction period and was also affected by low temperature. Thus, we propose that PtHB13 binds to PtFLC promoter to regulate its activity during the citrus floral induction process.

AtHB23 participates in the gene regulatory network controlling root branching, and reveals differences between secondary and tertiary roots

In Arabidopsis, lateral root (LR) development is mainly controlled by several known auxin-regulated transcription factors (TFs). Here, we show that AtHB23 (a homeodomain-leucine zipper I TF) participates in this intricate network. Our study of the expression pattern of AtHB23 revealed that it is transcriptionally activated in the early stages of secondary LR primordium (LRP). We found that AtHB23 directly limits the expression of LBD16, a key factor in LR initiation, and also directly induces the auxin transporter gene LAX3. We propose that this HD-Zip I mediates the regulation of LAX3 by ARF7/19. Furthermore, AtHB23 plays distinct roles during the formation of secondary and tertiary roots, exhibiting differential expression patterns. ATHB23 is expressed throughout the tertiary root primordium, whereas it is restricted to early stages in secondary primordia, likely later repressing LBD16 in tertiary LR development and further inhibiting root emergence. Our results suggest that different genetic programs govern the formation of LRP from the main or secondary roots, thereby shaping the global dynamic architecture of the root system.

A Co-Expression Network in Hexaploid Wheat Reveals Mostly Balanced Expression and Lack of Significant Gene Loss of Homeologous Meiotic Genes Upon Polyploidization

Polyploidization has played an important role in plant evolution. However, upon polyploidization, the process of meiosis must adapt to ensure the proper segregation of increased numbers of chromosomes to produce balanced gametes. It has been suggested that meiotic gene (MG) duplicates return to a single copy following whole genome duplication to stabilize the polyploid genome. Therefore, upon the polyploidization of wheat, a hexaploid species with three related (homeologous) genomes, the stabilization process may have involved rapid changes in content and expression of MGs on homeologous chromosomes (homeologs). To examine this hypothesis, sets of candidate MGs were identified in wheat using co-expression network analysis and orthology informed approaches. In total, 130 RNA-Seq samples from a range of tissues including wheat meiotic anthers were used to define co-expressed modules of genes. Three modules were significantly correlated with meiotic tissue samples but not with other tissue types. These modules were enriched for GO terms related to cell cycle, DNA replication, and chromatin modification and contained orthologs of known MGs. Overall, 74.4% of genes within these meiosis-related modules had three homeologous copies which was similar to other tissue-related modules. Amongst wheat MGs identified by orthology, rather than co-expression, the majority (93.7%) were either retained in hexaploid wheat at the same number of copies (78.4%) or increased in copy number (15.3%) compared to ancestral wheat species. Furthermore, genes within meiosis-related modules showed more balanced expression levels between homeologs than genes in non-meiosis-related modules. Taken together, our results do not support extensive gene loss nor changes in homeolog expression of MGs upon wheat polyploidization. The construction of the MG co-expression network allowed identification of hub genes and provided key targets for future studies.

Transcriptome Analyses of FY Mutants Reveal Its Role in mRNA Alternative Polyadenylation

A crucial step for mRNA polyadenylation is poly(A) signal recognition by trans-acting factors. The mammalian cleavage and polyadenylation specificity factor (CPSF) complex components CPSF30 and WD repeat-containing protein33 (WDR33) recognize the canonical AAUAAA for polyadenylation. In Arabidopsis (Arabidopsis thaliana), the flowering time regulator FY is the homolog of WDR33. However, its role in mRNA polyadenylation is poorly understood. Using poly(A) tag sequencing, we found that >50% of alternative polyadenylation (APA) events are altered in fy single mutants or double mutants with oxt6 (a null mutant of AtCPSF30), but mutation of the FY WD40-repeat has a stronger effect than deletion of the plant-unique Pro-Pro-Leu-Pro-Pro (PPLPP) domain. fy mutations disrupt AAUAAA or AAUAAA-like poly(A) signal recognition. Notably, A-rich signal usage is suppressed in the WD40-repeat mutation but promoted in PPLPP-domain deficiency. However, fy mutations do not aggravate the altered signal usage in oxt6 Furthermore, the WD40-repeat mutation shows a preference for 3' untranslated region shortening, but the PPLPP-domain deficiency shows a preference for lengthening. Interestingly, the WD40-repeat mutant exhibits shortened primary roots and late flowering with alteration of APA of related genes. Importantly, the long transcripts of two APA genes affected in fy are related to abiotic stress responses. These results reveal a conserved and specific role of FY in mRNA polyadenylation.

Co-expression networks for plant biology: why and how

Co-expression network analysis is one of the most powerful approaches for interpretation of large transcriptomic datasets. It enables characterization of modules of co-expressed genes that may share biological functional linkages. Such networks provide an initial way to explore functional associations from gene expression profiling and can be applied to various aspects of plant biology. This review presents the applications of co-expression network analysis in plant biology and addresses optimized strategies from the recent literature for performing co-expression analysis on plant biological systems. Additionally, we describe the combined interpretation of co-expression analysis with other genomic data to enhance the generation of biologically relevant information.

Physiological and transcriptome analyses of Opisthopappus taihangensis in response to drought stress

Background

Water scarcity is considered to be a severe environmental constraint to plant survival and productivity. Studies on drought-tolerant plants would definitely promote a better understanding of the regulatory mechanism lying behind the adaptive response of plants to drought. Opisthopappus taihangensis (ling) shih is a typical drought-tolerant perennial plant species endemically distributed across the Taihang Mountains in China, but the underlying mechanism for drought tolerance of this particular species remains elusive.

Results

To mimic natural drought stress, O. taihangensis plants were treated with two different concentrations (25% and 5%) of polyethylene glycol (PEG6000), which represent the H group (high salinity) and the L group (low salinity), respectively. The physiological characteristics of these two groups of plants, including relative water content maintenance (RWC), proline content and chlorophyll content were assessed and compared with plants in the control group (CK), which had normal irrigation. There was not a significant difference in RWC when comparing plants in the L group with the control group. Proline was accumulated to a higher level, and chlorophyll content was decreased slightly in plants under low drought stress. In plants from the H group, a lower RWC was observed. Proline was accumulated to an even higher level when compared with plants from the L group, and chlorophyll content was further reduced in plants under high drought stress. Transcriptomic analysis was carried out to look for genes that are differentially expressed (DEGs) in O. taihangensis plants coping adaptively with the two levels of drought stress. A total of 23,056 genes are differentially expressed between CK and L, among which 12,180 genes are up-regulated and 10,876 genes are down-regulated. Between H and L, 6182 genes are up-regulated and 1850 genes are down-regulated, which gives a total of 8032 genes. The highest number of genes, that are differentially expressed, was obtained when a comparison was made between CK and H. A total of 43,074 genes were found to be differentially expressed with 26,977 genes up-regulated and 16,097 genes down-regulated. Further analysis of these genes suggests that many of the up-regulated genes are enriched in pathways involved in amino acid metabolism. Besides, 39 transcription factors (TFs) were found to be continuously up-regulated with the increase of drought stress level.

Conclusion

Taken together, the results indicate that O. taihangensis plants are able to live adaptively under drought stress by responding physiologically and regulating the expression of a substantial number of drought-responsive genes and TFs to avoid adverse effects.

Electronic supplementary material

The online version of this article (10.1186/s13578-019-0318-7) contains supplementary material, which is available to authorized users.

Abscisic acid-dependent histone demethylation during postgermination growth arrest in Arabidopsis

After germination, seedlings undergo growth arrest in response to unfavourable conditions, a critical adaptation enabling plants to survive harsh environments. The plant hormone abscisic acid (ABA) plays a key role in this arrest. To arrest growth, ABA-dependent transcription factors change gene expression patterns in a flexible and reversible manner. Although the control of gene expression has important roles in growth arrest, the epigenetic mechanisms in the response to ABA are not fully understood. Here, we show that the histone demethylases JUMONJI-C domain-containing protein 30 (JMJ30) and JMJ32 control ABA-mediated growth arrest in Arabidopsis thaliana. During the postgermination stage (2-3 days after germination), the ABA-dependent transcription factor ABA-insensitive3 (ABI3) activates the expression of JMJ30 in response to ABA. JMJ30 then removes a repressive histone mark, H3 lysine 27 trimethylation (H3K27me3), from the SNF1-related protein kinase 2.8 (SnRK2.8) promoter, and hence activates SnRK2.8 expression. SnRK2.8 encodes a kinase that activates ABI3 and is responsible for JMJ30- and JMJ32-mediated growth arrest. A feed-forward loop involving the ABI3 transcription factor, JMJ histone demethylases, and the SnRK2.8 kinase fine-tunes ABA-dependent growth arrest in the postgermination phase. Our findings highlight the importance of the histone demethylases in mediating adaptation of plants to the environment.

Characterization of LRL5 as a key regulator of root hair growth in maize

Root hair, a special type of tubular-shaped cell, outgrows from the root epidermal cell and plays important roles in the acquisition of nutrients and water, as well as interactions with biotic and abiotic stresses. Studies in the model plant Arabidopsis have revealed that root-hair initiation and elongation are hierarchically regulated by a group of basic helix-loop-helix (bHLH) transcription factors (TFs). However, knowledge regarding the regulatory pathways of these bHLH TFs in controlling root hair growth remains limited. In this study, RNA-seq analysis was conducted to profile the transcriptome in the elongating maize root hair and >1000 genes with preferential expression in root hair were identified. A consensus cis-element previously featured as the potential bHLH-TF binding sites was present in the regulatory regions for the majority of the root hair-preferentially expressed genes. In addition, an individual change in ZmLRL5, the highest-expressed bHLH-TF in maize root hair resulted in a dramatic reduction in the elongation of root hair, and rendered the growth of root hair hypersensitive to translational inhibition. Moreover, RNA-seq, yeast-one-hybrid and ribosome profile analysis suggested that ZmLRL5 may function as a key player in orchestrating the translational process by directly regulating the expression of translational processes/ribosomal genes during maize root hair growth.

Genome-wide dynamic network analysis reveals a critical transition state of flower development in Arabidopsis

BACKGROUND

The flowering transition which is controlled by a complex and intricate gene regulatory network plays an important role in the reproduction for offspring of plants. It is a challenge to identify the critical transition state as well as the genes that control the transition of flower development. With the emergence of massively parallel sequencing, a great number of time-course transcriptome data greatly facilitate the exploration of the developmental phase transition in plants. Although some network-based bioinformatics analyses attempted to identify the genes that control the phase transition, they generally overlooked the dynamics of regulation and resulted in unreliable results. In addition, the results of these methods cannot be self-explained.

RESULTS

In this work, to reveal a critical transition state and identify the transition-specific genes of flower development, we implemented a genome-wide dynamic network analysis on temporal gene expression data in Arabidopsis by dynamic network biomarker (DNB) method. In the analysis, DNB model which can exploit collective fluctuations and correlations of different metabolites at a network level was used to detect the imminent critical transition state or the tipping point. The genes that control the phase transition can be identified by the difference of weighted correlations between the genes interested and the other genes in the global network. To construct the gene regulatory network controlling the flowering transition, we applied NARROMI algorithm which can reduce the noisy, redundant and indirect regulations on the expression data of the transition-specific genes. In the results, the critical transition state detected during the formation of flowers corresponded to the development of flowering on the 7th to 9th day in Arabidopsis. Among of 233 genes identified to be highly fluctuated at the transition state, a high percentage of genes with maximum expression in pollen was detected, and 24 genes were validated to participate in stress reaction process, as well as other floral-related pathways. Composed of three major subnetworks, a gene regulatory network with 150 nodes and 225 edges was found to be highly correlated with flowering transition. The gene ontology (GO) annotation of pathway enrichment analysis revealed that the identified genes are enriched in the catalytic activity, metabolic process and cellular process.

CONCLUSIONS

This study provides a novel insight to identify the real causality of the phase transition with genome-wide dynamic network analysis.

Comparative transcriptome study of hairy and hairless tea plant (Camellia sinensis) shoots

Trichome (also referred to as 'háo' in tea) is a key feature in both tea products and tea plant (Camellia sinensis) selection breeding. Although trichomes are used as a model for studying cell differentiation and have been well studied in many plant species, the regulation of trichome formation at the molecular level is poorly understood in tea plants. In the present study, the hairy and hairless tea plant cultivars Fudingdabaicha (FDDB) and Rongchunzao (RCZ), respectively, were used to study this mechanism. We characterised tea plant trichomes as unicellular and unbranched structures. High-throughput Illumina sequencing yielded approximately 277.0 million high-quality clean reads from the FDDB and RCZ cultivars. After de novo assembly, 161,444 unigenes were generated, with an average length of 937 bp. Among these unigenes, 81,425 were annotated using public databases, and 55,201 coding sequences and 4004 transcription factors (TFs) were identified. In total, 21,599 differentially expressed genes (DEGs) were identified between RCZ and FDDB, of which 10,785 DEGs were up-regulated and 10,814 DEGs were down-regulated. Genes involved in the DNA replication pathway were significantly enriched. Furthermore, between FDDB and RCZ, DEGs related to TFs, phytohormone signals, and cellulose synthesis were identified, suggesting that certain genes involved in these pathways are crucial for trichome initiation in tea plants. Together, the results of this study provide novel data to improve our understanding of the potential molecular mechanisms of trichome formation and lay a foundation for additional trichome studies in tea plants.

PlaD: A Transcriptomics Database for Plant Defense Responses to Pathogens, Providing New Insights into Plant Immune System

High-throughput transcriptomics technologies have been widely used to study plant transcriptional reprogramming during the process of plant defense responses, and a large quantity of gene expression data have been accumulated in public repositories. However, utilization of these data is often hampered by the lack of standard metadata annotation. In this study, we curated 2444 public pathogenesis-related gene expression samples from the model plant Arabidopsis and three major crops (maize, rice, and wheat). We organized the data into a user-friendly database termed as PlaD. Currently, PlaD contains three key features. First, it provides large-scale curated data related to plant defense responses, including gene expression and gene functional annotation data. Second, it provides the visualization of condition-specific expression profiles. Third, it allows users to search co-regulated genes under the infections of various pathogens. Using PlaD, we conducted a large-scale transcriptome analysis to explore the global landscape of gene expression in the curated data. We found that only a small fraction of genes were differentially expressed under multiple conditions, which might be explained by their tendency of having more network connections and shorter network distances in gene networks. Collectively, we hope that PlaD can serve as an important and comprehensive knowledgebase to the community of plant sciences, providing insightful clues to better understand the molecular mechanisms underlying plant immune responses. PlaD is freely available at http://systbio.cau.edu.cn/plad/index.php or http://zzdlab.com/plad/index.php.

CIPK9 is involved in seed oil regulation in Brassica napus L. and Arabidopsis thaliana (L.) Heynh

BACKGROUND

Accumulation of storage compounds during seed development plays an important role in the life cycle of oilseed plants; these compounds provide carbon and energy resources to support the establishment of seedlings.

RESULTS

In this study, we show that BnCIPK9 has a broad expression pattern in Brassica napus L. tissues and that wounding stress strongly induces its expression. The overexpression of BnCIPK9 during seed development reduced oil synthesis in transgenic B. napus compared to that observed in wild-type (WT) plants. Functional analysis revealed that seed oil content (OC) of complementation lines was similar to that of WT plants, whereas OC in Arabidopsis thaliana (L.) Heynh. Atcipk9 knockout mutants (cipk9) was higher than that of WT plants. Seedling of cipk9 mutants failed to establish roots on a sugar-free medium, but root establishment could be rescued by supplementation of sucrose or glucose. The phenotype of complementation transgenic lines was similar to that of WT plants when grown on sugar-free medium. Mutants, cipk9, cbl2, and cbl3 presented similar phenotypes, suggesting that CIPK9, CBL2, and CBL3 might work together and play similar roles in root establishment under sugar-free condition.

CONCLUSION

This study showed that BnCIPK9 and AtCIPK9 encode a protein kinase that is involved in sugar-related response and plays important roles in the regulation of energy reserves. Our results suggest that AtCIPK9 negatively regulates lipid accumulation and has a significant effect on early seedling establishment in A. thaliana. The functional characterization of CIPK9 provides insights into the regulation of OC, and might be used for improving OC in B. napus. We believe that our study makes a significant contribution to the literature because it provides information on how CIPKs coordinate stress regulation and energy signaling.

Seed Dormancy Involves a Transcriptional Program That Supports Early Plastid Functionality during Imbibition

Red rice fully dormant seeds do not germinate even under favorable germination conditions. In several species, including rice, seed dormancy can be removed by dry-afterripening (warm storage); thus, dormant and non-dormant seeds can be compared for the same genotype. A weedy (red) rice genotype with strong dormancy was used for mRNA expression profiling, by RNA-Seq, of dormant and non-dormant dehulled caryopses (here addressed as seeds) at two temperatures (30 °C and 10 °C) and two durations of incubation in water (8 h and 8 days). Aim of the study was to highlight the differences in the transcriptome of dormant and non-dormant imbibed seeds. Transcript data suggested important differences between these seeds (at least, as inferred by expression-based metabolism reconstruction): dry-afterripening seems to impose a respiratory impairment onto non-dormant seeds, thus glycolysis is deduced to be preferentially directed to alcoholic fermentation in non-dormant seeds but to alanine production in dormant ones; phosphoenolpyruvate carboxykinase, pyruvate phosphate dikinase and alanine aminotransferase pathways appear to have an important gluconeogenetic role associated with the restoration of plastid functions in the dormant seed following imbibition; correspondingly, co-expression analysis pointed out a commitment to guarantee plastid functionality in dormant seeds. At 8 h of imbibition, as inferred by gene expression, dormant seeds appear to preferentially use carbon and nitrogen resources for biosynthetic processes in the plastid, including starch and proanthocyanidins accumulation. Chromatin modification appears to be a possible mechanism involved in the transition from dormancy to germination. Non-dormant seeds show higher expression of genes related to cell wall modification, suggesting they prepare for acrospire/radicle elongation.

Melonet-DB, a Grand RNA-Seq Gene Expression Atlas in Melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important Cucurbitaceae crop produced worldwide, exhibiting wide genetic variations and comprising both climacteric and non-climacteric fruit types. The muskmelon cultivar "'Earl's favorite Harukei-3 (Harukei-3)"' known for its sweetness and rich aroma is used for breeding of high-grade muskmelon in Japan. We conducted RNA sequencing (RNA-seq) transcriptome studies in 30 different tissues of the 'Harukei-3' melon. These included root, stems, leaves, flowers, regenerating callus and ovaries, in addition to the flesh and peel sampled at seven stages of fruit development. The expression patterns of 20,752 genes were determined with fragments per kilobase of transcript per million fragments sequenced (FPKM) >1 in at least one tissue. Principal component analysis distinguished 30 melon tissues based on the global gene expression profile and, further, the weighted gene correlation network analysis classified melon genes into 45 distinct coexpression groups. Some coexpression groups exhibited tissue-specific gene expression. Furthermore, we developed and published web application tools designated "'Gene expression map viewer"' and "'Coexpression viewer"' on our website Melonet-DB (http://melonet-db.agbi.tsukuba.ac.jp/) to promote functional genomics research in melon. By using both tools, we analyzed melon homologs of tomato fruit ripening regulators such as E8, RIPENING-INHIBITOR (RIN) and NON-RIPENING (NOR). The "'Coexpression viewer"' clearly distinguished fruit ripening-associated melon RIN/NOR/CNR homologs from those expressed in other tissues. In addition, several other MADS-box, NAM/ATAF/CUC (NAC) and homeobox transcription factor genes were identified as fruit ripening-associated genes. Our tools provide useful information for research not only on melon but also on other fleshy fruit plants.

Metabolite profiling and associated gene expression reveal two metabolic shifts during the seed-to-seedling transition in Arabidopsis thaliana

Metabolic and transcriptomic correlation analysis identified two distinctive profiles involved in the metabolic preparation for seed germination and seedling establishment, respectively. Transcripts were identified that may control metabolic fluxes. The transition from a quiescent metabolic state (dry seed) to the active state of a vigorous seedling is crucial in the plant's life cycle. We analysed this complex physiological trait by measuring the changes in primary metabolism that occur during the transition in order to determine which metabolic networks are operational. The transition involves several developmental stages from seed germination to seedling establishment, i.e. between imbibition of the mature dry seed and opening of the cotyledons, the final stage of seedling establishment. We hypothesized that the advancement of growth is associated with certain signature metabolite profiles. Metabolite-metabolite correlation analysis underlined two specific profiles which appear to be involved in the metabolic preparation for seed germination and efficient seedling establishment, respectively. Metabolite profiles were also compared to transcript profiles and although transcriptional changes did not always equate to a proportional metabolic response, in depth correlation analysis identified several transcripts that may directly influence the flux through metabolic pathways during the seed-to-seedling transition. This correlation analysis also pinpointed metabolic pathways which are significant for the seed-to-seedling transition, and metabolite contents that appeared to be controlled directly by transcript abundance. This global view of the transcriptional and metabolic changes during the seed-to-seedling transition in Arabidopsis opens up new perspectives for understanding the complex regulatory mechanism underlying this transition.

Extensive transcriptomic and epigenomic remodelling occurs during Arabidopsis thaliana germination

BACKGROUND

Seed germination involves progression from complete metabolic dormancy to a highly active, growing seedling. Many factors regulate germination and these interact extensively, forming a complex network of inputs that control the seed-to-seedling transition. Our understanding of the direct regulation of gene expression and the dynamic changes in the epigenome and small RNAs during germination is limited. The interactions between genome, transcriptome and epigenome must be revealed in order to identify the regulatory mechanisms that control seed germination.

RESULTS

We present an integrated analysis of high-resolution RNA sequencing, small RNA sequencing and MethylC sequencing over ten developmental time points in Arabidopsis thaliana seeds, finding extensive transcriptomic and epigenomic transformations associated with seed germination. We identify previously unannotated loci from which messenger RNAs are expressed transiently during germination and find widespread alternative splicing and divergent isoform abundance of genes involved in RNA processing and splicing. We generate the first dynamic transcription factor network model of germination, identifying known and novel regulatory factors. Expression of both microRNA and short interfering RNA loci changes significantly during germination, particularly between the seed and the post-germinative seedling. These are associated with changes in gene expression and large-scale demethylation observed towards the end of germination, as the epigenome transitions from an embryo-like to a vegetative seedling state.

CONCLUSIONS

This study reveals the complex dynamics and interactions of the transcriptome and epigenome during seed germination, including the extensive remodelling of the seed DNA methylome from an embryo-like to vegetative-like state during the seed-to-seedling transition. Data are available for exploration in a user-friendly browser at https://jbrowse.latrobe.edu.au/germination_epigenome .

microRNAs participate in gene expression regulation and phytohormone cross-talk in barley embryo during seed development and germination

BACKGROUND

Small RNA and degradome sequencing have identified a large number of miRNA-target pairs in plant seeds. However, detailed spatial and temporal studies of miRNA-mediated regulation, which can reflect links between seed development and germination are still lacking.

RESULTS

In this study, we extended our investigation on miRNAs-involved gene regulation by a combined analysis of seed maturation and germination in barley. Through bioinformatics analysis of small RNA sequencing data, a total of 1324 known miRNA families and 448 novel miRNA candidates were identified. Of those, 16 known miRNAs with 40 target genes, and three novel miRNAs with four target genes were confirmed based on degradome sequencing data. Conserved miRNA families such as miR156, miR168, miR166, miR167, and miR894 were highly expressed in embryos of developing and germinating seeds. A barley-specific miRNA, miR5071, which was predicted to target an OsMLA10-like gene, accumulated at a high level, suggesting its involvement in defence response during these two developmental stages. Based on target prediction and Kyoto Encyclopedia of Genes and Genomes analysis of putative targets, nine highly expressed miRNAs were found to be related to phytohormone signalling and hormone cross-talk. Northern blot and qRT-PCR analysis showed that these miRNAs displayed differential expression patterns during seed development and germination, indicating their different roles in hormone signalling pathways. In addition, we showed that miR393 affected seed development through targeting two genes encoding the auxin receptors TIR1/AFBs in barley, as over-expression of miR393 led to an increased length-width ratio of seeds, whereas target mimic (MIM393)-mediated inhibition of its activity decreased the 1000-grain weight of seeds. Furthermore, the expression of auxin-responsive genes, abscisic acid- and gibberellic acid-related genes was altered in miR393 misexpression lines during germination and early seedling growth.

CONCLUSIONS

Our work indicates that miRNA-target pairs participate in gene expression regulation and hormone interaction in barley embryo and provides evidence that miR393-mediated auxin response regulation affects grain development and influences gibberellic acid and abscisic acid homeostasis during germination.

SeqEnrich: A tool to predict transcription factor networks from co-expressed Arabidopsis and Brassica napus gene sets

Transcription factors and their associated DNA binding sites are key regulatory elements of cellular differentiation, development, and environmental response. New tools that predict transcriptional regulation of biological processes are valuable to researchers studying both model and emerging-model plant systems. SeqEnrich predicts transcription factor networks from co-expressed Arabidopsis or Brassica napus gene sets. The networks produced by SeqEnrich are supported by existing literature and predicted transcription factor–DNA interactions that can be functionally validated at the laboratory bench. The program functions with gene sets of varying sizes and derived from diverse tissues and environmental treatments. SeqEnrich presents as a powerful predictive framework for the analysis of Arabidopsis and Brassica napus co-expression data, and is designed so that researchers at all levels can easily access and interpret predicted transcriptional circuits. The program outperformed its ancestral program ChipEnrich, and produced detailed transcription factor networks from Arabidopsis and Brassica napus gene expression data. The SeqEnrich program is ideal for generating new hypotheses and distilling biological information from large-scale expression data.

Extensive translational regulation during seed germination revealed by polysomal profiling

This work investigates the extent of translational regulation during seed germination. The polysome occupancy of each gene is determined by genome-wide profiling of total mRNA and polysome-associated mRNA. This reveals extensive translational regulation during Arabidopsis thaliana seed germination. The polysome occupancy of thousands of individual mRNAs changes to a large extent during the germination process. Intriguingly, these changes are restricted to two temporal phases (shifts) during germination, seed hydration and germination. Sequence features, such as upstream open reading frame number, transcript length, mRNA stability, secondary structures, and the presence and location of specific motifs correlated with this translational regulation. These features differed significantly between the two shifts, indicating that independent mechanisms regulate translation during seed germination. This study reveals substantial translational dynamics during seed germination and identifies development-dependent sequence features and cis elements that correlate with the translation control, uncovering a novel and important layer of gene regulation during seed germination.

Plant transcription factors from the homeodomain-leucine zipper family I. Role in development and stress responses

In front of stressful conditions plants display adaptation mechanisms leading to changes in their morphology, physiology, development and molecular composition. Transcription factors (TFs) play crucial roles in these complex adaptation processes. This work is focused in the homeodomain-leucine zipper I (HD-Zip I) family of TFs, unique to plants. First discovered in 1991, they were identified and isolated from monocotyledonous and dicotyledonous plants showing high structural similarity and diversified functions. These TFs have, besides the homeodomain and leucine zipper, conserved motifs in their carboxy-termini allowing the interaction with the basal machinery and with other regulatory proteins. The model dicotyledonous plant Arabidopsis thaliana has 17 HD-Zip I members; most of them regulated by external stimuli and hormones. These TFs are involved in key developmental processes like root and stem elongation, rosette leaves morphology determination, inflorescence stem branching, flowering and pollen hydration. Moreover, they are key players in responses to environmental stresses and illumination conditions. Several HD-Zip I encoding genes from different species were protected in patents because their overexpression or mutation generates improved agronomical phenotypes. Here we discuss many aspects about these TFs including structural features, biological functions and their utilization as biotechnological tools to improve crops. © 2017 IUBMB Life, 69(5):280-289, 2017.

Cross-kingdom comparison of the developmental hourglass

The developmental hourglass model has its foundations in classic anatomical studies by von Baer and Haeckel. In this context, even the conservation of animal body plans has been explained by evolutionary constraints acting on mid-embryogenic development. Recent studies have shown that developmental hourglass patterns also exist on the transcriptomic level, mirroring the corresponding morphological patterns. The identification of similar patterns in embryonic, post-embryonic, and life cycle spanning transcriptomes in plant and fungus development, however, contradict the notion of a direct coupling between morphological and molecular patterns. To explain the existence of hourglass patterns across kingdoms and developmental processes, we propose the organizational checkpoint model that integrates the developmental hourglass model into a framework of transcriptome switches.

JUNGBRUNNEN1 Confers Drought Tolerance Downstream of the HD-Zip I Transcription Factor AtHB13

Low water availability is the major environmental factor limiting growth and productivity of plants and crops and is therefore considered of high importance for agriculture affected by climate change. Identifying regulatory components controlling the response and tolerance to drought stress is thus of major importance. The NAC transcription factor (TF) JUNGBRUNNEN1 (JUB1) from Arabidopsis thaliana extends leaf longevity under non-stress growth conditions, lowers cellular hydrogen peroxide (H₂O₂) level, and enhances tolerance against heat stress and salinity. Here, we additionally find that JUB1 strongly increases tolerance to drought stress in Arabidopsis when expressed from both, a constitutive (CaMV 35S) and an abiotic stress-induced (RD29A) promoter. Employing a yeast one-hybrid screen we identified HD-Zip class I TF AtHB13 as an upstream regulator of JUB1. AtHB13 has previously been reported to act as a positive regulator of drought tolerance. AtHB13 and JUB1 thereby establish a joint drought stress control module.

A matter of quantity: Common features in the drought response of transgenic plants overexpressing HD-Zip I transcription factors

Plant responses to water deficit involve complex molecular mechanisms in which transcription factors have key roles. Previous reports ectopically overexpressed a few members of the homeodomain-leucine zipper I (HD-Zip I) family of transcription factors from different species, and the obtained transgenic plants exhibited drought tolerance which extent depended on the level of overexpression, triggering diverse molecular and physiological pathways. Here we show that most HD-Zip I genes are regulated by drought in the vegetative and/or reproductive stages. Moreover, uncharacterized members of this family were expressed as transgenes both in Col-0 and rdr6-12 backgrounds and were able to enhance drought tolerance in host plants. The extent of such tolerance depended on the expression level of the transgene and was significantly higher in transgenic rdr6-12 than in Col-0. Comparative transcriptome analyses of Arabidopsis thaliana plants overexpressing HD-Zip I proteins indicated that many members have common targets. Moreover, the water deficit tolerance exhibited by these plants is likely due to the induction and repression of certain of these common HD-Zip I-regulated genes. However, each HD-Zip I member regulates other pathways, which, in some cases, generate differential and potentially undesirable traits in addition to drought tolerance. In conclusion, only a few members of this family could become valuable tools to improve drought-tolerance.

Functional networks inference from rule-based machine learning models

BACKGROUND

Functional networks play an important role in the analysis of biological processes and systems. The inference of these networks from high-throughput (-omics) data is an area of intense research. So far, the similarity-based inference paradigm (e.g. gene co-expression) has been the most popular approach. It assumes a functional relationship between genes which are expressed at similar levels across different samples. An alternative to this paradigm is the inference of relationships from the structure of machine learning models. These models are able to capture complex relationships between variables, that often are different/complementary to the similarity-based methods.

RESULTS

We propose a protocol to infer functional networks from machine learning models, called FuNeL. It assumes, that genes used together within a rule-based machine learning model to classify the samples, might also be functionally related at a biological level. The protocol is first tested on synthetic datasets and then evaluated on a test suite of 8 real-world datasets related to human cancer. The networks inferred from the real-world data are compared against gene co-expression networks of equal size, generated with 3 different methods. The comparison is performed from two different points of view. We analyse the enriched biological terms in the set of network nodes and the relationships between known disease-associated genes in a context of the network topology. The comparison confirms both the biological relevance and the complementary character of the knowledge captured by the FuNeL networks in relation to similarity-based methods and demonstrates its potential to identify known disease associations as core elements of the network. Finally, using a prostate cancer dataset as a case study, we confirm that the biological knowledge captured by our method is relevant to the disease and consistent with the specialised literature and with an independent dataset not used in the inference process.

AVAILABILITY

The implementation of our network inference protocol is available at: http://ico2s.org/software/funel.html.