Skip to main content
Download PDF

CNVs into the wild: screening the genomes of conifer trees (Picea spp.) reveals fewer gene copy number variations in hybrids and links to adaptation

BMC Genomics volume 18, Article number: 97 (2017) Cite this article

Abstract

Background

Copy number variations (CNVs) have been linked to different phenotypes in human, including many diseases. A genome-scale understanding of CNVs is available in a few plants but none are wild species, leaving a knowledge gap regarding their genome biology and evolutionary role. We developed a reliable CNV detection method for species lacking contiguous reference genome. We selected multiple probes within 14,078 gene sequences and developed comparative genome hybridization on arrays. Gene CNVs were assessed in three full-sib families from species with 20 Gb genomes, i.e., white and black spruce, and interior spruce - a natural hybrid.

Results

We discovered hundreds of gene CNVs in each species, 3612 in total, which were enriched in functions related to stress and defense responses and narrow expression profiles, indicating a potential role in adaptation. The number of shared CNVs was in accordance with the degree of relatedness between individuals and species. The genetically mapped subset of these genes showed a wide distribution across the genome, implying numerous structural variations. The hybrid family presented significantly fewer CNVs, suggesting that the admixture of two species within one genome reduces the occurrence of CNVs.

Conclusions

The approach we developed is of particular interest in non-model species lacking a reference genome. Our findings point to a role for CNVs in adaptation. Their reduced abundance in the hybrid may limit genetic variability and evolvability of hybrids. We propose that CNVs make a qualitatively distinct contribution to adaptation which could be important for short term change.

Background

Studies of genetic variation have long focused on SNPs but over the last decade, novel insights into the genomic bases of phenotypic variation have come from analysing other types of polymorphisms including structural variations (SV). Genes with variable copy numbers have long been known in humans (e.g., [1]) but screening the entire genome for quantitatively variable DNA segments is more recent (e.g., [2, 3]). DNA copy number variations (CNVs) are a class of SV which may vary in size and are being intensively studied in the human as they are associated with a variety of phenotypes, mostly diseases [4].

Although, there is no strict consensus regarding the definition of CNVs, a few common features can be found. CNVs are typically described as DNA fragments exceeding 1 kb [3, 5, 6] with a minimum sequence identity of 90% [5] and most authors indicate that CNVs should not include transposons nor arise from transposon activity [7]. We refer to presence/absence variation (PAV) when zero copies are present in an individual genome [8], representing a subset of CNVs.

CNVs can either be inherited from the previous generations or appear de novo in an individual. The per-locus rate for genomic rearrangements was estimated between 2 × 10e−5 and 1.25 × 10e−4 in humans [4]. Mechanisms that have been identified as resulting in CNVs (for review, see [9]) include non-allelic homologous recombination [10], non-homologous end joining [11], break induced replication [12], single strand annealing [13], breakage-fusion-bridge [14] and all replication slippages of the DNA-polymerase during DNA replication and reparation [15, 16]. Despite this variety of mechanisms, the majority of CNVs are linked to the occurrence of low-copy-repeats (LCRs) that are DNA segments larger than 1 Kb, dispersed throughout the genome and share a sequence identity higher than 95%. LCRs have been reported to favor genomic rearrangements in flanking regions and described as being responsible for CNV ‘hotspots’ [17].

Gene CNVs in human have been linked to various phenotypes such as the ability to digest starch [18] and to the predisposition to diseases such breast cancer [19], among others. In plants, a SV affecting three genes was linked to pathogen resistance in soybean [20] and a gene duplication was shown to cause higher grain yield in rice [21]. Studying CNV at the genome level can even reveal evolutionary footprints, as within the primate lineage [22], which implies that inherited CNVs are subject to the same evolutionary forces as all other DNA polymorphisms. Taken together, these examples illustrate the importance of CNVs for phenotypic variation and the evolution. However, few studies have investigated CNVs and associated genes at the genome-scale in plants except for a few model species including maize [8, 23], Arabidopsis [24], soybean [25], barley [26], rice [27] and poplar [28].

A widely used approach for the discovery and analysis of CNVs at the genome scale involves comparative hybridizations on arrays of BACs or oligonucleotide sequences [2, 6, 8, 25, 29]. The method is known as aCGH (array Comparative Genome Hybridization). Oligonucleotide arrays afford rapid synthesis [29] and design flexibility, which can be advantageous to overcome the more variable signal of oligonucleotide probes compared to BACs, e.g., by using a minimum of three probes per genomic region to reduce variance [29]. As the technology has developed, statistical approaches to detect probe sets with abnormal intensities have also been implemented, most often looking for breakpoints in intensity signal between neighboring probes [29–31]. A disadvantage of these methods is that they require a high quality reference genome for the species under investigation which is often lacking for non-model organisms.

To date, gene CNVs have not been investigated in genomes as large as those encountered in most conifers (18 to 34 Gb). In this report, we developed an aCGH for white spruce (~20Gb; Picea glauca) and investigated North-American spruce species. White spruce and black spruce (Picea mariana) are largely sympatric and span from the Atlantic to the Pacific coasts in Canada. Picea glauca can hybridize with Picea Engelmanni (Additional file 1: Figure S1) to produce a hybrid known as interior spruce which is widely distributed in the Rocky Mountains [32]. Genome draft assemblies have been obtained for the diploid genomes of white spruce [33, 34] and Norway spruce [35] but are highly fragmented. However, considerable genome resources are available for white spruce including a gene catalog [36], sequence capture datasets [37], and a gene expression database [38].

In this report, we took advantage of this genomic knowledge to pursue three objectives: 1) to develop a CGH array for spruce species targeting thousands of genes and a custom analysis pipeline centered on gene sequences; 2) to determine in spruce the frequency and characteristics of CNVs and genes involved in by analysing pedigreed families; and 3) to compare results from white spruce, black spruce and interior spruce. We used our CGH array in offspring-parent comparisons and were able to detect thousands of CNVs within and across spruce families. Our findings indicated a role for CNVs in adaptation and suggested that the hybrid genome of interior spruce is less dynamic in regards to CNVs.

Results

Detection of CNV genes using CGH array and control for false discovery rate

