Transcriptome profiling of sorted endoreduplicated nuclei from tomato fruits: how the global shift in expression ascribed to DNA ploidy influences RNA-Seq data normalization and interpretation

CRISPR/Cas9 edited SlGT30 improved both drought resistance and fruit yield through endoreduplication

There is often a trade-off effect between different agronomic traits due to gene pleiotropy, leading to a negative correlation between yield and resistance. Consequently, using gene-editing techniques to develop superior traits becomes challenging. Genetic resources that defy this constraint are scarce but hold great potential as targets for improvement through the utilisation of CRISPR. Transcription factors are critical in modulating numerous gene expressions across diverse biological processes. Here, we found that the trihelix transcription factor SlGT30 plays a role in drought resistance and tomato fruit development. We edited the SlGT30 gene with CRISPR/Cas9 technology and found that the knockout lines showed decreased stomata density in the leaves and large fruits. Subsequent examination revealed that cell ploidy was impacted in the leaves and fruits of SlGT30 knockout lines. SlGT30 knockout affected cell size through the endoreduplication pathway, manifested in decreased stomata density and reduced water loss. Consequently, this resulted in an enhancement of drought resistance. For the fruit, both cell size and cell number increased in the fruit pericarp of knockout lines, improving the fruit size and weight accordingly. Therefore, SlGT30 represents a promising candidate gene for gene editing in breeding practice.

A stable combination of non-stable genes outperforms standard reference genes for RT-qPCR data normalization

Gene expression profiling is of key importance in all domains of life sciences, as medicine, environment, and plants, for both basic and applied research. Despite the emergence of microarrays and high-throughput sequencing, qPCR remains a standard method for gene expression analyses, with its data normalization step being crucial for ensuring accuracy. Currently, the most widely used normalization method is based on the use of reference genes, assumed to be stably expressed across all experimental conditions. In the present study, we show that finding a stable combination of genes, regardless of their individual stability, outperforms standard reference genes for RT-qPCR data normalization. A stable combination of genes consists of a fixed number of genes whose individual expression balance each other all along experimental conditions of interest. Moreover, the present study shows that such an optimal combination of genes can be found using a comprehensive database of RNA-Seq data. Indeed, assuming that such a comprehensive database contains accurate gene expression profiles, we can extract in silico, by the way of the mathematical variance calculation, a stable combination of genes that reflects in vivo stability. As a case study, this new method was developed using the tomato model plant, with corresponding RNA-Seq data from the TomExpress database. However, the method is potentially applicable to other organisms with available RNA-seq data. Our results demonstrate the superiority of the reported method over commonly used housekeeping genes or other stably expressed genes. We therefore recommend the use of our new method together with classic ones in order to always obtain the best reference genes for a given experimental design.

Homoeologs in Allopolyploids: Navigating Redundancy as Both an Evolutionary Opportunity and a Technical Challenge-A Transcriptomics Perspective

Allopolyploidy in plants involves the merging of two or more distinct parental genomes into a single nucleus, a significant evolutionary process in the plant kingdom. Transcriptomic analysis provides invaluable insights into allopolyploid plants by elucidating the fate of duplicated genes, revealing evolutionary novelties and uncovering their environmental adaptations. By examining gene expression profiles, scientists can discern how duplicated genes have evolved to acquire new functions or regulatory roles. This process often leads to the development of novel traits and adaptive strategies that allopolyploid plants leverage to thrive in diverse ecological niches. Understanding these molecular mechanisms not only enhances our appreciation of the genetic complexity underlying allopolyploidy but also underscores their importance in agriculture and ecosystem resilience. However, transcriptome profiling is challenging due to genomic redundancy, which is further complicated by the presence of multiple chromosomes sets and the variations among homoeologs and allelic genes. Prior to transcriptome analysis, sub-genome phasing and homoeology inference are essential for obtaining a comprehensive view of gene expression. This review aims to clarify the terminology in this field, identify the most challenging aspects of transcriptome analysis, explain their inherent difficulties, and suggest reliable analytic strategies. Furthermore, bulk RNA-seq is highlighted as a primary method for studying allopolyploid gene expression, focusing on critical steps like read mapping and normalization in differential gene expression analysis. This approach effectively captures gene expression from both parental genomes, facilitating a comprehensive analysis of their combined profiles. Its sensitivity in detecting low-abundance transcripts allows for subtle differences between parental genomes to be identified, crucial for understanding regulatory dynamics and gene expression balance in allopolyploids.

