The Norway spruce genome sequence and conifer genome evolution.

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding.

Partial correlation analysis of transcriptomes helps detangle the growth and defense network in spruce.

In plants, there can be a trade-off between resource allocations to growth vs defense. Here, we use partial correlation analysis of gene expression to make inferences about the nature of this interaction. We studied segregating progenies of Interior spruce subject to weevil attack. In a controlled experiment, we measured pre-attack plant growth and post-attack damage with several morphological measures, and profiled transcriptomes of 188 progeny. We used partial correlations of individual transcripts (expressed sequence tags, ESTs) with pairs of growth/defense traits to identify important nodes and edges in the inferred underlying gene network, for example, those pairs of growth/defense traits with high mutual correlation with a single EST transcript. We give a method to identify such ESTs. A terpenoid ABC transporter gene showed strongest correlations (P = 0.019); its transcript represented a hub within the compact 166-member gene-gene interaction network (P = 0.004) of the negative genetic correlations between growth and subsequent pest attack. A small 21-member interaction network (P = 0.004) represented the uncovered positive correlations. Our study demonstrates partial correlation analysis identifies important gene networks underlying growth and susceptibility to the weevil in spruce. In particular, we found transcripts that strongly modify the trade-off between growth and defense, and allow identification of networks more central to the trade-off.

Insights into conifer giga-genomes.

Insights from sequenced genomes of major land plant lineages have advanced research in almost every aspect of plant biology. Until recently, however, assembled genome sequences of gymnosperms have been missing from this picture. Conifers of the pine family (Pinaceae) are a group of gymnosperms that dominate large parts of the world's forests. Despite their ecological and economic importance, conifers seemed long out of reach for complete genome sequencing, due in part to their enormous genome size (20-30 Gb) and the highly repetitive nature of their genomes. Technological advances in genome sequencing and assembly enabled the recent publication of three conifer genomes: white spruce (Picea glauca), Norway spruce (Picea abies), and loblolly pine (Pinus taeda). These genome sequences revealed distinctive features compared with other plant genomes and may represent a window into the past of seed plant genomes. This Update highlights recent advances, remaining challenges, and opportunities in light of the publication of the first conifer and gymnosperm genomes.

Contrasting patterns of genetic diversity across the ranges of Pinus monticola and P. strobus: a comparison between eastern and western North American postglacial colonization histories.

UNLABELLED: • Premises of the study: Understanding the influence of recent glacial and postglacial periods on species' distributions is key for predicting the effects of future environmental changes. We investigated the influence of two physiographic landscapes on population structure and postglacial colonization of two white pine species of contrasting habitats: P. monticola, which occurs in the highly mountainous region of western North America, and P. strobus, which occurs in a much less mountainous area in eastern North America.• METHODS: To characterize the patterns of genetic diversity and population structure across the ranges of both species, 158 and 153 single nucleotide polymorphism (SNP) markers derived from expressed genes were genotyped on range-wide samples of 61 P. monticola and 133 P. strobus populations, respectively.• KEY RESULTS: In P. monticola, a steep latitudinal decrease in genetic diversity likely resulted from postglacial colonization involving rare long-distance dispersal (LDD) events. In contrast, no geographic patterns of diversity were detected in P. strobus, suggesting recolonization via a gradually advancing front or frequent LDD events. For each species, structure analyses identified two distinct southern and northern genetic groups that likely originated from two different glacial lineages. At a finer scale, and for the two species, smaller subgroups were detected that could be remnants of cryptic refugia.• CONCLUSION: During postglacial colonization, the western and eastern North American landscapes had different impacts on genetic signatures in P. monticola compared with P. strobus. We discuss the importance of our findings for conservation programs and predictions of species' response to climate change.

Reduced Genetic Diversity and Increased Dispersal in Guigna (Leopardus guigna) in Chilean Fragmented Landscapes.