A two-step approach was used to identify suitable aCGH probes for assessing relative gene copy number variations in comparisons between test and reference genomes (Fig. 1). We designed probes starting from a genomic assembly of the white spruce gene space which was developed from sequence-capture data and represented ~23,000 distinct genes [37] in the white spruce gene catalogue [36]. First, a set of 1 million sense and antisense probes tilled along the gene space assembly was designed by Genotypic Inc. (http://www.genotypic.co.in/) and an array was produced with the entire set of 1 million probes (Fig.1b). Second, preliminary comparative hybridizations were used to select a reduced set of genes meeting probe coverage and sequence length criteria (for details, see methods). A final set ~180,000 effective probes was used to design an array targeting 14,078 gene sequences of 500 bp or more (average 2869 bp) with a minimum of six probes (average of 12.6 probes) evenly distributed along the gene sequence; the maximum number of probes reached 114 for a single 17,259 bp gene sequence (Fig.1c).

Fig. 1
figure 1

Overview of the experimental approach deployed to investigate the occurrence of CNVs within spruce pedigrees. a) Full-sibs are produced between two mature individuals by means of controlled cross-pollination in nursery. b) The ‘female’ parent genomic DNA is employed as the reference genome while the genomic DNA of descendants are used as test genomes. c) Distribution of probes number per gene. d) Distribution of probes and corrected intensity ratios along one gene sequence (GQ027111_B03) in CNV when comparing one descendant to the #77111 parent

Full size image

The 180 k probe aCGH design was effective in discovering inherited CNVs in four different spruce families by comparing the DNA of full-sib individuals (i.e., the test genomes) to one parent (i.e., the reference genome; Fig.1a-b). Here, two families were from white spruce (WS1 and WS2) and the two other families were from black spruce (BS) and interior spruce (IS, i.e., the species hybrid), respectively; we analysed 19 individuals per family in a total of 76 test-reference hybridizations. A log ratio of hybridization fluorescence intensities was calculated for each probe following the test-reference comparisons by using R package (“limma”, [39]) and in-house Python programs taking into account the dye and GC gene content biases. The vast majority of probes gave a log intensity ratio close to 0 (Fig.2a) as expected when the reference and test genomes have the same number of copies. This expected number at a given locus is two since spruces are diploid but a higher numbers of copies cannot be dismissed. Over all comparisons, the number of positive and negative ratios were very similar.

Fig. 2
figure 2

Hybridizations and CNV detection results. a) Distribution of the log2 ratios for all probes over all families. b) Distributions of genes found in CNV among spruce families

Full size image

Several measures were taken to develop a reliable approach to identify genes displaying CNVs. The experimental design involved repeatedly assessing intensity ratios for a given gene by using several individuals from a pedigree and 12 probes per gene in average. Self-self hybridizations in which a single genome was hybridized in duplicate on the same array but with different dyes were used in addition to the test-reference comparisons. Self-self hybridizations were carried out in all three species (total of 9 self-self aCGH) and served to define robust CNV detection parameters, i.e., two major criteria for CNV detection, which resulted in a False Discovery Rate < < 1% in all comparisons over all species. First, an individual probe criterion was defined as: the absolute value of the log2 ratio must be superior to 0.42 (= |log2(4/3)| or |log2(3/4)|). Second, a gene level criterion was defined based on cross-probe repeatability such that a minimum of 83% of probes per gene must meet the log ratio criterion which translates in an average of 11 significant probes per gene declared as representing a CNV.

Inherited CNVs were found for many different genes but with low frequencies

We detected CNVs affecting a total of 3612 distinct genes when considering all comparisons and species (Fig. 2b, Table 1). As the DNA was extracted from a sample comprising several needles (folia tissue), hence large number of cells, the CNVs most likely affected the entire tested organism and were mostly representative of a variation inherited from the previous generation. CNVs with an absolute log ratio superior to 3 made up only 2.8% to 4.0% of the gene CNVs and represented the upper limit above which large copy number differences are undistinguishable from presence/absence variations (PAVs) given the maximum intensities that can be reached using aCGH. The number of individuals scored positively for any given gene CNV ranged from 1 to 19 per family (out of 19 tested descendants). This led us to classify gene CNVs as common (i.e., in 2 individuals or more) and infrequent (i.e., in only one individual per family), the latter accounted for 64% of gene CNVs on average (Table 1). The average number of gene CNVs detected within an individual was 109.8 over all families. Two individuals in each of the white spruce families and one individual in the black spruce family accounted for nearly half of the infrequent gene CNVs in each family and represented outliers within those families (Fig. 3).

Table 1 Numbers of gene CNVs within each family and percentages of frequent, infrequent and presence/absence variations (PAV)
Full size table
Fig. 3
figure 3

Numbers of gene CNVs within each individual for the four pedigrees. Despite some variation, the hybrid descendants presented fewer CNV genes than the pure species descendants

Full size image

CNVs were shared among species and less abundant in hybrid

The number of CNVs detected represented on average 10% of the tested genes but varied between families, i.e., 1338 to 1479 in white spruce, 1763 in black spruce and only 881 in interior spruce. The reduced number of gene CNVs in the hybrid (interior spruce) was found both at the individual level and the family level compared to the other spruce families (Table 1, Fig. 3). This difference was statistically tested and the number of gene CNVs was significantly lower in the interior spruce family than in white spruce and black spruce families (Tuckey HSD test, P < 1 × 10−7) while the difference between the three latter families was not significant (P > 0.69).

We examined the nature of the differences between species by considering the overlap in the gene CNV sets between the pedigrees. The different families shared from 289 to 646 gene CNVs, including a core set of 155 gene CNVs found in all spruces (Fig.2b). The number of shared gene CNVs among families was in accordance with the relative divergence between the spruce species. The two white spruce families had the most shared gene CNVs (646), white spruce and the interior spruce families had fewer (609) and the black spruce family shared the smallest numbers of gene CNVs with other pedigrees (547; Fig.2b).

Gene CNVs are enriched for adaptation-related functions and narrow expression profiles

The potential biological implications of gene CNVs were evaluated using enrichment tests for terms in Gene Ontology annotations. In white spruce, the comparison of gene CNVs to the entire set of tested genes on the CGH array revealed a significant enrichment in a number of functions related to 16 biological processes (Table 2) most significantly including the ‘response to stress’ term. Using the recently updated PFAM annotations of the white spruce gene catalog [35], enrichment was found for oxygenase, synthetase, hydrolase and cytochrome-P450 activities. In black spruce, GO-term enrichment was found in eight biological processes including most significantly ‘response to stress’ and ‘defense response’ (Table 2).

Table 2 Enrichment test in GO annotations for white and black spruce species, results from FatiGO in Babelomics 4.3 (Medina et al. 2010)
Full size table

The different numbers of gene CNVs detected between species did not translate into functional annotation differences. The core set of 155 gene CNVs found in all of the families is significantly enriched in functions related to defense response. The hybrid family sharing gene CNVs with the white spruce families may be expected by inheritance of genetic material from the white spruce species but did not explain gene CNVs also shared with black spruce. Hence, such CNVs sharing between all families is more likely indicative of ancestral genetic variations that remain variable in several North-American Picea species. In addition to defense responses-related functions, all spruce species and the hybrid were enriched in the stress responses while annotations of the 1242 genes specifically found in black spruce were not significantly enriched in any biological process compared to white spruce specific genes. On the other hand, the 199 genes specifically found in the hybrid were not enriched in functions related to defense nor stress responses but enriched in functions related to responses to hormonal, organic and chemical substances.

