1
|
Ma XG, Ren YB, Sun H. Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China. Cladistics 2024; 40:357-373. [PMID: 38197450 DOI: 10.1111/cla.12570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 11/27/2023] [Accepted: 12/09/2023] [Indexed: 01/11/2024] Open
Abstract
Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.
Collapse
Affiliation(s)
- Xiang-Guang Ma
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, 650201, China
| | - Yue-Bo Ren
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, 650201, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Hang Sun
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, 650201, China
| |
Collapse
|
2
|
Wagner S, Seguin-Orlando A, Leplé JC, Leroy T, Lalanne C, Labadie K, Aury JM, Poirier S, Wincker P, Plomion C, Kremer A, Orlando L. Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood. Mol Ecol 2024; 33:e16859. [PMID: 36748324 PMCID: PMC7615563 DOI: 10.1111/mec.16859] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2022] [Revised: 12/08/2022] [Accepted: 01/16/2023] [Indexed: 02/08/2023]
Abstract
Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.
Collapse
Affiliation(s)
- Stefanie Wagner
- Plant Genomic Resources Center (CNRGV), INRAE, Castanet-Tolosan, France
- Centre for Anthropobiology and Genomics of Toulouse (CAGT), CNRS UMR 5288, Université Paul Sabatier, Toulouse, France
| | - Andaine Seguin-Orlando
- Centre for Anthropobiology and Genomics of Toulouse (CAGT), CNRS UMR 5288, Université Paul Sabatier, Toulouse, France
| | | | - Thibault Leroy
- IRHS UMR1345, Université d’Angers, INRAE, Institut Agro, SFR 4207 QuaSaV, Beaucouzé, France
| | | | - Karine Labadie
- Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, Evry, France
| | - Jean-Marc Aury
- Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, Evry, France
| | | | - Patrick Wincker
- Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, Evry, France
| | | | | | - Ludovic Orlando
- Centre for Anthropobiology and Genomics of Toulouse (CAGT), CNRS UMR 5288, Université Paul Sabatier, Toulouse, France
| |
Collapse
|
3
|
Schroeder H, Kersten B. A Small Set of Nuclear Markers for Reliable Differentiation of the Two Closely Related Oak Species Quercus Robur and Q. Petraea. PLANTS (BASEL, SWITZERLAND) 2023; 12:566. [PMID: 36771649 PMCID: PMC9919172 DOI: 10.3390/plants12030566] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Revised: 01/19/2023] [Accepted: 01/24/2023] [Indexed: 06/18/2023]
Abstract
Quercus robur and Q. petraea are, in addition to Fagus sylvatica, the main economically used deciduous tree species in Europe. Identification of these two species is crucial because they differ in their ecological demands. Because of a changing climate, foresters must know more than ever which species will perform better under given environmental conditions. The search for differentiating molecular markers between these two species has already lasted for decades. Until now, differentiation has only been possible in approaches with a combination of several molecular markers and a subsequent statistical analysis to calculate the probability of being one or the other species. Here, we used MiSeq Illumina data from pools of Q. robur and Q. petraea specimens and identified nuclear SNPs and small InDels versus the Q. robur reference genome. Selected sequence variants with 100% allele frequency difference between the two pools were further validated in an extended set of Q. robur and Q. petraea specimens, and then the number of markers was deliberately reduced to the smallest possible set for species differentiation. A combination of six markers from four nuclear regions is enough to identify Q. robur, Q. petraea or hybrids between these two species quite well and represents a marker set that is cost-efficient and useable in every laboratory.
Collapse
|
4
|
Genetic Variability and Admixture Zones in the Italian Populations of Turkey Oak ( Quercus cerris L.). LIFE (BASEL, SWITZERLAND) 2022; 13:life13010018. [PMID: 36675965 PMCID: PMC9863001 DOI: 10.3390/life13010018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/06/2022] [Revised: 12/02/2022] [Accepted: 12/16/2022] [Indexed: 12/24/2022]
Abstract
The Turkey oak (Quercus cerris L.) is widely distributed in Italy, where it is the ecologically dominant oak on sandy and acidic soil. In this work, we analysed 23 natural populations by means of eight SSR (microsatellite) markers, to obtain the first synthetic map of genetic variability for this species and to study its dispersion during the Holocene, due to the possibility that at least one refugium during the Last Glacial Maximum was in Italy. The analyses showed a good amount of genetic variability together with fair differentiation between populations, as indicated by FST = 0.059. A Bayesian analysis of the amount of admixture among populations revealed the presence of four putative gene pools of origin and a rough subdivision of the populations according to their geographic location, as confirmed by the spatial analysis. No evidence for the existence of putative refugial populations was found; however, this study paves the way for the planning of conservation strategies also with regard to the relationship between Turkey oak and other oak species in Italy.
Collapse
|
5
|
Mader M, Schroeder H, Schott T, Schöning-Stierand K, Leite Montalvão AP, Liesebach H, Liesebach M, Fussi B, Kersten B. Mitochondrial Genome of Fagus sylvatica L. as a Source for Taxonomic Marker Development in the Fagales. PLANTS (BASEL, SWITZERLAND) 2020; 9:E1274. [PMID: 32992588 PMCID: PMC7650814 DOI: 10.3390/plants9101274] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 09/23/2020] [Accepted: 09/24/2020] [Indexed: 11/16/2022]
Abstract
European beech, Fagus sylvatica L., is one of the most important and widespread deciduous tree species in Central Europe and is widely managed for its hard wood. The complete DNA sequence of the mitochondrial genome of Fagus sylvatica L. was assembled and annotated based on Illumina MiSeq reads and validated using long reads from nanopore MinION sequencing. The genome assembled into a single DNA sequence of 504,715 bp in length containing 58 genes with predicted function, including 35 protein-coding, 20 tRNA and three rRNA genes. Additionally, 23 putative protein-coding genes were predicted supported by RNA-Seq data. Aiming at the development of taxon-specific mitochondrial genetic markers, the tool SNPtax was developed and applied to select genic SNPs potentially specific for different taxa within the Fagales. Further validation of a small SNP set resulted in the development of four CAPS markers specific for Fagus, Fagaceae, or Fagales, respectively, when considering over 100 individuals from a total of 69 species of deciduous trees and conifers from up to 15 families included in the marker validation. The CAPS marker set is suitable to identify the genus Fagus in DNA samples from tree tissues or wood products, including wood composite products.
Collapse
Affiliation(s)
- Malte Mader
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
| | - Hilke Schroeder
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
| | - Thomas Schott
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
| | - Katrin Schöning-Stierand
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
- Center for Bioinformatics, Universität Hamburg, 20146 Hamburg, Germany
| | - Ana Paula Leite Montalvão
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
| | - Heike Liesebach
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
| | - Mirko Liesebach
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
| | - Barbara Fussi
- Bavarian Office for Forest Genetics, 83317 Teisendorf, Germany;
| | - Birgit Kersten
- Thünen Institute of Forest Genetics, D-22927 Grosshansdorf, Germany; (M.M.); (H.S.); (T.S.); (K.S.-S.); (A.P.L.M.); (H.L.); (M.L.)
| |
Collapse
|
6
|
Brousseau L, Fine PVA, Dreyer E, Vendramin GG, Scotti I. Genomic and phenotypic divergence unveil microgeographic adaptation in the Amazonian hyperdominant tree Eperua falcata Aubl. (Fabaceae). Mol Ecol 2020; 30:1136-1154. [PMID: 32786115 DOI: 10.1111/mec.15595] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2018] [Revised: 06/19/2020] [Accepted: 07/31/2020] [Indexed: 01/04/2023]
Abstract
Plant populations can undergo very localized adaptation, allowing widely distributed populations to adapt to divergent habitats in spite of recurrent gene flow. Neotropical trees-whose large and undisturbed populations often span a variety of environmental conditions and local habitats-are particularly good models to study this process. Here, we explore patterns of adaptive divergence from large (i.e., regional) to small (i.e., microgeographic) spatial scales in the hyperdominant Amazonian tree Eperua falcata Aubl. (Fabaceae) under a replicated design involving two microhabitats (~300 m apart) in two study sites (~300 km apart). A three-year reciprocal transplant illustrates that, beyond strong maternal effects and phenotypic plasticity, genetically driven divergence in seedling growth and leaf traits was detected both between seedlings originating from different regions, and between seedlings from different microhabitats. In parallel, a complementary genome scan for selection was carried out through whole-genome sequencing of tree population pools. A set of 290 divergence outlier SNPs was detected at the regional scale (between study sites), while 185 SNPs located in the vicinity of 106 protein-coding genes were detected as replicated outliers between microhabitats within regions. Outlier-surrounding genomic regions are involved in a variety of physiological processes, including plant responses to stress (e.g., oxidative stress, hypoxia and metal toxicity) and biotic interactions. Together with evidence of microgeographic divergence in functional traits, the discovery of genomic candidates for microgeographic adaptive divergence represents a promising advance in our understanding of local adaptation, which probably operates across multiple spatial scales and underpins divergence and diversification in Neotropical trees.
Collapse
Affiliation(s)
- Louise Brousseau
- UMR EcoFoG, AgroParisTech, CIRAD, CNRS, INRAE, Université de Guyane, Université des Antilles, Kourou Cedex, France.,AMAP, Univ. Montpellier, CIRAD, CNRS, INRAE, IRD, Montpellier, France
| | - Paul V A Fine
- Department of Integrative Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Erwin Dreyer
- Université de Lorraine, AgroParisTech, INRAE, Silva, Nancy, France
| | - Giovanni G Vendramin
- Institute of Biosciences and BioResources (IBBR-CNR), National Research Council, Division of Florence, Sesto Fiorentino, Italy
| | - Ivan Scotti
- UR629 Ecologie des Forêts Méditerranéennes (URFM), INRAE, Avignon, France
| |
Collapse
|