Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae)

Morphological Study on the Differentiation of Flower Buds and the Embryological Stages of Male and Female Floral Organs in Lespedeza davurica (Laxm.) Schindl. cv. JinNong (Fabaceae)

Lespedeza davurica (Laxm.) is a leguminous plant with significant ecological benefits, but its embryonic development mechanism remains unclear. We investigated the flower bud differentiation, megaspore and microspore formation, gametophyte development, and embryo and endosperm development in L. davurica. Our aim was to elucidate the relationship between the external morphology and internal development processes of male and female floral organs during growth, as well as the reproductive factors influencing fruiting. The results indicated that although the pistil develops later than the stamen during flower bud differentiation, both organs mature synchronously before flowering. L. davurica pollen exhibits three germination grooves, a reticulate outer wall, and papillary structures on the anther surface. In vivo pollination experiments revealed abnormal spiral growth of L. davurica pollen tubes within the style and the occurrence of callus plugs, which may reduce the seed setting rate. The anther wall development follows the dicotyledonous type, with tetrads formed through microspore meiosis exhibiting both left-right symmetry and tetrahedral arrangements. L. davurica has a single ovule, and the embryo sac develops in the monosporic polygonum type. After dormancy, the zygote undergoes multiple divisions, progressing through spherical, heart-shaped, and torpedo-shaped embryo stages, culminating in a mature embryo. A mature seed comprises cotyledons, hypocotyl, embryo, radicle, and seed coat. Phylogenetic tree analysis reveals a close genetic relationship between L. davurica and other leguminous plants from the genera Lespedeza and Medicago. This study provides valuable insights into the regulation of flowering and hybrid breeding in leguminous plants and offers a new perspective on the development of floral organs and seed setting rates.

Transcriptome analysis of Lantana camara flower petals reveals candidate anthocyanin biosynthesis genes mediating red flower color development

Flower color plays a crucial role in the appeal and selection of ornamental plants, directly influencing breeding strategies and the broader horticulture industry. Lantana camara, a widely favored flowering shrub, presents a rich palette of flower colors. Yet, the intricate molecular mechanisms governing this color variation in the species have remained largely unidentified. With the aim of filling this gap, this study embarked on a comprehensive de novo transcriptome assembly and differential gene expression analysis across 3 distinct lantana accessions, each showcasing a unique flower color. By harnessing the capabilities of both PacBio and Illumina sequencing platforms, a robust transcriptome assembly, encompassing 123,492 gene clusters and boasting 94.2% BUSCO completeness, was developed. The differential expression analysis unveiled 72,862 unique gene clusters that exhibited varied expression across different flower stages. A pronounced upregulation of 8 candidate core anthocyanin biosynthesis genes in the red-flowered accession was uncovered. This was further complemented by an upregulation of candidate MYB75 (PAP1) and bHLH42 (TT8) transcription factors. A candidate carotenoid cleavage dioxygenase (CCD4a) gene cluster also manifested a marked upregulation in white flowers. The study unveils the molecular groundwork of lantana's flower color variation, offering insights for future research and potential applications in breeding ornamental plants with desired color traits.

Comparative transcriptomics and gene expression divergence associated with homoploid hybrid speciation in Argyranthemum

Ecological isolation is increasingly thought to play an important role in speciation, especially for the origin and reproductive isolation of homoploid hybrid species. However, the extent to which divergent and/or transgressive gene expression changes are involved in speciation is not well studied. In this study, we employ comparative transcriptomics to investigate gene expression changes associated with the origin and evolution of two homoploid hybrid plant species, Argyranthemum sundingii and A. lemsii (Asteraceae). As there is no standard methodology for comparative transcriptomics, we examined five different pipelines for data assembly and analysing gene expression across the four species (two hybrid and two parental). We note biases and problems with all pipelines, and the approach used affected the biological interpretation of the data. Using the approach that we found to be optimal, we identify transcripts showing DE between the parental taxa and between the homoploid hybrid species and their parents; in several cases, putative functions of these DE transcripts have a plausible role in ecological adaptation and could be the cause or consequence of ecological speciation. Although independently derived, the homoploid hybrid species have converged on similar expression phenotypes, likely due to adaptation to similar habitats.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Transcriptome profiles of yellowish-white and fuchsia colored flowers in the Rheum palmatum complex reveal genes related to color polymorphism

