Usability of reference-free transcriptome assemblies for detection of differential expression: a case study on Aethionema arabicum dimorphic seeds

The dimorphic diaspore model Aethionema arabicum (Brassicaceae): Distinct molecular and morphological control of responses to parental and germination temperatures

Plants in habitats with unpredictable conditions often have diversified bet-hedging strategies that ensure fitness over a wider range of variable environmental factors. A striking example is the diaspore (seed and fruit) heteromorphism that evolved to maximize species survival in Aethionema arabicum (Brassicaceae) in which external and endogenous triggers allow the production of two distinct diaspores on the same plant. Using this dimorphic diaspore model, we identified contrasting molecular, biophysical, and ecophysiological mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained by pericarp (fruit coat) removal from IND fruits. Large-scale comparative transcriptome and hormone analyses of M+ seeds, IND fruits, and M- seeds provided comprehensive datasets for their distinct thermal responses. Morph-specific differences in co-expressed gene modules in seeds, as well as in seed and pericarp hormone contents, identified a role of the IND pericarp in imposing coat dormancy by generating hypoxia affecting abscisic acid (ABA) sensitivity. This involved expression of morph-specific transcription factors, hypoxia response, and cell wall remodeling genes, as well as altered ABA metabolism, transport, and signaling. Parental temperature affected ABA contents and ABA-related gene expression and altered IND pericarp biomechanical properties. Elucidating the molecular framework underlying the diaspore heteromorphism can provide insight into developmental responses to globally changing temperatures.

Aethionema arabicum dimorphic seed trait resetting during transition to seedlings

The transition from germinating seeds to emerging seedlings is one of the most vulnerable plant life cycle stages. Heteromorphic diaspores (seed and fruit dispersal units) are an adaptive bet-hedging strategy to cope with spatiotemporally variable environments. While the roles and mechanisms of seedling traits have been studied in monomorphic species, which produce one type of diaspore, very little is known about seedlings in heteromorphic species. Using the dimorphic diaspore model Aethionema arabicum (Brassicaceae), we identified contrasting mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained from IND fruits by pericarp (fruit coat) removal. What follows the completion of germination is the pre-emergence seedling growth phase, which we investigated by comparative growth assays of early seedlings derived from the M+ seeds, bare M- seeds, and IND fruits. The dimorphic seedlings derived from M+ and M- seeds did not differ in their responses to ambient temperature and water potential. The phenotype of seedlings derived from IND fruits differed in that they had bent hypocotyls and their shoot and root growth was slower, but the biomechanical hypocotyl properties of 15-day-old seedlings did not differ between seedlings derived from germinated M+ seeds, M- seeds, or IND fruits. Comparison of the transcriptomes of the natural dimorphic diaspores, M+ seeds and IND fruits, identified 2,682 differentially expressed genes (DEGs) during late germination. During the subsequent 3 days of seedling pre-emergence growth, the number of DEGs was reduced 10-fold to 277 root DEGs and 16-fold to 164 shoot DEGs. Among the DEGs in early seedlings were hormonal regulators, in particular for auxin, ethylene, and gibberellins. Furthermore, DEGs were identified for water and ion transporters, nitrate transporter and assimilation enzymes, and cell wall remodeling protein genes encoding enzymes targeting xyloglucan and pectin. We conclude that the transcriptomes of seedlings derived from the dimorphic diaspores, M+ seeds and IND fruits, undergo transcriptional resetting during the post-germination pre-emergence growth transition phase from germinated diaspores to growing seedlings.

Single seeds exhibit transcriptional heterogeneity during secondary dormancy induction

Seeds are highly resilient to the external environment, which allows plants to persist in unpredictable and unfavorable conditions. Some plant species have adopted a bet-hedging strategy to germinate a variable fraction of seeds in any given condition, and this could be explained by population-based threshold models. Here, in the model plant Arabidopsis (Arabidopsis thaliana), we induced secondary dormancy (SD) to address the transcriptional heterogeneity among seeds that leads to binary germination/nongermination outcomes. We developed a single-seed RNA-seq strategy that allowed us to observe a reduction in seed transcriptional heterogeneity as seeds enter stress conditions, followed by an increase during recovery. We identified groups of genes whose expression showed a specific pattern through a time course and used these groups to position the individual seeds along the transcriptional gradient of germination competence. In agreement, transcriptomes of dormancy-deficient seeds (mutant of DELAY OF GERMINATION 1) showed a shift toward higher values of the germination competence index. Interestingly, a significant fraction of genes with variable expression encoded translation-related factors. In summary, interrogating hundreds of single-seed transcriptomes during SD-inducing treatment revealed variability among the transcriptomes that could result from the distribution of population-based sensitivity thresholds. Our results also showed that single-seed RNA-seq is the method of choice for analyzing seed bet-hedging-related phenomena.

