Transcriptome analysis of a petal anthocyanin polymorphism in the arctic mustard, Parrya nudicaulis

Genetic factors explaining anthocyanin pigmentation differences

BACKGROUND

Anthocyanins are important contributors to coloration across a wide phylogenetic range of plants. Biological functions of anthocyanins span from reproduction to protection against biotic and abiotic stressors. Owing to a clearly visible phenotype of mutants, the anthocyanin biosynthesis and its sophisticated regulation have been studied in numerous plant species. Genes encoding the anthocyanin biosynthesis enzymes are regulated by a transcription factor complex comprising MYB, bHLH and WD40 proteins.

RESULTS

A systematic comparison of anthocyanin-pigmented vs. non-pigmented varieties was performed within numerous plant species covering the taxonomic diversity of flowering plants. The literature was screened for cases in which genetic factors causing anthocyanin loss were reported. Additionally, transcriptomic data sets from four previous studies were reanalyzed to determine the genes possibly responsible for color variation based on their expression pattern. The contribution of different structural and regulatory genes to the intraspecific pigmentation differences was quantified. Differences concerning transcription factors are by far the most frequent explanation for pigmentation differences observed between two varieties of the same species. Among the transcription factors in the analyzed cases, MYB genes are significantly more prone to account for pigmentation differences compared to bHLH or WD40 genes. Among the structural genes, DFR genes are most often associated with anthocyanin loss.

CONCLUSIONS

These findings support previous assumptions about the susceptibility of transcriptional regulation to evolutionary changes and its importance for the evolution of novel coloration phenotypes. Our findings underline the particular significance of MYBs and their apparent prevalent role in the specificity of the MBW complex.

Transcriptomic and metabolomic analyses reveal the altitude adaptability and evolution of different-colored flowers in alpine Rhododendron species

Understanding the molecular mechanisms and evolutionary process of plant adaptation to the heterogeneous environment caused by altitude gradients in plateau mountain ecosystems can provide novel insight into species' responses to global changes. Flower color is the most conspicuous and highly diverse trait in nature. Herein, the gene expression patterns, evolutionary adaptation and metabolites changes of different-colored flowers of alpine Rhododendron L. species along altitude gradients were investigated based on a combined analysis of transcriptomics and metabolomics. Differentially expressed genes were found to be related to the biosynthesis of carbohydrates, fatty acids, amino acids and flavonoids, suggesting their important roles in the altitude adaptability of Rhododendron species. The evolution rate of high-altitude species was faster than that of low-altitude species. Genes related to DNA repair, mitogen-activated protein kinase and ABA signal transduction, and lipoic acid and propanoate metabolism were positively selected in the flowers of high-altitude Rhododendron species and those associated with carotenoid biosynthesis pathway, ABA signal transduction and ethylene signal transduction were positively selected in low-altitude species. These results indicated that the genes with differentiated expressions or functions exhibit varying evolution during the adaptive divergence of heterogeneous environment caused by altitude gradients. Flower-color variation might be attributed to the significant differences in gene expression or metabolites related to sucrose, flavonoids and carotenoids at the transcription or metabolism levels of Rhododendron species. This work suggests that Rhododendron species have multiple molecular mechanisms in their adaptation to changing environments caused by altitude gradients.

Phenotypic selection on floral traits in the arctic plant Parrya nudicaulis (Brassicaceae)

The evolution of floral traits is often attributed to pollinator‐mediated selection; however, the importance of pollinators as selective agents in arctic environments is poorly resolved. In arctic and subarctic regions that are thought to be pollen limited, selection is expected to either favor floral traits that increase pollinator attraction or promote reproductive assurance through selfing. We quantified phenotypic selection on floral traits in two arctic and two subarctic populations of the self‐compatible, but largely pollinator‐dependent, Parrya nudicaulis. Additionally, we measured selection in plants in both open pollination and pollen augmentation treatments to estimate selection imposed by pollinators in one population. Seed production was found to be limited by pollen availability and strong directional selection on flower number was observed. We did not detect consistently greater magnitudes of selection on floral traits in the arctic relative to the subarctic populations. Directional selection for more pigmented flowers in one arctic population was observed, however. In some populations, selection on flower color was found to interact with other traits. We did not detect consistently stronger selection gradients across all traits for plants exposed to pollinator selection relative to those in the pollen augmentation treatment; however, directional selection tended to be higher for some floral traits in open‐pollinated plants.

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.

Transcriptomics analysis of Psidium cattleyanum Sabine (Myrtaceae) unveil potential genes involved in fruit pigmentation

