Comparative Analysis of Developmental Transcriptome Maps of Arabidopsis thaliana and Solanum lycopersicum

A comparative analysis of stably expressed genes across diverse angiosperms exposes flexibility in underlying promoter architecture

Promoters regulate both the amplitude and pattern of gene expression-key factors needed for optimization of many synthetic biology applications. Previous work in Arabidopsis found that promoters that contain a TATA-box element tend to be expressed only under specific conditions or in particular tissues, while promoters which lack any known promoter elements, thus designated as Coreless, tend to be expressed more ubiquitously. To test whether this trend represents a conserved promoter design rule, we identified stably expressed genes across multiple angiosperm species using publicly available RNA-seq data. Comparisons between core promoter architectures and gene expression stability revealed differences in core promoter usage in monocots and eudicots. Furthermore, when tracing the evolution of a given promoter across species, we found that core promoter type was not a strong predictor of expression stability. Our analysis suggests that core promoter types are correlative rather than causative in promoter expression patterns and highlights the challenges in finding or building constitutive promoters that will work across diverse plant species.

Toward kingdom-wide analyses of gene expression

Gene expression data for Archaeplastida are accumulating exponentially, with more than 300 000 RNA-sequencing (RNA-seq) experiments available for hundreds of species. The gene expression data stem from thousands of experiments that capture gene expression in various organs, tissues, cell types, (a)biotic perturbations, and genotypes. Advances in software tools make it possible to process all these data in a matter of weeks on modern office computers, giving us the possibility to study gene expression in a kingdom-wide manner for the first time. We discuss how the expression data can be accessed and processed and outline analyses that take advantage of cross-species analyses, allowing us to generate powerful and robust hypotheses about gene function and evolution.

Interspecific comparison of gene expression profiles using machine learning

Interspecific gene comparisons are the keystones for many areas of biological research and are especially important for the translation of knowledge from model organisms to economically important species. Currently they are hampered by the low resolution of methods based on sequence analysis and by the complex evolutionary history of eukaryotic genes. This is especially critical for plants, whose genomes are shaped by multiple whole genome duplications and subsequent gene loss. This requires the development of new methods for comparing the functions of genes in different species. Here, we report ISEEML (Interspecific Similarity of Expression Evaluated using Machine Learning)-a novel machine learning-based algorithm for interspecific gene classification. In contrast to previous studies focused on sequence similarity, our algorithm focuses on functional similarity inferred from the comparison of gene expression profiles. We propose novel metrics for expression pattern similarity-expression score (ES)-that is suitable for species with differing morphologies. As a proof of concept, we compare detailed transcriptome maps of Arabidopsis thaliana, the model species, Zea mays (maize) and Fagopyrum esculentum (common buckwheat), which are species that represent distant clades within flowering plants. The classifier resulted in an AUC of 0.91; under the ES threshold of 0.5, the specificity was 94%, and sensitivity was 72%.

Function Analysis of the ERF and DREB Subfamilies in Tomato Fruit Development and Ripening

APETALA2/ethylene responsive factors (AP2/ERF) are unique regulators in the plant kingdom and are involved in the whole life activity processes such as development, ripening, and biotic and abiotic stresses. In tomato (Solanum lycopersicum), there are 140 AP2/ERF genes; however, their functionality remains poorly understood. In this work, the 14th and 19th amino acid differences in the AP2 domain were used to distinguish DREB and ERF subfamily members. Even when the AP2 domain of 68 ERF proteins from 20 plant species and motifs in tomato DREB and ERF proteins were compared, the binding ability of DREB and ERF proteins with DRE/CRT and/or GCC boxes remained unknown. During fruit development and ripening, the expressions of 13 DREB and 19 ERF subfamily genes showed some regular changes, and the promoters of most genes had ARF, DRE/CRT, and/or GCC boxes. This suggests that these genes directly or indirectly respond to IAA and/or ethylene (ET) signals during fruit development and ripening. Moreover, some of these may feedback regulate IAA or ET biosynthesis. In addition, 16 EAR motif-containing ERF genes in tomato were expressed in many organs and their total transcripts per million (TPM) values exceeded those of other ERF genes in most organs. To determine whether the EAR motif in EAR motif-containing ERF proteins has repression function, their EAR motifs were retained or deleted in a yeast one-hybrid (YIH) assay. The results indicate that most of EAR motif-containing ERF proteins lost repression activity after deleting the EAR motif. Moreover, some of these were expressed during ripening. Thus, these EAR motif-containing ERF proteins play vital roles in balancing the regulatory functions of other ERF proteins by completing the DRE/CRT and/or GCC box sites of target genes to ensure normal growth and development in tomato.

