Comparison of normalization methods for differential gene expression analysis in RNA-Seq experiments: A matter of relative size of studied transcriptomes

In recent years, RNA-Seq technologies became a powerful tool for transcriptome studies. However, computational methods dedicated to the analysis of high-throughput sequencing data are yet to be standardized. In particular, it is known that the choice of a normalization procedure leads to a great variability in results of differential gene expression analysis. The present study compares the most widespread normalization procedures and proposes a novel one aiming at removing an inherent bias of studied transcriptomes related to their relative size. Comparisons of the normalization procedures are performed on real and simulated data sets. Real RNA-Seq data sets analyses, performed with all the different normalization methods, show that only 50% of significantly differentially expressed genes are common. This result highlights the influence of the normalization step on the differential expression analysis. Real and simulated data sets analyses give similar results showing 3 different groups of procedures having the same behavior. The group including the novel method named “Median Ratio Normalization” (MR N) gives the lower number of false discoveries. Within this group the MR N method is less sensitive to the modification of parameters related to the relative size of transcriptomes such as the number of down- and upregulated genes and the gene expression levels. The newly proposed MR N method efficiently deals with intrinsic bias resulting from relative size of studied transcriptomes. Validation with real and simulated data sets confirmed that MR N is more consistent and robust than existing methods

Recent Developments on the Role of Ethylene in the Ripening of Climacteric Fruit

It has long been recognised that ethylene plays a major role in the ripening process of climacteric fruit. A more thorough analysis, however, has revealed that a number of biochemical and molecular processes associated with climacteric fruit ripening are ethylene-independent. One of the crucial steps of the onset of ripening is the induction of autocatalytic ethylene production. In ethylene-suppressed melons, ACC synthase activity is induced at the same time as in control melons, indicating that ACC biosynthesis during the early stages of ripening seems to be a developmentally-regulated (ethylene-independent) process. The various ripening events exhibit differential sensitivity to ethylene. For instance, the threshold level for degreening of the rind is 1ppm, while 2.5 ppm are required to trigger some components of the softening process. The saturating level of ethylene producing maximum effects is less than 5 ppm, which is by far lower than the internal ethylene concentrations found in the fruit at the climacteric peak (over 100 ppm). In many fruit chilling temperatures hasten ethylene production and ripening and in some late season pear varieties, exposure to chilling temperatures is even absolutely required for the attainment of the capacity to synthesize autocatalytic ethylene. This is correlated with the stimulation of expression of ACC oxidase and of members of the ACC synthase gene family. Ethylene operates via a perception and transduction pathway to induce the expression of genes responsible for the biochemical and physiological changes observed during ripening. However, only a few genes induced via the ethylene transduction pathway have been described so far. We have used a differential display method to isolate novel ethylene-reponsive (ER) cDNA clones of tomato that potentially play a role in propagating the ethylene response and in regulating fruit ripening. Collectively, these data permit a general scheme of the molecular mechanisms of fruit ripening to be proposed

TomExpress, a unified tomato RNA-Seq platform for visualization of expression data, clustering and correlation networks

The TomExpress platform was developed to provide the tomato research community with a browser and integrated web tools for public RNA-Seq data visualization and data mining. To avoid major biases that can result from the use of different mapping and statistical processing methods, RNA-Seq raw sequence data available in public databases were mapped de novo on a unique tomato reference genome sequence and post-processed using the same pipeline with accurate parameters. Following the calculation of the number of counts per gene in each RNA-Seq sample, a communal global normalization method was applied to all expression values. This unifies the whole set of expression data and makes them comparable. A database was designed where each expression value is associated with corresponding experimental annotations. Sample details were manually curated to be easily understandable by biologists. To make the data easily searchable, a user-friendly web interface was developed that provides versatile data mining web tools via on-the-fly generation of output graphics, such as expression bar plots, comprehensive in planta representations and heatmaps of hierarchically clustered expression data. In addition, it allows for the identification of co-expressed genes and the visualization of correlation networks of co-regulated gene groups. TomExpress provides one of the most complete free resources of publicly available tomato RNA-Seq data, and allows for the immediate interrogation of transcriptional programs that regulate vegetative and reproductive development in tomato under diverse conditions. The design of the pipeline developed in this project enables easy updating of the database with newly published RNA-Seq data, thereby allowing for continuous enrichment of the resource

Regulatory Features Underlying Pollination-Dependent and -Independent Tomato Fruit Set Revealed by Transcript and Primary Metabolite Profiling