Flower color variation is ubiquitous in many plant species, and several studies have been conducted to elucidate the underlying molecular mechanism. There are two flower color variants (yellowish-white and fuchsia) in the Rheum palmatum complex, however, few studies have investigated this phenomenon. Here, we used transcriptome sequencing of the two color variants to shed light on the molecular and biochemical basis for these color morphs. Comparison of the two transcriptomes identified 9641 differentially expressed unigenes (DEGs), including 6477 up-regulated and 3163 down-regulated genes. Functional analyses indicated that several DEGs were related to the anthocyanin biosynthesis pathway, and the expression profiles of these DEGs were coincident with the qRT-PCR validation results, indicating that expression levels of structural genes have a profound effect on the color variation in the R. palmatum complex. Our results suggested that the interaction of transcription factors (MYB, bHLH and WRKY) also regulated the anthocyanin biosynthesis in the R. palmatum complex. Estimation of selection pressures using the dN/dS ratio showed that 1106 pairs of orthologous genes have undergone positive selection. Of these positively selected genes, 21 were involved in the anthocyanin biosynthetic pathway, indicating that they may encode the proteins for structural alteration and affect flower color in the R. palmatum complex.

Co-expression network analysis of genes and networks associated with wheat pistillody

Crop male sterility has great value in theoretical research and breeding application. HTS-1, whose stamens transformed into pistils or pistil-like structures, is an important male sterility material selecting from Chinese Spring three-pistil (CSTP) wheat. However the molecular mechanism of pistillody development in HTS-1 remains a mystery. RNA-seq data of 11 wheat tissues were obtained from the National Center for Biotechnology Information (NCBI), including the stamens of CSTP and the pistils and pistillodic stamen of HTS-1. The Salmon program was utilized to quantify the gene expression levels of the 11 wheat tissues; and gene quantification results were normalized by transcripts per million (TPM). In total, 58,576 genes were used to construct block-wise network by co-expression networks analysis (WGCNA) R package. We obtained all of modules significantly associated with the 11 wheat tissues. AgriGO V2.0 was used to do Gene Ontology (GO) enrichment analysis; and genes and transcription factors (TFs) in these significant modules about wheat pistillody development were identified from GO enrichment results. Basic local alignment search tool (BLAST) was used to align HTS-1 proteins with the published pistillody-related proteins and TFs. Genes about wheat pistillody development were analyzed and validated by qRT-PCR. The MEturquoise, MEsaddlebrown, MEplum, MEcoral1, MElightsteelblue1, and MEdarkslateblue modules were significantly corelated to pistillodic stamen (correlation p < 0.05). Moreover, 206 genes related to carpel development (GO:0048440) or gynoecium development (GO:0048467) were identified only in the MEturquoise module by Gene Ontology (GO) analysis, and 42 of 206 genes were hub genes in MEturquoise module. qRT-PCR results showed that 38 of the 42 hub genes had highly expressed in pistils and pistillodic stamens than in stamens. A total of 15 pistillody development-related proteins were validated by BLAST. Transcription factors (TFs) were also analyzed in the MEturquoise module, and 618 TFs were identified. In total, 56 TFs from 11 families were considered to regulate the development of pistillodic stamen. The co-expression network showed that six of HB and three of BES1 genes were identified in 42 hub genes. This indicated that TFs played important roles in wheat pistillody development. In addition, there were 11 of ethylene-related genes connected with TFs or hub genes, suggesting the important roles of ethylene-related genes in pistillody development. These results provide important insights into the molecular interactions underlying pistillody development.

Testing candidate genes linked to corolla shape variation of a pollinator shift in Rhytidophyllum (Gesneriaceae)

