Evidence of Coevolution Between Cronartium harknessii Lineages and Their Corresponding Hosts, Lodgepole Pine and Jack Pine

Variation in rate of infection and susceptibility of Pinus spp. to the fungus Cronartium harknessii (syn. Endocronartium harknessii), the causative agent of western gall rust, has been well documented. To test the hypothesis that there is a coevolutionary relationship between C. harknessii and its hosts, we examined genetic structure and virulence of C. harknessii associated with lodgepole pine (P. contorta var. latifolia), jack pine (P. banksiana), and their hybrids. A secondary objective was to improve assessment and diagnosis of infection in hosts. Using 18 microsatellites, we assessed genetic structure of C. harknessii from 90 sites within the ranges of lodgepole pine and jack pine. We identified two lineages (East and West, F_ST = 0.677) associated with host genetic structure (r = 0.81, P = 0.001), with East comprising three sublineages. In parallel, we conducted a factorial experiment in which lodgepole pine, jack pine, and hybrid seedlings were inoculated with spores from the two primary genetic lineages. With this experiment, we refined the phenotypic categories associated with infection and demonstrated that stem width can be used as a quantitative measure of host response to infection. Overall, each host responded differentially to the fungal lineages, with jack pine exhibiting more resiliency to infection than lodgepole pine and hybrids exhibiting intermediate resiliency. Taken together, the shared genetic structure between fungus and host species, and the differential interaction of the fungal species with the hosts, supports a coevolutionary relationship between host and pathogen.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Linking genotype to phenotype to identify genetic variation relating to host susceptibility in the mountain pine beetle system

Identifying genetic variants responsible for phenotypic variation under selective pressure has the potential to enable productive gains in natural resource conservation and management. Despite this potential, identifying adaptive candidate loci is not trivial, and linking genotype to phenotype is a major challenge in contemporary genetics. Many of the population genetic approaches commonly used to identify adaptive candidates will simultaneously detect false positives, particularly in nonmodel species, where experimental evidence is seldom provided for putative roles of the adaptive candidates identified by outlier approaches. In this study, we use outcomes from population genetics, phenotype association, and gene expression analyses as multiple lines of evidence to validate candidate genes. Using lodgepole and jack pine as our nonmodel study species, we analyzed 17 adaptive candidate loci together with 78 putatively neutral loci at 58 locations across Canada (N > 800) to determine whether relationships could be established between these candidate loci and phenotype related to mountain pine beetle susceptibility. We identified two candidate loci that were significant across all population genetic tests, and demonstrated significant changes in transcript abundance in trees subjected to wounding or inoculation with the mountain pine beetle fungal associate Grosmannia clavigera. Both candidates are involved in central physiological processes that are likely to be invoked in a trees response to stress. One of these two candidate loci showed a significant association with mountain pine beetle attack status in lodgepole pine. The spatial distribution of the attack-associated allele further coincides with other indicators of susceptibility in lodgepole pine. These analyses, in which population genetics was combined with laboratory and field experimental validation approaches, represent first steps toward linking genetic variation to the phenotype of mountain pine beetle susceptibility in lodgepole and jack pine, and provide a roadmap for more comprehensive analyses.

A novel application of RNase H2-dependent quantitative PCR for detection and quantification of Grosmannia clavigera, a mountain pine beetle fungal symbiont, in environmental samples

