Selective histories of poplar protease inhibitors: elevated polymorphism, purifying selection, and positive selection driving divergence of recent duplicates

Poplar protease inhibitor expression differs in an herbivore specific manner

BACKGROUND

Protease inhibitors are defense proteins widely distributed in the plant kingdom. By reducing the activity of digestive enzymes in insect guts, they reduce the availability of nutrients and thus impair the growth and development of the attacking herbivore. One well-characterized class of protease inhibitors are Kunitz-type trypsin inhibitors (KTIs), which have been described in various plant species, including Populus spp. Long-lived woody perennials like poplar trees encounter a huge diversity of herbivores, but the specificity of tree defenses towards different herbivore species is hardly studied. We therefore aimed to investigate the induction of KTIs in black poplar (P. nigra) leaves upon herbivory by three different chewing herbivores, Lymantria dispar and Amata mogadorensis caterpillars, and Phratora vulgatissima beetles.

RESULTS

We identified and generated full-length cDNA sequences of 17 KTIs that are upregulated upon herbivory in black poplar leaves, and analyzed the expression patterns of the eight most up-regulated KTIs via qRT-PCR. We found that beetles elicited higher transcriptional induction of KTIs than caterpillars, and that both caterpillar species induced similar KTI expression levels. Furthermore, KTI expression strongly correlated with the trypsin-inhibiting activity in the herbivore-damaged leaves, but was not dependent on damage severity, i.e. leaf area loss, for most of the genes.

CONCLUSIONS

We conclude that the induction of KTIs in black poplar is controlled at the transcriptional level in a threshold-based manner and is strongly influenced by the species identity of the herbivore. However, the underlying molecular mechanisms and ecological consequences of these patterns remain to be investigated.

Recombinant cystatins in plants

Dozens of studies have assessed the practical value of plant cystatins as ectopic inhibitors of Cys proteases in biological systems. The potential of these proteins in crop protection to control herbivorous pests and pathogens has been documented extensively over the past 25 years. Their usefulness to regulate endogenous Cys proteases in planta has also been considered recently, notably to implement novel traits of agronomic relevance in crops or to generate protease activity-depleted environments in plants or plant cells used as bioreactors for recombinant proteins. After a brief update on the basic structural characteristics of plant cystatins, we summarize recent advances on the use of these proteins in plant biotechnology. Attention is also paid to the molecular improvement of their structural properties for the improvement of their protease inhibitory effects or the fine-tuning of their biological target range.

Idesia polycarpa (Salicaceae) leaf constituents and their toxic effect on Cerura vinula and Lymantria dispar (Lepidoptera) larvae

Phytochemical investigation of Idesia polycarpa (Salicaceae) resulted in the structure elucidation of nine previously undescribed phenolic natural products along with six known compounds. The compounds are structurally related to salicinoids that are known defense compounds from Salix and Populus species. The I. polycarpa diet was toxic, as shown in feeding experiments with larvae of Lymantria dispar, an herbivorous broadleaf tree generalist insect, and with larvae of Cerura vinula, a specialist adapted to poplar. The survival rate and mass gain of larvae was significantly lower when they fed on I. polycarpa leaves, compared to larvae fed on Populus nigra leaves. Potential reasons for the poor performance of both herbivores on I. polycarpa leaves are discussed.

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.

Prolonged Adaptive Evolution of a Defensive Gene in the Solanaceae

