Transcriptome and genome size analysis of the Venus flytrap

Omics Approaches in Uncovering Molecular Evolution and Physiology of Botanical Carnivory

Systems biology has been increasingly applied with multiple omics for a holistic comprehension of complex biological systems beyond the reductionist approach that focuses on individual molecules. Different high-throughput omics approaches, including genomics, transcriptomics, metagenomics, proteomics, and metabolomics have been implemented to study the molecular mechanisms of botanical carnivory. This covers almost all orders of carnivorous plants, namely Caryophyllales, Ericales, Lamiales, and Oxalidales, except Poales. Studies using single-omics or integrated multi-omics elucidate the compositional changes in nucleic acids, proteins, and metabolites. The omics studies on carnivorous plants have led to insights into the carnivory origin and evolution, such as prey capture and digestion as well as the physiological adaptations of trap organ formation. Our understandings of botanical carnivory are further enhanced by the discoveries of digestive enzymes and transporter proteins that aid in efficient nutrient sequestration alongside dynamic molecular responses to prey. Metagenomics studies revealed the mutualistic relationships between microbes and carnivorous plants. Lastly, in silico analysis accelerated the functional characterization of new molecules from carnivorous plants. These studies have provided invaluable molecular data for systems understanding of carnivorous plants. More studies are needed to cover the diverse species with convergent evolution of botanical carnivory.

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.

Calcium dynamics during trap closure visualized in transgenic Venus flytrap

The leaves of the carnivorous plant Venus flytrap, Dionaea muscipula (Dionaea) close rapidly to capture insect prey. The closure response usually requires two successive mechanical stimuli to sensory hairs on the leaf blade within approximately 30 s (refs. ^1-4). An unknown biological system in Dionaea is thought to memorize the first stimulus and transduce the signal from the sensory hair to the leaf blade². Here, we link signal memory to calcium dynamics using transgenic Dionaea expressing a Ca²⁺ sensor. Stimulation of a sensory hair caused an increase in cytosolic Ca²⁺ concentration ([Ca²⁺]_cyt) starting in the sensory hair and spreading to the leaf blade. A second stimulus increased [Ca²⁺]_cyt to an even higher level, meeting a threshold that is correlated to the leaf blade closure. Because [Ca²⁺]_cyt gradually decreased after the first stimulus, the [Ca²⁺]_cyt increase induced by the second stimulus was insufficient to meet the putative threshold for movement after about 30 s. The Ca²⁺ wave triggered by mechanical stimulation moved an order of magnitude faster than that induced by wounding in petioles of Arabidopsis thaliana⁵ and Dionaea. The capacity for rapid movement has evolved repeatedly in flowering plants. This study opens a path to investigate the role of Ca²⁺ in plant movement mechanisms and their evolution.

Annotated genome sequences of the carnivorous plant Roridula gorgonias and a non-carnivorous relative, Clethra arborea

OBJECTIVE

Plant carnivory is distributed across the tree of life and has evolved at least six times independently, but sequenced and annotated nuclear genomes of carnivorous plants are currently lacking. We have sequenced and structurally annotated the nuclear genome of the carnivorous Roridula gorgonias and that of a non-carnivorous relative, Madeira's lily-of-the-valley-tree, Clethra arborea, both within the Ericales. This data adds an important resource to study the evolutionary genetics of plant carnivory across angiosperm lineages and also for functional and systematic aspects of plants within the Ericales.

RESULTS

Our assemblies have total lengths of 284 Mbp (R. gorgonias) and 511 Mbp (C. arborea) and show high BUSCO scores of 84.2% and 89.5%, respectively. We used their predicted genes together with publicly available data from other Ericales' genomes and transcriptomes to assemble a phylogenomic data set for the inference of a species tree. However, groups of orthologs showed a marked absence of species represented by a transcriptome. We discuss possible reasons and caution against combining predicted genes from genome- and transriptome-based assemblies.

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.

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.

Molecular characterization and evolution of carnivorous sundew (Drosera rotundifolia L.) class V β-1,3-glucanase

MAIN CONCLUSION

A gene for β-1,3-glucanase was isolated from carnivorous sundew. It is active in leaves and roots, but not in digestive glands. Analyses in transgenic tobacco suggest its function in germination. Ancestral plant β-1,3-glucanases (EC 3.2.1.39) played a role in cell division and cell wall remodelling, but divergent evolution has extended their roles in plant defense against stresses to decomposition of prey in carnivorous plants. As available gene sequences from carnivorous plants are rare, we isolated a glucanase gene from roundleaf sundew (Drosera rotundifolia L.) by a genome walking approach. Computational predictions recognized typical gene features and protein motifs described for other plant β-1,3-glucanases. Phylogenetic reconstructions suggest strong support for evolutionary relatedness to class V β-1,3-glucanases, including homologs that are active in the traps of related carnivorous species. The gene is expressed in sundew vegetative tissues but not in flowers and digestive glands, and encodes for a functional enzyme when expressed in transgenic tobacco. Detailed analyses of the supposed promoter both in silico and in transgenic tobacco suggest that this glucanase plays a role in development. Specific spatiotemporal activity was observed during transgenic seed germination. Later during growth, the sundew promoter was active in marginal and sub-marginal areas of apical true leaf meristems of young tobacco plants. These results suggest that the isolated glucanase gene is regulated endogenously, possibly by auxin. This is the first report on a nuclear gene study from sundew.

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.

The Venus Flytrap Dionaea muscipula Counts Prey-Induced Action Potentials to Induce Sodium Uptake

Carnivorous plants, such as the Venus flytrap (Dionaea muscipula), depend on an animal diet when grown in nutrient-poor soils. When an insect visits the trap and tilts the mechanosensors on the inner surface, action potentials (APs) are fired. After a moving object elicits two APs, the trap snaps shut, encaging the victim. Panicking preys repeatedly touch the trigger hairs over the subsequent hours, leading to a hermetically closed trap, which via the gland-based endocrine system is flooded by a prey-decomposing acidic enzyme cocktail. Here, we asked the question as to how many times trigger hairs have to be stimulated (e.g., now many APs are required) for the flytrap to recognize an encaged object as potential food, thus making it worthwhile activating the glands. By applying a series of trigger-hair stimulations, we found that the touch hormone jasmonic acid (JA) signaling pathway is activated after the second stimulus, while more than three APs are required to trigger an expression of genes encoding prey-degrading hydrolases, and that this expression is proportional to the number of mechanical stimulations. A decomposing animal contains a sodium load, and we have found that these sodium ions enter the capture organ via glands. We identified a flytrap sodium channel DmHKT1 as responsible for this sodium acquisition, with the number of transcripts expressed being dependent on the number of mechano-electric stimulations. Hence, the number of APs a victim triggers while trying to break out of the trap identifies the moving prey as a struggling Na⁺-rich animal and nutrition for the plant.

Video Abstract

•

Carnivorous Dionaea muscipula captures and processes nutrient- and sodium-rich prey

•

Via mechano-sensor stimulation, an animal meal is recognized, captured, and processed

•

Mechano-electrical waves induce JA signaling pathways that trigger prey digestion

•

Number of stimulations controls the production of digesting enzymes and uptake modules