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

Molecular Phylogenomics Reveals the Deep Evolutionary History of Carnivory across Land Plants

Plastid molecular phylogenies that broadly sampled angiosperm lineages imply that carnivorous plants evolved at least 11 times independently in 13 families and 6 orders. Within and between these clades, the different prey capture strategies involving flypaper and pitfall structures arose in parallel with the subsequent evolution of snap traps and suction bladders. Attempts to discern the deep ontological history of carnivorous structures using multigene phylogenies have provided a plastid-level picture of sister relationships at the family level. Here, we present a molecular phylogeny of the angiosperms based on nuclear target sequence capture data (Angiosperms-353 probe set), assembled by the Kew Plant Trees of Life initiative, which aims to complete the tree of life for plants. This phylogeny encompasses all carnivorous and protocarnivorous families, although certain genera such as Philcoxia (Plantaginaceae) are excluded. This study offers a novel nuclear gene-based overview of relationships within and between carnivorous families and genera. Consistent with previous broadly sampled studies, we found that most carnivorous families are not affiliated with any single family. Instead, they emerge as sister groups to large clades comprising multiple non-carnivorous families. Additionally, we explore recent genomic studies across various carnivorous clades that examine the evolution of the carnivorous syndrome in relation to whole-genome duplication, subgenome dominance, small-scale gene duplication, and convergent evolution. Furthermore, we discuss insights into genome size evolution through the lens of carnivorous plant genomes.

Spatially resolved detection of small molecules from press-dried plant tissue using MALDI imaging

Premise

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) is a chemical imaging method that can visualize spatial distributions of particular molecules. Plant tissue imaging has so far mostly used cryosectioning, which can be impractical for the preparation of large-area imaging samples, such as full flower petals. Imaging unsectioned plant tissue presents its own difficulties in extracting metabolites to the surface due to the waxy cuticle.

Methods

We address this by using established delipidation techniques combined with a solvent vapor extraction prior to applying the matrix with many low-concentration sprays.

Results

Using this procedure, we imaged tissue from three different plant species (two flowers and one carnivorous plant leaf). Material factorization analysis of the resulting data reveals a wide range of plant-specific small molecules with varying degrees of localization to specific portions of the tissue samples, while facilitating detection and removal of signal from background sources.

Conclusions

This work demonstrates applicability of MALDI-MSI to press-dried plant samples without freezing or cryosectioning, setting the stage for spatially resolved molecule identification. Increased mass resolution and inclusion of tandem mass spectrometry are necessary next steps to allow more specific and reliable compound identification.

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.

Snatching Sundews-Analysis of Tentacle Movement in Two Species of Drosera in Terms of Response Rate, Response Time, and Speed of Movement

Drosera, Droseraceae, catch prey with sticky tentacles. Both Australian Drosera allantostigma and widespread D. rotundifolia show three types of anatomically different tentacles: short, peripheral, and snap-tentacles. The latter two are capable of fast movement. This motion was analysed after mechanical, chemical, and electrical stimulation with respect to response rate, response time, and angular velocity of bending. Compared to D. rotundifolia, D. allantostigma responds more frequently and faster; the tentacles bend with higher angular velocity. Snap-tentacles have a lower response rate, shorter response time, and faster angular velocity. The response rates for chemical and electrical stimuli are similar, and higher than the rates for mechanical stimulus. The response time is not dependent on stimulus type. The higher motility in D. allantostigma indicates increased dependence on mechanical prey capture, and a reduced role of adhesive mucilage. The same tentacle types are present in both species and show similar motility patterns. The lower response rate of snap-tentacles might be a safety measure against accidental triggering, since the motion of snap-tentacles is irreversible and tissue destructive. Furthermore, tentacles seem to discern stimuli and respond specifically. The established model of stereotypical tentacle movement may not fully explain these observations.

Molecular and in vivo studies of a glutamate-class prolyl-endopeptidase for coeliac disease therapy

