Comparative transcriptomics of Venus flytrap (Dionaea muscipula) across stages of prey capture and digestion.

The Venus flytrap, Dionaea muscipula, is perhaps the world's best-known botanical carnivore. The act of prey capture and digestion along with its rapidly closing, charismatic traps make this species a compelling model for studying the evolution and fundamental biology of carnivorous plants. There is a growing body of research on the genome, transcriptome, and digestome of Dionaea muscipula, but surprisingly limited information on changes in trap transcript abundance over time since feeding. Here we present the results of a comparative transcriptomics project exploring the transcriptomic changes across seven timepoints in a 72-hour time series of prey digestion and three timepoints directly comparing triggered traps with and without prey items. We document a dynamic response to prey capture including changes in abundance of transcripts with Gene Ontology (GO) annotations related to digestion and nutrient uptake. Comparisons of traps with and without prey documented 174 significantly differentially expressed genes at 1 hour after triggering and 151 genes with significantly different abundances at 24 hours. Approximately 50% of annotated protein-coding genes in Venus flytrap genome exhibit change (10041 of 21135) in transcript abundance following prey capture. Whereas peak abundance for most of these genes was observed within 3 hours, an expression cluster of 3009 genes exhibited continuously increasing abundance over the 72-hour sampling period, and transcript for these genes with GO annotation terms including both catabolism and nutrient transport may continue to accumulate beyond 72 hours.

Distinctive plastome evolution in carnivorous angiosperms.

BACKGROUND: Independent origins of carnivory in multiple angiosperm families are fabulous examples of convergent evolution using a diverse array of life forms and habitats. Previous studies have indicated that carnivorous plants have distinct evolutionary trajectories of plastid genome (plastome) compared to their non-carnivorous relatives, yet the extent and general characteristics remain elusive. RESULTS: We compared plastomes from 9 out of 13 carnivorous families and their non-carnivorous relatives to assess carnivory-associated evolutionary patterns. We identified inversions in all sampled Droseraceae species and four species of Utricularia, Pinguicula, Darlingtonia and Triphyophyllum. A few carnivores showed distinct shifts in inverted repeat boundaries and the overall repeat contents. Many ndh genes, along with some other genes, were independently lost in several carnivorous lineages. We detected significant substitution rate variations in most sampled carnivorous lineages. A significant overall substitution rate acceleration characterizes the two largest carnivorous lineages of Droseraceae and Lentibulariaceae. We also observe moderate substitution rates acceleration in many genes of Cephalotus follicularis, Roridula gorgonias, and Drosophyllum lusitanicum. However, only a few genes exhibit significant relaxed selection. CONCLUSION: Our results indicate that the carnivory of plants have different effects on plastome evolution across carnivorous lineages. The complex mechanism under carnivorous habitats may have resulted in distinctive plastome evolution with conserved plastome in the Brocchinia hechtioides to strongly reconfigured plastomes structures in Droseraceae. Organic carbon obtained from prey and the efficiency of utilizing prey-derived nutrients might constitute possible explanation.

Dramatic difference in rate of chromosome number evolution among sundew (Drosera L., Droseraceae) lineages.

Chromosome number change is a driver of speciation in eukaryotic organisms. Carnivorous sundews in the plant genus Drosera L. exhibit single chromosome number variation both among and within species, especially in the Australian Drosera subg. Ergaleium D.C., potentially linked to atypical centromeres that span much of the length of the chromosomes. We critically reviewed the literature on chromosome counts in Drosera, verified the taxonomy and quality of the original counts, and reconstructed dated phylogenies. We used the BiChrom model to test whether rates of single chromosome number increase and decrease, and chromosome number doubling differed between D. subg. Ergaleium and the other subgenera and between self-compatible and self-incompatible lineages. The best model for chromosome evolution among subgenera had equal rates of chromosome number doubling but higher rates of single chromosome number change in D. subg. Ergaleium than in the other subgenera. Contrary to expectation, self-incompatible lineages had a significantly higher rate of single chromosome loss than self-compatible lineages. We found no evidence for an association between differences in single chromosome number changes and diploidization after polyploidy or centromere type. This study presents an exemplar for critically examining published cytological data and rigorously testing factors that may impact the rates of chromosome number evolution.

