The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change

Species are the unit of analysis in many global change and conservation biology studies; however, species are not uniform entities but are composed of different, sometimes locally adapted, populations differing in plasticity. We examined how intraspecific variation in thermal niches and phenotypic plasticity will affect species distributions in a warming climate. We first developed a conceptual model linking plasticity and niche breadth, providing five alternative intraspecific scenarios that are consistent with existing literature. Secondly, we used ecological niche-modeling techniques to quantify the impact of each intraspecific scenario on the distribution of a virtual species across a geographically realistic setting. Finally, we performed an analogous modeling exercise using real data on the climatic niches of different tree provenances. We show that when population differentiation is accounted for and dispersal is restricted, forecasts of species range shifts under climate change are even more pessimistic than those using the conventional assumption of homogeneously high plasticity across a species' range. Suitable population-level data are not available for most species so identifying general patterns of population differentiation could fill this gap. However, the literature review revealed contrasting patterns among species, urging greater levels of integration among empirical, modeling and theoretical research on intraspecific phenotypic variation.

Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes

BACKGROUND

Molecular evolutionary studies share the common goal of elucidating historical relationships, and the common challenge of adequately sampling taxa and characters. Particularly at low taxonomic levels, recent divergence, rapid radiations, and conservative genome evolution yield limited sequence variation, and dense taxon sampling is often desirable. Recent advances in massively parallel sequencing make it possible to rapidly obtain large amounts of sequence data, and multiplexing makes extensive sampling of megabase sequences feasible. Is it possible to efficiently apply massively parallel sequencing to increase phylogenetic resolution at low taxonomic levels?

RESULTS

We reconstruct the infrageneric phylogeny of Pinus from 37 nearly-complete chloroplast genomes (average 109 kilobases each of an approximately 120 kilobase genome) generated using multiplexed massively parallel sequencing. 30/33 ingroup nodes resolved with > or = 95% bootstrap support; this is a substantial improvement relative to prior studies, and shows massively parallel sequencing-based strategies can produce sufficient high quality sequence to reach support levels originally proposed for the phylogenetic bootstrap. Resampling simulations show that at least the entire plastome is necessary to fully resolve Pinus, particularly in rapidly radiating clades. Meta-analysis of 99 published infrageneric phylogenies shows that whole plastome analysis should provide similar gains across a range of plant genera. A disproportionate amount of phylogenetic information resides in two loci (ycf1, ycf2), highlighting their unusual evolutionary properties.

CONCLUSION

Plastome sequencing is now an efficient option for increasing phylogenetic resolution at lower taxonomic levels in plant phylogenetic and population genetic analyses. With continuing improvements in sequencing capacity, the strategies herein should revolutionize efforts requiring dense taxon and character sampling, such as phylogeographic analyses and species-level DNA barcoding.

Characterisation of a pine MYB that regulates lignification

A member of the R2R3-MYB family of transcription factors was cloned from a cDNA library constructed from RNA isolated from differentiating pine xylem. This MYB, Pinus taeda MYB4 (PtMYB4), is expressed in cells undergoing lignification, as revealed by in situ RT-PCR. Electrophoretic mobility shift assays (EMSAs) showed that recombinant PtMYB4 protein is able to bind to DNA motifs known as AC elements. AC elements are ubiquitous in the promoters of genes encoding lignin biosynthetic enzymes. Transcriptional activation assays using yeast showed that PtMYB4 could activate transcription in an AC-element-dependent fashion. Overexpression of PtMYB4 in transgenic tobacco plants altered the accumulation of transcripts corresponding to genes encoding lignin biosynthetic enzymes. Lignin deposition increased in transgenic tobacco plants that overexpressed PtMYB4, and extended to cell types that do not normally lignify. Taken together, these findings are consistent with the hypothesis that PtMYB4 is sufficient to induce lignification, and that it may play this role during wood formation in pine.

TIR-X and TIR-NBS proteins: two new families related to disease resistance TIR-NBS-LRR proteins encoded in Arabidopsis and other plant genomes