Mountain pine beetle (Dendroctonus ponderosae Hopkins; MPB) is an economically and ecologically important pest of pine species in western North America. Mountain pine beetles form complex multipartite relationships with microbial partners, including the ophiostomoid fungi Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield, Ophiostoma montium (Rumbold) von Arx, Grosmannia aurea (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield, Leptographium longiclavatum (Lee, Kim, and Breuil) and Leptographium terebrantis (Barras and Perry). These fungi are vectored by MPB to new pine hosts, where the fungi overcome host defenses to grow into the sapwood. A tree's relative susceptibility to these fungi is conventionally assessed by measuring lesions that develop in response to fungal inoculation. However, these lesions represent a symptom of infection, representing both fungal growth and tree defense capacity. In order to more objectively assess fungal virulence and host tree susceptibility in studies of host-pathogen interactions, a reliable, consistent, sensitive method is required to accurately identify and quantify MPB-associated fungal symbionts in planta. We have adapted RNase H2-dependent PCR, a technique originally designed for rare allele discrimination, to develop a novel RNase H2-dependent quantitative PCR (rh-qPCR) assay that shows greater specificity and sensitivity than previously published PCR-based methods to quantify MPB fungal symbionts in pine xylem and MPB whole beetles. Two sets of assay probes were designed: one that amplifies a broad range of ophiostomoid species, and a second that amplifies G. clavigera but not other MPB-associated ophiostomoid species. Using these primers to quantify G. clavigera in pine stems, we provide evidence that lesion length does not accurately reflect the extent of fungal colonization along the stem nor the quantity of fungal growth within this colonized portion of stem. The sensitivity, specificity, reproducibility, cost effectiveness and high-throughput potential of the rh-qPCR assay makes the technology suitable for identification and quantification of a wide array of pathogenic and beneficial microbes that form associations with plants and other organisms, even when the microbial partner is present in low abundance.

Genetic and genomic evidence of niche partitioning and adaptive radiation in mountain pine beetle fungal symbionts

Bark beetles form multipartite symbiotic associations with blue stain fungi (Ophiostomatales, Ascomycota). These fungal symbionts play an important role during the beetle's life cycle by providing nutritional supplementation, overcoming tree defences and modifying host tissues to favour brood development. The maintenance of stable multipartite symbioses with seemingly less competitive symbionts in similar habitats is of fundamental interest to ecology and evolution. We tested the hypothesis that the coexistence of three fungal species associated with the mountain pine beetle is the result of niche partitioning and adaptive radiation using SNP genotyping coupled with genotype-environment association analysis and phenotypic characterization of growth rate under different temperatures. We found that genetic variation and population structure within each species is best explained by distinct spatial and environmental variables. We observed both common (temperature seasonality and the host species) and distinct (drought, cold stress, precipitation) environmental and spatial factors that shaped the genomes of these fungi resulting in contrasting outcomes. Phenotypic intraspecific variations in Grosmannia clavigera and Leptographium longiclavatum, together with high heritability, suggest potential for adaptive selection in these species. By contrast, Ophiostoma montium displayed narrower intraspecific variation but greater tolerance to extreme high temperatures. Our study highlights unique phenotypic and genotypic characteristics in these symbionts that are consistent with our hypothesis. By maintaining this multipartite relationship, the bark beetles have a greater likelihood of obtaining the benefits afforded by the fungi and reduce the risk of being left aposymbiotic. Complementarity among species could facilitate colonization of new habitats and survival under adverse conditions.

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land-use change have affected contemporary range-wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high-throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled 'Genomics and Forest Tree Genetics' was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome-enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.

Change in soil fungal community structure driven by a decline in ectomycorrhizal fungi following a mountain pine beetle (Dendroctonus ponderosae) outbreak

Western North American landscapes are rapidly being transformed by forest die-off caused by mountain pine beetle (Dendroctonus ponderosae), with implications for plant and soil communities. The mechanisms that drive changes in soil community structure, particularly for the highly prevalent ectomycorrhizal fungi in pine forests, are complex and intertwined. Critical to enhancing understanding will be disentangling the relative importance of host tree mortality from changes in soil chemistry following tree death. Here, we used a recent bark beetle outbreak in lodgepole pine (Pinus contorta) forests of western Canada to test whether the effects of tree mortality altered the richness and composition of belowground fungal communities, including ectomycorrhizal and saprotrophic fungi. We also determined the effects of environmental factors (i.e. soil nutrients, moisture, and phenolics) and geographical distance, both of which can influence the richness and composition of soil fungi. The richness of both groups of soil fungi declined and the overall composition was altered by beetle-induced tree mortality. Soil nutrients, soil phenolics and geographical distance influenced the community structure of soil fungi; however, the relative importance of these factors differed between ectomycorrhizal and saprotrophic fungi. The independent effects of tree mortality, soil phenolics and geographical distance influenced the community composition of ectomycorrhizal fungi, while the community composition of saprotrophic fungi was weakly but significantly correlated with the geographical distance of plots. Taken together, our results indicate that both deterministic and stochastic processes structure soil fungal communities following landscape-scale insect outbreaks and reflect the independent roles tree mortality, soil chemistry and geographical distance play in regulating the community composition of soil fungi.

The joint influence of photoperiod and temperature during growth cessation and development of dormancy in white spruce (Picea glauca)

Timely responses to environmental cues enable the synchronization of phenological life-history transitions essential for the health and survival of north-temperate and boreal tree species. While photoperiodic cues will remain persistent under climate change, temperature cues may vary, contributing to possible asynchrony in signals influencing developmental and physiological transitions essential to forest health. Understanding the relative contribution of photoperiod and temperature as determinants of the transition from active growth to dormancy is important for informing adaptive forest management decisions that consider future climates. Using a combination of photoperiod (long = 20 h or short = 8 h day lengths) and temperature (warm = 22 °C/16 °C and cool = 8 °C/4 °C day/night, respectively) treatments, we used microscopy, physiology and modeling to comprehensively examine hallmark traits of the growth-dormancy transition-including bud formation, growth cessation, cold hardiness and gas exchange-within two provenances of white spruce [Picea glauca (Moench) Voss] spanning a broad latitude in Alberta, Canada. Following exposure to experimental treatments, seedlings were transferred to favorable conditions, and the depth of dormancy was assessed by determining the timing and ability of spruce seedlings to resume growth. Short photoperiods promoted bud development and growth cessation, whereas longer photoperiods extended the growing season through the induction of lammas growth. In contrast, cool temperatures under both photoperiodic conditions delayed bud development. Photoperiod strongly predicted the development of cold hardiness, whereas temperature predicted photosynthetic rates associated with active growth. White spruce was capable of attaining endodormancy, but its release was environmentally determined. Dormancy depth varied substantially across experimental treatments suggesting that environmental cues experienced within one season could affect growth in the following season, which is particularly important for a determinate species such as white spruce. The joint influence of these environmental cues points toward the importance of including local constant photoperiod and shifting temperature cues into predictive models that consider how climate change may affect northern forests.

Hydraulic conductivity and aquaporin transcription in roots of trembling aspen (Populus tremuloides) seedlings colonized by Laccaria bicolor

Ectomycorrhizal fungi have been reported to increase root hydraulic conductivity (L pr) by altering apoplastic and plasma membrane intrinsic protein (PIP)-mediated cell-to-cell water transport pathways in associated roots, or to have little effect on root water transport, depending on the interacting species and imposed stresses. In this study, we investigated the water transport properties and PIP transcription in roots of aspen (Populus tremuloides) seedlings colonized by the wild-type strain of Laccaria bicolor and by strains overexpressing a major fungal water-transporting aquaporin JQ585595. Inoculation of aspen seedlings with L. bicolor resulted in about 30 % colonization rate of root tips, which developed dense mantle and the Hartig net that was restricted in the modified root epidermis. Transcript abundance of the aspen aquaporins PIP1;2, PIP2;1, and PIP2;2 decreased in colonized root tips. Root colonization by JQ585595-overexpressing strains had no significant impact on seedling shoot water potentials, gas exchange, or dry mass; however, it led to further decrease in transcript abundance of PIP1;2 and PIP2;3 and the significantly lower L pr than in non-inoculated roots. These results, taken together with our previous study that showed enhanced root water hydraulics of L. bicolor-colonized white spruce (Picea glauca), suggest that the impact of L. bicolor on root hydraulics varies by the ectomycorrhiza-associated tree species.

Responses of hybrid aspen over-expressing a PIP2;5 aquaporin to low root temperature

Aquaporins mediate the movement of water across cell membranes. Plasma membrane intrinsic protein 2;5 from Populus trichocarpa×deltoides (PtdPIP2;5) was previously demonstrated to be a functionally important water conducting aquaporin. To study the relevance of aquaporin-mediated root water transport at low temperatures, we generated transgenic Populus tremula×alba over-expressing PtdPIP2;5 under control of the maize ubiquitin promoter, and compared the physiological responses and water transport properties of the PtdPIP2;5 over-expressing lines (PtdPIP2;5ox) with wild-type plants. We hypothesized that over-expression of PtdPIP2;5 would reduce temperature sensitivity of root water transport and gas exchange. Decreasing root temperatures to 10 and 5°C significantly decreased hydraulic conductivities (Lp) in wild-type plants, but had no significant effect on Lp in PtdPIP2;5ox plants. Recovery of Lp in the transgenic lines returned to 20°C from 5°C was faster than in the wild-type plants. Low root temperature did not induce major changes in transcript levels for other PIPs. When roots were exposed to 5°C in solution culture and shoots were exposed to 20°C, wild-type plants had significantly lower net photosynthetic and transpiration rates compared to PtdPIP2;5ox plants. Taken together, our results demonstrate that over-expression of PtdPIP2;5 in P. tremula×alba was effective in alleviating the effects of low root temperature on Lp and gas exchange.

Transcript profiling of aquaporins during basidiocarp development in Laccaria bicolor ectomycorrhizal with Picea glauca

Sporocarp formation is part of the reproductive stage in the life cycle of many mycorrhizal macrofungi. Sporocarp formation is accompanied by a transcriptomic switch and profound changes in regulation of the gene families that play crucial roles in the sporocarp initiation and maturation. Since sporocarp growth requires efficient water delivery, in the present study, we investigated changes in transcript abundance of six fungal aquaporin genes that could be cloned from the ectomycorrhizal fungus Laccaria bicolor strain UAMH8232, during the initiation and development of its basidiocarp. Aquaporins are intrinsic membrane proteins facilitating the transmembrane transport of water and other small neutral molecules. In controlled-environment experiments, we induced basidiocarp formation in L. bicolor, which formed ectomycorrhizal associations with white spruce (Picea glauca) seedlings. We profiled transcript abundance corresponding to six fungal aquaporin genes at six different developmental stages of basidiocarp growth and development. We also compared physiological parameters of non-inoculated to mycorrhizal seedlings with and without the presence of basidiocarps. Two L. bicolor aquaporins--JQ585592, a functional channel for CO2, NO and H2O2, and JQ585595, a functional water channel--showed the greatest degree of upregulation during development of the basidiocarp. Our findings point to the importance of aquaporin-mediated transmembrane water and CO2 transport during distinct stages of basidiocarp development.

Overexpression of Laccaria bicolor aquaporin JQ585595 alters root water transport properties in ectomycorrhizal white spruce (Picea glauca) seedlings

The contribution of hyphae to water transport in ectomycorrhizal (ECM) white spruce (Picea glauca) seedlings was examined by altering expression of a major water-transporting aquaporin in Laccaria bicolor. Picea glauca was inoculated with wild-type (WT), mock transgenic or L. bicolor aquaporin JQ585595-overexpressing (OE) strains and exposed to root temperatures ranging from 5 to 20°C to examine the root water transport properties, physiological responses and plasma membrane intrinsic protein (PIP) expression in colonized plants. Mycorrhization increased shoot water potential, transpiration, net photosynthetic rates, root hydraulic conductivity and root cortical cell hydraulic conductivity in seedlings. At 20°C, OE plants had higher root hydraulic conductivity compared with WT plants and the increases were accompanied by higher expression of P. glauca PIP GQ03401_M18.1 in roots. In contrast to WT L. bicolor, the effects of OE fungi on root and root cortical cell hydraulic conductivities were abolished at 10 and 5°C in the absence of major changes in the examined transcript levels of P. glauca root PIPs. The results provide evidence for the importance of fungal aquaporins in root water transport of mycorrhizal plants. They also demonstrate links between hyphal water transport, root aquaporin expression and root water transport in ECM plants.

Influence of water deficit on the molecular responses of Pinus contorta × Pinus banksiana mature trees to infection by the mountain pine beetle fungal associate, Grosmannia clavigera

Conifers exhibit a number of constitutive and induced mechanisms to defend against attack by pests and pathogens such as mountain pine beetle (Dendroctonus ponderosae Hopkins) and their fungal associates. Ecological studies have demonstrated that stressed trees are more susceptible to attack by mountain pine beetle than their healthy counterparts. In this study, we tested the hypothesis that water deficit affects constitutive and induced responses of mature lodgepole pine × jack pine hybrids (Pinus contorta Dougl. ex Loud. var. latifolia Engelm. ex S. Wats. × Pinus banksiana Lamb.) to inoculation with the mountain pine beetle fungal associate Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield. The degree of stress induced by the imposed water-deficit treatment was sufficient to reduce photosynthesis. Grosmannia clavigera-induced lesions exhibited significantly reduced dimensions in water-deficit trees relative to well-watered trees at 5 weeks after inoculation. Treatment-associated cellular-level changes in secondary phloem were also observed. Quantitative RT-PCR was used to analyze transcript abundance profiles of 18 genes belonging to four families classically associated with biotic and abiotic stress responses: aquaporins (AQPs), dehydration-responsive element binding (DREB), terpene synthases (TPSs) and chitinases (CHIs). Transcript abundance profiles of a TIP2 AQP and a TINY-like DREB decreased significantly in fungus-inoculated trees, but not in response to water deficit. One TPS, Pcb(+)-3-carene synthase, and the Class II CHIs PcbCHI2.1 and PcbCHI2.2 showed increased expression under water-deficit conditions in the absence of fungal inoculation, while another TPS, Pcb(E)-β-farnesene synthase-like, and two CHIs, PcbCHI1.1 and PcbCHI4.1, showed attenuated expression under water-deficit conditions in the presence of fungal inoculation. The effects were observed both locally and systemically. These results demonstrate that both constitutive and induced carbon- and nitrogen-based defenses are affected by water deficit, suggesting potential consequences for mountain pine beetle dynamics, particularly in novel environments.

Cross-species outlier detection reveals different evolutionary pressures between sister species

Lodgepole pine (Pinus contorta var. latifolia) and jack pine (Pinus banksiana) hybridize in western Canada, an area of recent mountain pine beetle range expansion. Given the heterogeneity of the environment, and indications of local adaptation, there are many unknowns regarding the response of these forests to future outbreaks. To better understand this we aim to identify genetic regions that have adaptive potential. We used data collected on 472 single nucleotide polymorphism (SNP) loci from 576 tree samples collected across 13 lodgepole pine-dominated sites and four jack pine-dominated sites. We looked at the relationship of genetic diversity with the environment, and we identified candidate loci using both frequency-based (arlequin and bayescan) and correlation-based (matsam and bayenv) methods. We found contrasting relationships between environmental variation and genetic diversity for the species. While we identified a number of candidate outliers (34 in lodgepole pine, 25 in jack pine, and 43 interspecific loci), we did not find any loci in common between lodgepole and jack pine. Many of the outlier loci identified were correlated with environmental variation. Using rigorous criteria we have been able to identify potential outlier SNPs. We have also found evidence of contrasting environmental adaptations between lodgepole and jack pine which could have implications for beetle spread risk.

Population structure of mountain pine beetle symbiont Leptographium longiclavatum and the implication on the multipartite beetle-fungi relationships

Over 18 million ha of forests have been destroyed in the past decade in Canada by the mountain pine beetle (MPB) and its fungal symbionts. Understanding their population dynamics is critical to improving modeling of beetle epidemics and providing potential clues to predict population expansion. Leptographium longiclavatum and Grosmannia clavigera are fungal symbionts of MPB that aid the beetle to colonize and kill their pine hosts. We investigated the genetic structure and demographic expansion of L. longiclavatum in populations established within the historic distribution range and in the newly colonized regions. We identified three genetic clusters/populations that coincide with independent geographic locations. The genetic profiles of the recently established populations in northern British Columbia (BC) and Alberta suggest that they originated from central and southern BC. Approximate Bayesian Computation supports the scenario that this recent expansion represents an admixture of individuals originating from BC and the Rocky Mountains. Highly significant correlations were found among genetic distance matrices of L. longiclavatum, G. clavigera, and MPB. This highlights the concordance of demographic processes in these interacting organisms sharing a highly specialized niche and supports the hypothesis of long-term multipartite beetle-fungus co-evolutionary history and mutualistic relationships.

Decline of ectomycorrhizal fungi following a mountain pine beetle epidemic

Forest die-off caused by mountain pine beetle (MPB; Dendroctonus ponderosa) is rapidly transforming western North American landscapes. The rapid and widespread death of lodgepole pine (Pinus contorta) will likely have cascading effects on biodiversity. One group particularly prone to such declines associated with MPB are ectomycorrhizal fungi, symbiotic organisms that can depend on pine for their survival, and are critical for stand regeneration. We evaluated the indirect effects of MPB on above- (community composition of epigeous sporocarps) and belowground (hyphal abundance) occurrences of ectomycorrhizal fungi across 11 forest stands. Along a gradient of mortality (0-82% pine killed), macromycete community composition changed; this shift was driven by a decrease in the species richness of ectomycorrhizal fungi. Both the proportion of species that were ectomycorrhizal and hyphal length in the soil declined with increased MPB-caused pine mortality; < 10% of sporocarp species were ectomycorrhizal in stands with high pine mortality compared with > 70% in stands without MPB attacks. The rapid range expansion of a native insect results not only in the widespread mortality of an ecologically and economically important pine species, but the effect of MPB may also be exacerbated by the concomitant decline of fungi crucial for recovery of these forests.

How the mountain pine beetle (Dendroctonus ponderosae) breached the Canadian Rocky Mountains

The mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins), a major pine forest pest native to western North America, has extended its range north and eastward during an ongoing outbreak. Determining how the MPB has expanded its range to breach putative barriers, whether physical (nonforested prairie and high elevation of the Rocky Mountains) or climatic (extreme continental climate where temperatures can be below −40 °C), may contribute to our general understanding of range changes as well as management of the current epidemic. Here, we use a panel of 1,536 single nucleotide polymorphisms (SNPs) to assess population genetic structure, connectivity, and signals of selection within this MPB range expansion. Biallelic SNPs in MPB from southwestern Canada revealed higher genetic differentiation and lower genetic connectivity than in the northern part of its range. A total of 208 unique SNPs were identified using different outlier detection tests, of which 32 returned annotations for products with putative functions in cholesterol synthesis, actin filament contraction, and membrane transport. We suggest that MPB has been able to spread beyond its previous range by adjusting its cellular and metabolic functions, with genome scale differentiation enabling populations to better withstand cooler climates and facilitate longer dispersal distances. Our study is the first to assess landscape-wide selective adaptation in an insect. We have shown that interrogation of genomic resources can identify shifts in genetic diversity and putative adaptive signals in this forest pest species.