Mutational analysis of mechanosensitive ion channels in the carnivorous Venus flytrap plant.

How the Venus flytrap (Dionaea muscipula) evolved the remarkable ability to sense, capture, and digest animal prey for nutrients has long puzzled the scientific community.1 Recent genome and transcriptome sequencing studies have provided clues to the genes thought to play a role in these tasks.2,3,4,5 However, proving a causal link between these and any aspect of the plant's hunting behavior has been challenging due to the genetic intractability of this non-model organism. Here, we use CRISPR-Cas9 methods to generate targeted modifications in the Venus flytrap genome. The plant detects prey using touch-sensitive trigger hairs located on its bilobed leaves.6 Upon bending, these hairs convert mechanical touch signals into changes in the membrane potential of sensory cells, leading to rapid closure of the leaf lobes to ensnare the animal.7 Here, we generate mutations in trigger-hair-expressed MscS-like (MSL)-family mechanosensitive ion channel genes FLYCATCHER1 (FLYC1) and FLYCATCHER2 (FLYC2)5 and find that double-mutant plants have a reduced leaf-closing response to mechanical ultrasound stimulation. While we cannot exclude off-target effects of the CRISPR-Cas9 system, our genetic analysis is consistent with these and other functionally redundant mechanosensitive ion channels acting together to generate the sensory system necessary for prey detection.

Response of physiological parameters in Dionaea muscipula J. Ellis teratomas transformed with rolB oncogene.

BACKGROUND: Plant transformation with rol oncogenes derived from wild strains of Rhizobium rhizogenes is a popular biotechnology tool. Transformation effects depend on the type of rol gene, expression level, and the number of gene copies incorporated into the plant's genomic DNA. Although rol oncogenes are known as inducers of plant secondary metabolism, little is known about the physiological response of plants subjected to transformation. RESULTS: In this study, the physiological consequences of rolB oncogene incorporation into the DNA of Dionaea muscipula J. Ellis was evaluated at the level of primary and secondary metabolism. Examination of the teratoma (transformed shoots) cultures of two different clones (K and L) showed two different strategies for dealing with the presence of the rolB gene. Clone K showed an increased ratio of free fatty acids to lipids, superoxide dismutase activity, synthesis of the oxidised form of glutathione, and total pool of glutathione and carotenoids, in comparison to non-transformed plants (control). Clone L was characterised by increased accumulation of malondialdehyde, proline, activity of superoxide dismutase and catalase, total pool of glutathione, ratio of reduced form of glutathione to oxidised form, and accumulation of selected phenolic acids. Moreover, clone L had an enhanced ratio of total triglycerides to lipids and accumulated saccharose, fructose, glucose, and tyrosine. CONCLUSIONS: This study showed that plant transformation with the rolB oncogene derived from R. rhizogenes induces a pleiotropic effect in plant tissue after transformation. Examination of D. muscipula plant in the context of transformation with wild strains of R. rhizogenes can be a new source of knowledge about primary and secondary metabolites in transgenic organisms.

Stretch-activated ion channels identified in the touch-sensitive structures of carnivorous Droseraceae plants.

In response to touch, some carnivorous plants such as the Venus flytrap have evolved spectacular movements to capture animals for nutrient acquisition. However, the molecules that confer this sensitivity remain unknown. We used comparative transcriptomics to show that expression of three genes encoding homologs of the MscS-Like (MSL) and OSCA/TMEM63 family of mechanosensitive ion channels are localized to touch-sensitive trigger hairs of Venus flytrap. We focus here on the candidate with the most enriched expression in trigger hairs, the MSL homolog FLYCATCHER1 (FLYC1). We show that FLYC1 transcripts are localized to mechanosensory cells within the trigger hair, transfecting FLYC1 induces chloride-permeable stretch-activated currents in naïve cells, and transcripts coding for FLYC1 homologs are expressed in touch-sensing cells of Cape sundew, a related carnivorous plant of the Droseraceae family. Our data suggest that the mechanism of prey recognition in carnivorous Droseraceae evolved by co-opting ancestral mechanosensitive ion channels to sense touch.

