Expansion of the circadian transcriptome in Brassica rapa and genome-wide diversification of paralog expression patterns

Integrative transcriptomics reveals association of abscisic acid and lignin pathways with cassava whitefly resistance

BACKGROUND

Whiteflies are a global threat to crop yields, including the African subsistence crop cassava (Manihot esculenta). Outbreaks of superabundant whitefly populations throughout Eastern and Central Africa in recent years have dramatically increased the pressures of whitefly feeding and virus transmission on cassava. Whitefly-transmitted viral diseases threaten the food security of hundreds of millions of African farmers, highlighting the need for developing and deploying whitefly-resistant cassava. However, plant resistance to whiteflies remains largely poorly characterized at the genetic and molecular levels. Knowledge of cassava-defense programs also remains incomplete, limiting characterization of whitefly-resistance mechanisms. To better understand the genetic basis of whitefly resistance in cassava, we define the defense hormone- and Aleurotrachelus socialis (whitefly)-responsive transcriptome of whitefly-susceptible (COL2246) and whitefly-resistant (ECU72) cassava using RNA-seq. For broader comparison, hormone-responsive transcriptomes of Arabidopsis thaliana were also generated.

RESULTS

Whitefly infestation, salicylic acid (SA), jasmonic acid (JA), ethylene (ET), and abscisic acid (ABA) transcriptome responses of ECU72 and COL2246 were defined and analyzed. Strikingly, SA responses were largely reciprocal between the two cassava genotypes and we suggest candidate regulators. While susceptibility was associated with SA in COL2246, resistance to whitefly in ECU72 was associated with ABA, with SA-ABA antagonism observed. This was evidenced by expression of genes within the SA and ABA pathways and hormone levels during A. socialis infestation. Gene-enrichment analyses of whitefly- and hormone-responsive genes suggest the importance of fast-acting cell wall defenses (e.g., elicitor recognition, lignin biosynthesis) during early infestation stages in whitefly-resistant ECU72. A surge of ineffective immune and SA responses characterized the whitefly-susceptible COL2246's response to late-stage nymphs. Lastly, in comparison with the model plant Arabidopsis, cassava's hormone-responsive genes showed striking divergence in expression.

CONCLUSIONS

This study provides the first characterization of cassava's global transcriptome responses to whitefly infestation and defense hormone treatment. Our analyses of ECU72 and COL2246 uncovered possible whitefly resistance/susceptibility mechanisms in cassava. Comparative analysis of cassava and Arabidopsis demonstrated that defense programs in Arabidopsis may not always mirror those in crop species. More broadly, our hormone-responsive transcriptomes will also provide a baseline for the cassava community to better understand global responses to other yield-limiting pests/pathogens.

Diurnal oscillations of epigenetic modifications are associated with variation in rhythmic expression of homoeologous genes in Brassica napus

BACKGROUND

Epigenetic modifications that exhibit circadian oscillations also promote circadian oscillations of gene expression. Brassica napus is a heterozygous polyploid species that has undergone distant hybridization and genome doubling events and has a young and distinct species origin. Studies incorporating circadian rhythm analysis of epigenetic modifications can offer new insights into differences in diurnal oscillation behavior among subgenomes and the regulation of diverse expressions of homologous gene rhythms in biological clocks.

RESULTS

In this study, we created a high-resolution and multioscillatory gene expression dataset, active histone modification (H3K4me3, H3K9ac), and RNAPII recruitment in Brassica napus. We also conducted the pioneering characterization of the diurnal rhythm of transcription and epigenetic modifications in an allopolyploid species. We compared the evolution of diurnal rhythms between subgenomes and observed that the Cn subgenome had higher diurnal oscillation activity in both transcription and active histone modifications than the An subgenome. Compared to the A subgenome in Brassica rapa, the An subgenome of Brassica napus displayed significant changes in diurnal oscillation characteristics of transcription. Homologous gene pairs exhibited a higher proportion of diurnal oscillation in transcription than subgenome-specific genes, attributed to higher chromatin accessibility and abundance of active epigenetic modification types. We found that the diurnal expression of homologous genes displayed diversity, and the redundancy of the circadian system resulted in extensive changes in the diurnal rhythm characteristics of clock genes after distant hybridization and genome duplication events. Epigenetic modifications influenced the differences in the diurnal rhythm of homologous gene expression, and the diurnal oscillation of homologous gene expression was affected by the combination of multiple histone modifications.

