Five-leaf generalizations of the D-statistic reveal the directionality of admixture.

Over the past 15 years, the D-statistic, a four-taxon test for organismal admixture (hybridization, or introgression) which incorporates single nucleotide polymorphism data with allelic patterns ABBA and BABA, has seen considerable use. This statistic seeks to discern significant deviation from either a given species tree assumption, or from the balanced incomplete lineage sorting that could otherwise defy this species tree. However, while the D-statistic can successfully discriminate admixture from incomplete lineage sorting, it is not a simple matter to determine the directionality of admixture using only four-leaf tree models. As such, methods have been developed that use 5 leaves to evaluate admixture. Among these, the DFOIL method, which tests allelic patterns on the "symmetric" tree S = (((1,2),(3,4)),5), succeeds in finding admixture direction for many five-taxon examples. However, DFOIL does not make full use of all symmetry, nor can DFOIL function properly when ancient samples are included because of the reliance on singleton patterns (such as BAAAA and ABAAA). Here, we take inspiration from DFOIL to develop a new and completely general family of five-leaf admixture tests, dubbed Δ-statistics, that can either incorporate or exclude the singleton allelic patterns depending on individual taxon and age sampling choices. We describe two new shapes that are also fully testable, namely the "asymmetric" tree A = ((((1,2),3),4),5) and the "quasisymmetric" tree Q = (((1,2),3),(4,5)), which can considerably supplement the "symmetric" S = (((1,2),(3,4)),5) model used by DFOIL. We demonstrate the consistency of Δ-statistics under various simulated scenarios, and provide empirical examples using data from black, brown and polar bears, the latter also including two ancient polar bear samples from previous studies. Recently DFOIL and one of these ancient samples was used to argue for a dominant polar bear → brown bear introgression direction. However, we find, using both this ancient polar bear and our own, that by far the strongest signal using both DFOIL and Δ-statistics on tree S is actually bidirectional gene flow of indistinguishable direction. Further experiments on trees A and Q instead highlight what were likely two phases of admixture: one with stronger brown bear → polar bear introgression in ancient times, and a more recent phase with predominant polar bear → brown bear directionality. Code and documentation available at https://github.com/KalleLeppala/Delta-statistics.

Insights into bear evolution from a Pleistocene polar bear genome.

The polar bear (Ursus maritimus) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears (Ursus arctos). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear's lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.

The digestive systems of carnivorous plants.

To survive in the nutrient-poor habitats, carnivorous plants capture small organisms comprising complex substances not suitable for immediate reuse. The traps of carnivorous plants, which are analogous to the digestive systems of animals, are equipped with mechanisms for the breakdown and absorption of nutrients. Such capabilities have been acquired convergently over the past tens of millions of years in multiple angiosperm lineages by modifying plant-specific organs including leaves. The epidermis of carnivorous trap leaves bears groups of specialized cells called glands, which acquire substances from their prey via digestion and absorption. The digestive glands of carnivorous plants secrete mucilage, pitcher fluids, acids, and proteins, including digestive enzymes. The same (or morphologically distinct) glands then absorb the released compounds via various membrane transport proteins or endocytosis. Thus, these glands function in a manner similar to animal cells that are physiologically important in the digestive system, such as the parietal cells of the stomach and intestinal epithelial cells. Yet, carnivorous plants are equipped with strategies that deal with or incorporate plant-specific features, such as cell walls, epidermal cuticles, and phytohormones. In this review, we provide a systematic perspective on the digestive and absorptive capacity of convergently evolved carnivorous plants, with an emphasis on the forms and functions of glands.

The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation.