Cold-induced secondary dormancy and its regulatory mechanisms in Beta vulgaris

The dynamic behaviour of seeds in soil seed banks depends on their ability to act as sophisticated environmental sensors to adjust their sensitivity thresholds for germination by dormancy mechanisms. Here we show that prolonged incubation of sugar beet fruits at low temperature (chilling at 5°C, generally known to release seed dormancy of many species) can induce secondary nondeep physiological dormancy of an apparently nondormant crop species. The physiological and biophysical mechanisms underpinning this cold-induced secondary dormancy include the chilling-induced accumulation of abscisic acid in the seeds, a reduction in the embryo growth potential and a block in weakening of the endosperm covering the embryonic root. Transcriptome analysis revealed distinct gene expression patterns in the different temperature regimes and upon secondary dormancy induction and maintenance. The chilling caused reduced expression of cell wall remodelling protein genes required for embryo cell elongation growth and endosperm weakening, as well as increased expression of seed maturation genes, such as for late embryogenesis abundant proteins. A model integrating the hormonal signalling and master regulator expression with the temperature-control of seed dormancy and maturation programmes is proposed. The revealed mechanisms of the cold-induced secondary dormancy are important for climate-smart agriculture and food security.

Xyloglucan remodelling enzymes and the mechanics of plant seed and fruit biology

This article comments on: Di Marzo M, Ebeling Viana V, Banfi C, Cassina V, Corti R, Herrera-Ubaldo H, Babolin N, Guazzotti A, Kiegle E, Gregis V, de Folter S, Sampedro J, Mantegazza F, Colombo L, Ezquer I. 2022. Cell wall modifications by α-XYLOSIDASE1 are required for the control of seed and fruit size. Journal of Experimental Botany 73, 1499–1515.

Temporal Control of Seed Development in Dicots: Molecular Bases, Ecological Impact and Possible Evolutionary Ramifications

In flowering plants, seeds serve as organs of both propagation and dispersal. The developing seed passes through several consecutive stages, following a conserved general outline. The overall time needed for a seed to develop, however, may vary both within and between plant species, and these temporal developmental properties remain poorly understood. In the present paper, we summarize the existing data for seed development alterations in dicot plants. For genetic mutations, the reported cases were grouped in respect of the key processes distorted in the mutant specimens. Similar phenotypes arising from the environmental influence, either biotic or abiotic, were also considered. Based on these data, we suggest several general trends of timing alterations and how respective mechanisms might add to the ecological plasticity of the families considered. We also propose that the developmental timing alterations may be perceived as an evolutionary substrate for heterochronic events. Given the current lack of plausible models describing timing control in plant seeds, the presented suggestions might provide certain insights for future studies in this field.

A tale of two morphs: developmental patterns and mechanisms of seed coat differentiation in the dimorphic diaspore model Aethionema arabicum (Brassicaceae)

The developmental transition from a fertilized ovule to a dispersed diaspore (seed or fruit) involves complex differentiation processes of the ovule's integuments leading to the diversity in mature seed coat structures in angiosperms. In this study, comparative imaging and transcriptome analysis were combined to investigate the morph-specific developmental differences during outer seed coat differentiation and mucilage production in Aethionema arabicum, the Brassicaceae model for diaspore dimorphism. One of the intriguing adaptations of this species is the production and dispersal of morphologically distinct, mucilaginous and non-mucilaginous diaspores from the same plant (dimorphism). The dehiscent fruit morph programme producing multiple mucilaginous seed diaspores was used as the default trait combination, similar to Arabidopsis thaliana, and was compared with the indehiscent fruit morph programme leading to non-mucilaginous diaspores. Synchrotron-based radiation X-ray tomographic microscopy revealed a co-ordinated framework of morph-specific early changes in internal anatomy of developing A. arabicum gynoecia including seed abortion in the indehiscent programme and mucilage production by the mucilaginous seed coat. The associated comparative analysis of the gene expression patterns revealed that the unique seed coat dimorphism of Ae. arabicum provides an excellent model system for comparative study of the control of epidermal cell differentiation and mucilage biosynthesis by the mucilage transcription factor cascade and their downstream cell wall and mucilage remodelling genes. Elucidating the underlying molecular framework of the dimorphic diaspore syndrome is key to understanding differential regulation of bet-hedging survival strategies in challenging environments, timely in the face of global climatic change.

The EpiDiverse Plant Epigenome-Wide Association Studies (EWAS) Pipeline

Bisulfite sequencing is a widely used technique for determining DNA methylation and its relationship with epigenetics, genetics, and environmental parameters. Various techniques were implemented for epigenome-wide association studies (EWAS) to reveal meaningful associations; however, there are only very few plant studies available to date. Here, we developed the EpiDiverse EWAS pipeline and tested it using two plant datasets, from P. abies (Norway spruce) and Q. lobata (valley oak). Hence, we present an EWAS implementation tested for non-model plant species and describe its use.