Multilocus species identification and fungal DNA barcoding: insights from blue stain fungal symbionts of the mountain pine beetle

There is strong community-wide interest in applying molecular techniques to fungal species delimitation and identification, but selection of a standardized region or regions of the genome has not been finalized. A single marker, the ribosomal DNA internal transcribed spacer region, has frequently been suggested as the standard for fungi. We used a group of closely related blue stain fungi associated with the mountain pine beetle (Dendroctonus ponderosae Hopkins) to examine the success of such single-locus species identification, comparing the internal transcribed spacer with four other nuclear markers. We demonstrate that single loci varied in their utility for identifying the six fungal species examined, while use of multiple loci was consistently successful. In a literature survey of 21 similar studies, individual loci were also highly variable in their ability to provide consistent species identifications and were less successful than multilocus diagnostics. Accurate species identification is the essence of any molecular diagnostic system, and this consideration should be central to locus selection. Moreover, our study and the literature survey demonstrate the value of using closely related species as the proving ground for developing a molecular identification system. We advocate use of a multilocus barcode approach that is similar to the practice employed by the plant barcode community, rather than reliance on a single locus.

The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms

In tree species native to temperate and boreal regions, the activity-dormancy cycle is an important adaptive trait both for survival and growth. We discuss recent research on mechanisms controlling the overlapping developmental processes that define the activity-dormancy cycle, including cessation of apical growth, bud development, induction, maintenance and release of dormancy, and bud burst. The cycle involves an extensive reconfiguration of metabolism. Environmental control of the activity-dormancy cycle is based on perception of photoperiodic and temperature signals, reflecting adaptation to prevailing climatic conditions. Several molecular actors for control of growth cessation have been identified, with the CO/FT regulatory network and circadian clock having important coordinating roles in control of growth and dormancy. Other candidate regulators of bud set, dormancy and bud burst have been identified, such as dormancy-associated MADS-box factors, but their exact roles remain to be discovered. Epigenetic mechanisms also appear to factor in control of the activity-dormancy cycle. Despite evidence for gibberellins as negative regulators in growth cessation, and ABA and ethylene in bud formation, understanding of the roles that plant growth regulators play in controlling the activity-dormancy cycle is still very fragmentary. Finally, some of the challenges for further research in bud dormancy are discussed.

