Skip to main content


  • Poster presentation
  • Open Access

Genetic divergence in Cork Oak based on cpDNA sequence data

  • 1,
  • 2Email author,
  • 3,
  • 2,
  • 4,
  • 4,
  • 5,
  • 2,
  • 1 and
  • 6
BMC Proceedings20115 (Suppl 7) :P13

  • Published:


  • Introgressed Lineage
  • Intergenic Space Region
  • cpDNA Haplotype
  • Western Mediterranean Basin
  • cpDNA Region

Background and objectives

Cork oak (Quercus suber L.) is one of the dominant broadleaved woody species in the western Mediterranean Basin, defining unique open woods. These woodlands have an outstanding economical and ecological value in this region, particularly in Portugal, where they sustain a strong cork industry. In the context of a prospective management of these sustainable ecosystems and renowned reservoirs of biodiversity, it is vital to better understand how the genetic variation of Q. suber natural populations is spatially organized so reasonable guidelines for conservation can be provided. On the other hand, knowledge of how past climate fluctuations influenced the patterns and dynamics of Q. suber, shaping the ranges of the species in the Mediterranean peninsulas, is of the utmost importance for our perception of what can happen in the future. Although a great deal of details on the genetic divergence of the Mediterranean cork oak populations has been uncovered and several hypotheses have been advanced concerning its evolutionary history, there is still much to unravel. For instance, Portuguese natural population sampling included so far in previous studies has been very deficient. To achieve this goal a different and complementary analysis of cork oak´s genetic diversity was initiated under a phylogeographical framework based on chloroplastidial DNA sequences. This study is the starting up of a project aiming at assessing the genetic diversity and differentiation of natural cork oak populations from the entire Mediterranean distribution, with the intent of understanding patterns of biodiversity, gene flow and population admixture, as well as to infer possible evolutionary events.

Materials and methods

We used 3-5 samples from 25 populations collected across the range of distribution of cork oak, in a total of 115 individuals. Three chloroplastidial DNA (cpDNA) intergenic spacer regions (TrnL-F [1], TrnS-PsbC [2] and TrnH-psbA [3]) were amplified and sequenced. Fourteen individuals of Q. rotundifolia, five of Q. coccifera, two of Q. ilex, and one of Q. robur, Q. faginea, Q. pyrenaica, Q. lusitanicus, Q. rubra and Q. canariensis were also analyzed as comparative references. Castanea crenata was used as outgroup. The alignment was developed in Clustal X 2.0.12 [4] and manually refined. The phylogenetic trees were made with PAUP 4.0d99 [5], using a Maximum Parsimony analysis method.


Regarding the cpDNA regions under study, two main lineages of cork oak haplotypes were found: one lineage closely related to Quercus rotundifolia and Q. coccifera, (“introgressed lineage”), mainly present in Iberia and Morocco; and a second lineage that contrastingly seems “pure”, which is the most common and does not appear to be shared with any other Quercus species (“pure lineage”). Three distinct sublineages are shown in the “pure” lineage, corresponding to the Eastern populations, Sicilly and the Western populations, although separated by minor differences. Some populations present cpDNA haplotypes belonging to both lineages, while others are specific for one of them. An absence of introgression in the easternmost populations of cork oak seems also apparent.


Although these preliminary results seem to confirm previous ones, some major differences are suggested, such as too few differences between the western and eastern groups to be explained by a Tertiary divergence pattern, as previously suggested [6], but rather more consistent with a more recent expansion from few refugia; and widespread and multiple introgression events from Q. rotundifolia and Q. coccifera, common in peripheral western populations.



Fundação para a Ciência e a Tecnologia (FCT) is acknowledged for financial support, through project PTDC/AGR-GPL/104966/2008.

Authors’ Affiliations

Faculdade de Ciências da Universidade de Lisboa/Centro de Biologia Ambiental (CBA), Computational Biology and Population Genomics Group, Edifício C2, Campus da Faculdade de Ciências, 1759-016 Lisboa, Portugal
Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa (ITQB-UNL)/Instituto de Biologia Experimental e Tecnológica (IBET), Av. da República, 2780-501 Oeiras, Portugal
Centro de Estudos Florestais, Tapada da Ajuda, Universidade Técnica de Lisboa, Instituto Superior de Agronomia (ISA), 1349-017 Lisboa, Portugal
Instituto Nacional de Recursos Biológicos, I.P. (INRB), L-INIA - Pólo do Lumiar, Unidade Investigação Recursos Genéticos, Ecofisiologia e Melhoramento de Plantas, Grupo Biologia Molecular, Estrada Paço Lumiar 22, Ed.S 1°, 1649-038 Lisboa, Portugal
INRB, I.P., L-INIA - Pólo de Oeiras, Unidade Recursos Genéticos, Ecofisiologia e Melhoramento de Plantas, Quinta do Marquês, 2784-505 Oeiras, Portugal
Instituto de Investigação Científica Tropical IICT/CIFC, Quinta do Marquês, 2784-505 Oeiras, Portugal


  1. Taberlet P, Gielly L, Pautou G, Bouvet J: Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Mol Biol. 1991, 17: 1105-1109. 10.1007/BF00037152.View ArticlePubMedGoogle Scholar
  2. Nishizawa T, Watano Y: Primer pairs suitable for PCR-SSCP analysis of chloroplast DNA in angiosperms. J Phytogeo Taxon. 2000, 48 (1): 67-70.Google Scholar
  3. Kress WJ, Wurdack KJ, Zimmer EA, Weigt LA, Janzen DA: Use of DNA barcodes to identify flowering plants. PNAS. 2005, 102 (23): 8369-8374. 10.1073/pnas.0503123102.PubMed CentralView ArticlePubMedGoogle Scholar
  4. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
  5. Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 2002, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
  6. Magri D, Fineschi S, Bellarosa R, Buonamici A, Sebastiani F, Schirone B, Simeone MC, Vendramin GG: The distribution of Quercus suber chloroplast haplotypes matches the palaeogeographical history of the western Mediterranean. Mol Ecol. 2007, 16: 5259-5266. 10.1111/j.1365-294X.2007.03587.x.View ArticlePubMedGoogle Scholar


© Costa et al; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.