Psidium cattleyanum Sabine is an Atlantic Forest native species that presents some populations with red fruits and others with yellow fruits. This variation in fruit pigmentation in this species is an intriguing character that could be related to species evolution but still needs to be further explored. Our goal was to provide genomic information for these morphotypes to understand the molecular mechanisms of differences in fruit colour in this species. In this study, we performed a comparative transcriptome analysis of red and yellow morphotypes of P. cattleyanum, considering two stages of fruit ripening. The transcriptomic analysis performed encompassing leaves, unripe and ripe fruits, in triplicate for each morphotype. The transcriptome consensus from each morphotype showed 301,058 and 298,310 contigs from plants with yellow and red fruits, respectively. The differential expression revealed important genes that were involved in anthocyanins biosynthesis, such as the anthocyanidin synthase (ANS) and UDP-glucose:flavonoid-o-glucosyltransferase (UFGT) that were differentially regulated during fruit ripening. This study reveals stimulating data for the understanding of the pathways and mechanisms involved in the maturation and colouring of P. cattleyanum fruits and suggests that the ANS and UFGT genes are key factors involved in the synthase and pigmentation accumulation in red fruits.

MaCDSP32 From Mulberry Enhances Resilience Post-drought by Regulating Antioxidant Activity and the Osmotic Content in Transgenic Tobacco

Desiccation tolerance is a complex phenomenon that depends on the regulated expression of numerous genes during dehydration and subsequent rehydration. Our previous study identified a chloroplast drought-induced stress protein (MaCDSP32) in mulberry, a thioredoxin (Trx) that is upregulated under drought conditions and is likely to confer drought tolerance to transgenic plants. Mulberry (Morus spp.) is an ecologically and economically important perennial woody plant that is widely used in forest management to combat desertification. However, its stress tolerance physiology is not well understood. In this study, the functions of MaCDSP32 gene were investigated. The expression of MaCDSP32 exhibited a circadian rhythm and was induced by mild and severe water deficits. Under abiotic stress, MaCDSP32-overexpressing plants exhibited increased stress sensitivity with lower water retention capacity and more severe lipid peroxidation than the wild-type (WT) plants. Furthermore, the activity of superoxide dismutase (SOD), the contents of proline and soluble sugars and the expression of stress-related transcription factors were lower in the MaCDSP32-overexpressing plants than in the WT plants. However, the MaCDSP32-overexpressing lines exhibited stronger recovery capability after rewatering post-drought. Moreover, the SOD enzyme activity, proline content, and soluble sugar content were higher in the transgenic plants after rewatering than in the WT plants. The production of the reactive oxygen species (ROS) H₂O₂ and O₂ ^- was significantly lower in the transgenic plants than in the WT plants. In addition, under abiotic stress, the MaCDSP32-overexpressing lines exhibited improved seed germination and seedling growth, these effects were regulated by a positive redox reaction involving MaCDSP32 and one of its targets. In summary, this study indicated that MaCDSP32 from mulberry regulates plant drought tolerance and ROS homeostasis mainly by controlling SOD enzyme activity and proline and soluble sugar concentrations and that this control might trigger the stress response during seed germination and plant growth. Overall, MaCDSP32 exerts pleiotropic effects on the stress response and stress recovery in plants.

Unexpected predicted length variation for the coding sequence of the sleep related gene, BHLHE41 in gorilla amidst strong purifying selection across mammals

There is a molecular basis for many sleep patterns and disorders involving circadian clock genes. In humans, "short-sleeper" behavior has been linked to specific amino acid substitutions in BHLHE41 (DEC2), yet little is known about variation at these sites and across this gene in mammals. We compare BHLHE41 coding sequences for 27 mammals. Approximately half of the coding sequence was invariable at the nucleotide level and close to three-quarters of the amino acid alignment was identical. No other mammals had the same "short-sleeper" amino acid substitutions previously described from humans. Phylogenetic analyses based on the nucleotides of the coding sequence alignment are consistent with established mammalian relationships confirming orthology among the sampled sequences. Significant purifying selection was detected in about two-thirds of the variable codons and no codons exhibited significant signs of positive selection. Unexpectedly, the gorilla BHLHE41 sequence has a 318 bp insertion at the 5' end of the coding sequence and a deletion of 195 bp near the 3' end of the coding sequence (including the two short sleeper variable sites). Given the strong signal of purifying selection across this gene, phylogenetic congruence with expected relationships and generally conserved function among mammals investigated thus far, we suggest the indels predicted in the gorilla BHLHE41 may represent an annotation error and warrant experimental validation.

Comparative Transcriptomics Provides Insight into Floral Color Polymorphism in a Pleione limprichtii Orchid Population