The avocado, Persea americana, is a fruit crop of immense importance to Mexican agriculture with an increasing demand worldwide. Avocado lies in the anciently diverged magnoliid clade of angiosperms, which has a controversial phylogenetic position relative to eudicots and monocots. We sequenced the nuclear genomes of the Mexican avocado race, P. americana var. drymifolia, and the most commercially popular hybrid cultivar, Hass, and anchored the latter to chromosomes using a genetic map. Resequencing of Guatemalan and West Indian varieties revealed that ∼39% of the Hass genome represents Guatemalan source regions introgressed into a Mexican race background. Some introgressed blocks are extremely large, consistent with the recent origin of the cultivar. The avocado lineage experienced 2 lineage-specific polyploidy events during its evolutionary history. Although gene-tree/species-tree phylogenomic results are inconclusive, syntenic ortholog distances to other species place avocado as sister to the enormous monocot and eudicot lineages combined. Duplicate genes descending from polyploidy augmented the transcription factor diversity of avocado, while tandem duplicates enhanced the secondary metabolism of the species. Phenylpropanoid biosynthesis, known to be elicited by Colletotrichum (anthracnose) pathogen infection in avocado, is one enriched function among tandems. Furthermore, transcriptome data show that tandem duplicates are significantly up- and down-regulated in response to anthracnose infection, whereas polyploid duplicates are not, supporting the general view that collections of tandem duplicates contribute evolutionarily recent "tuning knobs" in the genome adaptive landscapes of given species.

TCP and MADS-Box Transcription Factor Networks Regulate Heteromorphic Flower Type Identity in Gerbera hybrida.

The large sunflower family, Asteraceae, is characterized by compressed, flower-like inflorescences that may bear phenotypically distinct flower types. The CYCLOIDEA (CYC)/TEOSINTE BRANCHED1-like transcription factors (TFs) belonging to the TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) protein family are known to regulate bilateral symmetry in single flowers. In Asteraceae, they function at the inflorescence level, and were recruited to define differential flower type identities. Here, we identified upstream regulators of GhCYC3, a gene that specifies ray flower identity at the flower head margin in the model plant Gerbera hybrida We discovered a previously unidentified expression domain and functional role for the paralogous CINCINNATA-like TCP proteins. They function upstream of GhCYC3 and affect the developmental delay of marginal ray primordia during their early ontogeny. At the level of single flowers, the Asteraceae CYC genes show a unique function in regulating the elongation of showy ventral ligules that play a major role in pollinator attraction. We discovered that during ligule development, the E class MADS-box TF GRCD5 activates GhCYC3 expression. We propose that the C class MADS-box TF GAGA1 contributes to stamen development upstream of GhCYC3 Our data demonstrate how interactions among and between the conserved floral regulators, TCP and MADS-box TFs, contribute to the evolution of the elaborate inflorescence architecture of Asteraceae.

Tandem gene duplications drive divergent evolution of caffeine and crocin biosynthetic pathways in plants.

BACKGROUND: Plants have evolved a panoply of specialized metabolites that increase their environmental fitness. Two examples are caffeine, a purine psychotropic alkaloid, and crocins, a group of glycosylated apocarotenoid pigments. Both classes of compounds are found in a handful of distantly related plant genera (Coffea, Camellia, Paullinia, and Ilex for caffeine; Crocus, Buddleja, and Gardenia for crocins) wherein they presumably evolved through convergent evolution. The closely related Coffea and Gardenia genera belong to the Rubiaceae family and synthesize, respectively, caffeine and crocins in their fruits. RESULTS: Here, we report a chromosomal-level genome assembly of Gardenia jasminoides, a crocin-producing species, obtained using Oxford Nanopore sequencing and Hi-C technology. Through genomic and functional assays, we completely deciphered for the first time in any plant the dedicated pathway of crocin biosynthesis. Through comparative analyses with Coffea canephora and other eudicot genomes, we show that Coffea caffeine synthases and the first dedicated gene in the Gardenia crocin pathway, GjCCD4a, evolved through recent tandem gene duplications in the two different genera, respectively. In contrast, genes encoding later steps of the Gardenia crocin pathway, ALDH and UGT, evolved through more ancient gene duplications and were presumably recruited into the crocin biosynthetic pathway only after the evolution of the GjCCD4a gene. CONCLUSIONS: This study shows duplication-based divergent evolution within the coffee family (Rubiaceae) of two characteristic secondary metabolic pathways, caffeine and crocin biosynthesis, from a common ancestor that possessed neither complete pathway. These findings provide significant insights on the role of tandem duplications in the evolution of plant specialized metabolism.

Long-read sequencing uncovers the adaptive topography of a carnivorous plant genome.

