Impact of recombination on polymorphism of genes encoding Kunitz-type protease inhibitors in the genus Solanum

Canonical or noncanonical? Structural plasticity of serine protease-binding loops in Kunitz-STI protease inhibitors

The Kunitz-Soybean Trypsin Inhibitor (Kunitz-STI) family is a large family of proteins with most of its members being protease inhibitors. The versatility of the inhibitory profile and the structural plasticity of these proteins, make this family a promising scaffold for designing new multifunctional proteins. Historically, Kunitz-STI inhibitors have been classified as canonical serine protease inhibitors, but new inhibitors with novel inhibition mechanisms have been described in recent years. Different inhibition mechanisms could be the result of different evolutionary pathways. In the present work, we performed a structural analysis of all the crystallographic structures available for Kunitz-STI inhibitors to characterize serine protease-binding loop structural features and locations. Our study suggests a relationship between the conformation of serine protease-binding loops and the inhibition mechanism, their location in the β-trefoil fold, and the plant source of the inhibitors. The classical canonical inhibitors of this family are restricted to plants from the Fabales order and bind their targets via the β4-β5 loop, whereas serine protease-binding loops in inhibitors from other plants lie mainly in the β5-β6 and β9-β10 loops. In addition, we found that the β5-β6 loop is used to inhibit two different families of serine proteases through a steric blockade inhibition mechanism. This work will help to change the general perception that all Kunitz-STI inhibitors are canonical inhibitors and proteins with protease-binding loops adopting noncanonical conformations are exceptions. Additionally, our results will help in the identification of protease-binding loops in uncharacterized or newly discovered inhibitors, and in the design of multifunctional proteins.

ProteinEvolverABC: coestimation of recombination and substitution rates in protein sequences by approximate Bayesian computation

MOTIVATION

The evolutionary processes of mutation and recombination, upon which selection operates, are fundamental to understand the observed molecular diversity. Unlike nucleotide sequences, the estimation of the recombination rate in protein sequences has been little explored, neither implemented in evolutionary frameworks, despite protein sequencing methods are largely used.

RESULTS

In order to accommodate this need, here I present a computational framework, called ProteinEvolverABC, to jointly estimate recombination and substitution rates from alignments of protein sequences. The framework implements the approximate Bayesian computation approach, with and without regression adjustments and includes a variety of substitution models of protein evolution, demographics and longitudinal sampling. It also implements several nuisance parameters such as heterogeneous amino acid frequencies and rate of change among sites and, proportion of invariable sites. The framework produces accurate coestimation of recombination and substitution rates under diverse evolutionary scenarios. As illustrative examples of usage, I applied it to several viral protein families, including coronaviruses, showing heterogeneous substitution and recombination rates.

AVAILABILITY AND IMPLEMENTATION

ProteinEvolverABC is freely available from https://github.com/miguelarenas/proteinevolverabc, includes a graphical user interface for helping the specification of the input settings, extensive documentation and ready-to-use examples. Conveniently, the simulations can run in parallel on multicore machines.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Phylogenetic Analysis of Different Ploidy Saccharum spontaneum Based on rDNA-ITS Sequences

Saccharum spontaneum L. is a crucial wild parent of modern sugarcane cultivars whose ploidy clones have been utilized successfully in improving the stress resistance and yield related traits of sugarcane cultivars. To establish knowledge regarding the genetic variances and evolutional relationships of ploidy clones of Saccharum spontaneum collected in China, the rDNA-ITS sequences of 62 ploidy clones including octaploid clones (2n = 64), nonaploid clones (2n = 72), decaploid clones (2n = 80), and dodecaploid clones (2n = 96), were obtained and analyzed. The rDNA-ITS sequences of four species from Saccharum and Sorghum bicolor selected as controls. The results showed that decaploid clones (2n = 80) possess the most abundant variances with 58 variable sites and 20 parsim-informative sites in ITS sequences, which were then followed by octaploid clones with 43 variable sites and 17 parsim-informative sites. In haplotype diversity, all four population exhibited high diversity, especially nonaploid and decaploid populations. By comparing the genetic distances among four ploidy populations, the dodecaploid population exhibited the closest relationship with the nonaploid population, and then the relationship strength decreased successively for the decaploid population and then for the octaploid population. Population differentiation analysis showed that the phenomena of population differentiation were not found among different ploidy populations, and low coefficient of gene differentiation(Gst) and high gene flow(Nm) occur among these populations possessing close genetic relationship. These results mentioned above will contribute to the understanding of the evolution of different ploidy populations of Saccharum spontaneum and provide vital knowledge for their utilization in sugarcane breeding and innovation.

Positive selection of digestive Cys proteases in herbivorous Coleoptera