The Toll/interleukin-1 receptor (TIR) domain is found in one of the two large families of homologues of plant disease resistance proteins (R proteins) in Arabidopsis and other dicotyledonous plants. In addition to these TIR-NBS-LRR (TNL) R proteins, we identified two families of TIR-containing proteins encoded in the Arabidopsis Col-0 genome. The TIR-X (TX) family of proteins lacks both the nucleotide-binding site (NBS) and the leucine rich repeats (LRRs) that are characteristic of the R proteins, while the TIR-NBS (TN) proteins contain much of the NBS, but lack the LRR. In Col-0, the TX family is encoded by 27 genes and three pseudogenes; the TN family is encoded by 20 genes and one pseudogene. Using massively parallel signature sequencing (MPSS), expression was detected at low levels for approximately 85% of the TN-encoding genes. Expression was detected for only approximately 40% of the TX-encoding genes, again at low levels. Physical map data and phylogenetic analysis indicated that multiple genomic duplication events have increased the numbers of TX and TN genes in Arabidopsis. Genes encoding TX, TN and TNL proteins were demonstrated in conifers; TX and TN genes are present in very low numbers in grass genomes. The expression, prevalence, and diversity of TX and TN genes suggests that these genes encode functional proteins rather than resulting from degradation or deletions of TNL genes. These TX and TN proteins could be plant analogues of small TIR-adapter proteins that function in mammalian innate immune responses such as MyD88 and Mal.

Transcriptional regulation of lignin biosynthesis

Lignin is the second most abundant plant biopolymer mainly present in the secondary walls of tracheary elements and fibers in wood. Understanding how lignin is biosynthesized has long been an interest to plant biologists and will have a significant impact on tree biotechnology. Lignin is polymerized from monolignols that are synthesized through the lignin biosynthetic pathway. To make lignin, all the genes in the lignin biosynthetic pathway need to be coordinately turned on. It has been shown that a common cis-element, namely the AC element, is present in the majority of the lignin biosynthetic genes and required for their expression in lignifying cells. Important progress has been made in the identification of transcription factors that bind to the AC elements and are potentially involved in the coordinated regulation of lignin biosynthesis. The Arabidopsis MYB58 and MYB63 as well as their poplar ortholog PtrMYB28 are transcriptional activators of the lignin biosynthetic pathway, whereas the eucalyptus EgMYB2 and pine PtMYB4 transcription factors are likely Arabidopsis MYB46 orthologs involved in the regulation of the entire secondary wall biosynthetic program. It was found that the transcriptional regulation of lignin biosynthesis is under the control of the same transcriptional network regulating the biosynthesis of other secondary wall components, including cellulose and xylan. The identification of transcription factors directly activating lignin biosynthetic genes provides unprecedented tools to potentially manipulate the amount of lignin in wood and other plant products based on our needs.

Evolution of genome size in pines (Pinus) and its life-history correlates: supertree analyses

Genome size has been suggested to be a fundamental biological attribute in determining life-history traits in many groups of organisms. We examined the relationships between pine genome sizes and pine phylogeny, environmental factors (latitude, elevation, annual rainfall), and biological traits (latitudinal and elevational ranges, seed mass, minimum generation time, interval between large seed crops, seed dispersal mode, relative growth rate, measures of potential and actual invasiveness, and level of rarity). Genome sizes were determined for 60 pine taxa and then combined with published values to make a dataset encompassing 85 species, or 70% of species in the genus. Supertrees were constructed using 20 published source phylogenies. Ancestral genome size was estimated as 32 pg. Genome size has apparently remained stable or increased over evolutionary time in subgenus Strobus, while it has decreased in most subsections in subgenus Pinus. We analyzed relationships between genome size and life-history variables using cross-species correlations and phylogenetically independent contrasts derived from supertree constructions. The generally assumed positive relation between genome size and minimum generation time could not be confirmed in phylogenetically controlled analyses. We found that the strongest correlation was between genome size and seed mass. Because the growth quantities specific leaf area and leaf area ratio (and to a lesser extent relative growth rate) are strongly negatively related to seed mass, they were also negatively correlated with genome size. Northern latitudinal limit was negatively correlated with genome size. Invasiveness, particularly of wind-dispersed species, was negatively associated with both genome size and seed mass. Seed mass and its relationships with seed number, dispersal mode, and growth rate contribute greatly to the differences in life-history strategies of pines. Many life-history patterns are therefore indirectly, but consistently, associated with genome size.