Utricularia gibba, the humped bladderwort, is a carnivorous plant that retains a tiny nuclear genome despite at least two rounds of whole genome duplication (WGD) since common ancestry with grapevine and other species. We used a third-generation genome assembly with several complete chromosomes to reconstruct the two most recent lineage-specific ancestral genomes that led to the modern U. gibba genome structure. Patterns of subgenome dominance in the most recent WGD, both architectural and transcriptional, are suggestive of allopolyploidization, which may have generated genomic novelty and led to instantaneous speciation. Syntenic duplicates retained in polyploid blocks are enriched for transcription factor functions, whereas gene copies derived from ongoing tandem duplication events are enriched in metabolic functions potentially important for a carnivorous plant. Among these are tandem arrays of cysteine protease genes with trap-specific expression that evolved within a protein family known to be useful in the digestion of animal prey. Further enriched functions among tandem duplicates (also with trap-enhanced expression) include peptide transport (intercellular movement of broken-down prey proteins), ATPase activities (bladder-trap acidification and transmembrane nutrient transport), hydrolase and chitinase activities (breakdown of prey polysaccharides), and cell-wall dynamic components possibly associated with active bladder movements. Whereas independently polyploid Arabidopsis syntenic gene duplicates are similarly enriched for transcriptional regulatory activities, Arabidopsis tandems are distinct from those of U. gibba, while still metabolic and likely reflecting unique adaptations of that species. Taken together, these findings highlight the special importance of tandem duplications in the adaptive landscapes of a carnivorous plant genome.

Architecture and evolution of a minute plant genome.

It has been argued that the evolution of plant genome size is principally unidirectional and increasing owing to the varied action of whole-genome duplications (WGDs) and mobile element proliferation. However, extreme genome size reductions have been reported in the angiosperm family tree. Here we report the sequence of the 82-megabase genome of the carnivorous bladderwort plant Utricularia gibba. Despite its tiny size, the U. gibba genome accommodates a typical number of genes for a plant, with the main difference from other plant genomes arising from a drastic reduction in non-genic DNA. Unexpectedly, we identified at least three rounds of WGD in U. gibba since common ancestry with tomato (Solanum) and grape (Vitis). The compressed architecture of the U. gibba genome indicates that a small fraction of intergenic DNA, with few or no active retrotransposons, is sufficient to regulate and integrate all the processes required for the development and reproduction of a complex organism.

Dissecting functions of SEPALLATA-like MADS box genes in patterning of the pseudanthial inflorescence of Gerbera hybrida.

The pseudanthial inflorescences of the sunflower family, Asteraceae, mimic a solitary flower but are composed of multiple flowers. Our studies in Gerbera hybrida indicate functional diversification for SEPALLATA (SEP)-like MADS box genes that often function redundantly in other core eudicots. We conducted phylogenetic and expression analysis for eight SEP-like GERBERA REGULATOR OF CAPITULUM DEVELOPMENT (GRCD) genes, including previously unstudied gene family members. Transgenic gerbera plants were used to infer gene functions. Adding to the previously identified stamen and carpel functions for GRCD1 and GRCD2, two partially redundant genes, GRCD4 and GRCD5, were found to be indispensable for petal development. Stepwise conversion of floral organs into leaves in the most severe RNA interference lines suggest redundant and additive GRCD activities in organ identity regulation. We show conserved and redundant functions for several GRCD genes in regulation of flower meristem maintenance, while functional diversification for three SEP1/2/4 clade genes in regulation of inflorescence meristem patterning was observed. GRCD genes show both specialized and pleiotropic functions contributing to organ differentiation and flower meristem fate, and uniquely, to patterning of the inflorescence meristem. Altogether, we provide an example of how plant reproductive evolution has used conserved genetic modules for regulating the elaborate inflorescence architecture in Asteraceae.

Evolutionary Co-Option of Floral Meristem Identity Genes for Patterning of the Flower-Like Asteraceae Inflorescence.