Metabolic model guided CRISPRi identifies a central role for phosphoglycerate mutase in Chlamydia trachomatis persistence

Upon nutrient starvation, Chlamydia trachomatis serovar L2 (CTL) shifts from its normal growth to a non-replicating form, termed persistence. It is unclear if persistence reflects an adaptive response or a lack thereof. To understand this, transcriptomics data were collected for CTL grown under nutrient-replete and nutrient-starved conditions. Applying K-means clustering on transcriptomics data revealed a global transcriptomic rewiring of CTL under stress conditions in the absence of any canonical global stress regulator. This is consistent with previous data that suggested that CTL's stress response is due to a lack of an adaptive response mechanism. To investigate the impact of this on CTL metabolism, we reconstructed a genome-scale metabolic model of CTL (iCTL278) and contextualized it with the collected transcriptomics data. Using the metabolic bottleneck analysis on contextualized iCTL278, we observed that phosphoglycerate mutase (pgm) regulates the entry of CTL to the persistence state. Our data indicate that pgm has the highest thermodynamics driving force and lowest enzymatic cost. Furthermore, CRISPRi-driven knockdown of pgm in the presence or absence of tryptophan revealed the importance of this gene in modulating persistence. Hence, this work, for the first time, introduces thermodynamics and enzyme cost as tools to gain a deeper understanding on CTL persistence.

IMPORTANCE

This study uses a metabolic model to investigate factors that contribute to the persistence of Chlamydia trachomatis serovar L2 (CTL) under tryptophan and iron starvation conditions. As CTL lacks many canonical transcriptional regulators, the model was used to assess two prevailing hypotheses on persistence-that the chlamydial response to nutrient starvation represents a passive response due to the lack of regulators or that it is an active response by the bacterium. K-means clustering of stress-induced transcriptomics data revealed striking evidence in favor of the lack of adaptive (i.e., a passive) response. To find the metabolic signature of this, metabolic modeling pin-pointed pgm as a potential regulator of persistence. Thermodynamic driving force, enzyme cost, and CRISPRi knockdown of pgm supported this finding. Overall, this work introduces thermodynamic driving force and enzyme cost as a tool to understand chlamydial persistence, demonstrating how systems biology-guided CRISPRi can unravel complex bacterial phenomena.

SlMYC2 promotes SlLBD40-mediated cell expansion in tomato fruit development

As a plant-specific transcription factor, lateral organ boundaries domain (LBD) protein was reported to regulate plant growth and stress response, but the functional research of subfamily II genes is limited. SlMYC2, a master regulator of Jasmonic acid response, has been found to exhibit high expression levels in fruit and has been implicated in the regulation of fruit ripening and resistance to Botrytis. However, its role in fruit expansion remains unknown. In this study, we present evidence that a subfamily II member of LBD, namely SlLBD40, collaborates with SlMYC2 in the regulation of fruit expansion. Overexpression of SlLBD40 significantly promoted fruit growth by promoting mesocarp cell expansion, while knockout of SlLBD40 showed the opposite result. Similarly, SlMYC2 knockout resulted in a significant decrease in cell expansion within the fruit. Genetic analysis indicated that SlLBD40-mediated cell expansion depends on the expression of SlMYC2. SlLBD40 bound to the promoter of SlEXPA5, an expansin gene, but did not activate its expression directly. While, the co-expression of SlMYC2 and SlLBD40 significantly stimulated the activation of SlEXPA5, leading to an increase in fruit size. SlLBD40 interacted with SlMYC2 and enhanced the stability and abundance of SlMYC2. Furthermore, SlMYC2 directly targeted and activated the expression of SlLBD40, which is essential for SlLBD40-mediated fruit expansion. In summary, our research elucidates the role of the interaction between SlLBD40 and SlMYC2 in promoting cell expansion in tomato fruits, thus providing novel insights into the molecular genetics underlying fruit growth.

A normalization method that controls for total RNA abundance affects the identification of differentially expressed genes, revealing bias toward morning-expressed responses