Floral adaptations to specific pollinators like corolla shape variation often result in reproductive isolation and thus speciation. But despite their ecological importance, the genetic bases of corolla shape transitions are still poorly understood, especially outside model species. Hence, our goal was to identify candidate genes potentially involved in corolla shape variation between two closely related species of the Rhytidophyllum genus (Gesneriaceae family) from the Antilles with contrasting pollination strategies. Rhytidophyllum rupincola has a tubular corolla and is strictly pollinated by hummingbirds, whereas R. auriculatum has more open flowers and is pollinated by hummingbirds, bats, and insects. We surveyed the literature and used a comparative transcriptome sequence analysis of synonymous and non-synonymous nucleotide substitutions to obtain a list of genes that could explain floral variation between R. auriculatum and R. rupincola. We then tested their association with corolla shape variation using QTL mapping in a F₂ hybrid population. Out of 28 genes tested, three were found to be good candidates because of a strong association with corolla shape: RADIALIS, GLOBOSA, and JAGGED. Although the role of these genes in Rhytidophyllum corolla shape variation remains to be confirmed, these findings are a first step towards identifying the genes that have been under selection by pollinators and thus involved in reproductive isolation and speciation in this genus.

Loss of anthocyanidin synthase gene is associated with white flowers of Salvia miltiorrhiza Bge. f. alba, a natural variant of S. miltiorrhiza

SmANS deletion leads to white flower mutation in Salvia miltiorrhiza. SmANS deletion leads to white flower mutation in Salvia miltiorrhiza. Abstract Salvia miltiorrhiza is an essential traditional Chinese medicine (TCM) with purple flowers, and S. miltiorrhiza Bge. f. alba is a unique intraspecific variation with white flowers. The molecular mechanism of flower color formation in S. miltiorrhiza will provide vital information for the variation and evolution. Here, we performed HPLC, transcriptomic, and re-sequencing analyses of purple-flowered S. miltiorrhiza line 'Zihua105' (ZH105) and white-flowered S. miltiorrhiza Bge. f. alba line 'Baihua18' (BH18). Delphinidin was the most anthocyanidin in ZH105, which become the main different between ZH105 vs. BH18 flowers. Transcriptome analysis revealed 299 differentially expressed genes (DEGs). SmANS, the anthocyanidin synthase gene in the down-stream anthocyanin biosynthesis pathway, was significantly expressed in ZH105 corollas, suggesting it might play a key role in white petal formation. Whole-genome re-sequencing revealed that a 6.75 kb segment located on chromosome 5, which contains the complete sequence of the SmANS genes, was lost in BH18 and another S. miltiorrhiza Bge. f. alba line. In contrast, the other five purple-flowered S. miltiorrhiza lines both possessed this segment. Further molecular marker identification also confirmed that wild S. miltiorrhiza Bge. f. alba lines lost regions that contained a complete or important part of SmANS sequences. Subsequently, the research showed that the deletion mutant of SmANS genes resulted in the natural white flower color variant of S. miltiorrhiza.

An Integrated Analysis of Transcriptome and miRNA Sequencing Provides Insights into the Dynamic Regulations during Flower Morphogenesis in Petunia

Published genome sequences can facilitate multiple genome sequencing studies of flower development, which can serve as the basis for later analysis of variation in flower phenotypes. To identify potential regulators related to flower morphology, we captured dynamic expression patterns under five different developmental stages of petunia flowers, a popular bedding plant, using transcriptome and miRNA sequencing. The significant transcription factor (TF) families, including MYB, MADS, and bHLH, were elucidated. MADS-box genes exhibited co-expression patterns with BBR-BPC, GATA, and Dof genes in different modules according to a weighted gene co-expression network analysis. Through miRNA sequencing, a total of 45 conserved and 26 novel miRNAs were identified. According to GO and KEGG enrichment analysis, the carbohydrate metabolic process, photosynthesis, and phenylalanine metabolism were significant at the transcriptomic level, while the response to hormone pathways was significantly enriched by DEmiR-targeted genes. Finally, an miRNA–RNA network was constructed, which suggested the possibility of novel miRNA-mediated regulation pathways being activated during flower development. Overall, the expression data in the present study provide novel insights into the developmental gene regulatory network facilitated by TFs, miRNA, and their target genes.

Genetic engineering of novel yellow color african violet (Saintpaulia ionantha) produced by accumulation of Aureusidin 6-O-glucoside