The evolutionary success of Asteraceae, the largest family of flowering plants, has been attributed to the unique inflorescence architecture of the family, which superficially resembles an individual flower. Here, we show that Asteraceae inflorescences (flower heads, or capitula) resemble solitary flowers not only morphologically but also at the molecular level. By conducting functional analyses for orthologs of the flower meristem identity genes LEAFY (LFY) and UNUSUAL FLORAL ORGANS (UFO) in Gerbera hybrida, we show that GhUFO is the master regulator of flower meristem identity, while GhLFY has evolved a novel, homeotic function during the evolution of head-like inflorescences. Resembling LFY expression in a single flower meristem, uniform expression of GhLFY in the inflorescence meristem defines the capitulum as a determinate structure that can assume floral fate upon ectopic GhUFO expression. We also show that GhLFY uniquely regulates the ontogeny of outer, expanded ray flowers but not inner, compact disc flowers, indicating that the distinction of different flower types in Asteraceae is connected with their independent evolutionary origins from separate branching systems.

Genotyping-by-sequencing provides the first well-resolved phylogeny for coffee (Coffea) and insights into the evolution of caffeine content in its species: GBS coffee phylogeny and the evolution of caffeine content.

A comprehensive and meaningful phylogenetic hypothesis for the commercially important coffee genus (Coffea) has long been a key objective for coffee researchers. For molecular studies, progress has been limited by low levels of sequence divergence, leading to insufficient topological resolution and statistical support in phylogenetic trees, particularly for the major lineages and for the numerous species occurring in Madagascar. We report here the first almost fully resolved, broadly sampled phylogenetic hypothesis for coffee, the result of combining genotyping-by-sequencing (GBS) technology with a newly developed, lab-based workflow to integrate short read next-generation sequencing for low numbers of additional samples. Biogeographic patterns indicate either Africa or Asia (or possibly the Arabian Peninsula) as the most likely ancestral locality for the origin of the coffee genus, with independent radiations across Africa, Asia, and the Western Indian Ocean Islands (including Madagascar and Mauritius). The evolution of caffeine, an important trait for commerce and society, was evaluated in light of our phylogeny. High and consistent caffeine content is found only in species from the equatorial, fully humid environments of West and Central Africa, possibly as an adaptive response to increased levels of pest predation. Moderate caffeine production, however, evolved at least one additional time recently (between 2 and 4Mya) in a Madagascan lineage, which suggests that either the biosynthetic pathway was already in place during the early evolutionary history of coffee, or that caffeine synthesis within the genus is subject to convergent evolution, as is also the case for caffeine synthesis in coffee versus tea and chocolate.

Evolution of RLSB, a nuclear-encoded S1 domain RNA binding protein associated with post-transcriptional regulation of plastid-encoded rbcL mRNA in vascular plants.

BACKGROUND: RLSB, an S-1 domain RNA binding protein of Arabidopsis, selectively binds rbcL mRNA and co-localizes with Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) within chloroplasts of C3 and C4 plants. Previous studies using both Arabidopsis (C3) and maize (C4) suggest RLSB homologs are post-transcriptional regulators of plastid-encoded rbcL mRNA. While RLSB accumulates in all Arabidopsis leaf chlorenchyma cells, in C4 leaves RLSB-like proteins accumulate only within Rubisco-containing bundle sheath chloroplasts of Kranz-type species, and only within central compartment chloroplasts in the single cell C4 plant Bienertia. Our recent evidence implicates this mRNA binding protein as a primary determinant of rbcL expression, cellular localization/compartmentalization, and photosynthetic function in all multicellular green plants. This study addresses the hypothesis that RLSB is a highly conserved Rubisco regulatory factor that occurs in the chloroplasts all higher plants. RESULTS: Phylogenetic analysis has identified RLSB orthologs and paralogs in all major plant groups, from ancient liverworts to recent angiosperms. RLSB homologs were also identified in algae of the division Charophyta, a lineage closely related to land plants. RLSB-like sequences were not identified in any other algae, suggesting that it may be specific to the evolutionary line leading to land plants. The RLSB family occurs in single copy across most angiosperms, although a few species with two copies were identified, seemingly randomly distributed throughout the various taxa, although perhaps correlating in some cases with known ancient whole genome duplications. Monocots of the order Poales (Poaceae and Cyperaceae) were found to contain two copies, designated here as RLSB-a and RLSB-b, with only RLSB-a implicated in the regulation of rbcL across the maize developmental gradient. Analysis of microsynteny in angiosperms revealed high levels of conservation across eudicot species and for both paralogs in grasses, highlighting the possible importance of maintaining this gene and its surrounding genomic regions. CONCLUSIONS: Findings presented here indicate that the RLSB family originated as a unique gene in land plant evolution, perhaps in the common ancestor of charophytes and higher plants. Purifying selection has maintained this as a highly conserved single- or two-copy gene across most extant species, with several conserved gene duplications. Together with previous findings, this study suggests that RLSB has been sustained as an important regulatory protein throughout the course of land plant evolution. While only RLSB-a has been directly implicated in rbcL regulation in maize, RLSB-b could have an overlapping function in the co-regulation of rbcL, or may have diverged as a regulator of one or more other plastid-encoded mRNAs. This analysis confirms that RLSB is an important and unique photosynthetic regulatory protein that has been continuously expressed in land plants as they emerged and diversified from their ancient common ancestor.