Aethionema arabicum genome annotation using PacBio full-length transcripts provides a valuable resource for seed dormancy and Brassicaceae evolution research

Aethionema arabicum is an important model plant for Brassicaceae trait evolution, particularly of seed (development, regulation, germination, dormancy) and fruit (development, dehiscence mechanisms) characters. Its genome assembly was recently improved but the gene annotation was not updated. Here, we improved the Ae. arabicum gene annotation using 294 RNA-seq libraries and 136 307 full-length PacBio Iso-seq transcripts, increasing BUSCO completeness by 11.6% and featuring 5606 additional genes. Analysis of orthologs showed a lower number of genes in Ae. arabicum than in other Brassicaceae, which could be partially explained by loss of homeologs derived from the At-α polyploidization event and by a lower occurrence of tandem duplications after divergence of Aethionema from the other Brassicaceae. Benchmarking of MADS-box genes identified orthologs of FUL and AGL79 not found in previous versions. Analysis of full-length transcripts related to ABA-mediated seed dormancy discovered a conserved isoform of PIF6-β and antisense transcripts in ABI3, ABI4 and DOG1, among other cases found of different alternative splicing between Turkey and Cyprus ecotypes. The presented data allow alternative splicing mining and proposition of numerous hypotheses to research evolution and functional genomics. Annotation data and sequences are available at the Ae. arabicum DB (https://plantcode.online.uni-marburg.de/aetar_db).

Ecological, (epi)genetic and physiological aspects of bet-hedging in angiosperms

KEY MESSAGE

Bet-hedging is a complex evolutionary strategy involving morphological, eco-physiological, (epi)genetic and population dynamics aspects. We review these aspects in flowering plants and propose further research needed for this topic. Bet-hedging is an evolutionary strategy that reduces the temporal variance in fitness at the expense of a lowered arithmetic mean fitness. It has evolved in organisms subjected to variable cues from the external environment, be they abiotic or biotic stresses such as irregular rainfall or predation. In flowering plants, bet-hedging is exhibited by hundreds of species and is mainly exerted by reproductive organs, in particular seeds but also embryos and fruits. The main example of bet-hedging in angiosperms is diaspore heteromorphism in which the same individual produces different seed/fruit morphs in terms of morphology, dormancy, eco-physiology and/or tolerance to biotic and abiotic stresses in order to 'hedge its bets' in unpredictable environments. The objective of this review is to provide a comprehensive overview of the ecological, genetic, epigenetic and physiological aspects involved in shaping bet-hedging strategies, and how these can affect population dynamics. We identify several open research questions about bet-hedging strategies in plants: 1) understanding ecological trade-offs among different traits; 2) producing more comprehensive phylogenetic analyses to understand the diffusion and evolutionary implications of this strategy; 3) clarifying epigenetic mechanisms related to bet-hedging and plant responses to environmental cues; and 4) applying multi-omics approaches to study bet-hedging at different levels of detail. Clarifying those aspects of bet-hedging will deepen our understanding of this fascinating evolutionary strategy.

The intergenerational toxic effects on offspring of medaka fish Oryzias melastigma from parental benzo[a]pyrene exposure via interference of the circadian rhythm

Benzo[a]pyrene (BaP), a widely existed polycyclic aromatic hydrocarbon pollutant in aquatic environment, has toxic effects on marine animals and their generations, but the intergenerational immunotoxic mechanism underlying has not been clearly understood. In the study, the offspring of marine medaka (oryzias melastigma) which were exposed to 0.5 μg L^-1 BaP suffered from circadian rhythm oscillation disorders and severe DNA damage. Many clock-associated genes like per1 were significantly modulated in offspring, both per1 and p53 were significantly inhibited that altered the progression of cell cycle and inhibited DNA repair, which possibly resulted in the increased mortality of offspring. The hypermethylation of the per1 promotor and abnormal levels of N⁶-methyladenosine (m⁶A) suggested that the underlying mechanism was probably related to the epigenetic modification. Moreover, the offspring from paternal BaP exposure had more severe DNA damage and a higher degree of hypermethylation than those from maternal exposure. F1 larvae from BaP-exposed parents were more sensitive to BaP exposure, showing that the expression of immune and metabolism-related genes were significantly up-regulated. Taken together, the parental toxicity induced by BaP could be passed to F1 generation and the mechanism underlying was probably associated with a characteristic circadian rhythm disorder.

Differential Gene Expression with an Emphasis on Floral Organ Size Differences in Natural and Synthetic Polyploids of Nicotiana tabacum (Solanaceae)

