1
|
Lo T, Coombe L, Gagalova KK, Marr A, Warren RL, Kirk H, Pandoh P, Zhao Y, Moore RA, Mungall AJ, Ritland C, Pavy N, Jones SJM, Bohlmann J, Bousquet J, Birol I, Thomson A. Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response. G3 (BETHESDA, MD.) 2023; 14:jkad247. [PMID: 37875130 PMCID: PMC10755193 DOI: 10.1093/g3journal/jkad247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Revised: 05/17/2023] [Accepted: 10/09/2023] [Indexed: 10/26/2023]
Abstract
Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.
Collapse
Affiliation(s)
- Theodora Lo
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Lauren Coombe
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Kristina K Gagalova
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Alex Marr
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - René L Warren
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Heather Kirk
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Pawan Pandoh
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Yongjun Zhao
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Richard A Moore
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Andrew J Mungall
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Carol Ritland
- Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
| | - Nathalie Pavy
- Canada Research Chair in Forest Genomics, Laval University, Quebec City, QC G1V 0A6, Canada
| | - Steven J M Jones
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Joerg Bohlmann
- Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
- Department of Botany, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
| | - Jean Bousquet
- Canada Research Chair in Forest Genomics, Laval University, Quebec City, QC G1V 0A6, Canada
| | - Inanç Birol
- Canada’s Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC V5Z 4S6, Canada
| | - Ashley Thomson
- Faculty of Natural Resources Management, Lakehead University, Thunder Bay, ON P7B 5E1, Canada
| |
Collapse
|
2
|
Girardin MP, Isabel N, Guo XJ, Lamothe M, Duchesne I, Lenz P. Annual aboveground carbon uptake enhancements from assisted gene flow in boreal black spruce forests are not long-lasting. Nat Commun 2021; 12:1169. [PMID: 33608515 PMCID: PMC7895975 DOI: 10.1038/s41467-021-21222-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 01/12/2021] [Indexed: 01/31/2023] Open
Abstract
Assisted gene flow between populations has been proposed as an adaptive forest management strategy that could contribute to the sequestration of carbon. Here we provide an assessment of the mitigation potential of assisted gene flow in 46 populations of the widespread boreal conifer Picea mariana, grown in two 42-year-old common garden experiments and established in contrasting Canadian boreal regions. We use a dendroecological approach taking into account phylogeographic structure to retrospectively analyse population phenotypic variability in annual aboveground net primary productivity (NPP). We compare population NPP phenotypes to detect signals of adaptive variation and/or the presence of phenotypic clines across tree lifespans, and assess genotype-by-environment interactions by evaluating climate and NPP relationships. Our results show a positive effect of assisted gene flow for a period of approximately 15 years following planting, after which there was little to no effect. Although not long lasting, well-informed assisted gene flow could accelerate the transition from carbon source to carbon sink after disturbance.
Collapse
Affiliation(s)
- Martin P. Girardin
- grid.146611.50000 0001 0775 5922Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Québec, QC Canada ,grid.38678.320000 0001 2181 0211Centre d’étude de la forêt, Université du Québec à Montréal, Montréal, QC Canada
| | - Nathalie Isabel
- grid.146611.50000 0001 0775 5922Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Québec, QC Canada ,grid.23856.3a0000 0004 1936 8390Canada Research Chair in Forest Genomics, Faculté de Foresterie, de Géographie et de Géomatique, Université Laval, Québec, QC Canada
| | - Xiao Jing Guo
- grid.146611.50000 0001 0775 5922Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Québec, QC Canada
| | - Manuel Lamothe
- grid.146611.50000 0001 0775 5922Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Québec, QC Canada
| | - Isabelle Duchesne
- grid.202033.00000 0001 2295 5236Natural Resources Canada, Canadian Wood Fibre Centre, Québec, QC Canada
| | - Patrick Lenz
- grid.23856.3a0000 0004 1936 8390Canada Research Chair in Forest Genomics, Faculté de Foresterie, de Géographie et de Géomatique, Université Laval, Québec, QC Canada ,grid.202033.00000 0001 2295 5236Natural Resources Canada, Canadian Wood Fibre Centre, Québec, QC Canada
| |
Collapse
|
3
|
Uckele KA, Adams RP, Schwarzbach AE, Parchman TL. Genome-wide RAD sequencing resolves the evolutionary history of serrate leaf Juniperus and reveals discordance with chloroplast phylogeny. Mol Phylogenet Evol 2020; 156:107022. [PMID: 33242585 DOI: 10.1016/j.ympev.2020.107022] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2020] [Revised: 10/06/2020] [Accepted: 11/17/2020] [Indexed: 12/22/2022]
Abstract
Juniper (Juniperus) is an ecologically important conifer genus of the Northern Hemisphere, the members of which are often foundational tree species of arid regions. The serrate leaf margin clade is native to topologically variable regions in North America, where hybridization has likely played a prominent role in their diversification. Here we use a reduced-representation sequencing approach (ddRADseq) to generate a phylogenomic data set for 68 accessions representing all 22 species in the serrate leaf margin clade, as well as a number of close and distant relatives, to improve understanding of diversification in this group. Phylogenetic analyses using three methods (SVDquartets, maximum likelihood, and Bayesian) yielded highly congruent and well-resolved topologies. These phylogenies provided improved resolution relative to past analyses based on Sanger sequencing of nuclear and chloroplast DNA, and were largely consistent with taxonomic expectations based on geography and morphology. Calibration of a Bayesian phylogeny with fossil evidence produced divergence time estimates for the clade consistent with a late Oligocene origin in North America, followed by a period of elevated diversification between 12 and 5 Mya. Comparison of the ddRADseq phylogenies with a phylogeny based on Sanger-sequenced chloroplast DNA revealed five instances of pronounced discordance, illustrating the potential for chloroplast introgression, chloroplast transfer, or incomplete lineage sorting to influence organellar phylogeny. Our results improve understanding of the pattern and tempo of diversification in Juniperus, and highlight the utility of reduced-representation sequencing for resolving phylogenetic relationships in non-model organisms with reticulation and recent divergence.
Collapse
Affiliation(s)
- Kathryn A Uckele
- Department of Biology, MS 314, University of Nevada, Reno, Max Fleischmann Agriculture Building, 1664 N Virginia St., Reno, NV 89557, USA.
| | - Robert P Adams
- Baylor University, Utah Lab, 201 N 5500 W, Hurricane, UT 84790, USA.
| | - Andrea E Schwarzbach
- Department of Health and Biomedical Sciences, University of Texas - Rio Grande Valley, 1 W University Drive, Brownsville, TX 78520, USA.
| | - Thomas L Parchman
- Department of Biology, MS 314, University of Nevada, Reno, Max Fleischmann Agriculture Building, 1664 N Virginia St., Reno, NV 89557, USA.
| |
Collapse
|
4
|
Napier JD, de Lafontaine G, Hu FS. Exploring genomic variation associated with drought stress in Picea mariana populations. Ecol Evol 2020; 10:9271-9282. [PMID: 32953060 PMCID: PMC7487243 DOI: 10.1002/ece3.6614] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Revised: 06/01/2020] [Accepted: 07/08/2020] [Indexed: 01/05/2023] Open
Abstract
Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype-by-environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.
Collapse
Affiliation(s)
- Joseph D. Napier
- Department of Plant BiologyUniversity of IllinoisUrbanaILUSA
- Department of Integrative BiologyThe University of Texas at AustinAustinTXUSA
| | - Guillaume de Lafontaine
- Department of Plant BiologyUniversity of IllinoisUrbanaILUSA
- Canada Research Chair in Integrative Biology of Northern FloraUniversité du Québec à RimouskiRimouskiQCCanada
| | - Feng Sheng Hu
- Department of GeologyUniversity of IllinoisUrbanaILUSA
- Department of BiologyWashington UniversitySt. LouisMOUSA
- Department of Earth and Planetary SciencesWashington UniversitySt. LouisMOUSA
| |
Collapse
|
5
|
Napier JD, Fernandez MC, de Lafontaine G, Hu FS. Ice-age persistence and genetic isolation of the disjunct distribution of larch in Alaska. Ecol Evol 2020; 10:1692-1702. [PMID: 32076544 PMCID: PMC7029070 DOI: 10.1002/ece3.6031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 12/18/2019] [Accepted: 01/06/2020] [Indexed: 01/07/2023] Open
Abstract
Larix laricina (eastern larch, tamarack) is a transcontinental North American conifer with a prominent disjunction in the Yukon isolating the Alaskan distribution from the rest of its range. We investigate whether in situ persistence during the last glacial maximum (LGM) or long-distance postglacial migration from south of the ice sheets resulted in the modern-day Alaskan distribution. We analyzed variation in three chloroplast DNA regions of 840 trees from a total of 69 populations (24 new sampling sites situated on both sides of the Yukon range disjunction pooled with 45 populations from a published source) and conducted ensemble species distribution modeling (SDM) throughout Canada and United States to hindcast the potential range of L. laricina during the LGM. We uncovered the genetic signature of a long-term isolation of larch populations in Alaska, identifying three endemic chlorotypes and low levels of genetic diversity. Range-wide analysis across North America revealed the presence of a distinct Alaskan lineage. Postglacial gene flow across the Yukon divide was unidirectional, from Alaska toward previously glaciated Canadian regions, and with no evidence of immigration into Alaska. Hindcast SDM indicates one of the broadest areas of past climate suitability for L. laricina existed in central Alaska, suggesting possible in situ persistence of larch in Alaska during the LGM. Our results provide the first unambiguous evidence for the long-term isolation of L. laricina in Alaska that extends beyond the last glacial period and into the present interglacial period. The lack of gene flow into Alaska along with the overall probability of larch occurrence in Alaska being currently lower than during the LGM suggests that modern-day Alaskan larch populations are isolated climate relicts of broader glacial distributions, and so are particularly vulnerable to current warming trends.
Collapse
Affiliation(s)
- Joseph D. Napier
- Department of Plant BiologyUniversity of IllinoisUrbanaILUSA
- Department of Integrative BiologyThe University of Texas at AustinAustinTXUSA
| | | | - Guillaume de Lafontaine
- Canada Research Chair in Integrative Biology of Northern FloraUniversité du Québec à RimouskiRimouskiQCCanada
| | - Feng Sheng Hu
- Department of Plant BiologyUniversity of IllinoisUrbanaILUSA
- Department of GeologyUniversity of IllinoisUrbanaILUSA
| |
Collapse
|
6
|
Godbout J, Gros-Louis M, Lamothe M, Isabel N. Going with the flow: Intraspecific variation may act as a natural ally to counterbalance the impacts of global change for the riparian species Populus deltoides. Evol Appl 2020; 13:176-194. [PMID: 31892951 PMCID: PMC6935597 DOI: 10.1111/eva.12854] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2019] [Revised: 07/16/2019] [Accepted: 07/31/2019] [Indexed: 12/19/2022] Open
Abstract
The speed and magnitude of global change will have major impacts on riparian ecosystems, thereby leading to greater forest vulnerability. Assessing species' adaptive capacities to provide relevant information for vulnerability assessments remains challenging, especially for nonmodel species like the North American Populus deltoides W. Bartram ex Marshall. The objective of this study was to understand how genomic diversity of this foundation species was shaped by its environment (climate, soil, and biotic interactions) to gauge its adaptive capacity. We used two complementary approaches to get a full portrait of P. deltoides genetic diversity at both the species and whole-genome ranges. First, we used a set of 93 nuclear and three chloroplastic SNP markers in 946 individuals covering most of the species' natural distribution. Then, to measure the degree of intraspecific divergence at the whole-genome level and to support the outlier and genomic-environment association analyses, we used a sequence capture approach on DNA pools. Three distinct lineages for P. deltoides were detected, and their current distribution was associated with abiotic and biotic variations. The comparison between both cpDNA and ncDNA patterns showed that gene flow between the lineages is unbalanced. The southern and northeastern populations may benefit from the input, through river flow, of novel alleles located upstream to their local gene pools. These alleles could migrate from populations that are already adapted to conditions that fit the predicted climates in the receiving local populations, hotter at the northeastern limit and drier in the Central United States. These "preadapted" incoming alleles may help to cope with maladaptation in populations facing changing conditions.
Collapse
Affiliation(s)
- Julie Godbout
- Ministère des Forêts, de la Faune et des Parcs, Direction de la recherche forestièreQuébecQCCanada
- Canadian Forest Service, Laurentian Forestry CentreNatural Resources CanadaQuébecQCCanada
| | | | - Manuel Lamothe
- Canadian Forest Service, Laurentian Forestry CentreNatural Resources CanadaQuébecQCCanada
| | - Nathalie Isabel
- Canadian Forest Service, Laurentian Forestry CentreNatural Resources CanadaQuébecQCCanada
| |
Collapse
|
7
|
Lait LA, Hebert PDN. Phylogeographic structure in three North American tent caterpillar species (Lepidoptera: Lasiocampidae): Malacosoma americana, M. californica, and M. disstria. PeerJ 2018; 6:e4479. [PMID: 29576956 PMCID: PMC5863710 DOI: 10.7717/peerj.4479] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2017] [Accepted: 02/19/2018] [Indexed: 11/23/2022] Open
Abstract
While phylogeographic structure has been examined in many North American vertebrate species, insects have received much less attention despite their central ecological roles. The moth genus Malacosoma (Hübner, 1820), is an important group of forestry pests responsible for large-scale defoliation across much of the Nearctic and Palearctic. The present study uses sequence variation in the mitochondrial cytochrome c oxidase 1 (COI) gene to examine the population genetic structure of the three widespread Malacosoma species (M. americana, M. californica, and M. disstria). Populations of all three species showed highest diversity in the south, suggesting that modern populations derived from southern refugia with loss of variation as these lineages dispersed northwards. However, despite similar life histories and dispersal abilities, the extent of regional variation varied among the taxa. M. americana, a species restricted to eastern North America, showed much less genetic structure than the western M. californica or the widespread M. disstria. The regional differentiation in the latter reflects the likely derivation of modern lineages from several refugia, as well as taxonomic uncertainty in M. californica. In these respects, the three species of Malacosoma share phylogeographic patterns similar to those detected in vertebrates which are characterised by greater phylogeographic breaks in the western half of the continent and limited structure in the east.
Collapse
Affiliation(s)
- Linda A Lait
- Centre for Biodiversity Genomics, University of Guelph, Guelph, Ontario, Canada
| | - Paul D N Hebert
- Centre for Biodiversity Genomics, University of Guelph, Guelph, Ontario, Canada
| |
Collapse
|
8
|
Case AL, Finseth FR, Barr CM, Fishman L. Selfish evolution of cytonuclear hybrid incompatibility in Mimulus. Proc Biol Sci 2017; 283:rspb.2016.1493. [PMID: 27629037 DOI: 10.1098/rspb.2016.1493] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2016] [Accepted: 08/23/2016] [Indexed: 11/12/2022] Open
Abstract
Intraspecific coevolution between selfish elements and suppressors may promote interspecific hybrid incompatibility, but evidence of this process is rare. Here, we use genomic data to test alternative models for the evolution of cytonuclear hybrid male sterility in Mimulus In hybrids between Iron Mountain (IM) Mimulus guttatus × Mimulus nasutus, two tightly linked M. guttatus alleles (Rf1/Rf2) each restore male fertility by suppressing a local mitochondrial male-sterility gene (IM-CMS). Unlike neutral models for the evolution of hybrid incompatibility loci, selfish evolution predicts that the Rf alleles experienced strong selection in the presence of IM-CMS. Using whole-genome sequences, we compared patterns of population-genetic variation in Rf at IM to a neighbouring population that lacks IM-CMS. Consistent with local selection in the presence of IM-CMS, the Rf region shows elevated FST, high local linkage disequilibrium and a distinct haplotype structure at IM, but not at Cone Peak (CP), suggesting a recent sweep in the presence of IM-CMS. In both populations, Rf2 exhibited lower polymorphism than other regions, but the low-diversity outliers were different between CP and IM. Our results confirm theoretical predictions of ubiquitous cytonuclear conflict in plants and provide a population-genetic mechanism for the evolution of a common form of hybrid incompatibility.
Collapse
Affiliation(s)
- Andrea L Case
- Department of Biological Sciences, Kent State University, Kent, OH 44242, USA
| | - Findley R Finseth
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
| | - Camille M Barr
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
| | - Lila Fishman
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
| |
Collapse
|
9
|
Dohms KM, Graham BA, Burg TM. Multilocus genetic analyses and spatial modeling reveal complex population structure and history in a widespread resident North American passerine ( Perisoreus canadensis). Ecol Evol 2017; 7:9869-9889. [PMID: 29238522 PMCID: PMC5723591 DOI: 10.1002/ece3.3478] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2017] [Revised: 08/08/2017] [Accepted: 08/19/2017] [Indexed: 11/29/2022] Open
Abstract
An increasing body of studies of widely distributed, high latitude species shows a variety of refugial locations and population genetic patterns. We examined the effects of glaciations and dispersal barriers on the population genetic patterns of a widely distributed, high latitude, resident corvid, the gray jay (Perisoreus canadensis), using the highly variable mitochondrial DNA (mtDNA) control region and microsatellite markers combined with species distribution modeling. We sequenced 914 bp of mtDNA control region for 375 individuals from 37 populations and screened seven loci for 402 individuals from 27 populations across the gray jay range. We used species distribution modeling and a range of phylogeographic analyses (haplotype diversity, ΦST, SAMOVA, FST, Bayesian clustering analyses) to examine evolutionary history and population genetic structure. MtDNA and microsatellite markers revealed significant genetic differentiation among populations with high concordance between markers. Paleodistribution models supported at least five potential areas of suitable gray jay habitat during the last glacial maximum and revealed distributions similar to the gray jay's contemporary during the last interglacial. Colonization from and prolonged isolation in multiple refugia is evident. Historical climatic fluctuations, the presence of multiple dispersal barriers, and highly restricted gene flow appear to be responsible for strong genetic diversification and differentiation in gray jays.
Collapse
Affiliation(s)
- Kimberly M Dohms
- Department of Biological Sciences University of Lethbridge Lethbridge AB Canada.,Present address: Canadian Wildlife Service, Environment and Climate Change Canada Delta BC Canada
| | - Brendan A Graham
- Department of Biological Sciences University of Windsor Windsor ON Canada
| | - Theresa M Burg
- Department of Biological Sciences University of Lethbridge Lethbridge AB Canada
| |
Collapse
|
10
|
Lenz PRN, Beaulieu J, Mansfield SD, Clément S, Desponts M, Bousquet J. Factors affecting the accuracy of genomic selection for growth and wood quality traits in an advanced-breeding population of black spruce (Picea mariana). BMC Genomics 2017; 18:335. [PMID: 28454519 PMCID: PMC5410046 DOI: 10.1186/s12864-017-3715-5] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2016] [Accepted: 04/21/2017] [Indexed: 11/11/2022] Open
Abstract
Background Genomic selection (GS) uses information from genomic signatures consisting of thousands of genetic markers to predict complex traits. As such, GS represents a promising approach to accelerate tree breeding, which is especially relevant for the genetic improvement of boreal conifers characterized by long breeding cycles. In the present study, we tested GS in an advanced-breeding population of the boreal black spruce (Picea mariana [Mill.] BSP) for growth and wood quality traits, and concurrently examined factors affecting GS model accuracy. Results The study relied on 734 25-year-old trees belonging to 34 full-sib families derived from 27 parents and that were established on two contrasting sites. Genomic profiles were obtained from 4993 Single Nucleotide Polymorphisms (SNPs) representative of as many gene loci distributed among the 12 linkage groups common to spruce. GS models were obtained for four growth and wood traits. Validation using independent sets of trees showed that GS model accuracy was high, related to trait heritability and equivalent to that of conventional pedigree-based models. In forward selection, gains per unit of time were three times higher with the GS approach than with conventional selection. In addition, models were also accurate across sites, indicating little genotype-by-environment interaction in the area investigated. Using information from half-sibs instead of full-sibs led to a significant reduction in model accuracy, indicating that the inclusion of relatedness in the model contributed to its higher accuracies. About 500 to 1000 markers were sufficient to obtain GS model accuracy almost equivalent to that obtained with all markers, whether they were well spread across the genome or from a single linkage group, further confirming the implication of relatedness and potential long-range linkage disequilibrium (LD) in the high accuracy estimates obtained. Only slightly higher model accuracy was obtained when using marker subsets that were identified to carry large effects, indicating a minor role for short-range LD in this population. Conclusions This study supports the integration of GS models in advanced-generation tree breeding programs, given that high genomic prediction accuracy was obtained with a relatively small number of markers due to high relatedness and family structure in the population. In boreal spruce breeding programs and similar ones with long breeding cycles, much larger gain per unit of time can be obtained from genomic selection at an early age than by the conventional approach. GS thus appears highly profitable, especially in the context of forward selection in species which are amenable to mass vegetative propagation of selected stock, such as spruces.
Collapse
Affiliation(s)
- Patrick R N Lenz
- Canadian Wood Fibre Centre, Canadian Forest Service, Natural Resources Canada, Government of Canada, 1055 du PEPS, P.O. Box 10380, Québec, Québec, G1V 4C7, Canada. .,Canada Research Chair in Forest Genomics, Institute of Systems and Integrative Biology and Centre for Forest Research, Université Laval, 1030, Avenue de la Médecine, Québec, Québec, G1V 0A6, Canada.
| | - Jean Beaulieu
- Canadian Wood Fibre Centre, Canadian Forest Service, Natural Resources Canada, Government of Canada, 1055 du PEPS, P.O. Box 10380, Québec, Québec, G1V 4C7, Canada.,Canada Research Chair in Forest Genomics, Institute of Systems and Integrative Biology and Centre for Forest Research, Université Laval, 1030, Avenue de la Médecine, Québec, Québec, G1V 0A6, Canada
| | - Shawn D Mansfield
- Department of Wood Science, Forest Sciences Centre, University of British Columbia, Vancouver, British Columbia, V6T 1Z4, Canada
| | - Sébastien Clément
- Canadian Wood Fibre Centre, Canadian Forest Service, Natural Resources Canada, Government of Canada, 1055 du PEPS, P.O. Box 10380, Québec, Québec, G1V 4C7, Canada
| | - Mireille Desponts
- Ministère des Forêts, de la Faune et des Parcs, Gouvernement du Québec, Direction de la recherche forestière, 2700 rue Einstein, Québec, Québec, G1P 3W8, Canada
| | - Jean Bousquet
- Canada Research Chair in Forest Genomics, Institute of Systems and Integrative Biology and Centre for Forest Research, Université Laval, 1030, Avenue de la Médecine, Québec, Québec, G1V 0A6, Canada
| |
Collapse
|
11
|
Dering M, Kosiński P, Wyka TP, Pers-Kamczyc E, Boratyński A, Boratyńska K, Reich PB, Romo A, Zadworny M, Żytkowiak R, Oleksyn J. Tertiary remnants and Holocene colonizers: Genetic structure and phylogeography of Scots pine reveal higher genetic diversity in young boreal than in relict Mediterranean populations and a dual colonization of Fennoscandia. DIVERS DISTRIB 2017. [DOI: 10.1111/ddi.12546] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Monika Dering
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
| | - Piotr Kosiński
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
- Department of Botany; Poznań University of Life Sciences; Poznań Poland
| | - Tomasz P. Wyka
- General Botany Laboratory; Institute of Experimental Biology; Department of Biology; Adam Mickiewicz University; Poznań Poland
| | - Emilia Pers-Kamczyc
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
| | - Adam Boratyński
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
| | - Krystyna Boratyńska
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
| | - Peter B. Reich
- Department of Forest Resources; University of Minnesota; St. Paul MN USA
- Hawkesbury Institute for the Environment; Western Sydney University; Penrith NSW Australia
| | - Angel Romo
- Institute of Botany; Consejo Superior de Investigaciones Científicas; IBB-CSIC-ICUB; Barcelona Spain
| | - Marcin Zadworny
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
| | - Roma Żytkowiak
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
| | - Jacek Oleksyn
- Laboratory of Systematics and Geography; Institute of Dendrology; Polish Academy of Sciences; Kórnik Poland
- Department of Forest Resources; University of Minnesota; St. Paul MN USA
| |
Collapse
|
12
|
Jackman SD, Warren RL, Gibb EA, Vandervalk BP, Mohamadi H, Chu J, Raymond A, Pleasance S, Coope R, Wildung MR, Ritland CE, Bousquet J, Jones SJM, Bohlmann J, Birol I. Organellar Genomes of White Spruce (Picea glauca): Assembly and Annotation. Genome Biol Evol 2015; 8:29-41. [PMID: 26645680 PMCID: PMC4758241 DOI: 10.1093/gbe/evv244] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The genome sequences of the plastid and mitochondrion of white spruce (Picea glauca) were assembled from whole-genome shotgun sequencing data using ABySS. The sequencing data contained reads from both the nuclear and organellar genomes, and reads of the organellar genomes were abundant in the data as each cell harbors hundreds of mitochondria and plastids. Hence, assembly of the 123-kb plastid and 5.9-Mb mitochondrial genomes were accomplished by analyzing data sets primarily representing low coverage of the nuclear genome. The assembled organellar genomes were annotated for their coding genes, ribosomal RNA, and transfer RNA. Transcript abundances of the mitochondrial genes were quantified in three developmental tissues and five mature tissues using data from RNA-seq experiments. C-to-U RNA editing was observed in the majority of mitochondrial genes, and in four genes, editing events were noted to modify ACG codons to create cryptic AUG start codons. The informatics methodology presented in this study should prove useful to assemble organellar genomes of other plant species using whole-genome shotgun sequencing data.
Collapse
Affiliation(s)
- Shaun D Jackman
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - René L Warren
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Ewan A Gibb
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Benjamin P Vandervalk
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Hamid Mohamadi
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Justin Chu
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Anthony Raymond
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Stephen Pleasance
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Robin Coope
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada
| | - Mark R Wildung
- School of Molecular Biosciences, Washington State University
| | - Carol E Ritland
- Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, BC, Canada
| | - Jean Bousquet
- Department of Forest and Environmental Genomics, Université Laval, Québec, QC, Canada
| | - Steven J M Jones
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada Department of Medical Genetics, University of British Columbia, Vancouver, BC, Canada School of Computing Science, Simon Fraser University, Burnaby, BC, Canada
| | - Joerg Bohlmann
- Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, BC, Canada Michael Smith Laboratories, University of British Columbia, Vancouver, BC, Canada Department of Botany, University of British Columbia, Vancouver, BC, Canada
| | - Inanç Birol
- Canada's Michael Smith Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, BC, Canada Department of Medical Genetics, University of British Columbia, Vancouver, BC, Canada School of Computing Science, Simon Fraser University, Burnaby, BC, Canada Department of Computer Science, University of British Columbia, Vancouver, BC, Canada
| |
Collapse
|
13
|
Mitochondrial introgression and complex biogeographic history of the genus Picea. Mol Phylogenet Evol 2015; 93:63-76. [DOI: 10.1016/j.ympev.2015.07.020] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2015] [Revised: 07/18/2015] [Accepted: 07/25/2015] [Indexed: 11/18/2022]
|
14
|
Prunier J, Tessier G, Bousquet J, MacKay J. From genotypes to phenotypes: expression levels of genes encompassing adaptive SNPs in black spruce. PLANT CELL REPORTS 2015; 34:2111-2125. [PMID: 26260097 DOI: 10.1007/s00299-015-1855-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/16/2015] [Revised: 07/30/2015] [Accepted: 08/04/2015] [Indexed: 06/04/2023]
Abstract
Measuring transcript levels for adaptive genes revealed polymorphisms having cis -effect upon gene expression levels related to phenotype variation in a black spruce natural population. Trees growing in temperate and boreal regions must acclimate to changes in climatic factors such as low winter temperatures to survive to seasonal variations. Common garden studies have shown that genetic variation in quantitative traits helps species to survive and adapt to environmental changes and local conditions. Twenty-four genes carrying SNPs were previously associated with genetic adaptation in black spruce (Picea mariana [Mill.] BSP). The objectives of this study were to investigate the potential role of these genes in regulation of winter acclimation and adaptation by studying their patterns of expression as a function of the physiological stage during the annual growth cycle, tissue type, and their SNP genotypic class. Considerable variability in gene expression was observed between different vegetative tissues or organs, and between physiological stages. The genes were expressed predominantly in tissues that could be linked more directly to winter acclimation and adaptation. The expression levels of several of the genes were significantly related to variation in tree height growth or budset timing and expression level variation related to SNP genotypic classes was observed in four of the genes. An interaction between genotypic classes and physiological stages was also observed for some genes, indicating genotypes with different reaction norms in terms of gene expression.
Collapse
Affiliation(s)
- Julien Prunier
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for System and Integrative Biology, Université Laval, 1030 Avenue de la Médecine, Québec, G1V0A6, Canada.
| | - Guillaume Tessier
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for System and Integrative Biology, Université Laval, 1030 Avenue de la Médecine, Québec, G1V0A6, Canada
| | - Jean Bousquet
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for System and Integrative Biology, Université Laval, 1030 Avenue de la Médecine, Québec, G1V0A6, Canada
| | - John MacKay
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for System and Integrative Biology, Université Laval, 1030 Avenue de la Médecine, Québec, G1V0A6, Canada
- Department of Plant Sciences, University of Oxford, South Parks Road, Oxford, OX1 3RB, UK
| |
Collapse
|
15
|
de Lafontaine G, Prunier J, Gérardi S, Bousquet J. Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor-derivative spruces (Picea mariana×P. rubens). Mol Ecol 2015; 24:5229-47. [DOI: 10.1111/mec.13377] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2015] [Revised: 08/05/2015] [Accepted: 09/01/2015] [Indexed: 01/17/2023]
Affiliation(s)
- Guillaume de Lafontaine
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute of Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC G1V 0A6 Canada
| | - Julien Prunier
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute of Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC G1V 0A6 Canada
| | - Sébastien Gérardi
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute of Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC G1V 0A6 Canada
| | - Jean Bousquet
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute of Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC G1V 0A6 Canada
| |
Collapse
|
16
|
Cinget B, de Lafontaine G, Gérardi S, Bousquet J. Integrating phylogeography and paleoecology to investigate the origin and dynamics of hybrid zones: insights from two widespread North American firs. Mol Ecol 2015; 24:2856-70. [DOI: 10.1111/mec.13194] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2014] [Revised: 04/02/2015] [Accepted: 04/08/2015] [Indexed: 12/20/2022]
Affiliation(s)
- Benjamin Cinget
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute for Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC Canada G1V0A6
| | - Guillaume de Lafontaine
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute for Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC Canada G1V0A6
| | - Sébastien Gérardi
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute for Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC Canada G1V0A6
| | - Jean Bousquet
- Canada Research Chair in Forest and Environmental Genomics; Centre for Forest Research and Institute for Systems and Integrative Biology; Université Laval; 1030 Avenue de la Médecine Québec QC Canada G1V0A6
| |
Collapse
|
17
|
Cinget B, Gérardi S, Beaulieu J, Bousquet J. Less pollen-mediated gene flow for more signatures of glacial lineages: congruent evidence from balsam fir cpDNA and mtDNA for multiple refugia in eastern and central North America. PLoS One 2015; 10:e0122815. [PMID: 25849816 PMCID: PMC4388536 DOI: 10.1371/journal.pone.0122815] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2014] [Accepted: 02/23/2015] [Indexed: 11/18/2022] Open
Abstract
The phylogeographic structure and postglacial history of balsam fir (Abies balsamea), a transcontinental North American boreal conifer, was inferred using mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) markers. Genetic structure among 107 populations (mtDNA data) and 75 populations (cpDNA data) was analyzed using Bayesian and genetic distance approaches. Population differentiation was high for mtDNA (dispersed by seeds only), but also for cpDNA (dispersed by seeds and pollen), indicating that pollen gene flow is more restricted in balsam fir than in other boreal conifers. Low cpDNA gene flow in balsam fir may relate to low pollen production due to the inherent biology of the species and populations being decimated by recurrent spruce budworm epidemics, and/or to low dispersal of pollen grains due to their peculiar structural properties. Accordingly, a phylogeographic structure was detected using both mtDNA and cpDNA markers and population structure analyses supported the existence of at least five genetically distinct glacial lineages in central and eastern North America. Four of these would originate from glacial refugia located south of the Laurentide ice sheet, while the last one would have persisted in the northern Labrador region. As expected due to reduced pollen-mediated gene flow, congruence between the geographic distribution of mtDNA and cpDNA lineages was higher than in other North American conifers. However, concordance was not complete, reflecting that restricted but nonetheless detectable cpDNA gene flow among glacial lineages occurred during the Holocene. As a result, new cpDNA and mtDNA genome combinations indicative of cytoplasmic genome capture were observed.
Collapse
Affiliation(s)
- Benjamin Cinget
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for Systems and Integrative Biology, Wood and Forest Sciences, Laval University, Québec, Québec, Canada
| | - Sébastien Gérardi
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for Systems and Integrative Biology, Wood and Forest Sciences, Laval University, Québec, Québec, Canada
| | - Jean Beaulieu
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for Systems and Integrative Biology, Wood and Forest Sciences, Laval University, Québec, Québec, Canada
- Natural Resources Canada, Canadian Wood Fibre Centre, Canadian Wood Fibre Centre, Natural Resources Canada, Québec, Québec, Canada
| | - Jean Bousquet
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for Systems and Integrative Biology, Wood and Forest Sciences, Laval University, Québec, Québec, Canada
- * E-mail:
| |
Collapse
|
18
|
van Loo M, Hintsteiner W, Pötzelsberger E, Schüler S, Hasenauer H. Intervarietal and intravarietal genetic structure in Douglas-fir: nuclear SSRs bring novel insights into past population demographic processes, phylogeography, and intervarietal hybridization. Ecol Evol 2015; 5:1802-17. [PMID: 26140197 PMCID: PMC4485962 DOI: 10.1002/ece3.1435] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2015] [Revised: 10/23/2014] [Accepted: 01/21/2015] [Indexed: 01/22/2023] Open
Abstract
Douglas-fir (Pseudotsuga menziesii) is one of numerous wide-range forest tree species represented by subspecies/varieties, which hybridize in contact zones. This study examined the genetic structure of this North American conifer and its two hybridizing varieties, coastal and Rocky Mountain, at intervarietal and intravarietal level. The genetic structure was subsequently associated with the Pleistocene refugial history, postglacial migration and intervarietal hybridization/introgression. Thirty-eight populations from the USA and Canada were genotyped for 13 nuclear SSRs and analyzed with simulations and traditional population genetic structuring methods. Eight genetic clusters were identified. The coastal clusters embodied five refugial populations originating from five distinct refugia. Four coastal refugial populations, three from California and one from western Canada, diverged during the Pleistocene (56.9–40.1 ka). The three Rocky Mountain clusters reflected distinct refugial populations of three glacial refugia. For Canada, ice covered during the Last Glacial Maximum, we present the following three findings. (1) One refugial population of each variety was revealed in the north of the distribution range. Additional research including paleodata is required to support and determine whether both northern populations originated from cryptic refugia situated south or north of the ice-covered area. (2) An interplay between intravarietal gene flow of different refugial populations and intervarietal gene flow by hybridization and introgression was identified. (3) The Canadian hybrid zone displayed predominantly introgressants of the Rocky Mountain into the coastal variety. This study provides new insights into the complex Quaternary dynamics of this conifer essential for understanding its evolution (outside and inside the native range), adaptation to future climates and for forest management.
Collapse
Affiliation(s)
- Marcela van Loo
- Institute of Silviculture, University of Natural Resources and Life Sciences Peter Jordan Straße 82, 1190, Wien, Austria
| | - Wolfgang Hintsteiner
- Institute of Silviculture, University of Natural Resources and Life Sciences Peter Jordan Straße 82, 1190, Wien, Austria ; alpS-GmbH Grabenweg 68, 6020, Innsbruck, Austria
| | - Elisabeth Pötzelsberger
- Institute of Silviculture, University of Natural Resources and Life Sciences Peter Jordan Straße 82, 1190, Wien, Austria
| | - Silvio Schüler
- Department of Forest Genetics, Federal Research and Training Centre for Forests, Natural Hazards and Landscapes Hauptstr. 7, 1140, Vienna, Austria
| | - Hubert Hasenauer
- Institute of Silviculture, University of Natural Resources and Life Sciences Peter Jordan Straße 82, 1190, Wien, Austria
| |
Collapse
|
19
|
When east meets west: population structure of a high-latitude resident species, the boreal chickadee (Poecile hudsonicus). Heredity (Edinb) 2013; 111:321-9. [PMID: 23759728 DOI: 10.1038/hdy.2013.54] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2012] [Revised: 04/16/2013] [Accepted: 04/30/2013] [Indexed: 11/09/2022] Open
Abstract
The population genetic structure of northern boreal species has been strongly influenced both by the Quaternary glaciations and the presence of contemporary barriers, such as mountain ranges and rivers. We used a combination of mitochondrial DNA (mtDNA), nuclear microsatellites and spatial distribution modelling to study the population genetic structure of the boreal chickadee (Poecile hudsonicus), a resident passerine, and to investigate whether historical or contemporary barriers have influenced this northern species. MtDNA data showed evidence of eastern and western groups, with secondary admixture occurring in central Canada. This suggests that the boreal chickadee probably persisted in multiple glacial refugia, one in Beringia and at least one in the east. Palaeo-distribution modelling identified suitable habitat in Beringia (Alaska), Atlantic Canada and the southern United States, and correspond to divergence dates of 60-96 kya. Pairwise FST values for both mtDNA and microsatellites were significant for all comparisons involving Newfoundland, though mtDNA data suggest a more recent separation. Furthermore, unlike mtDNA data, nuclear data support population connectivity among the continental populations, possibly due to male-biased dispersal. Although both are significant, the isolation-by-distance signal is much stronger for mtDNA (r(2)=0.51) than for microsatellites (r(2)=0.05), supporting the hypothesis of male-biased dispersal. The population structure of the boreal chickadee was influenced by isolation in multiple refugia and contemporary barriers. In addition to geographical distance, physical barriers such as the Strait of Belle Isle and northern mountains in Alaska are restricting gene flow, whereas the Rocky Mountains in the west are a porous barrier.
Collapse
|
20
|
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. [PMID: 23724860 PMCID: PMC3674900 DOI: 10.1186/1471-2164-14-368] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2013] [Accepted: 05/24/2013] [Indexed: 12/27/2022] Open
Abstract
BACKGROUND The genomic architecture of adaptive traits remains poorly understood in non-model plants. Various approaches can be used to bridge this gap, including the mapping of quantitative trait loci (QTL) in pedigrees, and genetic association studies in non-structured populations. Here we present results on the genomic architecture of adaptive traits in black spruce, which is a widely distributed conifer of the North American boreal forest. As an alternative to the usual candidate gene approach, a candidate SNP approach was developed for association testing. RESULTS A genetic map containing 231 gene loci was used to identify QTL that were related to budset timing and to tree height assessed over multiple years and sites. Twenty-two unique genomic regions were identified, including 20 that were related to budset timing and 6 that were related to tree height. From results of outlier detection and bulk segregant analysis for adaptive traits using DNA pool sequencing of 434 genes, 52 candidate SNPs were identified and subsequently tested in genetic association studies for budset timing and tree height assessed over multiple years and sites. A total of 34 (65%) SNPs were significantly associated with budset timing, or tree height, or both. Although the percentages of explained variance (PVE) by individual SNPs were small, several significant SNPs were shared between sites and among years. CONCLUSIONS The sharing of genomic regions and significant SNPs between budset timing and tree height indicates pleiotropic effects. Significant QTLs and SNPs differed quite greatly among years, suggesting that different sets of genes for the same characters are involved at different stages in the tree's life history. The functional diversity of genes carrying significant SNPs and low observed PVE further indicated that a large number of polymorphisms are involved in adaptive genetic variation. Accordingly, for undomesticated species such as black spruce with natural populations of large effective size and low linkage disequilibrium, efficient marker systems that are predictive of adaptation should require the survey of large numbers of SNPs. Candidate SNP approaches like the one developed in the present study could contribute to reducing these numbers.
Collapse
Affiliation(s)
- Julien Prunier
- Centre for Forest Research, and Institute for Systems and Integrative Biology, Université Laval, Québec, Québec G1V 0A6, Canada.
| | | | | | | | | | | | | |
Collapse
|
21
|
Prunier J, Gérardi 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. [PMID: 22805595 DOI: 10.1111/j.1365-294x.2012.05691.x] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
In response to selective pressure, adaptation may follow different genetic pathways throughout the natural range of a species due to historical differentiation in standing genetic variation. Using 41 populations of black spruce (Picea mariana), the objectives of this study were to identify adaptive genetic polymorphisms related to temperature and precipitation variation across the transcontinental range of the species, and to evaluate the potential influence of historical events on their geographic distribution. Population structure was first inferred using 50 control nuclear markers. Then, 47 candidate gene SNPs identified in previous genome scans were tested for relationship with climatic factors using an F(ST) -based outlier method and regressions between allele frequencies and climatic variations. Two main intraspecific lineages related to glacial vicariance were detected at the transcontinental scale. Within-lineage analyses of allele frequencies allowed the identification of 23 candidate SNPs significantly related to precipitation and/or temperature variation, among which seven were common to both lineages, eight were specific to the eastern lineage and eight were specific to the western lineage. The implication of these candidate SNPs in adaptive processes was further supported by gene functional annotations. Multiple evidences indicated that the occurrence of lineage-specific adaptive SNPs was better explained by selection acting on historically differentiated gene pools rather than differential selection due to heterogeneity of interacting environmental factors and pleiotropic effects. Taken together, these findings suggest that standing genetic variation of potentially adaptive nature has been modified by historical events, hence affecting the outcome of recent selection and leading to different adaptive routes between intraspecific lineages.
Collapse
Affiliation(s)
- Julien Prunier
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research, Université Laval, Québec, Québec, Canada G1V 0A6
| | | | | | | | | |
Collapse
|
22
|
Godbout J, Yeh FC, Bousquet J. Large-scale asymmetric introgression of cytoplasmic DNA reveals Holocene range displacement in a North American boreal pine complex. Ecol Evol 2012; 2:1853-66. [PMID: 22957188 PMCID: PMC3433990 DOI: 10.1002/ece3.294] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2012] [Revised: 05/08/2012] [Accepted: 05/09/2012] [Indexed: 11/11/2022] Open
Abstract
Jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta var. latifolia) are two North American boreal hard pines that hybridize in their zone of contact in western Canada. The main objective of this study was to characterize their patterns of introgression resulting from past and recent gene flow, using cytoplasmic markers having maternal or paternal inheritance. Mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) diversity was assessed in allopatric populations of each species and in stands from the current zone of contact containing morphological hybrids. Cluster analyses were used to identify genetic discontinuities among groups of populations. A canonical analysis was also conducted to detect putative associations among cytoplasmic DNA variation, tree morphology, and site ecological features. MtDNA introgression was extensive and asymmetric: it was detected in P. banksiana populations from the hybrid zone and from allopatric areas, but not in P. contorta populations. Very weak cpDNA introgression was observed, and only in P. banksiana populations. The mtDNA introgression pattern indicated that central Canada was first colonized by migrants from a P. contorta glacial population located west of the Rocky Mountains, before being replaced by P. banksiana migrating westward during the Holocene. In contrast, extensive pollen gene flow would have erased the cpDNA traces of this ancient presence of P. contorta. Additional evidence for this process was provided by the results of canonical analysis, which indicated that the current cpDNA background of trees reflected recent pollen gene flow from the surrounding dominant species rather than historical events that took place during the postglacial colonization.
Collapse
Affiliation(s)
- Julie Godbout
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for Systems and Integrative Biology, Université Laval1030 avenue de la Médecine, Québec, Québec, Canada, G1V 0A6
| | - Francis C Yeh
- Department of Renewable Resources, Faculty of Agricultural, Life and Environmental Sciences, University of AlbertaEdmonton, Alberta, Canada, T6G 2P5
| | - Jean Bousquet
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for Systems and Integrative Biology, Université Laval1030 avenue de la Médecine, Québec, Québec, Canada, G1V 0A6
| |
Collapse
|
23
|
Prunier J, Laroche J, Beaulieu J, Bousquet J. Scanning the genome for gene SNPs related to climate adaptation and estimating selection at the molecular level in boreal black spruce. Mol Ecol 2011; 20:1702-16. [PMID: 21375634 DOI: 10.1111/j.1365-294x.2011.05045.x] [Citation(s) in RCA: 91] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Outlier detection methods were used to scan the genome of the boreal conifer black spruce (Picea mariana [Mill.] B.S.P.) for gene single-nucleotide polymorphisms (SNPs) potentially involved in adaptations to temperature and precipitation variations. The scan involved 583 SNPs from 313 genes potentially playing adaptive roles. Differentiation estimates among population groups defined following variation in temperature and precipitation were moderately high for adaptive quantitative characters such as the timing of budset or tree height (Q(ST) = 0.189-0.314). Average differentiation estimates for gene SNPs were null, with F(ST) values of 0.005 and 0.006, respectively, among temperature and precipitation population groups. Using two detection approaches, a total of 26 SNPs from 25 genes distributed among 11 of the 12 linkage groups of black spruce were detected as outliers with F(ST) as high as 0.078. Nearly half of the outlier SNPs were located in exons and half of those were nonsynonymous. The functional annotations of genes carrying outlier SNPs and regression analyses between the frequencies of these SNPs and climatic variables supported their implication in adaptive processes. Several genes carrying outlier SNPs belonged to gene families previously found to harbour outlier SNPs in a reproductively isolated but largely sympatric congeneric species, suggesting differential subfunctionalization of gene duplicates. Selection coefficient estimates (S) were moderate but well above the magnitude of drift (>>1/N(e)), indicating that the signature of natural selection could be detected at the nucleotide level despite the recent establishment of these populations during the Holocene.
Collapse
Affiliation(s)
- Julien Prunier
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research 1030 avenue de Médecine, Université Laval, QC G1V 0A6, Canada
| | | | | | | |
Collapse
|
24
|
Godbout J, Beaulieu J, Bousquet J. Phylogeographic structure of jack pine (Pinus banksiana; Pinaceae) supports the existence of a coastal glacial refugium in northeastern North America. AMERICAN JOURNAL OF BOTANY 2010; 97:1903-1912. [PMID: 21616829 DOI: 10.3732/ajb.1000148] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
PREMISE OF THE STUDY The genetic structure of jack pine (Pinus banksiana Lamb.), a North American boreal conifer with large longitudinal distribution, was investigated to test for the possible existence of a genetically distinct lineage in the Maritimes region in northeastern North America, which could be indicative of a mid-latitude coastal refuge during the last glaciation. • METHODS One maternally inherited mitochondrial DNA (mtDNA) minisatellite marker and four paternally inherited chloroplast DNA (cpDNA) microsatellite markers were used to assess the range-wide geographical structure of jack pine populations with particular focus on northeastern North America. • KEY RESULTS The populations from the Maritimes region presented a unique mtDNA background characterized by very low diversity and the preponderance of a distinctive mitotype. The distribution of cpDNA diversity was not spatially structured, though three chlorotypes were restricted to the east. • CONCLUSIONS MtDNA data suggest that populations from the Maritimes region derive from a genetically depauperated north-coastal refugium. Contrastingly, the much higher geographical uniformity observed for cpDNA variation indicates that gene flow by pollen had been much more effective than seed gene flow at homogenizing population structure.
Collapse
Affiliation(s)
- Julie Godbout
- Canada Research Chair in Forest and Environmental Genomics, Centre for Forest Research and Institute for Systems and Integrative Biology, 1030 Avenue de la Médecine, Université Laval, Quebec, Quebec G1V 0A6 Canada
| | | | | |
Collapse
|