1
|
Neale DB, Zimin AV, Meltzer A, Bhattarai A, Amee M, Figueroa Corona L, Allen BJ, Puiu D, Wright J, De La Torre AR, McGuire PE, Timp W, Salzberg SL, Wegrzyn JL. A genome sequence for the threatened whitebark pine. G3 (Bethesda) 2024; 14:jkae061. [PMID: 38526344 PMCID: PMC11075562 DOI: 10.1093/g3journal/jkae061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 02/29/2024] [Accepted: 03/12/2024] [Indexed: 03/26/2024]
Abstract
Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.
Collapse
Affiliation(s)
- David B Neale
- Department of Plant Sciences, University of California, Davis, CA 95616, USA
- Whitebark Pine Ecosystem Foundation, Missoula, MT 59808, USA
| | - Aleksey V Zimin
- Department of Biomedical Engineering and Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Amy Meltzer
- Department of Biomedical Engineering and Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Akriti Bhattarai
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269, USA
| | - Maurice Amee
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269, USA
| | | | - Brian J Allen
- Department of Plant Sciences, University of California, Davis, CA 95616, USA
- University of California Cooperative Extension, Central Sierra, Jackson, CA 95642, USA
| | - Daniela Puiu
- Department of Biomedical Engineering and Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Jessica Wright
- USDA Forest Service, Pacific Southwest Research Station, Davis, CA 95618, USA
| | | | - Patrick E McGuire
- Department of Plant Sciences, University of California, Davis, CA 95616, USA
| | - Winston Timp
- Department of Biomedical Engineering and Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Steven L Salzberg
- Department of Biomedical Engineering and Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21218, USA
- Departments of Computer Science and Biostatistics, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Jill L Wegrzyn
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269, USA
- Institute for Systems Genomics, University of Connecticut, Storrs, CT 06269, USA
| |
Collapse
|
2
|
Neale DB, Zimin AV, Meltzer A, Bhattarai A, Amee M, Corona LF, Allen BJ, Puiu D, Wright J, Torre ARDL, McGuire PE, Timp W, Salzberg SL, Wegrzyn JL. A Genome Sequence for the Threatened Whitebark Pine. bioRxiv 2023:2023.11.16.567420. [PMID: 38014212 PMCID: PMC10680812 DOI: 10.1101/2023.11.16.567420] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Whitebark pine (WBP, Pinus albicaulis ) is a white pine of subalpine regions in western contiguous US and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola ) and additional threats from mountain pine beetle ( Dendroctonus ponderosae ), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short-reads of haploid megametophyte tissue and Oxford Nanopore long-reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gbp of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gbp). Approximately 87.2% (24.0 Gbp) of total sequence was placed on the twelve WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich-repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the three subclasses of NLRs (TNL, CNL, RNL). Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo assembled transcriptomes.
Collapse
|
3
|
Neale DB, Zimin AV, Zaman S, Scott AD, Shrestha B, Workman RE, Puiu D, Allen BJ, Moore ZJ, Sekhwal MK, De La Torre AR, McGuire PE, Burns E, Timp W, Wegrzyn JL, Salzberg SL. Assembled and annotated 26.5 Gbp coast redwood genome: a resource for estimating evolutionary adaptive potential and investigating hexaploid origin. G3 (Bethesda) 2022; 12:6460957. [PMID: 35100403 PMCID: PMC8728005 DOI: 10.1093/g3journal/jkab380] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/26/2021] [Accepted: 10/25/2021] [Indexed: 12/15/2022]
Abstract
Sequencing, assembly, and annotation of the 26.5 Gbp hexaploid genome of coast redwood (Sequoia sempervirens) was completed leading toward discovery of genes related to climate adaptation and investigation of the origin of the hexaploid genome. Deep-coverage short-read Illumina sequencing data from haploid tissue from a single seed were combined with long-read Oxford Nanopore Technologies sequencing data from diploid needle tissue to create an initial assembly, which was then scaffolded using proximity ligation data to produce a highly contiguous final assembly, SESE 2.1, with a scaffold N50 size of 44.9 Mbp. The assembly included several scaffolds that span entire chromosome arms, confirmed by the presence of telomere and centromere sequences on the ends of the scaffolds. The structural annotation produced 118,906 genes with 113 containing introns that exceed 500 Kbp in length and one reaching 2 Mb. Nearly 19 Gbp of the genome represented repetitive content with the vast majority characterized as long terminal repeats, with a 2.9:1 ratio of Copia to Gypsy elements that may aid in gene expression control. Comparison of coast redwood to other conifers revealed species-specific expansions for a plethora of abiotic and biotic stress response genes, including those involved in fungal disease resistance, detoxification, and physical injury/structural remodeling and others supporting flavonoid biosynthesis. Analysis of multiple genes that exist in triplicate in coast redwood but only once in its diploid relative, giant sequoia, supports a previous hypothesis that the hexaploidy is the result of autopolyploidy rather than any hybridizations with separate but closely related conifer species.
Collapse
Affiliation(s)
- David B Neale
- Department of Plant Sciences, University of California, Davis, Davis, CA 95616, USA
| | - Aleksey V Zimin
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA.,Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21211, USA
| | - Sumaira Zaman
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269, USA.,Department of Computer Science & Engineering, University of Connecticut, Storrs, CT 06269, USA
| | - Alison D Scott
- Department of Plant Sciences, University of California, Davis, Davis, CA 95616, USA
| | - Bikash Shrestha
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269, USA
| | - Rachael E Workman
- Department of Molecular Biology and Genetics, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Daniela Puiu
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA.,Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21211, USA
| | - Brian J Allen
- Department of Plant Sciences, University of California, Davis, Davis, CA 95616, USA
| | - Zane J Moore
- Department of Plant Sciences, University of California, Davis, Davis, CA 95616, USA
| | - Manoj K Sekhwal
- School of Forestry, Northern Arizona University, Flagstaff, AZ 86011, USA
| | | | - Patrick E McGuire
- Department of Plant Sciences, University of California, Davis, Davis, CA 95616, USA
| | - Emily Burns
- Save the Redwoods League, San Francisco, CA 94104, USA
| | - Winston Timp
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA.,Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21211, USA.,Department of Molecular Biology and Genetics, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Jill L Wegrzyn
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269, USA.,Institute for Systems Genomics, University of Connecticut, Storrs, CT 06269, USA
| | - Steven L Salzberg
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA.,Center for Computational Biology, Johns Hopkins University, Baltimore, MD 21211, USA.,Department of Computer Science, Johns Hopkins University, Baltimore, MD 21218, USA.,Department of Biostatistics, Johns Hopkins University, Baltimore, MD 21205, USA
| |
Collapse
|
4
|
De La Torre AR, Sekhwal MK, Puiu D, Salzberg SL, Scott AD, Allen B, Neale DB, Chin ARO, Buckley TN. Genome-wide association identifies candidate genes for drought tolerance in coast redwood and giant sequoia. Plant J 2022; 109:7-22. [PMID: 34800071 PMCID: PMC10773529 DOI: 10.1111/tpj.15592] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Revised: 11/05/2021] [Accepted: 11/16/2021] [Indexed: 06/13/2023]
Abstract
Drought is a major limitation for survival and growth in plants. With more frequent and severe drought episodes occurring due to climate change, it is imperative to understand the genomic and physiological basis of drought tolerance to be able to predict how species will respond in the future. In this study, univariate and multitrait multivariate genome-wide association study methods were used to identify candidate genes in two iconic and ecosystem-dominating species of the western USA, coast redwood and giant sequoia, using 10 drought-related physiological and anatomical traits and genome-wide sequence-capture single nucleotide polymorphisms. Population-level phenotypic variation was found in carbon isotope discrimination, osmotic pressure at full turgor, xylem hydraulic diameter, and total area of transporting fibers in both species. Our study identified new 78 new marker × trait associations in coast redwood and six in giant sequoia, with genes involved in a range of metabolic, stress, and signaling pathways, among other functions. This study contributes to a better understanding of the genomic basis of drought tolerance in long-generation conifers and helps guide current and future conservation efforts in the species.
Collapse
Affiliation(s)
- Amanda R. De La Torre
- School of Forestry, Northern Arizona University, 200 E. Pine Knoll, Flagstaff, AZ 86011, USA
| | - Manoj K. Sekhwal
- School of Forestry, Northern Arizona University, 200 E. Pine Knoll, Flagstaff, AZ 86011, USA
| | - Daniela Puiu
- Department of Biomedical Engineering, Computer Science and Biostatistics & Center for Computational Biology, John Hopkins University, 3100 Wyman Park Dr, Wyman Park Building, Room S220, Baltimore, MD 21211, USA
| | - Steven L. Salzberg
- Department of Biomedical Engineering, Computer Science and Biostatistics & Center for Computational Biology, John Hopkins University, 3100 Wyman Park Dr, Wyman Park Building, Room S220, Baltimore, MD 21211, USA
| | - Alison D. Scott
- Department of Plant Sciences, University of California, Davis, One Shields Avenue, Davis, CA 95616, USA
| | - Brian Allen
- Department of Plant Sciences, University of California, Davis, One Shields Avenue, Davis, CA 95616, USA
| | - David B. Neale
- Department of Plant Sciences, University of California, Davis, One Shields Avenue, Davis, CA 95616, USA
| | - Alana R. O. Chin
- Department of Plant Sciences, University of California, Davis, One Shields Avenue, Davis, CA 95616, USA
| | - Thomas N. Buckley
- Department of Plant Sciences, University of California, Davis, One Shields Avenue, Davis, CA 95616, USA
| |
Collapse
|
5
|
De La Torre AR, Sekhwal MK, Neale DB. Selective Sweeps and Polygenic Adaptation Drive Local Adaptation along Moisture and Temperature Gradients in Natural Populations of Coast Redwood and Giant Sequoia. Genes (Basel) 2021; 12:1826. [PMID: 34828432 PMCID: PMC8621000 DOI: 10.3390/genes12111826] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2021] [Revised: 11/18/2021] [Accepted: 11/18/2021] [Indexed: 12/26/2022] Open
Abstract
Dissecting the genomic basis of local adaptation is a major goal in evolutionary biology and conservation science. Rapid changes in the climate pose significant challenges to the survival of natural populations, and the genomic basis of long-generation plant species is still poorly understood. Here, we investigated genome-wide climate adaptation in giant sequoia and coast redwood, two iconic and ecologically important tree species. We used a combination of univariate and multivariate genotype-environment association methods and a selective sweep analysis using non-overlapping sliding windows. We identified genomic regions of potential adaptive importance, showing strong associations to moisture variables and mean annual temperature. Our results found a complex architecture of climate adaptation in the species, with genomic regions showing signatures of selective sweeps, polygenic adaptation, or a combination of both, suggesting recent or ongoing climate adaptation along moisture and temperature gradients in giant sequoia and coast redwood. The results of this study provide a first step toward identifying genomic regions of adaptive significance in the species and will provide information to guide management and conservation strategies that seek to maximize adaptive potential in the face of climate change.
Collapse
Affiliation(s)
- Amanda R. De La Torre
- School of Forestry, Northern Arizona University, 200 E. Pine Knoll, Flagstaff, AZ 86011, USA;
| | - Manoj K. Sekhwal
- School of Forestry, Northern Arizona University, 200 E. Pine Knoll, Flagstaff, AZ 86011, USA;
| | - David B. Neale
- Department of Plant Sciences, University of California-Davis, One Shields Avenue, Davis, CA 95616, USA
| |
Collapse
|
6
|
De La Torre AR, Wilhite B, Puiu D, St. Clair JB, Crepeau MW, Salzberg SL, Langley CH, Allen B, Neale DB. Dissecting the Polygenic Basis of Cold Adaptation Using Genome-Wide Association of Traits and Environmental Data in Douglas-fir. Genes (Basel) 2021; 12:110. [PMID: 33477542 PMCID: PMC7831106 DOI: 10.3390/genes12010110] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2020] [Revised: 01/14/2021] [Accepted: 01/15/2021] [Indexed: 02/06/2023] Open
Abstract
Understanding the genomic and environmental basis of cold adaptation is key to understand how plants survive and adapt to different environmental conditions across their natural range. Univariate and multivariate genome-wide association (GWAS) and genotype-environment association (GEA) analyses were used to test associations among genome-wide SNPs obtained from whole-genome resequencing, measures of growth, phenology, emergence, cold hardiness, and range-wide environmental variation in coastal Douglas-fir (Pseudotsuga menziesii). Results suggest a complex genomic architecture of cold adaptation, in which traits are either highly polygenic or controlled by both large and small effect genes. Newly discovered associations for cold adaptation in Douglas-fir included 130 genes involved in many important biological functions such as primary and secondary metabolism, growth and reproductive development, transcription regulation, stress and signaling, and DNA processes. These genes were related to growth, phenology and cold hardiness and strongly depend on variation in environmental variables such degree days below 0c, precipitation, elevation and distance from the coast. This study is a step forward in our understanding of the complex interconnection between environment and genomics and their role in cold-associated trait variation in boreal tree species, providing a baseline for the species' predictions under climate change.
Collapse
Affiliation(s)
- Amanda R. De La Torre
- School of Forestry, Northern Arizona University, 200 E. Pine Knoll, Flagstaff, AZ 86011, USA;
| | - Benjamin Wilhite
- School of Forestry, Northern Arizona University, 200 E. Pine Knoll, Flagstaff, AZ 86011, USA;
| | - Daniela Puiu
- Center for Computational Biology, Department of Biomedical Engineering, Computer Science and Biostatistics, John Hopkins University, 3100 Wyman Park Dr, Wyman Park Building, Room S220, Baltimore, MD 21211, USA; (D.P.); (S.L.S.)
| | - John Bradley St. Clair
- USDA Forest Service, Pacific Northwest Research Station, 3200 SW Jefferson Way, Corvallis, OR 97331, USA;
| | - Marc W. Crepeau
- Department of Evolution and Ecology, University of California-Davis, One Shields Avenue, Davis, CA 95616, USA; (M.W.C.); (C.H.L.)
| | - Steven L. Salzberg
- Center for Computational Biology, Department of Biomedical Engineering, Computer Science and Biostatistics, John Hopkins University, 3100 Wyman Park Dr, Wyman Park Building, Room S220, Baltimore, MD 21211, USA; (D.P.); (S.L.S.)
| | - Charles H. Langley
- Department of Evolution and Ecology, University of California-Davis, One Shields Avenue, Davis, CA 95616, USA; (M.W.C.); (C.H.L.)
| | - Brian Allen
- Department of Plant Sciences, University of California-Davis, One Shields Avenue, Davis, CA 95616, USA; (B.A.); (D.B.N.)
| | - David B. Neale
- Department of Plant Sciences, University of California-Davis, One Shields Avenue, Davis, CA 95616, USA; (B.A.); (D.B.N.)
| |
Collapse
|
7
|
Weiss M, Sniezko RA, Puiu D, Crepeau MW, Stevens K, Salzberg SL, Langley CH, Neale DB, De La Torre AR. Genomic basis of white pine blister rust quantitative disease resistance and its relationship with qualitative resistance. Plant J 2020; 104:365-376. [PMID: 32654344 PMCID: PMC10773528 DOI: 10.1111/tpj.14928] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Revised: 06/17/2020] [Accepted: 07/01/2020] [Indexed: 06/11/2023]
Abstract
The genomic architecture and molecular mechanisms controlling variation in quantitative disease resistance loci are not well understood in plant species and have been barely studied in long-generation trees. Quantitative trait loci mapping and genome-wide association studies were combined to test a large single nucleotide polymorphism (SNP) set for association with quantitative and qualitative white pine blister rust resistance in sugar pine. In the absence of a chromosome-scale reference genome, a high-density consensus linkage map was generated to obtain locations for associated SNPs. Newly discovered associations for white pine blister rust quantitative disease resistance included 453 SNPs involved in wide biological functions, including genes associated with disease resistance and others involved in morphological and developmental processes. In addition, NBS-LRR pathogen recognition genes were found to be involved in quantitative disease resistance, suggesting these newly reported genes are qualitative genes with partial resistance, they are the result of defeated qualitative resistance due to avirulent races, or they have epistatic effects on qualitative disease resistance genes. This study is a step forward in our understanding of the complex genomic architecture of quantitative disease resistance in long-generation trees, and constitutes the first step towards marker-assisted disease resistance breeding in white pine species.
Collapse
Affiliation(s)
- Matthew Weiss
- School of Forestry, Northern Arizona University, 200 E.
Pine Knoll, Flagstaff, AZ 86011
| | - Richard A. Sniezko
- Dorena Genetic Resource Center, USDA Forest Service,
Cottage-Grove, OR 97424
| | - Daniela Puiu
- Department of Biomedical Engineering, Computer Science and
Biostatistics and Center for Computational Biology, Johns Hopkins University, 3100
Wyman Park Dr., Wyman Park Building Room S220, Baltimore, MD 21211
| | - Marc W. Crepeau
- Department of Evolution and Ecology, University of
California-Davis, One Shields Avenue, Davis, CA 95616
| | - Kristian Stevens
- Department of Evolution and Ecology, University of
California-Davis, One Shields Avenue, Davis, CA 95616
| | - Steven L. Salzberg
- Department of Biomedical Engineering, Computer Science and
Biostatistics and Center for Computational Biology, Johns Hopkins University, 3100
Wyman Park Dr., Wyman Park Building Room S220, Baltimore, MD 21211
- Departments of Computer Science and Biostatistics, Johns
Hopkins University, Baltimore, MD 21218
| | - Charles H. Langley
- Department of Evolution and Ecology, University of
California-Davis, One Shields Avenue, Davis, CA 95616
| | - David B. Neale
- Department of Plant Sciences, University of
California-Davis, One Shields Avenue, Davis, CA 95616, USA
| | - Amanda R. De La Torre
- School of Forestry, Northern Arizona University, 200 E.
Pine Knoll, Flagstaff, AZ 86011
| |
Collapse
|
8
|
De La Torre AR, Wilhite B, Neale DB. Environmental Genome-Wide Association Reveals Climate Adaptation Is Shaped by Subtle to Moderate Allele Frequency Shifts in Loblolly Pine. Genome Biol Evol 2020; 11:2976-2989. [PMID: 31599932 PMCID: PMC6821164 DOI: 10.1093/gbe/evz220] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/07/2019] [Indexed: 01/21/2023] Open
Abstract
Understanding the genomic basis of local adaptation is crucial to determine the potential of long-lived woody species to withstand changes in their natural environment. In the past, efforts to dissect the genomic architecture in gymnosperms species have been limited due to the absence of reference genomes. Recently, the genomes of some commercially important conifers, such as loblolly pine, have become available, allowing whole-genome studies of these species. In this study, we test for associations between 87k SNPs, obtained from whole-genome resequencing of loblolly pine individuals, and 270 environmental variables and combinations of them. We determine the geographic location of significant loci and identify their genomic location using our newly constructed ultradense 26k SNP linkage map. We found that water availability is the main climatic variable shaping local adaptation of the species, and found 821 SNPs showing significant associations with climatic variables or combinations of them based on the consistent results of three different genotype–environment association methods. Our results suggest that adaptation to climate in the species might have occurred by many changes in the frequency of alleles with moderate to small effect sizes, and by the smaller contribution of large effect alleles in genes related to moisture deficit, temperature and precipitation. Genomic regions of low recombination and high population differentiation harbored SNPs associated with groups of environmental variables, suggesting climate adaptation might have evolved as a result of different selection pressures acting on groups of genes associated with an aspect of climate rather than on individual environmental variables.
Collapse
Affiliation(s)
| | | | - David B Neale
- Department of Plant Sciences, University of California-Davis
| |
Collapse
|
9
|
Sun YQ, Zhao W, Xu CQ, Xu Y, El-Kassaby YA, De La Torre AR, Mao JF. Genetic Variation Related to High Elevation Adaptation Revealed by Common Garden Experiments in Pinus yunnanensis. Front Genet 2020; 10:1405. [PMID: 32117429 PMCID: PMC7027398 DOI: 10.3389/fgene.2019.01405] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Accepted: 12/23/2019] [Indexed: 12/30/2022] Open
Abstract
Local adaptation, adaptation to specialized niches and environmental clines have been extensively reported for forest trees. Investigation of the adaptive genetic variation is crucial for forest resource management and breeding, especially in the context of global climate change. Here, we utilized a Pinus yunnanensis common garden experiments established at high and low elevation sites to assess the differences in growth and survival among populations and between the two common garden sites. The studied traits showed significant variation between the two test sites and among populations, suggesting adaptive divergence. To detect genetic variation related to environment, we captured 103,608 high quality SNPs based on RNA sequencing, and used them to assess the genetic diversity and population structure. We identified 321 outlier SNPs from 131 genes showing significant divergence in allelic frequency between survival populations of two sites. Functional categories associated with adaptation to high elevation were found to be related to flavonoid biosynthesis, response to UV, DNA repair, response to reactive oxygen species, and membrane lipid metabolic process. Further investigation of the outlier genes showed overrepresentation of the flavonoid biosynthesis pathway, suggesting that this pathway may play a key role in P. yunnanensis adaptation to high elevation environments. The outlier genes identified, and their variants, provide a basic reference for advanced investigations.
Collapse
Affiliation(s)
- Yan-Qiang Sun
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, National Engineering Laboratory for Tree Breeding, Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| | - Wei Zhao
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, National Engineering Laboratory for Tree Breeding, Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| | - Chao-Qun Xu
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, National Engineering Laboratory for Tree Breeding, Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| | - Yulan Xu
- Key Laboratory for Forest Resources Conservation and Utilization in the Southwest Mountains of China, Southwest Forestry University, Kunming, China
| | - Yousry A. El-Kassaby
- Department of Forest and Conservation Sciences, Faculty of Forestry, The University of British Columbia, Vancouver, BC, Canada
| | | | - Jian-Feng Mao
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, National Engineering Laboratory for Tree Breeding, Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| |
Collapse
|
10
|
De La Torre AR, Piot A, Liu B, Wilhite B, Weiss M, Porth I. Functional and morphological evolution in gymnosperms: A portrait of implicated gene families. Evol Appl 2020; 13:210-227. [PMID: 31892953 PMCID: PMC6935586 DOI: 10.1111/eva.12839] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Revised: 04/25/2019] [Accepted: 07/02/2019] [Indexed: 12/12/2022] Open
Abstract
Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.
Collapse
Affiliation(s)
| | - Anthony Piot
- Department of Wood and Forest SciencesLaval UniversityQuebec CityQuebecCanada
- Institute for System and Integrated Biology (IBIS)Laval UniversityQuebec CityQuebecCanada
- Centre for Forest Research (CEF)Laval UniversityQuebec CityQuebecCanada
| | - Bobin Liu
- School of ForestryNorthern Arizona UniversityFlagstaffAZUSA
- College of ForestryFujian Agricultural and Forestry UniversityFuzhouFujianChina
| | | | - Matthew Weiss
- School of ForestryNorthern Arizona UniversityFlagstaffAZUSA
| | - Ilga Porth
- Department of Wood and Forest SciencesLaval UniversityQuebec CityQuebecCanada
- Institute for System and Integrated Biology (IBIS)Laval UniversityQuebec CityQuebecCanada
- Centre for Forest Research (CEF)Laval UniversityQuebec CityQuebecCanada
| |
Collapse
|
11
|
De La Torre AR, Puiu D, Crepeau MW, Stevens K, Salzberg SL, Langley CH, Neale DB. Genomic architecture of complex traits in loblolly pine. New Phytol 2019; 221:1789-1801. [PMID: 30318590 DOI: 10.1111/nph.15535] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/25/2018] [Accepted: 10/06/2018] [Indexed: 05/02/2023]
Abstract
Dissecting the genetic and genomic architecture of complex traits is essential to understand the forces maintaining the variation in phenotypic traits of ecological and economical importance. Whole-genome resequencing data were used to generate high-resolution polymorphic single nucleotide polymorphism (SNP) markers and genotype individuals from common gardens across the loblolly pine (Pinus taeda) natural range. Genome-wide associations were tested with a large phenotypic dataset comprising 409 variables including morphological traits (height, diameter, carbon isotope discrimination, pitch canker resistance), and molecular traits such as metabolites and expression of xylem development genes. Our study identified 2335 new SNP × trait associations for the species, with many SNPs located in physical clusters in the genome of the species; and the genomic location of hotspots for metabolic × genotype associations. We found a highly polygenic basis of quantitative inheritance, with significant differences in number, effects size, genomic location and frequency of alleles contributing to variation in phenotypes in the different traits. While mutation-selection balance might be shaping the genetic variation in metabolic traits, balancing selection is more likely to shape the variation in expression of xylem development genes. Our work contributes to the study of complex traits in nonmodel plant species by identifying associations at a whole-genome level.
Collapse
Affiliation(s)
- Amanda R De La Torre
- School of Forestry, Northern Arizona University, 200 E. Pine Knoll Drive, Flagstaff, AZ, 86011, USA
- Department of Plant Sciences, University of California-Davis, One Shields Avenue, Davis, CA, 95616, USA
| | - Daniela Puiu
- Center for Computational Biology, Johns Hopkins University, 1900 E. Monument St, Baltimore, MD, 21205, USA
| | - Marc W Crepeau
- Department of Evolution and Ecology, University of California-Davis, One Shields Avenue, Davis, CA, 95616, USA
| | - Kristian Stevens
- Department of Evolution and Ecology, University of California-Davis, One Shields Avenue, Davis, CA, 95616, USA
| | - Steven L Salzberg
- Center for Computational Biology, Johns Hopkins University, 1900 E. Monument St, Baltimore, MD, 21205, USA
- Department of Biomedical Engineering, Computer Science and Biostatistics, Johns Hopkins University, Baltimore, MD, 21205, USA
| | - Charles H Langley
- Department of Evolution and Ecology, University of California-Davis, One Shields Avenue, Davis, CA, 95616, USA
| | - David B Neale
- Department of Plant Sciences, University of California-Davis, One Shields Avenue, Davis, CA, 95616, USA
| |
Collapse
|
12
|
Li Z, De La Torre AR, Sterck L, Cánovas FM, Avila C, Merino I, Cabezas JA, Cervera MT, Ingvarsson PK, Van de Peer Y. Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants. Genome Biol Evol 2017; 9:1130-1147. [PMID: 28460034 PMCID: PMC5414570 DOI: 10.1093/gbe/evx070] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/26/2017] [Indexed: 01/02/2023] Open
Abstract
Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.
Collapse
Affiliation(s)
- Zhen Li
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.,Center for Plant Systems Biology, VIB, Ghent, Belgium.,Bioinformatics Institute Ghent, Ghent, Belgium
| | - Amanda R De La Torre
- Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden.,Department of Plant Sciences, University of California-Davis, Davis, CA
| | - Lieven Sterck
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.,Center for Plant Systems Biology, VIB, Ghent, Belgium.,Bioinformatics Institute Ghent, Ghent, Belgium
| | - Francisco M Cánovas
- Departamento de Biología Molecular y Bioquímica, Facultad de Ciencias, Universidad de Málaga, Campus Universitario de Teatinos s/n, Málaga, Spain
| | - Concepción Avila
- Departamento de Biología Molecular y Bioquímica, Facultad de Ciencias, Universidad de Málaga, Campus Universitario de Teatinos s/n, Málaga, Spain
| | - Irene Merino
- Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | | | | | - Pär K Ingvarsson
- Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden.,Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Yves Van de Peer
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.,Center for Plant Systems Biology, VIB, Ghent, Belgium.,Bioinformatics Institute Ghent, Ghent, Belgium.,Genomics Research Institute, University of Pretoria, Hatfield Campus, Pretoria, South Africa
| |
Collapse
|
13
|
De La Torre AR, Li Z, Van de Peer Y, Ingvarsson PK. Contrasting Rates of Molecular Evolution and Patterns of Selection among Gymnosperms and Flowering Plants. Mol Biol Evol 2017; 34:1363-1377. [PMID: 28333233 PMCID: PMC5435085 DOI: 10.1093/molbev/msx069] [Citation(s) in RCA: 103] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The majority of variation in rates of molecular evolution among seed plants remains both unexplored and unexplained. Although some attention has been given to flowering plants, reports of molecular evolutionary rates for their sister plant clade (gymnosperms) are scarce, and to our knowledge differences in molecular evolution among seed plant clades have never been tested in a phylogenetic framework. Angiosperms and gymnosperms differ in a number of features, of which contrasting reproductive biology, life spans, and population sizes are the most prominent. The highly conserved morphology of gymnosperms evidenced by similarity of extant species to fossil records and the high levels of macrosynteny at the genomic level have led scientists to believe that gymnosperms are slow-evolving plants, although some studies have offered contradictory results. Here, we used 31,968 nucleotide sites obtained from orthologous genes across a wide taxonomic sampling that includes representatives of most conifers, cycads, ginkgo, and many angiosperms with a sequenced genome. Our results suggest that angiosperms and gymnosperms differ considerably in their rates of molecular evolution per unit time, with gymnosperm rates being, on average, seven times lower than angiosperm species. Longer generation times and larger genome sizes are some of the factors explaining the slow rates of molecular evolution found in gymnosperms. In contrast to their slow rates of molecular evolution, gymnosperms possess higher substitution rate ratios than angiosperm taxa. Finally, our study suggests stronger and more efficient purifying and diversifying selection in gymnosperm than in angiosperm species, probably in relation to larger effective population sizes.
Collapse
Affiliation(s)
- Amanda R De La Torre
- Department of Plant Sciences, University of California-Davis, Davis, CA.,Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden
| | - Zhen Li
- Department of Plant Systems Biology, VIB, Ghent, Belgium.,Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
| | - Yves Van de Peer
- Department of Plant Systems Biology, VIB, Ghent, Belgium.,Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.,Genomics Research Institute, University of Pretoria, Hatfield Campus, Pretoria, South Africa
| | - Pär K Ingvarsson
- Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden.,Department of Plant Biology, Uppsala Biocenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
14
|
Neale DB, Martínez-García PJ, De La Torre AR, Montanari S, Wei XX. Novel Insights into Tree Biology and Genome Evolution as Revealed Through Genomics. Annu Rev Plant Biol 2017; 68:457-483. [PMID: 28226237 DOI: 10.1146/annurev-arplant-042916-041049] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Reference genome sequences are the key to the discovery of genes and gene families that determine traits of interest. Recent progress in sequencing technologies has enabled a rapid increase in genome sequencing of tree species, allowing the dissection of complex characters of economic importance, such as fruit and wood quality and resistance to biotic and abiotic stresses. Although the number of reference genome sequences for trees lags behind those for other plant species, it is not too early to gain insight into the unique features that distinguish trees from nontree plants. Our review of the published data suggests that, although many gene families are conserved among herbaceous and tree species, some gene families, such as those involved in resistance to biotic and abiotic stresses and in the synthesis and transport of sugars, are often expanded in tree genomes. As the genomes of more tree species are sequenced, comparative genomics will further elucidate the complexity of tree genomes and how this relates to traits unique to trees.
Collapse
Affiliation(s)
- David B Neale
- Department of Plant Sciences, University of California, Davis, California 95616;
| | | | - Amanda R De La Torre
- Department of Plant Sciences, University of California, Davis, California 95616;
| | - Sara Montanari
- Department of Plant Sciences, University of California, Davis, California 95616;
| | - Xiao-Xin Wei
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
| |
Collapse
|
15
|
De La Torre AR, Lin YC, Van de Peer Y, Ingvarsson PK. Genome-wide analysis reveals diverged patterns of codon bias, gene expression, and rates of sequence evolution in picea gene families. Genome Biol Evol 2015; 7:1002-15. [PMID: 25747252 PMCID: PMC4419791 DOI: 10.1093/gbe/evv044] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The recent sequencing of several gymnosperm genomes has greatly facilitated studying the evolution of their genes and gene families. In this study, we examine the evidence for expression-mediated selection in the first two fully sequenced representatives of the gymnosperm plant clade (Picea abies and Picea glauca). We use genome-wide estimates of gene expression (>50,000 expressed genes) to study the relationship between gene expression, codon bias, rates of sequence divergence, protein length, and gene duplication. We found that gene expression is correlated with rates of sequence divergence and codon bias, suggesting that natural selection is acting on Picea protein-coding genes for translational efficiency. Gene expression, rates of sequence divergence, and codon bias are correlated with the size of gene families, with large multicopy gene families having, on average, a lower expression level and breadth, lower codon bias, and higher rates of sequence divergence than single-copy gene families. Tissue-specific patterns of gene expression were more common in large gene families with large gene expression divergence than in single-copy families. Recent family expansions combined with large gene expression variation in paralogs and increased rates of sequence evolution suggest that some Picea gene families are rapidly evolving to cope with biotic and abiotic stress. Our study highlights the importance of gene expression and natural selection in shaping the evolution of protein-coding genes in Picea species, and sets the ground for further studies investigating the evolution of individual gene families in gymnosperms.
Collapse
Affiliation(s)
| | - Yao-Cheng Lin
- Department of Plant Systems Biology, VIB, and Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
| | - Yves Van de Peer
- Department of Plant Systems Biology, VIB, and Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium Genomics Research Institute, University of Pretoria, South Africa
| | - Pär K Ingvarsson
- Department of Ecology and Environmental Science, Umeå University, Sweden Umeå Plant Science Centre, Umeå, Sweden
| |
Collapse
|
16
|
De La Torre AR, Birol I, Bousquet J, Ingvarsson PK, Jansson S, Jones SJM, Keeling CI, MacKay J, Nilsson O, Ritland K, Street N, Yanchuk A, Zerbe P, Bohlmann J. Insights into conifer giga-genomes. Plant Physiol 2014; 166:1724-32. [PMID: 25349325 PMCID: PMC4256843 DOI: 10.1104/pp.114.248708] [Citation(s) in RCA: 86] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Insights from sequenced genomes of major land plant lineages have advanced research in almost every aspect of plant biology. Until recently, however, assembled genome sequences of gymnosperms have been missing from this picture. Conifers of the pine family (Pinaceae) are a group of gymnosperms that dominate large parts of the world's forests. Despite their ecological and economic importance, conifers seemed long out of reach for complete genome sequencing, due in part to their enormous genome size (20-30 Gb) and the highly repetitive nature of their genomes. Technological advances in genome sequencing and assembly enabled the recent publication of three conifer genomes: white spruce (Picea glauca), Norway spruce (Picea abies), and loblolly pine (Pinus taeda). These genome sequences revealed distinctive features compared with other plant genomes and may represent a window into the past of seed plant genomes. This Update highlights recent advances, remaining challenges, and opportunities in light of the publication of the first conifer and gymnosperm genomes.
Collapse
Affiliation(s)
- Amanda R De La Torre
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Inanc Birol
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Jean Bousquet
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Pär K Ingvarsson
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Stefan Jansson
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Steven J M Jones
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Christopher I Keeling
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - John MacKay
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Ove Nilsson
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Kermit Ritland
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Nathaniel Street
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Alvin Yanchuk
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Philipp Zerbe
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| | - Jörg Bohlmann
- Department of Ecology and Environmental Sciences (A.R.D.L.T., P.K.I.) and Umeå Plant Science Center, Department of Plant Physiology (P.K.I., S.J., O.N., N.S.), Umeå University, SE-901 87 Umea, Sweden;Genome Sciences Centre, British Columbia Cancer Agency, Vancouver, British Columbia, Canada V5Z 4S6 (I.B., S.J.M.J.);Canada Research Chair in Forest and Environmental Genomics (J.Bou.) and Center for Forest Research and Institute for Systems and Integrative Biology (J.Bou., J.M.), Université Laval, Quebec, Quebec, Canada G1V 0A6;Michael Smith Laboratories (C.I.K., P.Z., J.Boh.) and Department of Forest and Conservation Sciences (K.R., J.Boh.), University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z4; andBritish Columbia Ministry of Forests, Lands, and Natural Resource Operations, Victoria, British Columbia, Canada V8W 9C2 (A.Y.)
| |
Collapse
|
17
|
De La Torre AR, Roberts DR, Aitken SN. Genome-wide admixture and ecological niche modelling reveal the maintenance of species boundaries despite long history of interspecific gene flow. Mol Ecol 2014; 23:2046-59. [PMID: 24597663 PMCID: PMC4228761 DOI: 10.1111/mec.12710] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2013] [Revised: 02/27/2014] [Accepted: 02/28/2014] [Indexed: 01/12/2023]
Abstract
The maintenance of species boundaries despite interspecific gene flow has been a continuous source of interest in evolutionary biology. Many hybridizing species have porous genomes with regions impermeable to introgression, conferring reproductive barriers between species. We used ecological niche modelling to study the glacial and postglacial recolonization patterns between the widely hybridizing spruce species Picea glauca and P. engelmannii in western North America. Genome-wide estimates of admixture based on a panel of 311 candidate gene single nucleotide polymorphisms (SNP) from 290 genes were used to assess levels of admixture and introgression and to identify loci putatively involved in adaptive differences or reproductive barriers between species. Our palaeoclimatic modelling suggests that these two closely related species have a long history of hybridization and introgression, dating to at least 21,000 years ago, yet species integrity is maintained by a combination of strong environmental selection and reduced current interspecific gene flow. Twenty loci showed evidence of divergent selection, including six loci that were both Fst outliers and associated with climatic gradients, and fourteen loci that were either outliers or showed associations with climate. These included genes responsible for carbohydrate metabolism, signal transduction and transcription factors.
Collapse
Affiliation(s)
- Amanda R De La Torre
- Centre for Forest Conservation Genetics, Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, V6T1Z4, Canada; Department of Ecology and Environmental Science, Umeå University, Linneaus väg 6, SE-901 87, Umeå, Sweden
| | | | | |
Collapse
|
18
|
De La Torre AR, Wang T, Jaquish B, Aitken SN. Adaptation and exogenous selection in a Picea glauca × Picea engelmannii hybrid zone: implications for forest management under climate change. New Phytol 2014; 201:687-699. [PMID: 24200028 PMCID: PMC4285121 DOI: 10.1111/nph.12540] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/13/2013] [Accepted: 08/30/2013] [Indexed: 05/04/2023]
Abstract
The nature of selection responsible for the maintenance of the economically and ecologically important Picea glauca × Picea engelmannii hybrid zone was investigated. Genomic, phenotypic and climatic data were used to test assumptions of hybrid zone maintenance and to model future scenarios under climate change. Genome-wide estimates of admixture based on a panel of 86 candidate gene single nucleotide polymorphisms were combined with long-term quantitative data on growth and survival (over 20 yr), as well as one-time assessments of bud burst and bud set phenology, and cold hardiness traits. A total of 15,498 individuals were phenotyped for growth and survival. Our results suggest that the P. glauca × P. engelmannii hybrid zone is maintained by local adaptation to growing season length and snowpack (exogenous selection). Hybrids appeared to be fitter than pure species in intermediate environments, which fits expectations of the bounded hybrid superiority model of hybrid zone maintenance. Adaptive introgression from parental species has probably contributed to increased hybrid fitness in intermediate habitats. While P. engelmannii ancestry is higher than P. glauca ancestry in hybrid populations, on average, selective breeding in managed hybrid populations is shifting genomic composition towards P. glauca, potentially pre-adapting managed populations to warmer climates.
Collapse
Affiliation(s)
- Amanda R De La Torre
- Centre for Forest Conservation Genetics, Department of Forest and Conservation Sciences, University of British Columbia3041–2424 Main Mall, Vancouver, BC, V6T1Z4, Canada
- Department of Ecology and Environmental Sciences, Umeå UniversityLinneaus väg 6, SE-901 87, Umeå, Sweden
- Author for correspondence:Amanda R. De La Torre, Tel: +46 090 7865475,
| | - Tongli Wang
- Centre for Forest Conservation Genetics, Department of Forest and Conservation Sciences, University of British Columbia3041–2424 Main Mall, Vancouver, BC, V6T1Z4, Canada
| | - Barry Jaquish
- Kalamalka Forestry Centre, Tree Improvement Branch, BC Ministry of Forests, Lands & Natural Resource Operations3401 Reservoir Rd, Vernon, BC, V1B2C7, Canada
| | - Sally N Aitken
- Centre for Forest Conservation Genetics, Department of Forest and Conservation Sciences, University of British Columbia3041–2424 Main Mall, Vancouver, BC, V6T1Z4, Canada
| |
Collapse
|