Characterizing the physical and genetic structure of the lodgepole pine × jack pine hybrid zone: mosaic structure and differential introgression

Understanding the physical and genetic structure of hybrid zones can illuminate factors affecting their formation and stability. In north-central Alberta, lodgepole pine (Pinus contorta Dougl. ex Loud. var. latifolia) and jack pine (Pinus banksiana Lamb) form a complex and poorly defined hybrid zone. Better knowledge of this zone is relevant, given the recent host expansion of mountain pine beetle into jack pine. We characterized the zone by genotyping 1998 lodgepole, jack pine, and hybrids from British Columbia, Alberta, Saskatchewan, Ontario, and Minnesota at 11 microsatellites. Using Bayesian algorithms, we calculated genetic ancestry and used this to model the relationship between species occurrence and environment. In addition, we analyzed the ancestry of hybrids to calculate the genetic contribution of lodgepole and jack pine. Finally, we measured the amount of gene flow between the pure species. We found the distribution of the pine classes is explained by environmental variables, and these distributions differ from classic distribution maps. Hybrid ancestry was biased toward lodgepole pine; however, gene flow between the two species was equal. The results of this study suggest that the hybrid zone is complex and influenced by environmental constraints. As a result of this analysis, range limits should be redefined.

Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal

Environmental change has a wide range of ecological consequences, including species extinction and range expansion. Many studies have shown that insect species respond rapidly to climatic change. A mountain pine beetle epidemic of record size in North America has led to unprecedented mortality of lodgepole pine, and a significant range expansion to the northeast of its historic range. Our goal was to determine the spatial genetic variation found among outbreak population from which genetic structure, and dispersal patterns may be inferred. Beetles from 49 sampling locations throughout the outbreak area in western Canada were analysed at 13 microsatellite loci. We found significant north-south population structure as evidenced by: (i) Bayesian-based analyses, (ii) north-south genetic relationships and diversity gradients; and (iii) a lack of isolation-by-distance in the northernmost cluster. The north-south structure is proposed to have arisen from the processes of postglacial colonization as well as recent climate-driven changes in population dynamics. Our data support the hypothesis of multiple sources of origin for the outbreak and point to the need for population specific information to improve our understanding and management of outbreaks. The recent range expansion across the Rocky Mountains into the jack/lodgepole hybrid and pure jack pine zones of northern Alberta is consistent with a northern British Columbia origin. We detected no loss of genetic variability in these populations, indicating that the evolutionary potential of mountain pine beetle to adapt has not been reduced by founder events. This study illustrates a rapid range-wide response to the removal of climatic constraints, and the potential for range expansion of a regional population.

Root herbivory: molecular analysis of the maize transcriptome upon infestation by Southern corn rootworm, Diabrotica undecimpunctata howardi