Positive selection is thought to contribute to the functional diversification of insect-inducible protease inhibitors in plants in response to selective pressures exerted by the digestive proteases of their herbivorous enemies. Here we assessed whether a reciprocal evolutionary process takes place on the insect side, and whether ingestion of a positively selected plant inhibitor may translate into a measurable rebalancing of midgut proteases in vivo. Midgut Cys proteases of herbivorous Coleoptera, including the major pest Colorado potato beetle (Leptinotarsa decemlineata), were first compared using a codon-based evolutionary model to look for the occurrence of hypervariable, positively selected amino acid sites among the tested sequences. Hypervariable sites were found, distributed within -or close to- amino acid regions interacting with Cys-type inhibitors of the plant cystatin protein family. A close examination of L. decemlineata sequences indicated a link between their assignment to protease functional families and amino acid identity at positively selected sites. A function-diversifying role for positive selection was further suggested empirically by in vitro protease assays and a shotgun proteomic analysis of L. decemlineata Cys proteases showing a differential rebalancing of protease functional family complements in larvae fed single variants of a model cystatin mutated at positively selected amino acid sites. These data confirm overall the occurrence of hypervariable, positively selected amino acid sites in herbivorous Coleoptera digestive Cys proteases. They also support the idea of an adaptive role for positive selection, useful to generate functionally diverse proteases in insect herbivores ingesting functionally diverse, rapidly evolving dietary cystatins.

Novel in vitro inhibitory functions of potato tuber proteinaceous inhibitors

Plant protease inhibitors are a structurally highly diverse and ubiquitous class of small proteins, which play various roles in plant development and defense against pests and pathogens. Particular isoforms inhibit in vitro proteases and other enzymes that are not their natural substrates, for example proteases that have roles in human diseases. Mature potato tubers are a rich source of several protease inhibitor families. Different cultivars have different inhibitor profiles. With the objective to explore the functional diversity of the natural diversity of potato protease inhibitors, we randomly selected and sequenced 9,600 cDNA clones originated from mature tubers of ten potato cultivars. Among these, 120 unique inhibitor cDNA clones were identified by homology searches. Eighty-eight inhibitors represented novel sequence variants of known plant protease inhibitor families. Most frequent were Kunitz-type inhibitors (KTI), potato protease inhibitors I and II (PIN), pectin methylesterase inhibitors, metallocarboxypeptidase inhibitors and defensins. Twenty-three inhibitors were functionally characterized after heterologous expression in the yeast Pichia pastoris. The purified recombinant proteins were tested for inhibitory activity on trypsin, eleven pharmacological relevant proteases and the non-proteolytic enzyme 5-lipoxygenase. Members of the KTI and PIN families inhibited pig pancreas elastase, β-Secretase, Cathepsin K, HIV-1 protease and potato 5-lipoxygenase. Our results demonstrate in vitro inhibitory diversity of small potato tuber proteins commonly known as protease inhibitors, which might have biotechnological or medical applications.

Defense response in non-genomic model species: methyl jasmonate exposure reveals the passion fruit leaves' ability to assemble a cocktail of functionally diversified Kunitz-type trypsin inhibitors and recruit two of them against papain

Multiplicity of protease inhibitors induced by predators may increase the understanding of a plant's intelligent behavior toward environmental challenges. Information about defense mechanisms of non-genomic model plant passion fruit (Passiflora edulis Sims) in response to predator attack is still limited. Here, via biochemical approaches, we showed its flexibility to build-up a broad repertoire of potent Kunitz-type trypsin inhibitors (KTIs) in response to methyl jasmonate. Seven inhibitors (20-25 kDa) were purified from exposed leaves by chromatographic techniques. Interestingly, the KTIs possessed truncated Kunitz motif in their N-terminus and some of them also presented non-consensus residues. Gelatin-Native-PAGE established multiple isoforms for each inhibitor. Significant differences regarding inhibitors' activity toward trypsin and chymotrypsin were observed, indicating functional polymorphism. Despite its rarity, two of them also inhibited papain, and such bifunctionality suggests a recruiting process onto another mechanistic class of target protease (cysteine-type). All inhibitors acted strongly on midgut proteases from sugarcane borer, Diatraea saccharalis (a lepidopteran insect) while in vivo assays supported their insecticide properties. Moreover, the bifunctional inhibitors displayed activity toward midgut proteases from cowpea weevil, Callosobruchus maculatus (a coleopteran insect). Unexpectedly, all inhibitors were highly effective against midgut proteases from Aedes aegypti a dipteran insect (vector of neglected tropical diseases) opening new avenues for plant-derived PIs for vector control-oriented research. Our results reflect the KTIs' complexities in passion fruit which could be wisely exploited by influencing plant defense conditions. Therefore, the potential of passion fruit as source of bioactive compounds with diversified biotechnological application was strengthened.