RNA-Sequencing is widely used to investigate changes in gene expression at the transcription level in plants. Most plant RNA-Seq analysis pipelines base the normalization approaches on the assumption that total transcript levels do not vary between samples. However, this assumption has not been demonstrated. In fact, many common experimental treatments and genetic alterations affect transcription efficiency or RNA stability, resulting in unequal transcript abundance. The addition of synthetic RNA controls is a simple correction that controls for variation in total mRNA levels. However, adding spike-ins appropriately is challenging with complex plant tissue, and carefully considering how they are added is essential to their successful use. We demonstrate that adding external RNA spike-ins as a normalization control produces differences in RNA-Seq analysis compared to traditional normalization methods, even between two times of day in untreated plants. We illustrate the use of RNA spike-ins with 3' RNA-Seq and present a normalization pipeline that accounts for differences in total transcriptional levels. We evaluate the effect of normalization methods on identifying differentially expressed genes in the context of identifying the effect of the time of day on gene expression and response to chilling stress in sorghum.

Ploidy-specific transcriptomes shed light on the heterogeneous identity and metabolism of developing tomato pericarp cells

Endoreduplication, during which cells increase their DNA content through successive rounds of full genome replication without cell division, is the major source of endopolyploidy in higher plants. Endoreduplication plays pivotal roles in plant growth and development and is associated with the activation of specific transcriptional programmes that are characteristic of each cell type, thereby defining their identity. In plants, endoreduplication is found in numerous organs and cell types, especially in agronomically valuable ones, such as the fleshy fruit (pericarp) of tomato presenting high ploidy levels. We used the tomato pericarp tissue as a model system to explore the transcriptomes associated with endoreduplication progression during fruit growth. We confirmed that expression globally scales with ploidy level and identified sets of differentially expressed genes presenting only developmental-specific, only ploidy-specific expression patterns or profiles resulting from an additive effect of ploidy and development. When comparing ploidy levels at a specific developmental stage, we found that non-endoreduplicated cells are defined by cell division state and cuticle synthesis while endoreduplicated cells are mainly defined by their metabolic activity changing rapidly over time. By combining this dataset with publicly available spatiotemporal pericarp expression data, we proposed a map describing the distribution of ploidy levels within the pericarp. These transcriptome-based predictions were validated by quantifying ploidy levels within the pericarp tissue. This in situ ploidy quantification revealed the dynamic progression of endoreduplication and its cell layer specificity during early fruit development. In summary, the study sheds light on the complex relationship between endoreduplication, cell differentiation and gene expression patterns in the tomato pericarp.

Functional consequences of somatic polyploidy in development

Polyploid cells contain multiple genome copies and arise in many animal tissues as a regulated part of development. However, polyploid cells can also arise due to cell division failure, DNA damage or tissue damage. Although polyploidization is crucial for the integrity and function of many tissues, the cellular and tissue-wide consequences of polyploidy can be very diverse. Nonetheless, many polyploid cell types and tissues share a remarkable similarity in function, providing important information about the possible contribution of polyploidy to cell and tissue function. Here, we review studies on polyploid cells in development, underlining parallel functions between different polyploid cell types, as well as differences between developmentally-programmed and stress-induced polyploidy.

Endoreduplication in plant organogenesis: a means to boost fruit growth

Endoreduplication is the major source of somatic endopolyploidy in higher plants, and leads to variation in cell ploidy levels due to iterative rounds of DNA synthesis in the absence of mitosis. Despite its ubiquitous occurrence in many plant organs, tissues, and cells, the physiological meaning of endoreduplication is not fully understood, although several roles during plant development have been proposed, mostly related to cell growth, differentiation, and specialization via transcriptional and metabolic reprogramming. Here, we review recent advances in our knowledge of the molecular mechanisms and cellular characteristics of endoreduplicated cells, and provide an overview of the multi-scale effects of endoreduplication on supporting growth in plant development. In addition, the effects of endoreduplication in fruit development are discussed, since it is highly prominent during fruit organogenesis where it acts as a morphogenetic factor supporting rapid fruit growth, as illustrated by case of the model fleshy fruit, tomato (Solanum lycopersicum).

Beyond Transcript Concentrations: Quantifying Polyploid Expression Responses per Biomass, per Genome, and per Cell with RNA-Seq