The Venus flytrap trigger hair-specific potassium channel KDM1 can reestablish the K+ gradient required for hapto-electric signaling.

The carnivorous plant Dionaea muscipula harbors multicellular trigger hairs designed to sense mechanical stimuli upon contact with animal prey. At the base of the trigger hair, mechanosensation is transduced into an all-or-nothing action potential (AP) that spreads all over the trap, ultimately leading to trap closure and prey capture. To reveal the molecular basis for the unique functional repertoire of this mechanoresponsive plant structure, we determined the transcriptome of D. muscipula's trigger hair. Among the genes that were found to be highly specific to the trigger hair, the Shaker-type channel KDM1 was electrophysiologically characterized as a hyperpolarization- and acid-activated K+-selective channel, thus allowing the reuptake of K+ ions into the trigger hair's sensory cells during the hyperpolarization phase of the AP. During trap development, the increased electrical excitability of the trigger hair is associated with the transcriptional induction of KDM1. Conversely, when KDM1 is blocked by Cs+ in adult traps, the initiation of APs in response to trigger hair deflection is reduced, and trap closure is suppressed. KDM1 thus plays a dominant role in K+ homeostasis in the context of AP and turgor formation underlying the mechanosensation of trigger hair cells and thus D. muscipula's hapto-electric signaling.

A Chemometry of Aldrovanda vesiculosa L. (Waterwheel, Droseraceae) Populations.

The genus Aldrovanda is a Palaeogene element containing a single extant species, Aldrovanda vesiculosa L. This aquatic carnivorous herb has a very wide range of distribution, natively covering four continents; however, it is a critically endangered aquatic plant species worldwide. Previous studies revealed that A. vesiculosa had an extremely low genetic variation. The main aim of the present paper is to explore, using chemometric tools, the diversity of 16 A. vesiculosa populations from various sites from four continents (Eurasia, Africa, Australia). Using chemometric data as markers for genetic diversity, we show the relationships of 16 A. vesiculosa populations from various sites, including four continents. Phytochemical markers allowed the identification of five well-supported (bootstrap > 90%) groups among the 16 populations sampled. The principal component analysis data support the idea that the strongly related African (Botswana) and Australian (Kimberley, NT, NW Australia) populations are the most distant ones, separated from the European and Asian ones. However, considering the five Australian populations sampled, three are nested within the Eurasian group. The chemometric data are correlated positively with the geographical distances between the samples, which suggests a tendency toward isolation for the most distant populations.

Genomes of the Venus Flytrap and Close Relatives Unveil the Roots of Plant Carnivory.

Most plants grow and develop by taking up nutrients from the soil while continuously under threat from foraging animals. Carnivorous plants have turned the tables by capturing and consuming nutrient-rich animal prey, enabling them to thrive in nutrient-poor soil. To better understand the evolution of botanical carnivory, we compared the draft genome of the Venus flytrap (Dionaea muscipula) with that of its aquatic sister, the waterwheel plant Aldrovanda vesiculosa, and the sundew Drosera spatulata. We identified an early whole-genome duplication in the family as source for carnivory-associated genes. Recruitment of genes to the trap from the root especially was a major mechanism in the evolution of carnivory, supported by family-specific duplications. Still, these genomes belong to the gene poorest land plants sequenced thus far, suggesting reduction of selective pressure on different processes, including non-carnivorous nutrient acquisition. Our results show how non-carnivorous plants evolved into the most skillful green hunters on the planet.

Plastome-Wide Rearrangements and Gene Losses in Carnivorous Droseraceae.

