Pseudogenization and Resurrection of a Speciation Gene

Complete chloroplast genomes and phylogenetic relationships of Pedicularis chinensis and Pedicularis kansuensis

The complete cp genomes of Pedicularis chinensis (GenBank accession number: OQ587614) and Pedicularis kansuensis (GenBank accession number: OQ587613) were sequenced, assembled, and annotated. Their chloroplast (cp) genome lengths were 146,452 bp, and 146,852 bp, respectively; 120 and 116 genes were identified, comprising 75 and 72 protein-coding genes (PCGs), 37 and 36 transfer RNA (tRNA) genes, and 8 and 8 ribosomal RNA (rRNA) genes, for P. chinensis and P. kansuensis, respectively. A simple sequence repeat (SSR) analysis revealed that the repetitive sequences were mainly composed of mononucleotide repeats (A/T motif) and dinucleotide repeats (AT/TA motif). Comparative genomics identified several variant genes (rpl22, rps19, rpl12, ycf1, trnH, psbA, and ndhH) and variant regions (trnS-GGA, trnV-UAC, ndhJ-trnV, ycf4-cemA, ndhE-nhdG, and rpl32-trnL) with a high Pi, indicating the potential to serve as deoxyribo nucleic acid (DNA) barcodes for Pedicularis species identification. The results show that the cp genomes of P. chinensis and P. kansuensis were the same as those of other plants in Pedicularis, with different degrees of AT preference for codons. Large differences in the number of SSRs and the expansion of the inverted repeat (IR) region showed strong variability and interspecific differentiation between these two species and other species represented in the genus Pedicularis. A phylogenetic analysis showed that P. kansuensis had the closest relationship with P. oliveriana, and P. chinensis had the closest relationship with P. aschistorhyncha. These results will facilitate the study of the phylogenetic classification and interspecific evolution of Pedicularis plants.

The physical and evolutionary energy landscapes of devolved protein sequences corresponding to pseudogenes

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners.

The road less taken: Dihydroflavonol 4-reductase inactivation and delphinidin anthocyanin loss underpins a natural intraspecific flower colour variation

Visual cues are of critical importance for the attraction of animal pollinators, however, little is known about the molecular mechanisms underpinning intraspecific floral colour variation. Here, we combined comparative spectral analysis, targeted metabolite profiling, multi-tissue transcriptomics, differential gene expression, sequence analysis and functional analysis to investigate a bee-pollinated orchid species, Glossodia major with common purple- and infrequent white-flowered morphs. We found uncommon and previously unreported delphinidin-based anthocyanins responsible for the conspicuous and pollinator-perceivable colour of the purple morph and three genetic changes underpinning the loss of colour in the white morph - (1) a loss-of-function (LOF; frameshift) mutation affecting dihydroflavonol 4-reductase (DFR1) coding sequence due to a unique 4-bp insertion, (2) specific downregulation of functional DFR1 expression and (3) the unexpected discovery of chimeric Gypsy transposable element (TE)-gene (DFR) transcripts with potential consequences to the genomic stability and post-transcriptional or epigenetic regulation of DFR. This is one of few known cases where regulatory changes and LOF mutation in an anthocyanin structural gene, rather than transcription factors, are important. Furthermore, if TEs prove to be a frequent source of mutation, the interplay between environmental stress-induced TE evolution and pollinator-mediated selection for adaptive colour variation may be an overlooked mechanism maintaining floral colour polymorphism in nature.

Shades of white: The Petunia long corolla tube clade evolutionary history

Delimiting species is challenging in recently diverged species, and adaptive radiation is fundamental to understanding the evolutionary processes because it requires multiple ecological opportunities associated with adaptation to biotic and abiotic environments. The young Petunia genus (Solanaceae) is an excellent opportunity to study speciation because of its association with pollinators and unique microenvironments. This study evaluated the phylogenetic relationships among a Petunia clade species with different floral syndromes that inhabit several environments. We based our work on multiple individuals per lineage and employed nuclear and plastid phylogenetic markers and nuclear microsatellites. The phylogenetic tree revealed two main groups regarding the elevation of the distribution range, whereas microsatellites showed high polymorphism-sharing splitting lineages into three clusters. Isolation by distance, migration followed by new environment colonization, and shifts in floral syndrome were the motors for lineage differentiation, including infraspecific structuring, which suggests the need for taxonomic revision in the genus.

Cell layer-specific expression of the homeotic MADS-box transcription factor PhDEF contributes to modular petal morphogenesis in petunia

Floral homeotic MADS-box transcription factors ensure the correct morphogenesis of floral organs, which are organized in different cell layers deriving from distinct meristematic layers. How cells from these distinct layers acquire their respective identities and coordinate their growth to ensure normal floral organ morphogenesis is unresolved. Here, we studied petunia (Petunia × hybrida) petals that form a limb and tube through congenital fusion. We identified petunia mutants (periclinal chimeras) expressing the B-class MADS-box gene DEFICIENS in the petal epidermis or in the petal mesophyll, called wico and star, respectively. Strikingly, wico flowers form a strongly reduced tube while their limbs are almost normal, while star flowers form a normal tube but greatly reduced and unpigmented limbs, showing that petunia petal morphogenesis is highly modular. These mutants highlight the layer-specific roles of PhDEF during petal development. We explored the link between PhDEF and petal pigmentation, a well-characterized limb epidermal trait. The anthocyanin biosynthesis pathway was strongly downregulated in star petals, including its major regulator ANTHOCYANIN2 (AN2). We established that PhDEF directly binds to the AN2 terminator in vitro and in vivo, suggesting that PhDEF might regulate AN2 expression and therefore petal epidermis pigmentation. Altogether, we show that cell layer-specific homeotic activity in petunia petals differently impacts tube and limb development, revealing the relative importance of the different cell layers in the modular architecture of petunia petals.