Floral color polymorphism can provide great insight into species evolution from a genetic and ecological standpoint. Color variations between species are often mediated by pollinators and are fixed characteristics, indicating their relevance to adaptive evolution, especially between plants within a single population or between similar species. The orchid genus Pleione has a wide variety of flower colors, from violet, rose-purple, pink, to white, but their color formation and its evolutionary mechanism are unclear. Here, we selected the P. limprichtii population in Huanglong, Sichuan Province, China, which displayed three color variations: Rose-purple, pink, and white, providing ideal material for exploring color variations with regard to species evolution. We investigated the distribution pattern of the different color morphs. The ratio of rose-purple:pink:white-flowered individuals was close to 6:3:1. We inferred that the distribution pattern may serve as a reproductive strategy to maintain the population size. Metabolome analysis was used to reveal that cyanindin derivatives and delphidin are the main color pigments involved. RNA sequencing was used to characterize anthocyanin biosynthetic pathway-related genes and reveal different color formation pathways and transcription factors in order to identify differentially-expressed genes and explore their relationship with color formation. In addition, qRT-PCR was used to validate the expression patterns of some of the genes. The results show that PlFLS serves as a crucial gene that contributes to white color formation and that PlANS and PlUFGT are related to the accumulation of anthocyanin which is responsible for color intensity, especially in pigmented flowers. Phylogenetic and co-expression analyses also identified a R2R3-MYB gene PlMYB10, which is predicted to combine with PlbHLH20 or PlbHLH26 along with PlWD40-1 to form an MBW protein complex (MYB, bHLH, and WDR) that regulates PlFLS expression and may serve as a repressor of anthocyanin accumulation-controlled color variations. Our results not only explain the molecular mechanism of color variation in P. limprichtii, but also contribute to the exploration of a flower color evolutionary model in Pleione, as well as other flowering plants.

Global gene expression and pigment analysis of two contrasting flower color cultivars of Canna

Development of flower color in plants is a complex process. Among others, it is an important trait for ornamental flowering plants. Canna is a flowering ornamental plant of family Cannaceae. To understand the molecular mechanism of flower color development in Canna, RNA sequencing from flower tissues of two contrasting flower color cultivars, Red President (RP) and Tropical Sunrise (TS) was performed. More than 27.0 million and 19.0 million clean reads were obtained from RP and TS, respectively. The combined clean reads were assembled into 147,295 unigenes. The Canna unigenes showed maximum homology with Populus trichocarpa (26.79%). A total of 2702 unigenes expressed differentially between the two cultivars of which 1972 were up-regulated and 730 were down-regulated in RP. Phenylpropanoid and flavonoid biosynthetic processes were the significant processes in RP. Expression of a vast number of transcription factors including MYB, bHLH, ARF, and WRKY were higher in RP than TS. The expression analysis of RNA sequencing data was validated by qRT-PCR analysis. Further, concentration of measured anthocyanidins and flavonols were very low or absent in TS, corroborating largely with our transcriptome data. These findings may help in understanding flower color development in Canna and in future crop breeding program.

Flower colour polymorphism in the Mediterranean Basin: occurrence, maintenance and implications for speciation

Flower colour polymorphism (FCP) is the occurrence of at least two discrete flower colour variants in the same population. Despite a vast body of research concerning the maintenance and evolutionary consequences of FCP, only recently has the spatial variation in morph frequencies among populations been explored. Here we summarise the biochemical and genetic basis of FCP, the factors that have been proposed to explain their maintenance, and the importance of FCP and its geographic variation in the speciation process. We also review the incidence of FCP in the environmentally heterogeneous Mediterranean Basin. Nearly 88% of Mediterranean FCP species showed anthocyanin-based polymorphisms. Concerning the evolutionary mechanisms that contribute to maintain FCP, selection by pollinators is suggested in some species, but in others, selection by non-pollinator agents, genetic drift or gene flow are also found; in some cases different processes interact in the maintenance of FCP. We emphasise the role of both autonomous selfing and clonal reproduction in FCP maintenance. Mediterranean polymorphic species show mainly monomorphic populations with only a few polymorphic ones, which generate clinal or mosaic patterns of variation in FCP. No cases of species with only polymorphic populations were found. We posit that different evolutionary processes maintaining polymorphism the Mediterranean Basin will result in a continuum of geographic patterns in morph compositions and relative frequencies of FCP species.

Transcriptome analysis of terpene chemotypes of Melaleuca alternifolia across different tissues