While many studies have characterized changes to the transcriptome of plants attacked by shoot-eating insect pests, few have examined transcriptome-level effects of root pests. Maize (Zea mays) seedlings were subjected to infestation for approximately 2 weeks by the root herbivore southern corn rootworm (SCR) Diabrotica undecimpunctata howardi, and changes in transcript abundance within both roots and shoots were analyzed using a 57K element microarray. A total of 541 genes showed statistically significant changes in transcript abundance in infested roots, including genes encoding many pathogenesis-related proteins such as chitinases, proteinase inhibitors, peroxidases and β-1,3-glucanases. Several WRKY transcription factors--often associated with biotic responses--exhibited increased transcript abundance upon SCR feeding. Differentially expressed (DE) genes were also detected in shoots of infested vs control plants. Quantitative Reverse Transcriptase Polymerase Chain Reaction (RT-PCR) was used to confirm patterns of transcript abundance for several significant DE genes using an independent experiment with a 2-6 day period of SCR infestation. Because of the well-documented roles that jasmonic acid (JA) or salicylic acid (SA) play in herbivory responses, the effect of exogenous JA or SA application on transcript abundance corresponding to the same subset of SCR-responsive genes was assessed. The response of these genes at the level of transcript abundance to SA and JA differed between roots and shoots and also differed among the genes that were examined. These data suggested that SA- and JA-dependent and independent signals contributed to the transcriptome-level changes in maize roots and shoots in response to SCR infestation.

Integrated transcriptomic and proteomic profiling of white spruce stems during the transition from active growth to dormancy

In the autumn, stems of woody perennials such as forest trees undergo a transition from active growth to dormancy. We used microarray transcriptomic profiling in combination with a proteomics analysis to elucidate processes that occur during this growth-to-dormancy transition in a conifer, white spruce (Picea glauca[Moench] Voss). Several differentially expressed genes were likely associated with the developmental transition that occurs during growth cessation in the cambial zone and the concomitant completion of cell maturation in vascular tissues. Genes encoding for cell wall and membrane biosynthetic enzymes showed transcript abundance patterns consistent with completion of cell maturation, and also of cell wall and membrane modifications potentially enabling cells to withstand the harsh conditions of winter. Several differentially expressed genes were identified that encoded putative regulators of cambial activity, cell development and of the photoperiodic pathway. Reconfiguration of carbon allocation figured centrally in the tree's overwintering preparations. For example, genes associated with carbon-based defences such as terpenoids were down-regulated, while many genes associated with protein-based defences and other stress mitigation mechanisms were up-regulated. Several of these correspond to proteins that were accumulated during the growth-to-dormancy transition, emphasizing the importance of stress protection in the tree's adaptive response to overwintering.

Population structure and migration pattern of a conifer pathogen, Grosmannia clavigera, as influenced by its symbiont, the mountain pine beetle

We investigated the population structure of Grosmannia clavigera (Gc), a fungal symbiont of the mountain pine beetle (MPB) that plays a crucial role in the establishment and reproductive success of this pathogen. This insect-fungal complex has destroyed over 16 million ha of lodgepole pine forests in Canada, the largest MPB epidemic in recorded history. During this current epidemic, MPB has expanded its range beyond historically recorded boundaries, both northward and eastward, and has now reached the jack pine of Alberta, potentially threatening the Canadian boreal forest. To better understand the dynamics between the beetle and its fungal symbiont, we sampled 19 populations in western North America and genotyped individuals from these populations with eight microsatellite markers. The fungus displayed high haplotype diversity, with over 250 unique haplotypes observed in 335 single spore isolates. Linkage equilibria in 13 of the 19 populations suggested that the fungus reproduces sexually. Bayesian clustering and distance analyses identified four genetic clusters that corresponded to four major geographical regions, which suggested that the epidemic arose from multiple geographical sources. A genetic cluster north of the Rocky Mountains, where the MPB has recently become established, experienced a population bottleneck, probably as a result of the recent range expansion. The two genetic clusters located north and west of the Rocky Mountains contained many fungal isolates admixed from all populations, possibly due to the massive movement of MPB during the epidemic. The general agreement in north-south differentiation of MPB and G. clavigera populations points to the fungal pathogen's dependence on the movement of its insect vector. In addition, the patterns of diversity and the individual assignment tests of the fungal associate suggest that migration across the Rocky Mountains occurred via a northeastern corridor, in accordance with meteorological patterns and observation of MPB movement data. Our results highlight the potential of this pathogen for both expansion and sexual reproduction, and also identify some possible barriers to gene flow. Understanding the ecological and evolutionary dynamics of this fungus-beetle association is important for the modelling and prediction of MPB epidemics.

A white spruce gene catalog for conifer genome analyses

Several angiosperm plant genomes, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), poplar (Populus trichocarpa), and grapevine (Vitis vinifera), have been sequenced, but the lack of reference genomes in gymnosperm phyla reduces our understanding of plant evolution and restricts the potential impacts of genomics research. A gene catalog was developed for the conifer tree Picea glauca (white spruce) through large-scale expressed sequence tag sequencing and full-length cDNA sequencing to facilitate genome characterizations, comparative genomics, and gene mapping. The resource incorporates new and publicly available sequences into 27,720 cDNA clusters, 23,589 of which are represented by full-length insert cDNAs. Expressed sequence tags, mate-pair cDNA clone analysis, and custom sequencing were integrated through an iterative process to improve the accuracy of clustering outcomes. The entire catalog spans 30 Mb of unique transcribed sequence. We estimated that the P. glauca nuclear genome contains up to 32,520 transcribed genes owing to incomplete, partially sequenced, and unsampled transcripts and that its transcriptome could span up to 47 Mb. These estimates are in the same range as the Arabidopsis and rice transcriptomes. Next-generation methods confirmed and enhanced the catalog by providing deeper coverage for rare transcripts, by extending many incomplete clusters, and by augmenting the overall transcriptome coverage to 38 Mb of unique sequence. Genomic sample sequencing at 8.5% of the 19.8-Gb P. glauca genome identified 1,495 clusters representing highly repeated sequences among the cDNA clusters. With a conifer transcriptome in full view, functional and protein domain annotations clearly highlighted the divergences between conifers and angiosperms, likely reflecting their respective evolutionary paths.

Effect of water stress and fungal inoculation on monoterpene emission from an historical and a new pine host of the mountain pine beetle

The mountain pine beetle (Dendroctonus ponderosae, MPB) has killed millions of lodgepole pine (Pinus contorta) trees in Western Canada, and recent range expansion has resulted in attack of jack pine (Pinus banksiana) in Alberta. Establishment of MPB in the Boreal forest will require use of jack pine under a suite of environmental conditions different from those it typically encounters in its native range. Lodgepole and jack pine seedlings were grown under controlled environment conditions and subjected to either water deficit or well watered conditions and inoculated with Grosmannia clavigera, a MPB fungal associate. Soil water content, photosynthesis, stomatal conductance, and emission of volatile organic compounds (VOCs) were monitored over the duration of the six-week study. Monoterpene content of bark and needle tissue was measured at the end of the experiment. β-Phellandrene, the major monoterpene in lodgepole pine, was almost completely lacking in the volatile emission profile of jack pine. The major compound in jack pine was α-pinene. The emission of both compounds was positively correlated with stomatal conductance. 3-Carene was emitted at a high concentration from jack pine seedlings, which is in contrast to monoterpene profiles of jack pine from more southern and eastern parts of its range. Fungal inoculation caused a significant increase in total monoterpene emission in water deficit lodgepole pine seedlings right after its application. By 4 weeks into the experiment, water deficit seedlings of both species released significantly lower levels of total monoterpenes than well watered seedlings. Needle tissue contained lower total monoterpene content than bark. Generally, monoterpene tissue content increased over time independent from any treatment. The results suggest that monoterpenes that play a role in pine-MPB interactions differ between lodgepole and jack pine, and also that they are affected by water availability.