BACKGROUND

Flower color is one of the main characteristics of ornamental plants. Aurones are light yellow flavonoids produced in the petals of a limited number of plant species including snapdragon (Antirrhinum majus). As a commercially-recognized species, African violet can be found in various colors except yellow. This research, aiming at changing the petals' color of African violet from white to yellow, was conducted using the simultaneous expressions of chalcone 4'-O-glucosyltransferase (4'CGT) and aureusidin synthase (AS1) genes without the need for silencing anthocyanin biosynthesis pathway genes via both transient and stable transfer methods.

RESULTS

The transient gene transfer among transgenic plants led to a clear change of petals' color from white to light yellow. This occurs while no change was observed in non-transgenic (Wild type) petals. In total, 15 positive transgenic plants, produced via stable gene transfer, were detected. Moreover, since their flower color was yellow, both genes were present. Meanwhile, the corresponding transformation yield was determined 20-30%. The transformation, expression and integration of genes among T0 transgenic plants were verified using the PCR, qRT-PCR and Southern blotting techniques, respectively. Furthermore, the probable color change of petals' cross-section and existence of Aureusidin 6-O-glucoside (AOG) compound were determined using a light microscope and HPLC-DAD-MSn analysis, correspondingly.

CONCLUSIONS

Generally, the creation of aurones biosynthesis pathway is only viable through the simultaneous expression of genes which leads to color change of African violet's petal from white to yellow. This conclusion can lead to an effective strategy to produce yellow color in ornamental plant species.

A Survey of Enhanced Cold Tolerance and Low-Temperature-Induced Anthocyanin Accumulation in a Novel Zoysia japonica Biotype

Zoysia japonica is a warm-season turfgrass that is extensively used in landscaping, sports fields, and golf courses worldwide. Uncovering the low-temperature response mechanism of Z. japonica can help to accelerate the development of new cold-tolerant cultivars, which could be used to prolong the ornamental and usage duration of turf. A novel Z. japonica biotype, YueNong-9 (YN-9), was collected from northeastern China for this study. Phenotypic measurements, cold-tolerance investigation, and whole-transcriptome surveys were performed on YN-9 and LanYin-3 (LY-3), the most popular Z. japonica cultivar in Southern China. The results indicated the following: YN-9 has longer second and third leaves than LY-3; when exposed to the natural low temperature during winter in Guangzhou, YN-9 accumulated 4.74 times more anthocyanin than LY-3; after cold acclimation and freezing treatment, 83.25 ± 9.55% of YN-9 survived while all LY-3 leaves died, and the dark green color index (DGCI) value of YN-9 was 1.78 times that of LY-3; in YN-9, there was a unique up-regulation of Phenylalanine ammonia-lyase (PAL), Homeobox-leucine Zipper IV (HD-ZIP), and ATP-Binding Cassette transporter B8 (ABCB8) expressions, as well as a unique down-regulation of zinc-regulated transporters and iron-regulated transporter-like proteins (ZIPs) expression, which may promote anthocyanin biosynthesis, transport, and accumulation. In conclusion, YN-9 exhibited enhanced cold tolerance and is thus an excellent candidate for breeding cold-tolerant Z. japonica variety, and its unique low-temperature-induced anthocyanin accumulation and gene responses provide ideas and candidate genes for the study of low-temperature tolerance mechanisms and genetic engineering breeding.

Species and population specific gene expression in blood transcriptomes of marine turtles

BACKGROUND

Transcriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations.

RESULTS

We generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance.

CONCLUSIONS

Obtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.

Transgenerational Plasticity in Flower Color Induced by Caterpillars

Variation in flower color due to transgenerational plasticity could stem directly from abiotic or biotic environmental conditions. Finding a link between biotic ecological interactions across generations and plasticity in flower color would indicate that transgenerational effects of ecological interactions, such as herbivory, might be involved in flower color evolution. We conducted controlled experiments across four generations of wild radish (Raphanus sativus, Brassicaceae) plants to explore whether flower color is influenced by herbivory, and to determine whether flower color is associated with transgenerational chromatin modifications. We found transgenerational effects of herbivory on flower color, partly related to chromatin modifications. Given the presence of herbivory in plant populations worldwide, our results are of broad significance and contribute to our understanding of flower color evolution.

Changes at a Critical Branchpoint in the Anthocyanin Biosynthetic Pathway Underlie the Blue to Orange Flower Color Transition in Lysimachia arvensis

Anthocyanins are the primary pigments contributing to the variety of flower colors among angiosperms and are considered essential for survival and reproduction. Anthocyanins are members of the flavonoids, a broader class of secondary metabolites, of which there are numerous structural genes and regulators thereof. In western European populations of Lysimachia arvensis, there are blue- and orange-petaled individuals. The proportion of blue-flowered plants increases with temperature and daylength yet decreases with precipitation. Here, we performed a transcriptome analysis to characterize the coding sequences of a large group of flavonoid biosynthetic genes, examine their expression and compare our results to flavonoid biochemical analysis for blue and orange petals. Among a set of 140 structural and regulatory genes broadly representing the flavonoid biosynthetic pathway, we found 39 genes with significant differential expression including some that have previously been reported to be involved in similar flower color transitions. In particular, F3'5'H and DFR, two genes at a critical branchpoint in the ABP for determining flower color, showed differential expression. The expression results were complemented by careful examination of the SNPs that differentiate the two color types for these two critical genes. The decreased expression of F3'5'H in orange petals and differential expression of two distinct copies of DFR, which also exhibit amino acid changes in the color-determining substrate specificity region, strongly correlate with the blue to orange transition. Our biochemical analysis was consistent with the transcriptome data indicating that the shift from blue to orange petals is caused by a change from primarily malvidin to largely pelargonidin forms of anthocyanins. Overall, we have identified several flavonoid biosynthetic pathway loci likely involved in the shift in flower color in L. arvensis and even more loci that may represent the complex network of genetic and physiological consequences of this flower color polymorphism.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

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.

Performance of gene expression analyses using de novo assembled transcripts in polyploid species

MOTIVATION

Quality of gene expression analyses using de novo assembled transcripts in species that experienced recent polyploidization remains unexplored.

RESULTS

Differential gene expression (DGE) analyses using putative genes inferred by Trinity, Corset and Grouper performed slightly differently across five plant species that experienced various polyploidy histories. In species that lack recent polyploidy events that occurred in the past several millions of years, DGE analyses using de novo assembled transcriptomes identified 54-82% of the differentially expressed genes recovered by mapping reads to the reference genes. However, in species that experienced more recent polyploidy events, the percentage decreased to 21-65%. Gene co-expression network analyses using de novo assemblies versus mapping to the reference genes recovered the same module that significantly correlated with treatment in one species that lacks recent polyploidization.

AVAILABILITY AND IMPLEMENTATION

Commands and scripts used in this study are available at https://bitbucket.org/lychen83/chen_et_al_2018_benchmark_dge/; Analysis files are available at Dryad doi: 10.5061/dryad.4p6n481.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Co-expression clustering across flower development identifies modules for diverse floral forms in Achimenes (Gesneriaceae)

BACKGROUND

Genetic pathways involved with flower color and shape are thought to play an important role in the development of flowers associated with different pollination syndromes, such as those associated with bee, butterfly, or hummingbird pollination. Because pollination syndromes are complex traits that are orchestrated by multiple genes and pathways, the gene regulatory networks have not been explored. Gene co-expression networks provide a systems level approach to identify important contributors to floral diversification.

METHODS

RNA-sequencing was used to assay gene expression across two stages of flower development (an early bud and an intermediate stage) in 10 species of Achimenes (Gesneriaceae). Two stage-specific co-expression networks were created from 9,503 orthologs and analyzed to identify module hubs and the network periphery. Module association with bee, butterfly, and hummingbird pollination syndromes was tested using phylogenetic mixed models. The relationship between network connectivity and evolutionary rates (d N/d S) was tested using linear models.

RESULTS

Networks contained 65 and 62 modules that were largely preserved between developmental stages and contained few stage-specific modules. Over a third of the modules in both networks were associated with flower color, shape, and pollination syndrome. Within these modules, several hub nodes were identified that related to the production of anthocyanin and carotenoid pigments and the development of flower shape. Evolutionary rates were decreased in highly connected genes and elevated in peripheral genes.

DISCUSSION

This study aids in the understanding of the genetic architecture and network properties underlying the development of floral form and provides valuable candidate modules and genes for future studies.

Temporal Distinction between Male and Female Floral Organ Development in Nicotiana tabacum cv. Xanthi (Solanaceae)

Early floral developmental investigations provide crucial evidence for phylogenetic and molecular studies of plants. The developmental and evolutionary mechanisms underlying the variations in floral organs are critical for a thorough understanding of the diversification of flowers. Ontogenetic comparisons between anthers and pistil within single flowers were characterized over time in Nicotiana tabacum cv. Xanthi. The ages of 42 tobacco flower or flower primordia were estimated using corolla growth analysis. Results showed that the protodermal layer in carpel primordia contributes to carpel development by both anticlinal and periclinal divisions. Periclinal divisions in the hypodermal layer of the placenta were observed around 4.8 ± 1.3 days after the formation of early carpel primordia (ECP) and ovule initiation occurred 10.0 ± 0.5 days after ECP. Meiosis in anthers and ovules began about 8.9 ± 1.1 days and 14.4 ± 1.3 days after ECP, respectively. Results showed an evident temporal distinction between megasporogenesis and microsporogenesis. Flower ages spanned a 17-day interval, starting with flower primordia containing the ECP and anther primordia to the tetrad stage of meiosis in megasporocytes and the bicellular stage in pollen grains. These results establish a solid foundation for future studies in order to identify the developmental and molecular mechanisms responsible for the mating system in tobacco.

Chemical Basis of Floral Color Signals in Gesneriaceae: The Effect of Alternative Anthocyanin Pathways

Changes in floral pigmentation can have dramatic effects on angiosperm evolution by making flowers either attractive or inconspicuous to different pollinator groups. Flower color largely depends on the type and abundance of pigments produced in the petals, but it is still unclear whether similar color signals rely on same biosynthetic pathways and to which extent the activation of certain pathways influences the course of floral color evolution. To address these questions, we investigated the physical and chemical aspects of floral color in the Neotropical Gesnerioideae (ca. 1,200 spp.), in which two types of anthocyanins, hydroxyanthocyanins, and deoxyanthocyanins, have been recorded as floral pigments. Using spectrophotometry, we measured flower reflectance for over 150 species representing different clades and pollination syndromes. We analyzed these reflectance data to estimate how the Gesnerioideae flowers are perceived by bees and hummingbirds using the visual system models of these pollinators. Floral anthocyanins were further identified using high performance liquid chromatography coupled to mass spectrometry. We found that orange/red floral colors in Gesnerioideae are produced either by deoxyanthocyanins (e.g., apigenidin, luteolinidin) or hydroxyanthocyanins (e.g., pelargonidin). The presence of deoxyanthocyanins in several lineages suggests that the activation of the deoxyanthocyanin pathway has evolved multiple times in the Gesnerioideae. The hydroxyanthocyanin-producing flowers span a wide range of colors, which enables them to be discriminated by hummingbirds or bees. By contrast, color diversity among the deoxyanthocyanin-producing species is lower and mainly represented at longer wavelengths, which is in line with the hue discrimination optima for hummingbirds. These results indicate that Gesnerioideae have evolved two different biochemical mechanisms to generate orange/red flowers, which is associated with hummingbird pollination. Our findings also suggest that the activation of the deoxyanthocyanin pathway has restricted flower color diversification to orange/red hues, supporting the potential constraining role of this alternative biosynthetic pathway on the evolutionary outcome of phenotypical and ecological diversification.

In situ transcriptomic and metabolomic study of the loss of photosynthesis in the leaves of mixotrophic plants exploiting fungi

Mycoheterotrophic plants have lost photosynthesis and obtain carbon through mycorrhizal fungi colonizing their roots. They are likely to have evolved from mixotrophic ancestors, which rely on both photosynthesis and fungal carbon for their development. Whereas our understanding of the ecological and genomic changes associated with the evolutionary shift to mycoheterotrophy is deepening, little information is known about the specific metabolic and physiological features driving this evolution. We investigated this issue in naturally occurring achlorophyllous variants of temperate mixotrophic orchids. We carried out an integrated transcriptomic and metabolomic analysis of the response to achlorophylly in the leaves of three mixotrophic species sampled in natura. Achlorophyllous leaves showed major impairment of their photosynthetic and mineral nutrition functions, strong accumulation of free amino acids, overexpression of enzymes and transporters related to sugars, amino acids and fatty acid catabolism, as well as induction of some autophagy-related and biotic stress genes. Such changes were reminiscent of these reported for variegated leaves and appeared to be symptomatic of a carbon starvation response. Rather than decisive metabolic innovations, we suggest that the evolution towards mycoheterotrophy in orchids is more likely to be reliant on the versatility of plant metabolism and an ability to exploit fungal organic resources, especially amino acids, to replace missing photosynthates.

Comparative transcriptomics sheds light on differential adaptation and species diversification between two Melastoma species and their F1 hybrid

Variation in gene expression has been shown to promote adaptive divergence, and can lead to speciation. The plant genus Melastoma, thought to have diversified through adaptive radiation, provides an excellent model for the study of gene expressional changes during adaptive differentiation and following interspecific hybridization. In this study, we performed RNA-seq on M. candidum, M. sanguineum and their F1 hybrid, to investigate the role of gene expression in species diversification within the genus. Reference transcriptomes were assembled using combined data from both parental species, resulting in 50 519 and 48 120 transcripts for the leaf and flower petal, after removing redundancy. Differential expression analysis uncovered 3793 and 2116 differentially expressed (DE) transcripts, most of which are between M. candidum and M. sanguineum. Differential expression was observed for genes related to light responses, as well as genes that regulate the development of leaf trichomes, a trait that among others is thought to protect plants against sunlight, suggesting the differential adaptation of the species to sunlight intensity. The analysis of positively selected genes between the two species also revealed possible differential adaptation to other abiotic stresses such as drought and temperature. In the hybrid, almost all possible modes of expression were observed at the DE transcripts, although at most transcripts, the expression levels were similar to that of either parent instead of being intermediate. A small number of transgressively expressed transcripts that matched genes known to promote plant growth and adaptation to stresses in new environments were also found, possibly explaining the vigour observed in the hybrid. The findings in this study provided insights into the role of gene expression in the diversification of Melastoma, which we believe is an important example for more cross-taxa comparisons in the future.

Unraveling the Developmental and Genetic Mechanisms Underpinning Floral Architecture in Proteaceae

Proteaceae are a basal eudicot family with a highly conserved floral groundplan but which displays considerable variation in other aspects of floral and inflorescence morphology. Their morphological diversity and phylogenetic position make them good candidates for understanding the evolution of floral architecture, in particular the question of the homology of the undifferentiated perianth with the differentiated perianth of core eudicots, and the mechanisms underlying the repeated evolution of zygomorphy. In this paper, we combine a morphological approach to explore floral ontogenesis and a transcriptomic approach to access the genes involved in floral organ identity and development, focusing on Grevillea juniperina, a species from subfamily Grevilleoideae. We present developmental data for Grevillea juniperina and three additional species that differ in their floral symmetry using stereomicroscopy, SEM and High Resolution X-Ray Computed Tomography. We find that the adnation of stamens to tepals takes place at early developmental stages, and that the establishment of bilateral symmetry coincides with the asymmetrical growth of the single carpel. To set a framework for understanding the genetic basis of floral development in Proteaceae, we generated and annotated de novo a reference leaf/flower transcriptome from Grevillea juniperina. We found Grevillea homologs of all lineages of MADS-box genes involved in floral organ identity. Using Arabidopsis thaliana gene expression data as a reference, we found homologs of other genes involved in floral development in the transcriptome of G. juniperina. We also found at least 21 class I and class II TCP genes, a gene family involved in the regulation of growth processes, including floral symmetry. The expression patterns of a set of floral genes obtained from the transcriptome were characterized during floral development to assess their organ specificity and asymmetry of expression.

Lessons from Cleomaceae, the Sister of Crucifers

Cleomaceae is a diverse group well-suited to addressing fundamental genomic and evolutionary questions as the sister group to Brassicaceae, facilitating transfer of knowledge from the model Arabidopsis thaliana. Phylogenetic and taxonomic revisions provide a framework for examining the evolution of substantive morphological and physiology diversity in Cleomaceae, but not necessarily in Brassicaceae. The investigation of both nested and contrasting whole-genome duplications (WGDs) between Cleomaceae and Brassicaceae allows comparisons of independently duplicated genes and investigation of whether they may be drivers of the observed innovations. Further, a wealth of outstanding genetic research has provided insight into how the important alternative carbon fixation pathway, C₄ photosynthesis, has evolved via differential expression of a suite of genes, of which the underlying mechanisms are being elucidated.

Early responses to dehydration in contrasting wild Arachis species

Wild peanut relatives (Arachis spp.) are genetically diverse and were selected throughout evolution to a range of environments constituting, therefore, an important source of allelic diversity for abiotic stress tolerance. In particular, A. duranensis and A. stenosperma, the parents of the reference Arachis A-genome genetic map, show contrasting transpiration behavior under limited water conditions. This study aimed to build a comprehensive gene expression profile of these two wild species under dehydration stress caused by the withdrawal of hydroponic nutrient solution. For this purpose, roots of both genotypes were collected at seven time-points during the early stages of dehydration and used to construct cDNA paired-end libraries. Physiological analyses indicated initial differences in gas exchange parameters between the drought-tolerant genotype of A. duranensis and the drought-sensitive genotype of A. stenosperma. High-quality Illumina reads were mapped against the A. duranensis reference genome and resulted in the identification of 1,235 and 799 Differentially Expressed Genes (DEGs) that responded to the stress treatment in roots of A. duranensis and A. stenosperma, respectively. Further analysis, including functional annotation and identification of biological pathways represented by these DEGs confirmed the distinct gene expression behavior of the two contrasting Arachis species genotypes under dehydration stress. Some species-exclusive and common DEGs were then selected for qRT-PCR analysis, which corroborated the in silico expression profiling. These included genes coding for regulators and effectors involved in drought tolerance responses, such as activation of osmosensing molecular cascades, control of hormone and osmolyte content, and protection of macromolecules. This dataset of transcripts induced during the dehydration process in two wild Arachis genotypes constitute new tools for the understanding of the distinct gene regulation processes in these closely related species but with contrasting drought responsiveness. In addition, our findings provide insights into the nature of drought tolerance in wild germoplasm, which might be explored as novel sources of diversity and useful wild alleles to develop climate-resilient crop varieties.

Phylogenomic analyses reveal extensive gene flow within the magic flowers (Achimenes)

PREMISE OF THE STUDY

The Neotropical Gesneriaceae is a lineage known for its colorful and diverse flowers, as well as an extensive history of intra- and intergeneric hybridization, particularly among Achimenes (the magic flowers) and other members of subtribe Gloxiniinae. Despite numerous studies seeking to elucidate the evolutionary relationships of these lineages, relatively few have sought to infer specific patterns of gene flow despite evidence of widespread hybridization.

METHODS

To explore the utility of phylogenomic data for reassessing phylogenetic relationships and inferring patterns of gene flow among species of Achimenes, we sequenced 12 transcriptomes. We used a variety of methods to infer the species tree, examine gene tree discordance, and infer patterns of gene flow.

KEY RESULTS

Phylogenomic analyses resolve clade relationships at the crown of the lineage with strong support. In contrast to previous analyses, we recovered strong support for several new relationships despite a significant amount of gene tree discordance. We present evidence for at least two introgression events between two species pairs that share pollinators, and suggest that the species status of Achimenes admirabilis be reexamined.

CONCLUSIONS

Our study demonstrates the utility of transcriptome data for phylogenomic analyses, and inferring patterns of gene flow despite gene tree discordance. Moreover, these data provide another example of prevalent interspecific gene flow among Neotropical plants that share pollinators.