Niche evolution in a northern temperate tree lineage: biogeographical legacies in cork oaks (Quercus section Cerris).

BACKGROUND AND AIMS: Cork oaks (Quercus section Cerris) comprise 15 extant species in Eurasia. Despite being a small clade, they display a range of leaf morphologies comparable to the largest sections (>100 spp.) in Quercus. Their fossil record extends back to the Eocene. Here, we explore how cork oaks achieved their modern ranges and how legacy effects might explain niche evolution in modern species of section Cerris and its sister section Ilex, the holly oaks. METHODS: We inferred a dated phylogeny for cork and holly oaks using a reduced-representation next-generation sequencing method, restriction site-associated DNA sequencing (RAD-seq), and used D-statistics to investigate gene flow hypotheses. We estimated divergence times using a fossilized birth-death model calibrated with 47 fossils. We used Köppen profiles, selected bioclimatic parameters and forest biomes occupied by modern species to infer ancestral climatic and biotic niches. KEY RESULTS: East Asian and Western Eurasian cork oaks diverged initially in the Eocene. Subsequently, four Western Eurasian lineages (subsections) differentiated during the Oligocene and Miocene. Evolution of leaf size, form and texture was correlated, in part, with multiple transitions from ancestral humid temperate climates to mediterranean, arid and continental climates. Distantly related but ecologically similar species converged on similar leaf traits in the process. CONCLUSIONS: Originating in temperate (frost-free) biomes, Eocene to Oligocene ranges of the primarily deciduous cork oaks were restricted to higher latitudes (Siberia to north of Paratethys). Members of the evergreen holly oaks (section Ilex) also originated in temperate biomes but migrated southwards and south-westwards into then-(sub)tropical southern China and south-eastern Tibet during the Eocene, then westwards along existing pre-Himalayan mountain ranges. Divergent biogeographical histories and deep-time phylogenetic legacies (in cold and drought tolerance, nutrient storage and fire resistance) thus account for the modern species mosaic of Western Eurasian oak communities, which are composed of oaks belonging to four sections.

Corrigendum: ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids.

[This corrects the article DOI: 10.3389/fpls.2020.564414.].

ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids.

Hybridization and its relevance is a hot topic in ecology and evolutionary biology. Interspecific gene flow may play a key role in species adaptation to environmental change, as well as in the survival of endangered populations. Despite the fact that hybridization is quite common in plants, many hybridizing species, such as Quercus spp., maintain their integrity, while precise determination of genomic boundaries between species remains elusive. Novel high throughput sequencing techniques have opened up new perspectives in the comparative analysis of genomes and in the study of historical and current interspecific gene flow. In this work, we applied ddRADseq technique and developed an ad hoc bioinformatics pipeline for the study of ongoing hybridization between two relevant Mediterranean oaks, Q. ilex and Q. suber. We adopted a local scale approach, analyzing adult hybrids (sensu lato) identified in a mixed stand and their open-pollinated progenies. We have identified up to 9,435 markers across the genome and have estimated individual introgression levels in adults and seedlings. Estimated contribution of Q. suber to the genome is higher, on average, in hybrid progenies than in hybrid adults, suggesting preferential backcrossing with this parental species, maybe followed by selection during juvenile stages against individuals with higher Q. suber genomic contribution. Most discriminating markers seem to be scattered throughout the genome, suggesting that a large number of small genomic regions underlie boundaries between these species. A noticeable proportion of the markers (26%) showed allelic frequencies in adult hybrids very similar to one of the parental species, and very different from the other; a finding that seems relevant for understanding the hybridization process and the occurrence of adaptive introgression. Candidate marker databases developed in this study constitute a valuable resource to design large scale re-sequencing experiments in Mediterranean sclerophyllous oak species and could provide insight in species boundaries and on adaptive introgression between Q. suber and Q. ilex.

Genomic landscape of the global oak phylogeny.

The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny. We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 individuals representing nearly 250 Quercus species to infer a time-calibrated phylogeny of the world's oaks. We use a reversible-jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among-clade sampling biases. We then map the > 20 000 RAD-seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny. Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity. The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.

Comment on "Eocene Fagaceae from Patagonia and Gondwanan legacy in Asian rainforests".

Wilf et al (Research Articles, 7 June 2019, eaaw5139) claim that Castanopsis evolved in the Southern Hemisphere from where it spread to its modern distribution in Southeast Asia. However, extensive paleobotanical records of Antarctica and Australia lack evidence of any Fagaceae, and molecular patterns indicate shared biogeographic histories of Castanopsis, Castanea, Lithocarpus, and Quercus subgenus Cerris, making the southern route unlikely.