Spatial community structure of mountain pine beetle fungal symbionts across a latitudinal gradient

Symbiont redundancy in obligate insect-fungal systems is thought to buffer the insect host against symbiont loss and to extend the environmental conditions under which the insect can persist. The mountain pine beetle is associated with at least three well-known and putatively obligate ophiostomatoid fungal symbionts that vary in their environmental tolerances. To better understand the spatial variation in beetle-fungal symbiotic associations, we examined the community composition of ophiostomatoid fungi associated with the mountain pine beetle as a function of latitude and elevation. The region investigated represents the leading edge of a recent outbreak of mountain pine beetle in western Canada. Using regression and principal components analysis, we identified significant spatial patterns in fungal species abundances that indicate symmetrical replacement between two of the three fungi along a latitudinal gradient and little variation in response to elevation. We also identified significant variation in the prevalence of pair-wise species combinations that occur within beetle galleries. Frequencies of pair-wise combinations were significantly different from what was expected given overall species abundances. These results suggest that complex processes of competitive exclusion and coexistence help determine fungal community composition and that the consequences of these processes vary spatially. The presence of three fungal symbionts in different proportions and combinations across a wide range of environmental conditions may help explain the success of mountain pine beetle attacks across a broad geographic range.

Mountain pine beetle host-range expansion threatens the boreal forest

The current epidemic of the mountain pine beetle (MPB), an indigenous pest of western North American pine, has resulted in significant losses of lodgepole pine. The leading edge has reached Alberta where forest composition shifts from lodgepole to jack pine through a hybrid zone. The susceptibility of jack pine to MPB is a major concern, but there has been no evidence of host-range expansion, in part due to the difficulty in distinguishing the parentals and their hybrids. We tested the utility of a panel of microsatellite loci optimized for both species to classify lodgepole pine, jack pine and their hybrids using simulated data. We were able to accurately classify simulated individuals, and hence applied these markers to identify the ancestry of attacked trees. Here we show for the first time successful MPB attack in natural jack pine stands at the leading edge of the epidemic. This once unsuitable habitat is now a novel environment for MPB to exploit, a potential risk which could be exacerbated by further climate change. The consequences of host-range expansion for the vast boreal ecosystem could be significant.

Molecular events of apical bud formation in white spruce, Picea glauca

Bud formation is an adaptive trait that temperate forest trees have acquired to facilitate seasonal synchronization. We have characterized transcriptome-level changes that occur during bud formation of white spruce [Picea glauca (Moench) Voss], a primarily determinate species in which preformed stem units contained within the apical bud constitute most of next season's growth. Microarray analysis identified 4460 differentially expressed sequences in shoot tips during short day-induced bud formation. Cluster analysis revealed distinct temporal patterns of expression, and functional classification of genes in these clusters implied molecular processes that coincide with anatomical changes occurring in the developing bud. Comparing expression profiles in developing buds under long day and short day conditions identified possible photoperiod-responsive genes that may not be essential for bud development. Several genes putatively associated with hormone signalling were identified, and hormone quantification revealed distinct profiles for abscisic acid (ABA), cytokinins, auxin and their metabolites that can be related to morphological changes to the bud. Comparison of gene expression profiles during bud formation in different tissues revealed 108 genes that are differentially expressed only in developing buds and show greater transcript abundance in developing buds than other tissues. These findings provide a temporal roadmap of bud formation in white spruce.

Gene family structure, expression and functional analysis of HD-Zip III genes in angiosperm and gymnosperm forest trees

BACKGROUND

Class III Homeodomain Leucine Zipper (HD-Zip III) proteins have been implicated in the regulation of cambium identity, as well as primary and secondary vascular differentiation and patterning in herbaceous plants. They have been proposed to regulate wood formation but relatively little evidence is available to validate such a role. We characterised and compared HD-Zip III gene family in an angiosperm tree, Populus spp. (poplar), and the gymnosperm Picea glauca (white spruce), representing two highly evolutionarily divergent groups.

RESULTS

Full-length cDNA sequences were isolated from poplar and white spruce. Phylogenetic reconstruction indicated that some of the gymnosperm sequences were derived from lineages that diverged earlier than angiosperm sequences, and seem to have been lost in angiosperm lineages. Transcript accumulation profiles were assessed by RT-qPCR on tissue panels from both species and in poplar trees in response to an inhibitor of polar auxin transport. The overall transcript profiles HD-Zip III complexes in white spruce and poplar exhibited substantial differences, reflecting their evolutionary history. Furthermore, two poplar sequences homologous to HD-Zip III genes involved in xylem development in Arabidopsis and Zinnia were over-expressed in poplar plants. PtaHB1 over-expression produced noticeable effects on petiole and primary shoot fibre development, suggesting that PtaHB1 is involved in primary xylem development. We also obtained evidence indicating that expression of PtaHB1 affected the transcriptome by altering the accumulation of 48 distinct transcripts, many of which are predicted to be involved in growth and cell wall synthesis. Most of them were down-regulated, as was the case for several of the poplar HD-Zip III sequences. No visible physiological effect of over-expression was observed on PtaHB7 transgenic trees, suggesting that PtaHB1 and PtaHB7 likely have distinct roles in tree development, which is in agreement with the functions that have been assigned to close homologs in herbaceous plants.

CONCLUSIONS

This study provides an overview of HD-zip III genes related to woody plant development and identifies sequences putatively involved in secondary vascular growth in angiosperms and in gymnosperms. These gene sequences are candidate regulators of wood formation and could be a source of molecular markers for tree breeding related to wood properties.

Functional characterization of drought-responsive aquaporins in Populus balsamifera and Populus simonii×balsamifera clones with different drought resistance strategies

We have characterized poplar aquaporins (AQPs) to investigate their possible functions in differential drought responses of Populus balsamifera and Populus simonii×balsamifera leaves. Plants were exposed to mild and severe levels of drought stress and to drought stress recovery treatment, and their responses were compared with well-watered controls. Compared with P. balsamifera, P. simonii×balsamifera used drought avoidance as the main drought resistance strategy, and rapidly reduced stomatal conductance in response to stress. This strategy is correlated with growth rate reductions. Eleven AQPs were transcriptionally profiled in leaves from these experiments and five were functionally characterized for water channel activity. PIP1;3 and PIP2;5 were among the most highly expressed leaf AQPs that were responsive to drought. Expression of PIP1;3 and five other AQPs increased in response to drought in the leaves of P. simonii×balsamifera but not in P. balsamifera, suggesting a possible role of these AQPs in water redistribution in the leaf tissues. PIP2;5 was upregulated in P. balsamifera, but not in P. simonii×balsamifera, suggesting that this AQP supports the transpiration-driven water flow. Functional characterization of five drought-responsive plasma membrane intrinsic proteins (PIPs) demonstrated that three PIP2 AQPs (PIP2;2, PIP2;5, PIP2;7) functioned as water transporters in Xenopus laevis oocytes, while the two PIP1 AQPs (PIP1;2 and PIP1;3) did not, consistent with the notion that they may be functional only as heterotetramers.