RNA-seq has been used extensively to study expression responses to polyploidy. Most current methods for normalizing RNA-seq data yield estimates of transcript concentrations (transcripts per transcriptome). The implicit assumption of these normalization methods is that transcriptome size is equivalent between the samples being compared such that transcript concentrations are equivalent to transcripts per cell. In recent years, however, evidence has mounted that transcriptome size can vary dramatically in response to a range of factors including polyploidy and that such variation is ubiquitous. Where such variation exists, transcript concentration is often a poor or even misleading proxy for expression responses at other biologically relevant scales (e.g., expression per cell). Thus, it is important that transcriptomic studies of polyploids move beyond simply comparing transcript concentrations if we are to gain a complete understanding of how genome multiplication affects gene expression. I discuss this issue in more detail and summarize a suite of approaches that can leverage RNA-seq to quantify expression responses per genome, per cell, and per unit of biomass.

Characterization of Transcriptome Dynamics during Early Fruit Development in Olive (Olea europaea L.)

In the olive (Olea europaea L.), an economically leading oil crop worldwide, fruit size and yield are determined by the early stages of fruit development. However, few detailed analyses of this stage of fruit development are available. This study offers an extensive characterization of the various processes involved in early olive fruit growth (cell division, cell cycle regulation, and cell expansion). For this, cytological, hormonal, and transcriptional changes characterizing the phases of early fruit development were analyzed in olive fruit of the cv. 'Picual'. First, the surface area and mitotic activity (by flow cytometry) of fruit cells were investigated during early olive fruit development, from 0 to 42 days post-anthesis (DPA). The results demonstrate that the cell division phase extends up to 21 DPA, during which the maximal proportion of 4C cells in olive fruits was reached at 14 DPA, indicating that intensive cell division was activated in olive fruits at that time. Subsequently, fruit cell expansion lasted as long as 3 weeks more before endocarp lignification. Finally, the molecular mechanisms controlling the early fruit development were investigated by analyzing the transcriptome of olive flowers at anthesis (fruit set) as well as olive fruits at 14 DPA (cell division phase) and at 28 DPA (cell expansion phase). Sequential induction of the cell cycle regulating genes is associated with the upregulation of genes involved in cell wall remodeling and ion fluxes, and with a shift in plant hormone metabolism and signaling genes during early olive fruit development. This occurs together with transcriptional activity of subtilisin-like protease proteins together with transcription factors potentially involved in early fruit growth signaling. This gene expression profile, together with hormonal regulators, offers new insights for understanding the processes that regulate cell division and expansion, and ultimately fruit yield and olive size.

Genome size and endoreplication in two pairs of cytogenetically contrasting species of Pulmonaria (Boraginaceae) in Central Europe

Genome size is species-specific feature and commonly constant in an organism. In various plants, DNA content in cell nucleus is commonly increased in process of endoreplication, cellular-specific multiplication of DNA content without mitosis. This leads to the endopolyploidy, the presence of multiplied chromosome sets in a subset of cells. The relationship of endopolyploidy to species-specific genome size is rarely analysed and is not fully understood. While negative correlation between genome size and endopolyploidy level is supposed, this is species- and lineage-specific. In the present study, we shed light on this topic, exploring both genome size and endoreplication-induced DNA content variation in two pairs of morphologically similar species of Pulmonaria, P. obscura-P. officinalis and P. mollis-P. murinii. We aim (i) to characterize genome size and chromosome numbers in these species using cytogenetic, root-tip squashing and flow cytometry (FCM) techniques; (ii) to investigate the degree of endopolyploidy in various plant organs, including the root, stem, leaf, calyx and corolla using FCM; and (iii) to comprehensively characterize and compare the level of endopolyploidy and DNA content in various organs of all four species in relation to species systematic relationships and genome size variation. We have confirmed the diploid-dysploid nature of chromosome complements, and divergent genome sizes for Pulmonaria species: P. murinii with 2n = 2x = 14, 2.31 pg/2C, P. obscura 2n = 2x = 14, 2.69 pg/2C, P. officinalis 2n = 2x = 16, 2.96 pg/2C and P. mollis 2n = 2x = 18, 3.18 pg/2C. Endopolyploidy varies between species and organs, and we have documented 4C-8C in all four organs and up to 32C (64C) endopolyploid nuclei in stems at least in some species. Two species with lower genome sizes tend to have higher endopolyploidy levels than their closest relatives. Endoreplication-generated tissue-specific mean DNA content is increased and more balanced among species in all four organs compared to genome size. Our results argue for the narrow relationship between genome size and endopolyploidy in the present plant group within the genus Pulmonaria, and endopolyploidization seems to play a compensatory developmental role in organs of related morphologically similar species.