Transcriptome atlas of Phalaenopsis equestris

The vast diversity of Orchidaceae together with sophisticated adaptations to pollinators and other unique features make this family an attractive model for evolutionary and functional studies. The sequenced genome of Phalaenopsis equestris facilitates Orchidaceae research. Here, we present an RNA-seq-based transcriptome map of P. equestris that covers 19 organs of the plant, including leaves, roots, floral organs and the shoot apical meristem. We demonstrated the high quality of the data and showed the similarity of the P. equestris transcriptome map with the gene expression atlases of other plants. The transcriptome map can be easily accessed through our database Transcriptome Variation Analysis (TraVA) for visualizing gene expression profiles. As an example of the application, we analyzed the expression of Phalaenopsis “orphan” genes–those that do not have recognizable similarity with the genes of other plants. We found that approximately half of these genes were not expressed; the ones that were expressed were predominantly expressed in reproductive structures.

Transcriptome profiling and differential gene expression analysis provides insights into Lr24-based resistance in wheat against Puccinia triticina

Leaf rust caused by Puccinia triticina is an important disease of wheat and Lr24 gene confers resistance to all known pathotypes of P. triticina in India. Transcripts associated with the Lr24 mediated resistance were identified through transcriptome sequencing and further expression analysis of differentially regulated genes was performed using qPCR technique. De novo transcriptome assembly showed 66,415 and 68,688 transcripts in resistant and susceptible genotypes, respectively. The study revealed that 5873 genes unique to resistant; 6782 genes unique to susceptible, while 10,841 genes were common to both. Gene Ontology distribution statistics showed 1030 and 1068 CDS in biological processes; 1234 and 1326 CDS in cellular processes; 1321 and 1352 CDS in molecular functions, respectively. A total of 659 genes were found to be differentially expressed, of which 349 were upregulated and 310 were downregulated in resistant genotype. Pathway analysis of transcripts appeared in resistant genotype revealed that 279 transcripts had homology with genes involved in signal transduction, 18 transcripts in membrane transport, one transcript in signaling molecules. Real-time PCR study showed that most of the up-regulated defense related genes expressed in early hours indicating that a cascade of defense starts early in Lr24 mediated resistance, which successfully inhibited pathogen establishment.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02972-9.

High-Resolution Transcriptome Atlas and Improved Genome Assembly of Common Buckwheat, Fagopyrum esculentum

Common buckwheat (Fagopyrum esculentum) is an important non-cereal grain crop and a prospective component of functional food. Despite this, the genomic resources for this species and for the whole family Polygonaceae, to which it belongs, are scarce. Here, we report the assembly of the buckwheat genome using long-read technology and a high-resolution expression atlas including 46 organs and developmental stages. We found that the buckwheat genome has an extremely high content of transposable elements, including several classes of recently (0.5-1 Mya) multiplied TEs ("transposon burst") and gradually accumulated TEs. The difference in TE content is a major factor contributing to the three-fold increase in the genome size of F. esculentum compared with its sister species F. tataricum. Moreover, we detected the differences in TE content between the wild ancestral subspecies F. esculentum ssp. ancestrale and buckwheat cultivars, suggesting that TE activity accompanied buckwheat domestication. Expression profiling allowed us to test a hypothesis about the genetic control of petaloidy of tepals in buckwheat. We showed that it is not mediated by B-class gene activity, in contrast to the prediction from the ABC model. Based on a survey of expression profiles and phylogenetic analysis, we identified the MYB family transcription factor gene tr_18111 as a potential candidate for the determination of conical cells in buckwheat petaloid tepals. The information on expression patterns has been integrated into the publicly available database TraVA: http://travadb.org/browse/Species=Fesc/. The improved genome assembly and transcriptomic resources will enable research on buckwheat, including practical applications.