Expression levels for 81.5% of the gene CNVs were previously assessed in various organs at different stages of development in white spruce for a total of eight tissues [38]. Gene CNVs appeared expressed in wide range of tissues and organs, including vegetative buds, foliage, primary and secondary xylem, phelloderm, roots, megagametophyte and embryogenic cells with no significant enrichment for any of them. However, the set of gene CNVs is significantly enriched in genes expressed in only one or two tissues, with no or little variation between tissues.

CNVs were distributed across the entire genome

A dense genetic map has been continuously developed since 2006 for white spruce [40–42], which has allowed to locate the position of 1766 genes, observe a high synteny between spruce species and dissect the genomic architecture for several traits of adaptive and economic interest in spruce species [43, 44]. A total of 210 (5.8%) of the CNV genes were located on this spruce genetic map. The CNVs were distributed over the entire genome; each linkage group included from 10 (LG1) to 24 (LG10) gene CNVs and had an average of 17.5 CNVs (Fig. 4). The average distances between two mapped gene CNVs were 27.26, 18.75, 16.28 and 45.26 cM for the WS1, WS2, BS and IS families respectively, and was inversely proportional to the number of mapped gene CNVs in each family. No significant clustering of gene CNVs was observed within linkage groups although a total of 10 cases of two or more genes within a 1 cM region were found across the three spruce species (Additional file 2).

Fig. 4
figure 4

Distribution over the genome and frequency of genes found in CNV also positioned upon the spruce genetic map [40]. Each point represents the position of a gene found in CNV in one CGH comparison; the piling is indicative of the number of individuals presenting the CNV

Full size image

Discussion

We analysed thousands of gene CNVs and found that they are less abundant in a hybrid compared to two pure spruce species. The biological functions of the gene CNVs indicate a role in stress and defense responses. We discuss the development of CNV detection in wild species and the insights gained from our findings into genome biology, evolution and adaptation.

Detecting gene CNVs within pedigrees

We developed gene-based CNV detection procedures to analyse wild species, even those with large genomes. A low threshold was used for CNV detection (test to reference genome log2 ratios > 0.42) compared to other recent studies (e.g., [8, 23, 26]). This apparently low value was permitted because of an array design that included a relatively large number of probes which presented such significant signal for the same gene. Previous reports have recommended a minimum of three probes per genomic segment in order to effectively estimate the relative number of copies when using oligonucleotide CGH (e.g., [45]). When analysing genes, some authors set a requirement that all of the three to four probes tested be significant for positive scoring of CNVs (e.g., [23]). We used many more probes per gene sequence (12.6 probes on average) which provided a strong support for the estimated ratios and robust detection parameters with very low FDR. In these conditions, the number of CNVs detected was inflated substantially when we lowered the proportion of significant probes (results not shown) and was only slightly affected when we lowered the significance threshold for intensity ratios.

We did not detect a bias towards negative or positive log ratios in the overall data set nor when comparing individual offspring to a parental reference genome. This result appeared to be at odds with those from other organisms showing that losses often dominate over gains [8]; however, the detection of more losses than gains could be due to a reference genome bias, at least in some cases [8, 23]. This type of bias is not expected in our experiment as we used parent-offspring comparisons which reduced sequence divergences between the test and reference genomes, both being similarly distant from the genome used for probe design. It has also been argued that a higher number of positive CNVs (i.e., copy gains) may be expected because genomes are likely more tolerant to duplications than deletion possibly resulting in loss of function [46].

Our approach, in which we sampled multiple offspring, was designed to uncover and identify inherited pre-existing gene CNVs shared across a pedigree rather than de novo CNVs arising from recent structural variations. However, many of the CNVs in spruces (50 to 66%) were found in only one individual within a pedigree, as observed in other plants [26, 46, 47]. These infrequent CNVs may represent three different cases: i) random effects due to the small sample size (19 individuals per family), which is supported by the detection of these same CNVs in at least one other family (275 cases); ii) false-positives, despite the very small number of CNVs in self-self hybridizations; iii) de novo structural variations which occurred in or during reproduction of the previous generation. In the subset of CNVs found only once over the entire dataset (N = 1919; including the four pedigrees) the proportion of negative and positive log ratios was balanced, suggesting that there was no prevalence of de novo copy gains or losses.

Impact and distribution of CNVs across the genome

The proportion of genes positively scored for CNVs was 15.5% in white spruce (considering both families) which is similar or slightly higher in magnitude than found in whole genome analyses in model species and crop plants ([3] (human); [48] (mouse); [49] (rat); [50] (human and chimpanzee); [8, 23] (maize), [26] (barley)). The high proportion found in spruce can be either attributed to the detection of relatively short CNV regions enabled by the high probe density in our array design, or the occurrence of more than one gene located within the same long DNA sequence in CNV. In that case, some CNVs may represent only one and the same DNA segment that has been duplicated or lost. This could not be verified for most of the CNVs as their genomic positions were unknown. Hence, the rate of CNVs may not precisely reflect the rate of actual structural variations but rather provided an estimate of the order of magnitude.

In total, 9% of the genes surveyed (i.e., 1,284 out of 14,078) using aCGH and 6% of the genes found in CNV had previously been positioned on a genetic linkage map in Picea glauca [40]. This lower proportion of mapped CNVs than expected by chance alone is not surprising because the map was mainly developed using SNP from large high-throughput genotyping arrays [40] which likely selected against genes that display CNVs. Localizations on the genetic map showed that CNV genes were widely distributed at a low density over the entire genome (Fig. 4) with a minimum of 10 CNV genes per linkage group and average distance of 27 cM between successive CNVs within species. In addition, 10 cases of 2 or 3 CNV genes located within 1-cM of each other suggested the occurrence of structural variations involving chromosomal segments with more than one gene. Even though the number of mapped CNVs was too small to precisely assess the frequency of such events, it would represent an estimated rate of SVs close to 4.3%. Multi-gene SVs have already been described in soybean [20] and a rate of 5–10% was reported in maize/teosinte by using aCGH [23]. In light of these reports, more SVs affecting more than one gene in spruce could be expected. However, conifer genes are often separated by very long non-coding regions [33–35, 51] such that multi-gene SVs may be less common in this phylum. Furthermore, the discovery of CNVs in close proximity to each other raises the question of their origin. The first hypothesis may be that these genes are among the few genes actually separated by relatively small non-coding sequences. Another explanation would be that these genes are widely separated but structural variations that are typically even wider in conifers, hence allowing to encompass several coding sequences. Finally, one possible explanation is that these genes may otherwise originate from different events of SV within the same genomic region; such hotspots for structural variations have been found in other model organisms [52–54]. Further characterization of these genetically linked CNVs genes and a finer dissection of structural variations could be achieved by analysing larger populations (~100 segregating individuals), hence allowing to directly map these genetic variations and infer the occurrence of one or more events of structural variation.