High nitrogen fertilization and stem leaning have overlapping effects on wood formation in poplar but invoke largely distinct molecular pathways

Previous studies indicated that high nitrogen fertilization may impact secondary xylem development and alter fibre anatomy and composition. The resulting wood shares some resemblance with tension wood, which has much thicker cell walls than normal wood due to the deposition of an additional layer known as the G-layer. This report compares the short-term effects of high nitrogen fertilization and tree leaning to induce tension wood, either alone or in combination, upon wood formation in young trees of Populus trichocarpa (Torr. & Gray) × P. deltoides Bartr. ex Marsh. Fibre anatomy, chemical composition and transcript profiles were examined in newly formed secondary xylem. Each of the treatments resulted in thicker cell walls relative to the controls. High nitrogen and tree leaning had overlapping effects on chemical composition based on Fourier transform infrared analysis, specifically indicating that secondary cell wall composition was shifted in favour of cellulose and hemicelluloses relative to lignin content. In contrast, the high-nitrogen trees had shorter fibres, whilst the leaning trees had longer fibres that the controls. Microarray transcript profiling carried out after 28 days of treatment identified 180 transcripts that accumulated differentially in one or more treatments. Only 10% of differentially expressed transcripts were affected in all treatments relative to the controls. Several of the affected transcripts were related to carbohydrate metabolism, secondary cell wall formation, nitrogen metabolism and osmotic stress. RT-qPCR analyses at 1, 7 and 28 days showed that several transcripts followed very different accumulation profiles in terms of rate and level of accumulation, depending on the treatment. Our findings suggest that high nitrogen fertilization and tension wood induction elicit largely distinct and molecular pathways with partial overlap. When combined, the two types of environmental cue yielded additive effects.

Influence of nitrogen fertilization on xylem traits and aquaporin expression in stems of hybrid poplar

We studied the influence of nitrogen (N) on hydraulic traits and aquaporin (AQP) expression in the stem xylem of hybrid poplar saplings (Populus trichocarpa (Torr. & Gray) x deltoides Bartr. ex Marsh clone H11-11). Plants were grown in a controlled environment and were kept well watered throughout the experiments. Hydraulic measurements were done on basal and distal stem segments of plants receiving high N fertilization (high N plants) versus plants receiving only adequate N fertilization (adequate N plants). High N plants grew faster and exhibited more leaf area than adequate N controls. These morphological differences were paralleled by wider vessels and higher specific conductivities (K(S)) in high N plants. However, stems of high N plants were more vulnerable to xylem cavitation, at least in one of two experiments, and showed lower wood densities than stems of adequate N plants. Leaf area was strongly correlated with cross-sectional xylem area in both plant groups. Since higher K(S) in high N plants was accompanied by concomitant increases in leaf area, leaf-specific conductivities were similar in both plant groups. Influences of N on hydraulic traits were paralleled by changes in AQP expression. Seven AQPs were upregulated in the stem xylem of high N plants, five of which have been identified recently as water transporters. The enhanced growth of secondary xylem of high N plants has been shown to result from both increased cambial activity as well as increased cell size. We suggest that some of these water-transporting AQPs could play a role in xylogenesis, facilitating the influx of water into the zone of differentiating and maturing cells in secondary xylem, including expanding vessels.

Aquaporin gene expression and apoplastic water flow in bur oak (Quercus macrocarpa) leaves in relation to the light response of leaf hydraulic conductance

It has previously been shown that hydraulic conductance in bur oak leaves (Quercus macrocarpa Michx.), measured with the high pressure flow meter technique (HPFM), can significantly increase within 30 min following exposure to high irradiance. The present study investigated whether this increase could be explained by an increase in the cell-to-cell pathway and whether the response is linked to changes in the transcript level corresponding to aquaporin genes. Four cDNA sequences showing high similarity to members of the aquaporin gene family from other plant species were characterized from bur oak leaves and the expression levels of these cDNA sequences were examined in leaves by quantitative real-time PCR (QRT-PCR). No change was found in the relative transcript abundance corresponding to these four putative aquaporin genes in leaves with light-induced high hydraulic conductance (exposed to high irradiance) compared to leaves with low hydraulic conductance (exposed to low irradiance). However, in sun leaves that were exposed to different light levels prior to leaf collection (full sunlight, shade, and covered with aluminium foil for 16 h), the relative transcript levels of two of the putative aquaporin genes increased several-fold in shaded leaves compared to the sun-exposed or covered leaves. When the leaves were pressure-infiltrated with the apoplastic tracer dye trisodium 3-hydroxy-5,8,10-pyrenetrisulphonate (PTS(3), 0.02%), there was no change in the PTS(3) concentration of leaf exudates collected in ambient light or in high irradiance, but there was a small apoplastic acidification. There was also no change in PTS(3) concentration between the leaves infiltrated under high irradiance with 0.02% PTS(3) or with 0.1 mM HgCl(2) in 0.02% PTS(3). The results suggest that the putative aquaporin genes that were identified in the present study probably do not play a role in the light responses of hydraulic conductance at the transcript level, but they may function in regulating water homeostasis in leaves adapted to different light conditions. In addition, it is shown that high irradiance induced changes in the pH of the apoplast and that there does not appear to be a significant shift to the cell-to-cell mediated water transport in bur oak leaves exposed to high irradiance as measured by the apoplastic tracer dye.

Involvement of Pinus taeda MYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis

The involvement of two R2R3-MYB genes from Pinus taeda L., PtMYB1 and PtMYB8, in phenylpropanoid metabolism and secondary cell wall biogenesis was investigated in planta. These pine MYBs were constitutively overexpressed (OE) in Picea glauca (Moench) Voss, used as a heterologous conifer expression system. Morphological, histological, chemical (lignin and soluble phenols), and transcriptional analyses, i.e. microarray and reverse transcription quantitative PCR (RT-qPCR) were used for extensive phenotyping of MYB-overexpressing spruce plantlets. Upon germination of somatic embryos, root growth was reduced in both transgenics. Enhanced lignin deposition was also a common feature but ectopic secondary cell wall deposition was more strongly associated with PtMYB8-OE. Microarray and RT-qPCR data showed that overexpression of each MYB led to an overlapping up-regulation of many genes encoding phenylpropanoid enzymes involved in lignin monomer synthesis, while misregulation of several cell wall-related genes and other MYB transcription factors was specifically associated with PtMYB8-OE. Together, the results suggest that MYB1 and MYB8 may be part of a conserved transcriptional network involved in secondary cell wall deposition in conifers.

Potato, Solanum tuberosum, defense against Colorado potato beetle, Leptinotarsa decemlineata (Say): microarray gene expression profiling of potato by Colorado potato beetle regurgitant treatment of wounded leaves

