1
|
Roberts M, Josephs EB. Previously unmeasured genetic diversity explains part of Lewontin's paradox in a k-mer-based meta-analysis of 112 plant species. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.17.594778. [PMID: 38798362 PMCID: PMC11118579 DOI: 10.1101/2024.05.17.594778] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
At the molecular level, most evolution is expected to be neutral. A key prediction of this expectation is that the level of genetic diversity in a population should scale with population size. However, as was noted by Richard Lewontin in 1974 and reaffirmed by later studies, the slope of the population size-diversity relationship in nature is much weaker than expected under neutral theory. We hypothesize that one contributor to this paradox is that current methods relying on single nucleotide polymorphisms (SNPs) called from aligning short reads to a reference genome underestimate levels of genetic diversity in many species. To test this idea, we calculated nucleotide diversity ( π ) and k-mer-based metrics of genetic diversity across 112 plant species, amounting to over 205 terabases of DNA sequencing data from 27,488 individual plants. We then compared how these different metrics correlated with proxies of population size that account for both range size and population density variation across species. We found that our population size proxies scaled anywhere from about 3 to over 20 times faster with k-mer diversity than nucleotide diversity after adjusting for evolutionary history, mating system, life cycle habit, cultivation status, and invasiveness. The relationship between k-mer diversity and population size proxies also remains significant after correcting for genome size, whereas the analogous relationship for nucleotide diversity does not. These results suggest that variation not captured by common SNP-based analyses explains part of Lewontin's paradox in plants.
Collapse
|
2
|
Shi T, Zhang X, Hou Y, Jia C, Dan X, Zhang Y, Jiang Y, Lai Q, Feng J, Feng J, Ma T, Wu J, Liu S, Zhang L, Long Z, Chen L, Street NR, Ingvarsson PK, Liu J, Yin T, Wang J. The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees. MOLECULAR PLANT 2024; 17:725-746. [PMID: 38486452 DOI: 10.1016/j.molp.2024.03.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Revised: 02/22/2024] [Accepted: 03/11/2024] [Indexed: 04/05/2024]
Abstract
Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.
Collapse
Affiliation(s)
- Tingting Shi
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Xinxin Zhang
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Yukang Hou
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Changfu Jia
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Xuming Dan
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Yulin Zhang
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Yuanzhong Jiang
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Qiang Lai
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Jiajun Feng
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Jianju Feng
- College of Horticulture and Forestry, Tarim University, Alar 843300, China
| | - Tao Ma
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Jiali Wu
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Shuyu Liu
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Lei Zhang
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Zhiqin Long
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Liyang Chen
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Nathaniel R Street
- Umeå Plant Science Centre, Department of Plant Physiology, Umeå University, Umeå, Västerbotten, Sweden
| | - Pär K Ingvarsson
- Linnean Centre for Plant Biology, Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Jianquan Liu
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China.
| | - Tongming Yin
- The Key Laboratory of Tree Genetics and Biotechnology of Jiangsu Province and Education Department of China, Nanjing Forestry University, Nanjing, Jiangsu, China.
| | - Jing Wang
- Key Laboratory for Bio-Resources and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China.
| |
Collapse
|
3
|
Zhao W, Gao J, Hall D, Andersson BA, Bruxaux J, Tomlinson KW, Drouzas AD, Suyama Y, Wang XR. Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow. THE NEW PHYTOLOGIST 2024. [PMID: 38515228 DOI: 10.1111/nph.19694] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Accepted: 03/04/2024] [Indexed: 03/23/2024]
Abstract
Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers withc . $$ c. $$ 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.
Collapse
Affiliation(s)
- Wei Zhao
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
| | - Jie Gao
- CAS Key Laboratory of Tropical Forest Ecology, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Yunnan, 666303, China
| | - David Hall
- Forestry Research Institute of Sweden (Skogforsk), Sävar, SE-91833, Sweden
| | - Bea Angelica Andersson
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
| | - Jade Bruxaux
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
| | - Kyle W Tomlinson
- Center for Integrative Conservation & Yunnan Key Laboratory for Conservation of Tropical Rainforests and Asian Elephant, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Yunnan, 666303, China
| | - Andreas D Drouzas
- Laboratory of Systematic Botany and Phytogeography, School of Biology, Aristotle University of Thessaloniki, Thessaloniki, 54124, Greece
| | - Yoshihisa Suyama
- Graduate School of Agricultural Science, Tohoku University, Miyagi, 989-6711, Japan
| | - Xiao-Ru Wang
- Department of Ecology and Environmental Science, Umeå Plant Science Center, Umeå University, Umeå, SE-90187, Sweden
- National Engineering Research Center of Tree Breeding and Ecological Restoration, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, 100083, China
| |
Collapse
|
4
|
Riehl JFL, Cole CT, Morrow CJ, Barker HL, Bernhardsson C, Rubert‐Nason K, Ingvarsson PK, Lindroth RL. Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen ( Populus tremuloides Michx.). Ecol Evol 2023; 13:e10541. [PMID: 37780087 PMCID: PMC10534199 DOI: 10.1002/ece3.10541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Revised: 08/30/2023] [Accepted: 09/04/2023] [Indexed: 10/03/2023] Open
Abstract
Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.
Collapse
Affiliation(s)
| | | | - Clay J. Morrow
- Department of Forest and Wildlife EcologyUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
| | - Hilary L. Barker
- Department of EntomologyUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
- Present address:
Office of Student SuccessWisconsin Technical College SystemMadisonWisconsinUSA
| | - Carolina Bernhardsson
- Department of Ecology and Environmental ScienceUmeå UniversityUmeåSweden
- Present address:
Department of Organismal Biology, Center for Evolutionary BiologyUppsala UniversityUppsalaSweden
| | - Kennedy Rubert‐Nason
- Department of EntomologyUniversity of Wisconsin‐MadisonMadisonWisconsinUSA
- Present address:
Division of Natural SciencesUniversity of Maine at Fort KentFort KentMaineUSA
| | - Pär K. Ingvarsson
- Department of Plant BiologySwedish University of Agricultural Sciences, Uppsala BioCenterUppsalaSweden
| | | |
Collapse
|
5
|
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
|
6
|
Yan W, Wang Z, Zhou B. Population evolution of seagrasses returning to the ocean. Heliyon 2023; 9:e20231. [PMID: 37809433 PMCID: PMC10559988 DOI: 10.1016/j.heliyon.2023.e20231] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 08/05/2023] [Accepted: 09/14/2023] [Indexed: 10/10/2023] Open
Abstract
Seagrasses are higher flowering plants that live entirely in marine environments, with the greatest habitat variation occurring from land to sea. Genetic structure or population differentiation history is a hot topic in evolutionary biology, which is of great significance for understanding speciation. Genetic information is obtained from geographically distributed subpopulations, different subspecies, or strains of the same species using next-generation sequencing techniques. Genetic variation is identified by comparison with reference genomes. Genetic diversity is explored using population structure, principal component analysis (PCA), and phylogenetic relationships. Patterns of population genetic differentiation are elucidated by combining the isolation by distance (IBD) model, linkage disequilibrium levels, and genetic statistical analysis. Demographic history is simulated using effective population size, divergence time, and site frequency spectrum (SFS). Through various population genetic analyses, the genetic structure and historical population dynamics of seagrass can be clarified, and their evolutionary processes can be further explored at the molecular level to understand how evolutionary processes contributed to the formation of early ecological species and provide data support for seagrass conservation.
Collapse
Affiliation(s)
- Wenjie Yan
- Key Laboratory of Mariculture (Ocean University of China), Ministry of Education, Qingdao, 266003, China
| | - Zhaohua Wang
- First Institute of Oceanography, MNR, Qingdao, 266061, China
| | - Bin Zhou
- College of Marine Life Science, Ocean University of China, Qingdao, 266003, China
| |
Collapse
|
7
|
Zhang H, Zhang X, Wu G, Dong C, Liu J, Li M. Genomic divergence and introgression among three Populus species. Mol Phylogenet Evol 2023; 180:107686. [PMID: 36586545 DOI: 10.1016/j.ympev.2022.107686] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Revised: 12/15/2022] [Accepted: 12/22/2022] [Indexed: 12/29/2022]
Abstract
Genomic divergence with gene flow is very common in both plants and animals. However, divergence and gene flow are two counteracting factors during speciation. Identifying the types of genes that are likely to be introgressed and what genetic factors restrict further effective reproduction of interspecific hybrids is of great interest to biologists. We aimed to address these issues using three related tree species, Populus alba (Pa), P. tremula (Pt), and P. tremuloides (Ps), and the interspecific hybrid of the former two species, P. × canescens (Pc). We collected 105 genomes for these four poplar lineages, including 28 Pa, 38Pt, 21 Ps, and 18 Pc individuals, to reconstruct their evolutionary histories. Our coalescence-based simulations indicated that Pa diverged earliest from Ps and Pt, and asymmetrical gene flow existed between any two lineages, with especially large ancient gene flow occurring between Pa and Pt. The genomic landscape of divergence between pairs of the three species are highly heterogeneous, which may have arisen through both divergent sorting of ancient polymorphisms and ongoing gene flow. We found that extant regions of the genome with introgressed ancestry reduced genetic divergence but elevated recombination rates and accounted for 5.76 % of the total genome. Introgressed genes were functionally associated with stress resistance, including innate immune response, anti-adversity response, and programmed cell death. However, candidate genes underlying postmating barriers of Pc were homozygous and resistant to introgression due to the incompatibility of alleles between loci after hybridization and were associated with endosperm and gamete formation and disease resistance. Our study revealed genomic dynamics during speciation with gene flow and identified regions of the genome that were likely introgressed and adaptive as well as candidate loci responsible for hybrid incompatibility that resulted in the formation of postmating barriers after hybridization.
Collapse
Affiliation(s)
- Han Zhang
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, China
| | - Xu Zhang
- Internet Medical and System Applications of National Engineering Laboratory, Zhengzhou University First Affiliated Hospital, Zhengzhou 450000, China
| | - Guili Wu
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, China
| | - Congcong Dong
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, China
| | - Jianquan Liu
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, China; Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu 610065, China
| | - Minjie Li
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, China.
| |
Collapse
|
8
|
Abeyratne CR, Macaya-Sanz D, Zhou R, Barry KW, Daum C, Haiby K, Lipzen A, Stanton B, Yoshinaga Y, Zane M, Tuskan GA, DiFazio SP. High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa. G3 (BETHESDA, MD.) 2023; 13:jkac269. [PMID: 36250890 PMCID: PMC9836356 DOI: 10.1093/g3journal/jkac269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Accepted: 09/28/2022] [Indexed: 12/14/2022]
Abstract
Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.
Collapse
Affiliation(s)
| | - David Macaya-Sanz
- Department of Forest Ecology & Genetics, CIFOR-INIA, CSIC, Madrid 28040, Spain
| | - Ran Zhou
- Warnell School of Forestry and Natural Resources, Department of Genetics, and Department of Plant Biology, University of Georgia, Athens, GA 30602, USA
| | - Kerrie W Barry
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | - Christopher Daum
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | | | - Anna Lipzen
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | | | - Yuko Yoshinaga
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | - Matthew Zane
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | - Gerald A Tuskan
- Biosciences Division, Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | - Stephen P DiFazio
- Department of Biology, West Virginia University, Morgantown, WV 26506, USA
| |
Collapse
|
9
|
Moreira LR, Klicka J, Smith BT. Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age. Mol Ecol 2023; 32:1739-1759. [PMID: 36617622 DOI: 10.1111/mec.16841] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Revised: 12/13/2022] [Accepted: 12/20/2022] [Indexed: 01/10/2023]
Abstract
The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.
Collapse
Affiliation(s)
- Lucas R Moreira
- Department of Ecology, Evolution, and Environmental Biology, Columbia University, New York, New York, USA.,Department of Ornithology, American Museum of Natural History, New York City, New York, USA.,Program in Bioinformatics and Integrative Biology, University of Massachusetts Chan Medical School, Worcester, Massachusetts, USA
| | - John Klicka
- Burke Museum of Natural History and Culture and Department of Biology, University of Washington, Seattle, Washington, USA
| | - Brian Tilston Smith
- Department of Ornithology, American Museum of Natural History, New York City, New York, USA
| |
Collapse
|
10
|
Rodríguez BD, Kloth KJ, Albrectsen BR. Effects of condensed tannins on behavior and performance of a specialist aphid on aspen. Ecol Evol 2022; 12:e9229. [PMID: 36016819 PMCID: PMC9396707 DOI: 10.1002/ece3.9229] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2021] [Revised: 07/05/2022] [Accepted: 08/01/2022] [Indexed: 11/18/2022] Open
Abstract
Genes involved in plant defences against herbivores and pathogens are often highly polymorphic. This is a putative sign that balancing selection may have operated reciprocally on the hosts and their herbivores. Spatial and temporal variations (for example, in soil nutrients and the plants' ontogenetic development) may also modulate resistance traits, and thus selection pressures, but have been largely overlooked in theories of plant defences. Important elements of defences in Populus tremula (hereafter aspen) are phenolic compounds, including condensed tannins (CTs). Concentrations of CTs vary considerably with both variations in external factors and time, but they are also believed to provide genotype-dependent resistance, mainly against chewing herbivores and pathogens. However, evidence of their contributions to resistance is sparse. Detailed studies of co-evolved plant-herbivore associations could provide valuable insights into these contributions. Therefore, we examined correlations between CT levels in aspen leaves and both the feeding behavior and reproduction of the specialist aspen leaf aphid (Chaitophorus tremulae) in varied conditions. We found that xylem sap intake and probing difficulties were higher on genotypes with high-CT concentrations. However, aphids engaged in more nonprobing activities on low-CT genotypes, indicating that CTs were not the only defence traits involved. Thus, high-CT genotypes were not necessarily more resistant than low-CT genotypes, but aphid reproduction was generally negatively correlated with local CT accumulation. Genotype-specific resistance ranking also depended on the experimental conditions. These results support the hypothesis that growth conditions may affect selection pressures mediated by aphids in accordance with balancing selection theory.
Collapse
Affiliation(s)
- Bárbara Díez Rodríguez
- Department of Plant Ecology and GeobotanyPhilipps‐University MarburgMarburgGermany
- Department of Plant PhysiologyUmeå Plant Science CentreUmeåSweden
| | - Karen J. Kloth
- Laboratory of EntomologyWageningen University and ResearchWageningenThe Netherlands
| | | |
Collapse
|
11
|
Sandercock AM, Westbrook JW, Zhang Q, Johnson HA, Saielli TM, Scrivani JA, Fitzsimmons SF, Collins K, Perkins MT, Craddock JH, Schmutz J, Grimwood J, Holliday JA. Frozen in time: rangewide genomic diversity, structure, and demographic history of relict American chestnut populations. Mol Ecol 2022; 31:4640-4655. [PMID: 35880415 DOI: 10.1111/mec.16629] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Revised: 07/15/2022] [Accepted: 07/18/2022] [Indexed: 11/30/2022]
Abstract
American chestnut (Castanea dentata) was once the most economically and ecologically important hardwood species in the eastern United States. In the first half of the 20th century, an exotic fungal pathogen - Cryphonectria parasitica - decimated the species, killing billions of chestnut trees. Two approaches to developing blight resistant American chestnut populations show promise, but both will require introduction of adaptive genomic diversity from wild germplasm to produce diverse, locally adapted restoration populations. Here we characterize population structure, demographic history, and genomic diversity in a range-wide sample of 384 wild American chestnuts to inform conservation and breeding with blight resistant varieties. Population structure analyses suggest that the chestnut range can be roughly divided into northeast, central, and southwest populations. Within-population genomic diversity estimates revealed a clinal pattern with the highest diversity in the southwest, which likely reflects bottleneck events associated with Quaternary glaciation. Finally, we identified genomic regions under positive selection within each population, which suggests that defense against fungal pathogens is a common target of selection across all populations. Taken together, these results show that American chestnut underwent a postglacial expansion from the southern portion of its range leading to three extant genetic populations. These populations will serve as management units for breeding adaptive genetic variation into the blight-resistant tree populations for targeted reintroduction efforts.
Collapse
Affiliation(s)
| | | | - Qian Zhang
- Department of Forest Resources and Environmental Conservation, Virginia Tech, Blacksburg, VA, USA
| | - Hayley A Johnson
- Department of Forest Resources and Environmental Conservation, Virginia Tech, Blacksburg, VA, USA
| | | | | | | | | | - M Taylor Perkins
- Department of Biology, Geology, and Environmental Science, The University of Tennessee at Chattanooga, Chattanooga, TN, USA
| | - J Hill Craddock
- Department of Biology, Geology, and Environmental Science, The University of Tennessee at Chattanooga, Chattanooga, TN, USA
| | - Jeremy Schmutz
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Jane Grimwood
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Jason A Holliday
- Department of Forest Resources and Environmental Conservation, Virginia Tech, Blacksburg, VA, USA
| |
Collapse
|
12
|
Booker TR, Payseur BA, Tigano A. Background selection under evolving recombination rates. Proc Biol Sci 2022; 289:20220782. [PMID: 35730151 PMCID: PMC9233929 DOI: 10.1098/rspb.2022.0782] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Background selection (BGS), the effect that purifying selection exerts on sites linked to deleterious alleles, is expected to be ubiquitous across eukaryotic genomes. The effects of BGS reflect the interplay of the rates and fitness effects of deleterious mutations with recombination. A fundamental assumption of BGS models is that recombination rates are invariant over time. However, in some lineages, recombination rates evolve rapidly, violating this central assumption. Here, we investigate how recombination rate evolution affects genetic variation under BGS. We show that recombination rate evolution modifies the effects of BGS in a manner similar to a localized change in the effective population size, potentially leading to underestimation or overestimation of the genome-wide effects of selection. Furthermore, we find evidence that recombination rate evolution in the ancestors of modern house mice may have impacted inferences of the genome-wide effects of selection in that species.
Collapse
Affiliation(s)
- Tom R. Booker
- Department of Zoology, University of British Columbia, Vancouver Campus, Vancouver, BC, Canada
| | - Bret A. Payseur
- Laboratory of Genetics, University of Wisconsin - Madison, Madison, WI, USA
| | - Anna Tigano
- Department of Biology, University of British Columbia, Okanagan Campus, Kelowna, BC, Canada
| |
Collapse
|
13
|
Li Z, Zhao W, Zhang J, Pan Z, Bai S, Tong C. A Novel Strategy to Reveal the Landscape of Crossovers in an F1 Hybrid Population of Populus deltoides and Populus simonii. PLANTS 2022; 11:plants11081046. [PMID: 35448774 PMCID: PMC9025136 DOI: 10.3390/plants11081046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Revised: 04/03/2022] [Accepted: 04/07/2022] [Indexed: 11/16/2022]
Abstract
Although the crossover (CO) patterns of different species have been extensively investigated, little is known about the landscape of CO patterns in Populus because of its high heterozygosity and long-time generation. A novel strategy was proposed to reveal the difference of CO rate and interference between Populus deltoides and Populus simonii using their F1 hybrid population. We chose restriction site-associated DNA (RAD) tags that contained two SNPs, one only receiving the CO information from the female P. deltoides and the other from the male P. simonii. These RAD tags allowed us to investigate the CO patterns between the two outbred species, instead of using the traditional backcross populations in inbred lines. We found that the CO rate in P. deltoides was generally greater than that in P. simonii, and that the CO interference was a common phenomenon across the two genomes. The COs landscape of the different Populus species facilitates not only to understand the evolutionary mechanism for adaptability but also to rebuild the statistical model for precisely constructing genetic linkage maps that are critical in genome assembly in Populus. Additionally, the novel strategy could be applied in other outbred species for investigating the CO patterns.
Collapse
Affiliation(s)
| | | | | | | | | | - Chunfa Tong
- Correspondence: ; Tel.: +86-025-85428817 (ext. 815)
| |
Collapse
|
14
|
Liang YY, Shi Y, Yuan S, Zhou BF, Chen XY, An QQ, Ingvarsson PK, Plomion C, Wang B. Linked selection shapes the landscape of genomic variation in three oak species. THE NEW PHYTOLOGIST 2022; 233:555-568. [PMID: 34637540 DOI: 10.1111/nph.17793] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2021] [Accepted: 09/27/2021] [Indexed: 06/13/2023]
Abstract
Natural selection shapes genome-wide patterns of diversity within species and divergence between species. However, quantifying the efficacy of selection and elucidating the relative importance of different types of selection in shaping genomic variation remain challenging. We sequenced whole genomes of 101 individuals of three closely related oak species to track the divergence history, and to dissect the impacts of selective sweeps and background selection on patterns of genomic variation. We estimated that the three species diverged around the late Neogene and experienced a bottleneck during the Pleistocene. We detected genomic regions with elevated relative differentiation ('FST -islands'). Population genetic inferences from the site frequency spectrum and ancestral recombination graph indicated that FST -islands were formed by selective sweeps. We also found extensive positive selection; the fixation of adaptive mutations and reduction neutral diversity around substitutions generated a signature of selective sweeps. Prevalent negative selection and background selection have reduced genetic diversity in both genic and intergenic regions, and contributed substantially to the baseline variation in genetic diversity. Our results demonstrate the importance of linked selection in shaping genomic variation, and illustrate how the extent and strength of different selection models vary across the genome.
Collapse
Affiliation(s)
- Yi-Ye Liang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China
- University of the Chinese Academy of Sciences, Beijing, 100049, China
| | - Yong Shi
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China
| | - Shuai Yuan
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China
| | - Biao-Feng Zhou
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China
| | - Xue-Yan Chen
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China
| | - Qing-Qing An
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China
| | - Pär K Ingvarsson
- Department of Plant Biology, Linnean Center for Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Sciences, Uppsala, SE-75007, Sweden
| | | | - Baosheng Wang
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China
- Center of Conservation Biology, Core Botanical Gardens, Chinese Academy of Sciences, Guangzhou, 510650, China
| |
Collapse
|
15
|
Population structure and adaptive variation of Helichrysum italicum (Roth) G. Don along eastern Adriatic temperature and precipitation gradient. Sci Rep 2021; 11:24333. [PMID: 34934087 PMCID: PMC8692458 DOI: 10.1038/s41598-021-03548-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Accepted: 12/06/2021] [Indexed: 01/07/2023] Open
Abstract
Immortelle (Helichrysum italicum (Roth) G. Don; Asteraceae) is a perennial plant species native to the Mediterranean region, known for many properties with wide application mainly in perfume and cosmetic industry. A total of 18 wild H. italicum populations systematically sampled along the eastern Adriatic environmental gradient were studied using AFLP markers to determine genetic diversity and structure and to identify loci potentially responsible for adaptive divergence. Results showed higher levels of intrapopulation diversity than interpopulation diversity. Genetic differentiation among populations was significant but low, indicating extensive gene flow between populations. Bayesian analysis of population structure revealed the existence of two genetic clusters. Combining the results of FST - outlier analysis (Mcheza and BayeScan) and genome-environment association analysis (Samβada, LFMM) four AFLP loci strongly associated with the bioclimatic variables Bio03 Isothermality, Bio08 Mean temperature of the wettest quarter, Bio15 Precipitation seasonality, and Bio17 Precipitation of driest quarter were found to be the main variables driving potential adaptive genetic variation in H. italicum along the eastern Adriatic environmental gradient. Redundancy analysis revealed that the partitioning of genetic variation was mainly associated with the adaptation to temperature oscillations. The results of the research may contribute to a clearer understanding of the importance of local adaptations for the genetic differentiation of Mediterranean plants and allow the planning of appropriate conservation strategies. However, considering that the identified outlier loci may be linked to genes under selection rather than being the target of natural selection, future studies must aim at their additional analysis.
Collapse
|
16
|
Lucena-Perez M, Kleinman-Ruiz D, Marmesat E, Saveljev AP, Schmidt K, Godoy JA. Bottleneck-associated changes in the genomic landscape of genetic diversity in wild lynx populations. Evol Appl 2021; 14:2664-2679. [PMID: 34815746 PMCID: PMC8591332 DOI: 10.1111/eva.13302] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Revised: 08/17/2021] [Accepted: 09/08/2021] [Indexed: 01/06/2023] Open
Abstract
Demographic bottlenecks generally reduce genetic diversity through more intense genetic drift, but their net effect may vary along the genome due to the random nature of genetic drift and to local effects of recombination, mutation, and selection. Here, we analyzed the changes in genetic diversity following a bottleneck by comparing whole-genome diversity patterns in populations with and without severe recent documented declines of Iberian (Lynx pardinus, n = 31) and Eurasian lynx (Lynx lynx, n = 29). As expected, overall genomic diversity correlated negatively with bottleneck intensity and/or duration. Correlations of genetic diversity with divergence, chromosome size, gene or functional site content, GC content, or recombination were observed in nonbottlenecked populations, but were weaker in bottlenecked populations. Also, functional features under intense purifying selection and the X chromosome showed an increase in the observed density of variants, even resulting in higher θ W diversity than in nonbottlenecked populations. Increased diversity seems to be related to both a higher mutational input in those regions creating a large collection of low-frequency variants, a few of which increase in frequency during the bottleneck to the point they become detectable with our limited sample, and the reduced efficacy of purifying selection, which affects not only protein structure and function but also the regulation of gene expression. The results of this study alert to the possible reduction of fitness and adaptive potential associated with the genomic erosion in regulatory elements. Further, the detection of a gain of diversity in ultra-conserved elements can be used as a sensitive and easy-to-apply signature of genetic erosion in wild populations.
Collapse
Affiliation(s)
- Maria Lucena-Perez
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
| | - Daniel Kleinman-Ruiz
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
- Departamento de Genética Facultad de Biología Universidad Complutense Madrid Spain
| | - Elena Marmesat
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
| | - Alexander P Saveljev
- Department of Animal Ecology Russian Research Institute of Game Management and Fur Farming Kirov Russia
| | - Krzysztof Schmidt
- Mammal Research Institute Polish Academy of Sciences Białowieża Poland
| | - José A Godoy
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
| |
Collapse
|
17
|
Janzen T, Miró Pina V. Estimating the time since admixture from phased and unphased molecular data. Mol Ecol Resour 2021; 22:908-926. [PMID: 34599646 PMCID: PMC9291888 DOI: 10.1111/1755-0998.13519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Revised: 09/20/2021] [Accepted: 09/22/2021] [Indexed: 11/26/2022]
Abstract
After admixture, recombination breaks down genomic blocks of contiguous ancestry. The breakdown of these blocks forms a new “molecular clock” that ticks at a much faster rate than the mutation clock, enabling accurate dating of admixture events in the recent past. However, existing theory on the breakdown of these blocks, or the accumulation of delineations between blocks, so‐called “junctions”, has mostly been limited to using regularly spaced markers on phased data. Here, we present an extension to the theory of junctions using the ancestral recombination graph that describes the expected number of junctions for any distribution of markers along the genome. Furthermore, we provide a new framework to infer the time since admixture using unphased data. We demonstrate both the phased and unphased methods on simulated data and show that our new extensions have improved accuracy with respect to previous methods, especially for smaller population sizes and more ancient admixture times. Lastly, we demonstrate the applicability of our method on three empirical data sets, including labcrosses of yeast (Saccharomyces cerevisae) and two case studies of hybridization in swordtail fish and Populus trees.
Collapse
Affiliation(s)
- Thijs Janzen
- Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, The Netherlands.,Carl von Ossietzky University, Oldenburg, Germany
| | - Verónica Miró Pina
- Instituto de Investigaciones en Matemáticas Aplicadas y Sistemas (IIMAS), Universidad Nacional Autónoma de México (UNAM), México City, México.,Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Barcelona, Spain.,Universitat Pompeu Fabra (UPF), Barcelona, Spain
| |
Collapse
|
18
|
Danguy des Déserts A, Bouchet S, Sourdille P, Servin B. Evolution of Recombination Landscapes in Diverging Populations of Bread Wheat. Genome Biol Evol 2021; 13:evab152. [PMID: 34185074 PMCID: PMC8350361 DOI: 10.1093/gbe/evab152] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/24/2021] [Indexed: 11/16/2022] Open
Abstract
Reciprocal exchanges of DNA (crossovers) that occur during meiosis are mandatory to ensure the production of fertile gametes in sexually reproducing species. They also contribute to shuffle parental alleles into new combinations thereby fueling genetic variation and evolution. However, due to biological constraints, the recombination landscape is highly heterogeneous along the genome which limits the range of allelic combinations and the adaptability of populations. An approach to better understand the constraints on the recombination process is to study how it evolved in the past. In this work, we tackled this question by constructing recombination profiles in four diverging bread wheat (Triticum aestivum L.) populations established from 371 landraces genotyped at 200,062 SNPs. We used linkage disequilibrium (LD) patterns to estimate in each population the past distribution of recombination along the genome and characterize its fine-scale heterogeneity. At the megabase scale, recombination rates derived from LD patterns were consistent with family-based estimates obtained from a population of 406 recombinant inbred lines. Among the four populations, recombination landscapes were positively correlated between each other and shared a statistically significant proportion of highly recombinant intervals. However, this comparison also highlighted that the similarity in recombination landscapes between populations was significantly decreasing with their genetic differentiation in most regions of the genome. This observation was found to be robust to SNPs ascertainment and demography and suggests a relatively rapid evolution of factors determining the fine-scale localization of recombination in bread wheat.
Collapse
Affiliation(s)
- Alice Danguy des Déserts
- INRAE-Université Clermont-Auvergne, UMR1095, Génétique Diversité Ecophysiologie des Céréales, Clermont-Ferrand, France
| | - Sophie Bouchet
- INRAE-Université Clermont-Auvergne, UMR1095, Génétique Diversité Ecophysiologie des Céréales, Clermont-Ferrand, France
| | - Pierre Sourdille
- INRAE-Université Clermont-Auvergne, UMR1095, Génétique Diversité Ecophysiologie des Céréales, Clermont-Ferrand, France
| | - Bertrand Servin
- INRAE, Université de Toulouse, GenPhySE, Castanet-Tolosan, France
| |
Collapse
|
19
|
Apuli RP, Richards T, Rendón-Anaya M, Karacic A, Rönnberg-Wästljung AC, Ingvarsson PK. The genetic basis of adaptation in phenology in an introduced population of Black Cottonwood (Populus trichocarpa, Torr. & Gray). BMC PLANT BIOLOGY 2021; 21:317. [PMID: 34215191 PMCID: PMC8252265 DOI: 10.1186/s12870-021-03103-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Accepted: 06/16/2021] [Indexed: 06/13/2023]
Abstract
BACKGROUND Entering and exiting winter dormancy present important trade-offs between growth and survival at northern latitudes. Many forest trees display local adaptation across latitude in traits associated with these phenology transitions. Transfers of a species outside its native range introduce the species to novel combinations of environmental conditions potentially requiring different combinations of alleles to optimize growth and survival. In this study, we performed genome wide association analyses and a selection scan in a P. trichocarpa mapping population derived from crossings between clones collected across the native range and introduced into Sweden. GWAS analyses were performed using phenotypic data collected across two field seasons and in a controlled phytotron experiment. RESULTS We uncovered 584 putative candidate genes associated with spring and autumn phenology traits as well as with growth. Many regions harboring variation significantly associated with the initiation of leaf shed and leaf autumn coloring appeared to have been evolving under positive selection in the native environments of P. trichocarpa. A comparison between the candidate genes identified with results from earlier GWAS analyses performed in the native environment found a smaller overlap for spring phenology traits than for autumn phenology traits, aligning well with earlier observations that spring phenology transitions have a more complex genetic basis than autumn phenology transitions. CONCLUSIONS In a small and structured introduced population of P. trichocarpa, we find complex genetic architectures underlying all phenology and growth traits, and identify multiple putative candidate genes despite the limitations of the study population.
Collapse
Affiliation(s)
- Rami-Petteri Apuli
- Linnean Centre for Plant Biology, Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Science, Uppsala, Sweden
| | - Thomas Richards
- Linnean Centre for Plant Biology, Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Science, Uppsala, Sweden
- Plant Ecology and Evolution, Department of Ecology and Genetics, Evolutionary Biology Centre, Uppsala University, Uppsala, Sweden
| | - Martha Rendón-Anaya
- Linnean Centre for Plant Biology, Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Science, Uppsala, Sweden
| | - Almir Karacic
- Institute for Crop Production Ecology, Swedish University of Agricultural Science, Uppsala, Sweden
| | - Ann-Christin Rönnberg-Wästljung
- Linnean Centre for Plant Biology, Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Science, Uppsala, Sweden
| | - Pär K Ingvarsson
- Linnean Centre for Plant Biology, Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Science, Uppsala, Sweden.
| |
Collapse
|
20
|
Schield DR, Pasquesi GIM, Perry BW, Adams RH, Nikolakis ZL, Westfall AK, Orton RW, Meik JM, Mackessy SP, Castoe TA. Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9. Mol Biol Evol 2021; 37:1272-1294. [PMID: 31926008 DOI: 10.1093/molbev/msaa003] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.
Collapse
Affiliation(s)
- Drew R Schield
- Department of Biology, University of Texas at Arlington, Arlington, TX
| | | | - Blair W Perry
- Department of Biology, University of Texas at Arlington, Arlington, TX
| | - Richard H Adams
- Department of Biology, University of Texas at Arlington, Arlington, TX.,Department of Computer & Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, FL
| | | | | | - Richard W Orton
- Department of Biology, University of Texas at Arlington, Arlington, TX
| | - Jesse M Meik
- Department of Biological Sciences, Tarleton State University, Stephenville, TX
| | - Stephen P Mackessy
- School of Biological Sciences, University of Northern Colorado, Greeley, CO
| | - Todd A Castoe
- Department of Biology, University of Texas at Arlington, Arlington, TX
| |
Collapse
|
21
|
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
|
22
|
Sanderson BJ, DiFazio SP, Cronk QCB, Ma T, Olson MS. A targeted sequence capture array for phylogenetics and population genomics in the Salicaceae. APPLICATIONS IN PLANT SCIENCES 2020; 8:e11394. [PMID: 33163293 PMCID: PMC7598885 DOI: 10.1002/aps3.11394] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Accepted: 08/12/2020] [Indexed: 05/03/2023]
Abstract
PREMISE The family Salicaceae has proved taxonomically challenging, especially in the genus Salix, which is speciose and features frequent hybridization and polyploidy. Past efforts to reconstruct the phylogeny with molecular barcodes have failed to resolve the species relationships of many sections of the genus. METHODS We used the wealth of sequence data in the family to design sequence capture probes to target regions of 300-1200 bp of exonic regions of 972 genes. RESULTS We recovered sequence data for nearly all of the targeted genes in three species of Populus and three species of Salix. We present a species tree, discuss concordance among gene trees, and present population genomic summary statistics for these loci. CONCLUSIONS Our sequence capture array has extremely high capture efficiency within the genera Populus and Salix, resulting in abundant phylogenetic information. Additionally, these loci show promise for population genomic studies.
Collapse
Affiliation(s)
- Brian J. Sanderson
- Department of Biological SciencesTexas Tech UniversityLubbockTexas79409‐3131USA
- Department of BiologyWest Virginia UniversityMorgantownWest Virginia26506USA
| | - Stephen P. DiFazio
- Department of BiologyWest Virginia UniversityMorgantownWest Virginia26506USA
| | - Quentin C. B. Cronk
- Department of BotanyUniversity of British ColumbiaVancouverBritish ColumbiaV6T 1Z4Canada
| | - Tao Ma
- Key Laboratory of Bio‐Resource and Eco‐Environment of Ministry of EducationCollege of Life SciencesSichuan UniversityChengdu610065People’s Republic of China
| | - Matthew S. Olson
- Department of Biological SciencesTexas Tech UniversityLubbockTexas79409‐3131USA
| |
Collapse
|
23
|
Zheng C, Tan L, Sang M, Ye M, Wu R. Genetic adaptation of Tibetan poplar ( Populus szechuanica var. tibetica) to high altitudes on the Qinghai-Tibetan Plateau. Ecol Evol 2020; 10:10974-10985. [PMID: 33144942 PMCID: PMC7593140 DOI: 10.1002/ece3.6508] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2019] [Revised: 05/14/2020] [Accepted: 05/28/2020] [Indexed: 12/26/2022] Open
Abstract
Plant adaptation to high altitudes has long been a substantial focus of ecological and evolutionary research. However, the genetic mechanisms underlying such adaptation remain poorly understood. Here, we address this issue by sampling, genotyping, and comparing populations of Tibetan poplar, Populus szechuanica var. tibetica, distributed from low (~2,000 m) to high altitudes (~3,000 m) of Sejila Mountain on the Qinghai-Tibet Plateau. Population structure analyses allow clear classification of two groups according to their altitudinal distributions. However, in contrast to the genetic variation within each population, differences between the two populations only explain a small portion of the total genetic variation (3.64%). We identified asymmetrical gene flow from high- to low-altitude populations. Integrating population genomic and landscape genomic analyses, we detected two hotspot regions, one containing four genes associated with altitudinal variation, and the other containing ten genes associated with response to solar radiation. These genes participate in abiotic stress resistance and regulation of reproductive processes. Our results provide insight into the genetic mechanisms underlying high-altitude adaptation in Tibetan poplar.
Collapse
Affiliation(s)
- Chenfei Zheng
- Beijing Advanced Innovation Center for Tree Breeding by Molecular DesignCenter for Computational BiologyCollege of Biological Sciences and TechnologyBeijing Forestry UniversityBeijingChina
| | - Lizhi Tan
- Beijing Advanced Innovation Center for Tree Breeding by Molecular DesignCenter for Computational BiologyCollege of Biological Sciences and TechnologyBeijing Forestry UniversityBeijingChina
| | - Mengmeng Sang
- Beijing Advanced Innovation Center for Tree Breeding by Molecular DesignCenter for Computational BiologyCollege of Biological Sciences and TechnologyBeijing Forestry UniversityBeijingChina
| | - Meixia Ye
- Beijing Advanced Innovation Center for Tree Breeding by Molecular DesignCenter for Computational BiologyCollege of Biological Sciences and TechnologyBeijing Forestry UniversityBeijingChina
| | - Rongling Wu
- Beijing Advanced Innovation Center for Tree Breeding by Molecular DesignCenter for Computational BiologyCollege of Biological Sciences and TechnologyBeijing Forestry UniversityBeijingChina
- Center for Statistical GeneticsPennsylvania State UniversityHersheyPAUSA
| |
Collapse
|
24
|
Khedikar Y, Clarke WE, Chen L, Higgins EE, Kagale S, Koh CS, Bennett R, Parkin IAP. Narrow genetic base shapes population structure and linkage disequilibrium in an industrial oilseed crop, Brassica carinata A. Braun. Sci Rep 2020; 10:12629. [PMID: 32724070 PMCID: PMC7387349 DOI: 10.1038/s41598-020-69255-w] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 07/09/2020] [Indexed: 12/16/2022] Open
Abstract
Ethiopian mustard (Brassica carinata A. Braun) is an emerging sustainable source of vegetable oil, in particular for the biofuel industry. The present study exploited genome assemblies of the Brassica diploids, Brassica nigra and Brassica oleracea, to discover over 10,000 genome-wide SNPs using genotype by sequencing of 620 B. carinata lines. The analyses revealed a SNP frequency of one every 91.7 kb, a heterozygosity level of 0.30, nucleotide diversity levels of 1.31 × 10-05, and the first five principal components captured only 13% molecular variation, indicating low levels of genetic diversity among the B. carinata collection. Genome bias was observed, with greater SNP density found on the B subgenome. The 620 lines clustered into two distinct sub-populations (SP1 and SP2) with the majority of accessions (88%) clustered in SP1 with those from Ethiopia, the presumed centre of origin. SP2 was distinguished by a collection of breeding lines, implicating targeted selection in creating population structure. Two selective sweep regions on B3 and B8 were detected, which harbour genes involved in fatty acid and aliphatic glucosinolate biosynthesis, respectively. The assessment of genetic diversity, population structure, and LD in the global B. carinata collection provides critical information to assist future crop improvement.
Collapse
Affiliation(s)
- Yogendra Khedikar
- Agriculture and Agri-Food Canada, 107 Science Place, Saskatoon, SK, Canada
| | - Wayne E Clarke
- Agriculture and Agri-Food Canada, 107 Science Place, Saskatoon, SK, Canada
| | - Lifeng Chen
- Agrisoma Biosciences Inc., 110 Gymnasium Place, Saskatoon, SK, Canada
| | - Erin E Higgins
- Agriculture and Agri-Food Canada, 107 Science Place, Saskatoon, SK, Canada
| | - Sateesh Kagale
- National Research Council Canada, 110 Gymnasium Place, Saskatoon, SK, Canada
| | - Chu Shin Koh
- Global Institute of Food Security, Saskatoon, SK, Canada
| | - Rick Bennett
- Agrisoma Biosciences Inc., 110 Gymnasium Place, Saskatoon, SK, Canada
| | - Isobel A P Parkin
- Agriculture and Agri-Food Canada, 107 Science Place, Saskatoon, SK, Canada.
| |
Collapse
|
25
|
Shang H, Hess J, Pickup M, Field DL, Ingvarsson PK, Liu J, Lexer C. Evolution of strong reproductive isolation in plants: broad-scale patterns and lessons from a perennial model group. Philos Trans R Soc Lond B Biol Sci 2020; 375:20190544. [PMID: 32654641 DOI: 10.1098/rstb.2019.0544] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Many recent studies have addressed the mechanisms operating during the early stages of speciation, but surprisingly few studies have tested theoretical predictions on the evolution of strong reproductive isolation (RI). To help address this gap, we first undertook a quantitative review of the hybrid zone literature for flowering plants in relation to reproductive barriers. Then, using Populus as an exemplary model group, we analysed genome-wide variation for phylogenetic tree topologies in both early- and late-stage speciation taxa to determine how these patterns may be related to the genomic architecture of RI. Our plant literature survey revealed variation in barrier complexity and an association between barrier number and introgressive gene flow. Focusing on Populus, our genome-wide analysis of tree topologies in speciating poplar taxa points to unusually complex genomic architectures of RI, consistent with earlier genome-wide association studies. These architectures appear to facilitate the 'escape' of introgressed genome segments from polygenic barriers even with strong RI, thus affecting their relationships with recombination rates. Placed within the context of the broader literature, our data illustrate how phylogenomic approaches hold great promise for addressing the evolution and temporary breakdown of RI during late stages of speciation. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.
Collapse
Affiliation(s)
- Huiying Shang
- Department of Botany and Biodiversity Research, University of Vienna, Rennweg 14, 1030 Vienna, Austria.,Vienna Graduate School of Population Genetics, Vienna, Austria
| | - Jaqueline Hess
- Department of Botany and Biodiversity Research, University of Vienna, Rennweg 14, 1030 Vienna, Austria.,Helmholtz Centre for Environmental Research, Halle (Saale), Germany
| | - Melinda Pickup
- Institute of Science and Technology (IST), Klosterneuburg, Austria
| | - David L Field
- Department of Botany and Biodiversity Research, University of Vienna, Rennweg 14, 1030 Vienna, Austria.,Edith Cowan University, Perth, Australia
| | - Pär K Ingvarsson
- Swedish University of Agricultural Sciences (SLU), Uppsala, Sweden
| | - Jianquan Liu
- Key Laboratory for Bio-resources and Eco-environment, College of Life Science, Sichuan University, Chengdu, People's Republic of China
| | - Christian Lexer
- Department of Botany and Biodiversity Research, University of Vienna, Rennweg 14, 1030 Vienna, Austria
| |
Collapse
|
26
|
Bagley JC, Heming NM, Gutiérrez EE, Devisetty UK, Mock KE, Eckert AJ, Strauss SH. Genotyping-by-sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen ( Populus tremuloides). Ecol Evol 2020; 10:4609-4629. [PMID: 32551047 PMCID: PMC7297775 DOI: 10.1002/ece3.6214] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2020] [Revised: 02/23/2020] [Accepted: 02/24/2020] [Indexed: 01/01/2023] Open
Abstract
Populus tremuloides is the widest-ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome-wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal-Cascades (cluster 1), east-slope Cascades-Sierra Nevadas-Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the "stable-edge hypothesis" that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited "trailing-edge" dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the "inland dispersal hypothesis" predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific-coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable-edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific-coastal genetic lineage of quaking aspen.
Collapse
Affiliation(s)
- Justin C. Bagley
- Plant Evolutionary Genomics LaboratoryDepartment of BiologyVirginia Commonwealth UniversityRichmondVAUSA
- Departamento de ZoologiaInstituto de Ciências BiológicasUniversidade de BrasíliaBrasíliaBrazil
| | - Neander M. Heming
- Departamento de ZoologiaInstituto de Ciências BiológicasUniversidade de BrasíliaBrasíliaBrazil
| | - Eliécer E. Gutiérrez
- Departamento de ZoologiaInstituto de Ciências BiológicasUniversidade de BrasíliaBrasíliaBrazil
- Programa de Pos‐Graduação em Biodiversidade AnimalCentro de Ciências Naturais e ExatasUniversidade Federal de Santa MariaSanta MariaBrazil
| | | | - Karen E. Mock
- Department of Wildland Resources and Ecology CenterUtah State UniversityLoganUTUSA
| | - Andrew J. Eckert
- Plant Evolutionary Genomics LaboratoryDepartment of BiologyVirginia Commonwealth UniversityRichmondVAUSA
| | - Steven H. Strauss
- Department of Forest Ecosystems and SocietyOregon State UniversityCorvallisORUSA
| |
Collapse
|
27
|
From molecules to populations: appreciating and estimating recombination rate variation. Nat Rev Genet 2020; 21:476-492. [DOI: 10.1038/s41576-020-0240-1] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/15/2020] [Indexed: 02/07/2023]
|
28
|
Wang J, Street NR, Park EJ, Liu J, Ingvarsson PK. Evidence for widespread selection in shaping the genomic landscape during speciation of Populus. Mol Ecol 2020; 29:1120-1136. [PMID: 32068935 DOI: 10.1111/mec.15388] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Revised: 02/13/2020] [Accepted: 02/14/2020] [Indexed: 12/13/2022]
Abstract
Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome-wide patterns of within- and between- species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole-genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long-term balancing selection have also been crucial components in shaping patterns of genome-wide variation during the speciation process.
Collapse
Affiliation(s)
- Jing Wang
- Key Laboratory for Bio-Resources and Eco-Environment, College of Life Science, Sichuan University, Chengdu, China
| | - Nathaniel R Street
- Department of Plant Physiology, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
| | - Eung-Jun Park
- Department of Bioresources, National Institute of Forest Science, Suwon, Korea
| | - Jianquan Liu
- Key Laboratory for Bio-Resources and Eco-Environment, College of Life Science, Sichuan University, Chengdu, China
| | - Pär K Ingvarsson
- Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
29
|
Inferring the Genomic Landscape of Recombination Rate Variation in European Aspen ( Populus tremula). G3-GENES GENOMES GENETICS 2020; 10:299-309. [PMID: 31744900 PMCID: PMC6945010 DOI: 10.1534/g3.119.400504] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
The rate of meiotic recombination is one of the central factors determining genome-wide levels of linkage disequilibrium which has important consequences for the efficiency of natural selection and for the dissection of quantitative traits. Here we present a new, high-resolution linkage map for Populus tremula that we use to anchor approximately two thirds of the P. tremula draft genome assembly on to the expected 19 chromosomes, providing us with the first chromosome-scale assembly for P. tremula (Table 2). We then use this resource to estimate variation in recombination rates across the P. tremula genome and compare these results to recombination rates based on linkage disequilibrium in a large number of unrelated individuals. We also assess how variation in recombination rates is associated with a number of genomic features, such as gene density, repeat density and methylation levels. We find that recombination rates obtained from the two methods largely agree, although the LD-based method identifies a number of genomic regions with very high recombination rates that the map-based method fails to detect. Linkage map and LD-based estimates of recombination rates are positively correlated and show similar correlations with other genomic features, showing that both methods can accurately infer recombination rate variation across the genome. Recombination rates are positively correlated with gene density and negatively correlated with repeat density and methylation levels, suggesting that recombination is largely directed toward gene regions in P. tremula.
Collapse
|
30
|
Peñalba JV, Deng Y, Fang Q, Joseph L, Moritz C, Cockburn A. Genome of an iconic Australian bird: High-quality assembly and linkage map of the superb fairy-wren (Malurus cyaneus). Mol Ecol Resour 2020; 20:560-578. [PMID: 31821695 DOI: 10.1111/1755-0998.13124] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Revised: 11/09/2019] [Accepted: 12/02/2019] [Indexed: 12/13/2022]
Abstract
The superb fairy-wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high-quality reference. To rectify that, we present the first high-quality genome assembly of a Meliphagides representative: the superb fairy-wren. We combined Illumina shotgun and mate-pair sequences, PacBio long-reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07-Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high-quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree-based bird genetic maps, we find that the fairy-wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high-quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy-wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.
Collapse
Affiliation(s)
- Joshua V Peñalba
- Division of Evolutionary Biology, Ludwig Maximilian University of Munich, Munich, Germany
| | | | - Qi Fang
- BGI-Shenzhen, Shenzhen, China
| | - Leo Joseph
- Australian National Wildlife Collection, CSIRO National Research Collections, Australia, Canberra, ACT, Australia
| | - Craig Moritz
- Centre for Biodiversity Analysis, Acton, ACT, Australia.,Division of Ecology and Evolution, Australian National University, Acton, ACT, Australia
| | - Andrew Cockburn
- Division of Ecology and Evolution, Australian National University, Acton, ACT, Australia
| |
Collapse
|
31
|
Guerra FP, Suren H, Holliday J, Richards JH, Fiehn O, Famula R, Stanton BJ, Shuren R, Sykes R, Davis MF, Neale DB. Exome resequencing and GWAS for growth, ecophysiology, and chemical and metabolomic composition of wood of Populus trichocarpa. BMC Genomics 2019; 20:875. [PMID: 31747881 PMCID: PMC6864938 DOI: 10.1186/s12864-019-6160-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2019] [Accepted: 10/09/2019] [Indexed: 12/26/2022] Open
Abstract
Background Populus trichocarpa is an important forest tree species for the generation of lignocellulosic ethanol. Understanding the genomic basis of biomass production and chemical composition of wood is fundamental in supporting genetic improvement programs. Considerable variation has been observed in this species for complex traits related to growth, phenology, ecophysiology and wood chemistry. Those traits are influenced by both polygenic control and environmental effects, and their genome architecture and regulation are only partially understood. Genome wide association studies (GWAS) represent an approach to advance that aim using thousands of single nucleotide polymorphisms (SNPs). Genotyping using exome capture methodologies represent an efficient approach to identify specific functional regions of genomes underlying phenotypic variation. Results We identified 813 K SNPs, which were utilized for genotyping 461 P. trichocarpa clones, representing 101 provenances collected from Oregon and Washington, and established in California. A GWAS performed on 20 traits, considering single SNP-marker tests identified a variable number of significant SNPs (p-value < 6.1479E-8) in association with diameter, height, leaf carbon and nitrogen contents, and δ15N. The number of significant SNPs ranged from 2 to 220 per trait. Additionally, multiple-marker analyses by sliding-windows tests detected between 6 and 192 significant windows for the analyzed traits. The significant SNPs resided within genes that encode proteins belonging to different functional classes as such protein synthesis, energy/metabolism and DNA/RNA metabolism, among others. Conclusions SNP-markers within genes associated with traits of importance for biomass production were detected. They contribute to characterize the genomic architecture of P. trichocarpa biomass required to support the development and application of marker breeding technologies.
Collapse
Affiliation(s)
- Fernando P Guerra
- Department of Plant Sciences, University of California at Davis, 262C Robbins Hall, Mail Stop 4, Davis, CA, 95616, USA.,Instituto de Ciencias Biológicas, Universidad de Talca, Talca, P.O. Box 747, 3460000, Chile
| | - Haktan Suren
- Department of Forest Resources and Environmental Conservation, Virginia Polytechnic Institute and State University, Blacksburg, VA, 24061, USA
| | - Jason Holliday
- Department of Forest Resources and Environmental Conservation, Virginia Polytechnic Institute and State University, Blacksburg, VA, 24061, USA
| | - James H Richards
- Department of Land, Air and Water Resources, University of California, Davis, CA, 95616, USA
| | - Oliver Fiehn
- Department of Molecular and Cellular Biology & Genome Center, University of California, Davis, CA, 95616, USA
| | - Randi Famula
- Department of Plant Sciences, University of California at Davis, 262C Robbins Hall, Mail Stop 4, Davis, CA, 95616, USA
| | - Brian J Stanton
- Biological Research Group, GreenWood Resources, Portland, OR, 97201, USA
| | - Richard Shuren
- Biological Research Group, GreenWood Resources, Portland, OR, 97201, USA
| | - Robert Sykes
- National Renewable Energy Laboratory, Golden, CO, 80401, USA
| | - Mark F Davis
- National Renewable Energy Laboratory, Golden, CO, 80401, USA
| | - David B Neale
- Department of Plant Sciences, University of California at Davis, 262C Robbins Hall, Mail Stop 4, Davis, CA, 95616, USA. .,Bioenergy Research Center, University of California at Davis, Davis, CA, 95616, USA.
| |
Collapse
|
32
|
V. Barroso G, Puzović N, Dutheil JY. Inference of recombination maps from a single pair of genomes and its application to ancient samples. PLoS Genet 2019; 15:e1008449. [PMID: 31725722 PMCID: PMC6879166 DOI: 10.1371/journal.pgen.1008449] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Revised: 11/26/2019] [Accepted: 09/30/2019] [Indexed: 12/11/2022] Open
Abstract
Understanding the causes and consequences of recombination landscape evolution is a fundamental goal in genetics that requires recombination maps from across the tree of life. Such maps can be obtained from population genomic datasets, but require large sample sizes. Alternative methods are therefore necessary to research organisms where such datasets cannot be generated easily, such as non-model or ancient species. Here we extend the sequentially Markovian coalescent model to jointly infer demography and the spatial variation in recombination rate. Using extensive simulations and sequence data from humans, fruit-flies and a fungal pathogen, we demonstrate that iSMC accurately infers recombination maps under a wide range of scenarios-remarkably, even from a single pair of unphased genomes. We exploit this possibility and reconstruct the recombination maps of ancient hominins. We report that the ancient and modern maps are correlated in a manner that reflects the established phylogeny of Neanderthals, Denisovans, and modern human populations.
Collapse
Affiliation(s)
- Gustavo V. Barroso
- Max Planck Institute for Evolutionary Biology, Department of Evolutionary Genetics, August-Thienemann-Straße , Plön–GERMANY
| | - Nataša Puzović
- Max Planck Institute for Evolutionary Biology, Department of Evolutionary Genetics, August-Thienemann-Straße , Plön–GERMANY
| | - Julien Y. Dutheil
- Max Planck Institute for Evolutionary Biology, Department of Evolutionary Genetics, August-Thienemann-Straße , Plön–GERMANY
| |
Collapse
|
33
|
Mattila TM, Laenen B, Horvath R, Hämälä T, Savolainen O, Slotte T. Impact of demography on linked selection in two outcrossing Brassicaceae species. Ecol Evol 2019; 9:9532-9545. [PMID: 31534673 PMCID: PMC6745670 DOI: 10.1002/ece3.5463] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 06/28/2019] [Accepted: 07/02/2019] [Indexed: 12/13/2022] Open
Abstract
Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.
Collapse
Affiliation(s)
- Tiina M. Mattila
- Department of Ecology and GeneticsUniversity of OuluOuluFinland
- Present address:
Department of Organismal BiologyUppsala UniversityUppsalaSweden
| | - Benjamin Laenen
- Science for Life Laboratory, Department of Ecology, Environment, and Plant SciencesStockholm UniversityStockholmSweden
| | - Robert Horvath
- Science for Life Laboratory, Department of Ecology, Environment, and Plant SciencesStockholm UniversityStockholmSweden
| | - Tuomas Hämälä
- Department of Ecology and GeneticsUniversity of OuluOuluFinland
- Biocenter OuluUniversity of OuluOuluFinland
- Present address:
Department of Plant and Microbial BiologyUniversity of Minnesota Twin CitiesSt. PaulMNUSA
| | - Outi Savolainen
- Department of Ecology and GeneticsUniversity of OuluOuluFinland
- Biocenter OuluUniversity of OuluOuluFinland
| | - Tanja Slotte
- Science for Life Laboratory, Department of Ecology, Environment, and Plant SciencesStockholm UniversityStockholmSweden
| |
Collapse
|
34
|
Shen C, Wang N, Huang C, Wang M, Zhang X, Lin Z. Population genomics reveals a fine-scale recombination landscape for genetic improvement of cotton. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2019; 99:494-505. [PMID: 31002209 DOI: 10.1111/tpj.14339] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2019] [Revised: 03/17/2019] [Accepted: 04/01/2019] [Indexed: 05/28/2023]
Abstract
Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population-scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population-scaled recombination rate between Gossypium hirsutum and G. arboreum and sub-genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome-wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high-recombination regions than were those of yield and early maturity traits. These results provide insights into the population-scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.
Collapse
Affiliation(s)
- Chao Shen
- National Key Laboratory of Crop Genetic Improvement, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei, China
| | - Nian Wang
- National Key Laboratory of Crop Genetic Improvement, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei, China
| | - Cong Huang
- National Key Laboratory of Crop Genetic Improvement, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei, China
| | - Maojun Wang
- National Key Laboratory of Crop Genetic Improvement, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei, China
| | - Xianlong Zhang
- National Key Laboratory of Crop Genetic Improvement, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei, China
| | - Zhongxu Lin
- National Key Laboratory of Crop Genetic Improvement, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei, China
| |
Collapse
|
35
|
Barker HL, Riehl JF, Bernhardsson C, Rubert-Nason KF, Holeski LM, Ingvarsson PK, Lindroth RL. Linking plant genes to insect communities: Identifying the genetic bases of plant traits and community composition. Mol Ecol 2019; 28:4404-4421. [PMID: 31233634 DOI: 10.1111/mec.15158] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 05/31/2019] [Accepted: 06/03/2019] [Indexed: 12/30/2022]
Abstract
Community genetics aims to understand the effects of intraspecific genetic variation on community composition and diversity, thereby connecting community ecology with evolutionary biology. Thus far, research has shown that plant genetics can underlie variation in the composition of associated communities (e.g., insects, lichen and endophytes), and those communities can therefore be considered as extended phenotypes. This work, however, has been conducted primarily at the plant genotype level and has not identified the key underlying genes. To address this gap, we used genome-wide association mapping with a population of 445 aspen (Populus tremuloides) genets to identify the genes governing variation in plant traits (defence chemistry, bud phenology, leaf morphology, growth) and insect community composition. We found 49 significant SNP associations in 13 Populus genes that are correlated with chemical defence compounds and insect community traits. Most notably, we identified an early nodulin-like protein that was associated with insect community diversity and the abundance of interacting foundation species (ants and aphids). These findings support the concept that particular plant traits are the mechanistic link between plant genes and the composition of associated insect communities. In putting the "genes" into "genes to ecosystems ecology", this work enhances understanding of the molecular genetic mechanisms that underlie plant-insect associations and the consequences thereof for the structure of ecological communities.
Collapse
Affiliation(s)
- Hilary L Barker
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, WI, USA
| | - Jennifer F Riehl
- Department of Entomology, University of Wisconsin-Madison, Madison, WI, USA
| | | | | | - Liza M Holeski
- Department of Biological Sciences, Northern Arizona University, Flagstaff, AZ, USA
| | - Pär K Ingvarsson
- Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Richard L Lindroth
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, WI, USA.,Department of Entomology, University of Wisconsin-Madison, Madison, WI, USA
| |
Collapse
|
36
|
Zhang B, Zhu W, Diao S, Wu X, Lu J, Ding C, Su X. The poplar pangenome provides insights into the evolutionary history of the genus. Commun Biol 2019; 2:215. [PMID: 31240253 PMCID: PMC6581948 DOI: 10.1038/s42003-019-0474-7] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2018] [Accepted: 05/23/2019] [Indexed: 12/20/2022] Open
Abstract
The genus Populus comprises a complex amalgam of ancient and modern species that has become a prime model for evolutionary and taxonomic studies. Here we sequenced the genomes of 10 species from five sections of the genus Populus, identified 71 million genomic variations, and observed new correlations between the single-nucleotide polymorphism-structural variation (SNP-SV) density and indel-SV density to complement the SNP-indel density correlation reported in mammals. Disease resistance genes (R genes) with heterozygous loss-of-function (LOF) were significantly enriched in the 10 species, which increased the diversity of poplar R genes during evolution. Heterozygous LOF mutations in the self-incompatibility genes were closely related to the self-fertilization of poplar, suggestive of genomic control of self-fertilization in dioecious plants. The phylogenetic genome-wide SNPs tree also showed possible ancient hybridization among species in sections Tacamahaca, Aigeiros, and Leucoides. The pangenome resource also provided information for poplar genetics and breeding.
Collapse
Affiliation(s)
- Bingyu Zhang
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
- Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
| | - Wenxu Zhu
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
- Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
| | - Shu Diao
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
- Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
| | - Xiaojuan Wu
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
- Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
| | - Junqian Lu
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
- Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
| | - ChangJun Ding
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
- Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
| | - Xiaohua Su
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
- Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, 100091 Beijing, China
| |
Collapse
|
37
|
Acosta JJ, Fahrenkrog AM, Neves LG, Resende MFR, Dervinis C, Davis JM, Holliday JA, Kirst M. Exome Resequencing Reveals Evolutionary History, Genomic Diversity, and Targets of Selection in the Conifers Pinus taeda and Pinus elliottii. Genome Biol Evol 2019; 11:508-520. [PMID: 30689841 PMCID: PMC6385631 DOI: 10.1093/gbe/evz016] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/22/2019] [Indexed: 12/22/2022] Open
Abstract
Loblolly pine (Pinus taeda) and slash pine (Pinus elliottii) are ecologically and economically important pine species that dominate many forest ecosystems in the southern United States, but like all conifers, the study of their genetic diversity and demographic history has been hampered by their large genome size. A small number of studies mainly based on candidate-gene sequencing have been reported for P. taeda to date, whereas none are available for P. elliottii. Targeted exome resequencing has recently enabled population genomics studies for conifers, approach used here to assess genomic diversity, signatures of selection, population structure, and demographic history of P. elliottii and P. taeda. Extensive similarities were revealed between these species: both species feature rapid linkage disequilibrium decay and high levels of genetic diversity. Moreover, genome-wide positive correlations for measures of genetic diversity between the species were also observed, likely due to shared structural genomic constraints. Also, positive selection appears to be targeting a common set of genes in both pines. Demographic history differs between both species, with only P. taeda being affected by a dramatic bottleneck during the last glacial period. The ability of P. taeda to recover from a dramatic reduction in population size while still retaining high levels of genetic diversity shows promise for other pines facing environmental stressors associated with climate change, indicating that these too may be able to adapt successfully to new future conditions even after a drastic population size contraction.
Collapse
Affiliation(s)
- Juan J Acosta
- School of Forest Resources and Conservation, University of Florida.,University of Florida Genetics Institute, University of Florida.,Camcore, Department of Forestry and Environmental Resources, North Carolina State University, Raleigh, NC
| | - Annette M Fahrenkrog
- School of Forest Resources and Conservation, University of Florida.,Plant Molecular and Cellular Biology Graduate Program, University of Florida
| | - Leandro G Neves
- School of Forest Resources and Conservation, University of Florida.,Plant Molecular and Cellular Biology Graduate Program, University of Florida.,RAPiD Genomics, Gainesville, FL
| | | | | | - John M Davis
- School of Forest Resources and Conservation, University of Florida
| | - Jason A Holliday
- Department of Forest Resources and Environmental Conservation, Virginia Polytechnic Institute and State University
| | - Matias Kirst
- School of Forest Resources and Conservation, University of Florida.,Plant Molecular and Cellular Biology Graduate Program, University of Florida.,University of Florida Genetics Institute, University of Florida
| |
Collapse
|
38
|
Silva-Junior OB, Grattapaglia D, Novaes E, Collevatti RG. Design and evaluation of a sequence capture system for genome-wide SNP genotyping in highly heterozygous plant genomes: a case study with a keystone Neotropical hardwood tree genome. DNA Res 2019; 25:535-545. [PMID: 30020434 PMCID: PMC6191306 DOI: 10.1093/dnares/dsy023] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2017] [Accepted: 06/22/2018] [Indexed: 12/12/2022] Open
Abstract
Targeted sequence capture coupled to high-throughput sequencing has become a powerful method for the study of genome-wide sequence variation. Following our recent development of a genome assembly for the Pink Ipê tree (Handroanthus impetiginosus), a widely distributed Neotropical timber species, we now report the development of a set of 24,751 capture probes for single-nucleotide polymorphisms (SNPs) characterization and genotyping across 18,216 distinct loci, sampling more than 10 Mbp of the species genome. This system identifies nearly 200,000 SNPs located inside or in close proximity to almost 14,000 annotated protein-coding genes, generating quality genotypic data in populations spanning wide geographic distances across the species native range. To provide recommendations for future developments of similar systems for highly heterozygous plant genomes we investigated issues such as probe design, sequencing coverage and bioinformatics, including the evaluation of the capture efficiency and a reassessment of the technical reproducibility of the assay for SNPs recall and genotyping precision. Our results highlight the value of a detailed probe screening on a preliminary genome assembly to produce reliable data for downstream genetic studies. This work should inspire and assist the development of similar genomic resources for other orphan crops and forest trees with highly heterozygous genomes.
Collapse
Affiliation(s)
- Orzenil Bonfim Silva-Junior
- EMBRAPA Recursos Genéticos e Biotecnologia, EPqB, Brasília, DF, Brazil.,Programa de Ciências Genômicas e Biotecnologia, Universidade Católica de Brasília, SGAN 916 Modulo B, Brasilia, DF, Brazil
| | - Dario Grattapaglia
- EMBRAPA Recursos Genéticos e Biotecnologia, EPqB, Brasília, DF, Brazil.,Programa de Ciências Genômicas e Biotecnologia, Universidade Católica de Brasília, SGAN 916 Modulo B, Brasilia, DF, Brazil
| | - Evandro Novaes
- Departamento de Biologia, Universidade Federal de Lavras, Lavras, MG, Brazil
| | - Rosane G Collevatti
- Laboratório de Genética & Biodiversidade, Instituto de Ciências Biológicas, Universidade Federal de Goiás, Goiânia, GO, Brazil
| |
Collapse
|
39
|
Du Q, Lu W, Quan M, Xiao L, Song F, Li P, Zhou D, Xie J, Wang L, Zhang D. Genome-Wide Association Studies to Improve Wood Properties: Challenges and Prospects. FRONTIERS IN PLANT SCIENCE 2018; 9:1912. [PMID: 30622554 PMCID: PMC6309013 DOI: 10.3389/fpls.2018.01912] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Accepted: 12/10/2018] [Indexed: 05/02/2023]
Abstract
Wood formation is an excellent model system for quantitative trait analysis due to the strong associations between the transcriptional and metabolic traits that contribute to this complex process. Investigating the genetic architecture and regulatory mechanisms underlying wood formation will enhance our understanding of the quantitative genetics and genomics of complex phenotypic variation. Genome-wide association studies (GWASs) represent an ideal statistical strategy for dissecting the genetic basis of complex quantitative traits. However, elucidating the molecular mechanisms underlying many favorable loci that contribute to wood formation and optimizing GWAS design remain challenging in this omics era. In this review, we summarize the recent progress in GWAS-based functional genomics of wood property traits in major timber species such as Eucalyptus, Populus, and various coniferous species. We discuss several appropriate experimental designs for extensive GWAS in a given undomesticated tree population, such as omics-wide association studies and high-throughput phenotyping technologies. We also explain why more attention should be paid to rare allelic and major structural variation. Finally, we explore the potential use of GWAS for the molecular breeding of trees. Such studies will help provide an integrated understanding of complex quantitative traits and should enable the molecular design of new cultivars.
Collapse
Affiliation(s)
- Qingzhang Du
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Wenjie Lu
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Mingyang Quan
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Liang Xiao
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Fangyuan Song
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Peng Li
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Daling Zhou
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Jianbo Xie
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Longxin Wang
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
| | - Deqiang Zhang
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 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
|
40
|
Corcoran P, Gossmann TI, Barton HJ, Slate J, Zeng K. Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species. Genome Biol Evol 2018; 9:2987-3007. [PMID: 29045655 PMCID: PMC5714183 DOI: 10.1093/gbe/evx213] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/16/2017] [Indexed: 02/06/2023] Open
Abstract
Population genetic theory predicts that selection should be more effective when the effective population size (Ne) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger Ne. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.
Collapse
Affiliation(s)
- Pádraic Corcoran
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | - Toni I Gossmann
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | - Henry J Barton
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | | | - Jon Slate
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| | - Kai Zeng
- Department of Animal and Plant Sciences, University of Sheffield, South Yorkshire, United Kingdom
| |
Collapse
|
41
|
Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen. Proc Natl Acad Sci U S A 2018; 115:E10970-E10978. [PMID: 30373829 PMCID: PMC6243237 DOI: 10.1073/pnas.1801437115] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We performed de novo, full-genome sequence analysis of two Populus species, North American quaking and Eurasian trembling aspen, that contain striking levels of genetic variation. Our results showed that positive and negative selection broadly affects patterns of genomic variation, but to varying degrees across coding and noncoding regions. The strength of selection and rates of sequence divergence were strongly related to differences in gene expression and coexpression network connectivity. These results highlight the importance of both positive and negative selection in shaping genome-wide levels of genetic variation in an obligately outcrossing, perennial plant. The resources we present establish aspens as a powerful study system enabling future studies for understanding the genomic determinants of adaptive evolution. The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).
Collapse
|
42
|
Chhatre VE, Evans LM, DiFazio SP, Keller SR. Adaptive introgression and maintenance of a trispecies hybrid complex in range‐edge populations of
Populus. Mol Ecol 2018; 27:4820-4838. [DOI: 10.1111/mec.14820] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2017] [Revised: 07/19/2018] [Accepted: 07/20/2018] [Indexed: 12/16/2022]
Affiliation(s)
- Vikram E. Chhatre
- Department of Plant Biology University of Vermont Burlington Vermont
| | - Luke M. Evans
- Department of Ecology and Evolutionary Biology Institute of Behavioral Genetics University of Colorado Boulder Colorado
| | | | - Stephen R. Keller
- Department of Plant Biology University of Vermont Burlington Vermont
| |
Collapse
|
43
|
Jia Y, Zhu J, Wu Y, Fan WB, Zhao GF, Li ZH. Effects of Geological and Environmental Events on the Diversity and Genetic Divergence of Four Closely Related Pines: Pinus koraiensis, P. armandii, P. griffithii, and P. pumila. FRONTIERS IN PLANT SCIENCE 2018; 9:1264. [PMID: 30210523 PMCID: PMC6121107 DOI: 10.3389/fpls.2018.01264] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Accepted: 08/10/2018] [Indexed: 06/08/2023]
Abstract
The effects of mountain uplift and environmental oscillations on nucleotide variability and species divergence remain largely unknown in East Asia. In this study, based on multiple nuclear DNA markers, we investigated the levels and patterns of nucleotide diversity and interspecific divergence in four closely related pines in China, i.e., Pinus koraiensis, P. armandii, P. griffithii, and P. pumila. The four pine taxa shared low levels of nucleotide polymorphisms at the species level. P. pumila had the highest silent nucleotide diversity (πsil = 0.00661) whereas P. griffithii had the lowest (πsil = 0.00175), while the levels of genetic polymorphism in P. armandii (πsil = 0.00508) and P. koraiensis (πsil = 0.00652) were intermediate between the other two species. Population genetic structure analysis showed that variations primarily existed within populations of the four pine species, presumably due to habitat fragmentation or the island-like distributions of Pinus species. Population divergence (FST) analysis showed that the genetic divergence between P. griffithii and P. koraiensis was much greater than that between P. koraiensis and the other two pines species. Isolation-with-migration analysis suggested that asymmetric gene flow had occurred between any two pairs of pine species. Phylogenetic analyses indicated that the four allied species split into two groups about 1.37 million years ago, where P. armandii and P. pumila were closer and clustered as sister species, whereas P. koraiensis and P. griffithii were clustered on another branch. Our results and those obtained in previous studies suggest that mountain uplift and geological climate oscillations may have led to the patterns of genetic divergence and nucleotide variations in these four pine species.
Collapse
Affiliation(s)
| | | | | | | | | | - Zhong-Hu Li
- Key Laboratory of Resource Biology and Biotechnology in Western China, Ministry of Education, College of Life Sciences, Northwest University, Xi’an, China
| |
Collapse
|
44
|
Cole CT, Ingvarsson PK. Pathway position constrains the evolution of an ecologically important pathway in aspens (Populus tremula L.). Mol Ecol 2018; 27:3317-3330. [PMID: 29972878 DOI: 10.1111/mec.14785] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Revised: 01/30/2018] [Accepted: 02/20/2018] [Indexed: 12/22/2022]
Abstract
Many ecological interactions of aspens and their relatives (Populus spp.) are affected by products of the phenylpropanoid pathway synthesizing condensed tannins (CTs), whose production involves trade-offs with other ecologically important compounds and with growth. Genes of this pathway are candidates for investigating the role of selection on ecologically important, polygenic traits. We analysed sequences from 25 genes representing 10 steps of the CT synthesis pathway, which produces CTs used in defence and lignins used for growth, in 12 individuals of European aspen (Populus tremula). We compared these to homologs from P. trichocarpa, to a control set of 77 P. tremula genes, to genome-wide resequencing data and to RNA-seq expression levels, in order to identify signatures of selection distinct from those of demography. In Populus, pathway position exerts a strong influence on the evolution of these genes. Nonsynonymous diversity, divergence and allele frequency shifts (Tajima's D) were much lower than for synonymous measures. Expression levels were higher, and the direction of selection more negative, for upstream genes than for those downstream. Selective constraints act with increasing intensity on upstream genes, despite the presence of multiple paralogs in most gene families. Pleiotropy, expression level, flux control and codon bias appear to interact in determining levels and patterns of variation in genes of this pathway, whose products mediate a wide array of ecological interactions for this widely distributed species.
Collapse
Affiliation(s)
- Christopher T Cole
- Division of Science and Mathematics, University of Minnesota, Morris, Morris, Minnesota
| | - Pär K Ingvarsson
- Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden
| |
Collapse
|
45
|
Wang J, Ding J, Tan B, Robinson KM, Michelson IH, Johansson A, Nystedt B, Scofield DG, Nilsson O, Jansson S, Street NR, Ingvarsson PK. A major locus controls local adaptation and adaptive life history variation in a perennial plant. Genome Biol 2018; 19:72. [PMID: 29866176 PMCID: PMC5985590 DOI: 10.1186/s13059-018-1444-y] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2017] [Accepted: 05/03/2018] [Indexed: 02/17/2023] Open
Abstract
BACKGROUND The initiation of growth cessation and dormancy represent critical life-history trade-offs between survival and growth and have important fitness effects in perennial plants. Such adaptive life-history traits often show strong local adaptation along environmental gradients but, despite their importance, the genetic architecture of these traits remains poorly understood. RESULTS We integrate whole genome re-sequencing with environmental and phenotypic data from common garden experiments to investigate the genomic basis of local adaptation across a latitudinal gradient in European aspen (Populus tremula). A single genomic region containing the PtFT2 gene mediates local adaptation in the timing of bud set and explains 65% of the observed genetic variation in bud set. This locus is the likely target of a recent selective sweep that originated right before or during colonization of northern Scandinavia following the last glaciation. Field and greenhouse experiments confirm that variation in PtFT2 gene expression affects the phenotypic variation in bud set that we observe in wild natural populations. CONCLUSIONS Our results reveal a major effect locus that determines the timing of bud set and that has facilitated rapid adaptation to shorter growing seasons and colder climates in European aspen. The discovery of a single locus explaining a substantial fraction of the variation in a key life-history trait is remarkable, given that such traits are generally considered to be highly polygenic. These findings provide a dramatic illustration of how loci of large-effect for adaptive traits can arise and be maintained over large geographical scales in natural populations.
Collapse
Affiliation(s)
- Jing Wang
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, 90187, Umeå, Sweden.
- Centre for Integrative Genetics, Department of Animal and Aquacultural Sciences, Faculty of Life Sciences, Norwegian University of Life Sciences, PO Box 5003, Ås, Norway.
| | - Jihua Ding
- Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, 901 83, Umeå, Sweden
| | - Biyue Tan
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, 90187, Umeå, Sweden
- Stora Enso Biomaterials, 13104, Nacka, Sweden
| | - Kathryn M Robinson
- Umeå Plant Science Centre, Department of Plant Physiology, Umeå University, 90187, Umeå, Sweden
| | - Ingrid H Michelson
- Umeå Plant Science Centre, Department of Plant Physiology, Umeå University, 90187, Umeå, Sweden
| | - Anna Johansson
- Wallenberg Advanced Bioinformatics Infrastructure, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Björn Nystedt
- Wallenberg Advanced Bioinformatics Infrastructure, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Douglas G Scofield
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, 90187, Umeå, Sweden
- Department of Ecology and Genetics, Evolutionary Biology, Uppsala University, Uppsala, Sweden
- Uppsala Multidisciplinary Center for Advanced Computational Science, Uppsala University, Uppsala, Sweden
| | - Ove Nilsson
- Umeå Plant Science Centre, Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, 901 83, Umeå, Sweden
| | - Stefan Jansson
- Umeå Plant Science Centre, Department of Plant Physiology, Umeå University, 90187, Umeå, Sweden
| | - Nathaniel R Street
- Umeå Plant Science Centre, Department of Plant Physiology, Umeå University, 90187, Umeå, Sweden
| | - Pär K Ingvarsson
- Umeå Plant Science Centre, Department of Ecology and Environmental Science, Umeå University, 90187, Umeå, Sweden.
- Present address: Department of Plant Biology, Uppsala BioCenter, Swedish University of Agricultural Sciences, PO Box 7080, 750 07, Uppsala, Sweden.
| |
Collapse
|
46
|
Morata J, Tormo M, Alexiou KG, Vives C, Ramos-Onsins SE, Garcia-Mas J, Casacuberta JM. The Evolutionary Consequences of Transposon-Related Pericentromer Expansion in Melon. Genome Biol Evol 2018; 10:1584-1595. [PMID: 29901717 PMCID: PMC6009578 DOI: 10.1093/gbe/evy115] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/04/2018] [Indexed: 12/16/2022] Open
Abstract
Transposable elements (TEs) are a major driver of plant genome evolution. A part from being a rich source of new genes and regulatory sequences, TEs can also affect plant genome evolution by modifying genome size and shaping chromosome structure. TEs tend to concentrate in heterochromatic pericentromeric regions and their proliferation may expand these regions. Here, we show that after the split of melon and cucumber, TEs have expanded the pericentromeric regions of melon chromosomes that, probably as a consequence, show a very low recombination frequency. In contrast, TEs have not proliferated to a high extent in cucumber, which has small TE-dense pericentromeric regions and shows a relatively constant recombination rate along chromosomes. These differences in chromosome structure also translate in differences in gene nucleotide diversity. Although gene nucleotide diversity is essentially constant along cucumber chromosomes, melon chromosomes show a bimodal pattern of genetic variability, with a gene-poor region where variability is negatively correlated with gene density. Interestingly, genes are not homogeneously distributed in melon, and the high variable low-recombining pericentromeric regions show a higher concentration of melon-specific genes whereas genes shared with cucumber and other plants are essentially found in gene-rich chromosomal arms. The results presented here suggest that melon pericentromeric regions may allow gene sequences to evolve more freely than in other chromosomal compartments which may allow new ORFs to arise and eventually be selected. These results show that TEs can drastically change the structure of chromosomes creating different chromosomal compartments imposing different constraints for gene evolution.
Collapse
Affiliation(s)
- Jordi Morata
- Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Marc Tormo
- Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Konstantinos G Alexiou
- Institut de Recerca i Tecnologia Agroalimentàries, Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Cristina Vives
- Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Sebastián E Ramos-Onsins
- Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Jordi Garcia-Mas
- Institut de Recerca i Tecnologia Agroalimentàries, Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Josep M Casacuberta
- Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| |
Collapse
|
47
|
Marand AP, Jansky SH, Zhao H, Leisner CP, Zhu X, Zeng Z, Crisovan E, Newton L, Hamernik AJ, Veilleux RE, Buell CR, Jiang J. Meiotic crossovers are associated with open chromatin and enriched with Stowaway transposons in potato. Genome Biol 2017; 18:203. [PMID: 29084572 PMCID: PMC5663088 DOI: 10.1186/s13059-017-1326-8] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2017] [Accepted: 09/25/2017] [Indexed: 12/25/2022] Open
Abstract
Background Meiotic recombination is the foundation for genetic variation in natural and artificial populations of eukaryotes. Although genetic maps have been developed for numerous plant species since the late 1980s, few of these maps have provided the necessary resolution needed to investigate the genomic and epigenomic features underlying meiotic crossovers. Results Using a whole genome sequencing-based approach, we developed two high-density reference-based haplotype maps using diploid potato clones as parents. The vast majority (81%) of meiotic crossovers were mapped to less than 5 kb. The fine-scale accuracy of crossover detection was validated by Sanger sequencing for a subset of ten crossover events. We demonstrate that crossovers reside in genomic regions of “open chromatin”, which were identified based on hypersensitivity to DNase I digestion and association with H3K4me3-modified nucleosomes. The genomic regions spanning crossovers were significantly enriched with the Stowaway family of miniature inverted-repeat transposable elements (MITEs). The occupancy of Stowaway elements in gene promoters is concomitant with an increase in recombination rate. A generalized linear model identified the presence of Stowaway elements as the third most important genomic or chromatin feature behind genes and open chromatin for predicting crossover formation over 10-kb windows. Conclusions Collectively, our results suggest that meiotic crossovers in potato are largely determined by the local chromatin status, marked by accessible chromatin, H3K4me3-modified nucleosomes, and the presence of Stowaway transposons. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1326-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Alexandre P Marand
- Department of Horticulture, University of Wisconsin-Madison, Madison, Wisconsin, 53706, USA
| | - Shelley H Jansky
- Department of Horticulture, University of Wisconsin-Madison, Madison, Wisconsin, 53706, USA. .,USDA-ARS, Vegetable Crops Research Unit, Madison, Wisconsin, 53706, USA.
| | - Hainan Zhao
- Department of Horticulture, University of Wisconsin-Madison, Madison, Wisconsin, 53706, USA
| | - Courtney P Leisner
- Department of Plant Biology, Michigan State University, East Lansing, Michigan, 48824, USA
| | - Xiaobiao Zhu
- Department of Horticulture, University of Wisconsin-Madison, Madison, Wisconsin, 53706, USA
| | - Zixian Zeng
- Department of Horticulture, University of Wisconsin-Madison, Madison, Wisconsin, 53706, USA
| | - Emily Crisovan
- Department of Plant Biology, Michigan State University, East Lansing, Michigan, 48824, USA
| | - Linsey Newton
- Department of Plant Biology, Michigan State University, East Lansing, Michigan, 48824, USA
| | - Andy J Hamernik
- Department of Horticulture, University of Wisconsin-Madison, Madison, Wisconsin, 53706, USA.,USDA-ARS, Vegetable Crops Research Unit, Madison, Wisconsin, 53706, USA
| | | | - C Robin Buell
- Department of Plant Biology, Michigan State University, East Lansing, Michigan, 48824, USA
| | - Jiming Jiang
- Department of Horticulture, University of Wisconsin-Madison, Madison, Wisconsin, 53706, USA. .,Current address: Departments of Plant Biology and Horticulture, Michigan State University, East Lansing, Michigan, 48824, USA.
| |
Collapse
|
48
|
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: 111] [Impact Index Per Article: 15.9] [Reference Citation Analysis] [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
|
49
|
Fahrenkrog AM, Neves LG, Resende MFR, Dervinis C, Davenport R, Barbazuk WB, Kirst M. Population genomics of the eastern cottonwood ( Populus deltoides). Ecol Evol 2017; 7:9426-9440. [PMID: 29187979 PMCID: PMC5696417 DOI: 10.1002/ece3.3466] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2017] [Revised: 08/14/2017] [Accepted: 08/16/2017] [Indexed: 12/30/2022] Open
Abstract
Despite its economic importance as a bioenergy crop and key role in riparian ecosystems, little is known about genetic diversity and adaptation of the eastern cottonwood (Populus deltoides). Here, we report the first population genomics study for this species, conducted on a sample of 425 unrelated individuals collected in 13 states of the southeastern United States. The trees were genotyped by targeted resequencing of 18,153 genes and 23,835 intergenic regions, followed by the identification of single nucleotide polymorphisms (SNPs). This natural P. deltoides population showed low levels of subpopulation differentiation (FST = 0.022–0.106), high genetic diversity (θW = 0.00100, π = 0.00170), a large effective population size (Ne ≈ 32,900), and low to moderate levels of linkage disequilibrium. Additionally, genomewide scans for selection (Tajima's D), subpopulation differentiation (XTX), and environmental association analyses with eleven climate variables carried out with two different methods (LFMM and BAYENV2) identified genes putatively involved in local adaptation. Interestingly, many of these genes were also identified as adaptation candidates in another poplar species, Populus trichocarpa, indicating possible convergent evolution. This study constitutes the first assessment of genetic diversity and local adaptation in P. deltoides throughout the southern part of its range, information we expect to be of use to guide management and breeding strategies for this species in future, especially in the face of climate change.
Collapse
Affiliation(s)
- Annette M Fahrenkrog
- School of Forest Resources and Conservation University of Florida Gainesville FL USA.,Plant Molecular and Cellular Biology Graduate Program University of Florida Gainesville FL USA
| | - Leandro G Neves
- School of Forest Resources and Conservation University of Florida Gainesville FL USA.,Plant Molecular and Cellular Biology Graduate Program University of Florida Gainesville FL USA.,Present address: RAPiD Genomics LLC756 2nd Avenue Gainesville FL 32601 USA
| | - Márcio F R Resende
- Horticultural Sciences Department University of Florida Gainesville FL USA
| | - Christopher Dervinis
- School of Forest Resources and Conservation University of Florida Gainesville FL USA
| | - Ruth Davenport
- Biology Department University of Florida Gainesville FL USA
| | - W Brad Barbazuk
- Biology Department University of Florida Gainesville FL USA.,University of Florida Genetics Institute University of Florida Gainesville FL USA
| | - Matias Kirst
- School of Forest Resources and Conservation University of Florida Gainesville FL USA.,University of Florida Genetics Institute University of Florida Gainesville FL USA
| |
Collapse
|
50
|
Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch. Nat Genet 2017; 49:904-912. [PMID: 28481341 DOI: 10.1038/ng.3862] [Citation(s) in RCA: 125] [Impact Index Per Article: 17.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2017] [Accepted: 04/12/2017] [Indexed: 12/26/2022]
Abstract
Silver birch (Betula pendula) is a pioneer boreal tree that can be induced to flower within 1 year. Its rapid life cycle, small (440-Mb) genome, and advanced germplasm resources make birch an attractive model for forest biotechnology. We assembled and chromosomally anchored the nuclear genome of an inbred B. pendula individual. Gene duplicates from the paleohexaploid event were enriched for transcriptional regulation, whereas tandem duplicates were overrepresented by environmental responses. Population resequencing of 80 individuals showed effective population size crashes at major points of climatic upheaval. Selective sweeps were enriched among polyploid duplicates encoding key developmental and physiological triggering functions, suggesting that local adaptation has tuned the timing of and cross-talk between fundamental plant processes. Variation around the tightly-linked light response genes PHYC and FRS10 correlated with latitude and longitude and temperature, and with precipitation for PHYC. Similar associations characterized the growth-promoting cytokinin response regulator ARR1, and the wood development genes KAK and MED5A.
Collapse
|