The digestion of gluten generates toxic peptides, among which a highly immunogenic proline-rich 33-mer from wheat α-gliadin, that trigger coeliac disease. Neprosin from the pitcher plant is a reported prolyl endopeptidase. Here, we produce recombinant neprosin and its mutants, and find that full-length neprosin is a zymogen, which is self-activated at gastric pH by the release of an all-β pro-domain via a pH-switch mechanism featuring a lysine plug. The catalytic domain is an atypical 7+8-stranded β-sandwich with an extended active-site cleft containing an unprecedented pair of catalytic glutamates. Neprosin efficiently degrades both gliadin and the 33-mer in vitro under gastric conditions and is reversibly inactivated at pH > 5. Moreover, co-administration of gliadin and the neprosin zymogen at the ratio 500:1 reduces the abundance of the 33-mer in the small intestine of mice by up to 90%. Neprosin therefore founds a family of eukaryotic glutamate endopeptidases that fulfils requisites for a therapeutic glutenase.

Celiac disease is characterized by intolerance to gluten, a cereal protein. Here, the authors show that neprosin, a glutamate peptidase from the pitcher plant, efficiently cleaves gluten components under physiological conditions in vitro and in the gut of mice.

The genome of the recretohalophyte Limonium bicolor provides insights into salt gland development and salinity adaptation during terrestrial evolution

Halophytes have evolved specialized strategies to cope with high salinity. The extreme halophyte sea lavender (Limonium bicolor) lacks trichomes but possesses salt glands on its epidermis that can excrete harmful ions, such as sodium, to avoid salt damage. Here, we report a high-quality, 2.92-Gb, chromosome-scale L. bicolor genome assembly based on a combination of Illumina short reads, single-molecule, real-time long reads, chromosome conformation capture (Hi-C) data, and Bionano genome maps, greatly enriching the genomic information on recretohalophytes with multicellular salt glands. Although the L. bicolor genome contains genes that show similarity to trichome fate genes from Arabidopsis thaliana, it lacks homologs of the decision fate genes GLABRA3, ENHANCER OF GLABRA3, GLABRA2, TRANSPARENT TESTA GLABRA2, and SIAMESE, providing a molecular explanation for the absence of trichomes in this species. We identified key genes (LbHLH and LbTTG1) controlling salt gland development among classical trichome homologous genes and confirmed their roles by showing that their mutations markedly disrupted salt gland initiation, salt secretion, and salt tolerance, thus offering genetic support for the long-standing hypothesis that salt glands and trichomes may share a common origin. In addition, a whole-genome duplication event occurred in the L. bicolor genome after its divergence from Tartary buckwheat and may have contributed to its adaptation to high salinity. The L. bicolor genome resource and genetic evidence reported in this study provide profound insights into plant salt tolerance mechanisms that may facilitate the engineering of salt-tolerant crops.

Genome-Wide Analyses of Aspartic Proteases on Potato Genome (Solanum tuberosum): Generating New Tools to Improve the Resistance of Plants to Abiotic Stress

Aspartic proteases are proteolytic enzymes widely distributed in living organisms and viruses. Although they have been extensively studied in many plant species, they are poorly described in potatoes. The present study aimed to identify and characterize S. tuberosum aspartic proteases. Gene structure, chromosome and protein domain organization, phylogeny, and subcellular predicted localization were analyzed and integrated with RNAseq data from different tissues, organs, and conditions focused on abiotic stress. Sixty-two aspartic protease genes were retrieved from the potato genome, distributed in 12 chromosomes. A high number of intronless genes and segmental and tandem duplications were detected. Phylogenetic analysis revealed eight StAP groups, named from StAPI to StAPVIII, that were differentiated into typical (StAPI), nucellin-like (StAPIIIa), and atypical aspartic proteases (StAPII, StAPIIIb to StAPVIII). RNAseq data analyses showed that gene expression was consistent with the presence of cis-acting regulatory elements on StAP promoter regions related to water deficit. The study presents the first identification and characterization of 62 aspartic protease genes and proteins on the potato genome and provides the baseline material for functional gene determinations and potato breeding programs, including gene editing mediated by CRISPR.

Signaling and transport processes related to the carnivorous lifestyle of plants living on nutrient-poor soil