Although plants and their natural enemies may coevolve for prolonged periods, little is known about how long individual plant defensive genes are involved in the coevolutionary process. We address this issue by examining patterns of selection on the defensive gene threonine deaminase (TD). Tomato (Solanum lycopersicum) has two copies of this gene. One performs the canonical housekeeping function in amino acid metabolism of catalyzing the first reaction in the conversion of threonine to isoleucine. The second copy functions as an antinutritive defense against lepidopteran herbivores by depleting threonine in the insect gut. Wild tobacco (Nicotiana attenuata) also contains a defensive copy. We show that a single copy of TD underwent two or three duplications near the base of the Solanaceae. One copy retains the housekeeping function, whereas a second copy evolved defensive functions. Positive selection occurred on the branch of the TD2 gene tree subtending the common ancestor of the Nicotianoideae and Solanoideae. It also occurred within the Solanoideae clade but not within the Nicotianoideae clade. Finally, it occurred on most branches leading from the common ancestor to S. lycopersicum. Based on recent calibrations of the Solanaceae phylogeny, TD2 experienced adaptive substitutions for a period of 30-50 My. We suggest that the most likely explanation for this result is fluctuating herbivore abundances: When herbivores are rare, relaxed selection increases the likelihood that slightly disadvantageous mutations will be fixed by drift; when herbivores are common, increased selection causes the evolution of compensatory adaptive mutations. Alternative explanations are also discussed.

Influence of Genotype, Environment, and Gypsy Moth Herbivory on Local and Systemic Chemical Defenses in Trembling Aspen (Populus tremuloides)

Numerous studies have explored the impacts of intraspecific genetic variation and environment on the induction of plant chemical defenses by herbivory. Relatively few, however, have considered how those factors affect within-plant distribution of induced defenses. This work examined the impacts of plant genotype and soil nutrients on the local and systemic phytochemical responses of trembling aspen (Populus tremuloides) to defoliation by gypsy moth (Lymantria dispar). We deployed larvae onto foliage on individual tree branches for 15 days and then measured chemistry in leaves from: 1) branches receiving damage, 2) undamaged branches of insect-damaged trees, and 3) branches of undamaged control trees. The relationship between post-herbivory phytochemical variation and insect performance also was examined. Plant genotype, soil nutrients, and damage all influenced phytochemistry, with genotype and soil nutrients being stronger determinants than damage. Generally, insect damage decreased foliar nitrogen, increased levels of salicinoids and condensed tannins, but had little effect on levels of a Kunitz trypsin inhibitor, TI3. The largest damage-mediated tannin increases occurred in leaves on branches receiving damage, whereas the largest salicinoid increases occurred in leaves of adjacent, undamaged branches. Foliar nitrogen and the salicinoid tremulacin had the strongest positive and negative relationships, respectively, with insect growth. Overall, plant genetics and environment concomitantly influenced both local and systemic phytochemical responses to herbivory. These findings suggest that herbivory can contribute to phytochemical heterogeneity in aspen foliage, which may in turn influence future patterns of herbivory and nutrient cycling over larger spatial scales.

Physiological and proteomic responses of different willow clones (Salix fragilis x alba) exposed to dredged sediment contaminated by heavy metals

High biomass producing species are considered as tools for remediation of contaminated soils. Willows (Salix spp.) are prominent study subjects in this regard. In this study, different willow clones (Salix fragilis x alba) were planted on heavy-metal polluted dredging sludge. A first objective was assessment of the biomass production for these clones. Using a Gupta statistic, four clones were identified as high biomass producers (HBP). For comparison, a group of four clones with lowest biomass production were selected (LBP). A second objective was to compare metal uptake as well as the physiological and proteomic responses of these two groups. All these complementary data's allow us to have a better picture of the health of the clones that would be used in phytoremediation programs. Cd, Zn, and Ni total uptake was higher in the HBPs but Pb total uptake was higher in LBPs. Our proteomic and physiological results showed that the LBPs were able to maintain cellular activity as much as the HBPs although the oxidative stress response was more pronounced in the LBPs. This could be due to the high Pb content found in this group although a combined effect of the other metals cannot be excluded.

Gene duplication as a mechanism of genomic adaptation to a changing environment