Landscape fragmentation is often a major cause of species extinction as it can affect a wide variety of ecological processes. The impact of fragmentation varies among species depending on many factors, including their life-history traits and dispersal abilities. Felids are one of the groups most threatened by fragmented landscapes because of their large home ranges, territorial behavior, and low population densities. Here, we model the impacts of habitat fragmentation on patterns of genetic diversity in the guigna (Leopardus guigna), a small felid that is closely associated with the heavily human-impacted temperate rainforests of southern South America. We assessed genetic variation in 1798 base pairs of mitochondrial DNA sequences, 15 microsatellite loci, and 2 sex chromosome genes and estimated genetic diversity, kinship, inbreeding, and dispersal in 38 individuals from landscapes with differing degrees of fragmentation on Chiloé Island in southern Chile. Increased fragmentation was associated with reduced genetic diversity, but not with increased kinship or inbreeding. However, in fragmented landscapes, there was a weaker negative correlation between pairwise kinship and geographic distance, suggesting increased dispersal distances. These results highlight the importance of biological corridors to maximize connectivity in fragmented landscapes and contribute to our understanding of the broader genetic consequences of habitat fragmentation, especially for forest-specialist carnivores.

Evolution of gene structure in the conifer Picea glauca: a comparative analysis of the impact of intron size.

BACKGROUND: A positive relationship between genome size and intron length is observed across eukaryotes including Angiosperms plants, indicating a co-evolution of genome size and gene structure. Conifers have very large genomes and longer introns on average than most plants, but impacts of their large genome and longer introns on gene structure has not be described. RESULTS: Gene structure was analyzed for 35 genes of Picea glauca obtained from BAC sequencing and genome assembly, including comparisons with A. thaliana, P. trichocarpa and Z. mays. We aimed to develop an understanding of impact of long introns on the structure of individual genes. The number and length of exons was well conserved among the species compared but on average, P. glauca introns were longer and genes had four times more intronic sequence than Arabidopsis, and 2 times more than poplar and maize. However, pairwise comparisons of individual genes gave variable results and not all contrasts were statistically significant. Genes generally accumulated one or a few longer introns in species with larger genomes but the position of long introns was variable between plant lineages. In P. glauca, highly expressed genes generally had more intronic sequence than tissue preferential genes. Comparisons with the Pinus taeda BACs and genome scaffolds showed a high conservation for position of long introns and for sequence of short introns. A survey of 1836 P. glauca genes obtained by sequence capture mostly containing introns <1 Kbp showed that repeated sequences were 10× more abundant in introns than in exons. CONCLUSION: Conifers have large amounts of intronic sequence per gene for seed plants due to the presence of few long introns and repetitive element sequences are ubiquitous in their introns. Results indicate a complex landscape of intron sizes and distribution across taxa and between genes with different expression profiles.

Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data.