Genome evolution in plants and the origins of innovation

Plant evolution has been characterised by a series of major novelties in their vegetative and reproductive traits that have led to greater complexity. Underpinning this diversification has been the evolution of the genome. When viewed at the scale of the plant kingdom, plant genome evolution has been punctuated by conspicuous instances of gene and whole-genome duplication, horizontal gene transfer and extensive gene loss. The periods of dynamic genome evolution often coincide with the evolution of key traits, demonstrating the coevolution of plant genomes and phenotypes at a macroevolutionary scale. Conventionally, plant complexity and diversity have been considered through the lens of gene duplication and the role of gene loss in plant evolution remains comparatively unexplored. However, in light of reductive evolution across multiple plant lineages, the association between gene loss and plant phenotypic diversity warrants greater attention.

Reversibility of digit loss revisited: Limb diversification in Bachia lizards (gymnophthalmidae)

Strict interpretations of the Dollo's Law lead to postulation that trait loss is irreversible and organisms never recover ancestral phenotypes. Dollo, however, admitted the possibility of reversals in trait loss when predicted differences between reversed (derived) and ancestral forms. Phenotypic signatures from reversals are expected, as the historical context of a reversal in trait loss differs from the initial setting where the trait originally evolved. This article combines morphological and molecular information for Bachia scolecoides to discuss phenotypic and genetic patterns established during processes that reversed digit loss in Gymnophthalmidae (also termed microteiid lizards). Results suggest that pathways leading to the derived tetradactyl state of B. scolecoides comprise particularities in their origin and associated processes. Autopodial bones of B. scolecoides lack digit identity, and muscle anatomy is very similar between manus and pes. Gymnophthalmidae sequence patterns in the limb-specific sonic hedgehog enhancer (ZRS) suggest that regulation of shh expression did not degenerate in Bachia, given the prediction of similar motifs despite mutations specific to Bachia. Persistence of developmental mechanisms might explain intermittent character expression leading to reversals of digit loss, as ZRS signaling pathways remain active during the development of at least one pair of appendices in Bachia, especially if some precursors persisted at early stages. Patterns of ZRS sequences suggest that irreversibility of trait loss might be lineage-specific (restricted to Gymnophthalmini) and contingent to the type of signature established. These results provide insights regarding possible mechanisms that may allow reactivation of developmental programs in specific regions of the embryo.

Unravelling the due importance of pseudogenes and their resurrection in plants

The complexities of a genome are underpinned to the vast expanses of the intergenic region, which constitutes ∼97-98% of the genome. This region is essentially composed of what is colloquially referred to as the "junk DNA" and is composed of various elements like transposons, repeats, pseudogenes, etc. The latter have long been considered as dead elements merely contributing to transcriptional noise in the genome. Many studies now describe the previously unknown regulatory functions of these genes. Recent advances in the Next-generation sequencing (NGS) technologies have allowed unprecedented access to these regions. With the availability of whole genome sequences of more than 788 different plant species in past 20 years, genome annotation has become feasible like never before. Different bioinformatic pipelines are available for the identification of pseudogenes. However, still little is known about their biological functions. The functional validation of these genes remains challenging and research in this area is still in infancy, particularly in plants. CRISPR/Cas-based genome editing could provide solutions to understand the biological roles of these genes by allowing creation of precise edits within these genes. The possibility of pseudogene reactivation or resurrection as has been demonstrated in a few studies might open new avenues of genetic manipulation to yield a desirable phenotype. This review aims at comprehensively summarizing the progress made with regards to the identification of pseudogenes and understanding their biological functions in plants.

New insights into the genetic manipulation of the R2R3-MYB and CHI gene families on anthocyanin pigmentation in Petunia hybrida

Several R2R3-MYB genes control anthocyanin pigmentation in petunia, and ANTHOCYANIN-2 (AN2) is treated as the main player in petal limbs. However, the actual roles of R2R3-MYBs in the coloration of different floral tissues in the so called "darkly-veined" petunias are still not clear. The genetic background and expression of AN2 paralogs from various petunias with different color patterns were identified. All "darkly-veined" genotypes have the identical mutation in the AN2 gene, but express a different functional paralog - ANTHOCYANIN-4 (AN4) - abundantly in flowers. Constitutive overexpression of PhAN4 in this petunia resulted not only in a fully colored flower but also in a clearly visible pigmentation in the green tissue and roots, which can be rapidly increased by stress conditions. Suppression of AN4 gene resulted in discolored petals and whitish anthers. Interestingly, when a similar white flower phenotype was achieved by knockout of an essential structural gene of anthocyanin biosynthesis - CHALCONE ISOMERASE-A (CHI-A) - the plant responded directly by upregulating of another paralogs - DEEP PURPLE (DPL) and PURPLE HAZE (PHZ). Moreover, we also found that CHI-B can partially substitute for CHI-A in anthers, but not in vegetative tissues. Further, no significant effects on the longevity of white or enhanced colored flowers were observed compared with the wild type. We concluded that endogenous up-regulation of AN4 leads to the restoration of petal color in the "darkly-veined" phenotypes as a result of the breeding process under human selection, and CHI-B is a backup for CHI-A acitvity in some floral tissues.

Towards understanding paleoclimate impacts on primate de novo genes