Histidine-Rich Defensins from the Solanaceae and Brasicaceae Are Antifungal and Metal Binding Proteins

Plant defensins are best known for their antifungal activity and contribution to the plant immune system. The defining feature of plant defensins is their three-dimensional structure known as the cysteine stabilized alpha-beta motif. This protein fold is remarkably tolerant to sequence variation with only the eight cysteines that contribute to the stabilizing disulfide bonds absolutely conserved across the family. Mature defensins are typically 46–50 amino acids in length and are enriched in lysine and/or arginine residues. Examination of a database of approximately 1200 defensin sequences revealed a subset of defensin sequences that were extended in length and were enriched in histidine residues leading to their classification as histidine-rich defensins (HRDs). Using these initial HRD sequences as a query, a search of the available sequence databases identified over 750 HRDs in solanaceous plants and 20 in brassicas. Histidine residues are known to contribute to metal binding functions in proteins leading to the hypothesis that HRDs would have metal binding properties. A selection of the HRD sequences were recombinantly expressed and purified and their antifungal and metal binding activity was characterized. Of the four HRDs that were successfully expressed all displayed some level of metal binding and two of four had antifungal activity. Structural characterization of the other HRDs identified a novel pattern of disulfide linkages in one of the HRDs that is predicted to also occur in HRDs with similar cysteine spacing. Metal binding by HRDs represents a specialization of the plant defensin fold outside of antifungal activity.

An Update to the TraVA Database: Time Series of Capsella bursa-pastoris Shoot Apical Meristems during Transition to Flowering

Transition to flowering is a crucial part of plant life directly affecting the fitness of a plant. Time series of transcriptomes is a useful tool for the investigation of process dynamics and can be used for the identification of novel genes and gene networks involved in the process. We present a detailed time series of polyploid Capsella bursa-pastoris shoot apical meristems created with RNA-seq. The time series covers transition to flowering and can be used for thorough analysis of the process. To make the data easy to access, we uploaded them in our database Transcriptome Variation Analysis (TraVA), which provides a convenient depiction of the gene expression profiles, the differential expression analysis between the homeologs and quick data extraction.

RNA sequencing by direct tagmentation of RNA/DNA hybrids

RNA sequencing is widely used to measure gene expression in biomedical research; therefore, improvements in the simplicity and accuracy of the technology are desirable. All existing RNA sequencing methods rely on the conversion of RNA into double-stranded DNA through reverse transcription followed by second-strand synthesis. The latter step requires additional enzymes and purification, and introduces sequence-dependent bias. Here, we show that Tn5 transposase, which randomly binds and cuts double-stranded DNA, can directly fragment and prime the RNA/DNA heteroduplexes generated by reverse transcription. The primed fragments are then subject to PCR amplification. This provides an approach for simple and accurate RNA characterization and quantification.

Transcriptome profiling by RNA sequencing (RNA-seq) has been widely used to characterize cellular status, but it relies on second-strand complementary DNA (cDNA) synthesis to generate initial material for library preparation. Here we use bacterial transposase Tn5, which has been increasingly used in various high-throughput DNA analyses, to construct RNA-seq libraries without second-strand synthesis. We show that Tn5 transposome can randomly bind RNA/DNA heteroduplexes and add sequencing adapters onto RNA directly after reverse transcription. This method, Sequencing HEteRo RNA-DNA-hYbrid (SHERRY), is versatile and scalable. SHERRY accepts a wide range of starting materials, from bulk RNA to single cells. SHERRY offers a greatly simplified protocol and produces results with higher reproducibility and GC uniformity compared with prevailing RNA-seq methods.

Expression of galactinol synthase from Ammopiptanthus nanus in tomato improves tolerance to cold stress