Numbers of shared CNVs agrees with the estimated time of species divergence

The proportion of shared gene CNVs across the different Picea species was in accordance with their phylogenetic distance, i.e., more gene CNVs were shared among more closely related families and species. The highest overlap was found between the two white spruce families, suggesting that the 646 shared gene CNVs are not limited to one parental lineage but more likely to be spread at a significant level in the population. The second highest overlaps were between interior spruce and white spruce families despite the lower number of gene CNVs in the hybrid. Engelmann and white spruces, the parental species of the hybrid interior spruce, diverged only 9.9 Myears ago, a relatively short timeframe compared to the 38.1 Myears separating the white/Engelmann spruce lineage from the black spruce lineage [55]. In line with their relative times of divergence, Engelmann and white spruce can hybridize while black spruce is reproductively isolated from both of the other species [56, 57] and also shared the smallest number of CNVs with the other spruce families.

Evolutionary forces have also been linked to CNV diversity in other plants. In a study comparing domesticated and wild barely accessions for instance, it was shown that strong directional selection associated with breeding in barley reduced CNV diversity [26]. Hence, CNV diversity in barley mirrored the behaviour of other DNA polymorphisms indicating that evolutionary forces likely have similar impacts on CNVs than on other types of DNA polymorphism. Similarly, CNV diversity assessed using a human aCGH allowed to detect gene CNVs shared among primates that likely resulted from gene duplication in common ancestors [22]. Therefore, gene CNVs shared between the different spruce species likely represent ancestral polymorphisms that have not been eliminated by genetic drift or directional selection during the course of evolution, or alternatively, have been maintained by balancing or divergent selection.

Role of CNVs in the evolution of forest trees

A wide range of gene CNVs have been linked to quantitative phenotypic variation in model organisms [18, 20, 23, 58]. In forest trees, quantitative genetic variations have been linked to adaptation and is characteristic of key adaptive traits such as growth phenology and cold tolerance. Genetic variation in adaptive traits has been studied in many conifers including Picea spp. and was often observed to follow gradual changes over wide ranging geographic areas and environmental conditions as opposed to drastic changes linked to specific conditions [59–61]. Here we examine the potential impact of CNVs on adaptation in trees and predict that they make a qualitatively distinct contribution compared to variations in coding sequences resulting in amino acid substitutions.

The two types of mutation have different impacts at the molecular level and are expected to have different evolutionary consequences. Changes in amino acid sequences may alter protein function and therefore are more often likely to be deleterious and purged, or to be recessive and only impact phenotype in homozygotes [62]. CNVs alter gene dosage and may affect expression which is akin to cis- acting regulatory variation. As such CNVs are more likely to alter phenotype in heterozygotes and thus be more directly visible to selection [63]. Gene dosage effects caused by CNVs may also be detrimental, for example by creating an imbalance in multi-protein complexes [62], but this is likely to be restricted to certain classes of genes and be more prevalent in house-keeping functions.

The relative effects of CNVs and amino acid changes will also depend on the function and expression of the genes product. Here, CNV genes were enriched for biological functions related to stress response, defense response and metabolic processes, which was also observed in barley and maize [23, 26] and is in line with a role in adaptation. The CNV genes were also enriched for narrow transcript accumulation profiles in white spruce [38], suggesting more specialized functions rather than house-keeping roles which reinforces the potential for CNVs to contribute to adaptive phenotypes. Considering these observations, CNVs affecting genes have the potential to make a unique contribution to adaptation in trees as a result of their qualitatively distinct impact on evolution and the types of genes affected. We further suggest that CNVs could be especially important for short term adaptability and thus in responding to the acceleration of environmental change.

Fewer CNVs in the hybrid species

Significantly fewer CNV genes were found for the interior spruce family (IS), a hybrid resulting from the natural introgression of white and Engelmann spruce. We examined whether the difference could be due to technical reasons such as IS producing fewer successful hybridizations as the probe set was designed with pure white spruce sequences. However, the data indicated that the numbers of successfully tested genes were nearly identical in all three species including the more distantly related black spruce. This is consistent with the expectation that the procedure we used was able to tolerate up to three mismatches per probe (Genotypic Inc. pers. comm.) and the average degree of polymorphism found among sequences of white, black and Engelmann spruces.

The occurrence of CNVs is underpinned by mechanisms that are likely to be conserved across all living organisms [9] and are predominantly associated with low-copy-repeats (LCRs), such as unequal interchromosomal recombinations, among others. LCRs are conserved (identity > 95%) DNA segments that occur a few times throughout the genome [9]. Interior spruce has a diploid genome that originally came from contributions of white spruce and Engelmann spruce. As a result, a higher level of heterozygosity is observed in the individuals of interior spruce compared to individuals of pure species. Such dual origin is expected to have introduced a higher level of diversity in LCRs found in hybrid than that found in individuals of pure species. Hence, a lower sequence identity between LCR sequences of homologous chromosomes is expected in the hybrid genome which would have entangled the mechanisms leading to gene CNVs. As a consequence, the interior spruce individuals has less numerous CNVs than white spruce or black spruce individuals; these CNVs may result from other types of breakage-reparation mechanisms that are not based upon sequence identity. These results suggest that hybrids have systemically reduced copy number variability indicative of a less dynamic genome which in turn would limit genome evolvability. Future studies investigating other cases of introgression, either from the wild or from artificial cross-breeding, would support or dismiss this hypothesis.

Conclusion

Investigations of non-model species have broadened the scope of genomic studies and contribute to the understanding of ecological diversity. More affordable sequencing and genotyping techniques as well as the advent of Next-generation sequencing have paved the way to these developments. This report illustrates how techniques of transcriptomics and gene-space characterization in a non-model system can in turn be used to design CGH arrays and study CNVs on a genome scale. Approaches using whole genome sequencing data from model organisms to detect usually small CNVs have been developed; however different detection methods and algorithms have produced different results from the same datasets [64]. More recently, continuous development of technologies to sequence and characterize larger DNA segments have allowed to investigate larger structural variations in model organisms using third generation PacBio sequencing [65] or genome mapping on nanochannel arrays [66]. These technologies can be applied to non-model organisms, although the 20-Gb conifer genomes would still remain challenging.

An advantage of the approach outlined here was the very low FDR which indicated robust results. Although it is time and resource consuming, testing several offspring from the same parents was also helpful as it confirmed the discovery of CNVs by repeated identification in related individuals. This approach also throws light onto inherited gene CNVs which are of interest for investigating the role of CNVs in adaptation and evolution as well as establishing a framework to study their relation to phenotype, e.g., through QTL mapping. As gene CNVs could result in dosage effects on phenotypic variation, they may play a key role in quantitative trait variations such as growth and phenology in plants.