UNLABELLED: White spruce (Picea glauca) is a dominant conifer of the boreal forests of North America, and providing genomics resources for this commercially valuable tree will help improve forest management and conservation efforts. Sequencing and assembling the large and highly repetitive spruce genome though pushes the boundaries of the current technology. Here, we describe a whole-genome shotgun sequencing strategy using two Illumina sequencing platforms and an assembly approach using the ABySS software. We report a 20.8 giga base pairs draft genome in 4.9 million scaffolds, with a scaffold N50 of 20,356 bp. We demonstrate how recent improvements in the sequencing technology, especially increasing read lengths and paired end reads from longer fragments have a major impact on the assembly contiguity. We also note that scalable bioinformatics tools are instrumental in providing rapid draft assemblies. AVAILABILITY: The Picea glauca genome sequencing and assembly data are available through NCBI (Accession#: ALWZ0100000000 PID: PRJNA83435). http://www.ncbi.nlm.nih.gov/bioproject/83435.

Relatedness coefficients and their applications for triplets and quartets of genetic markers.

Relatedness coefficients which seek the identity-by-descent of genetic markers are described. The markers are in groups of two, three or four, and if four, can consist of two pairs. It is essential to use cumulants (not moments) for four-marker-gene probabilities, as the covariance of homozygosity, used in four-marker applications, can only be described with cumulants. A covariance of homozygosity between pairs of markers arises when populations follow a mixture distribution. Also, the probability of four markers all identical-by-descent equals the normalized fourth cumulant. In this article, a "genetic marker" generally represents either a gene locus or an allele at a locus. Applications of three marker coefficients mainly involve conditional regression, and applications of four marker coefficients can involve identity disequilibrium. Estimation of relatedness using genetic marker data is discussed. However, three- and four-marker estimators suffer from statistical and numerical problems, including higher statistical variance, complexity of estimation formula, and singularity at some intermediate allele frequencies.

Aging in a long-lived clonal tree.

From bacteria to multicellular animals, most organisms exhibit declines in survivorship or reproductive performance with increasing age ("senescence"). Evidence for senescence in clonal plants, however, is scant. During asexual growth, we expect that somatic mutations, which negatively impact sexual fitness, should accumulate and contribute to senescence, especially among long-lived clonal plants. We tested whether older clones of Populus tremuloides (trembling aspen) from natural stands in British Columbia exhibited significantly reduced reproductive performance. Coupling molecular-based estimates of clone age with male fertility data, we observed a significant decline in the average number of viable pollen grains per catkin per ramet with increasing clone age in trembling aspen. We found that mutations reduced relative male fertility in clonal aspen populations by about 5.8 x 10(-5) to 1.6 x 10(-3) per year, leading to an 8% reduction in the number of viable pollen grains, on average, among the clones studied. The probability that an aspen lineage ultimately goes extinct rises as its male sexual fitness declines, suggesting that even long-lived clonal organisms are vulnerable to senescence.

Variation and Natural Heritability of Blueberry Floral Volatiles.

Vaccinium corymbosum (highbush blueberry) in British Columbia relies on honeybee pollination for a reliable fruit set. As floral volatiles may help explain pollinator preference for blueberry, we surveyed components of variation for volatiles using gas chromatography-mass spectrometry (GC/MS). Principal component analysis of GC chromatogram peaks revealed a grouping of cultivars by a biosynthetic pathway that also corresponds with their known pedigree. To identify genetic variance, we identified 34 chemicals with adequate sample sizes. We estimated "natural heritability" (using uncontrolled crosses in natural environments) in two ways: (1) as "clonal repeatability," which is equivalent to broad-sense heritability and is an upper bound for narrow-sense heritability, and (2) "marker-based heritability," which serves as a lower bound for narrow-sense heritability. Both methods indicate that heritability is relatively low, ca. 15%, and variable among traits. This is expected as the floral volatile release is changeable and dependent on environmental conditions. It might be possible to use highly heritable volatiles for breeding.

Slow but not low: genomic comparisons reveal slower evolutionary rate and higher dN/dS in conifers compared to angiosperms.

BACKGROUND: Comparative genomics can inform us about the processes of mutation and selection across diverse taxa. Among seed plants, gymnosperms have been lacking in genomic comparisons. Recent EST and full-length cDNA collections for two conifers, Sitka spruce (Picea sitchensis) and loblolly pine (Pinus taeda), together with full genome sequences for two angiosperms, Arabidopsis thaliana and poplar (Populus trichocarpa), offer an opportunity to infer the evolutionary processes underlying thousands of orthologous protein-coding genes in gymnosperms compared with an angiosperm orthologue set. RESULTS: Based upon pairwise comparisons of 3,723 spruce and pine orthologues, we found an average synonymous genetic distance (dS) of 0.191, and an average dN/dS ratio of 0.314. Using a fossil-established divergence time of 140 million years between spruce and pine, we extrapolated a nucleotide substitution rate of 0.68 × 10(-9) synonymous substitutions per site per year. When compared to angiosperms, this indicates a dramatically slower rate of nucleotide substitution rates in conifers: on average 15-fold. Coincidentally, we found a three-fold higher dN/dS for the spruce-pine lineage compared to the poplar-Arabidopsis lineage. This joint occurrence of a slower evolutionary rate in conifers with higher dN/dS, and possibly positive selection, showcases the uniqueness of conifer genome evolution. CONCLUSIONS: Our results are in line with documented reduced nucleotide diversity, conservative genome evolution and low rates of diversification in conifers on the one hand and numerous examples of local adaptation in conifers on the other hand. We propose that reduced levels of nucleotide mutation in large and long-lived conifer trees, coupled with large effective population size, were the main factors leading to slow substitution rates but retention of beneficial mutations.

Long-distance gene flow and adaptation of forest trees to rapid climate change.

Forest trees are the dominant species in many parts of the world and predicting how they might respond to climate change is a vital global concern. Trees are capable of long-distance gene flow, which can promote adaptive evolution in novel environments by increasing genetic variation for fitness. It is unclear, however, if this can compensate for maladaptive effects of gene flow and for the long-generation times of trees. We critically review data on the extent of long-distance gene flow and summarise theory that allows us to predict evolutionary responses of trees to climate change. Estimates of long-distance gene flow based both on direct observations and on genetic methods provide evidence that genes can move over spatial scales larger than habitat shifts predicted under climate change within one generation. Both theoretical and empirical data suggest that the positive effects of gene flow on adaptation may dominate in many instances. The balance of positive to negative consequences of gene flow may, however, differ for leading edge, core and rear sections of forest distributions. We propose future experimental and theoretical research that would better integrate dispersal biology with evolutionary quantitative genetics and improve predictions of tree responses to climate change.

Estimating selfing rates from reconstructed pedigrees using multilocus genotype data.

Several methods have been developed to estimate the selfing rate of a population from a sample of individuals genotyped for several marker loci. These methods can be based on homozygosity excess (or inbreeding), identity disequilibrium, progeny array (PA) segregation or population assignment incorporating partial selfing. Progeny array-based method is generally the best because it is not subject to some assumptions made by other methods (such as lack of misgenotyping, absence of biparental inbreeding and presence of inbreeding equilibrium), and it can reveal other facets of a mixed-mating system such as patterns of shared paternity. However, in practice, it is often difficult to obtain PAs, especially for animal species. In this study, we propose a method to reconstruct the pedigree of a sample of individuals taken from a monoecious diploid population practicing mixed mating, using multilocus genotypic data. Selfing and outcrossing events are then detected when an individual derives from identical parents and from two distinct parents, respectively. Selfing rate is estimated by the proportion of selfed offspring in the reconstructed pedigree of a sample of individuals. The method enjoys many advantages of the PA method, but without the need of a priori family structure, although such information, if available, can be utilized to improve the inference. Furthermore, the new method accommodates genotyping errors, estimates allele frequencies jointly and is robust to the presence of biparental inbreeding and inbreeding disequilibrium. Both simulated and empirical data were analysed by the new and previous methods to compare their statistical properties and accuracies.

Defense mechanisms against herbivory in Picea: sequence evolution and expression regulation of gene family members in the phenylpropanoid pathway.

BACKGROUND: In trees, a substantial amount of carbon is directed towards production of phenolics for development and defense. This metabolic pathway is also a major factor in resistance to insect pathogens in spruce. In such gene families, environmental stimuli may have an important effect on the evolutionary fate of duplicated genes, and different expression patterns may indicate functional diversification. RESULTS: Gene families in spruce (Picea) have expanded to superfamilies, including O-methyltransferases, cytochrome-P450, and dirigents/classIII-peroxidases. Neo-functionalization of superfamily members from different clades is reflected in expression diversification. Genetical genomics can provide new insights into the genetic basis and evolution of insect resistance in plants. Adopting this approach, we merged genotype data (252 SNPs in a segregating pedigree), gene expression levels (for 428 phenylpropanoid-related genes) and measures of susceptibility to Pissodes stobi, using a partial-diallel crossing-design with white spruce (Picea glauca). Thirty-eight expressed phenylpropanoid-related genes co-segregated with weevil susceptibility, indicating either causative or reactive effects of these genes to weevil resistance. We identified eight regulatory genomic regions with extensive overlap of quantitative trait loci from susceptibility and growth phenotypes (pQTLs) and expression QTL (eQTL) hotspots. In particular, SNPs within two different CCoAOMT loci regulate phenotypic variation from a common set of 24 genes and three resistance traits. CONCLUSIONS: Pest resistance was associated with individual candidate genes as well as with trans-regulatory hotspots along the spruce genome. Our results showed that specific genes within the phenylpropanoid pathway have been duplicated and diversified in the conifer in a process fundamentally different from short-lived angiosperm species. These findings add to the information about the role of the phenylpropanoid pathway in the evolution of plant defense mechanisms against insect pests and provide substantial potential for the functional characterization of several not yet resolved alternative pathways in plant defenses.

QTL mapping in white spruce: gene maps and genomic regions underlying adaptive traits across pedigrees, years and environments.

BACKGROUND: The genomic architecture of bud phenology and height growth remains poorly known in most forest trees. In non model species, QTL studies have shown limited application because most often QTL data could not be validated from one experiment to another. The aim of our study was to overcome this limitation by basing QTL detection on the construction of genetic maps highly-enriched in gene markers, and by assessing QTLs across pedigrees, years, and environments. RESULTS: Four saturated individual linkage maps representing two unrelated mapping populations of 260 and 500 clonally replicated progeny were assembled from 471 to 570 markers, including from 283 to 451 gene SNPs obtained using a multiplexed genotyping assay. Thence, a composite linkage map was assembled with 836 gene markers.For individual linkage maps, a total of 33 distinct quantitative trait loci (QTLs) were observed for bud flush, 52 for bud set, and 52 for height growth. For the composite map, the corresponding numbers of QTL clusters were 11, 13, and 10. About 20% of QTLs were replicated between the two mapping populations and nearly 50% revealed spatial and/or temporal stability. Three to four occurrences of overlapping QTLs between characters were noted, indicating regions with potential pleiotropic effects. Moreover, some of the genes involved in the QTLs were also underlined by recent genome scans or expression profile studies.Overall, the proportion of phenotypic variance explained by each QTL ranged from 3.0 to 16.4% for bud flush, from 2.7 to 22.2% for bud set, and from 2.5 to 10.5% for height growth. Up to 70% of the total character variance could be accounted for by QTLs for bud flush or bud set, and up to 59% for height growth. CONCLUSIONS: This study provides a basic understanding of the genomic architecture related to bud flush, bud set, and height growth in a conifer species, and a useful indicator to compare with Angiosperms. It will serve as a basic reference to functional and association genetic studies of adaptation and growth in Picea taxa. The putative QTNs identified will be tested for associations in natural populations, with potential applications in molecular breeding and gene conservation programs. QTLs mapping consistently across years and environments could also be the most important targets for breeding, because they represent genomic regions that may be least affected by G × E interactions.

Evolutionary potential in the wild: more than meets the eye.

The genus Aquilegia consists of 60-70 perennial plant species widely distributed throughout the northern hemisphere. Its flowers have a delicate and ornamental appearance that makes them a favourite of gardeners. In this genus, adaptive radiations for both floral and vegetative traits have occurred. These adaptive radiations, and the key phylogenetic placement of Aquilegia between Arabidopsis and rice, make this genus a 'model system' for plant evolution (Kramer 2009). In this issue, Castellanos et al. (2011) use a marker-based method to infer heritability for floral and vegetative traits in two Aquilegia species. Layered on top of this are estimates of the strength of natural selection. This novel joint estimation of heritability and selection in the wild showed that vegetative traits, compared to floral traits, have the highest evolutionarily potential. Evolutionary potential is the most important quantity to measure in wild populations. It combines inheritance and strength of selection and predicts the potential for populations to adapt to changing environments. The combination of molecular techniques with species in natural environments makes this work a model for molecular ecological investigations.

Performing parentage analysis for polysomic inheritances based on allelic phenotypes.

Polyploidy poses several problems for parentage analysis. We present a new polysomic inheritance model for parentage analysis based on genotypes or allelic phenotypes to solve these problems. The effects of five factors are simultaneously accommodated in this model: (1) double-reduction, (2) null alleles, (3) negative amplification, (4) genotyping errors and (5) self-fertilization. To solve genotyping ambiguity (unknown allele dosage), we developed a new method to establish the likelihood formulas for allelic phenotype data and to simultaneously include the effects of our five chosen factors. We then evaluated and compared the performance of our new method with three established methods by using both simulated data and empirical data from the cultivated blueberry (Vaccinium corymbosum). We also developed and compared the performance of two additional estimators to estimate the genotyping error rate and the sample rate. We make our new methods freely available in the software package polygene, at http://github.com/huangkang1987/polygene.

The correlation of evolutionary rate with pathway position in plant terpenoid biosynthesis.

Genes are expected to face stronger selective constraint and to evolve more slowly if they encode enzymes upstream as opposed to downstream in metabolic pathways, because upstream genes are more pleiotropic, being required for a wider range of end products. However, few clear examples of this trend in evolutionary rate variation exist. We examined whether genes involved in plant terpenoid biosynthesis exhibit such a pattern, using data for 40 genes from five fully sequenced angiosperms, Oryza, Vitis, Arabidopsis, Populus, and Ricinus. Our results show that d(N)/d(S) does in fact correlate with pathway position along pathways converting glucose to the terpenoid phytohormones abscissic acid, gibberellic acid (GA), and brassinosteroids. Upstream versus downstream rate variation is particularly strong in the GA pathway. In contrast, we found no or little apparent variation in d(N)/d(S) with gene copy number. We also introduce a new measure of pathway position, the Pathway Pleiotropy Index (PPI), which counts groups of enzymes between pathway branch points. We found that this measure is superior to pathway position in explaining variation in d(N)/d(S) along each pathway. Although at least 8 of the 40 genes showed evidence of positive selection, correlations of d(N)/d(S) with PPI remain significant when these genes are removed. Therefore, our results are consistent with the prediction that selective constraint is progressively relaxed along metabolic pathways.

Widespread, ecologically relevant genetic markers developed from association mapping of climate-related traits in Sitka spruce (Picea sitchensis).

• Genecological studies in widespread tree species have revealed steep genetic clines along environmental gradients for climate-related traits. In a changing climate, the ecological and economic importance of conifers necessitates an appraisal of how molecular genetic variation shapes quantitative trait variation, and one of the most promising approaches to answer this question is association mapping. • We phenotyped a wide collection of 410 individuals of the widely distributed conifer Sitka spruce rangewide (Picea sitchensis) for budset timing and autumn cold hardiness, and genotyped these individuals for a panel of 768 single nucleotide polymorphisms (SNPs) representing > 200 expressed nuclear genes. • After correcting for population structure, associations were detected in 28 of the candidate genes, which cumulatively explained 28 and 34% of the phenotypic variance in cold hardiness and budset, respectively. Most notable among the associations were five genes putatively involved in light signal transduction, the key pathway regulating autumn growth cessation in perennials. Many SNPs with phenotypic associations were also correlated with at least one climate variable. • This study represents a significant step toward the goal of characterizing the genomic basis of adaptation to local climate in conifers, and provides an important resource for breeding and conservation genetics in a changing climate.

Quantitative iTRAQ proteome and comparative transcriptome analysis of elicitor-induced Norway spruce (Picea abies) cells reveals elements of calcium signaling in the early conifer defense response.

Long-lived conifer trees depend on both constitutive and induced defenses for resistance against a myriad of potential pathogens and herbivores. In species of spruce (Picea spp.), several of the late events of pathogen-, insect-, or elicitor-induced defense responses have previously been characterized at the anatomical, biochemical, transcriptome, and proteome levels in stems and needles. However, accurately measuring the early events of induced cellular responses in a conifer is technically challenging due to limitations in the precise timing of induction and tissue sampling from intact trees following insect or fungal treatment. In the present study, we used the advantages of Norway spruce (Picea abies) cell suspensions combined with chitosan elicitation to investigate the early proteome response in a conifer. A combination of iTRAQ labeling and a new design of iterative sample analysis employing data-dependent exclusion lists were used for proteome analysis. This approach improved the coverage of the spruce proteome beyond that achieved in any prior study in a conifer system. Comparison of elicitor-induced proteome and transcriptome responses in Norway spruce cells consistently identified features associated with calcium-mediated signaling and response to oxidative stress that have not previously been observed in the response of intact trees to fungal attack.