Spatially resolved transcriptomic analysis of the germinating barley grain.

Seeds are a vital source of calories for humans and a unique stage in the life cycle of flowering plants. During seed germination, the embryo undergoes major developmental transitions to become a seedling. Studying gene expression in individual seed cell types has been challenging due to the lack of spatial information or low throughput of existing methods. To overcome these limitations, a spatial transcriptomics workflow was developed for germinating barley grain. This approach enabled high-throughput analysis of spatial gene expression, revealing specific spatial expression patterns of various functional gene categories at a sub-tissue level. This study revealed over 14 000 genes differentially regulated during the first 24 h after imbibition. Individual genes, such as the aquaporin gene family, starch degradation, cell wall modification, transport processes, ribosomal proteins and transcription factors, were found to have specific spatial expression patterns over time. Using spatial autocorrelation algorithms, we identified auxin transport genes that had increasingly focused expression within subdomains of the embryo over time, suggesting their role in establishing the embryo axis. Overall, our study provides an unprecedented spatially resolved cellular map for barley germination and identifies specific functional genomics targets to better understand cellular restricted processes during germination. The data can be viewed at https://spatial.latrobe.edu.au/.

Coordinated regulation of the mitochondrial retrograde response by circadian clock regulators and ANAC017.

Mitochondrial retrograde signaling (MRS) supports photosynthetic function under a variety of conditions. Induction of mitochondrial dysfunction with myxothiazol (a specific inhibitor of the mitochondrial bc1 complex) or antimycin A (an inhibitor of the mitochondrial bc1 complex and cyclic electron transport in the chloroplast under light conditions) in the light and dark revealed diurnal control of MRS. This was evidenced by (1) significantly enhanced binding of ANAC017 to promoters in the light compared with the dark in Arabidopsis plants treated with myxothiazol (but not antimycin A), (2) overlap in the experimentally determined binding sites for ANAC017 and circadian clock regulators in the promoters of ANAC013 and AOX1a, (3) a diurnal expression pattern for ANAC017 and transcription factors it regulates, (4) altered expression of ANAC017-regulated genes in circadian clock mutants with and without myxothiazol treatment, and (5) a decrease in the magnitude of LHY and CCA1 expression in an ANAC017-overexpressing line and protein-protein interaction between ANAC017 and PIF4. This study also shows a large difference in transcriptome responses to antimycin A and myxothiazol in the dark: these responses are ANAC017 independent, observed in shoots and roots, similar to biotic challenge and salicylic acid responses, and involve ERF and ZAT transcription factors. This suggests that antimycin A treatment stimulates a second MRS pathway that is mediated or converges with salicylic acid signaling and provides a merging point with chloroplast retrograde signaling.

The retrograde signaling regulator ANAC017 recruits the MKK9-MPK3/6, ethylene, and auxin signaling pathways to balance mitochondrial dysfunction with growth.

In plant cells, mitochondria are ideally positioned to sense and balance changes in energy metabolism in response to changing environmental conditions. Retrograde signaling from mitochondria to the nucleus is crucial for adjusting the required transcriptional responses. We show that ANAC017, the master regulator of mitochondrial stress, directly recruits a signaling cascade involving the plant hormones ethylene and auxin as well as the MAP KINASE KINASE (MKK) 9-MAP KINASE (MPK) 3/6 pathway in Arabidopsis thaliana. Chromatin immunoprecipitation followed by sequencing and overexpression demonstrated that ANAC017 directly regulates several genes of the ethylene and auxin pathways, including MKK9, 1-AMINO-CYCLOPROPANE-1-CARBOXYLATE SYNTHASE 2, and YUCCA 5, in addition to genes encoding transcription factors regulating plant growth and stress responses such as BASIC REGION/LEUCINE ZIPPER MOTIF (bZIP) 60, bZIP53, ANAC081/ATAF2, and RADICAL-INDUCED CELL DEATH1. A time-resolved RNA-seq experiment established that ethylene signaling precedes the stimulation of auxin signaling in the mitochondrial stress response, with a large part of the transcriptional regulation dependent on ETHYLENE-INSENSITIVE 3. These results were confirmed by mutant analyses. Our findings identify the molecular components controlled by ANAC017, which integrates the primary stress responses to mitochondrial dysfunction with whole plant growth via the activation of regulatory and partly antagonistic feedback loops.