Here, we reported a lower number of CNVs within a hybrid pedigree than in families from pure species, suggesting that hybrids may systemically have a lower diversity in CNVs. Considering the gene annotations related to stress and defense responses, this reduced genetic diversity may act on phenotypic variability and evolution of hybrids in response to natural selection. Future investigations to assess the impact of CNVs on the evolution of species could focus on CNVs at the population level, looking for signatures of selection across geographical ranges and varying environmental conditions. Applied to a variety of species, such analyses may reveal new insights into adaptation and speciation.

Methods

Pedigrees & plant material

Full-sib trees issued from four controlled crosses were sampled in common garden sites. Two pedigrees (2516: 77111 × 2388; 2856: 80109 × 80112) were sampled in the Quebec province (East Canada, supporting Additional file 1: Figure S1) and represent pure white spruce trees that were already investigated and genotyped for genetic mapping and QTL mapping of adaptive traits and growth [41, 43]. The third pedigree (Pg29 × Pg79) was sampled in British-Columbia (West Canada, Additional file 1: Figure S1) and represents an introgressed progeny of white spruce and Engelmann spruce. Finally, the fourth pedigree (Pm1442 × Pm1425) was sampled in the North of the species distribution in the Quebec province (East Canada, Additional file 1: Figure S1), remotely from the sympatric zone between Picea mariana and Picea rubens, hence represent pure black spruce trees.

DNA extraction