In Eukaryotes, long-distance and rapid signal transmission is required in order to be able to react fast and flexibly to external stimuli. This long-distance signal transmission cannot take place by diffusion of signal molecules from the site of perception to the target tissue, as their speed is insufficient. Therefore, for adequate stimulus transmission, plants as well as animals make use of electrical signal transmission, as this can quickly cover long distances. This update summarises the most important advances in plant electrical signal transduction with a focus on the carnivorous Venus flytrap. It highlights the different types of electrical signals, examines their underlying ion fluxes and summarises the carnivorous processes downstream of the electrical signals.

Recent ecophysiological, biochemical and evolutional insights into plant carnivory

BACKGROUND

Carnivorous plants are an ecological group of approx. 810 vascular species which capture and digest animal prey, absorb prey-derived nutrients and utilize them to enhance their growth and development. Extant carnivorous plants have evolved in at least ten independent lineages, and their adaptive traits represent an example of structural and functional convergence. Plant carnivory is a result of complex adaptations to mostly nutrient-poor, wet and sunny habitats when the benefits of carnivory exceed the costs. With a boost in interest and extensive research in recent years, many aspects of these adaptations have been clarified (at least partly), but many remain unknown.

SCOPE

We provide some of the most recent insights into substantial ecophysiological, biochemical and evolutional particulars of plant carnivory from the functional viewpoint. We focus on those processes and traits in carnivorous plants associated with their ecological characterization, mineral nutrition, cost-benefit relationships, functioning of digestive enzymes and regulation of the hunting cycle in traps. We elucidate mechanisms by which uptake of prey-derived nutrients leads to stimulation of photosynthesis and root nutrient uptake.

CONCLUSIONS

Utilization of prey-derived mineral (mainly N and P) and organic nutrients is highly beneficial for plants and increases the photosynthetic rate in leaves as a prerequisite for faster plant growth. Whole-genome and tandem gene duplications brought gene material for diversification into carnivorous functions and enabled recruitment of defence-related genes. Possible mechanisms for the evolution of digestive enzymes are summarized, and a comprehensive picture on the biochemistry and regulation of prey decomposition and prey-derived nutrient uptake is provided.

Research advances in and prospects of ornamental plant genomics

The term 'ornamental plant' refers to all plants with ornamental value, which generally have beautiful flowers or special plant architectures. China is rich in ornamental plant resources and known as the "mother of gardens". Genomics is the science of studying genomes and is useful for carrying out research on genome evolution, genomic variations, gene regulation, and important biological mechanisms based on detailed genome sequence information. Due to the diversity of ornamental plants and high sequencing costs, the progress of genome research on ornamental plants has been slow for a long time. With the emergence of new sequencing technologies and a reduction in costs since the whole-genome sequencing of the first ornamental plant (Prunus mume) was completed in 2012, whole-genome sequencing of more than 69 ornamental plants has been completed in <10 years. In this review, whole-genome sequencing and resequencing of ornamental plants will be discussed. We provide analysis with regard to basic data from whole-genome studies of important ornamental plants, the regulation of important ornamental traits, and application prospects.

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.

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.

The Droserasin 1 PSI: A Membrane-Interacting Antimicrobial Peptide from the Carnivorous Plant Drosera capensis

The Droserasins, aspartic proteases from the carnivorous plant Drosera capensis, contain a 100-residue plant-specific insert (PSI) that is post-translationally cleaved and independently acts as an antimicrobial peptide. PSIs are of interest not only for their inhibition of microbial growth, but also because they modify the size of lipid vesicles and strongly interact with biological membranes. PSIs may therefore be useful for modulating lipid systems in NMR studies of membrane proteins. Here we present the expression and biophysical characterization of the Droserasin 1 PSI (D1 PSI.) This peptide is monomeric in solution and maintains its primarily α-helical secondary structure over a wide range of temperatures and pH values, even under conditions where its three disulfide bonds are reduced. Vesicle fusion assays indicate that the D1 PSI strongly interacts with bacterial and fungal lipids at pH 5 and lower, consistent with the physiological pH of D. capensis mucilage. It binds lipids with a variety of head groups, highlighting its versatility as a potential stabilizer for lipid nanodiscs. Solid-state NMR spectra collected at a field strength of 36 T, using a unique series-connected hybrid magnet, indicate that the peptide is folded and strongly bound to the membrane. Molecular dynamics simulations indicate that the peptide is stable as either a monomer or a dimer in a lipid bilayer. Both the monomer and the dimer allow the passage of water through the membrane, albeit at different rates.