Apparent homology of expressed genes from wood-forming tissues of loblolly pine (Pinus taeda L.) with Arabidopsis thaliana

Pinus taeda L. (loblolly pine) and Arabidopsis thaliana differ greatly in form, ecological niche, evolutionary history, and genome size. Arabidopsis is a small, herbaceous, annual dicotyledon, whereas pines are large, long-lived, coniferous forest trees. Such diverse plants might be expected to differ in a large number of functional genes. We have obtained and analyzed 59,797 expressed sequence tags (ESTs) from wood-forming tissues of loblolly pine and compared them to the gene sequences inferred from the complete sequence of the Arabidopsis genome. Approximately 50% of pine ESTs have no apparent homologs in Arabidopsis or any other angiosperm in public databases. When evaluated by using contigs containing long, high-quality sequences, we find a higher level of apparent homology between the inferred genes of these two species. For those contigs 1,100 bp or longer, approximately 90% have an apparent Arabidopsis homolog (E value < 10-10). Pines and Arabidopsis last shared a common ancestor approximately 300 million years ago. Few genes would be expected to retain high sequence similarity for this time if they did not have essential functions. These observations suggest substantial conservation of gene sequence in seed plants.

Plant expansins are a complex multigene family with an ancient evolutionary origin

Expansins are a group of extracellular proteins that directly modify the mechanical properties of plant cell walls, leading to turgor-driven cell extension. Within the completely sequenced Arabidopsis genome, we identified 38 expansin sequences that fall into three discrete subfamilies. Based on phylogenetic analysis and shared intron patterns, we propose a new, systematic nomenclature of Arabidopsis expansins. Further phylogenetic analysis, including expansin sequences found here in monocots, pine (Pinus radiata, Pinus taeda), fern (Regnellidium diphyllum, Marsilea quadrifolia), and moss (Physcomitrella patens) indicate that the three plant expansin subfamilies arose and began diversifying very early in, if not before, colonization of land by plants. Closely related "expansin-like" sequences were also identified in the social amoeba, Dictyostelium discoidium, suggesting that these wall-modifying proteins have a very deep evolutionary origin.

Phylogeny, historical biogeography, and patterns of diversification for Pinus (Pinaceae): Phylogenetic tests of fossil-based hypotheses

Pines comprise one of the largest coniferous genera, are distributed throughout the Northern Hemisphere, and have an abundant fossil record. Distributions of fossils have been used to derive a three-step hypothesis of early pine evolution, which postulates a Mesozoic origin for the genus, east-west expansions across Laurasia, and retraction into Eocene refugia. Here, we present phylogenetic tests of this hypothesis using chloroplast sequence data from four loci for 83 pine species. We used the fossil-based hypothesis to derive null expectations concerning monophyly of taxonomic groups, dates of cladogenesis, and patterns of diversification. Phylogenetic analyses using several algorithms subsequently provided rigorous tests of these expectations. Our inferred phylogenies illustrated broad congruence with taxonomic groups, but highlighted consistent problems within subgenus Strobus. Estimated minimum dates of divergence derived from relaxed clock methods were largely consistent with the fossil record and yielded a date for the ingroup node of Pinus of 128+/-4 mya, depending upon the calibration used for subgenus Pinus. Ancestral area reconstructions showed Pinus to have most likely originated in Eurasia. Major clades differed in biogeographic patterns, but were consistent with the fossil-based hypothesis. We found weak support, however, for a change in diversification rate in the Eocene as interpretations of fossil distributions would have predicted.

Phylogeography of maritime pine inferred with organelle markers having contrasted inheritance