Whole shoot tips composed of needles, buds and apex were collected during early spring and kept at minus 80 °C until further use. Tissues were ground to a fine powder using a MixerMill MM 300 (Retsch, http://www.retsch.com/). Genomic DNA (gDNA) was extracted with DNeasy Plant Mini Kit (QIAgen, http://www.qiagen.com/). DNA concentration was determined using a NanoDrop 1000 (Thermo Scientific, http://www.thermoscientific.com/). Quality was assessed by agarose gel electrophoresis (1%, TAE buffer) of 100 ng gDNA per sample.

Marking, hybridization and visualizing

Array CGH hybridizations were performed with the SureTag DNA Labeling Kit (Agilent, protocol version 7.2), following manufacturer’s instructions using the manual configuration for the 1 × 1 M array and the Tecan configuration for the 4 × 180 K array. HS400 Pro hybridization stations (Tecan, http://www.tecan.com/) were used with the Tecan configuration.

Probe and array designs from genome capture

A sequence capture assembly is available for white spruce since 2014 [37]. For this work, custom Nimblegen probes were designed from a gene catalog for white spruce [36] and included in a chip targeting 23,864 genes. The captured DNA was sequenced using 454 technology and sequences were assembled using the gsAssembler module of Newbler (v2.53). A total of 83,362 contigs (55.87 Mb total length) were obtained, only those corresponding to known Picea glauca trancripts were considered, representing ~23,000 genes and including mostly exons and introns < 1Kbp.

Probe design was realized by Genotypic Inc. (http://www.genotypic.co.in/) using these genomic sequence data after filtering out highly repeated sequences. Probes were designed in sense and antisense orientation based on GC content, Tm and PolyX, with an average probe tiling of 30 bp. Probes were validated according to these criteria: single hit against the target, alignment length of 60 bp, less than 3 allowed mismatches and less than 2 allowed gaps. After validation, cross-hybridizing probes were removed. In the end, 3,661 contigs were removed from sequence data because probes did not pass criteria. 971,042 probes were kept for the final design, representing 22,219 Picea glauca genes [36].

Final design was sent to Agilent (Agilent Technologies, http://www.agilent.com/) through their SureDesign platform (https://earray.chem.agilent.com/suredesign/) for printing on SurePrint G3 Custom CGH Microarrays 1x1M. This first CGH microarray was tested 5 times using the individuals from the different crosses, comparing the gDNA from offspring to the gDNA of its parent.

These comparisons were analyzed by Genotypic who selected high quality probes according to the following criteria. Raw intensity had to be twice background intensity in at least 6 datasets out of 10 (5 hybs × 2 channels). Sense or antisense probe targeting the same position were removed based on lowest average intensity across the 10 datasets. There had to be a minimum 6 probes per gene. The length covered (i.e., distance between start of first probe and end of last probe) had to be more than 500 bp, otherwise the gene was not considered. Probe distribution across contigs was adjusted with minimum spacing of 200 bp. However if number of probes per transcript was less than 16, spacing was reduced to 100 bp.

The final number of probes on the array is 177 223, representing 14,078 transcripts. Final design was sent to Agilent through the SureDesign platform for printing on SurePrint G3 Custom CGH Microarrays 4x180K.

Data treatment

The fluorescence intensities were analyzed using an in-house pipeline of R (v3.1.3) and Pyton (v2.7.9) scripts (available in https://bitbucket.org/jprunier/git-CNVgene-discovery) in order to correct for dye bias, take into account the GC content in the intensity variations and finally compare the ‘test’ genome (a descendant) to the ‘reference’ genome (a parent).

To take into account the technical variability among slides and cytochromes, a LOESS correction was performed using the ‘limma’ package in the R distribution [38]. Intensity ratios were normalized applying the robust LOESS regression on the M-A plot, where M = log2(I1/I2) and A = log2[(I1 _ I2)1/2], with I1 and I2 being the background subtracted intensities of the spot in the two images.

In addition, given the well-known effect of GC content on the intensities, also called “genome waves” [67], the effect of the GC content of each contig on the signal intensity was tested and correction was made to take into account such effect. In line with this, the effect of the bait (=probe) composition was also highlighted as a possible biasing factor in arrays. Hence, the composition effect on signal intensities was also tested and the intensities ratio was adjusted accordingly for each probe.

Our dataset presented an originality since it presented an average number of 12.6 probes per gene, in other words 12.6 probes per DNA regions. Thus, two parameters could be explored, i.e., the minimum absolute intensity ratio to consider a probe significant and the proportion of probes presenting such ratio to consider the gene in significantly different copy numbers between the test and reference samples. In order to delineate a set of parameters allowing to robustly identify gene CNVs, we tested series of detection thresholds and proportion of significant probes to analyze self-self hybridizations and assess the False Discovery Rate (FDR). These self-self hybridizations were performed for the ‘reference’ genomes of WS1, BS and IS families, for a total of 9 self-self aCGH tests thus covering all species and hybrid. An absolute value of the log2 ratio of 0.42 corresponding to a variation of three over four copies of the gene (|log2(3/4)| = |log2(4/3)|) as a significant threshold for each probe and a proportion of 0.85 significant probes per genes resulted in a very low FDR (<<1%).

All tested genes (14,078) were annotated according to PFAM and the Gene Ontology for homologous sequence in TAIR [36] and all enrichment tests for genes found in CNV were performed using FatiGO in Babelomics 4.3 (http://v4.babelomics.org/functional.html) [68]. Gene genomic locations were obtained from the spruce genetic map [40] and the density distributions over the genetic map for gene CNVs and all mapped genes were compared by means of Kolmogorov-Smirnov test using the ‘ks.test’ function in R (v3.1.3).

Abbreviations

aCGH:

Comparative Genomic Hybridization on array

CNV:

Copy Number Variation

PAV:

Presence/Absence Variation

QTL:

Quantitative Trait Loci

SV:

Structural Variation

References

  1. Nathans J, Piantanida TP, Eddy RL, Shows TB, Hogness DS. Molecular genetics of inherited variation in human color vision. Science. 1986;232:203–10.

    Article  CAS  PubMed  Google Scholar 

  2. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C. Detection of large-scale variation in the human genome. Nat Genet. 2004;36:949–51.

    Article  CAS  PubMed  Google Scholar 

  3. Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, et al. Global variation in copy number in the human genome. Nature. 2006;444:444–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Lupski JR. Structural variation in the human genome. N Engl J Med. 2007;356:1169–71.

    Article  CAS  PubMed  Google Scholar 

  5. Feuk L, Carson AR, Scherer SW. Structural variation in the human genome. Nat Rev Genet. 2006;7:85–97.

    Article  CAS  PubMed  Google Scholar 

  6. Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Maner S, Massa H, Walker M, Chi M, et al. Large-scale copy number polymorphism in the human genome. Science. 2004;305:525–8.

    Article  CAS  PubMed  Google Scholar 

  7. Freeman JL, Perry GH, Feuk L, Redon R, McCarroll SA, Altshuler DM, Aburatani H, Jones KW, Tyler-Smith C, Hurles ME, et al. Copy number variation: new insights in genome diversity. Genome Res. 2006;16:949–61.

    Article  CAS  PubMed  Google Scholar 

  8. Springer NM, Ying K, Fu Y, Ji T, Yeh CT, Jia Y, Wu W, Richmond T, Kitzman J, Rosenbaum H, et al. Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content. PLoS Genet. 2009;5:e1000734.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Hastings PJ, Lupski JR, Rosenberg SM, Ira G. Mechanisms of change in gene copy number. Nat Rev Genet. 2009;10:551–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lupski JR. Genomic disorders: structural features of the genome can lead to DNA rearrangements and human disease traits. Trends Genet. 1998;14:417–22.

    Article  CAS  PubMed  Google Scholar 

  11. Daley JM, Palmbos PL, Wu D, Wilson TE. Nonhomologous end joining in yeast. Annu Rev Genet. 2005;39:431–51.

    Article  CAS  PubMed  Google Scholar 

  12. Lam KW, Jeffreys AJ. Processes of de novo duplication of human alpha-globin genes. Proc Natl Acad Sci U S A. 2007;104:10950–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Sweigert SE, Carroll D. Repair and recombination of X-irradiated plasmids in Xenopus laevis oocytes. Mol Cell Biol. 1990;10:5849–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Tanaka H, Yao MC. Palindromic gene amplification-an evolutionarily conserved role for DNA inverted repeats in the genome. Nat Rev Cancer. 2009;9:216–24.

    Article  CAS  PubMed  Google Scholar 

  15. Lovett ST. Encoded errors: mutations and rearrangements mediated by misalignment at repetitive DNA sequences. Mol Microbiol. 2004;52:1243–53.

    Article  CAS  PubMed  Google Scholar 

  16. Slack A, Thornton PC, Magner DB, Rosenberg SM, Hastings PJ. On the mechanism of gene amplification induced under stress in Escherichia coli. PLoS Genet. 2006;2:e48.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Lupski JR. Genomic rearrangements and sporadic disease. Nat Genet. 2007;39:S43–7.

    Article  CAS  PubMed  Google Scholar 

  18. Perry GH, Dominy NJ, Claw KG, Lee AS, Fiegler H, Redon R, Werner J, Villanea FA, Mountain JL, Misra R, et al. Diet and the evolution of human amylase gene copy number variation. Nat Genet. 2007;39:1256–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Frank B, Bermejo JL, Hemminki K, Sutter C, Wappenschmidt B, Meindl A, Kiechle-Bahat M, Bugert P, Schmutzler RK, Bartram CR, Burwinkel B. Copy number variant in the candidate tumor suppressor gene MTUS1 and familial breast cancer risk. Carcinogenesis. 2007;28:1442–5.

    Article  CAS  PubMed  Google Scholar 

  20. Cook DE, Lee TG, Guo X, Melito S, Wang K, Bayless AM, Wang J, Hughes TJ, Willis DK, Clemente TE, et al. Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean. Science. 2012;338:1206–9.

    Article  CAS  PubMed  Google Scholar 

  21. Wang L, Xu L, Liu X, Zhang T, Li N, Hayel H, Zhang Y, Yan H, Zhao K, Liu GE, Zhang L. Copy number variation-based genome wide association study reveals additional variants contributing to meat quality in Swine. Sci Rep. 2015;5:12535.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Dumas L, Kim YH, Karimpour-Fard A, Cox M, Hopkins J, Pollack JR, Sikela JM. Gene copy number variation spanning 60 million years of human and primate evolution. Genome Res. 2007;17:1266–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Swanson-Wagner RA, Eichten SR, Kumari S, Tiffin P, Stein JC, Ware D, Springer NM. Pervasive gene content variation and copy number variation in maize and its undomesticated progenitor. Genome Res. 2010;20:1689–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. DeBolt S. Copy number variation shapes genome diversity in Arabidopsis over immediate family generational scales. Genome Biol Evol. 2010;2:441–53.

    Article  PubMed  PubMed Central  Google Scholar 

  25. McHale LK, Haun WJ, Xu WW, Bhaskar PB, Anderson JE, Hyten DL, Gerhardt DJ, Jeddeloh JA, Stupar RM. Structural variants in the soybean genome localize to clusters of biotic stress-response genes. Plant Physiol. 2012;159:1295–308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Munoz-Amatriain M, Eichten SR, Wicker T, Richmond TA, Mascher M, Steuernagel B, Scholz U, Ariyadasa R, Spannagl M, Nussbaumer T, et al. Distribution, functional impact, and origin mechanisms of copy number variation in the barley genome. Genome Biol. 2013;14:R58.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Bai Z, Chen J, Liao Y, Wang M, Liu R, Ge S, Wing RA, Chen M. The impact and origin of copy number variations in the Oryza species. BMC Genomics. 2016;17:261.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Pinosio S, Giacomello S, Faivre-Rampant P, Taylor G, Jorge V, Le Paslier MC, Zaina G, Bastien C, Cattonaro F, Marroni F, et al. Characterization of the poplar pan-genome by genome-wide identification of structural variation. Mol Biol Evol. 2016;33(10):2706–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Carter NP. Methods and strategies for analyzing copy number variation using DNA microarrays. Nat Genet. 2007;39:S16–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Olshen AB, Venkatraman ES, Lucito R, Wigler M. Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004;5:557–72.

    Article  PubMed  Google Scholar 

  31. Price TS, Regan R, Mott R, Hedman A, Honey B, Daniels RJ, Smith L, Greenfield A, Tiganescu A, Buckle V, et al. SW-ARRAY: a dynamic programming solution for the identification of copy-number changes in genomic DNA using array comparative genome hybridization data. Nucleic Acids Res. 2005;33:3455–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. De La Torre AR, Wang T, Jaquish B, Aitken SN. Adaptation and exogenous selection in a Picea glauca x Picea engelmannii hybrid zone: implications for forest management under climate change. New Phytol. 2014;201:687–99.

    Article  Google Scholar 

  33. Birol I, Raymond A, Jackman SD, Pleasance S, Coope R, Taylor GA, Yuen MM, Keeling CI, Brand D, Vandervalk BP, et al. Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data. Bioinformatics. 2013;29:1492–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Warren RL, Keeling CI, Yuen MM, Raymond A, Taylor GA, Vandervalk BP, Mohamadi H, Paulino D, Chiu R, Jackman SD, et al. Improved white spruce (Picea glauca) genome assemblies and annotation of large gene families of conifer terpenoid and phenolic defense metabolism. Plant J. 2015;83:189–212.

    Article  CAS  PubMed  Google Scholar 

  35. Nystedt B, Street NR, Wetterbom A, Zuccolo A, Lin YC, Scofield DG, Vezzi F, Delhomme N, Giacomello S, Alexeyenko A, et al. The Norway spruce genome sequence and conifer genome evolution. Nature. 2013;497:579–84.

    Article  CAS  PubMed  Google Scholar 

  36. Rigault P, Boyle B, Lepage P, Cooke JEK, Bousquet J, MacKay JJ. A white spruce gene catalog for conifer genome analyses. Plant Physiol. 2011;157:14–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Stival Sena J, Giguere I, Boyle B, Rigault P, Birol I, Zuccolo A, Ritland K, Ritland C, Bohlmann J, Jones S, et al. Evolution of gene structure in the conifer Picea glauca: a comparative analysis of the impact of intron size. BMC Plant Biol. 2014;14:95.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Raherison E, Rigault P, Caron S, Poulin PL, Boyle B, Verta JP, Giguere I, Bomal C, Bohlmann J, MacKay J. Transcriptome profiling in conifers and the PiceaGenExpress database show patterns of diversification within gene families and interspecific conservation in vascular gene expression. BMC Genomics. 2012;13:434.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Smyth GK: Limma: linear models for microarray data. In Bioinformatics and computational biology solutions using R and Bioconductor. New York: Springer; 2005. p. 397–420.

  40. Pavy N, Pelgas B, Laroche J, Rigault P, Isabel N, Bousquet J. A spruce gene map infers ancient plant genome reshuffling and subsequent slow evolution in the gymnosperm lineage leading to extant conifers. BMC Biol. 2012;10:1.

    Article  Google Scholar 

  41. Pelgas B, Beauseigle S, Achere V, Jeandroz S, Bousquet J, Isabel N. Comparative genome mapping among Picea glauca, P-mariana x P-rubens and P-abies, and correspondence with other Pinaceae. Theor Appl Genet. 2006;113:1371–93.

    Article  CAS  PubMed  Google Scholar 

  42. Pavy N, Pelgas B, Beauseigle S, Blais S, Gagnon F, Gosselin I, Lamothe M, Isabel N, Bousquet J. 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. BMC Genomics. 2008;9:21.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Pelgas B, Bousquet J, Meirmans PG, Ritland K, Isabel N. QTL mapping in white spruce: gene maps and genomic regions underlying adaptive traits across pedigrees, years and environments. BMC Genomics. 2011;12:145.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Prunier J, Pelgas B, Gagnon F, Desponts M, Isabel N, Beaulieu J, Bousquet J. The genomic architecture and association genetics of adaptive characters using a candidate SNP approach in boreal black spruce. BMC Genomics. 2013;14:368.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Carvalho B, Ouwerkerk E, Meijer GA, Ylstra B. High resolution microarray comparative genomic hybridisation analysis using spotted oligonucleotides. J Clin Pathol. 2004;57:644–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Brewer C, Holloway S, Zawalnyski P, Schinzel A, FitzPatrick D. A chromosomal duplication map of malformations: regions of suspected haplo-and triplolethality-and tolerance of segmental aneuploidy-in humans. Am J Hum Genet. 1999;64:1702–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Santuari L, Pradervand S, Amiguet-Vercher A-M, Thomas J, Dorcey E, Harshman K, Xenarios I, Juenger TE, Hardtke CS. Substantial deletion overlap among divergent Arabidopsis genomes revealed by intersection of short reads and tiling arrays. Genome Biol. 2010;11:1.

    Article  Google Scholar 

  48. Graubert TA, Cahan P, Edwin D, Selzer RR, Richmond TA, Eis PS, Shannon WD, Li X, McLeod HL, Cheverud JM. A high-resolution map of segmental DNA copy number variation in the mouse genome. PLoS Genet. 2007;3:e3.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Guryev V, Saar K, Adamovic T, Verheul M, Van Heesch SA, Cook S, Pravenec M, Aitman T, Jacob H, Shull JD. Distribution and functional impact of DNA copy number variation in the rat. Nat Genet. 2008;40:538–45.

    Article  CAS  PubMed  Google Scholar 

  50. Perry GH, Yang F, Marques-Bonet T, Murphy C, Fitzgerald T, Lee AS, Hyland C, Stone AC, Hurles ME, Tyler-Smith C. Copy number variation and evolution in humans and chimpanzees. Genome Res. 2008;18:1698–710.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Neale DB, Wegrzyn JL, Stevens KA, Zimin AV, Puiu D, Crepeau MW, Cardeno C, Koriabine M, Holtz-Morris AE, Liechty JD. Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies. Genome Biol. 2014;15:1.

    Article  Google Scholar 

  52. Perry GH, Tchinda J, McGrath SD, Zhang J, Picker SR, Cáceres AM, Iafrate AJ, Tyler-Smith C, Scherer SW, Eichler EE. Hotspots for copy number variation in chimpanzees and humans. Proc Natl Acad Sci. 2006;103:8006–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. De Smith AJ, Tsalenko A, Sampas N, Scheffer A, Yamada NA, Tsang P, Ben-Dor A, Yakhini Z, Ellis RJ, Bruhn L. Array CGH analysis of copy number variation identifies 1284 new genes variant in healthy white males: implications for association studies of complex diseases. Hum Mol Genet. 2007;16:2783–94.

    Article  PubMed  Google Scholar 

  54. Clop A, Vidal O, Amills M. Copy number variation in the genomes of domestic animals. Anim Genet. 2012;43:503–17.

    Article  CAS  PubMed  Google Scholar 

  55. He T, Pausas JG, Belcher CM, Schwilk DW, Lamont BB. Fire‐adapted traits of Pinus arose in the fiery Cretaceous. New Phytol. 2012;194:751–9.

    Article  PubMed  Google Scholar 

  56. Nienstaedt H, Zasada JC. Picea glauca (Moench) Voss, White Spruce. In Silvics of North America. Volume 1: USDA. 1990. p. 227–37.

    Google Scholar 

  57. Viereck L, Johnston W. Picea mariana (Mill.) BSP, Black spruce. In Silvics of North America. Volume 1: USDA. 1990. p. 227–37.

    Google Scholar 

  58. Knox AK, Dhillon T, Cheng H, Tondelli A, Pecchioni N, Stockinger EJ. CBF gene copy number variation at Frost Resistance-2 is associated with levels of freezing tolerance in temperate-climate cereals. Theor Appl Genet. 2010;121:21–35.

    Article  PubMed  Google Scholar 

  59. Morgenstern EK. Genetic variation in seedlings of picea mariana (Mill.) BSP. : I. Correlation with ecological factors. Silvae Genetica. 1969. p. 18.

    Google Scholar 

  60. Nienstaedt H. Inheritance and correlations of frost injury, growth, flowering, and cone characteristics in white spruce, Picea glauca (Moench) Voss. Can J For Res. 1984;15:498–504.

    Article  Google Scholar 

  61. Prunier J, Gerardi S, Laroche J, Beaulieu J, Bousquet J. Parallel and lineage-specific molecular adaptation to climate in boreal black spruce. Mol Ecol. 2012;21:4270–86.

    Article  CAS  PubMed  Google Scholar 

  62. Lynch M, Conery JS. The origins of genome complexity. Science. 2003;302:1401–4.

    Article  CAS  PubMed  Google Scholar 

  63. Wray GA. The evolutionary significance of cis-regulatory mutations. Nat Rev Genet. 2007;8:206–16.

    Article  CAS  PubMed  Google Scholar 

  64. Alkan C, Coe BP, Eichler EE. Genome structural variation discovery and genotyping. Nat Rev Genetics. 2011;12(5):363–76.

    Article  CAS  PubMed  Google Scholar 

  65. Dong J, Feng Y, Kumar D, Zhang W, Zhu T, Luo MC, Messing J. Analysis of tandem gene copies in maize chromosomal regions reconstructed from long sequence reads. Proc Natl Acad Sci U S A. 2016;2016:201608775.

    Google Scholar 

  66. Mak AC, Lai YY, Lam ET, Kwok TP, Leung AK, Poon A, Mostovoy Y, Hastie AR, Stedman W, Anantharaman T, et al. Genome-wide structural variation detection by genome mapping on nanochannel arrays. Genetics. 2016;202(1):351–62.

    Article  CAS  PubMed  Google Scholar 

  67. Marioni JC, Thorne NP, Valsesia A, Fitzgerald T, Redon R, Fiegler H, Andrews TD, Stranger BE, Lynch AG, Dermitzakis ET. Breaking the waves: improved detection of copy number variation from microarray-based comparative genomic hybridization. Genome Biol. 2007;8:1.

    Article  Google Scholar 

  68. Medina I, Carbonell J, Pulido L, Madeira SC, Goetz S, Conesa A, Tárraga J, Pascual-Montano A, Nogales-Cadenas R, Santoyo J. Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling. Nucleic Acids Res. 2010;38:W210–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We want to thank Nathalie Isabel (Canadian Forest Service, Canada Natural Resources), Barry Jaquish (Ministry of Forests, Lands and Natural Resource Operations, British Columbia, Canada) and Mireille Desponts (Ministère des Forêts, Faune et Parcs, Québec province, Canada) cooperation in supplying the samples coming from the various spruce pedigrees. We also appreciate the insights from MacKay’s lab team during the designing of the experiment and analysis, and Jérôme Laroche and Philippe Rigault during the bioinformatics pipeline development. We finally thank Ilga Porth from Université Laval (Québec city, QC, Canada) for useful commentaries on the first version of the manuscript and two anonymous reviewers for their beneficial comments.

Funding

This work was financially supported by the Canadian SmartForests and the GenAC Projects lead by Prof. John MacKay.

Author’s contributions

JP, SC and JM conceived the experiment. JM obtained the funding allowing to perform the experiment and scientific work. SC performed and supervised all laboratory work with help of JP. JP developed the bioinformatics pipeline and performed bioinformatics, statistical and enrichment analyses. JP wrote most of the manuscript with help from JM and SC. All co-authors read and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Data access

The raw datasets supporting the conclusions of this article were compressed in one file available in the GEO repository, at the accession no. GSE92329. All R and Python scripts and codes as well as a demonstration input file can be found at https://bitbucket.org/jprunier/git-CNVgene-discovery/.

Author information

Authors and Affiliations

  1. Institute for System and Integrative Biology (IBIS), Université Laval, Quebec, QC, G1V 0A6, Canada

    Julien Prunier & Sébastien Caron

  2. Centre for Forest Research, Université Laval, Quebec, QC, G1V 0A6, Canada

    Julien Prunier & John MacKay

  3. Department of Plant Sciences, University of Oxford, Oxford, OX1 3RB, UK

    John MacKay

Authors
  1. Julien Prunier
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Sébastien Caron
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. John MacKay
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Julien Prunier.

Additional files

Additional file 1: Figure S1.

Map of parent provenances for the four pedigrees. (DOCX 240 kb)

Additional file 2: Table S1.

Distribution of the gene CNVs along the spruce genetic map with position in cM. (XLSX 15 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Prunier, J., Caron, S. & MacKay, J. CNVs into the wild: screening the genomes of conifer trees (Picea spp.) reveals fewer gene copy number variations in hybrids and links to adaptation. BMC Genomics 18, 97 (2017). https://doi.org/10.1186/s12864-016-3458-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-016-3458-8

Keywords

  • Copy number variation
  • Comparative genomic hybridization
  • Non-model species
  • Genomic structural variation
  • Genome architecture
  • Species hybridization
  • Picea glauca
  • Picea mariana
  • Picea Engelmanni
  • conifers

BMC Genomics

ISSN: 1471-2164

Contact us