De novo genes are genes that emerge as new genes in some species, such as primate de novo genes that emerge in certain primate species. Over the past decade, a great deal of research has been conducted regarding their emergence, origins, functions, and various attributes in different species, some of which have involved estimating the ages of de novo genes. However, limited by the number of species available for whole-genome sequencing, relatively few studies have focused specifically on the emergence time of primate de novo genes. Among those, even fewer investigate the association between primate gene emergence with environmental factors, such as paleoclimate (ancient climate) conditions. This study investigates the relationship between paleoclimate and human gene emergence at primate species divergence. Based on 32 available primate genome sequences, this study has revealed possible associations between temperature changes and the emergence of de novo primate genes. Overall, findings in this study are that de novo genes tended to emerge in the recent 13 MY when the temperature continues cooling, which is consistent with past findings. Furthermore, in the context of an overall trend of cooling temperature, new primate genes were more likely to emerge during local warming periods, where the warm temperature more closely resembled the environmental condition that preceded the cooling trend. Results also indicate that both primate de novo genes and human cancer-associated genes have later origins in comparison to random human genes. Future studies can be in-depth on understanding human de novo gene emergence from an environmental perspective as well as understanding species divergence from a gene emergence perspective.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

Genetic architecture of a pollinator shift and its fate in secondary hybrid zones of two Petunia species

BACKGROUND

Theory suggests that the genetic architecture of traits under divergent natural selection influences how easily reproductive barriers evolve and are maintained between species. Divergently selected traits with a simple genetic architecture (few loci with major phenotypic effects) should facilitate the establishment and maintenance of reproductive isolation between species that are still connected by some gene flow. While empirical support for this idea appears to be mixed, most studies test the influence of trait architectures on reproductive isolation only indirectly. Petunia plant species are, in part, reproductively isolated by their different pollinators. To investigate the genetic causes and consequences of this ecological isolation, we deciphered the genetic architecture of three floral pollination syndrome traits in naturally occurring hybrids between the widespread Petunia axillaris and the highly endemic and endangered P. exserta.

RESULTS

Using population genetics, Bayesian linear mixed modelling and genome-wide association studies, we found that the three pollination syndrome traits vary in genetic architecture. Few genome regions explain a majority of the variation in flavonol content (defining UV floral colour) and strongly predict the trait value in hybrids irrespective of interspecific admixture in the rest of their genomes. In contrast, variation in pistil exsertion and anthocyanin content (defining visible floral colour) is controlled by many genome-wide loci. Opposite to flavonol content, the genome-wide proportion of admixture between the two species predicts trait values in their hybrids. Finally, the genome regions strongly associated with the traits do not show extreme divergence between individuals representing the two species, suggesting that divergent selection on these genome regions is relatively weak within their contact zones.

CONCLUSIONS

Among the traits analysed, those with a more complex genetic architecture are best maintained in association with the species upon their secondary contact. We propose that this maintained genotype-phenotype association is a coincidental consequence of the complex genetic architectures of these traits: some of their many underlying small-effect loci are likely to be coincidentally linked with the actual barrier loci keeping these species partially isolated upon secondary contact. Hence, the genetic architecture of a trait seems to matter for the outcome of hybridization not only then when the trait itself is under selection.

Tight genetic linkage of genes causing hybrid necrosis and pollinator isolation between young species

The mechanisms of reproductive isolation that cause phenotypic diversification and eventually speciation are a major topic of evolutionary research. Hybrid necrosis is a post-zygotic isolation mechanism in which cell death develops in the absence of pathogens. It is often due to the incompatibility between proteins from two parents. Here we describe a unique case of hybrid necrosis due to an incompatibility between loci on chromosomes 2 and 7 between two pollinator-isolated Petunia species. Typical immune responses as well as endoplasmic reticulum stress responses are induced in the necrotic line. The locus on chromosome 2 encodes ChiA1, a bifunctional GH18 chitinase/lysozyme. The enzymatic activity of ChiA1 is dispensable for the development of necrosis. We propose that the extremely high expression of ChiA1 involves a positive feedback loop between the loci on chromosomes 2 and 7. ChiA1 is tightly linked to major genes involved in the adaptation to different pollinators, a form of pre-zygotic isolation. This linkage of pre- and post-zygotic barriers strengthens reproductive isolation and probably contributes to rapid diversification and speciation.

Physiological changes besides the enhancement of pigmentation in Petunia hybrida caused by overexpression of PhAN2, an R2R3-MYB transcription factor

Ectopic expression of PhAN2 in vegetative tissue can improve regeneration and adventitious rooting but inhibit axillary bud outgrowth of petunia, while overexpression specifically in flowers could shorten longevity. Anthocyanin 2 has been only treated as a critical positive regulation factor of anthocyanin biosynthesis in petunia flowers. To determine if this gene had other functions in plant growth, we overexpressed this gene in an an2 mutant petunia cultivar driven by promoters with different strengths or tissue specificity. Various physiological processes of transformants in different growth stages and environments were analyzed. Besides the expected pigmentation improvement in different tissues, the results also showed that ectopic expression of AN2 could improve the regeneration skill but inhibit the axillary bud germination of in vitro plants. Moreover, the rooting ability of shoot tips of transformants was significantly improved, while some transgenic lines' flower longevity was shortened. Gene expression analysis showed that the transcripts level of AN2, partner genes anthocyanin 1 (AN1), anthocyanin 11 (AN11), and target gene dihydroflavonol 4-reductase (DFR) was altered in the different transgenic lines. In addition, ethylene biosynthesis-related genes 1-aminocyclopropane-1-carboxylic acid synthase (ACS1) and ACC oxidase (ACO1) were upregulated in rooting and flower senescence processes but at different time points. Overall, our data demonstrate that the critical role of this AN2 gene in plant growth physiology may extend beyond that of a single activator of anthocyanin biosynthesis.