High Gene Family Turnover Rates and Gene Space Adaptation in the Compact Genome of the Carnivorous Plant Utricularia gibba.

Utricularia gibba is an aquatic carnivorous plant with highly specialized morphology, featuring fibrous floating networks of branches and leaf-like organs, no recognizable roots, and bladder traps that capture and digest prey. We recently described the compressed genome of U. gibba as sufficient to control the development and reproduction of a complex organism. We hypothesized intense deletion pressure as a mechanism whereby most noncoding DNA was deleted, despite evidence for three independent whole-genome duplications (WGDs). Here, we explore the impact of intense genome fractionation in the evolutionary dynamics of U. gibba's functional gene space. We analyze U. gibba gene family turnover by modeling gene gain/death rates under a maximum-likelihood statistical framework. In accord with our deletion pressure hypothesis, we show that the U. gibba gene death rate is significantly higher than those of four other eudicot species. Interestingly, the gene gain rate is also significantly higher, likely reflecting the occurrence of multiple WGDs and possibly also small-scale genome duplications. Gene ontology enrichment analyses of U. gibba-specific two-gene orthogroups, multigene orthogroups, and singletons highlight functions that may represent adaptations in an aquatic carnivorous plant. We further discuss two homeodomain transcription factor gene families (WOX and HDG/HDZIP-IV) showing conspicuous differential expansions and contractions in U. gibba. Our results 1) reconcile the compactness of the U. gibba genome with its accommodation of a typical number of genes for a plant genome, and 2) highlight the role of high gene family turnover in the evolutionary diversification of U. gibba's functional gene space and adaptations to its unique lifestyle and highly specialized body plan.

The Fragaria vesca homolog of suppressor of overexpression of constans1 represses flowering and promotes vegetative growth.

In the annual long-day plant Arabidopsis thaliana, suppressor of overexpression of constans1 (SOC1) integrates endogenous and environmental signals to promote flowering. We analyzed the function and regulation of the SOC1 homolog (Fragaria vesca [Fv] SOC1) in the perennial short-day plant woodland strawberry (Fragaria vesca). We found that Fv SOC1 overexpression represses flower initiation under inductive short days, whereas its silencing causes continuous flowering in both short days and noninductive long days, similar to mutants in the floral repressor Fv terminal flower1 (Fv TFL1). Molecular analysis of these transgenic lines revealed that Fv SOC1 activates Fv TFL1 in the shoot apex, leading to the repression of flowering in strawberry. In parallel, Fv SOC1 regulates the differentiation of axillary buds to runners or axillary leaf rosettes, probably through the activation of gibberellin biosynthetic genes. We also demonstrated that Fv SOC1 is regulated by photoperiod and Fv flowering locus T1, suggesting that it plays a central role in the photoperiodic control of both generative and vegetative growth in strawberry. In conclusion, we propose that Fv SOC1 is a signaling hub that regulates yearly cycles of vegetative and generative development through separate genetic pathways.

Functional diversification of duplicated CYC2 clade genes in regulation of inflorescence development in Gerbera hybrida (Asteraceae).