The plastid genomes of four related carnivorous plants (Drosera regia, Drosera erythrorhiza, Aldrovanda vesiculosa, and Dionaea muscipula) were sequenced to examine changes potentially induced by the transition to carnivory. The plastid genomes of the Droseraceae show multiple rearrangements, gene losses, and large expansions or contractions of the inverted repeat. All the ndh genes are lost or nonfunctional, as well as in some of the species, clpP1, ycf1, ycf2 and some tRNA genes. Uniquely, among land plants, the trnK gene has no intron. Carnivory in the Droseraceae coincides with changes in plastid gene content similar to those induced by parasitism and mycoheterotrophy, suggesting parallel changes in chloroplast function due to the similar switch from autotrophy to (mixo-) heterotrophy. A molecular phylogeny of the taxa based on all shared plastid genes indicates that the "snap-traps" of Aldrovanda and Dionaea have a common origin.

Digestive glands extraction and precise pigment analysis support the exclusion of the carnivorous plant Dionaea muscipula Ellis from the Caryophyllales order.

In the order Caryophyllales, plants synthesize betalains instead of anthocyanins, with only two exceptions, the Caryophyllaceae and Molluginaceae. Dionaea muscipula Ellis was included in the Caryophyllales order but recent research based on genetic studies proposed the consideration of the Droseraceae family into the Nepenthales order. In this work we face the dilemma of the phylogenetic classification of Dionaea from a phytochemical point of view. Dionaea's pigments were analyzed by using techniques of structural analysis. Extracts from the leaves, mature stem and flowers of different specimens of Dionaea were analyzed, to find possible differences in the types of pigments or in their proportion in different parts of the plant. These extracts were analyzed by spectrophotometry, HPLC co-elution and ESI-MS/MS. In addition, digestive glands were extracted from the snap trap with minor sample manipulation and by reducing the non-pigmented plant tissue. Considering only the digestive glands instead of whole snap traps, the analyses allowed to quantitate and elucidate the structure of the compounds responsible for the red coloration: delphinidin-3-O-glucoside (myrtillin), cyanidin-3-O-glucoside (kuromanin) and a third compound, the aglycone cyanidin, detected in the species for the first time. The unambiguous results of the present work support the exclusion of Dionaea from the Caryophyllales.

Venus Flytrap: How an Excitable, Carnivorous Plant Works.

The carnivorous plant Dionaea possesses very sensitive mechanoreceptors. Upon contact with prey an action potential is triggered which, via an electrical network - comparable to the nervous system of vertebrates - rapidly closes its bivalved trap. The 'hunting cycle' comprises a constitutively activated mechanism for the rapid capture of prey, followed by a well-orchestrated sequence of activation of genes responsible for tight trap closure, digestion of the prey, and uptake of nutrients. Decisions on the step-by-step activation are based on 'counting' the number of stimulations of sensory organs. These remarkable animal-like skills in the carnivore are achieved not by taking over genes from its prey but by modifying and rearranging the functions of genes that are ubiquitous in plants.

Evolution of genome size and genomic GC content in carnivorous holokinetics (Droseraceae).

BACKGROUND AND AIMS: Studies in the carnivorous family Lentibulariaceae in the last years resulted in the discovery of the smallest plant genomes and an unusual pattern of genomic GC content evolution. However, scarcity of genomic data in other carnivorous clades still prevents a generalization of the observed patterns. Here the aim was to fill this gap by mapping genome evolution in the second largest carnivorous family, Droseraceae, where this evolution may be affected by chromosomal holokinetism in Drosera METHODS: The genome size and genomic GC content of 71 Droseraceae species were measured by flow cytometry. A dated phylogeny was constructed, and the evolution of both genomic parameters and their relationship to species climatic niches were tested using phylogeny-based statistics. KEY RESULTS: The 2C genome size of Droseraceae varied between 488 and 10 927 Mbp, and the GC content ranged between 37·1 and 44·7 %. The genome sizes and genomic GC content of carnivorous and holocentric species did not differ from those of their non-carnivorous and monocentric relatives. The genomic GC content positively correlated with genome size and annual temperature fluctuations. The genome size and chromosome numbers were inversely correlated in the Australian clade of Drosera CONCLUSIONS: Our results indicate that neither carnivory (nutrient scarcity) nor the holokinetism have a prominent effect on size and DNA base composition of Droseraceae genomes. However, the holokinetic drive seems to affect karyotype evolution in one of the major clades of Drosera Our survey confirmed that the evolution of GC content is tightly connected with the evolution of genome size and also with environmental conditions.