The origin and loss of interferon regulatory factor 10 (IRF10) in different lineages of vertebrates

The vertebrate IFN regulatory factor (IRF) family consists of 11 members that exert distinct roles in a variety of biological processes, including antiviral defense, regulation of cell proliferation, differentiation and apoptosis. Of these, IRF10 is widely present in different vertebrate lineages, but appears to have been lost in primates and rodents. To understand the evolutionary occurrence of IRF10, we performed comparative analyses of currently available genomic data in a taxonomically diverse set of vertebrates, and found that IRF10 originated after the divergence of chondrichthyans from gnathostomes. Phylogenetically, vertebrate IRF10 is much more closely related to IRF4 than to IRF8 or IRF9, although these four IRFs may have a common ancestor. In addition, the loss of IRF10 in Euarchontoglires might be resulted from mutation accumulation, and the rate of mutations in rodents appears to be higher than in the primate lineage. In primates, the gene-disruptive mutations may have occurred at least prior to the separation of new world monkey and old world primates, roughly 40 million years ago. Overall, we propose a detailed evolutionary scenario for IRF10, which may help us understand the evolutionary mechanisms in the expansion and contraction of the IRF family.

Flower Development in the Solanaceae

Flower development is the process leading from a reproductive meristem to a mature flower with fully developed floral organs. This multi-step process is complex and involves thousands of genes in intertwined regulatory pathways; navigating through the FLOR-ID website will give an impression of this complexity and of the astonishing amount of work that has been carried on the topic (Bouché et al., Nucleic Acids Res 44:D1167-D1171, 2016). Our understanding of flower development mostly comes from the model species Arabidopsis thaliana, but numerous other studies outside of Brassicaceae have helped apprehend the conservation of these mechanisms in a large evolutionary context (Moyroud and Glover, Curr Biol 27:R941-R951, 2017; Smyth, New Phytol 220:70-86, 2018; Soltis et al., Ann Bot 100:155-163, 2007). Integrating additional species and families to the research on this topic can only advance our understanding of flower development and its evolution.In this chapter, we review the contribution that the Solanaceae family has made to the comprehension of flower development. While many of the general features of flower development (i.e., the key molecular players involved in flower meristem identity, inflorescence architecture or floral organ development) are similar to Arabidopsis, our main objective in this chapter is to highlight the points of divergence and emphasize specificities of the Solanaceae. We will not discuss the large topics of flowering time regulation, inflorescence architecture and fruit development, and we will restrict ourselves to the mechanisms included in a time window after the floral transition and before the fertilization. Moreover, this review will not be exhaustive of the large amount of work carried on the topic, and the choices that we made to describe in large details some stories from the literature are based on the soundness of the functional work performed, and surely as well on our own preferences and expertise.First, we will give a brief overview of the Solanaceae family and some of its specificities. Then, our focus will be on the molecular mechanisms controlling floral organ identity, for which extended functional work in petunia led to substantial revisions to the famous ABC model. Finally, after reviewing some studies on floral organ initiation and growth, we will discuss floral organ maturation, using the examples of the inflated calyx of the Chinese lantern Physalis and petunia petal pigmentation.

Single gene mutation in a plant MYB transcription factor causes a major shift in pollinator preference

Understanding the molecular basis of reproductive isolation and speciation is a key goal of evolutionary genetics. In the South American genus Petunia, the R2R3-MYB transcription factor MYB-FL regulates the biosynthesis of UV-absorbing flavonol pigments, a major determinant of pollinator preference. MYB-FL is highly expressed in the hawkmoth-pollinated P. axillaris, but independent losses of its activity in sister taxa P. secreta and P. exserta led to UV-reflective flowers and associated pollinator shifts in each lineage (bees and hummingbirds, respectively). We created a myb-fl CRISPR mutant in P. axillaris and studied the effect of this single gene on innate pollinator preference. The mutation strongly reduced the expression of the two key flavonol-related biosynthetic genes but only affected the expression of few other genes. The mutant flowers were UV reflective as expected but additionally contained low levels of visible anthocyanin pigments. Hawkmoths strongly preferred the wild-type P. axillaris over the myb-fl mutant, whereas both social and solitary bee preference depended on the level of visible color of the mutants. MYB-FL, with its specific expression pattern, small number of target genes, and key position at the nexus of flavonol and anthocyanin biosynthetic pathways, provides a striking example of evolution by single mutations of large phenotypic effect.

Phylogenetic Analyses of Some Key Genes Provide Information on Pollinator Attraction in Solanaceae

Floral syndromes are known by the conserved morphological traits in flowers associated with pollinator attraction, such as corolla shape and color, aroma emission and composition, and rewards, especially the nectar volume and sugar concentration. Here, we employed a phylogenetic approach to investigate sequences of genes enrolled in the biosynthetic pathways responsible for some phenotypes that are attractive to pollinators in Solanaceae genomes. We included genes involved in visible color, UV-light response, scent emission, and nectar production to test the hypothesis that these essential genes have evolved by convergence under pollinator selection. Our results refuted this hypothesis as all four studied genes recovered the species' phylogenetic relationships, even though some sites were positively selected. We found differences in protein motifs among genera in Solanaceae that were not necessarily associated with the same floral syndrome. Although it has had a crucial role in plant diversification, the plant-pollinator interaction is complex and still needs further investigation, with genes evolving not only under the influence of pollinators, but by the sum of several evolutionary forces along the speciation process in Solanaceae.

Plant phylogenetics: The never-ending cycle of evolutionary gains and losses

The Zygnematophyceae is the sister clade to the land plants, but their biology remains mysterious. In a new study, a resolved phylogeny and a scenario for the evolution of multicellularity in that clade are proposed.