The complex inflorescences (capitula) of Asteraceae consist of different types of flowers. In Gerbera hybrida (gerbera), the peripheral ray flowers are bilaterally symmetrical and lack functional stamens while the central disc flowers are more radially symmetrical and hermaphroditic. Proteins of the CYC2 subclade of the CYC/TB1-like TCP domain transcription factors have been recruited several times independently for parallel evolution of bilaterally symmetrical flowers in various angiosperm plant lineages, and have also been shown to regulate flower-type identity in Asteraceae. The CYC2 subclade genes in gerbera show largely overlapping gene expression patterns. At the level of single flowers, their expression domain in petals shows a spatial shift from the dorsal pattern known so far in species with bilaterally symmetrical flowers, suggesting that this change in expression may have evolved after the origin of Asteraceae. Functional analysis indicates that GhCYC2, GhCYC3 and GhCYC4 mediate positional information at the proximal-distal axis of the inflorescence, leading to differentiation of ray flowers, but that they also regulate ray flower petal growth by affecting cell proliferation until the final size and shape of the petals is reached. Moreover, our data show functional diversification for the GhCYC5 gene. Ectopic activation of GhCYC5 increases flower density in the inflorescence, suggesting that GhCYC5 may promote the flower initiation rate during expansion of the capitulum. Our data thus indicate that modification of the ancestral network of TCP factors has, through gene duplications, led to the establishment of new expression domains and to functional diversification.

Syntenic block overlap multiplicities with a panel of reference genomes provide a signature of ancient polyploidization events.

BACKGROUND: Following whole genome duplication (WGD), there is a compact distribution of gene similarities within the genome reflecting duplicate pairs of all the genes in the genome. With time, the distribution broadens and loses volume due to variable decay of duplicate gene similarity and to the process of duplicate gene loss. If there are two WGD, the older one becomes so reduced and broad that it merges with the tail of the distributions resulting from more recent events, and it becomes difficult to distinguish them. The goal of this paper is to advance statistical methods of identifying, or at least counting, the WGD events in the lineage of a given genome. METHODS: For a set of 15 angiosperm genomes, we analyze all 15 × 14 = 210 ordered pairs of target genome versus reference genome, using SynMap to find syntenic blocks. We consider all sets of B ≥ 2 syntenic blocks in the target genome that overlap in the reference genome as evidence of WGD activity in the target, whether it be one event or several. We hypothesize that in fitting an exponential function to the tail of the empirical distribution f (B) of block multiplicities, the size of the exponent will reflect the amount of WGD in the history of the target genome. RESULTS: By amalgamating the results from all reference genomes, a range of values of SynMap parameters, and alternative cutoff points for the tail, we find a clear pattern whereby multiple-WGD core eudicots have the smallest (negative) exponents, followed by core eudicots with only the single "γ" triplication in their history, followed by a non-core eudicot with a single WGD, followed by the monocots, with a basal angiosperm, the WGD-free Amborella having the largest exponent. CONCLUSION: The hypothesis that the exponent of the fit to the tail of the multiplicity distribution is a signature of the amount of WGD is verified, but there is also a clear complicating factor in the monocot clade, where a history of multiple WGD is not reflected in a small exponent.

Deep sequencing of the Mexican avocado transcriptome, an ancient angiosperm with a high content of fatty acids.

BACKGROUND: Avocado (Persea americana) is an economically important tropical fruit considered to be a good source of fatty acids. Despite its importance, the molecular and cellular characterization of biochemical and developmental processes in avocado is limited due to the lack of transcriptome and genomic information. RESULTS: The transcriptomes of seeds, roots, stems, leaves, aerial buds and flowers were determined using different sequencing platforms. Additionally, the transcriptomes of three different stages of fruit ripening (pre-climacteric, climacteric and post-climacteric) were also analyzed. The analysis of the RNAseqatlas presented here reveals strong differences in gene expression patterns between different organs, especially between root and flower, but also reveals similarities among the gene expression patterns in other organs, such as stem, leaves and aerial buds (vegetative organs) or seed and fruit (storage organs). Important regulators, functional categories, and differentially expressed genes involved in avocado fruit ripening were identified. Additionally, to demonstrate the utility of the avocado gene expression atlas, we investigated the expression patterns of genes implicated in fatty acid metabolism and fruit ripening. CONCLUSIONS: A description of transcriptomic changes occurring during fruit ripening was obtained in Mexican avocado, contributing to a dynamic view of the expression patterns of genes involved in fatty acid biosynthesis and the fruit ripening process.