Current status and trends in forest genomics

Forests are not only the most predominant of the Earth's terrestrial ecosystems, but are also the core supply for essential products for human use. However, global climate change and ongoing population explosion severely threatens the health of the forest ecosystem and aggravtes the deforestation and forest degradation. Forest genomics has great potential of increasing forest productivity and adaptation to the changing climate. In the last two decades, the field of forest genomics has advanced quickly owing to the advent of multiple high-throughput sequencing technologies, single cell RNA-seq, clustered regularly interspaced short palindromic repeats (CRISPR)-mediated genome editing, and spatial transcriptomes, as well as bioinformatics analysis technologies, which have led to the generation of multidimensional, multilayered, and spatiotemporal gene expression data. These technologies, together with basic technologies routinely used in plant biotechnology, enable us to tackle many important or unique issues in forest biology, and provide a panoramic view and an integrative elucidation of molecular regulatory mechanisms underlying phenotypic changes and variations. In this review, we recapitulated the advancement and current status of 12 research branches of forest genomics, and then provided future research directions and focuses for each area. Evidently, a shift from simple biotechnology-based research to advanced and integrative genomics research, and a setup for investigation and interpretation of many spatiotemporal development and differentiation issues in forest genomics have just begun to emerge.

Plant Cell Identity in the Era of Single-Cell Transcriptomics

High-throughput single-cell transcriptomic approaches have revolutionized our view of gene expression at the level of individual cells, providing new insights into their heterogeneity, identities, and functions. Recently, technical challenges to the application of single-cell transcriptomics to plants have been overcome, and many plant organs and tissues have now been subjected to analyses at single-cell resolution. In this review, we describe these studies and their impact on our understanding of the diversity, differentiation, and activities of plant cells. We particularly highlight their impact on plant cell identity, including unprecedented views of cell transitions and definitions of rare and novel cell types. We also point out current challenges and future opportunities for the application and analyses of single-cell transcriptomics in plants. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution

The scarcity of accessible sites that are dynamic or cell type-specific in plants may be due in part to tissue heterogeneity in bulk studies. To assess the effects of tissue heterogeneity, we apply single-cell ATAC-seq to Arabidopsis thaliana roots and identify thousands of differentially accessible sites, sufficient to resolve all major cell types of the root. We find that the entirety of a cell's regulatory landscape and its transcriptome independently capture cell type identity. We leverage this shared information on cell identity to integrate accessibility and transcriptome data to characterize developmental progression, endoreduplication and cell division. We further use the combined data to characterize cell type-specific motif enrichments of transcription factor families and link the expression of family members to changing accessibility at specific loci, resolving direct and indirect effects that shape expression. Our approach provides an analytical framework to infer the gene regulatory networks that execute plant development.

The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution

The scarcity of accessible sites that are dynamic or cell type-specific in plants may be due in part to tissue heterogeneity in bulk studies. To assess the effects of tissue heterogeneity, we apply single-cell ATAC-seq to Arabidopsis thaliana roots and identify thousands of differentially accessible sites, sufficient to resolve all major cell types of the root. We find that the entirety of a cell's regulatory landscape and its transcriptome independently capture cell type identity. We leverage this shared information on cell identity to integrate accessibility and transcriptome data to characterize developmental progression, endoreduplication and cell division. We further use the combined data to characterize cell type-specific motif enrichments of transcription factor families and link the expression of family members to changing accessibility at specific loci, resolving direct and indirect effects that shape expression. Our approach provides an analytical framework to infer the gene regulatory networks that execute plant development.

Single-nucleus RNA and ATAC sequencing reveals the impact of chromatin accessibility on gene expression in Arabidopsis roots at the single-cell level