RNA-seq analysis of laser microdissected Arabidopsis thaliana leaf epidermis, mesophyll and vasculature defines tissue-specific transcriptional responses to multiple stress treatments.

Acclimation of plants to adverse conditions requires the coordination of gene expression and signalling pathways between tissues and cell types. As the energy and carbon capturing organs, leaves are significantly affected by abiotic and biotic stresses. However, tissue- or cell type-specific analyses of stress responses have focussed on the Arabidopsis root. Here, we comparatively explore the transcriptomes of three leaf tissues (epidermis, mesophyll, vasculature) after induction of diverse stress pathways by chemical stimuli (antimycin A, 3-amino-1,2,4-triazole, methyl viologen, salicylic acid) and ultraviolet light in Arabidopsis using laser capture microdissection followed by RNA sequencing. Stimulation of stress pathways caused an overall reduction in the number of genes expressed in a tissue-specific manner, though a small subset gained or changed their tissue specificity. We find no evidence of a common stress response, with only a few genes consistently responsive to two or more treatments in the analysed tissues. However, differentially expressed genes overlap between tissues for individual treatments. A focussed analysis provided evidence for an interaction of auxin and ethylene that mediates retrograde signalling during mitochondrial dysfunction specifically in the epidermis, and a gene regulatory network defined the hierarchy of interactions. Taken together, we have generated an extensive reference dataset that will be valuable for future experiments analysing transcriptional responses on a tissue or single-cell level. Our results will enable the tailoring of the tissue-specific engineering of stress-tolerant plants.

Laser-Capture Microdissection RNA-Sequencing for Spatial and Temporal Tissue-Specific Gene Expression Analysis in Plants.

The development of a complex multicellular organism is governed by distinct cell types that have different transcriptional profiles. To identify transcriptional regulatory networks that govern developmental processes it is necessary to measure the spatial and temporal gene expression profiles of these individual cell types. Therefore, insight into the spatio-temporal control of gene expression is essential to gain understanding of how biological and developmental processes are regulated. Here, we describe a laser-capture microdissection (LCM) method to isolate small number of cells from three barley embryo organs over a time-course during germination followed by transcript profiling. The method consists of tissue fixation, tissue processing, paraffin embedding, sectioning, LCM and RNA extraction followed by real-time PCR or RNA-seq. This method has enabled us to obtain spatial and temporal profiles of seed organ transcriptomes from varying numbers of cells (tens to hundreds), providing much greater tissue-specificity than typical bulk-tissue analyses. From these data we were able to define and compare transcriptional regulatory networks as well as predict candidate regulatory transcription factors for individual tissues. The method should be applicable to other plant tissues with minimal optimization.

Temporal tissue-specific regulation of transcriptomes during barley (Hordeum vulgare) seed germination.

The distinct functions of individual cell types require cells to express specific sets of genes. The germinating seed is an excellent model to study genome regulation between cell types since the majority of the transcriptome is differentially expressed in a short period, beginning from a uniform, metabolically inactive state. In this study, we applied laser-capture microdissection RNA-sequencing to small numbers of cells from the plumule, radicle tip and scutellum of germinating barley seeds every 8 h, over a 48 h time course. Tissue-specific gene expression was notably common; 25% (910) of differentially expressed transcripts in plumule, 34% (1876) in radicle tip and 41% (2562) in scutellum were exclusive to that organ. We also determined that tissue-specific storage of transcripts occurs during seed development and maturation. Co-expression of genes had strong spatiotemporal structure, with most co-expression occurring within one organ and at a subset of specific time points during germination. Overlapping and distinct enrichment of functional categories were observed in the tissue-specific profiles. We identified candidate transcription factors amongst these that may be regulators of spatiotemporal gene expression programs. Our findings contribute to the broader goal of generating an integrative model that describes the structure and function of individual cells within seeds during germination.

A novel role of the soybean clock gene LUX ARRHYTHMO in male reproductive development.