Roles of Plant-Specific Inserts in Plant Defense

Ubiquitously expressed in plants, the plant-specific insert (PSI) of typical plant aspartic proteases (tpAPs) has been associated with plant development, stress response, and defense processes against invading pathogens. Despite sharing high sequence identity, structural studies revealed possible different mechanisms of action among species. The PSI induces signaling pathways of defense hormones in vivo and demonstrates broad-spectrum activity against phytopathogens in vitro. Recent characterization of the PSI-tpAP relationship uncovered novel, nonconventional intracellular protein transport pathways and improved tpAP production yields for industrial applications. In spite of research to date, relatively little is known about the structure-function relationships of PSIs. A comprehensive understanding of their biological roles may benefit plant protection strategies against virulent phytopathogens.

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.

Atypical and nucellin-like aspartic proteases: emerging players in plant developmental processes and stress responses

Members of the pepsin-like family (A1) of aspartic proteases (APs) are widely distributed in plants. A large number of genes encoding putative A1 APs are found in different plant genomes, the vast majority of which exhibit distinct features when compared with the so-called typical APs (and, therefore, grouped as atypical and nucellin-like APs). These features include the absence of the plant-specific insert; an unusually high number of cysteine residues; the nature of the amino acids preceding the first catalytic aspartate; and unexpected localizations. The over-representation of atypical and nucellin-like APs in plants is suggestive of greater diversification of protein functions and a more regulatory role for these APs, as compared with the housekeeping function generally attributed to typical APs. New functions have been uncovered for non-typical APs, with proposed roles in biotic and abiotic stress responses, chloroplast metabolism, and reproductive development, clearly suggesting functional specialization and tight regulation of activity. Furthermore, unusual enzymatic properties have also been documented for some of these proteases. Here, we give an overview of the current knowledge on the distinctive features and functions of both atypical and nucellin-like APs, and discuss this emerging pattern of functional complexity and specialization among plant pepsin-like proteases.

Evolutionary and structural analyses uncover a role for solvent interactions in the diversification of cocoonases in butterflies

Multi-omic approaches promise to supply the power to detect genes underlying disease and fitness-related phenotypes. Optimal use of the resulting profusion of data requires detailed investigation of individual candidate genes, a challenging proposition. Here, we combine transcriptomic and genomic data with molecular modelling of candidate enzymes to characterize the evolutionary history and function of the serine protease cocoonase. Heliconius butterflies possess the unique ability to feed on pollen; recent work has identified cocoonase as a candidate gene in pollen digestion. Cocoonase was first described in moths, where it aids in eclosure from the cocoon and is present as a single copy gene. In heliconiine butterflies it is duplicated and highly expressed in the mouthparts of adults. At least six copies of cocoonase are present in Heliconius melpomene and copy number varies across H. melpomene sub-populations. Most cocoonase genes are under purifying selection, however branch-site analyses suggest cocoonase 3 genes may have evolved under episodic diversifying selection. Molecular modelling of cocoonase proteins and examination of their predicted structures revealed that the active site region of each type has a similar structure to trypsin, with the same predicted substrate specificity across types. Variation among heliconiine cocoonases instead lies in the outward-facing residues involved in solvent interaction. Thus, the neofunctionalization of cocoonase duplicates appears to have resulted from the need for these serine proteases to operate in diverse biochemical environments. We suggest that cocoonase may have played a buffering role in feeding during the diversification of Heliconius across the neotropics by enabling these butterflies to digest protein from a range of biochemical milieux.

Protein structure networks provide insight into active site flexibility in esterase/lipases from the carnivorous plant Drosera capensis