Similar to other complex organisms, plants consist of diverse and specialized cell types. The gain of unique biological functions of these different cell types is the consequence of the establishment of cell-type-specific transcriptional programs. As a necessary step in gaining a deeper understanding of the regulatory mechanisms controlling plant gene expression, we report the use of single-nucleus RNA sequencing (sNucRNA-seq) and single-nucleus assay for transposase accessible chromatin sequencing (sNucATAC-seq) technologies on Arabidopsis roots. The comparison of our single-nucleus transcriptomes to the published protoplast transcriptomes validated the use of nuclei as biological entities to establish plant cell-type-specific transcriptomes. Furthermore, our sNucRNA-seq results uncovered the transcriptomes of additional cell subtypes not identified by single-cell RNA-seq. Similar to our transcriptomic approach, the sNucATAC-seq approach led to the distribution of the Arabidopsis nuclei into distinct clusters, suggesting the differential accessibility of chromatin between groups of cells according to their identity. To reveal the impact of chromatin accessibility on gene expression, we integrated sNucRNA-seq and sNucATAC-seq data and demonstrated that cell-type-specific marker genes display cell-type-specific patterns of chromatin accessibility. Our data suggest that the differential chromatin accessibility is a critical mechanism to regulate gene activity at the cell-type level.

Single-Cell Transcriptome Analysis in Plants: Advances and Challenges

The rapid and enthusiastic adoption of single-cell RNA sequencing (scRNA-seq) has demonstrated that this technology is far more than just another way to perform transcriptome analysis. It is not an exaggeration to say that the advent of scRNA-seq is revolutionizing the details of whole-transcriptome snapshots from a tissue to a cell. With this disruptive technology, it is now possible to mine heterogeneity between tissue types and within cells like never before. This enables more rapid identification of rare and novel cell types, simultaneous characterization of multiple different cell types and states, more accurate and integrated understanding of their roles in life processes, and more. However, we are only at the beginning of unlocking the full potential of scRNA-seq applications. This is particularly true for plant sciences, where single-cell transcriptome profiling is in its early stage and has many exciting challenges to overcome. In this review, we compare and evaluate recent pioneering studies using the Arabidopsis root model, which has established new paradigms for scRNA-seq studies in plants. We also explore several new and promising single-cell analysis tools that are available to those wishing to study plant development and physiology at unprecedented resolution and scale. In addition, we propose some future directions on the use of scRNA-seq technology to tackle some of the critical challenges in plant research and breeding.

Single-cell RNA-seq analysis reveals ploidy-dependent and cell-specific transcriptome changes in Arabidopsis female gametophytes

BACKGROUND

Polyploidy provides new genetic material that facilitates evolutionary novelty, species adaptation, and crop domestication. Polyploidy often leads to an increase in cell or organism size, which may affect transcript abundance or transcriptome size, but the relationship between polyploidy and transcriptome changes remains poorly understood. Plant cells often undergo endoreduplication, confounding the polyploid effect.

RESULTS

To mitigate these effects, we select female gametic cells that are developmentally stable and void of endoreduplication. Using single-cell RNA sequencing (scRNA-seq) in Arabidopsis thaliana tetraploid lines and isogenic diploids, we show that transcriptome abundance doubles in the egg cell and increases approximately 1.6-fold in the central cell, consistent with cell size changes. In the central cell of tetraploid plants, DEMETER (DME) is upregulated, which can activate PRC2 family members FIS2 and MEA, and may suppress the expression of other genes. Upregulation of cell size regulators in tetraploids, including TOR and OSR2, may increase the size of reproductive cells. In diploids, the order of transcriptome abundance is central cell, synergid cell, and egg cell, consistent with their cell size variation. Remarkably, we uncover new sets of female gametophytic cell-specific transcripts with predicted biological roles; the most abundant transcripts encode families of cysteine-rich peptides, implying roles in cell-cell recognition during double fertilization.

CONCLUSIONS

Transcriptome in single cells doubles in tetraploid plants compared to diploid, while the degree of change and relationship to the cell size depends on cell types. These scRNA-seq resources are free of cross-contamination and are uniquely valuable for advancing plant hybridization, reproductive biology, and polyploid genomics.

Gene Balance Predicts Transcriptional Responses Immediately Following Ploidy Change in Arabidopsis thaliana