The evening complex of ELF4-ELF3-LUX proteins is an integral component of a plant circadian clock. LUX ARRHYTHMO (LUX) is one of the key components of the evening complex, and that play a key role in circadian rhythms and flowering. Here, we report that diverged soybean LUX has the additional role in male reproductive development. We studied diurnal and circadian rhythms of soybean LUX (GmLUXa, GmLUXb, and GmLUXc) using qRT-PCR, and show its nuclear localisation by particle bombardment. Yeast-two hybrid (Y2H) studies indicate that both GmLUXb and GmLUXc form an evening complex with GmELF4b and GmELF3a, respectively. Ectopic expression of GmLUXb in Arabidopsis lux mutants can complement functions of AtLUX, whereas GmLUXc generates novel phenotypes of serrated leaves, stunted plants, shortened anther filament, and low seed set. Overall, our results suggest that the LUX gene has diverged in soybean where GmLUXb and GmLUXc share the role to control flowering time, but GmLUXc has evolved to regulate anther filament growth and seed set by regulating the Gibberellin hormone biosynthesis pathway.

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 .

EARLY FLOWERING3 Redundancy Fine-Tunes Photoperiod Sensitivity.

Three pea (Pisum sativum) loci controlling photoperiod sensitivity, HIGH RESPONSE (HR), DIE NEUTRALIS (DNE), and STERILE NODES (SN), have recently been shown to correspond to orthologs of Arabidopsis (Arabidopsis thaliana) circadian clock genes EARLY FLOWERING3 (ELF3), ELF4, and LUX ARRHYTHMO, respectively. A fourth pea locus, PHOTOPERIOD (PPD), also contributes to the photoperiod response in a similar manner to SN and DNE, and recessive ppd mutants on a spring-flowering hr mutant background show early, photoperiod-insensitive flowering. However, the molecular identity of PPD has so far remained elusive. Here, we show that the PPD locus also has a role in maintenance of diurnal and circadian gene expression rhythms and identify PPD as an ELF3 co-ortholog, termed ELF3b Genetic interactions between pea ELF3 genes suggest that loss of PPD function does not affect flowering time in the presence of functional HR, whereas PPD can compensate only partially for the lack of HR These results provide an illustration of how gene duplication and divergence can generate potential for the emergence of more subtle variations in phenotype that may be adaptively significant.

Unique and conserved features of floral evocation in legumes.

Legumes, with their unique ability to fix atmospheric nitrogen, play a vital role in ensuring future food security and mitigating the effects of climate change because they use less fossil energy and produce less greenhouse gases compared with N-fertilized systems. Grain legumes are second only to cereal crops as a source of human and animal food, and they contribute approximately one third of the protein consumed by the human population. The productivity of seed crops, such as grain legumes, is dependent on flowering. Despite the genetic variation and importance of flowering in legume production, studies of the molecular pathways that control flowering in legumes are limited. Recent advances in genomics have revealed that legume flowering pathways are divergent from those of such model species as Arabidopsis thaliana. Here, we discuss the current understanding of flowering time regulation in legumes and highlight the unique and conserved features of floral evocation in legumes.

The Pea Photoperiod Response Gene STERILE NODES Is an Ortholog of LUX ARRHYTHMO.

The STERILE NODES (SN) locus in pea (Pisum sativum) was one of the first photoperiod response genes to be described and provided early evidence for the genetic control of long-distance signaling in flowering-time regulation. Lines homozygous for recessive sn mutations are early flowering and photoperiod insensitive, with an increased ability to promote flowering across a graft union in short-day conditions. Here, we show that SN controls developmental regulation of genes in the FT family and rhythmic regulation of genes related to circadian clock function. Using a positional and functional candidate approach, we identify SN as the pea ortholog of LUX ARRHYTHMO, a GARP transcription factor from Arabidopsis (Arabidopsis thaliana) with an important role in circadian clock function. In addition to induced mutants, sequence analysis demonstrates the presence of at least three other independent, naturally occurring loss-of-function mutations among known sn cultivars. Examination of genetic and regulatory interactions between SN and two other circadian clock genes, HIGH RESPONSE TO PHOTOPERIOD (HR) and DIE NEUTRALIS (DNE), suggests a complex relationship in which HR regulates expression of SN and the role of DNE and HR in control of flowering is dependent on SN. These results extend previous work to show that pea orthologs of all three Arabidopsis evening complex genes regulate clock function and photoperiod-responsive flowering and suggest that the function of these genes may be widely conserved.