Uncovering anthocyanin diversity in potato landraces (Solanum tuberosum L. Phureja) using RNA-seq

Potato (Solanum tuberosum L.) is the third largest source of antioxidants in the human diet, after maize and tomato. Potato landraces have particularly diverse contents of antioxidant compounds such as anthocyanins. We used this diversity to study the evolutionary and genetic basis of anthocyanin pigmentation. Specifically, we analyzed the transcriptomes and anthocyanin content of tubers from 37 landraces with different colorations. We conducted analyses of differential expression between potatoes with different colorations and used weighted correlation network analysis to identify genes whose expression is correlated to anthocyanin content across landraces. A very significant fraction of the genes identified in these two analyses had annotations related to the flavonoid-anthocyanin biosynthetic pathway, including 18 enzymes and 5 transcription factors. Importantly, the causal genes at the D, P and R loci governing anthocyanin accumulation in potato cultivars also showed correlations to anthocyanin production in the landraces studied here. Furthermore, we found that 60% of the genes identified in our study were located within anthocyanin QTLs. Finally, we identified new candidate enzymes and transcription factors that could have driven the diversification of anthocyanins. Our results indicate that many anthocyanins biosynthetic genes were manipulated in ancestral potato breeding and can be used in future breeding programs.

Floral color is not as simple as it once seemed

A false start mutation produces reduced protein and flower color, highlighting the role of mutations affecting protein translation in phenotypic evolution and variation.

Lost in translation: Molecular basis of reduced flower coloration in a self-pollinated monkeyflower (Mimulus) species

Phenotypic evolution is usually attributed to changes in protein function or gene transcription. In principle, mutations that affect protein abundance through enhancing or attenuating protein translation also could be an important source for phenotypic evolution. However, these types of mutations remain largely unexplored in the studies of phenotypic variation in nature. Through fine-scale genetic mapping and functional interrogation, we identify a single nucleotide substitution in an anthocyanin-activating R2R3-MYB gene causing flower color variation between a pair of closely related monkeyflower (Mimulus) species, the hummingbird-pollinated Mimulus cardinalis, and self-pollinated Mimulus parishii. This causal mutation is located in the 5' untranslated region and generates an upstream ATG start codon, leading to attenuated protein translation and reduced flower coloration in the self-pollinated species. Together, our results provide empirical support for the role of mutations affecting protein translation, as opposed to protein function or transcript level, in natural phenotypic variation.

The mungbean VrP locus encoding MYB90, an R2R3-type MYB protein, regulates anthocyanin biosynthesis

Anthocyanins are water-soluble pigments present in several tissues/parts of plants. The pigments provide color and are wildly known for health benefits for human, insect attraction for plant pollination, and stress resistance in plants. Anthocyanin content variations in mungbean [Vigna radiata (L.) Wilczek] were first noticed a long time ago, but the genetic mechanism controlling the anthocyanins in mungbean remains unknown. An F₂ population derived from the cross between purple-hypocotyl (V2709) and green-hypocotyl (Sulv1) mungbeans was used to map the VrP locus controlling purple hypocotyl. The VrP locus was mapped to a 78.9-kb region on chromosome 4. Sequence comparison and gene expression analysis identified an R2R3-MYB gene VrMYB90 as the candidate gene for the VrP locus. Haplotype analysis using 124 mungbean accessions suggested that 10 single nucleotide polymorphisms (SNPs) in exon 3 may lead to an abolished expression of VrMYB90 and an absence of anthocyanin accumulation in the hypocotyl of Sulv1 and KPS2. The overexpression of VrMYB90 in mungbean hairy root, tobacco leaf, and Arabidopsis resulted in anthocyanin accumulation (purple color). Gene expression analysis demonstrated that VrMYB90 regulated anthocyanin accumulation in the hypocotyl, stem, petiole, and flowers, and the expression was sensitive to light. VrMYB90 protein may upregulate VrDFR encoding dihydroflavonol 4-reductase at the late biosynthesis step of anthocyanins in mungbeans. These results suggest that VrMYB90 is the dominator in the spatiotemporal regulation of anthocyanin biosynthesis. Our results provide insight into the biosynthesis mechanism of anthocyanin and a theoretical basis for breeding mungbeans.

Macroevolutionary Analyses Provide New Evidence of Phasmid Wings Evolution as a Reversible Process

The concept that complex ancestral traits can never be recovered after their loss is still widely accepted, despite phylogenetic and molecular approaches suggest instances where phenotypes may have been lost throughout the evolutionary history of a clade and subsequently reverted back in derived lineages. One of the first and most notable examples of such a process is wing evolution in phasmids; this polyneopteran order of insects, which comprises stick and leaf insects, has played a central role in initiating a long-standing debate on the topic. In this study, a novel and comprehensive time tree including over 300 Phasmatodea species is used as a framework for investigating wing evolutionary patterns in the clade. Despite accounting for several possible biases and sources of uncertainty, macroevolutionary analyses consistently revealed multiple reversals to winged states taking place after their loss, and reversibility is coupled with higher species diversification rates. Our findings support a loss of or reduction in wings that occurred in the lineage leading to the extant phasmid most recent common ancestor, and brachyptery is inferred to be an unstable state unless co-opted for nonaerodynamic adaptations. We also explored how different assumptions of wing reversals probability could impact their inference: we found that until reversals are assumed to be over 30 times more unlikely than losses, they are consistently inferred despite uncertainty in tree and model parameters. Our findings demonstrate that wing evolution is a reversible and dynamic process in phasmids and contribute to our understanding of complex trait evolution.