Range-wide variation of maritime pine was studied at maternally inherited and paternally inherited markers (mitochondrial DNA and chloroplast DNA). While chloroplast DNA exhibits the highest diversity, phylogeographic inferences from this marker are blurred by homoplasy and extensive pollen flow. In contrast, the only three mitochondrial haplotypes found provide a clear picture of nonoverlapping areas colonized from different refugia, with no single population having a mixed composition (GST = 1). Comparison of the genetic structure inferred from both organelle genomes allows the investigation of differential seed and pollen dispersal, pointing to pollen, but not seed, dispersal across the Strait of Gibraltar (from Morocco into Iberia). A comparison with already available genetic information, especially that of one of the maritime pine's most threatening insect pests, the bast scale Matsucoccus feytaudi, further completes the picture.

Fossil Calibration of Molecular Divergence Infers a Moderate Mutation Rate and Recent Radiations for Pinus

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.

An autoregulatory feedback loop involving PAP1 and TAS4 in response to sugars in Arabidopsis

miR828 in Arabidopsis triggers the cleavage of Trans-Acting SiRNA Gene 4 (TAS4) transcripts and production of small interfering RNAs (ta-siRNAs). One siRNA, TAS4-siRNA81(-), targets a set of MYB transcription factors including PAP1, PAP2, and MYB113 which regulate the anthocyanin biosynthesis pathway. Interestingly, miR828 also targets MYB113, suggesting a close relationship between these MYBs, miR828, and TAS4, but their evolutionary origins are unknown. We found that PAP1, PAP2, and TAS4 expression is induced specifically by exogenous treatment with sucrose and glucose in seedlings. The induction is attenuated in abscisic acid (ABA) pathway mutants, especially in abi3-1 and abi5-1 for PAP1 or PAP2, while no such effect is observed for TAS4. PAP1 is under regulation by TAS4, demonstrated by the accumulation of PAP1 transcripts and anthocyanin in ta-siRNA biogenesis pathway mutants. TAS4-siR81(-) expression is induced by physiological concentrations of Suc and Glc and in pap1-D, an activation-tagged line, indicating a feedback regulatory loop exists between PAP1 and TAS4. Bioinformatic analysis revealed MIR828 homologues in dicots and gymnosperms, but only in one basal monocot, whereas TAS4 is only found in dicots. Consistent with this observation, PAP1, PAP2, and MYB113 dicot paralogs show peptide and nucleotide footprints for the TAS4-siR81(-) binding site, providing evidence for purifying selection in contrast to monocots. Extended sequence similarities between MIR828, MYBs, and TAS4 support an inverted duplication model for the evolution of MIR828 from an ancestral gymnosperm MYB gene and subsequent formation of TAS4 by duplication of the miR828* arm. We obtained evidence by modified 5'-RACE for a MYB mRNA cleavage product guided by miR828 in Pinus resinosa. Taken together, our results suggest that regulation of anthocyanin biosynthesis by TAS4 and miR828 in higher plants is evolutionarily significant and consistent with the evolution of TAS4 since the dicot-monocot divergence.

Genome hypermethylation in Pinus silvestris of Chernobyl--a mechanism for radiation adaptation?

Adaptation is a complex process by which populations of organisms respond to long-term environmental stresses by permanent genetic change. Here we present data from the natural "open-field" radiation adaptation experiment after the Chernobyl accident and provide the first evidence of the involvement of epigenetic changes in adaptation of a eukaryote-Scots pine (Pinus silvestris), to chronic radiation exposure. We have evaluated global genome methylation of control and radiation-exposed pine trees using a method based on cleavage by a methylation-sensitive HpaII restriction endonuclease that leaves a 5' guanine overhang and subsequent single nucleotide extension with labeled [3H] dCTP. We have found that genomic DNA of exposed pine trees was considerably hypermethylated. Moreover, hypermethylation appeared to be dependent upon the radiation dose absorbed by the trees. Such hypermethylation may be viewed as a defense strategy of plants that prevents genome instability and reshuffling of the hereditary material, allowing survival in an extreme environment. Further studies are clearly needed to analyze in detail the involvement of DNA methylation and other epigenetic mechanisms in the complex process of radiation stress and adaptive response.