CONCLUSIONS

Herein, we presented, for the first time, a characterization of the diurnal rhythm characteristics of gene expression and its epigenetic modifications in an allopolyploid species. Our discoveries shed light on the epigenetic factors responsible for the diurnal oscillation activity imbalance between subgenomes and homologous genes' rhythmic expression differences. The comprehensive time-series dataset we generated for gene expression and epigenetic modifications provides a valuable resource for future investigations into the regulatory mechanisms of protein-coding genes in Brassica napus.

Genome-wide circadian gating of a cold temperature response in bread wheat

Circadian rhythms coordinate the responses of organisms with their daily fluctuating environments, by establishing a temporal program of gene expression. This schedules aspects of metabolism, physiology, development and behaviour according to the time of day. Circadian regulation in plants is extremely pervasive, and is important because it underpins both productivity and seasonal reproduction. Circadian regulation extends to the control of environmental responses through a regulatory process known as circadian gating. Circadian gating is the process whereby the circadian clock regulates the response to an environmental cue, such that the magnitude of response to an identical cue varies according to the time of day of the cue. Here, we show that there is genome-wide circadian gating of responses to cold temperatures in plants. By using bread wheat as an experimental model, we establish that circadian gating is crucial to the programs of gene expression that underlie the environmental responses of a crop of major socioeconomic importance. Furthermore, we identify that circadian gating of cold temperature responses are distributed unevenly across the three wheat subgenomes, which might reflect the geographical origins of the ancestors of modern wheat.

A novel CLAVATA1 mutation causes multilocularity in Brassica rapa

Locules are the seed-bearing structure of fruits. Multiple locules are associated with increased fruit size and seed set, and therefore, control of locule number is an important agronomic trait. Locule number is controlled in part by the CLAVATA-WUSCHEL pathway. Disruption of either the CLAVATA1 receptor-like kinase or its ligand CLAVATA3 can cause larger floral meristems and an increased number of locules. In an EMS mutagenized population of Brassica rapa, we identified a mutant allele that raises the number of locules from four to a range of from six to eight. Linkage mapping and genetic analysis support that the mutant phenotype is due to a missense mutation in a CLAVATA 1 (CLV1) homolog. In addition to increased locule number, additional internal gynoecia are formed in brclv1 individuals, suggesting a failure to terminate floral meristem development, which results in decreased seed production.

Functional characterization of C-TERMINALLY ENCODED PEPTIDE (CEP) family in Brassica rapa L

The small regulatory C-TERMINALLY ENCODED PEPTIDE (CEP) peptide family plays crucial roles in plant growth and stress response. However, little is known about this peptide family in Brassica species. Here, we performed a systematic analysis to identify the putative Brassica rapa L. CEP (BrCEP) gene family. In total, 27 BrCEP genes were identified and they were classified into four subgroups based on the CEP motifs similarity. BrCEP genes displayed distinct expression patterns in response to both developmental and several environmental signals, suggesting their broad roles during Brassica rapa development. Furthuremore, the synthetic BrCEP3 peptide accelerated Brassica rapa primary root growth in a hydrogen peroxide (H₂O₂) and Ca²⁺ dependent manner. In summary, our work will provide fundamental insights into the physiological function of CEP peptides during Brassica rapa development.

Time of Day Analysis over a Field Grown Developmental Time Course in Rice

Plants integrate time of day (TOD) information over an entire season to ensure optimal growth, flowering time, and grain fill. However, most TOD expression studies have focused on a limited number of combinations of daylength and temperature under laboratory conditions. Here, an Oryza sativa (rice) expression study that followed TOD expression in the field over an entire growing season was re-analyzed. Similar to Arabidopsis thaliana, almost all rice genes have a TOD-specific expression over the developmental time course. As has been suggested in other grasses, thermocycles were a stronger cue for TOD expression than the photocycles over the growing season. All the core circadian clock genes display consistent TOD expression over the season with the interesting exception that the two grass paralogs of EARLY FLOWERING 3 (ELF3) display a distinct phasing based on the interaction between thermo- and photo-cycles. The dataset also revealed how specific pathways are modulated to distinct TOD over the season consistent with the changing biology. The data presented here provide a resource for researchers to study how TOD expression changes under natural conditions over a developmental time course, which will guide approaches to engineer more resilient and prolific crops.