Transcription factors evolve faster than their structural gene targets in the flavonoid pigment pathway

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

Pairs of mutually compensatory frameshifting mutations contribute to protein evolution

Insertions and deletions of lengths not divisible by 3 in protein-coding sequences cause frameshifts that usually induce premature stop codons and may carry a high fitness cost. However, this cost can be partially offset by a second compensatory indel restoring the reading frame. The role of such pairs of compensatory frameshifting mutations (pCFMs) in evolution has not been studied systematically. Here, we use whole-genome alignments of protein-coding genes of 100 vertebrate species, and of 122 insect species, studying the prevalence of pCFMs in their divergence. We detect a total of 624 candidate pCFM genes; six of them pass stringent quality filtering, including three human genes: RAB36, ARHGAP6, and NCR3LG1. In some instances, amino acid substitutions closely predating or following pCFMs restored the biochemical similarity of the frameshifted segment to the ancestral amino acid sequence, possibly reducing or negating the fitness cost of the pCFM. Typically, however, the biochemical similarity of the frameshifted sequence to the ancestral one was not higher than the similarity of a random sequence of a protein-coding gene to its frameshifted version, indicating that pCFMs can uncover radically novel regions of protein space. In total, pCFMs represent an appreciable and previously overlooked source of novel variation in amino acid sequences.

Noninsect-Based Diet Leads to Structural and Functional Changes of Acidic Chitinase in Carnivora

Acidic chitinase (Chia) digests the chitin of insects in the omnivorous stomach and the chitinase activity in carnivorous Chia is significantly lower than that of the omnivorous enzyme. However, mechanistic and evolutionary insights into the functional changes in Chia remain unclear. Here we show that a noninsect-based diet has caused structural and functional changes in Chia during the course of evolution in Carnivora. By creating mouse-dog chimeric Chia proteins and modifying the amino acid sequences, we revealed that F214L and A216G substitutions led to the dog enzyme activation. In 31 Carnivora, Chia was present as a pseudogene with stop codons in the open reading frame (ORF) region. Importantly, the Chia proteins of skunk, meerkat, mongoose, and hyena, which are insect-eating species, showed high chitinolytic activity. The cat Chia pseudogene product was still inactive even after ORF restoration. However, the enzyme was activated by matching the number and position of Cys residues to an active form and by introducing five meerkat Chia residues. Mutations affecting the Chia conformation and activity after pseudogenization have accumulated in the common ancestor of Felidae due to functional constraints. Evolutionary analysis indicates that Chia genes are under relaxed selective constraint in species with noninsect-based diets except for Canidae. These results suggest that there are two types of inactivating processes in Carnivora and that dietary changes affect the structure and activity of Chia.

Incipient diploidization of the medicinal plant Perilla within 10,000 years

Perilla is a young allotetraploid Lamiaceae species widely used in East Asia as herb and oil plant. Here, we report the high-quality, chromosome-scale genomes of the tetraploid (Perilla frutescens) and the AA diploid progenitor (Perilla citriodora). Comparative analyses suggest post Neolithic allotetraploidization within 10,000 years, and nucleotide mutation in tetraploid is 10% more than in diploid, both of which are dominated by G:C → A:T transitions. Incipient diploidization is characterized by balanced swaps of homeologous segments, and subsequent homeologous exchanges are enriched towards telomeres, with excess of replacements of AA genes by fractionated BB homeologs. Population analyses suggest that the crispa lines are close to the nascent tetraploid, and involvement of acyl-CoA: lysophosphatidylcholine acyltransferase gene for high α-linolenic acid content of seed oil is revealed by GWAS. These resources and findings provide insights into incipient diploidization and basis for breeding improvement of this medicinal plant.

Repeated horizontal gene transfer of GALactose metabolism genes violates Dollo's law of irreversible loss

Dollo's law posits that evolutionary losses are irreversible, thereby narrowing the potential paths of evolutionary change. While phenotypic reversals to ancestral states have been observed, little is known about their underlying genetic causes. The genomes of budding yeasts have been shaped by extensive reductive evolution, such as reduced genome sizes and the losses of metabolic capabilities. However, the extent and mechanisms of trait reacquisition after gene loss in yeasts have not been thoroughly studied. Here, through phylogenomic analyses, we reconstructed the evolutionary history of the yeast galactose utilization pathway and observed widespread and repeated losses of the ability to utilize galactose, which occurred concurrently with the losses of GALactose (GAL) utilization genes. Unexpectedly, we detected multiple galactose-utilizing lineages that were deeply embedded within clades that underwent ancient losses of galactose utilization. We show that at least two, and possibly three, lineages reacquired the GAL pathway via yeast-to-yeast horizontal gene transfer. Our results show how trait reacquisition can occur tens of millions of years after an initial loss via horizontal gene transfer from distant relatives. These findings demonstrate that the losses of complex traits and even whole pathways are not always evolutionary dead-ends, highlighting how reversals to ancestral states can occur.

Complex evolution of novel red floral color in Petunia

Red flower color has arisen multiple times and is generally associated with hummingbird pollination. The majority of evolutionary transitions to red color proceeded from purple lineages and tend to be genetically simple, almost always involving a few loss-of-function mutations of major phenotypic effect. Here we report on the complex evolution of a novel red floral color in the hummingbird-pollinated Petunia exserta (Solanaceae) from a colorless ancestor. The presence of a red color is remarkable because the genus cannot synthesize red anthocyanins and P. exserta retains a nonfunctional copy of the key MYB transcription factor AN2. We show that moderate upregulation and a shift in tissue specificity of an AN2 paralog, DEEP PURPLE, restores anthocyanin biosynthesis in P. exserta. An essential shift in anthocyanin hydroxylation occurred through rebalancing the expression of three hydroxylating genes. Furthermore, the downregulation of an acyltransferase promotes reddish hues in typically purple pigments by preventing acyl group decoration of anthocyanins. This study presents a rare case of a genetically complex evolutionary transition toward the gain of a novel red color.