The gene balance hypothesis postulates that there is selection on gene copy number (gene dosage) to preserve the stoichiometric balance among interacting proteins. This presupposes that gene product abundance is governed by gene dosage and that gene dosage responses are consistent for interacting genes in a dosage-balance-sensitive network or complex. Gene dosage responses, however, have rarely been quantified, and the available data suggest that they are highly variable. We sequenced the transcriptomes of two synthetic autopolyploid accessions of Arabidopsis (Arabidopsis thaliana) and their diploid progenitors, as well as one natural tetraploid and its synthetic diploid produced via haploid induction, to estimate transcriptome size and dosage responses immediately following ploidy change. Similar to what has been observed in previous studies, overall transcriptome size does not exhibit a simple doubling in response to genome doubling, and individual gene dosage responses are highly variable in all three accessions, indicating that expression is not strictly coupled with gene dosage. Nonetheless, putatively dosage balance-sensitive gene groups (Gene Ontology terms, metabolic networks, gene families, and predicted interacting proteins) exhibit smaller and more coordinated dosage responses than do putatively dosage-insensitive gene groups, suggesting that constraints on dosage balance operate immediately following whole-genome duplication and that duplicate gene retention patterns are shaped by selection to preserve dosage balance.

Endoreplication - a means to an end in cell growth and stress response

Endoreplication, also called endoreduplication or endopolyploidization, is a cell cycle variant in which the genome is re-replicated in the absence of mitosis causing cellular polyploidization. Despite the common occurrence of endoreplication in plants and the tremendous extent in specific tissues and cell types such as the endosperm, the underlying molecular regulation and the physiological consequences have only now started to be understood. Endoreplication is often associated with cell differentiation and withdrawal from mitotic cycles. Recent studies have underlined the importance of endoreplication as a stress response and we summarize here this progress with particular focus on future perspectives offered by the recent advances in genomics and biotechnology.

Retracted Article: Investigation of yellow horn (Xanthoceras sorbifolia Bunge) transcriptome in response to different abiotic stresses: a comparative RNA-Seq study

Yellow horn (Xanthoceras sorbifolia Bunge) is a well-known oil-rich seed shrub which can grow well in barren and arid environments in the northern part of China. Yellow horn has received worldwide attention because of its excellent economic and environmental value. However, because of its limited genetic data, little information can be found regarding the molecular defense mechanisms of yellow horn exposed to various abiotic stresses. In view of this, the current study aims to investigate the impact of different abiotic stresses (i.e. NaCl, ABA and low temperature) on the transcriptome of yellow horn using RNA-Seq. Based on the transcriptome sequencing data, approximately 27% to 45% of stress-responsive genes were found highly expressed after stress treatment for 24 h. In addition, these genes were found to be still expressed after stress treatment for 48 h. However, many additional genes were stress-regulated after 48 h treatment compared with the 24 h treatment. GO enrichment analysis revealed that the expression patterns of the stress-responsive, type-specific terms were generally down-regulated. Most shared GO terms were primarily involved in protein folding, unfolding protein binding, protein transport and protein modification. Further, transcription factors (TFs), such as ERFs, bHLH, GRAS and NAC, were found to be enriched only in the low temperature treatment group, particularly the ERF TFs families. These combined results suggested that yellow horn may have developed specific molecular defense systems against diverse abiotic stresses.

Organ-wide and ploidy-dependent regulation both contribute to cell-size determination: evidence from a computational model of tomato fruit

The development of a new organ is the result of coordinated events of cell division and expansion, in strong interaction with each other. This study presents a dynamic model of tomato fruit development that includes cell division, endoreduplication, and expansion processes. The model is used to investigate the potential interactions among these developmental processes within the context of the neo-cellular theory. In particular, different control schemes (either cell-autonomous or organ-controlled) are tested and compared to experimental data from two contrasting genotypes. The model shows that a pure cell-autonomous control fails to reproduce the observed cell-size distribution, and that an organ-wide control is required in order to get realistic cell-size variations. The model also supports the role of endoreduplication as an important determinant of the final cell size and suggests that a direct effect of endoreduplication on cell expansion is needed in order to obtain a significant correlation between size and ploidy, as observed in real data.

Role of a cotton endoreduplication-related gene, GaTOP6B, in response to drought stress