Accelerating the domestication of trees using genomic selection: accuracy of prediction models across ages and environments

• Genomic selection is increasingly considered vital to accelerate genetic improvement. However, it is unknown how accurate genomic selection prediction models remain when used across environments and ages. This knowledge is critical for breeders to apply this strategy in genetic improvement. • Here, we evaluated the utility of genomic selection in a Pinus taeda population of c. 800 individuals clonally replicated and grown on four sites, and genotyped for 4825 single-nucleotide polymorphism (SNP) markers. Prediction models were estimated for diameter and height at multiple ages using genomic random regression best linear unbiased predictor (BLUP). • Accuracies of prediction models ranged from 0.65 to 0.75 for diameter, and 0.63 to 0.74 for height. The selection efficiency per unit time was estimated as 53-112% higher using genomic selection compared with phenotypic selection, assuming a reduction of 50% in the breeding cycle. Accuracies remained high across environments as long as they were used within the same breeding zone. However, models generated at early ages did not perform well to predict phenotypes at age 6 yr. • These results demonstrate the feasibility and remarkable gain that can be achieved by incorporating genomic selection in breeding programs, as long as models are used at the relevant selection age and within the breeding zone in which they were estimated.

Evolution of diterpene metabolism: Sitka spruce CYP720B4 catalyzes multiple oxidations in resin acid biosynthesis of conifer defense against insects

Diterpene resin acids (DRAs) are specialized (secondary) metabolites of the oleoresin defense of conifers produced by diterpene synthases and cytochrome P450s of the CYP720B family. The evolution of DRA metabolism shares common origins with the biosynthesis of ent-kaurenoic acid, which is highly conserved in general (primary) metabolism of gibberellin biosynthesis. Transcriptome mining in species of spruce (Picea) and pine (Pinus) revealed CYP720Bs of four distinct clades. We cloned a comprehensive set of 12 different Sitka spruce (Picea sitchensis) CYP720Bs as full-length cDNAs. Spatial expression profiles, methyl jasmonate induction, and transcript enrichment in terpenoid-producing resin ducts suggested a role of CYP720B4 in DRA biosynthesis. CYP720B4 was characterized as a multisubstrate, multifunctional enzyme by the formation of oxygenated diterpenoids in metabolically engineered yeast, yeast in vivo transformation of diterpene substrates, in vitro assays with CYP720B4 protein produced in Escherichia coli, and alteration of DRA profiles in RNA interference-suppressed spruce seedlings. CYP720B4 was active with 24 different diterpenoid substrates, catalyzing consecutive C-18 oxidations in the biosynthesis of an array of diterpene alcohols, aldehydes, and acids. CYP720B4 was most active in the formation of dehydroabietic acid, a compound associated with insect resistance of Sitka spruce. We identified patterns of convergent evolution of CYP720B4 in DRA metabolism and ent-kaurene oxidase CYP701 in gibberellin metabolism and revealed differences in the evolution of specialized and general diterpene metabolism in a gymnosperm. The genomic and functional characterization of the gymnosperm CYP720B family highlights that the evolution of specialized metabolism involves substantial diversification relative to conserved, general metabolism.

Reciprocal selection causes a coevolutionary arms race between crossbills and lodgepole pine

Few studies have shown both reciprocal selection and reciprocal adaptations for a coevolving system in the wild. The goal of our study was to determine whether the patterns of selection on Rocky Mountain lodgepole pine (Pinus contorta spp. latifolia) and red crossbills (Loxia curvirostra complex) were concordant with earlier published evidence of reciprocal adaptations in lodgepole pine and crossbills on isolated mountain ranges in the absence of red squirrels (Tamiasciurus hudsonicus). We found that selection (directional) by crossbills on lodgepole pine where Tamiasciurus are absent was divergent from the selection (directional) exerted by Tamiasciurus on lodgepole pine. This resulted in divergent selection between areas with and without Tamiasciurus that was congruent with the geographic patterns of cone variation. In the South Hills, Idaho, where Tamiasciurus are absent and red crossbills are thought to be coevolving with lodgepole pine, crossbills experienced stabilizing selection on bill size, with cone structure as the agent of selection. These results show that crossbills and lodgepole pine exhibit reciprocal adaptations in response to reciprocal selection, and they provide insight into the traits mediating and responding to selection in a coevolutionary arms race.