Indole Acetic Acid 9 (IAA9) is a negative auxin response regulator belonging to the Aux/IAA transcription factor gene family whose downregulation triggers fruit set before pollination, thus giving rise to parthenocarpy. In situ hybridization experiments revealed that a tissue-specific gradient of IAA9 expression is established during flower development, the release of which upon pollination triggers the initiation of fruit development. Comparative transcriptome and targeted metabolome analysis uncovered important features of the molecular events underlying pollination-induced and pollination-independent fruit set. Comprehensive transcriptomic profiling identified a high number of genescommonto both types of fruit set,amongwhich only a small subset are dependent on IAA9 regulation. The fine-tuning of Aux/IAA and ARF genes and the downregulation of TAG1 and TAGL6 MADS box genes are instrumental in triggering the fruit set program. Auxin and ethylene emerged as the most active signaling hormones involved in the flower-to-fruit transition. However, while these hormones affected only a small number of transcriptional events, dramatic shifts were observed at the metabolic and developmental levels. The activation of photosynthesis and sucrose metabolism-related genes is an integral regulatory component of fruit set process. The combined results allow a far greater comprehension of the regulatory and metabolic events controlling early fruit development both in the presence and absence of pollination/fertilization

Comprehensive Profiling of Ethylene Response Factor Expression Identifies Ripening-Associated ERF Genes and Their Link to Key Regulators of Fruit Ripening in Tomato

Our knowledge of the factors mediating ethylene-dependent ripening of climacteric fruit remains limited. The transcription of ethylene-regulated genes is mediated by ethylene response factors (ERFs), but mutants providing information on the specific role of the ERFs in fruit ripening are still lacking, likely due to functional redundancy among this large multigene family of transcription factors. We present here a comprehensive expression profiling of tomato (Solanum lycopersicum) ERFs in wild-type and tomato ripening-impaired tomato mutants (Never-ripe [Nr], ripening-inhibitor [rin], and non-ripening [nor]), indicating that out of the 77 ERFs present in the tomato genome, 27 show enhanced expression at the onset of ripening while 28 display a ripeningassociated decrease in expression, suggesting that different ERFs may have contrasting roles in fruit ripening. Among the 19 ERFs exhibiting the most consistent up-regulation during ripening, the expression of 11 ERFs is strongly down-regulated in rin, nor, and Nr tomato ripening mutants, while only three are consistently up-regulated. Members of subclass E, SlERF.E1, SlERF.E2, and SlERF.E4, show dramatic down-regulation in the ripening mutants, suggesting that their expression might be instrumental in fruit ripening. This study illustrates the high complexity of the regulatory network connecting RIN and ERFs and identifies subclass E members as the most active ERFs in ethylene- and RIN/NOR-dependent ripening

Interaction of two MADS-box genes leads to growth phenotype divergence of all-flesh type of tomatoes

[EN] All-flesh tomato cultivars are devoid of locular gel and exhibit enhanced firmness and improved postharvest storage. Here, we show that SlMBP3 is a master regulator of locular tissue in tomato fruit and that a deletion at the gene locus underpins the All-flesh trait. Intriguingly, All-flesh varieties lack the deleterious phenotypes reported previously for SlMBP3 under-expressing lines and which preclude any potential commercial use. We resolve the causal factor for this phenotypic divergence through the discovery of a natural mutation at the SlAGL11 locus, a close homolog of SlMBP3. Misexpressing SlMBP3 impairs locular gel formation through massive transcriptomic reprogramming at initial phases of fruit development. SlMBP3 influences locule gel formation by controlling cell cycle and cell expansion genes, indicating that important components of fruit softening are determined at early pre-ripening stages. Our findings define potential breeding targets for improved texture in tomato and possibly other fleshy fruits. The all-flesh type of tomato fruits is caused by mutation of the MBP3 gene, however, knocking down MBP3 in certain genotypes also affect plant and fruit development. Here, the authors show that a natural mutation of AGL11, a close homolog of MBP3, is responsible for the phenotypic divergence.The authors are grateful to L. Lemonnier and D. Saint-Martin for transformation and cultivation of tomato plants and GeT-PlaGe core facility (INRAe Toulouse) for ChIP deep sequencing. The authors also want to thank Dr. Christian Chevalier (INRAE et Univsersite de Bordeaux) for helping in analyzing genes related to cell cycle, cell division, and endoreduplication in tomato. This research was supported by the EU H2020 TomGEM 679796 and HARNESSTOM 101000716 projects.Huang, B.; Hu, G.; Wang, K.; Frasse, P.; Maza, E.; Djari, A.; Deng, W.... (2021). Interaction of two MADS-box genes leads to growth phenotype divergence of all-flesh type of tomatoes. Nature Communications. 12(1):1-14. https://doi.org/10.1038/s41467-021-27117-711412

Overexpression of the class D MADS-box gene Sl-AGL11 impacts fleshy tissue differentiation and structure in tomato fruits