A subject of extensive study in evolutionary theory has been the issue of how neutral, redundant copies can be maintained in the genome for long periods of time. Concurrently, examples of adaptive gene duplications to various environmental conditions in different species have been described. At this point, it is too early to tell whether or not a substantial fraction of gene copies have initially achieved fixation by positive selection for increased dosage. Nevertheless, enough examples have accumulated in the literature that such a possibility should be considered. Here, I review the recent examples of adaptive gene duplications and make an attempt to draw generalizations on what types of genes may be particularly prone to be selected for under certain environmental conditions. The identification of copy-number variation in ecological field studies of species adapting to stressful or novel environmental conditions may improve our understanding of gene duplications as a mechanism of adaptation and its relevance to the long-term persistence of gene duplications.

Geographical structure and adaptive population differentiation in herbivore defence genes in European aspen (Populus tremula L., Salicaceae)

When a phenotypic trait is subjected to spatially variable selection and local adaptation, the underlying genes controlling the trait are also expected to show strong patterns of genetic differentiation because alternative alleles are favoured in different geographical locations. Here, we study 71 single nucleotide polymorphisms (SNPs) from seven genes associated with inducible defence responses in a sample of Populus tremula collected from across Sweden. Four of these genes (PPO2, TI2, TI4 and TI5) show substantial population differentiation, and a principal component analyses conducted on the defence SNPs divides the Swedish population into three distinct clusters. Several defence SNPs show latitudinal clines, although these were not robust to multiple testing. However, five SNPs (located within TI4 and TI5) show strong longitudinal clines that remain significant after multiple test correction. Genetic geographical variation, supporting local adaptation, has earlier been confirmed in genes involved in the photoperiod pathway in P. tremula, but this is, to our knowledge, one of the first times that geographical variation has been found in genes involved in plant defence against antagonists.

Pleistocene speciation in the genus Populus (salicaceae)

The macroevolutionary consequences of recent climate change remain controversial, and there is little paleobotanical or morphological evidence that Pleistocene (1.8-0.12 Ma) glacial cycles acted as drivers of speciation, especially among lineages with long generation times, such as trees. We combined genetic and ecogeographic data from 2 closely related North American tree species, Populus balsamifera and P. trichocarpa (Salicacaeae), to determine if their divergence coincided with and was possibly caused by Pleistocene climatic events. We analyzed 32 nuclear loci from individuals of P. balsamifera and P. trichocarpa to produce coalescent-based estimates of the divergence time between the 2 species. We coupled the coalescent analyses with paleodistribution models to assess the influence of climate change on species' range. Furthermore, measures of niche overlap were used to investigate patterns of ecological differentiation between species. We estimated the divergence date of P. balsamifera and P. trichocarpa at approximately 75 Ka, which corresponds closely with the onset of Marine Isotope Stage 4 (∼76 Ka) and a rapid increase in global ice volume. Significance tests of niche overlap, in conjunction with genetic estimates of migration, suggested that speciation occurred in allopatry, possibly resulting from the environmental effects of Pleistocene glacial cycles. Our results indicate that the divergence of keystone tree species, which have shaped community diversity in northern North American ecosystems, was recent and may have been a consequence of Pleistocene-era glaciation and climate change.

Coexistence of trichome variation in a natural plant population: a combined study using ecological and candidate gene approaches

The coexistence of distinct phenotypes within populations has long been investigated in evolutionary ecology. Recent studies have identified the genetic basis of distinct phenotypes, but it is poorly understood how the variation in candidate loci is maintained in natural environments. In this study, we examined fitness consequences and genetic basis of variation in trichome production in a natural population of Arabidopsis halleri subsp. gemmifera. Half of the individuals in the study population produced trichomes while the other half were glabrous, and the leaf beetle Phaedon brassicae imposed intensive damage to both phenotypes. The fitness of hairy and glabrous plants showed no significant differences in the field during two years. A similar result was obtained when sibling hairy and glabrous plants were transplanted at the same field site, whereas a fitness cost of trichome production was detected under a weak herbivory condition. Thus, equivalent fitness of hairy and glabrous plants under natural herbivory allows their coexistence in the contemporary population. The pattern of polymorphism of the candidate trichome gene GLABROUS1 (GL1) showed no evidence of long-term maintenance of trichome variation within the population. Although balancing selection under fluctuating biotic environments is often proposed to explain the maintenance of defense variation, the lack of clear evidence of balancing selection in the study population suggests that other factors such as gene flow and neutral process may have played relatively large roles in shaping trichome variation at least for the single population level.