The evolution of chloroplast RNA editing

RNA editing alters the nucleotide sequence of an RNA molecule so that it deviates from the sequence of its DNA template. Different RNA-editing systems are found in the major eukaryotic lineages, and these systems are thought to have evolved independently. In this study, we provide a detailed analysis of data on C-to-U editing sites in land plant chloroplasts and propose a model for the evolution of RNA editing in land plants. First, our data suggest that the limited RNA-editing system of seed plants and the much more extensive systems found in hornworts and ferns are of monophyletic origin. Further, although some eukaryotic editing systems appear to have evolved to regulate gene expression, or at least are now involved in gene regulation, there is no evidence that RNA editing plays a role in gene regulation in land plant chloroplasts. Instead, our results suggest that land plant chloroplast C-to-U RNA editing originated as a mechanism to generate variation at the RNA level, which could complement variation at the DNA level. Under this model, many of the original sites, particularly in seed plants, have been subsequently lost due to mutation at the DNA level, and the function of extant sites is merely to conserve certain codons. This is the first comprehensive model for the evolution of the chloroplast RNA-editing system of land plants and may also be applicable to the evolution of RNA editing in plant mitochondria.

Suppression of 4-coumarate-CoA ligase in the coniferous gymnosperm Pinus radiata

Severe suppression of 4-coumarate-coenzyme A ligase (4CL) in the coniferous gymnosperm Pinus radiata substantially affected plant phenotype and resulted in dwarfed plants with a "bonsai tree-like" appearance. Microscopic analyses of stem sections from 2-year-old plants revealed substantial morphological changes in both wood and bark tissues. This included the formation of weakly lignified tracheids that displayed signs of collapse and the development of circumferential bands of axial parenchyma. Acetyl bromide-soluble lignin assays and proton nuclear magnetic resonance studies revealed lignin reductions of 36% to 50% in the most severely affected transgenic plants. Two-dimensional nuclear magnetic resonance and pyrolysis-gas chromatography-mass spectrometry studies indicated that lignin reductions were mainly due to depletion of guaiacyl but not p-hydroxyphenyl lignin. 4CL silencing also caused modifications in the lignin interunit linkage distribution, including elevated beta-aryl ether (beta-O-4 unit) and spirodienone (beta-1) levels, accompanied by lower phenylcoumaran (beta-5), resinol (beta-beta), and dibenzodioxocin (5-5/beta-O-4) levels. A sharp depletion in the level of saturated (dihydroconiferyl alcohol) end groups was also observed. Severe suppression of 4CL also affected carbohydrate metabolism. Most obvious was an up to approximately 2-fold increase in galactose content in wood from transgenic plants due to increased compression wood formation. The molecular, anatomical, and analytical data verified that the isolated 4CL clone is associated with lignin biosynthesis and illustrated that 4CL silencing leads to complex, often surprising, physiological and morphological changes in P. radiata.

Range-wide phylogeography and gene zones in Pinus pinaster Ait. revealed by chloroplast microsatellite markers