Novel proteases from the genome of the carnivorous plant Drosera capensis: Structural prediction and comparative analysis.

In his 1875 monograph on insectivorous plants, Darwin described the feeding reactions of Drosera flypaper traps and predicted that their secretions contained a "ferment" similar to mammalian pepsin, an aspartic protease. Here we report a high-quality draft genome sequence for the cape sundew, Drosera capensis, the first genome of a carnivorous plant from order Caryophyllales, which also includes the Venus flytrap (Dionaea) and the tropical pitcher plants (Nepenthes). This species was selected in part for its hardiness and ease of cultivation, making it an excellent model organism for further investigations of plant carnivory. Analysis of predicted protein sequences yields genes encoding proteases homologous to those found in other plants, some of which display sequence and structural features that suggest novel functionalities. Because the sequence similarity to proteins of known structure is in most cases too low for traditional homology modeling, 3D structures of representative proteases are predicted using comparative modeling with all-atom refinement. Although the overall folds and active residues for these proteins are conserved, we find structural and sequence differences consistent with a diversity of substrate recognition patterns. Finally, we predict differences in substrate specificities using in silico experiments, providing targets for structure/function studies of novel enzymes with biological and technological significance. Proteins 2016; 84:1517-1533. © 2016 Wiley Periodicals, Inc.

Venus flytrap carnivorous lifestyle builds on herbivore defense strategies.

Although the concept of botanical carnivory has been known since Darwin's time, the molecular mechanisms that allow animal feeding remain unknown, primarily due to a complete lack of genomic information. Here, we show that the transcriptomic landscape of the Dionaea trap is dramatically shifted toward signal transduction and nutrient transport upon insect feeding, with touch hormone signaling and protein secretion prevailing. At the same time, a massive induction of general defense responses is accompanied by the repression of cell death-related genes/processes. We hypothesize that the carnivory syndrome of Dionaea evolved by exaptation of ancient defense pathways, replacing cell death with nutrient acquisition.

Venus Flytrap HKT1-Type Channel Provides for Prey Sodium Uptake into Carnivorous Plant Without Conflicting with Electrical Excitability.

The animal diet of the carnivorous Venus flytrap, Dionaea muscipula, contains a sodium load that enters the capture organ via an HKT1-type sodium channel, expressed in special epithelia cells on the inner trap lobe surface. DmHKT1 expression and sodium uptake activity is induced upon prey contact. Here, we analyzed the HKT1 properties required for prey sodium osmolyte management of carnivorous Dionaea. Analyses were based on homology modeling, generation of model-derived point mutants, and their functional testing in Xenopus oocytes. We showed that the wild-type HKT1 and its Na(+)- and K(+)-permeable mutants function as ion channels rather than K(+) transporters driven by proton or sodium gradients. These structural and biophysical features of a high-capacity, Na(+)-selective ion channel enable Dionaea glands to manage prey-derived sodium loads without confounding the action potential-based information management of the flytrap.

Calcium sensor kinase activates potassium uptake systems in gland cells of Venus flytraps.