De novo sequencing and analysis of Lophophora williamsii transcriptome, and searching for putative genes involved in mescaline biosynthesis.

BACKGROUND: Lophophora williamsii (commonly named peyote) is a small, spineless cactus with psychoactive alkaloids, particularly mescaline. Peyote utilizes crassulacean acid metabolism (CAM), an alternative form of photosynthesis that exists in succulents such as cacti and other desert plants. Therefore, its transcriptome can be considered an important resource for future research focused on understanding how these plants make more efficient use of water in marginal environments and also for research focused on better understanding of the overall mechanisms leading to production of plant natural products and secondary metabolites. RESULTS: In this study, two cDNA libraries were generated from L. williamsii. These libraries, representing buttons (tops of stems) and roots were sequenced using different sequencing platforms (GS-FLX, GS-Junior and PGM, respectively). A total of 5,541,550 raw reads were generated, which were assembled into 63,704 unigenes with an average length of 564.04 bp. A total of 25,149 unigenes (62.19 %) was annotated using public databases. 681 unigenes were found to be differentially expressed when comparing the two libraries, where 400 were preferentially expressed in buttons and 281 in roots. Some of the major alkaloids, including mescaline, were identified by GC-MS and relevant metabolic pathways were reconstructed using the Kyoto encyclopedia of genes and genomes database (KEGG). Subsequently, the expression patterns of preferentially expressed genes putatively involved in mescaline production were examined and validated by qRT-PCR. CONCLUSIONS: High throughput transcriptome sequencing (RNA-seq) analysis allowed us to efficiently identify candidate genes involved in mescaline biosynthetic pathway in L. williamsii; these included tyrosine/DOPA decarboxylase, hydroxylases, and O-methyltransferases. This study sets the theoretical foundation for bioassay design directed at confirming the participation of these genes in mescaline production.

Molecular evolutionary mechanisms driving functional diversification of the HSP90A family of heat shock proteins in eukaryotes.

The ubiquitous and conserved cytosolic heat-shock proteins 90 (HSP90A) perform essential functions in the cell. To understand the evolutionary origin of HSP90A functional diversification, we analyzed the distribution of HSP90A family from 54 species representing the main eukaryotic lineages. Three independent HSP90A duplications led to the paralog subfamilies HSP90AA (heat-stress inducible) and HSP90AB (constitutive) and trace back to key time points during vertebrate, seed plant, and yeast evolution. HSP90AA and HSP90AB present divergent selection pressures, positive selection (PS), and signatures of functional divergence (FD) after duplication. The differential evolutionary patterns support different mechanisms for HSP90A functional diversification in vertebrates and seed plants. Mapping of PS and FD residues onto the HSP90A structure suggests the acquisition of novel and/or specialized client protein and/or cochaperone binding functions. We propose these residues as targets for further experimental studies of HSP90A proteins, reported to be capacitors of rapid evolutionary change, and targets for anticancer therapeutics.

Functional diversification of duplicated chalcone synthase genes in anthocyanin biosynthesis of Gerbera hybrida.

• Chalcone synthase (CHS) is the key enzyme in the first committed step of the flavonoid biosynthetic pathway and catalyzes the stepwise condensation of 4-coumaroyl-CoA and malonyl-CoA to naringenin chalcone. In plants, CHS is often encoded by a small family of genes that are temporally and spatially regulated. Our earlier studies have shown that GCHS4 is highly activated by ectopic expression of an MYB-type regulator GMYB10 in gerbera (Gerbera hybrida). • The tissue- and development-specific expression patterns of three gerbera CHS genes were examined. Virus-induced gene silencing (VIGS) was used to knock down GCHS1 and GCHS4 separately in gerbera inflorescences. • Our data show that GCHS4 is the only CHS encoding gene that is expressed in the cyanidin-pigmented vegetative tissues of gerbera cv Terraregina. GCHS3 expression is pronounced in the pappus bristles of the flowers. Expression of both GCHS1 and GCHS4 is high in the epidermal cells of gerbera petals, but only GCHS1 is contributing to flavonoid biosynthesis. • Gerbera contains a family of three CHS encoding genes showing different spatial and temporal regulation. GCHS4 expression in gerbera petals is regulated post-transcriptionally, at the level of either translation elongation or protein stability.