In plants, esterase/lipases perform transesterification reactions, playing an important role in the synthesis of useful molecules, such as those comprising the waxy coatings of leaf surfaces. Plant genomes and transcriptomes have provided a wealth of data about expression patterns and the circumstances under which these enzymes are upregulated, e.g. pathogen defense and response to drought; however, predicting their functional characteristics from genomic or transcriptome data is challenging due to weak sequence conservation among the diverse members of this group. Although functional sequence blocks mediating enzyme activity have been identified, progress to date has been hampered by the paucity of information on the structural relationships among these regions and how they affect substrate specificity. Here we present methodology for predicting overall protein flexibility and active site flexibility based on molecular modeling and analysis of protein structure networks (PSNs). We define two new types of specialized PSNs: sequence region networks (SRNs) and active site networks (ASNs), which provide parsimonious representations of molecular structure in reference to known features of interest. Our approach, intended as an aid to target selection for poorly characterized enzyme classes, is demonstrated for 26 previously uncharacterized esterase/lipases from the genome of the carnivorous plant Drosera capensis and validated using a case/control design. Analysis of the network relationships among functional blocks and among the chemical moieties making up the catalytic triad reveals potentially functionally significant differences that are not apparent from sequence analysis alone.

Discovery of digestive enzymes in carnivorous plants with focus on proteases

BACKGROUND

Carnivorous plants have been fascinating researchers with their unique characters and bioinspired applications. These include medicinal trait of some carnivorous plants with potentials for pharmaceutical industry.

METHODS

This review will cover recent progress based on current studies on digestive enzymes secreted by different genera of carnivorous plants: Drosera (sundews), Dionaea (Venus flytrap), Nepenthes (tropical pitcher plants), Sarracenia (North American pitcher plants), Cephalotus (Australian pitcher plants), Genlisea (corkscrew plants), and Utricularia (bladderworts).

RESULTS

Since the discovery of secreted protease nepenthesin in Nepenthes pitcher, digestive enzymes from carnivorous plants have been the focus of many studies. Recent genomics approaches have accelerated digestive enzyme discovery. Furthermore, the advancement in recombinant technology and protein purification helped in the identification and characterisation of enzymes in carnivorous plants.

DISCUSSION

These different aspects will be described and discussed in this review with focus on the role of secreted plant proteases and their potential industrial applications.

Evaluating the adaptive evolutionary convergence of carnivorous plant taxa through functional genomics

Carnivorous plants are striking examples of evolutionary convergence, displaying complex and often highly similar adaptations despite lack of shared ancestry. Using available carnivorous plant genomes along with non-carnivorous reference taxa, this study examines the convergence of functional overrepresentation of genes previously implicated in plant carnivory. Gene Ontology (GO) coding was used to quantitatively score functional representation in these taxa, in terms of proportion of carnivory-associated functions relative to all functional sequence. Statistical analysis revealed that, in carnivorous plants as a group, only two of the 24 functions tested showed a signal of substantial overrepresentation. However, when the four carnivorous taxa were analyzed individually, 11 functions were found to be significant in at least one taxon. Though carnivorous plants collectively may show overrepresentation in functions from the predicted set, the specific functions that are overrepresented vary substantially from taxon to taxon. While it is possible that some functions serve a similar practical purpose such that one taxon does not need to utilize both to achieve the same result, it appears that there are multiple approaches for the evolution of carnivorous function in plant genomes. Our approach could be applied to tests of functional convergence in other systems provided on the availability of genomes and annotation data for a group.

Expression, activation and processing of a novel plant milk-clotting aspartic protease in Pichia pastoris

Galium verum, also known as Lady's Bedstraw or Cheese Rennet, is an herbaceous perennial plant traditionally used in cheese-making. We used RACE PCR to isolate novel enzymes from Galium verum with the ability to clot milk. This approach generated two cDNA sequences (named preprogaline A and B) encoding proteins displaying the typical plant aspartic protease primary structure. Preprogaline B was expressed in the yeast Pichia pastoris, after deleting and replacing its original signal peptide with the yeast α-factor signal peptide from Saccharomyces cerevisiae. The secreted recombinant protein was obtained by growing P. pastoris in YPD medium and had the ability to clot milk. The mature form of progaline B is a heterodimeric glycosylated enzyme, with a molecular weight of approximately 48 kDa, that contains a heavy (30.7 kDa) and a light (13.5 kDa) polypeptide chains linked by disulfide bonds. Western blot analysis revealed that progaline B is activated by the acidification of the yeast culture medium and that enzymatic activation requires two steps. First the precursor protein is cleaved into two polypeptide chains by partial removal of the plant-specific insert (PSI) present in plant aspartic proteases; this is later followed by propeptide removal. By altering the pH of the P. pastoris culture medium, we were able to obtain either active or inactive forms of the enzyme. Recombinant progaline B displayed a κ-casein hydrolysis pattern analogous to those produced by the animal and microbial coagulants currently used in the dairy industry, but it exhibited a different digestion profile on α- and β-caseins. The plant protease progaline B displays milk-clotting activities suitable for the production of novel dairy products.