An RNA-seq transcriptome analysis of histone modifiers and RNA silencing genes in soybean during floral initiation process.

Epigenetics has been recognised to play vital roles in many plant developmental processes, including floral initiation through the epigenetic regulation of gene expression. The histone modifying proteins that mediate these modifications involve the SET domain-containing histone methyltransferases, JmjC domain-containing demethylase, acetylases and deacetylases. In addition, RNA interference (RNAi)-associated genes are also involved in epigenetic regulation via RNA-directed DNA methylation and post-transcriptional gene silencing. Soybean, a major crop legume, requires a short day to induce flowering. How histone modifications regulate the plant response to external cues that initiate flowering is still largely unknown. Here, we used RNA-seq to address the dynamics of transcripts that are potentially involved in the epigenetic programming and RNAi mediated gene silencing during the floral initiation of soybean. Soybean is a paleopolyploid that has been subjected to at least two rounds of whole genome duplication events. We report that the expanded genomic repertoire of histone modifiers and RNA silencing genes in soybean includes 14 histone acetyltransferases, 24 histone deacetylases, 47 histone methyltransferases, 15 protein arginine methyltransferases, 24 JmjC domain-containing demethylases and 47 RNAi-associated genes. To investigate the role of these histone modifiers and RNA silencing genes during floral initiation, we compared the transcriptional dynamics of the leaf and shoot apical meristem at different time points after a short-day treatment. Our data reveal that the extensive activation of genes that are usually involved in the epigenetic programming and RNAi gene silencing in the soybean shoot apical meristem are reprogrammed for floral development following an exposure to inductive conditions.

A conserved molecular basis for photoperiod adaptation in two temperate legumes.

Legumes were among the first plant species to be domesticated, and accompanied cereals in expansion of agriculture from the Fertile Crescent into diverse environments across the Mediterranean basin, Europe, Central Asia, and the Indian subcontinent. Although several recent studies have outlined the molecular basis for domestication and eco-geographic adaptation in the two main cereals from this region, wheat and barley, similar questions remain largely unexplored in their legume counterparts. Here we identify two major loci controlling differences in photoperiod response between wild and domesticated pea, and show that one of these, high response to photoperiod (HR), is an ortholog of early flowering 3 (ELF3), a gene involved in circadian clock function. We found that a significant proportion of flowering time variation in global pea germplasm is controlled by HR, with a single, widespread functional variant conferring altered circadian rhythms and the reduced photoperiod response associated with the spring habit. We also present evidence that ELF3 has a similar role in lentil, another major legume crop, with a distinct functional variant contributing to reduced photoperiod response in cultivars widely deployed in short-season environments. Our results identify the factor likely to have permitted the successful prehistoric expansion of legume cultivation to Northern Europe, and define a conserved genetic basis for major adaptive changes in flowering phenology and growth habit in an important crop group.

The pea GIGAS gene is a FLOWERING LOCUS T homolog necessary for graft-transmissible specification of flowering but not for responsiveness to photoperiod.

Garden pea (Pisum sativum) was prominent in early studies investigating the genetic control of flowering and the role of mobile flowering signals. In view of recent evidence that genes in the FLOWERING LOCUS T (FT) family play an important role in generating mobile flowering signals, we isolated the FT gene family in pea and examined the regulation and function of its members. Comparison with Medicago truncatula and soybean (Glycine max) provides evidence of three ancient subclades (FTa, FTb, and FTc) likely to be common to most crop and model legumes. Pea FT genes show distinctly different expression patterns with respect to developmental timing, tissue specificity, and response to photoperiod and differ in their activity in transgenic Arabidopsis thaliana, suggesting they may have different functions. We show that the pea FTa1 gene corresponds to the GIGAS locus, which is essential for flowering under long-day conditions and promotes flowering under short-day conditions but is not required for photoperiod responsiveness. Grafting, expression, and double mutant analyses show that GIGAS/FTa1 regulates a mobile flowering stimulus but also provide clear evidence for a second mobile flowering stimulus that is correlated with expression of FTb2 in leaf tissue. These results suggest that induction of flowering by photoperiod in pea results from interactions among several members of a diversified FT family.

