1
|
Woudstra Y, Tumas H, van Ghelder C, Hung TH, Ilska JJ, Girardi S, A’Hara S, McLean P, Cottrell J, Bohlmann J, Bousquet J, Birol I, Woolliams JA, MacKay JJ. Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome. Genome Biol Evol 2024; 16:evae113. [PMID: 38787537 PMCID: PMC11171428 DOI: 10.1093/gbe/evae113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Revised: 04/23/2024] [Accepted: 05/21/2024] [Indexed: 05/25/2024] Open
Abstract
Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.
Collapse
Affiliation(s)
| | - Hayley Tumas
- Department of Biology, University of Oxford, Oxford OX1 3RB, UK
| | - Cyril van Ghelder
- INRAE, Université Côte d’Azur, CNRS, ISA, Sophia Antipolis 06903, France
| | - Tin Hang Hung
- Department of Biology, University of Oxford, Oxford OX1 3RB, UK
| | - Joana J Ilska
- The Roslin Institute, Royal (Dick) School of Veterinary Science, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - Sebastien Girardi
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC, Canada G1V 0A6
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC, Canada GIV 0A6
| | - Stuart A’Hara
- Forest Research, Northern Research Station, Roslin, Midlothian EH25 9SY, UK
| | - Paul McLean
- Forest Research, Northern Research Station, Roslin, Midlothian EH25 9SY, UK
| | - Joan Cottrell
- Forest Research, Northern Research Station, Roslin, Midlothian EH25 9SY, UK
| | - Joerg Bohlmann
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC, Canada V6T 1Z4
| | - Jean Bousquet
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC, Canada G1V 0A6
| | - Inanc Birol
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, Canada V5Z 4S6
| | - John A Woolliams
- The Roslin Institute, Royal (Dick) School of Veterinary Science, University of Edinburgh, Easter Bush, Midlothian EH25 9RG, UK
| | - John J MacKay
- Department of Biology, University of Oxford, Oxford OX1 3RB, UK
| |
Collapse
|
2
|
Tumas H, Ilska JJ, Gérardi S, Laroche J, A’Hara S, Boyle B, Janes M, McLean P, Lopez G, Lee SJ, Cottrell J, Gorjanc G, Bousquet J, Woolliams JA, MacKay JJ. High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers. G3 (BETHESDA, MD.) 2024; 14:jkae020. [PMID: 38366548 PMCID: PMC10989875 DOI: 10.1093/g3journal/jkae020] [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: 12/04/2023] [Accepted: 01/03/2024] [Indexed: 02/18/2024]
Abstract
In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.
Collapse
Affiliation(s)
- Hayley Tumas
- Department of Biology, University of Oxford, Oxford OX1 3RB, UK
| | - Joana J Ilska
- The Roslin Institute, Royal (Dick) School of Veterinary Science, University of Edinburgh, Midlothian EH25 9RG, UK
| | - Sebastien Gérardi
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC GIV 0A6, Canada
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC GIV 0A6, Canada
| | - Jerome Laroche
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC GIV 0A6, Canada
| | - Stuart A’Hara
- Forest Research, Northern Research Station, Midlothian EH25 9SY, UK
| | - Brian Boyle
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC GIV 0A6, Canada
| | - Mateja Janes
- The Roslin Institute, Royal (Dick) School of Veterinary Science, University of Edinburgh, Midlothian EH25 9RG, UK
| | - Paul McLean
- Forest Research, Northern Research Station, Midlothian EH25 9SY, UK
| | - Gustavo Lopez
- Forest Research, Northern Research Station, Midlothian EH25 9SY, UK
| | - Steve J Lee
- Forest Research, Northern Research Station, Midlothian EH25 9SY, UK
| | - Joan Cottrell
- Forest Research, Northern Research Station, Midlothian EH25 9SY, UK
| | - Gregor Gorjanc
- The Roslin Institute, Royal (Dick) School of Veterinary Science, University of Edinburgh, Midlothian EH25 9RG, UK
| | - Jean Bousquet
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC GIV 0A6, Canada
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC GIV 0A6, Canada
| | - John A Woolliams
- The Roslin Institute, Royal (Dick) School of Veterinary Science, University of Edinburgh, Midlothian EH25 9RG, UK
| | - John J MacKay
- Department of Biology, University of Oxford, Oxford OX1 3RB, UK
| |
Collapse
|
3
|
Wang X, Ingvarsson PK. Quantifying adaptive evolution and the effects of natural selection across the Norway spruce genome. Mol Ecol 2023; 32:5288-5304. [PMID: 37622583 DOI: 10.1111/mec.17106] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Revised: 08/07/2023] [Accepted: 08/09/2023] [Indexed: 08/26/2023]
Abstract
Detecting natural selection is one of the major goals of evolutionary genomics. Here, we sequenced the whole genome of 25 Picea abies individuals and quantified the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we showed that both negative selection and the rate of positively selected substitutions are very limited in coding regions. We found a positive correlation between the rate of adaptive substitutions and recombination rate and a negative correlation between the rate of adaptive substitutions and gene density, suggesting a widespread influence from Hill-Robertson interference on the efficiency of protein adaptation in P. abies. Finally, the distinct population statistics between genomic regions under either positive or balancing selection with that under neutral regions indicated the impact of natural selection on the genomic architecture of Norway spruce. Further gene ontology enrichment analysis for genes located in regions identified as undergoing either positive or long-term balancing selection also highlighted the specific molecular functions and biological processes that appear to be targets of selection in Norway spruce.
Collapse
Affiliation(s)
- Xi Wang
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden
| | - Pär K Ingvarsson
- Linnean Centre for Plant Biology, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
4
|
Kakui H, Ujino-Ihara T, Hasegawa Y, Tsurisaki E, Futamura N, Iwai J, Higuchi Y, Fujino T, Suzuki Y, Kasahara M, Yamaguchi K, Shigenobu S, Otani M, Nakano M, Nameta M, Shibata S, Ueno S, Moriguchi Y. A single-nucleotide substitution of CjTKPR1 determines pollen production in the gymnosperm plant Cryptomeria japonica. PNAS NEXUS 2023; 2:pgad236. [PMID: 37559748 PMCID: PMC10408704 DOI: 10.1093/pnasnexus/pgad236] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 05/05/2023] [Accepted: 07/07/2023] [Indexed: 08/11/2023]
Abstract
Pollinosis, also known as pollen allergy or hay fever, is a global problem caused by pollen produced by various plant species. The wind-pollinated Japanese cedar (Cryptomeria japonica) is the largest contributor to severe pollinosis in Japan, where increasing proportions of people have been affected in recent decades. The MALE STERILITY 4 (MS4) locus of Japanese cedar controls pollen production, and its homozygous mutants (ms4/ms4) show abnormal pollen development after the tetrad stage and produce no mature pollen. In this study, we narrowed down the MS4 locus by fine mapping in Japanese cedar and found TETRAKETIDE α-PYRONE REDUCTASE 1 (TKPR1) gene in this region. Transformation experiments using Arabidopsis thaliana showed that single-nucleotide substitution ("T" to "C" at 244-nt position) of CjTKPR1 determines pollen production. Broad conservation of TKPR1 beyond plant division could lead to the creation of pollen-free plants not only for Japanese cedar but also for broader plant species.
Collapse
Affiliation(s)
- Hiroyuki Kakui
- Graduate School of Science and Technology, Niigata University, Niigata 950-2181, Japan
- Institute for Sustainable Agro-ecosystem Services, Graduate School of Agricultural and Life Science, University of Tokyo, Tokyo 188-0002, Japan
| | - Tokuko Ujino-Ihara
- Department of Forest Molecular Genetics and Biotechnology, Forestry and Forest Products Research Institute, Forest Research and Management Organization, Ibaraki 305-8687, Japan
| | - Yoichi Hasegawa
- Department of Forest Molecular Genetics and Biotechnology, Forestry and Forest Products Research Institute, Forest Research and Management Organization, Ibaraki 305-8687, Japan
| | - Eriko Tsurisaki
- Graduate School of Science and Technology, Niigata University, Niigata 950-2181, Japan
| | - Norihiro Futamura
- Department of Forest Molecular Genetics and Biotechnology, Forestry and Forest Products Research Institute, Forest Research and Management Organization, Ibaraki 305-8687, Japan
| | - Junji Iwai
- Forest and Forestry Technology Division, Niigata Prefectural Forest Research Institute, Niigata 958-0264, Japan
| | - Yuumi Higuchi
- Forest and Forestry Technology Division, Niigata Prefectural Forest Research Institute, Niigata 958-0264, Japan
| | - Takeshi Fujino
- Graduate School of Frontier Sciences, The University of Tokyo, Chiba 277-8561, Japan
| | - Yutaka Suzuki
- Graduate School of Frontier Sciences, The University of Tokyo, Chiba 277-8561, Japan
| | - Masahiro Kasahara
- Graduate School of Frontier Sciences, The University of Tokyo, Chiba 277-8561, Japan
| | - Katsushi Yamaguchi
- Trans-Scale Biology Center, National Institute for Basic Biology, Aichi 444-8585, Japan
| | - Shuji Shigenobu
- Trans-Scale Biology Center, National Institute for Basic Biology, Aichi 444-8585, Japan
| | - Masahiro Otani
- Faculty of Agriculture, Niigata University, Niigata 950-2181, Japan
| | - Masaru Nakano
- Faculty of Agriculture, Niigata University, Niigata 950-2181, Japan
| | - Masaaki Nameta
- Graduate School of Medical and Dental Sciences, Niigata University, Niigata 951-8122, Japan
| | - Shinsuke Shibata
- Graduate School of Medical and Dental Sciences, Niigata University, Niigata 951-8122, Japan
| | - Saneyoshi Ueno
- Department of Forest Molecular Genetics and Biotechnology, Forestry and Forest Products Research Institute, Forest Research and Management Organization, Ibaraki 305-8687, Japan
| | | |
Collapse
|
5
|
Tiret M, Olsson L, Grahn T, Karlsson B, Milesi P, Lascoux M, Lundqvist S, García‐Gil MR. Divergent selection predating the Last Glacial Maximum mainly acted on macro-phenotypes in Norway spruce. Evol Appl 2022; 16:163-172. [PMID: 36699125 PMCID: PMC9850012 DOI: 10.1111/eva.13519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 11/27/2022] [Accepted: 12/03/2022] [Indexed: 12/24/2022] Open
Abstract
The current distribution and population structure of many species were, to a large extent, shaped by cycles of isolation in glacial refugia and subsequent population expansions. Isolation in and postglacial expansion through heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a constant interplay between evolutionary and demographic processes. We investigated population differentiation and adaptation of Norway spruce for juvenile growth, diameter of the stem, wood density, and tracheid traits at breast height. Data from 4461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of Q ST vs. F ST. We show that the macroscopic resultant trait (stem diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection that predated the Last Glacial Maximum. Altogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly preadapted, mainly through growth-related traits.
Collapse
Affiliation(s)
- Mathieu Tiret
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala UniversityUppsalaSweden,Department of Forest Genetics and Plant PhysiologySLU, Umeå Plant Science Centre (UPSC)UmeåSweden,IGEPP, INRAE, Institut Agro, Université de RennesDomaine de la MotteLe RheuFrance
| | | | | | | | - Pascal Milesi
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala UniversityUppsalaSweden
| | - Martin Lascoux
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala UniversityUppsalaSweden
| | | | - Maria Rosario García‐Gil
- Department of Forest Genetics and Plant PhysiologySLU, Umeå Plant Science Centre (UPSC)UmeåSweden
| |
Collapse
|
6
|
Li L, Milesi P, Tiret M, Chen J, Sendrowski J, Baison J, Chen Z, Zhou L, Karlsson B, Berlin M, Westin J, Garcia‐Gil MR, Wu HX, Lascoux M. Teasing apart the joint effect of demography and natural selection in the birth of a contact zone. THE NEW PHYTOLOGIST 2022; 236:1976-1987. [PMID: 36093739 PMCID: PMC9828440 DOI: 10.1111/nph.18480] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 08/23/2022] [Indexed: 05/26/2023]
Abstract
Vast population movements induced by recurrent climatic cycles have shaped the genetic structure of plant species. During glacial periods species were confined to low-latitude refugia from which they recolonized higher latitudes as the climate improved. This multipronged recolonization led to many lineages that later met and formed large contact zones. We utilize genomic data from 5000 Picea abies trees to test for the presence of natural selection during recolonization and establishment of a contact zone in Scandinavia. Scandinavian P. abies is today made up of a southern genetic cluster originating from the Baltics, and a northern one originating from Northern Russia. The contact zone delineating them closely matches the limit between two major climatic regions. We show that natural selection contributed to its establishment and maintenance. First, an isolation-with-migration model with genome-wide linked selection fits the data better than a purely neutral one. Second, many loci show signatures of selection or are associated with environmental variables. These loci, regrouped in clusters on chromosomes, are often related to phenology. Altogether, our results illustrate how climatic cycles, recolonization and selection can establish strong local adaptation along contact zones and affect the genetic architecture of adaptive traits.
Collapse
Affiliation(s)
- Lili Li
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala University75236UppsalaSweden
| | - Pascal Milesi
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala University75236UppsalaSweden
| | - Mathieu Tiret
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala University75236UppsalaSweden
| | - Jun Chen
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala University75236UppsalaSweden
- College of Life SciencesZhejiang UniversityHangzhouZhejiang310058China
| | - Janek Sendrowski
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala University75236UppsalaSweden
| | - John Baison
- Department Forest Genetics and Plant Physiology, Umeå Plant Science CentreSwedish University of Agricultural SciencesUmeåSE‐90183Sweden
| | - Zhi‐qiang Chen
- Department Forest Genetics and Plant Physiology, Umeå Plant Science CentreSwedish University of Agricultural SciencesUmeåSE‐90183Sweden
| | - Linghua Zhou
- Department Forest Genetics and Plant Physiology, Umeå Plant Science CentreSwedish University of Agricultural SciencesUmeåSE‐90183Sweden
| | | | - Mats Berlin
- SkogforskUppsala Science Park751 83UppsalaSweden
| | - Johan Westin
- Unit for Field‐Based Forest ResearchSwedish University of Agricultural SciencesSE‐922 91VindelnSweden
| | - Maria Rosario Garcia‐Gil
- Department Forest Genetics and Plant Physiology, Umeå Plant Science CentreSwedish University of Agricultural SciencesUmeåSE‐90183Sweden
| | - Harry X. Wu
- Department Forest Genetics and Plant Physiology, Umeå Plant Science CentreSwedish University of Agricultural SciencesUmeåSE‐90183Sweden
- CSIRO National Collection Research AustraliaBlack Mountain LaboratoryCanberraACT2601Australia
| | - Martin Lascoux
- Program in Plant Ecology and Evolution, Department of Ecology and Genetics, EBC and SciLife LabUppsala University75236UppsalaSweden
| |
Collapse
|
7
|
Freeman JS, Slavov GT, Butler JB, Frickey T, Graham NJ, Klápště J, Lee J, Telfer EJ, Wilcox P, Dungey HS. High density linkage maps, genetic architecture, and genomic prediction of growth and wood properties in Pinus radiata. BMC Genomics 2022; 23:731. [PMID: 36307760 PMCID: PMC9617409 DOI: 10.1186/s12864-022-08950-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Accepted: 10/10/2022] [Indexed: 11/19/2022] Open
Abstract
Background The growing availability of genomic resources in radiata pine paves the way for significant advances in fundamental and applied genomic research. We constructed robust high-density linkage maps based on exome-capture genotyping in two F1 populations, and used these populations to perform quantitative trait locus (QTL) scans, genomic prediction and quantitative analyses of genetic architecture for key traits targeted by tree improvement programmes. Results Our mapping approach used probabilistic error correction of the marker data, followed by an iterative approach based on stringent parameters. This approach proved highly effective in producing high-density maps with robust marker orders and realistic map lengths (1285–4674 markers per map, with sizes ranging from c. 1643–2292 cM, and mean marker intervals of 0.7–2.1 cM). Colinearity was high between parental linkage maps, although there was evidence for a large chromosomal rearrangement (affecting ~ 90 cM) in one of the parental maps. In total, 28 QTL were detected for growth (stem diameter) and wood properties (wood density and fibre properties measured by Silviscan) in the QTL discovery population, with 1–3 QTL of small to moderate effect size detected per trait in each parental map. Four of these QTL were validated in a second, unrelated F1 population. Results from genomic prediction and analyses of genetic architecture were consistent with those from QTL scans, with wood properties generally having moderate to high genomic heritabilities and predictive abilities, as well as somewhat less complex genetic architectures, compared to growth traits. Conclusions Despite the economic importance of radiata pine as a plantation forest tree, robust high-density linkage maps constructed from reproducible, sequence-anchored markers have not been published to date. The maps produced in this study will be a valuable resource for several applications, including the selection of marker panels for genomic prediction and anchoring a recently completed de novo whole genome assembly. We also provide the first map-based evidence for a large genomic rearrangement in radiata pine. Finally, results from our QTL scans, genomic prediction, and genetic architecture analyses are informative about the genomic basis of variation in important phenotypic traits. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08950-6.
Collapse
|
8
|
Gagalova KK, Warren RL, Coombe L, Wong J, Nip KM, Yuen MMS, Whitehill JGA, Celedon JM, Ritland C, Taylor GA, Cheng D, Plettner P, Hammond SA, Mohamadi H, Zhao Y, Moore RA, Mungall AJ, Boyle B, Laroche J, Cottrell J, Mackay JJ, Lamothe M, Gérardi S, Isabel N, Pavy N, Jones SJM, Bohlmann J, Bousquet J, Birol I. Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2022; 111:1469-1485. [PMID: 35789009 DOI: 10.1111/tpj.15889] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Revised: 06/22/2022] [Accepted: 06/27/2022] [Indexed: 06/15/2023]
Abstract
Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.
Collapse
Affiliation(s)
- Kristina K Gagalova
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - René L Warren
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Lauren Coombe
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Johnathan Wong
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Ka Ming Nip
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Macaire Man Saint Yuen
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
| | - Justin G A Whitehill
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
| | - Jose M Celedon
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
| | - Carol Ritland
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
| | - Greg A Taylor
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Dean Cheng
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Patrick Plettner
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - S Austin Hammond
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
- Next-Generation Sequencing Facility, University of Saskatchewan, Saskatoon, SK, S7N 5E5, Canada
| | - Hamid Mohamadi
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Yongjun Zhao
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Richard A Moore
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Andrew J Mungall
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Brian Boyle
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC, GIV 0A6, Canada
| | - Jérôme Laroche
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC, GIV 0A6, Canada
| | - Joan Cottrell
- Forest Research, U.K. Forestry Commission, Northern Research Station, Roslin, EH25 9SY, Midlothian, UK
| | - John J Mackay
- Department of Plant Sciences, University of Oxford, Oxford, OX1 3RB, UK
| | - Manuel Lamothe
- Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Québec, QC, G1V 4C7, Canada
| | - Sébastien Gérardi
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC, GIV 0A6, Canada
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Nathalie Isabel
- Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Québec, QC, G1V 4C7, Canada
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Nathalie Pavy
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC, GIV 0A6, Canada
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Steven J M Jones
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| | - Joerg Bohlmann
- Michael Smith Laboratories, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
| | - Jean Bousquet
- Institute for Systems and Integrative Biology, Université Laval, Québec, QC, GIV 0A6, Canada
- Canada Research Chair in Forest Genomics, Forest Research Centre, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Inanc Birol
- Canada's Michael Smith Genome Sciences Centre, Vancouver, BC, V5Z 4S6, Canada
| |
Collapse
|
9
|
Jasper RJ, McDonald TK, Singh P, Lu M, Rougeux C, Lind BM, Yeaman S. Evaluating the accuracy of variant calling methods using the frequency of parent-offspring genotype mismatch. Mol Ecol Resour 2022; 22:2524-2533. [PMID: 35510784 PMCID: PMC9544674 DOI: 10.1111/1755-0998.13628] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Revised: 03/11/2022] [Accepted: 03/29/2022] [Indexed: 11/28/2022]
Abstract
The use of next‐generation sequencing (NGS) data sets has increased dramatically over the last decade, but there have been few systematic analyses quantifying the accuracy of the commonly used variant caller programs. Here we used a familial design consisting of diploid tissue from a single lodgepole pine (Pinus contorta) parent and the maternally derived haploid tissue from 106 full‐sibling offspring, where mismatches could only arise due to mutation or bioinformatic error. Given the rarity of mutation, we used the rate of mismatches between parent and offspring genotype calls to infer the single nucleotide polymorphism (SNP) genotyping error rates of FreeBayes, HaplotypeCaller, SAMtools, UnifiedGenotyper, and VarScan. With baseline filtering HaplotypeCaller and UnifiedGenotyper yielded more SNPs and higher error rates by one to two orders of magnitude, whereas FreeBayes, SAMtools and VarScan yielded lower numbers of SNPs and more modest error rates. To facilitate comparison between variant callers we standardized each SNP set to the same number of SNPs using additional filtering, where UnifiedGenotyper consistently produced the smallest proportion of genotype errors, followed by HaplotypeCaller, VarScan, SAMtools, and FreeBayes. Additionally, we found that error rates were minimized for SNPs called by more than one variant caller. Finally, we evaluated the performance of various commonly used filtering metrics on SNP calling. Our analysis provides a quantitative assessment of the accuracy of five widely used variant calling programs and offers valuable insights into both the choice of variant caller program and the choice of filtering metrics, especially for researchers using non‐model study systems.
Collapse
Affiliation(s)
- Russ J Jasper
- Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | | | - Pooja Singh
- Department of Biological Sciences, University of Calgary, Calgary, AB, Canada.,Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland.,3EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
| | - Mengmeng Lu
- Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Clément Rougeux
- Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Brandon M Lind
- Centre for Forest Conservation Genetics and Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, BC, Canada
| | - Sam Yeaman
- Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| |
Collapse
|
10
|
Redondo MA, Oliva J, Elfstrand M, Boberg J, Capador-Barreto HD, Karlsson B, Berlin A. Host genotype interacts with aerial spore communities and influences the needle mycobiome of Norway spruce. Environ Microbiol 2022; 24:3640-3654. [PMID: 35315253 PMCID: PMC9544151 DOI: 10.1111/1462-2920.15974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Revised: 03/07/2022] [Accepted: 03/08/2022] [Indexed: 11/28/2022]
Abstract
The factors shaping the composition of the tree mycobiome are still under investigation. We tested the effects of host genotype, site, host phenotypic traits, and air fungal spore communities on the assembly of the fungi inhabiting Norway spruce needles. We used Norway spruce clones and spore traps within the collection sites and characterized both needle and air mycobiome communities by high‐throughput sequencing of the ITS2 region. The composition of the needle mycobiome differed between Norway spruce clones, and clones with high genetic similarity had a more similar mycobiome. The needle mycobiome also varied across sites and was associated with the composition of the local air mycobiome and climate. Phenotypic traits such as diameter at breast height or crown health influenced the needle mycobiome to a lesser extent than host genotype and air mycobiome. Altogether, our results suggest that the needle mycobiome is mainly driven by the host genotype in combination with the composition of the local air spore communities. Our work highlights the role of host intraspecific variation in shaping the mycobiome of trees and provides new insights on the ecological processes structuring fungal communities inhabiting woody plants.
Collapse
Affiliation(s)
- Miguel A Redondo
- Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Box 7026, 750 07, Sweden
| | - Jonàs Oliva
- Department of Crop and Forest Sciences, University of Lleida, Alcalde Rovira Roure 191, Lleida, 25198, Spain.,Joint Research Unit CTFC-AGROTECNIO, Alcalde Rovira Roure 191, Lleida, 25198, Spain
| | - Malin Elfstrand
- Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Box 7026, 750 07, Sweden
| | - Johanna Boberg
- Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Box 7026, 750 07, Sweden
| | - Hernán D Capador-Barreto
- Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Box 7026, 750 07, Sweden
| | - Bo Karlsson
- Skogforsk, Svalöv, Ekebo 2250, 268 90, Sweden
| | - Anna Berlin
- Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Box 7026, 750 07, Sweden
| |
Collapse
|
11
|
Wang D, Yang L, Shi C, Li S, Tang H, He C, Cai N, Duan A, Gong H. QTL mapping for growth-related traits by constructing the first genetic linkage map in Simao pine. BMC PLANT BIOLOGY 2022; 22:48. [PMID: 35065611 PMCID: PMC8783431 DOI: 10.1186/s12870-022-03425-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2020] [Accepted: 01/04/2022] [Indexed: 05/31/2023]
Abstract
BACKGROUND Simao pine is one of the primary economic tree species for resin and timber production in southwest China. The exploitation and utilization of Simao pine are constrained by the relatively lacking of genetic information. Construction a fine genetic linkage map and detecting quantitative trait locis (QTLs) for growth-related traits is a prerequisite section of Simao Pine's molecular breeding program. RESULTS In our study, a high-resolution Simao pine genetic map employed specific locus amplified fragment sequencing (SLAF-seq) technology and based on an F1 pseudo-testcross population has been constructed. There were 11,544 SNPs assigned to 12 linkage groups (LGs), and the total length of the map was 2,062.85 cM with a mean distance of 0.37 cM between markers. According to the phenotypic variation analysis for three consecutive years, a total of seventeen QTLs for four traits were detected. Among 17 QTLs, there were six for plant height (Dh.16.1, Dh16.2, Dh17.1, Dh18.1-3), five for basal diameter (Dbd.17.1-5), four for needle length (Dnl17.1-3, Dnl18.1) and two for needle diameter (Dnd17.1 and Dnd18.1) respectively. These QTLs individually explained phenotypic variance from 11.0-16.3%, and the logarithm of odds (LOD) value ranged from 2.52 to 3.87. CONCLUSIONS In our study, a fine genetic map of Simao pine applied the technology of SLAF-seq has been constructed for the first time. Based on the map, a total of 17 QTLs for four growth-related traits were identified. It provides helpful information for genomic studies and marker-assisted selection (MAS) in Simao pine.
Collapse
Affiliation(s)
- Dawei Wang
- Key Laboratory for Forest Resource Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, Southwest Forestry University, Kunming, China
- Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming, China
| | - Lin Yang
- Key Laboratory for Forest Resource Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, Southwest Forestry University, Kunming, China
- Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming, China
| | - Chen Shi
- Key Laboratory for Forest Resource Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, Southwest Forestry University, Kunming, China
- Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming, China
| | - Siguang Li
- Yunnan Academy of Forestry, Kunming, China
| | - Hongyan Tang
- Puer City Institute of Forestry Sciences, Puer, China
| | - Chengzhong He
- Key Laboratory for Forest Resource Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, Southwest Forestry University, Kunming, China
- Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming, China
| | - Nianhui Cai
- Key Laboratory for Forest Resource Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, Southwest Forestry University, Kunming, China
- Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming, China
| | - Anan Duan
- Key Laboratory for Forest Resource Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, Southwest Forestry University, Kunming, China
- Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming, China
| | - Hede Gong
- School of Geography, Southwest Forestry University, Kunming, China.
| |
Collapse
|
12
|
Zhao J, Li H, Xu Y, Yin Y, Huang T, Zhang B, Wang Y, Li Y, Cao Y, An W. A consensus and saturated genetic map provides insight into genome anchoring, synteny of Solanaceae and leaf- and fruit-related QTLs in wolfberry (Lycium Linn.). BMC PLANT BIOLOGY 2021; 21:350. [PMID: 34303361 PMCID: PMC8306383 DOI: 10.1186/s12870-021-03115-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2021] [Accepted: 06/22/2021] [Indexed: 06/02/2023]
Abstract
BACKGROUND Lycium Linn. (Solanaceae) is a genus of economically important plants producing fruits and leaves with high nutritional value and medicinal benefits. However, genetic analysis of this plant and molecular breeding for quality improvement are limited by the lack of sufficient molecular markers. RESULTS In this study, two parental strains, 'Ningqi No. 1' (Lycium barbarum L.) and 'Yunnan Gouqi' (Lycium yunnanense Kuang et A.M. Lu), and 200 F1 hybrid individuals were resequenced for genetic analysis. In total, 8,507 well-selected SNPs were developed, and a high-density genetic map (NY map) was constructed with a total genetic distance of 2,122.24 cM. A consensus genetic map was established by integrating the NY map and a previously published genetic map (NC map) containing 15,240 SNPs, with a total genetic distance of 3,058.19 cM and an average map distance of 0.21 cM. The 12 pseudochromosomes of the Lycium reference genome were anchored using this consensus genetic map, with an anchoring rate of 64.3%. Moreover, weak collinearities between the consensus map and the pepper, potato, and tomato genomes were observed. Twenty-five stable QTLs were identified for leaf- and fruit-related phenotypes, including fruit weight, fruit longitude, leaf length, the fruit index, and the leaf index; these stable QTLs were mapped to four different linkage groups, with LOD scores ranging from 2.51 to 19.37 and amounts of phenotypic variance explained from 6.2% to 51.9%. Finally, 82 out of 188 predicted genes underlying stable QTLs for fruit-related traits were differentially expressed according to RNA-seq analysis. CONCLUSIONS A chromosome-level assembly can provide a foundation for further functional genomics research for wolfberry. The genomic regions of these stably expressed QTLs could be used as targets for further fine mapping and development of molecular markers for marker-assisted selection (MAS). The present study provided valuable information on saturated SNP markers and reliable QTLs for map-based cloning of functional genes related to yield and morphological traits in Lycium spp.
Collapse
Affiliation(s)
- Jianhua Zhao
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| | - Haoxia Li
- Desertification Control Research Institute, Ningxia Academy of Agriculture and Forestry Sciences, Yinchuan, 750002 China
| | - Yuhui Xu
- Adsen Biotechnology Co., Ltd, Urumchi, 830022 China
| | - Yue Yin
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| | - Ting Huang
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| | - Bo Zhang
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| | - Yajun Wang
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| | - Yanlong Li
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| | - Youlong Cao
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| | - Wei An
- Wolfberry Science Research Institute, Ningxia Academy of Agriculture and Forestry Sciences/National Wolfberry Engineering Research Center, Yinchuan, 750002 China
| |
Collapse
|
13
|
Capador-Barreto HD, Bernhardsson C, Milesi P, Vos I, Lundén K, Wu HX, Karlsson B, Ingvarsson PK, Stenlid J, Elfstrand M. Killing two enemies with one stone? Genomics of resistance to two sympatric pathogens in Norway spruce. Mol Ecol 2021; 30:4433-4447. [PMID: 34218489 DOI: 10.1111/mec.16058] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2021] [Revised: 06/22/2021] [Accepted: 06/29/2021] [Indexed: 12/31/2022]
Abstract
Trees must cope with the attack of multiple pathogens, often simultaneously during their long lifespan. Ironically, the genetic and molecular mechanisms controlling this process are poorly understood. The objective of this study was to compare the genetic component of resistance in Norway spruce to Heterobasidion annosum s.s. and its sympatric congener Heterobasidion parviporum. Heterobasidion root- and stem-rot is a major disease of Norway spruce caused by members of the Heterobasidion annosum species complex. Resistance to both pathogens was measured using artificial inoculations in half-sib families of Norway spruce trees originating from central to northern Europe. The genetic component of resistance was analysed using 63,760 genome-wide exome-capture sequenced SNPs and multitrait genome-wide associations. No correlation was found for resistance to the two pathogens; however, associations were found between genomic variants and resistance traits with synergic or antagonist pleiotropic effects to both pathogens. Additionally, a latitudinal cline in resistance in the bark to H. annosum s.s. was found; trees from southern latitudes, with a later bud-set and thicker stem diameter, allowed longer lesions, but this was not the case for H. parviporum. In summary, this study detects genomic variants with pleiotropic effects which explain multiple disease resistance from a genic level and could be useful for selection of resistant trees to both pathogens. Furthermore, it highlights the need for additional research to understand the evolution of resistance traits to multiple pathogens in trees.
Collapse
Affiliation(s)
- Hernán D Capador-Barreto
- Uppsala Biocentre, Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Carolina Bernhardsson
- Uppsala Biocentre, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Pascal Milesi
- Department of Ecology and Genetics, Evolutionary Biology Centre, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Ingrid Vos
- Forestry Research Institute of Sweden (Skogforsk), Ekebo, Sweden
| | - Karl Lundén
- Uppsala Biocentre, Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Harry X Wu
- Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, Umeå, Sweden
| | - Bo Karlsson
- Forestry Research Institute of Sweden (Skogforsk), Ekebo, Sweden
| | - Pär K Ingvarsson
- Uppsala Biocentre, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Jan Stenlid
- Uppsala Biocentre, Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Malin Elfstrand
- Uppsala Biocentre, Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
14
|
Wang X, Bernhardsson C, Ingvarsson PK. Demography and Natural Selection Have Shaped Genetic Variation in the Widely Distributed Conifer Norway Spruce (Picea abies). Genome Biol Evol 2020; 12:3803-3817. [PMID: 31958121 PMCID: PMC7046165 DOI: 10.1093/gbe/evaa005] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/14/2020] [Indexed: 12/13/2022] Open
Abstract
Under the neutral theory, species with larger effective population size are expected to harbor higher genetic diversity. However, across a wide variety of organisms, the range of genetic diversity is orders of magnitude more narrow than the range of effective population size. This observation has become known as Lewontin’s paradox and although aspects of this phenomenon have been extensively studied, the underlying causes for the paradox remain unclear. Norway spruce (Picea abies) is a widely distributed conifer species across the northern hemisphere, and it consequently plays a major role in European forestry. Here, we use whole-genome resequencing data from 35 individuals to perform population genomic analyses in P. abies in an effort to understand what drives genome-wide patterns of variation in this species. Despite having a very wide geographic distribution and an corresponding enormous current population size, our analyses find that genetic diversity of P. abies is low across a number of populations (π = 0.0049 in Central-Europe, π = 0.0063 in Sweden-Norway, π = 0.0063 in Finland). To assess the reasons for the low levels of genetic diversity, we infer the demographic history of the species and find that it is characterized by several reoccurring bottlenecks with concomitant decreases in effective population size can, at least partly, provide an explanation for low polymorphism we observe in P. abies. Further analyses suggest that recurrent natural selection, both purifying and positive selection, can also contribute to the loss of genetic diversity in Norway spruce by reducing genetic diversity at linked sites. Finally, the overall low mutation rates seen in conifers can also help explain the low genetic diversity maintained in Norway spruce.
Collapse
Affiliation(s)
- Xi Wang
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, Sweden.,Linnean Centre for Plant Biology, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Carolina Bernhardsson
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, Sweden.,Linnean Centre for Plant Biology, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden.,Department of Organismal Biology, Uppsala University, Uppsala, Sweden
| | - Pär K Ingvarsson
- Linnean Centre for Plant Biology, Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
15
|
Bernhardsson C, Zan Y, Chen Z, Ingvarsson PK, Wu HX. Development of a highly efficient 50K single nucleotide polymorphism genotyping array for the large and complex genome of Norway spruce (Picea abies L. Karst) by whole genome resequencing and its transferability to other spruce species. Mol Ecol Resour 2020; 21:880-896. [PMID: 33179386 PMCID: PMC7984398 DOI: 10.1111/1755-0998.13292] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Revised: 10/23/2020] [Accepted: 11/04/2020] [Indexed: 12/30/2022]
Abstract
Norway spruce (Picea abies L. Karst) is one of the most important forest tree species with significant economic and ecological impact in Europe. For decades, genomic and genetic studies on Norway spruce have been challenging due to the large and repetitive genome (19.6 Gb with more than 70% being repetitive). To accelerate genomic studies, including population genetics, genome‐wide association studies (GWAS) and genomic selection (GS), in Norway spruce and related species, we here report on the design and performance of a 50K single nucleotide polymorphism (SNP) genotyping array for Norway spruce. The array is developed based on whole genome resequencing (WGS), making it the first WGS‐based SNP array in any conifer species so far. After identifying SNPs using genome resequencing data from 29 trees collected in northern Europe, we adopted a two‐step approach to design the array. First, we built a 450K screening array and used this to genotype a population of 480 trees sampled from both natural and breeding populations across the Norway spruce distribution range. These samples were then used to select high‐confidence probes that were put on the final 50K array. The SNPs selected are distributed over 45,552 scaffolds from the P. abies version 1.0 genome assembly and target 19,954 unique gene models with an even coverage of the 12 linkage groups in Norway spruce. We show that the array has a 99.5% probe specificity, >98% Mendelian allelic inheritance concordance, an average sample call rate of 96.30% and an SNP call rate of 98.90% in family trios and haploid tissues. We also observed that 23,797 probes (50%) could be identified with high confidence in three other spruce species (white spruce [Picea glauca], black spruce [P. mariana] and Sitka spruce [P. sitchensis]). The high‐quality genotyping array will be a valuable resource for genetic and genomic studies in Norway spruce as well as in other conifer species of the same genus.
Collapse
Affiliation(s)
- Carolina Bernhardsson
- Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden.,Department of Organismal Biology, Uppsala University, Uppsala, Sweden
| | - Yanjun Zan
- Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Science, Umeå, Sweden
| | - Zhiqiang Chen
- Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Science, Umeå, Sweden
| | - Pär K Ingvarsson
- Linnean Centre for Plant Biology, Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Science, Uppsala, Sweden
| | - Harry X Wu
- Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Science, Umeå, Sweden.,Beijing Advanced Innovation Centre for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China.,Black Mountain Laboratory, CSIRO National Research Collection Australia, Canberra, ACT, Australia
| |
Collapse
|
16
|
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] [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
|
17
|
Elfstrand M, Zhou L, Baison J, Olson Å, Lundén K, Karlsson B, Wu HX, Stenlid J, García‐Gil MR. Genotypic variation in Norway spruce correlates to fungal communities in vegetative buds. Mol Ecol 2020; 29:199-213. [PMID: 31755612 PMCID: PMC7003977 DOI: 10.1111/mec.15314] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 10/31/2019] [Accepted: 11/20/2019] [Indexed: 12/19/2022]
Abstract
The taxonomically diverse phyllosphere fungi inhabit leaves of plants. Thus, apart from the fungi's dispersal capacities and environmental factors, the assembly of the phyllosphere community associated with a given host plant depends on factors encoded by the host's genome. The host genetic factors and their influence on the assembly of phyllosphere communities under natural conditions are poorly understood, especially in trees. Recent work indicates that Norway spruce (Picea abies) vegetative buds harbour active fungal communities, but these are hitherto largely uncharacterized. This study combines internal transcribed spacer sequencing of the fungal communities associated with dormant vegetative buds with a genome-wide association study (GWAS) in 478 unrelated Norway spruce trees. The aim was to detect host loci associated with variation in the fungal communities across the population, and to identify loci correlating with the presence of specific, latent, pathogens. The fungal communities were dominated by known Norway spruce phyllosphere endophytes and pathogens. We identified six quantitative trait loci (QTLs) associated with the relative abundance of the dominating taxa (i.e., top 1% most abundant taxa). Three additional QTLs associated with colonization by the spruce needle cast pathogen Lirula macrospora or the cherry spruce rust (Thekopsora areolata) in asymptomatic tissues were detected. The identification of the nine QTLs shows that the genetic variation in Norway spruce influences the fungal community in dormant buds and that mechanisms underlying the assembly of the communities and the colonization of latent pathogens in trees may be uncovered by combining molecular identification of fungi with GWAS.
Collapse
Affiliation(s)
- Malin Elfstrand
- Uppsala BiocentreDepartment of Forest Mycology and Plant PathologySwedish University of Agricultural SciencesUppsalaSweden
| | - Linghua Zhou
- Umeå Plant Science CentreDepartment of Forest Genetics and Plant PhysiologySwedish University of Agricultural SciencesUmeåSweden
| | - John Baison
- Umeå Plant Science CentreDepartment of Forest Genetics and Plant PhysiologySwedish University of Agricultural SciencesUmeåSweden
| | - Åke Olson
- Uppsala BiocentreDepartment of Forest Mycology and Plant PathologySwedish University of Agricultural SciencesUppsalaSweden
| | - Karl Lundén
- Uppsala BiocentreDepartment of Forest Mycology and Plant PathologySwedish University of Agricultural SciencesUppsalaSweden
| | | | - Harry X. Wu
- Umeå Plant Science CentreDepartment of Forest Genetics and Plant PhysiologySwedish University of Agricultural SciencesUmeåSweden
| | - Jan Stenlid
- Uppsala BiocentreDepartment of Forest Mycology and Plant PathologySwedish University of Agricultural SciencesUppsalaSweden
| | - M. Rosario García‐Gil
- Umeå Plant Science CentreDepartment of Forest Genetics and Plant PhysiologySwedish University of Agricultural SciencesUmeåSweden
| |
Collapse
|
18
|
Genome-Wide Variant Identification and High-Density Genetic Map Construction Using RADseq for Platycladus orientalis (Cupressaceae). G3-GENES GENOMES GENETICS 2019; 9:3663-3672. [PMID: 31506321 PMCID: PMC6829139 DOI: 10.1534/g3.119.400684] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Platycladus orientalis is an ecologically important native conifer in Northern China and exotic species in many parts of the world; however, knowledge about the species’ genetics and genome are very limited. The availability of well-developed battery of genetic markers, with large genome coverage, is a prerequisite for the species genetic dissection of adaptive attributes and efficient selective breeding. Here, we present a genome-wide genotyping method with double-digestion restriction site associated DNA sequencing (ddRAD-seq) that is effective in generating large number of Mendelian markers for genome mapping and other genetic applications. Using 139 megagametophytes collected from a single mother tree, we assembled 397,226 loci, of which 108,683 (27.4%) were polymorphic. After stringent filtering for 1:1 segregation ratio and missing rate of <20%, the remaining 23,926 loci (22% of the polymorphic loci) were ordered into 11 linkage groups (LGs) and distributed across 7,559 unique positions, with a total map length of 1,443 cM and an average spacing of 0.2 cM between adjacent unique positions. The 11 LGs correspond to the species’ 11 haploid genome chromosome number. This genetic map is among few high-density maps available for conifers to date, and represents the first genetic map for P. orientalis. The information generated serves as a solid foundation not only for marker-assisted breeding efforts, but also for comparative conifer genomic studies.
Collapse
|
19
|
Baison J, Vidalis A, Zhou L, Chen Z, Li Z, Sillanpää MJ, Bernhardsson C, Scofield D, Forsberg N, Grahn T, Olsson L, Karlsson B, Wu H, Ingvarsson PK, Lundqvist S, Niittylä T, García‐Gil MR. Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2019; 100:83-100. [PMID: 31166032 PMCID: PMC6852177 DOI: 10.1111/tpj.14429] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 04/16/2019] [Accepted: 05/20/2019] [Indexed: 05/26/2023]
Abstract
Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.
Collapse
Affiliation(s)
- John Baison
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
| | - Amaryllis Vidalis
- Section of Population Epigenetics and EpigenomicsCentre of Life and Food Sciences WeihenstephanTechnische Universität MünchenLichtenbergstr. 2aMünchen85748Germany
| | - Linghua Zhou
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
| | - Zhi‐Qiang Chen
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
| | - Zitong Li
- Ecological Genetics Research UnitDepartment of BiosciencesUniversity of HelsinkiP.O. Box 65FI‐00014HelsinkiFinland
| | - Mikko J. Sillanpää
- Department of Mathematical SciencesBiocenter OuluUniversity of OuluPentti Kaiteran katu 1OuluFinland
| | - Carolina Bernhardsson
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
- Department of Ecology and Environmental ScienceUmeå UniversityLinnaeus väg 4-6Umeå907 36Sweden
| | - Douglas Scofield
- Uppsala Multidisciplinary Centre for Advanced Computational ScienceUppsala UniversityLägerhyddsvägen 2Uppsala752 37Sweden
| | - Nils Forsberg
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
| | - Thomas Grahn
- RISE BioeconomyDrottning Kristinas väg 61SE‐114 86StockholmSweden
| | - Lars Olsson
- RISE BioeconomyDrottning Kristinas väg 61SE‐114 86StockholmSweden
| | | | - Harry Wu
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
| | - Pär K. Ingvarsson
- Department of Ecology and Environmental ScienceUmeå UniversityLinnaeus väg 4-6Umeå907 36Sweden
- Department of Ecology and Genetics: Evolutionary BiologyUppsala UniversityKåbovägen 4Uppsala752 36Sweden
| | - Sven‐Olof Lundqvist
- RISE BioeconomyDrottning Kristinas väg 61SE‐114 86StockholmSweden
- IICRosenlundsgatan 48BSE‐118 63StockholmSweden
| | - Totte Niittylä
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
| | - M Rosario García‐Gil
- Department of Forest Genetics and Plant PhysiologyUmeå Plant Science CentreSwedish University of Agricultural ScienceParallellvägen 21Umeå907 36Sweden
| |
Collapse
|
20
|
Chen J, Li L, Milesi P, Jansson G, Berlin M, Karlsson B, Aleksic J, Vendramin GG, Lascoux M. Genomic data provide new insights on the demographic history and the extent of recent material transfers in Norway spruce. Evol Appl 2019; 12:1539-1551. [PMID: 31462913 PMCID: PMC6708423 DOI: 10.1111/eva.12801] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2019] [Revised: 03/28/2019] [Accepted: 04/02/2019] [Indexed: 12/23/2022] Open
Abstract
Primeval forests are today exceedingly rare in Europe, and transfer of forest reproductive material for afforestation and improvement has been very common, especially over the last two centuries. This can be a serious impediment when inferring past population movements in response to past climate changes such as the last glacial maximum (LGM), some 18,000 years ago. In the present study, we genotyped 1,672 individuals from three Picea species (P. abies, P. obovata, and P. omorika) at 400K SNPs using exome capture to infer the past demographic history of Norway spruce (P. abies) and estimate the amount of recent introduction used to establish the Norway spruce breeding program in southern Sweden. Most of these trees belong to P. abies and originate from the base populations of the Swedish breeding program. Others originate from populations across the natural ranges of the three species. Of the 1,499 individuals stemming from the breeding program, a large proportion corresponds to recent introductions from mainland Europe. The split of P. omorika occurred 23 million years ago (mya), while the divergence between P. obovata and P. abies began 17.6 mya. Demographic inferences retrieved the same main clusters within P. abies than previous studies, that is, a vast northern domain ranging from Norway to central Russia, where the species is progressively replaced by Siberian spruce (P. obovata) and two smaller domains, an Alpine domain and a Carpathian one, but also revealed further subdivision and gene flow among clusters. The three main domains divergence was ancient (15 mya), and all three went through a bottleneck corresponding to the LGM. Approximately 17% of P. abies Nordic domain migrated from P. obovata ~103K years ago, when both species had much larger effective population sizes. Our analysis of genomewide polymorphism data thus revealed the complex demographic history of Picea genus in Western Europe and highlighted the importance of material transfer in Swedish breeding program.
Collapse
Affiliation(s)
- Jun Chen
- Department of Ecology and Genetics, Evolutionary Biology CentreUppsala UniversityUppsalaSweden
| | - Lili Li
- Department of Ecology and Genetics, Evolutionary Biology CentreUppsala UniversityUppsalaSweden
| | - Pascal Milesi
- Department of Ecology and Genetics, Evolutionary Biology CentreUppsala UniversityUppsalaSweden
| | - Gunnar Jansson
- Forestry Research Institute of Sweden (Skogforsk)UppsalaSweden
| | - Mats Berlin
- Forestry Research Institute of Sweden (Skogforsk)UppsalaSweden
| | - Bo Karlsson
- Forestry Research Institute of Sweden (Skogforsk)EkeboSweden
| | - Jelena Aleksic
- Institute of Molecular Genetics and Genetic EngineeringUniversity of BelgradeBelgradeSerbia
| | - Giovanni G. Vendramin
- Division of Florence, Institute of Biosciences and BioResourcesNational Research Council (IBBR‐CNR)Sesto FiorentinoItaly
| | - Martin Lascoux
- Department of Ecology and Genetics, Evolutionary Biology CentreUppsala UniversityUppsalaSweden
| |
Collapse
|