The wheat clock strikes a balance across subgenomes to regulate gene expression

Increasing crop yields is complicated by the polyploid nature of our major crops. This Primer explores the implications of a PLOS Biology study that provides a transcriptomic view of the influence of the circadian clock on agriculturally relevant traits in the polyploid bread wheat.

Growth-limiting drought stress induces time-of-day-dependent transcriptome and physiological responses in hybrid poplar

Drought stress negatively impacts the health of long-lived trees. Understanding the genetic mechanisms that underpin response to drought stress is requisite for selecting or enhancing climate change resilience. We aimed to determine how hybrid poplars respond to prolonged and uniform exposure to drought; how responses to moderate and more severe growth-limiting drought stresses differed; and how drought responses change throughout the day. We established hybrid poplar trees (Populus × 'Okanese') from unrooted stem cutting with abundant soil moisture for 6 weeks. We then withheld water to establish well-watered, moderate and severe growth-limiting drought conditions. These conditions were maintained for 3 weeks during which growth was monitored. We then measured photosynthetic rates and transcriptomes of leaves that had developed during the drought treatments at two times of day. The moderate and severe drought treatments elicited distinct changes in growth and development, photosynthetic rates and global transcriptome profiles. Notably, the time of day of sampling produced the strongest effect in the transcriptome data. The moderate drought treatment elicited global transcriptome changes that were intermediate to the severe and well-watered treatments in the early evening but did not elicit a strong drought response in the morning. Stable drought conditions that are sufficient to limit plant growth elicit distinct transcriptional profiles depending on the degree of water limitation and on the time of day at which they are measured. There appears to be a limited number of genes and functional gene categories that are responsive to all of the tested drought conditions in this study emphasizing the complex nature of drought regulation in long-lived trees.

Circadian regulation of the transcriptome in a complex polyploid crop

The circadian clock is a finely balanced timekeeping mechanism that coordinates programmes of gene expression. It is currently unknown how the clock regulates expression of homoeologous genes in polyploids. Here, we generate a high-resolution time-course dataset to investigate the circadian balance between sets of 3 homoeologous genes (triads) from hexaploid bread wheat. We find a large proportion of circadian triads exhibit imbalanced rhythmic expression patterns, with no specific subgenome favoured. In wheat, period lengths of rhythmic transcripts are found to be longer and have a higher level of variance than in other plant species. Expression of transcripts associated with circadian controlled biological processes is largely conserved between wheat and Arabidopsis; however, striking differences are seen in agriculturally critical processes such as starch metabolism. Together, this work highlights the ongoing selection for balance versus diversification in circadian homoeologs and identifies clock-controlled pathways that might provide important targets for future wheat breeding.

CAST-R: An application to visualize circadian and heat stress-responsive genes in plants

The circadian clock helps organisms to anticipate and coordinate gene regulatory responses to changes in environmental stimuli. Under stresses, both time of day and the circadian clock closely control the magnitude of plant responses. The identification of clock-regulated genes is, therefore, important when studying the influence of environmental factors. Here, we present CAST-R (Circadian And heat STress-Responsive), a "Shiny" application that allows users to identify and visualize circadian and heat stress-responsive genes in plants. More specifically, users can generate and export profiles and heatmaps representing transcript abundance of a single or of multiple Arabidopsis (Arabidopsis thaliana) genes over a 24-h time course, in response to heat stress and during recovery following the stress. The application also takes advantage of published Arabidopsis chromatin immunoprecipitation-sequencing datasets to visualize the connections between clock proteins and their targets in an interactive network. In addition, CAST-R offers the possibility to perform phase (i.e. timing of expression) enrichment analyses for rhythmic datasets from any species, within and beyond plants. This functionality combines statistical analyses and graphical representations to identify significantly over- and underrepresented phases within a subset of genes. Lastly, profiles of transcript abundance can be visualized from multiple circadian datasets generated in Arabidopsis, Brassica rapa, barley (Hordeum vulgare), and rice (Oryza sativa). In summary, CAST-R is a user-friendly interface that allows the rapid identification of circadian and stress-responsive genes through multiple modules of visualization. We anticipate that this tool will make it easier for users to obtain temporal and dynamic information on genes of interest that links plant responses to environmental signals.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.