The Darwin plant Dionaea muscipula is able to grow on mineral-poor soil, because it gains essential nutrients from captured animal prey. Given that no nutrients remain in the trap when it opens after the consumption of an animal meal, we here asked the question of how Dionaea sequesters prey-derived potassium. We show that prey capture triggers expression of a K(+) uptake system in the Venus flytrap. In search of K(+) transporters endowed with adequate properties for this role, we screened a Dionaea expressed sequence tag (EST) database and identified DmKT1 and DmHAK5 as candidates. On insect and touch hormone stimulation, the number of transcripts of these transporters increased in flytraps. After cRNA injection of K(+)-transporter genes into Xenopus oocytes, however, both putative K(+) transporters remained silent. Assuming that calcium sensor kinases are regulating Arabidopsis K(+) transporter 1 (AKT1), we coexpressed the putative K(+) transporters with a large set of kinases and identified the CBL9-CIPK23 pair as the major activating complex for both transporters in Dionaea K(+) uptake. DmKT1 was found to be a K(+)-selective channel of voltage-dependent high capacity and low affinity, whereas DmHAK5 was identified as the first, to our knowledge, proton-driven, high-affinity potassium transporter with weak selectivity. When the Venus flytrap is processing its prey, the gland cell membrane potential is maintained around -120 mV, and the apoplast is acidified to pH 3. These conditions in the green stomach formed by the closed flytrap allow DmKT1 and DmHAK5 to acquire prey-derived K(+), reducing its concentration from millimolar levels down to trace levels.

Transcriptome and genome size analysis of the Venus flytrap.

The insectivorous Venus flytrap (Dionaea muscipula) is renowned from Darwin's studies of plant carnivory and the origins of species. To provide tools to analyze the evolution and functional genomics of D. muscipula, we sequenced a normalized cDNA library synthesized from mRNA isolated from D. muscipula flowers and traps. Using the Oases transcriptome assembler 79,165,657 quality trimmed reads were assembled into 80,806 cDNA contigs, with an average length of 679 bp and an N50 length of 1,051 bp. A total of 17,047 unique proteins were identified, and assigned to Gene Ontology (GO) and classified into functional categories. A total of 15,547 full-length cDNA sequences were identified, from which open reading frames were detected in 10,941. Comparative GO analyses revealed that D. muscipula is highly represented in molecular functions related to catalytic, antioxidant, and electron carrier activities. Also, using a single copy sequence PCR-based method, we estimated that the genome size of D. muscipula is approx. 3 Gb. Our genome size estimate and transcriptome analyses will contribute to future research on this fascinating, monotypic species and its heterotrophic adaptations.

Structural and functional characteristics of S-like ribonucleases from carnivorous plants.

Although the S-like ribonucleases (RNases) share sequence homology with the S-RNases involved in the self-incompatibility mechanism in plants, they are not associated with this mechanism. They usually function in stress responses in non-carnivorous plants and in carnivory in carnivorous plants. In this study, we clarified the structures of the S-like RNases of Aldrovanda vesiculosa, Nepenthes bicalcarata and Sarracenia leucophylla, and compared them with those of other plants. At ten positions, amino acid residues are conserved or almost conserved only for carnivorous plants (six in total). In contrast, two positions are specific to non-carnivorous plants. A phylogenetic analysis revealed that the S-like RNases of the carnivorous plants form a group beyond the phylogenetic relationships of the plants. We also prepared and characterized recombinant S-like RNases of Dionaea muscipula, Cephalotus follicularis, A. vesiculosa, N. bicalcarata and S. leucophylla, and RNS1 of Arabidopsis thaliana. The recombinant carnivorous plant enzymes showed optimum activities at about pH 4.0. Generally, poly(C) was digested less efficiently than poly(A), poly(I) and poly(U). The kinetic parameters of the recombinant D. muscipula enzyme (DM-I) and A. thaliana enzyme RNS1 were similar. The k cat/K m of recombinant RNS1 was the highest among the enzymes, followed closely by that of recombinant DM-I. On the other hand, the k cat/K m of the recombinant S. leucophylla enzyme was the lowest, and was ~1/30 of that for recombinant RNS1. The magnitudes of the k cat/K m values or k cat values for carnivorous plant S-like RNases seem to correlate negatively with the dependency on symbionts for prey digestion.