Jean Bousquet

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.

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome‐wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce (Picea glauca) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations (FST) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates (HE = 0.270 and FST = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with FST at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes’ putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.

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

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. Contact: ibirol@bcgsc.ca Supplementary information: Supplementary data are available at Bioinformatics online.

RAPDs as an aid to evaluate the genetic integrity of somatic embryogenesis-derived populations of Picea mariana (Mill.) B.S.P.

SummaryThe usefulness of random amplified polymorphic DNA (RAPD) in assessing the genetic stability of somatic embryogenesis-derived populations of black spruce [Picea mariana (Mill.) B.S.P.] was evaluated. Three arbitrary 11-mer primers were successfully used to amplify DNA from both in-vivo and in-vitro material. Twenty-five embryogenic cell lines, additional zygotic embryos and megagametophytes from three controlled crosses involving four selected genotypes of black spruce were used for the segregation analysis of RAPD variants. Ten markers were genetically characterized and used to evaluate the genetic stability of somatic embryos derived from three embryogenic cell lines (one cell line per cross, 30 somatic embryos per cell line). No variation was detected within clones. The utilization of RAPD markers both for the assessment of genetic stability of clonal materials and to certify genetic stability throughout the process of somatic embryogenesis is discussed.

Variation in mitochondrial DNA reveals multiple distant glacial refugia in black spruce (Picea mariana), a transcontinental North American conifer

Range‐wide genetic variation of black spruce (Picea mariana) was studied using polymerase chain reaction‐random fragment length polymorphism markers of the mitochondrial genome. Four polymorphic mitochondrial DNA (mtDNA) loci were surveyed and two or three alleles were detected at each locus, resulting in 10 multilocus mtDNA types or mitotypes. A significant subdivision of population genetic diversity was detected (GST = 0.671; NST = 0.726), suggesting low levels of gene flow among populations. The distribution of mitotypes was not random (NST > GST; P < 0.05) and revealed four partially overlapping zones, presumably representative of different glacial populations. Comparison of the genetic structure derived from mtDNA markers and the colonization paths previously deduced from the fossil and pollen records allow us to infer at least three southern and one northeastern glacial populations for black spruce. The patterns revealed in this study suggest that black spruce shares its biogeographical history with other forest‐associated North American species.

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.

Enhancing genetic mapping of complex genomes through the design of highly-multiplexed SNP arrays: application to the large and unsequenced genomes of white spruce and black spruce

BackgroundTo explore the potential value of high-throughput genotyping assays in the analysis of large and complex genomes, we designed two highly multiplexed Illumina bead arrays using the GoldenGate SNP assay for gene mapping in white spruce (Picea glauca [Moench] Voss) and black spruce (Picea mariana [Mill.] B.S.P.).ResultsEach array included 768 SNPs, identified by resequencing genomic DNA from parents of each mapping population. For white spruce and black spruce, respectively, 69.2% and 77.1% of genotyped SNPs had valid GoldenGate assay scores and segregated in the mapping populations. For each of these successful SNPs, on average, valid genotyping scores were obtained for over 99% of progeny. SNP data were integrated to pre-existing ALFP, ESTP, and SSR markers to construct two individual linkage maps and a composite map for white spruce and black spruce genomes. The white spruce composite map contained 821 markers including 348 gene loci. Also, 835 markers including 328 gene loci were positioned on the black spruce composite map. In total, 215 anchor markers (mostly gene markers) were shared between the two species. Considering lineage divergence at least 10 Myr ago between the two spruces, interspecific comparison of homoeologous linkage groups revealed remarkable synteny and marker colinearity.ConclusionThe design of customized highly multiplexed Illumina SNP arrays appears as an efficient procedure to enhance the mapping of expressed genes and make linkage maps more informative and powerful in such species with poorly known genomes. This genotyping approach will open new avenues for co-localizing candidate genes and QTLs, partial genome sequencing, and comparative mapping across conifers.

A mitochondrial DNA minisatellite reveals the postglacial history of jack pine (Pinus banksiana), a broad‐range North American conifer

Jack pine (Pinus banksiana Lamb.) is a broadly distributed North American conifer and its current range was covered by the Laurentian ice sheet during the last glacial maximum. To infer about the history and postglacial colonization of this boreal species, range‐wide genetic variation was assessed using a new and highly variable minisatellite‐like marker of the mitochondrial genome. Among the 543 trees analysed, 14 distinct haplotypes were detected, which corresponded to different repeat numbers of the 32‐nucleotide minisatellite‐like motif. Several haplotypes were rare with limited distribution, suggesting recent mutation events during the Holocene. At the population level, an average of 2.6 haplotypes and a mean haplotype diversity (H) of 0.328 were estimated. Population subdivision of genetic diversity was quite high with GST and RST values of 0.569 and 0.472, respectively. Spatial analyses identified three relatively homogeneous groups of populations presumably representative of genetically distinct glacial populations, one west and one east of the Appalachian Mountains in the United States and a third one presumably on the unglaciated northeastern coastal area in Canada. These results indicate the significant role of the northern part of the US Appalachian Mountains as a factor of vicariance during the ice age. A fourth distinct group of populations was observed in central Québec where the continental glacier retreated last. It included populations harbouring haplotypes present into the three previous groups, and it had higher level of haplotype diversity per population (H = 0.548) and lower population differentiation (GST = 0.265), which indicates a zone of suture or secondary contact between the migration fronts of the three glacial populations. Introgression from Pinus contorta Dougl. var. latifolia Engelm. was apparent in one western population from Alberta. Altogether, these results indicate that the mitochondrial DNA variation of jack pine is geographically highly structured and it correlates well with large‐scale patterns emerging from recent phylogeographical studies of other tree boreal species in North America.

Direct evidence for biased gene diversity estimates from dominant random amplified polymorphic DNA (RAPD) fingerprints

The relevance of using dominant random amplified polymorphic DNA (RAPD) fingerprints for estimating population differentiation was investigated when typically small population sample sizes were used. Haploid sexual tissues were first used to determine genotypes at RAPD loci for 75 eastern white pines (Pinus strobus L.) representing five populations. Dominant RAPD fingerprints were then inferred from genotypic data for each individual at each locus, and gene diversity estimates from both sources of data were compared. Genotypic information at RAPD loci indicated little or no differentiation among populations, similar to allozyme loci. However, estimates of population differentiation derived from dominant RAPD fingerprints according to various common methods of analysis were generally inflated, especially when all fragments were considered. Simulations showed that an increase in loci sampling and population sample sizes did not significantly alleviate the biases observed.

Trans-species shared polymorphisms at orthologous nuclear gene loci among distant species in the conifer Picea (Pinaceae): implications for the long-term maintenance of genetic diversity in trees.

For each of three nuclear gene loci, intraspecific- as well as trans-specific shared polymorphisms were detected in DNA among three distantly related species in the genus Picea. Few fixed interspecific polymorphisms were observed. Allele genealogies did not match species phylogenies, and species lineages were not reciprocally monophyletic. Based on molecular clocks and morphological evidence from the fossil record, the divergence time between species was estimated at 13-20 million years (my), and a mutation rate of 2.23 × 10(-10) to 3.42 × 10(-10) per site per year was estimated. Large historical population sizes in excess of 100 000 were inferred, which would have delayed the fixation of polymorphisms. These numbers translated into allele coalescence times in the order of 10 to 18 my, which implies the sharing of polymorphisms since common ancestry. These results suggest that trans-species shared polymorphisms might be frequent at plant nuclear gene loci, leading to high allelic diversity. Such a trend is more likely in trees and plants characterized by ecological and life-history determinants favoring large population sizes such as an outcrossing mating system, wind pollination, and a dominant position in ecosystem. These polymorphisms also call for caution in estimating congeneric species phylogenies from nuclear gene sequences in such plant groups.