Current Strategies of Polyploid Plant Genome Sequence Assembly

Polyploidy or duplication of an entire genome occurs in the majority of angiosperms. The understanding of polyploid genomes is important for the improvement of those crops, which humans rely on for sustenance and basic nutrition. As climate change continues to pose a potential threat to agricultural production, there will increasingly be a demand for plant cultivars that can resist biotic and abiotic stresses and also provide needed and improved nutrition. In the past decade, Next Generation Sequencing (NGS) has fundamentally changed the genomics landscape by providing tools for the exploration of polyploid genomes. Here, we review the challenges of the assembly of polyploid plant genomes, and also present recent advances in genomic resources and functional tools in molecular genetics and breeding. As genomes of diploid and less heterozygous progenitor species are increasingly available, we discuss the lack of complexity of these currently available reference genomes as they relate to polyploid crops. Finally, we review recent approaches of haplotyping by phasing and the impact of third generation technologies on polyploid plant genome assembly.

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.

Structure prediction and network analysis of chitinases from the Cape sundew, Drosera capensis

BACKGROUND

Carnivorous plants possess diverse sets of enzymes with novel functionalities applicable to biotechnology, proteomics, and bioanalytical research. Chitinases constitute an important class of such enzymes, with future applications including human-safe antifungal agents and pesticides. Here, we compare chitinases from the genome of the carnivorous plant Drosera capensis to those from related carnivorous plants and model organisms.

METHODS

Using comparative modeling, in silico maturation, and molecular dynamics simulation, we produce models of the mature enzymes in aqueous solution. We utilize network analytic techniques to identify similarities and differences in chitinase topology.

RESULTS

Here, we report molecular models and functional predictions from protein structure networks for eleven new chitinases from D. capensis, including a novel class IV chitinase with two active domains. This architecture has previously been observed in microorganisms but not in plants. We use a combination of comparative and de novo structure prediction followed by molecular dynamics simulation to produce models of the mature forms of these proteins in aqueous solution. Protein structure network analysis of these and other plant chitinases reveal characteristic features of the two major chitinase families.

GENERAL SIGNIFICANCE

This work demonstrates how computational techniques can facilitate quickly moving from raw sequence data to refined structural models and comparative analysis, and to select promising candidates for subsequent biochemical characterization. This capability is increasingly important given the large and growing body of data from high-throughput genome sequencing, which makes experimental characterization of every target impractical.

Sequence comparison, molecular modeling, and network analysis predict structural diversity in cysteine proteases from the Cape sundew, Drosera capensis

Carnivorous plants represent a so far underexploited reservoir of novel proteases with potentially useful activities. Here we investigate 44 cysteine proteases from the Cape sundew, Drosera capensis, predicted from genomic DNA sequences. D. capensis has a large number of cysteine protease genes; analysis of their sequences reveals homologs of known plant proteases, some of which are predicted to have novel properties. Many functionally significant sequence and structural features are observed, including targeting signals and occluding loops. Several of the proteases contain a new type of granulin domain. Although active site residues are conserved, the sequence identity of these proteases to known proteins is moderate to low; therefore, comparative modeling with all-atom refinement and subsequent atomistic MD-simulation is used to predict their 3D structures. The structure prediction data, as well as analysis of protein structure networks, suggest multifarious variations on the papain-like cysteine protease structural theme. This in silico methodology provides a general framework for investigating a large pool of sequences that are potentially useful for biotechnology applications, enabling informed choices about which proteins to investigate in the laboratory.

•

44 new cysteine proteases from the carnivorous plant Drosera capensis are described.

•

Structure prediction and molecular dynamics simulation predict overall folds similar to papain.

•

Functionally significant sequence and structural features are observed, including targeting signals and occluding loops.

•

Several of the proteases contain a new type of granulin domain.

•

Protein structure networks reveal global differences in interactions among chemical groups.