Prospects of barcoding the Italian wild dendroflora: oaks reveal severe limitations to tracking species identity.

DNA barcoding may be particularly important in influencing ecology, economic issues, and the fundamental crisis facing biodiversity as a standardized, species-level identification tool for taxonomy assessment. Trees play important roles in the conservation of many land ecosystems, the wood trade, and the definition of biogeographical processes; nevertheless, peculiar biological, evolutionary and taxonomical features will probably constitute an intriguing challenge to barcoders. We examined whether four marker regions (trnh-psba, rbcL, rpoc1, matK) proposed by the Consortium for the Barcode of Life (CBOL) matched species taxonomy in a preliminary tree biodiversity survey of Italian forested land. Our objective was to provide a test of future in situ applications of DNA barcodes by evaluating the efficacy of species discrimination under the criteria of uniformity of methods and natural co-occurrence of the species in the main forest ecosystems. Fifty-two species were included in a floristic study. We obtained 73% total discrimination success, with trnH-psbA as the best performing marker and oaks as the least responsive plants to the markers used. A further taxon-based study of Quercus (thirty specimens, 12 species) revealed that this genus is refractory to barcoding (0% discrimination success), a probable consequence of low variation rate at the plastid genome level, hybridization, and the incidence of biogeography. We conclude that some species-rich tree genera in small geographical regions may prove exceptionally difficult to barcode. Until more efficient markers are developed, we recommend that improved and diversified sampling (multiple locations of sympatric and co-occurring congenerics) be embraced as a timely and important goal for the precise assessment of haplotype specificity to facilitate the productive application of barcoding in practice.

Molecular characterization of new waxy mutants identified in bread and durum wheat.

Recently, electrophoretic analyses of waxy proteins in several hexaploid and tetraploid wheat accessions from worldwide collections have permitted the identification of new variants at the waxy loci, including allelic forms with different mobilities and partial null types. In this paper, the molecular characterization of mutated waxy loci in four bread wheat cultivars (two lacking the Wx-B1 and two lacking the Wx-D1 protein, respectively) and in four durum wheat cultivars (one lacking Wx-A1 and the remainder with Wx-B1 proteins showing different electrophoretic mobilities) was conducted by means of PCR, Southern and DNA sequence analyses. Three primer pairs were developed that identified six of the above-mentioned mutations and allowed their molecular description, providing a useful tool for further germplasm screening or marker assisted progeny selection in breeding programs involving the newly identified material. We have found that a complete gene deletion is responsible for a null allele at the Wx-B1 locus in one bread wheat line, whereas sequencing of the corresponding fragments showed a 724 bp deletion in the Wx-D1 locus in one line of bread wheat and an insertion of 89 bp in the Wx-A1 locus in one line of durum wheat, respectively. In addition, nucleotide substitutions and various insertions/deletions ranging from 3 to 30 bp were detected in the PCR fragments of one durum wheat line with a Wx-B1 protein with a different electrophoretic mobility. A fourth primer set, specific for this mutation, was consequently derived.

Utility of ITS sequence data for phylogenetic reconstruction of Italian Quercus spp.

Nuclear ribosomal DNA sequences encoding the 5.8S RNA and the flanking internal transcribed spacers (ITS1 and ITS2) were used to test the phylogenetic relationships within 12 Italian Quercus taxa (Fagaceae). Hypotheses of sequence orthology are tested by detailed inspection of some basic features of oak ITS sequences (i.e., general patterns of conserved domains, thermodynamic stability and predicted conformation of the secondary structure of transcripts) that also allowed more accurate sequence alignment. Analysis of ITS variation supported three monophyletic groups, corresponding to subg. Cerris, Schlerophyllodrys (=Ilex group sensu Nixon) and Quercus, as proposed by Schwarz [Feddes Rep., Sonderbeih. D, 1-200]. A derivation of the "Cerris group" from the "Ilex group" is suggested, with Q. cerris sister to the rest of the "Cerris group." Quercus pubescens was found to be sister to the rest of the "Quercus group." The status of hybrispecies of Q. crenata (Q. cerrisxQ. suber) and Q. morisii (Q. ilexxQ. suber) was evaluated and discussed. Finally, the phylogenetic position of the Italian species in a broader context of the genus is presented. The utility of the ITS marker to assess the molecular systematics of oaks is therefore confirmed. The importance of Italy as a region with a high degree of diversity at the population and genetic level is discussed.