Discrete bHLH transcription factors play functionally overlapping roles in pigmentation patterning in flowers of Antirrhinum majus

Floral pigmentation patterning is important for pollinator attraction as well as aesthetic appeal. Patterning of anthocyanin accumulation is frequently associated with variation in activity of the Myb, bHLH and WDR transcription factor complex (MBW) that regulates anthocyanin biosynthesis. Investigation of two classic mutants in Antirrhinum majus, mutabilis and incolorata I, showed they affect a gene encoding a bHLH protein belonging to subclade bHLH-2. The previously characterised gene, Delila, which encodes a bHLH-1 protein, has a bicoloured mutant phenotype, with residual lobe-specific pigmentation conferred by Incolorata I. Both Incolorata I and Delila induce expression of the anthocyanin biosynthetic gene DFR. Rosea 1 (Myb) and WDR1 proteins compete for interaction with Delila, but interact positively to promote Incolorata I activity. Delila positively regulates Incolorata I and WDR1 expression. Hierarchical regulation can explain the bicoloured patterning of delila mutants, through effects on both regulatory gene expression and the activity of promoters of biosynthetic genes like DFR that mediate MBW regulation. bHLH-1 and bHLH-2 proteins contribute to establishing patterns of pigment distribution in A. majus flowers in two ways: through functional redundancy in regulating anthocyanin biosynthetic gene expression, and through differences between the proteins in their ability to regulate genes encoding transcription factors.

Phenotypic plasticity explains apparent reverse evolution of fat synthesis in parasitic wasps

Numerous cases of evolutionary trait loss and regain have been reported over the years. Here, we argue that such reverse evolution can also become apparent when trait expression is plastic in response to the environment. We tested this idea for the loss and regain of fat synthesis in parasitic wasps. We first show experimentally that the wasp Leptopilina heterotoma switches lipogenesis on in a fat-poor environment, and completely off in a fat-rich environment. Plasticity suggests that this species did not regain fat synthesis, but that it can be switched off in some environmental settings. We then compared DNA sequence variation and protein domains of several more distantly related parasitoid species thought to have lost lipogenesis, and found no evidence for non-functionality of key lipogenesis genes. This suggests that other parasitoids may also show plasticity of fat synthesis. Last, we used individual-based simulations to show that a switch for plastic expression can remain functional in the genome for thousands of generations, even if it is only used sporadically. The evolution of plasticity could thus also explain other examples of apparent reverse evolution.

How to become a crab: Phenotypic constraints on a recurring body plan

A fundamental question in biology is whether phenotypes can be predicted by ecological or genomic rules. At least five cases of convergent evolution of the crab-like body plan (with a wide and flattened shape, and a bent abdomen) are known in decapod crustaceans, and have, for over 140 years, been known as "carcinization." The repeated loss of this body plan has been identified as "decarcinization." In reviewing the field, we offer phylogenetic strategies to include poorly known groups, and direct evidence from fossils, that will resolve the history of crab evolution and the degree of phenotypic variation within crabs. Proposed ecological advantages of the crab body are summarized into a hypothesis of phenotypic integration suggesting correlated evolution of the carapace shape and abdomen. Our premise provides fertile ground for future studies of the genomic and developmental basis, and the predictability, of the crab-like body form.

Evidence for the loss and recovery of SLAMF9 during human evolution: implications on Dollo's law

Signaling lymphocyte activation molecule family member 9 (SLAMF9) is a cell surface protein of the CD2/SLAM family of leukocyte surface receptors. It is conserved throughout mammals and has roles in the initiation of inflammatory responses and regulation of plasmacytoid dendritic cell function. Through comparison of reference sequences encoding SLAMF9 in human, mouse, and primate sequences, we have determined that the SLAMF9 gene underwent successive mutation events, resulting in the loss of the protein and subsequent recovery of a less stable version. The mutations included a single base pair deletion in the second exon and a change in the splice acceptor site of that same exon. These changes would have had the effect of creating and later repairing a frameshift in the coding sequence. These events took place since the divergence of the human lineage from the chimpanzee-human last common ancestor and represent the first known case of the functional loss and recovery of a gene within the human lineage.

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.

Structure and contingency determine mutational hotspots for flower color evolution

Evolutionary genetic studies have uncovered abundant evidence for genomic hotspots of phenotypic evolution, as well as biased patterns of mutations at those loci. However, the theoretical basis for this concentration of particular types of mutations at particular loci remains largely unexplored. In addition, historical contingency is known to play a major role in evolutionary trajectories, but has not been reconciled with the existence of such hotspots. For example, do the appearance of hotspots and the fixation of different types of mutations at those loci depend on the starting state and/or on the nature and direction of selection? Here, we use a computational approach to examine these questions, focusing the anthocyanin pigmentation pathway, which has been extensively studied in the context of flower color transitions. We investigate two transitions that are common in nature, the transition from blue to purple pigmentation and from purple to red pigmentation. Both sets of simulated transitions occur with a small number of mutations at just four loci and show strikingly similar peaked shapes of evolutionary trajectories, with the mutations of the largest effect occurring early but not first. Nevertheless, the types of mutations (biochemical vs. regulatory) as well as their direction and magnitude are contingent on the particular transition. These simulated color transitions largely mirror findings from natural flower color transitions, which are known to occur via repeated changes at a few hotspot loci. Still, some types of mutations observed in our simulated color evolution are rarely observed in nature, suggesting that pleiotropic effects further limit the trajectories between color phenotypes. Overall, our results indicate that the branching structure of the pathway leads to a predictable concentration of evolutionary change at the hotspot loci, but the types of mutations at these loci and their order is contingent on the evolutionary context.