Soluble carbohydrates not only directly affect plant growth and development but also act as signal molecules in processes that enhance tolerance to cold stress. Raffinose family oligosaccharides (RFOs) are an example and play an important role in abiotic stress tolerance. This study aimed to determine whether galactinol, a key limiting factor in RFO biosynthesis, functions as a signal molecule in triggering cold tolerance. Exposure to low temperatures induces the expression of galactinol synthase (AnGolS1) in Ammopiptanthus nanus, a desert plant that survives temperatures between -30 °C to 47 °C. AnGolS1 has a greater catalytic activity than tomato galactinol synthase (SlGolS2). Moreover, SlGolS2 is expressed only at low levels. Expression of AnGolS1 in tomato enhanced cold tolerance and led to changes in the sugar composition of the seeds and seedlings. AnGolS1 transgenic tomato lines exhibited an enhanced capacity for ethylene (ET) signaling. The application of galactinol abolished the repression of the ET signaling pathway by 1-methylcyclopropene during seed germination. In addition, the expression of ERF transcription factors was increased. Galactinol may therefore act as a signal molecule affecting the ET pathway.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

Transcriptome Profiling Provides Insight into the Genes in Carotenoid Biosynthesis during the Mesocarp and Seed Developmental Stages of Avocado (Persea americana)

Avocado (Persea americana Mill.) is an economically important crop because of its high nutritional value. However, the absence of a sequenced avocado reference genome has hindered investigations of secondary metabolism. For next-generation high-throughput transcriptome sequencing, we obtained 365,615,152 and 348,623,402 clean reads as well as 109.13 and 104.10 Gb of sequencing data for avocado mesocarp and seed, respectively, during five developmental stages. High-quality reads were assembled into 100,837 unigenes with an average length of 847.40 bp (N50 = 1725 bp). Additionally, 16,903 differentially expressed genes (DEGs) were detected, 17 of which were related to carotenoid biosynthesis. The expression levels of most of these 17 DEGs were higher in the mesocarp than in the seed during five developmental stages. In this study, the avocado mesocarp and seed transcriptome were also sequenced using single-molecule long-read sequencing to acquired 25.79 and 17.67 Gb clean data, respectively. We identified 233,014 and 238,219 consensus isoforms in avocado mesocarp and seed, respectively. Furthermore, 104 and 59 isoforms were found to correspond to the putative 11 carotenoid biosynthetic-related genes in the avocado mesocarp and seed, respectively. The isoform numbers of 10 out of the putative 11 genes involved in the carotenoid biosynthetic pathway were higher in the mesocarp than those in the seed. Besides, alpha- and beta-carotene contents in the avocado mesocarp and seed during five developmental stages were also measured, and they were higher in the mesocarp than in the seed, which validated the results of transcriptome profiling. Gene expression changes and the associated variations in gene dosage could influence carotenoid biosynthesis. These results will help to further elucidate carotenoid biosynthesis in avocado.

Transcriptome Sequencing of Different Avocado Ecotypes: de novo Transcriptome Assembly, Annotation, Identification and Validation of EST-SSR Markers

Avocado (Persea americana Mill.) could be considered as an important tropical and subtropical woody oil crop with high economic and nutritional value. Despite the importance of this species, genomic information is currently unavailable for avocado and closely related congeners. In this study, we generated more than 216 million clean reads from different avocado ecotypes using Illumina HiSeq high-throughput sequencing technology. The high-quality reads were assembled into 154,310 unigenes with an average length of 922 bp. A total of 55,558 simple sequence repeat (SSR) loci detected among the 43,270 SSR-containing unigene sequences were used to develop 74,580 expressed sequence tag (EST)-SSR markers. From these markers, a subset of 100 EST-SSR markers was randomly chosen to identify polymorphic EST-SSR markers in 28 avocado accessions. Sixteen EST-SSR markers with moderate to high polymorphism levels were detected, with polymorphism information contents ranging from 0.33 to 0.84 and averaging 0.63. These 16 polymorphic EST-SSRs could clearly and effectively distinguish the 28 avocado accessions. In summary, our study is the first presentation of transcriptome data of different avocado ecotypes and comprehensive study on the development and analysis of a set of EST-SSR markers in avocado. The application of next-generation sequencing techniques for SSR development is a potentially powerful tool for genetic studies.