Some 1339 trees from 48 Pinus pinaster stands were characterized by five chloroplast microsatellites, detecting a total of 103 distinct haplotypes. Frequencies for the 16 most abundant haplotypes (p(k) > 0.01) were spatially interpolated over a lattice made by 430 grid points. Fitting of spatially interpolated values on raw haplotype frequencies at the same geographical location was tested by regression analysis. A range-wide 'diversity map' based on interpolated haplotype frequencies allowed the identification of one 'hotspot' of diversity in central and southeastern Spain, and two areas of low haplotypic diversity located in the western Iberian peninsula and Morocco. Principal component analysis (PCA) carried out on haplotypes frequency surfaces allowed the construction of a colour-based 'synthetic' map of the first three PC components, enabling the detection of the main range-scale genetic trends and the identification of three main 'gene pools' for the species: (i) a 'southeastern' gene pool, including southeastern France, Italy, Corsica, Sardinia, Pantelleria and northern Africa; (ii) an 'Atlantic' gene pool, including all the western areas of the Iberian peninsula; and (iii) a 'central' gene pool, located in southeastern Spain. Multivariate and AMOVA analyses carried out on interpolated grid point frequency values revealed the existence of eight major clusters ('gene zones'), whose genetic relationships were related with the history of the species. In addition, demographic models showed more ancient expansions in the eastern and southern ranges of maritime pine probably associated to early postglacial recolonization. The delineation of the gene zones provides a baseline for designing conservation areas in this key Mediterranean pine.

Fire-adapted traits of Pinus arose in the fiery Cretaceous

• The mapping of functional traits onto chronograms is an emerging approach for the identification of how agents of natural selection have shaped the evolution of organisms. Recent research has reported fire-dependent traits appearing among flowering plants from 60 million yr ago (Ma). Although there are many records of fossil charcoal in the Cretaceous (65-145 Ma), evidence of fire-dependent traits evolving in that period is lacking. • We link the evolutionary trajectories for five fire-adapted traits in Pinaceae with paleoatmospheric conditions over the last 250 million yr to determine the time at which fire originated as a selective force in trait evolution among seed plants. • Fire-protective thick bark originated in Pinus c. 126 Ma in association with low-intensity surface fires. More intense crown fires emerged c. 89 Ma coincident with thicker bark and branch shedding, or serotiny with branch retention as an alternative strategy. These innovations appeared at the same time as the Earth's paleoatmosphere experienced elevated oxygen levels that led to high burn probabilities during the mid-Cretaceous. • The fiery environments of the Cretaceous strongly influenced trait evolution in Pinus. Our evidence for a strong correlation between the evolution of fire-response strategies and changes in fire regime 90-125 Ma greatly backdates the key role that fire has played in the evolution of seed plants.

DNA from bird-dispersed seed and wind-disseminated pollen provides insights into postglacial colonization and population genetic structure of whitebark pine (Pinus albicaulis)

Uniparentally inherited mitochondrial (mt)DNA and chloroplast (cp)DNA microsatellites (cpSSRs) were used to examine population genetic structure and biogeographic patterns of bird-dispersed seed and wind-disseminated pollen of whitebark pine (Pinus albicaulis Engelm.). Sampling was conducted from 41 populations throughout the range of the species. Analyses provide evidence for an ancestral haplotype and two derived mtDNA haplotypes with distinct regional distributions. An abrupt contact zone between mtDNA haplotypes in the Cascade Range suggests postglacial biogeographic movements. Among three cpSSR loci, 42 haplotypes were detected within 28 cpSSR sample populations that were aggregated into six regions. Analysis of molecular variance (amova) was used to determine the hierarchical genetic structure of cpSSRs. amova and population pairwise comparisons (FST ) of cpSSR, and geographical distribution of mtDNA haplotypes provide insights into historical changes in biogeography. The genetic data suggest that whitebark pine has been intimately tied to climatic change and associated glaciation, which has led to range movements facilitated by seed dispersal by Clark's nutcracker (Nucifraga columbiana Wilson). The two hypotheses proposed to explain the genetic structure are: (i) a northward expansion into Canada and the northern Cascades in the early Holocene; and (ii) historical gene flow between Idaho and the Oregon Cascades when more continuous habitat existed in Central Oregon during the late Pleistocene. Genetic structure and insights gained from historical seed movements provide a basis on which to develop recovery plans for a species that is at risk from multiple threats.

Evolution of genome size and complexity in Pinus

BACKGROUND

Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood.

METHODOLOGY/PRINCIPAL FINDINGS

Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA.

CONCLUSIONS/SIGNIFICANCE

Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.

Characterisation of Pt MYB1, an R2R3-MYB from pine xylem