What could be the fate of secondary contact zones between closely related plant species?

Interspecific hybridization has been fundamental in plant evolution. Nevertheless, the fate of hybrid zones throughout the generations remains poorly addressed. We analyzed a pair of recently diverged, interfertile, and sympatric Petunia species to ask what fate the interspecific hybrid population has met over time. We analyzed the genetic diversity in two generations from two contact sites and evaluated the effect of introgression. To do this, we collected all adult plants from the contact zones, including canonicals and intermediary colored individuals, and compared them with purebred representatives of both species based on seven highly informative microsatellite loci. We compared the genetic diversity observed in the contact zones with what is seen in isolated populations of each species, considering two generations of these annual species. Our results have confirmed the genetic differentiation between the species and the hybrid origin of the majority of the intermediary colored individuals. We also observed a differentiation related to genetic variability and inbreeding levels among the populations. Over time, there were no significant differences per site related to genetic diversity or phenotype composition. We found two stable populations kept by high inbreeding and backcross rates that influence the genetic diversity of their parental species through introgression.

Pollen colour morphs take different paths to fitness

Colour phenotypes are often involved in communication and are thus under selection by species interactions. However, selection may also act on colour through correlated traits or alternative functions of biochemical pigments. Such forms of selection are instrumental in maintaining petal colour diversity in plants. Pollen colour also varies markedly, but the maintenance of this variation is little understood. In Campanula americana, pollen ranges from white to dark purple, with darker morphs garnering more pollinator visits and exhibiting elevated pollen performance under heat stress. Here, we generate an F2 population segregating for pollen colour and measure correlations with floral traits, pollen attributes and plant-level traits related to fitness. We determine the pigment biochemistry of colour variants and evaluate maternal and paternal fitness of light and dark morphs by crossing within and between morphs. Pollen colour was largely uncorrelated with floral traits (petal colour, size, nectar traits) suggesting it can evolve independently. Darker pollen grains were larger and had higher anthocyanin content (cyanidin and peonidin) which may explain why they outperform light pollen under heat stress. Overall, pollen-related fitness metrics were greater for dark pollen, and dark pollen sires generated seeds with higher germination potential. Conversely, light pollen plants produce 61% more flowers than dark, and 18% more seeds per fruit, suggesting a seed production advantage. Results indicate that light and dark morphs may achieve fitness through different means-dark morphs appear to have a pollen advantage whereas light morphs have an ovule advantage-helping to explain the maintenance of pollen colour variation.

Computational Modeling of Anthocyanin Pathway Evolution: Biases, Hotspots, and Trade-offs

The alteration of metabolic pathways is a common mechanism underlying the evolution of new phenotypes. Flower color is a striking example of the importance of metabolic evolution in a complex phenotype, wherein shifts in the activity of the underlying pathway lead to a wide range of pigments. Although experimental work has identified common classes of mutations responsible for transitions among colors, we lack a unifying model that relates pathway function and activity to the evolution of distinct pigment phenotypes. One challenge in creating such a model is the branching structure of pigment pathways, which may lead to evolutionary trade-offs due to competition for shared substrates. In order to predict the effects of shifts in enzyme function and activity on pigment production, we created a simple kinetic model of a major plant pigmentation pathway: the anthocyanin pathway. This model describes the production of the three classes of blue, purple, and red anthocyanin pigments, and accordingly, includes multiple branches and substrate competition. We first studied the general behavior of this model using a naïve set of parameters. We then stochastically evolved the pathway toward a defined optimum and analyzed the patterns of fixed mutations. This approach allowed us to quantify the probability density of trajectories through pathway state space and identify the types and number of changes. Finally, we examined whether our simulated results qualitatively align with experimental observations, i.e., the predominance of mutations which change color by altering the function of branching genes in the pathway. These analyses provide a theoretical framework that can be used to predict the consequences of new mutations in terms of both pigment phenotypes and pleiotropic effects.

Retracing the molecular basis and evolutionary history of the loss of benzaldehyde emission in the genus Capsella

The transition from pollinator-mediated outbreeding to selfing has occurred many times in angiosperms. This is generally accompanied by a reduction in traits attracting pollinators, including reduced emission of floral scent. In Capsella, emission of benzaldehyde as a main component of floral scent has been lost in selfing C. rubella by mutation of cinnamate-CoA ligase CNL1. However, the biochemical basis and evolutionary history of this loss remain unknown, as does the reason for the absence of benzaldehyde emission in the independently derived selfer Capsella orientalis. We used plant transformation, in vitro enzyme assays, population genetics and quantitative genetics to address these questions. CNL1 has been inactivated twice independently by point mutations in C. rubella, causing a loss of enzymatic activity. Both inactive haplotypes are found within and outside of Greece, the centre of origin of C. rubella, indicating that they arose before its geographical spread. By contrast, the loss of benzaldehyde emission in C. orientalis is not due to an inactivating mutation in CNL1. CNL1 represents a hotspot for mutations that eliminate benzaldehyde emission, potentially reflecting the limited pleiotropy and large effect of its inactivation. Nevertheless, even closely related species have followed different evolutionary routes in reducing floral scent.

Evolution: Flip-Flopping Flower Color Defies Dollo's Law

Even complex traits can re-evolve after being lost. A new study details the molecular mechanisms causing the regain of floral color pigment in a lineage that evolved white flowers.