DIE NEUTRALIS and LATE BLOOMER 1 contribute to regulation of the pea circadian clock.

The DIE NEUTRALIS (DNE) locus in garden pea (Pisum sativum) was previously shown to inhibit flowering under noninductive short-day conditions and to affect a graft-transmissible flowering signal. In this study, we establish that DNE has a role in diurnal and/or circadian regulation of several clock genes, including the pea GIGANTEA (GI) ortholog LATE BLOOMER 1 (LATE1) and orthologs of the Arabidopsis thaliana genes LATE ELONGATED HYPOCOTYL and TIMING OF CHLOROPHYLL A/B BINDING PROTEIN EXPRESSION 1. We also confirm that LATE1 participates in the clock and provide evidence that DNE is the ortholog of Arabidopsis EARLY FLOWERING4 (ELF4). Circadian rhythms of clock gene expression in wild-type plants under constant light were weaker in pea than in Arabidopsis, and a number of differences were also seen in the effects of both DNE/ELF4 and LATE1/GI on clock gene expression. Grafting studies suggest that DNE controls flowering at least in part through a LATE1-dependent mobile stimulus, and dne mutants show elevated expression of a FLOWERING LOCUS T homolog under short-day conditions. However, the early flowering of the dne mutant is not associated with altered expression of a previously described CONSTANS-like gene. Collectively, our results characterize the clock system and reveal its importance for photoperiod responsiveness in a model legume.

Update on the genetic control of flowering in garden pea.

The garden pea has been a model for the genetics of flowering for several decades and numerous flowering loci have been identified, but until recently little was known about the molecular nature of these loci. This paper presents an update on recent work on the molecular genetics of flowering in pea, outlining progress in gene and mutant isolation, expression analyses, grafting and other physiological studies, and candidate gene assessment. Work so far has led to the identification of the LATE1 and DNE loci as orthologues of Arabidopsis GIGANTEA and ELF4, respectively, and candidate genes for several other loci are being evaluated. Expression analysis of an expanded FT-like gene family suggests a more complex role for this group of genes. These results provide the first insight into the circadian clock, photoperiod response mechanism, and mobile signals in pea, and identify both conserved and divergent features in comparison with Arabidopsis.

Pea LATE BLOOMER1 is a GIGANTEA ortholog with roles in photoperiodic flowering, deetiolation, and transcriptional regulation of circadian clock gene homologs.

Genes controlling the transition to flowering have been studied in several species, including Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), but have not yet received much attention in legumes. Here, we describe a new allelic series of late-flowering, photoperiod-insensitive mutants in the pea (Pisum sativum) LATE BLOOMER1 (LATE1) gene and show that LATE1 is an ortholog of Arabidopsis GIGANTEA. Mutants display defects in phytochrome B-dependent deetiolation under red light and in the diurnal regulation of pea homologs of several Arabidopsis circadian clock genes, including TIMING OF CAB1, EARLY FLOWERING4, and CIRCADIAN CLOCK ASSOCIATED1/LATE ELONGATED HYPOCOTYL. LATE1 itself shows strongly rhythmic expression with a small but distinct acute peak following dark-to-light transfer. Mutations in LATE1 prevent the induction of a FLOWERING LOCUS T (FT) homolog FTL in long days but cause only minor alteration to the rhythmic expression pattern of the only known group Ia CONSTANS homolog COLa. The late-flowering phenotype of late1 mutants can be completely rescued by grafting to the wild type, but this rescue is not associated with a significant increase in FTL transcript level in shoot apices. Genetic interactions of late1 with the photoperiod-insensitive, early-flowering sterile nodes (sn) mutant and impairment of the LATE1 diurnal expression rhythm in sn plants suggest that SN may also affect the circadian clock. These results show that several functions of Arabidopsis GIGANTEA are conserved in its pea ortholog and demonstrate that genetic pathways for photoperiodic flowering are likely to be conserved between these two species. They also suggest that in addition to its role in the floral transition, LATE1 also acts throughout reproductive development.