Plant chemotypes or chemical polymorphisms are defined by discrete variation in secondary metabolites within a species. This variation can have consequences for ecological interactions or the human use of plants. Understanding the molecular basis of chemotypic variation can help to explain how variation of plant secondary metabolites is controlled. We explored the transcriptomes of the 3 cardinal terpene chemotypes of Melaleuca alternifolia in young leaves, mature leaves, and stem and compared transcript abundance to variation in the constitutive profile of terpenes. Leaves from chemotype 1 plants (dominated by terpinen-4-ol) show a similar pattern of gene expression when compared to chemotype 5 plants (dominated by 1,8-cineole). Only terpene synthases in young leaves were differentially expressed between these chemotypes, supporting the idea that terpenes are mainly synthetized in young tissue. Chemotype 2 plants (dominated by terpinolene) show a greater degree of differential gene expression compared to the other chemotypes, which might be related to the isolation of plant populations that exhibit this chemotype and the possibility that the terpinolene synthase gene in M. alternifolia was derived by introgression from a closely related species, Melaleuca trichostachya. By using multivariate analyses, we were able to associate terpenes with candidate terpene synthases.

Transcriptome and Biochemical Analysis of a Flower Color Polymorphism in Silene littorea (Caryophyllaceae)

Flower color polymorphisms are widely used as model traits from genetics to ecology, yet determining the biochemical and molecular basis can be challenging. Anthocyanin-based flower color variations can be caused by at least 12 structural and three regulatory genes in the anthocyanin biosynthetic pathway (ABP). We use mRNA-Seq to simultaneously sequence and estimate expression of these candidate genes in nine samples of Silene littorea representing three color morphs (dark pink, light pink and white) across three developmental stages in hopes of identifying the cause of flower color variation. We identified 29 putative paralogs for the 15 candidate genes in the ABP. We assembled complete coding sequences for 16 structural loci and nine of ten regulatory loci. Among these 29 putative paralogs, we identified 622 SNPs, yet only nine synonymous SNPs in Ans had allele frequencies that differentiated pigmented petals (dark pink and light pink) from white petals. These Ans allele frequency differences were further investigated with an expanded sequencing survey of 38 individuals, yet no SNPs consistently differentiated the color morphs. We also found one locus, F3h1, with strong differential expression between pigmented and white samples (>42x). This may be caused by decreased expression of Myb1a in white petal buds. Myb1a in S. littorea is a regulatory locus closely related to Subgroup 7 Mybs known to regulate F3h and other loci in the first half of the ABP in model species. We then compare the mRNA-Seq results with petal biochemistry which revealed cyanidin as the primary anthocyanin and five flavonoid intermediates. Concentrations of three of the flavonoid intermediates were significantly lower in white petals than in pigmented petals (rutin, quercetin and isovitexin). The biochemistry results for rutin, quercetin, luteolin and apigenin are consistent with the transcriptome results suggesting a blockage at F3h, possibly caused by downregulation of Myb1a.

The Phenotypic and Genetic Underpinnings of Flower Size in Polemoniaceae

Corolla length is a labile flower feature and has strong implications for pollinator success. However, the phenotypic and genetic bases of corolla elongation are not well known, largely due to a lack of good candidate genes for potential genetic exploration and functional work. We investigate both the cellular phenotypic differences in corolla length, as well as the genetic control of this trait, in Saltugilia (Polemoniaceae). Taxa in this clade exhibit a large range of flower sizes and differ dramatically in pollinator guilds. Flowers of each species were collected from multiple individuals during four stages of flower development to ascertain if cell number or cell size is more important in determining flower size. In Saltugilia, increased flower size during development appears to be driven more by cell size than cell number. Differences in flower size between species are governed by both cell size and cell number, with the large-flowered S. splendens subsp. grantii having nearly twice as many cells as the small-flowered species. Fully mature flowers of all taxa contain jigsaw cells similar to cells seen in sepals and leaves; however, these cells are not typically found in the developing flowers of most species. The proportion of this cell type in mature flowers appears to have substantial implications, comprising 17-68% of the overall flower size. To identify candidate genes responsible for differences in cell area and cell type, transcriptomes were generated for two individuals of the species with the smallest (S. australis) and largest (S. splendens subsp. grantii) flowers across the same four developmental stages visualized with confocal microscopy. Analyses identified genes associated with cell wall formation that are up-regulated in the mature flower stage compared to mid-stage flowers (75% of mature size). This developmental change is associated with the origin of jigsaw cells in the corolla tube of mature flowers. Further comparisons between mature flowers in the two species revealed 354 transcripts that are up-regulated in the large-flowered S. splendens subsp. grantii compared to the small-flowered S. australis. These results are likely broadly applicable to Polemoniaceae, a clade of nearly 400 species, with extensive variation in floral form and shape.