Adaptation to drought in two wild tomato species: the evolution of the Asr gene family

Wild tomato species are a valuable system in which to study local adaptation to drought: they grow in diverse environments ranging from mesic to extremely arid conditions. Here, we investigate the evolution of members of the Asr (ABA/water stress/ripening induced) gene family, which have been reported to be involved in the water stress response. We analysed molecular variation in the Asr gene family in populations of two closely related species, Solanum chilense and Solanum peruvianum. We concluded that Asr1 has evolved under strong purifying selection. In contrast to previous reports, we did not detect evidence for positive selection at Asr2. However, Asr4 shows patterns consistent with local adaptation in an S. chilense population that lives in an extremely dry environment. We also discovered a new member of the gene family, Asr5. Our results show that the Asr genes constitute a dynamic gene family and provide an excellent example of tandemly arrayed genes that are of importance in adaptation. Taking the potential distribution of the species into account, it appears that S. peruvianum can cope with a great variety of environmental conditions without undergoing local adaptation, whereas S. chilense undergoes local adaptation more frequently.

Local selection across a latitudinal gradient shapes nucleotide diversity in balsam poplar, Populus balsamifera L

Molecular studies of adaptive evolution often focus on detecting selective sweeps driven by positive selection on a species-wide scale; however, much adaptation is local, particularly of ecologically important traits. Here, we look for evidence of range-wide and local adaptation at candidate genes for adaptive phenology in balsam poplar, Populus balsamifera, a widespread forest tree whose range extends across environmental gradients of photoperiod and growing season length. We examined nucleotide diversity of 27 poplar homologs of the flowering-time network-a group of genes that control plant developmental phenology through interactions with environmental cues such as photoperiod and temperature. Only one gene, ZTL2, showed evidence of reduced diversity and an excess of fixed replacement sites, consistent with a species-wide selective sweep. Two other genes, LFY and FRI, harbored high levels of nucleotide diversity and exhibited elevated differentiation between northern and southern accessions, suggesting local adaptation along a latitudinal gradient. Interestingly, FRI has also been identified as a target of local selection between northern and southern accessions of Arabidopsis thaliana, indicating that this gene may be commonly involved in ecological adaptation in distantly related species. Our findings suggest an important role for local selection shaping molecular diversity and reveal limitations of inferring molecular adaptation from analyses designed only to detect species-wide selective sweeps.

Nucleotide diversity and linkage disequilibrium in balsam poplar (Populus balsamifera)

*Current perceptions that poplars have high levels of nucleotide variation, large effective population sizes, and rapid decay of linkage disequilibrium are based primarily on studies from one poplar species, Populus tremula. *We analysed 590 gene fragments (average length 565 bp) from each of 15 individuals from different populations from throughout the range of Populus balsamifera. *Nucleotide diversity (theta(total) = 0.0028, pi = 0.0027) was low compared with other trees and model agricultural systems. Patterns of nucleotide diversity and site frequency spectra were consistent with purifying selection on replacement and intron sites. When averaged across all loci we found no evidence for decay of linkage disequilibrium across 750 bp, consistent with the low estimates of the scaled recombination parameter, rho = 0.0092. *Compared with P. tremula, a well studied congener with a similar distribution, P. balsamifera has low diversity and low effective recombination, both of which indicate a lower effective population size in P. balsamifera. Patterns of diversity and linkage indicate that there is considerable variation in population genomic patterns among poplar species and unlike P. tremula, association mapping techniques in balsam poplar should consider sampling single nucleotide polymorphisms (SNPs) at well-spaced intervals.