Colorado potato beetle (CPB) is a leading pest of solanaceous plants. Despite the economic importance of this pest, surprisingly few studies have been carried out to characterize its molecular interaction with the potato plant. In particular, little is known about the effect of CPB elicitors on gene expression associated with the plant's defense response. In order to discover putative CPB elicitor-responsive genes, the TIGR 11,421 EST Solanaceae microarray was used to identify genes that are differentially expressed in response to the addition of CPB regurgitant to wounded potato leaves. By applying a cutoff corresponding to an adjusted P-value of <0.01 and a fold change of >1.5 or <0.67, we found that 73 of these genes are induced by regurgitant treatment of wounded leaves when compared to wounding alone, whereas 54 genes are repressed by this treatment. This gene set likely includes regurgitant-responsive genes as well as wounding-responsive genes whose expression patterns are further enhanced by the presence of regurgitant. Real-time polymerase chain reaction was used to validate differential expression by regurgitant treatment for five of these genes. In general, genes that encoded proteins involved in secondary metabolism and stress were induced by regurgitant; genes associated with photosynthesis were repressed. One induced gene that encodes aromatic amino acid decarboxylase is responsible for synthesis of the precursor of 2-phenylethanol. This is significant because 2-phenylethanol is recognized by the CPB predator Perillus bioculatis. In addition, three of the 16 type 1 and type 2 proteinase inhibitor clones present on the potato microarray were repressed by application of CPB regurgitant to wounded leaves. Given that proteinase inhibitors are known to interfere with digestion of proteins in the insect midgut, repression of these proteinase inhibitors by CPB may inhibit this component of the plant's defense arsenal. These data suggest that beyond the wound response, CPB elicitors play a role in mediating the plant/insect interaction.

Effects of increased nitrogen supply on the lignification of poplar wood

The short-term influence of adequate and high nitrogen fertilization on poplar lignification was investigated. The high nitrogen supply decreased lignin staining in the newly formed secondary xylem, indicating that lignin deposition was affected. Acetyl bromide determinations gave a 9-10% decrease in lignin content; however, Klason lignin content was unchanged. Thioacidolysis showed that elevated N supply affected lignin structure such that there was a reduced frequency of lignin units involved in beta-O-4 bonds, a reduced syringyl/guaiacyl ratio, an increased frequency of p-hydroxyphenyl lignin units, more guaiacyl units with free phenolic groups, and more p-hydroxybenzoic acid ester-linked to poplar lignins. These features suggest that lignins from poplars grown under high N bear structural similarities to lignins formed during early stages of wood development. The findings also indicate that a gravitational stimulus inducing the formation of tension wood and high N availability lead to similar and additive effects on lignin content and structure.

Short-term physiological and developmental responses to nitrogen availability in hybrid poplar

Nitrogen fertilization induces dramatic changes in the growth and development of plants, including forest trees. In this study we examined short-term responses of hybrid poplar, Populus balsamifera ssp. trichocarpa x deltoides, to N fertilization. Glasshouse-grown saplings subjected to limiting, intermediate, and luxuriant levels of ammonium nitrate over a 28 d time course demonstrated rapid changes to whole-plant architecture and biomass accumulation. Nitrogen-associated shifts in allocation occurred in temporally distinct stages. Nitrogen availability modulated parameters that affect carbon gain, including light-saturated net photosynthesis and leaf area. These parameters were affected by N-induced changes to leaf maturation and senescence. Leaf area was also affected by N-induced sylleptic branch development. Genes encoding vegetative storage proteins and starch biosynthetic enzymes exhibited contrasting patterns of expression under differential N availability. A gene encoding a previously uncharacterized putative pectin methylesterase inhibitor displayed expression patterns comparable to the starch biosynthetic genes. The results of this study illustrate the phenotypic plasticity that P. balsamifera ssp. trichocarpa x deltoides exhibits in response to differential N availability.

Nitrogen storage and seasonal nitrogen cycling in Populus: bridging molecular physiology and ecophysiology

While both annual and perennial plants store nitrogen resources during the growing season, seasonal N cycling is a hallmark of the perennial habit. In Populus the vegetative storage proteins BSP, WIN4 and PNI288 all play a role in N storage during active growth, whereas BSP is the major form of reduced N storage during winter dormancy. In this review we explore cellular and molecular events implicated in seasonal N cycling in Populus, as well as environmental cues that modulate both the phenology of seasonal N cycling, and the efficiency and proficiency of autumn N resorption. We highlight recent advances that have been made using Populus genomics resources to address processes germane to seasonal N cycling. Genetic and genetological studies are enabling us to connect our understanding of seasonal N cycling at molecular and cellular levels with that at ecophysiological levels. With the genomics resources and foundational knowledge that are now in place, Populus researchers are poised to build an integrative understanding of seasonal N cycling that spans from genomes to ecosystems.

A robust neural networks approach for spatial and intensity-dependent normalization of cDNA microarray data

MOTIVATION

Microarray experiments are affected by numerous sources of non-biological variation that contribute systematic bias to the resulting data. In a dual-label (two-color) cDNA or long-oligonucleotide microarray, these systematic biases are often manifested as an imbalance of measured fluorescent intensities corresponding to Sample A versus those corresponding to Sample B. Systematic biases also affect between-slide comparisons. Making effective corrections for these systematic biases is a requisite for detecting the underlying biological variation between samples. Effective data normalization is therefore an essential step in the confident identification of biologically relevant differences in gene expression profiles. Several normalization methods for the correction of systemic bias have been described. While many of these methods have addressed intensity-dependent bias, few have addressed both intensity-dependent and spatiality-dependent bias.

RESULTS

We present a neural network-based normalization method for correcting the intensity- and spatiality-dependent bias in cDNA microarray datasets. In this normalization method, the dependence of the log-intensity ratio (M) on the average log-intensity (A) as well as on the spatial coordinates (X,Y) of spots is approximated with a feed-forward neural network function. Resistance to outliers is provided by assigning weights to each spot based on how distant their M values is from the median over the spots whose A values are similar, as well as by using pseudospatial coordinates instead of spot row and column indices. A comparison of the robust neural network method with other published methods demonstrates its potential in reducing both intensity-dependent bias and spatial-dependent bias, which translates to more reliable identification of truly regulated genes.

Pathogen challenge, salicylic acid, and jasmonic acid regulate expression of chitinase gene homologs in pine

To better understand the molecular regulation of defense responses in members of the genus Pinus, we tested the expression of various chitinase homologs in response to pathogen-associated signals. PSCHI4, a putative extracellular class II chitinase, was secreted into liquid medium by pine cells and was also secreted by transgenic tobacco cells that ectopically expressed pschi4. Extracellular proteins of pine were separated by isoelectric focusing; PSCHI4 was not associated with fractions containing detectable beta-N-acetylglucosaminidase or lysozyme activities. However, other fractions contained enzyme activities that increased markedly after elicitor treatment. The pschi4 transcript and protein accumulated in pine seedlings challenged with the necrotrophic pathogen Fusarium subglutinans f. sp. pini, with the protein reaching detectable levels in susceptible seedlings concomitant with the onset of visible disease symptoms. Additional chitinase transcripts, assigned to classes I and IV based on primary sequence analysis, were also induced by pathogen challenge. Jasmonic acid induced class I and class IV but not class II chitinase, whereas salicylic acid induced all three classes of chitinase. These results show that multiple chitinase homologs are induced after challenge by a necrotrophic pathogen and by potential signaling molecules identified in angiosperms. This suggests the potential importance of de novo pathogenesis-related (PR) gene expression in pathogen defense responses of pine trees.