Cotton GaTOP6B is involved in cellular endoreduplication and a positive response to drought stress via promoting plant leaf and root growth. Drought is deemed as one of adverse conditions that could cause substantial reductions in crop yields worldwide. Since cotton exhibits a moderate-tolerant phenotype under water-deficit conditions, the plant could therefore be used to characterize potential new genes regulating drought tolerance in crop plants. In this work, GaTOP6B, encoding DNA topoisomerase VI subunit B, was identified in Asian cotton (Gossypium arboreum). Virus-induced gene silencing (VIGS) and overexpression (OE) were used to investigate the biological function of GaTOP6B in G. arboreum and Arabidopsis thaliana under drought stress. The GaTOP6B-silencing plants showed a reduced ploidy level, and displayed a compromised tolerance phenotype including lowered relative water content (RWC), decreased proline content and antioxidative enzyme activity, and an increased malondialdehyde (MDA) content under drought stress. GaTOP6B-overexpressing Arabidopsis lines, however, had increased ploidy levels, and were more tolerant to drought treatment, associated with improved RWC maintenance, higher proline accumulation, and reduced stomatal aperture under drought stress. Transcriptome analysis showed that genes involved in the processes like cell cycle, transcription and signal transduction, were substantially up-regulated in GaTOP6B-overexpressing Arabidopsis, promoting plant growth and development. More specifically, under drought stress, the genes involved in the biosynthesis of secondary metabolites such as phenylpropanoid, starch and sucrose were selectively enhanced to improve tolerance in plants. Taken together, the results demonstrated that GaTOP6B could coordinately regulate plant leaf and root growth via cellular endoreduplication, and positively respond to drought stress. Thus, GaTOP6B could be a competent candidate gene for improvement of drought tolerance in crop species.

Divergent gene expression levels between diploid and autotetraploid Tolmiea relative to the total transcriptome, the cell, and biomass

PREMISE OF THE STUDY

Studies of gene expression and polyploidy are typically restricted to characterizing differences in transcript concentration. Using diploid and autotetraploid Tolmiea, we present an integrated approach for cross-ploidy comparisons that account for differences in transcriptome size and cell density and make multiple comparisons of transcript abundance.

METHODS

We use RNA spike-in standards in concert with cell size and density to identify and correct for differences in transcriptome size and compare levels of gene expression across multiple scales: per transcriptome, per cell, and per biomass.

KEY RESULTS

In total, ~17% of all loci were identified as differentially expressed (DEGs) between the diploid and autopolyploid species. The per-transcriptome normalization, the method researchers typically use, captured the fewest DEGs (58% of total DEGs) and failed to detect any DEGs not found by the alternative normalizations. When transcript abundance was normalized per biomass and per cell, ~66% and ~82% of the total DEGs were recovered, respectively. The discrepancy between per-transcriptome and per-cell recovery of DEGs occurs because per-transcriptome normalizations are concentration-based and therefore blind to differences in transcriptome size.

CONCLUSIONS

While each normalization enables valid comparisons at biologically relevant scales, a holistic comparison of multiple normalizations provides additional explanatory power not available from any single approach. Notably, autotetraploid loci tend to conserve diploid-like transcript abundance per biomass through increased gene expression per cell, and these loci are enriched for photosynthesis-related functions.

Meeting report - INDEPTH kick-off meeting

The precise location of chromatin domains within the cell nucleus has seen growing recognition in the past decade as an additional mechanism of controlling gene expression in both plants and animals (Dekker et al., 2017). Consequently, international efforts are devoted to understanding the organising principle of this organelle in plants, and notably the nature and the role of functional compartments on gene expression (Graumann et al., 2013; Sotelo-Silveira et al., 2018). The European cooperation 'Impact of Nuclear Domains on Gene Expression and Plant Traits' (INDEPTH) brings together molecular cell biologists, plant physiologists, bioinformaticians, image analysts and computer scientists. They aim to address the question of how nuclear architecture, chromatin organisation and gene expression are connected in plants, particularly in relation to traits of interest such as biomass, reproduction and resistance to pathogens (https://www.brookes.ac.uk/indepth/). The kick-off meeting of the INDEPTH consortium took place in Clermont-Ferrand, France, on 12-14th March 2018, where more than 80 researchers set the agenda for the coming four years of research and collaboration.