MADS-box transcription factors are key elements of the genetic networks controlling flower and fruit development. Among these, the class D clade gathers AGAMOUS-like genes which are involved in seed, ovule, and funiculus development. The tomato genome comprises two class D genes, Sl-AGL11 and Sl-MBP3 , both displaying high expression levels in seeds and in central tissues of young fruits. The potential effects of Sl-AGL11 on fruit development were addressed through RNAi silencing and ectopic expression strategies. Sl-AGL11-down-regulated tomato lines failed to show obvious phenotypes except a slight reduction in seed size. In contrast, Sl-AGL11 overexpression triggered dramatic modifications of flower and fruit structure that include: the conversion of sepals into fleshy organs undergoing ethylene-dependent ripening, a placenta hypertrophy to the detriment of locular space, starch and sugar accumulation, and an extreme softening that occurs well before the onset of ripening. RNA-Seq transcriptomic profiling high-lighted substantial metabolic reprogramming occurring in sepals and fruits, with major impacts on cell wall-related genes. While several Sl-AGL11-related phenotypes are reminiscent of class C MADS-box genes (TAG1 and TAGL1), the modifications observed on the placenta and cell wall and the Sl-AGL11 expression pattern suggest an action of this class D MADS-box factor on early fleshy fruit development

High resolution synteny maps allowing direct comparisons between the coffee and tomato genomes

Tomato (Solanum lycopersicum) and coffee (Coffea canephora) belong to the sister families Solanaceae and Rubiaceae, respectively. We report herein the mapping of a common set of 257 Conserved Ortholog Set II genes in the genomes of both species. The mapped markers are well distributed across both genomes allowing the first syntenic comparison between species from these two families. The majority (75%) of the synteny blocks are short (<4 cM); however, some extend up to 50 cM. In an effort to further characterize the synteny between these two genomes, we took advantage of the available sequence for the tomato genome to show that tomato chromosome 7 is syntenic to half of the two coffee linkage groups E and F with the putative break point in tomato localized to the boundary of the heterochromatin and euchromatin on the long arm. In addition to the new insight on genome conservation and evolution between the plant families Solanaceae and Rubiaceae, the comparative maps presented herein provide a translational tool by which coffee researchers may take benefit of DNA sequence and genetic information from tomato and vice versa. It is thus expected that these comparative genome information will help to facilitate and expedite genetic and genomic research in coffee

Sl-IAA3, a tomato Aux/IAA at the crossroads of auxin and ethylene signalling involved in differential growth

Whereas the interplay of multiple hormones is essential for most plant developmental processes, the key integrating molecular players remain largely undiscovered or uncharacterized. It is shown here that a member of the tomato auxin/indole-3-acetic acid (Aux/IAA) gene family, Sl-IAA3, intersects the auxin and ethylene signal transduction pathways. Aux/IAA genes encode short-lived transcriptional regulators central to the control of auxin responses. Their functions have been defined primarily by dominant, gain-of-function mutant alleles in Arabidopsis. The Sl-IAA3 gene encodes a nuclear-targeted protein that can repress transcription from auxin-responsive promoters. Sl-IAA3 expression is auxin and ethylene dependent, is regulated on a tight tissue-specific basis, and is associated with tissues undergoing differential growth such as in epinastic petioles and apical hook. Antisense down-regulation of Sl-IAA3 results in auxin and ethylene-related phenotypes, including altered apical dominance, lower auxin sensitivity, exaggerated apical hook curvature in the dark and reduced petiole epinasty in the light. The results provide novel insights into the roles of Aux/IAAs and position the Sl-IAA3 protein at the crossroads of auxin and ethylene signalling in tomato

Characterization of the Tomato ARF Gene Family Uncovers a Multi-Levels Post-Transcriptional Regulation Including Alternative Splicing

Background: The phytohormone auxin is involved in a wide range of developmental processes and auxin signaling is known to modulate the expression of target genes via two types of transcriptional regulators, namely, Aux/IAA and Auxin Response Factors (ARF). ARFs play a major role in transcriptional activation or repression through direct binding to the promoter of auxin-responsive genes. The present study aims at gaining better insight on distinctive structural and functional features among ARF proteins. Results: Building on the most updated tomato (Solanum lycopersicon) reference genome sequence, a comprehensive set of ARF genes was identified, extending the total number of family members to 22. Upon correction of structural annotation inconsistencies, renaming the tomato ARF family members provided a consensus nomenclature for all ARF genes across plant species. In silico search predicted the presence of putative target site for small interfering RNAs within twelve Sl-ARFs while sequence analysis of the 59-leader sequences revealed the presence of potential small uORF regulatory elements. Functional characterization carried out by transactivation assay partitioned tomato ARFs into repressors and activators of auxin-dependent gene transcription. Expression studies identified tomato ARFs potentially involved in the fruit set process. Genome-wide expression profiling using RNA-seq revealed that at least one third of the gene family members display alternative splicing mode of regulation during the flower to fruit transition. Moreover, the regulation of several tomato ARF genes by both ethylene and auxin, suggests their potential contribution to the convergence mechanism between the signaling pathways of these two hormones. Conclusion: All together, the data bring new insight on the complexity of the expression control of Sl-ARF genes at the transcriptional and post-transcriptional levels supporting the hypothesis that these transcriptional mediators might represent one of the main components that enable auxin to regulate a wide range of physiological processes in a highly specific and coordinated manner