A cDNA encoding a member of the R2R3-MYB family of transcription factors was cloned from a library constructed from differentiating Pinus taeda (loblolly pine) xylem RNA. This MYB family member, Pinus taeda MYB1 (PtMYB1), was most abundantly expressed in differentiating xylem, as assessed by both ribonuclease protection assays, and by northern blot analysis with poly(A)-enriched RNA. Similar to other plant R2R3-MYB family members, recombinant Pt MYB1 protein was able to bind to AC elements in electrophoretic mobility shift assays (EMSAs). AC elements are DNA motifs rich in adenosine and cytosine that have been implicated in the xylem-localised regulation of genes encoding lignin biosynthetic enzymes. Pt MYB1 not only bound to AC elements, but was also able to induce AC-element-dependent shifts in the electrophoretic mobility of a plant promoter that contains three AC elements, the minimal PHENYLALANINE AMMONIA-LYASE 2 (PAL2) promoter from Phaseolus vulgaris. Transcriptional activation assays conducted using yeast showed that Pt MYB1 also activated transcription, and that it did so in an AC-element-dependent fashion. Pt MYB1 also activated transcription from the minimal PAL2 promoter in plant cells in an AC-element-dependent fashion, as revealed by transient transcriptional activation assays with microprojectile-bombarded tobacco NT-1 cells. Taken together, these finding are consistent with the hypothesis that Pt MYB1 may regulate transcription from cis -acting AC elements in pine xylem.

Limited genetic variability and phenotypic plasticity detected for cavitation resistance in a Mediterranean pine

Resistance to cavitation is a major determinant of plant survival under severe drought and can be used to quantify species adaptive potential. Interspecific variation in this key trait is well defined in woody species, but intraspecific variation (level and structure) resulting from standing genetic variation and phenotypic plasticity has never been determined. Combining for the first time in situ characterization of natural populations and two reciprocal common gardens in dry and wet sites, we estimated variance components (phenotypic, genetic, environmental, and genetic × environmental) of cavitation resistance based on 513 genotypes of a Mediterranean pine, Pinus pinaster. Despite the selected populations being climatically contrasted, phenotypic plasticity in resistance to cavitation remained low and was essentially attributed to family level. Between-population variation in cavitation resistance for both phenotypic and genetic variation was limited. These results strongly suggest that cavitation resistance is buffered against genetic and to a lesser extent environmental variation (canalization) in maritime pine. Consequently, in a drier world, the increasing drought tolerance of Pinus species might be severely constrained by the low level of cavitation resistance variation, resulting in a large-scale loss of productivity.

Exploring lignification in conifers by silencing hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyltransferase in Pinus radiata

The enzyme hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyltransferase (HCT) is involved in the production of methoxylated monolignols that are precursors to guaiacyl and syringyl lignin in angiosperm species. We identified and cloned a putative HCT gene from Pinus radiata, a coniferous gymnosperm that does not produce syringyl lignin. This gene was up-regulated during tracheary element (TE) formation in P. radiata cell cultures and showed 72.6% identity to the amino acid sequence of the Nicotiana tabacum HCT isolated earlier. RNAi-mediated silencing of the putative HCT gene had a strong impact on lignin content, monolignol composition, and interunit linkage distribution. AcBr assays revealed an up to 42% reduction in lignin content in TEs. Pyrolysis-GC/MS, thioacidolysis, and NMR detected substantial changes in lignin composition. Most notable was the rise of p-hydroxyphenyl units released by thioacidolysis, which increased from trace amounts in WT controls to up to 31% in transgenics. Two-dimensional 13C-1H correlative NMR confirmed the increase in p-hydroxyphenyl units in the transgenics and revealed structural differences, including an increase in resinols, a reduction in dibenzodioxocins, and the presence of glycerol end groups. The observed modifications in silenced transgenics validate the targeted gene as being associated with lignin biosynthesis in P. radiata and thus likely to encode HCT. This enzyme therefore represents the metabolic entry point leading to the biosynthesis of methoxylated phenylpropanoids in angiosperm species and coniferous gymnosperms such as P. radiata.