The adaptive nature of the plant circadian clock in natural environments

The plant circadian clock coordinates developmental, physiological, and metabolic processes with diel changes in light and temperature throughout the year. The balance between the persistence and plasticity of the clock in response to predictable and unpredictable environmental changes may be key to the clock's adaptive nature across temporal and spatial scales. Studies under controlled conditions have uncovered critical signaling pathways involved in light and temperature perception by the clock; however, they don't account for the natural lag of temperature behind photoperiod. Studies in natural environments provide key insights into the clock's adaptive advantage under more complex natural settings. Here, we discuss the role of the circadian clock in light and temperature perception and signaling, how the clock integrates these signals for a coordinated and adaptive response, and the adaptive advantage conferred by the clock across time and space in natural environments.

Keeping time in the dark: Potato diel and circadian rhythmic gene expression reveals tissue-specific circadian clocks

The circadian clock is an internal molecular oscillator and coordinates numerous physiological processes through regulation of molecular pathways. Tissue-specific clocks connected by mobile signals have previously been found to run at different speeds in Arabidopsis thaliana tissues. However, tissue variation in circadian clocks in crop species is unknown. In this study, leaf and tuber global gene expression in cultivated potato under cycling and constant environmental conditions was profiled. In addition, we used a circadian-regulated luciferase reporter construct to study tuber gene expression rhythms. Diel and circadian expression patterns were present among 17.9% and 5.6% of the expressed genes in the tuber. Over 500 genes displayed differential tissue specific diel phases. Intriguingly, few core circadian clock genes had circadian expression patterns, while all such genes were circadian rhythmic in cultivated tomato leaves. Furthermore, robust diel and circadian transcriptional rhythms were observed among detached tubers. Our results suggest alternative regulatory mechanisms and/or clock composition is present in potato, as well as the presence of tissue-specific independent circadian clocks. We have provided the first evidence of a functional circadian clock in below-ground storage organs, holding important implications for other storage root and tuberous crops.

The biology of time: dynamic responses of cell types to developmental, circadian and environmental cues

As sessile organisms, plants are finely tuned to respond dynamically to developmental, circadian and environmental cues. Genome-wide studies investigating these types of cues have uncovered the intrinsically different ways they can impact gene expression over time. Recent advances in single-cell sequencing and time-based bioinformatic algorithms are now beginning to reveal the dynamics of these time-based responses within individual cells and plant tissues. Here, we review what these techniques have revealed about the spatiotemporal nature of gene regulation, paying particular attention to the three distinct ways in which plant tissues are time sensitive. (i) First, we discuss how studying plant cell identity can reveal developmental trajectories hidden in pseudotime. (ii) Next, we present evidence that indicates that plant cell types keep their own local time through tissue-specific regulation of the circadian clock. (iii) Finally, we review what determines the speed of environmental signaling responses, and how they can be contingent on developmental and circadian time. By these means, this review sheds light on how these different scales of time-based responses can act with tissue and cell-type specificity to elicit changes in whole plant systems.

Circadian coordination of cellular processes and abiotic stress responses

Diel changes in the environment are perceived by the circadian clock which transmits temporal information throughout the plant cell to synchronize daily and seasonal environmental signals with internal biological processes. Dynamic modulations of diverse levels of clock gene regulation within the plant cell are impacted by stress. Recent insights into circadian control of cellular processes such as alternative splicing, polyadenylation, and noncoding RNAs are discussed. We highlight studies on the circadian regulation of reactive oxygen species, calcium signaling, and gating of temperature stress responses. Finally, we briefly summarize recent work on the translation-specific rhythmicity of cell cycle genes and the control of subcellular localization and relocalization of oscillator components. Together, this mini-review highlights these cellular events in the context of clock gene regulation and stress responses in Arabidopsis.

QTL Mapping and Diurnal Transcriptome Analysis Identify Candidate Genes Regulating Brassica napus Flowering Time