Floral organ size, especially the size of the corolla, plays an important role in plant reproduction by facilitating pollination efficiency. Previous studies have outlined a hypothesized organ size pathway. However, the expression and function of many of the genes in the pathway have only been investigated in model diploid species; therefore, it is unknown how these genes interact in polyploid species. Although correlations between ploidy and cell size have been shown in many systems, it is unclear whether there is a difference in cell size between naturally occurring and synthetic polyploids. To address these questions comparing floral organ size and cell size across ploidy, we use natural and synthetic polyploids of Nicotiana tabacum (Solanaceae) as well as their known diploid progenitors. We employ a comparative transcriptomics approach to perform analyses of differential gene expression, focusing on candidate genes that may be involved in floral organ size, both across developmental stages and across accessions. We see differential expression of several known floral organ candidate genes including ARF2, BIG BROTHER, and GASA/GAST1. Results from linear models show that ploidy, cell width, and cell number positively influence corolla tube circumference; however, the effect of cell width varies by ploidy, and diploids have a significantly steeper slope than both natural and synthetic polyploids. These results demonstrate that polyploids have wider cells and that polyploidy significantly increases corolla tube circumference.

Hybrid transcriptome sequencing approach improved assembly and gene annotation in Cynara cardunculus (L.)

Background

The investigation of transcriptome profiles using short reads in non-model organisms, which lack of well-annotated genomes, is limited by partial gene reconstruction and isoform detection. In contrast, long-reads sequencing techniques revealed their potential to generate complete transcript assemblies even when a reference genome is lacking. Cynara cardunculus var. altilis (DC) (cultivated cardoon) is a perennial hardy crop adapted to dry environments with many industrial and nutraceutical applications due to the richness of secondary metabolites mostly produced in flower heads. The investigation of this species benefited from the recent release of a draft genome, but the transcriptome profile during the capitula formation still remains unexplored. In the present study we show a transcriptome analysis of vegetative and inflorescence organs of cultivated cardoon through a novel hybrid RNA-seq assembly approach utilizing both long and short RNA-seq reads.

Results

The inclusion of a single Nanopore flow-cell output in a hybrid sequencing approach determined an increase of 15% complete assembled genes and 18% transcript isoforms respect to short reads alone. Among 25,463 assembled unigenes, we identified 578 new genes and updated 13,039 gene models, 11,169 of which were alternatively spliced isoforms. During capitulum development, 3424 genes were differentially expressed and approximately two-thirds were identified as transcription factors including bHLH, MYB, NAC, C2H2 and MADS-box which were highly expressed especially after capitulum opening. We also show the expression dynamics of key genes involved in the production of valuable secondary metabolites of which capitulum is rich such as phenylpropanoids, flavonoids and sesquiterpene lactones. Most of their biosynthetic genes were strongly transcribed in the flower heads with alternative isoforms exhibiting differentially expression levels across the tissues.

Conclusions

This novel hybrid sequencing approach allowed to improve the transcriptome assembly, to update more than half of annotated genes and to identify many novel genes and different alternatively spliced isoforms. This study provides new insights on the flowering cycle in an Asteraceae plant, a valuable resource for plant biology and breeding in Cynara and an effective method for improving gene annotation.

Genome Improvement and Genetic Map Construction for Aethionema arabicum, the First Divergent Branch in the Brassicaceae Family

The genus Aethionema is a sister-group to the core-group of the Brassicaceae family that includes Arabidopsis thaliana and the Brassica crops. Thus, Aethionema is phylogenetically well-placed for the investigation and understanding of genome and trait evolution across the family. We aimed to improve the quality of the reference genome draft version of the annual species Aethionema arabicum. Second, we constructed the first Ae. arabicum genetic map. The improved reference genome and genetic map enabled the development of each other. We started with the initially published genome (version 2.5). PacBio and MinION sequencing together with genetic map v2.5 were incorporated to produce the new reference genome v3.0. The improved genome contains 203 MB of sequence, with approximately 94% of the assembly made up of called (non-gap) bases, assembled into 2,883 scaffolds (with only 6% of the genome made up of non-called bases (Ns)). The N₅₀ (10.3 MB) represents an 80-fold increase over the initial genome release. We generated a Recombinant Inbred Line (RIL) population that was derived from two ecotypes: Cyprus and Turkey (the reference genotype. Using a Genotyping by Sequencing (GBS) approach, we generated a high-density genetic map with 749 (v2.5) and then 632 SNPs (v3.0) was generated. The genetic map and reference genome were integrated, thus greatly improving the scaffolding of the reference genome into 11 linkage groups. We show that long-read sequencing data and genetics are complementary, resulting in an improved genome assembly in Ae. arabicum. They will facilitate comparative genetic mapping work for the Brassicaceae family and are also valuable resources to investigate wide range of life history traits in Aethionema.