Timely flowering is important for seed formation and maximization of rapeseed (Brassica napus) yield. Here, we performed flowering-time quantitative trait loci (QTL) mapping using a double haploid (DH) population grown in three environments to study the genetic architecture. Brassica 60 K Illumina Infinium™ single nucleotide polymorphism (SNP) array and simple sequence repeat (SSR) markers were used for genotyping of the DH population, and a high-density genetic linkage map was constructed. QTL analysis of flowering time from the three environments revealed five consensus QTLs, including two major QTLs. A major QTL located on chromosome A03 was detected specifically in the semi-winter rapeseed growing region, and the one on chromosome C08 was detected in all environments. Ribonucleic acid sequencing (RNA-seq) was performed on the parents’ leaves at seven time-points in a day to determine differentially expressed genes (DEGs). The biological processes and pathways with significant enrichment of DEGs were obtained. The DEGs in the QTL intervals were analyzed, and four flowering time-related candidate genes were found. These results lay a foundation for the genetic regulation of rapeseed flowering time and create a rapeseed gene expression library for seven time-points in a day.

Abiotic stress through time

Throughout plant evolution the circadian clock has expanded into a complex signaling network, coordinating physiological and metabolic processes with the environment. Early land plants faced new environmental pressures that required energy-demanding stress responses. Integrating abiotic stress response into the circadian system provides control over daily energy expenditure. Here, we describe the evolution of the circadian clock in plants and the limited, yet compelling, evidence for conserved regulation of abiotic stress. The need to introduce abiotic stress tolerance into current crops has expanded research into wild accessions and revealed extensive variation in circadian clock parameters across monocot and eudicot species. We argue that research into the ancestral links between the clock and abiotic stress will benefit crop improvement efforts.

Closing the protein gap in plant chronobiology

Our modern understanding of diel cell regulation in plants stems from foundational work in the late 1990s that analysed the dynamics of selected genes and mutants in Arabidopsis thaliana. The subsequent rise of transcriptomics technologies such as microarrays and RNA sequencing has substantially increased our understanding of anticipatory (circadian) and reactive (light- or dark-triggered) diel events in plants. However, it is also becoming clear that gene expression data fail to capture critical events in diel regulation that can only be explained by studying protein-level dynamics. Over the past decade, mass spectrometry technologies and quantitative proteomic workflows have significantly advanced, finally allowing scientists to characterise diel protein regulation at high throughput. Initial proteomic investigations suggest that the diel transcriptome and proteome generally lack synchrony and that the timing of daily regulatory events in plants is impacted by multiple levels of protein regulation (e.g., post-translational modifications [PTMs] and protein-protein interactions [PPIs]). Here, we highlight and summarise how the use of quantitative proteomics to elucidate diel plant cell regulation has advanced our understanding of these processes. We argue that this new understanding, coupled with the extraordinary developments in mass spectrometry technologies, demands greater focus on protein-level regulation of, and by, the circadian clock. This includes hitherto unexplored diel dynamics of protein turnover, PTMs, protein subcellular localisation and PPIs that can be masked by simple transcript- and protein-level changes. Finally, we propose new directions for how the latest advancements in quantitative proteomics can be utilised to answer outstanding questions in plant chronobiology.

Expansion of the circadian transcriptome in Brassica rapa and genome-wide diversification of paralog expression patterns

An important challenge of crop improvement strategies is assigning function to paralogs in polyploid crops. Here we describe the circadian transcriptome in the polyploid crop Brassica rapa. Strikingly, almost three-quarters of the expressed genes exhibited circadian rhythmicity. Genetic redundancy resulting from whole genome duplication is thought to facilitate evolutionary change through sub- and neo-functionalization among paralogous gene pairs. We observed genome-wide expansion of the circadian expression phase among retained paralogous pairs. Using gene regulatory network models, we compared transcription factor targets between B. rapa and Arabidopsis circadian networks to reveal evidence for divergence between B. rapa paralogs that may be driven in part by variation in conserved non-coding sequences (CNS). Additionally, differential drought response among retained paralogous pairs suggests further functional diversification. These findings support the rapid expansion and divergence of the transcriptional network in a polyploid crop and offer a new approach for assessing paralog activity at the transcript level.