Allele discovery of ten candidate drought-response genes in Austrian oak using a systematically informatics approach based on 454 amplicon sequencing
© Homolka et al.; licensee BioMed Central Ltd. 2012
Received: 26 September 2011
Accepted: 3 April 2012
Published: 3 April 2012
Rise of temperatures and shortening of available water as result of predicted climate change will impose significant pressure on long-lived forest tree species. Discovering allelic variation present in drought related genes of two Austrian oak species can be the key to understand mechanisms of natural selection and provide forestry with key tools to cope with future challenges.
In the present study we have used Roche 454 sequencing and developed a bioinformatic pipeline to process multiplexed tagged amplicons in order to identify single nucleotide polymorphisms and allelic sequences of ten candidate genes related to drought/osmotic stress from sessile oak (Quercus robur) and sessile oak (Q. petraea) individuals. Out of these, eight genes of 336 oak individuals growing in Austria have been detected with a total number of 158 polymorphic sites. Allele numbers ranged from ten to 52 with observed heterozygosity ranging from 0.115 to 0.640. All loci deviated from Hardy-Weinberg equilibrium and linkage disequilibrium was found among six combinations of loci.
We have characterized 183 alleles of drought related genes from oak species and detected first evidences of natural selection. Beside the potential for marker development, we have created an expandable bioinformatic pipeline for the analysis of next generation sequencing data.
White oaks are native to Europe, Asia, North Africa and North America, and include sessile oak, pedunculate oak, pubescent oak (Q. pubescens) and bur oak (Q. macrocarpa) among their most prominent species . About two percent of the Austrian forests harbour oak trees which correspond to an area of 66,000 hectare. Q. petraea and Q. robur are the predominant species while Q. cerris (Turkey oak) and Q. pubescens play only a minor role . In Austria as well as in Europe oak species colonize huge areas with vastly differing climatic conditions.
Rise of temperatures and shortening of available water as result of predicted climate change will impose significant pressures on long-lived forest tree species like European white oak. According to the latest predictions, we expect the average global surface temperature to increase by a maximum of 6.4°C within the next 90 years  leading to a higher frequency of severe drought events. But it is well known that different species diverge in their ability to resist drought induced damages and even within a species there is tremendous variability [4, 5]. Inter- as well as intraspecific allelic diversity is the key element of a plants potential to adapt to a changing environment and tolerance towards drought stress. There is an increased demand for molecular tools helping to describe these variations and to prepare forestry for future challenges.
Molecular markers are the first choice for plant research and breeding . Using them as landmarks, genetic maps can be established and subsequently used for identification of traits controlled by different genes (quantitative trait loci). Marker assisted selection provides breeders with an efficient tool for identifying desired phenotypes in large populations. Beside their use in breeding, molecular markers are highly valuable for population genetics permitting evolutionary studies and population structure interference.
Single nucleotide polymorphisms (SNPs) are commonly used for functional diversity assessment. Although they are highly abundant in the human genome and every SNP locus could potentially serve as utile marker, there are still few studies dealing with a high number of SNPs in plants . Latest advances in human and animal genome analysis have created several sophisticated technologies which are capable to analyze millions of SNPs in reasonable time and with low costs  and can be easily transferred to applications in plants. SNP discovery technologies include bioinformatic mining of expressed sequence tag (EST) databases , array based methods , comparison of whole genomes  and the application of next generation sequencing (NGS) for amplicon resequencing . Recently developed sequencing technologies summarized as next generation sequencing have already replaced traditional methods for detecting polymorphisms in genomes .
Roche 454 sequencing technology  is based on single strand amplification with emulsion polymerase chain reaction followed by pyrosequencing. Average readlengths of around 400 bp and high achievable coverage makes this technology well suited for discovery of SNPs and even detection of rare alleles. Short oligonucleotide barcodes can be used to tag individual sequences and enable the parallel analysis of several targets which has been successfully demonstrated in different species [14–16]. Multiplexing capabilities enable large studies including several hundred individuals with a high coverage for each sample. Therefore extensive cloning procedures to identify both haplotypes of diploid individuals as used in Sanger sequencing or generation of inbred lines  can be avoided. Bioinformatic haplotype interference with parsimony  or maximum-likelihood methods  is no longer necessary because each haplotype will be covered by a sufficient number of reads. The present paper describes the discovery and characterization of alleles in ten drought stress related genes originating from two oak species growing in Austria based on multiplex 454 amplicon sequencing and the development of a bioinformatic analysis pipeline.
Processing of 454 sequencing data
Comparison between raw data and preprocessed reads
Number of reads
Base pairs (Mbp)
Number of reads
Base pairs (Mbp)
Recovery statistics and individual coverage
Polymorphic sites and allele identification
Several studies describe the in silico detection of SNPs from the growing amount of EST data available [22, 23]. However this approach is likely to be biased by the small number of individuals coming from a limited amount of populations . To overcome this limitation we made use of high throughput sequencing techniques and a bioinformatic pipeline. 454 Sequencing technology developed by Roche has been chosen to analyze coding regions of drought responsive candidate genes. The main of advantage of this technique is the high number of reads generated leading to adequate coverage of each allele. Major drawbacks include the preferential amplification of short sequences and the high abundance of reading errors after homopolymer stretches. To cope with these problems we have developed a highly customizable pipeline which integrates preprocessing, statistics and error correction for 454 amplicon sequencing data and can be used to resolve multiplexed samples and to detect variants in a large number of sequences.
Different studies based on molecular markers have been conducted in oak. But the majority was investigating the population history of Quercus spp.  and hybridisation between Q. robur and Q. petraea. Only Derory et al. and Quang et al. used SNPs to assess allelic diversity of candidate genes. So far, this is the first study which aimed at detecting alleles of drought stress related genes of oak. With the aid of the bioinformatic pipeline presented above were able to detect a total amount of 183 alleles within eight genes of 336 oak individuals growing in Austria. The amount of alleles is comparable to results from white spruce where 173 alleles were detected in six loci from 283 individuals growing in Alaska . The average number of 6.71 mutations per 100 bp is higher or at least in the same range than in other plant species. Frequencies of 2.96 to 3.70, 3.85 or 3.83 to 7.30 have been discovered in maize , black poplar  and eucalyptus  which shows that the number of mutations strongly varies between species and gene analyzed. A high level of heterozygosity is present within ERD8 and LEA14 which points towards a large amount of genetic variability. Regarding the large number of individuals sampled across the heterogenic climatic area of Austria, a high variability is expected especially in genes like ERD8 and LEA14 which are directly involved in processes regulating drought response.
The high number of discovered alleles will be used as valuable source for association studies between allele frequencies and environmental variables in connection with drought stress. SNPs discovered in this study can be used to describe the genetic diversity present in two Austrian oak species and to develop molecular markers for drought tolerance if natural selection can be proven for some loci. Significant deviation from Hardy-Weinberg equilibrium as well as lower values of observed heterozygosity than expected at all loci may be the first evidence for natural selection. The deviations could partly be a result of sampling, be explained by the presence of null alleles (Table 3) or might arise from selective pressure on the coding regions. Although varying among genes, low levels of observed heterozygosity and an excess of homozygotes at RD26 support the latter hypothesis as well as the presence of linkage disequilibrium between some pairs of loci. This might either occur due to epistatic natural selection where favourable combinations of genes are linked and function together as supergenes  or several other factors including gene flow which is well documented among several oak species [27, 33, 34]. Present linkage between pairs of loci is the basis of association studies which can help to identify alleles occurring more frequently in plants exposed to dry conditions. Linkage between several loci of drought reactive genes is generally expected as drought response is under control of a huge regulatory network [35, 36]. However, this approach is limited by the fact that rapid decay of LD is commonly observed in forest tree species [37, 38] increasing the marker density needed for successful association studies.
Considering these findings, it points out clearly that the available dataset should be subject of further analysis like decay of LD, estimation of genetic diversity or differentiation of populations. To evaluate if selective forces were acting, different neutrality test as well as detection of outlier loci will be necessary. To evaluate if these deviations were only generated by demographic processes the results should be related to neutral markers like microsatellites found in chloroplasts (cpSSRs). If the presence of loci under selection and their association with climatic variables like temperature and precipitation could be proven, development of molecular markers will be possible. Given the predicted climate change and the resulting pressure on sessile organisms, functional markers might be a valuable tool for forestry and a basis for marker assisted selection of drought tolerant genotypes. Additional benefit arises from the fact that two different oak species were included in the analysis. Although they are mainly sympatric, Q. petraea is not as susceptible to drought induced damages as Q. robur. If species specific alleles can be discovered and related to environmental conditions the genetic basis of this benefit may be revealed.
Given next generation sequencing data of ten drought stress related genes originating from different oak species growing in Austria, we were able to discover new SNPs and characterize 183 alleles. In order to obtain these results, we have developed a semiautomatic analysis pipeline based on freely available tools and scripts. This pipeline can be fully automated and provided with a graphical user interface to make it more valuable for the scientific community. First analysis of the genetic data provides evidence of natural selection acting on the genes which makes them a target for future evolutionary studies and a potential source for molecular marker development. The alleles discovered will definitely help to understand drought adaptation processes acting in forest tree species.
Availability of supporting data
The raw data supporting the results of this article are available in the Dryad repository,http://dx.doi.org/10.5061/dryad.83gf113b. The data sets supporting the results of this article are included within the article and its additional files.
Plant material and DNA isolation
Summary data of candidate genes
Auxin responsive factor 16
Heat shock cognate protein 80
Germin-like protein 3A
Desiccation protectant protein Lea14
Responsive to desiccation 26
Amplification of candidate genes
Sequencing adaptors and barcodes were attached to the gene of interest following a modified two-step approach used by Schuelke . In the first step, genespecific primer pairs (internal primers) with M13-tails were used to amplify regions of interest. Internal primers targeting regions with a minimum of 200 and a maximum of 375 bp were planned with Primer3  using default settings. M13Fw (5′-TGTAAAACGACGGCCAGT-3′) and M13Re (5′-CAGGAAACAGCTATGACC-3′) were synthesized to the 5′-end of the internal primers. Amplicons for 454 sequencing were generated from 20 ng of template DNA. PCR reactions were performed in 25 μl total volume using 5 μl 5 x HOT FIREPol® Blend Master Mix 12.5 mM MgCl2 without dye (Solis BioDyne) and 3 pmol of each primer. A three-step PCR program consisting of 15 min. initial denaturation at 95°C followed by 32 cycles denaturation at 95°C for 30 sec., annealing at a temperature of 68°C for 45 sec., extension at 72°C for one minute and a final extension at 72°C for 8 minutes was used.
Annealing time (sec.)
Primer concentration (pmol)
Multiplexed 454 amplicon sequencing
PCR products of forty-two individuals per gene were pooled resulting in a total number of 80 pools. Cleaning was performed after pooling with the QIAquick PCR purification kit (QIAGEN) according to the manufacturer’s instructions. Concentrations of these pools were measured on Nanodrop (Thermo Scientific) and adjusted to 30 ng/μl adding TE buffer pH 8.0. PCR products of the same individuals originating from different genes were combined and sent for sequencing which was carried out by GATC Biotech in an 8 gasket format run on the Genome Sequencer FLX system (454 Life Sciences) with Titanium chemistry.
Bioinformatic analysis pipeline
To cope with errors resulting from the misinterpretation of homopolymer stretches by 454 sequencing we made use of the assembly tool SCARF. Sanger sequences  of the amplicons were used as references and homopolymer correction was turned on for a minimum length of 2 bp. Therefore homopolymer errors longer than the reference sequence will be trimmed to reduce read errors. Minimum percent identity and minimum hit score were set to 80 and 100, respectively (./scarf -f input.fasta –r referrence.fasta -c T -l 2 -p 80 -s 100). Alignment of the assembled reads was done with MultAlin using the AltDNA symbol comparison table and the gap penalty at extremities parameter set to end (ma -c: altdna.tab –x:1 input.clusters). Reads with gaps larger than 10 bp were masked in the alignment using a bash script for detection and not considered for further analysis.
SNP and allele detection
We created an extensive script (FindAlleles) using BioPerl which is able to identify mutations and detect alleles. Msf alignment files produced by MultAlin serve as input and are read using the Bio::AlignIO module. In a next step three different consensus sequences are calculated using the Bio::SimpleAlign module: consensus_iupac()—consensus sequence using IUPAC ambiguity codes for DNA, consensus_string()—standard consensus sequence displaying bases which occur in plurality and consensus_string(40)—produces a standard consensus and marks positions with a lower percent-identity than 40 % with “?”. A column representation of all sequence reads in the alignment is created. Using the three consensus sequences created, each column is screened for possible allelic distribution or arbitrary nucleotide insertions or deletions created by 454 sequencing technology and both types are annotated with a tag. In a next step the script checks each column for an allelic distribution of a 40/60 ratio. This correction threshold value was identified by manually examining randomly chosen alignments. In the examined alignments reading errors were not present in more than ten percent of the reads. The 40/60 distribution might be created by two different nucleotides (nucleotide variation) or by the insertion or deletion of a nucleotide (InDel). If this ratio is found the position is considered as valid “allelic” distribution. If it is not found, it is considered as reading error and is automatically corrected with the nucleotide found in the Sanger reference if a previously set correction threshold of ten percent is reached. As SCARF only treats homopolymer errors longer than the length present in the reference, we implemented an additional routine to fill up homopolymers shorter than the reference. As we observed errors occurring already in homopolymers consisting of two nucleotides, the min homopolymer parameter which sets the amount of nucleotides in a row to trigger correction was set to one. Nucleotide positions exceeding the reference are ignored.
After the correction phase the allele identification procedure is started running again over each column. If a tag for allelic distribution is found, the process splits up the available reads into two subgroups and a recursive subroutine is started for each split alignment. Each subroutine receives the associated reads and continues with them. Reads shorter than 70 % of the reference were discarded. Then the tags are checked and split up again if an allelic distribution is found. After all subroutines are finished alleles are exported to a multiple fasta file as well as statistics about the number of reads accounting for their creation and a text file with the excluded reads are generated. A cluster for allele generation was considered as valid if at least five reads covered more than 70 % of the reference sequence derived from Sanger sequencing. Number of alleles, observed and expected heterozygosities, deviation from Hardy-Weinberg equilibrium and linkage disequilibrium were calculated with Genepop 4.0 . Frequency of null alleles was estimated with FreeNA  and effective number of alleles was calculated with GenAlEx 6  using default settings.
The project was funded by the Austrian Federal Ministry of Agriculture, Forestry, Environment and Water Management (LE 2109 GZ LE.3.2.3/0014-IV 2/2005), the Regional Office of Vienna , the Regional Office of Lower Austria (BD1 NU94-2005, LF2-FO-7059/001-05, LF4-A-25/009-2004), the Regional Office of Upper Austria, the Regional Office of Styria (A3-36E4-05/2) and the Regional Office of Burgenland. We thank Dr. Schüler and the Federal Forest Office Vienna for sharing genomic DNA with us.
- Liesebach M, Schüler S, Weißenbacher L, Franner T, Henninger F, Geburek T: Herkunftsversuch mit Eiche in Vorbereitung. Forstzeitung. 2006, 117: 40-41.Google Scholar
- Dickson RE, Tomlinson PT: Oak growth, development and carbon metabolism in response to water stress. Ann Sci For. 1996, 53: 181-196. 10.1051/forest:19960202.View ArticleGoogle Scholar
- Zhenzhu X, Guangsheng Z, Hideyuki S: Plant response to drought and rewatering. Plant Signal Behav. 2010, 5: 649-654. 10.4161/psb.5.6.11398.View ArticleGoogle Scholar
- Beikircher B, Mayr S: Intraspecific differences in drought tolerance and acclimation in hydraulics ofLigustrum vulgareandViburnum lantana. Tree Physiol. 2009, 29: 767-775.View ArticleGoogle Scholar
- Luo N, Liu J, Yu X, Jiang Y: Natural variation of drought response inBrachypodium distachyon. Physiol Plant. 2011, 141: 19-29. 10.1111/j.1399-3054.2010.01413.x.PubMedView ArticleGoogle Scholar
- Ganal MW, Altmann T, Röder MS: SNP identification in crop plants. Curr Opin Plant Biol. 2009, 12: 211-217. 10.1016/j.pbi.2008.12.009.PubMedView ArticleGoogle Scholar
- Rafalski A: Applications of single nucleotide polymorphism in crop genetics. Curr Opin Plant Biol. 2002, 5: 94-100. 10.1016/S1369-5266(02)00240-6.PubMedView ArticleGoogle Scholar
- Le Dantec L, Chagné D, Pot D, Cantin O, Garnier-Géré P, Bedon F, Frigerio JM, Chaumeil P, Léger P, Garcia V, Laigret F, De Daruvar A, Plomion C: Automated SNP detection in expressed sequence tags: statistical considerations and application to maritime pine sequences. Plant Mol Biol. 2004, 54: 461-470.View ArticleGoogle Scholar
- Becker A, Chao DY, Zhang X, Salt DE, Baxter I: Bulk segregant analysis using single nucleotide polymorphism microarrays. PLoS One. 2011, 6: e15993-10.1371/journal.pone.0015993.PubMedPubMed CentralView ArticleGoogle Scholar
- Yao J, Lin H, Doddapaneni H, Civerolo EL: nWayComp: a genome-wide sequence comparison tool for multiple strains/species of phylogenetically related microorganisms. In Silico Biol. 2007, 7: 195-200.PubMedGoogle Scholar
- Meyer M, Stenzel U, Myles S, Prüfer K, Hofreiter M: Targeted high-throughput sequencing of tagged nucleic acid samples. Nucleic Acids Res. 2007, 35: e97-10.1093/nar/gkm566.PubMedPubMed CentralView ArticleGoogle Scholar
- Koboldt DC, Chen K, Wylie T, Larson DE, McLellan MD, Mardis ER, Weinstock GM, Wilson RK, Ding L: VarScan: variant detection in massively parallel sequencing of individual and pooled samples. Bioinformatics. 2009, 25: 2283-2285. 10.1093/bioinformatics/btp373.PubMedPubMed CentralView ArticleGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.PubMedPubMed CentralGoogle Scholar
- Stiller M, Knapp M, Stenzel U, Hofreiter M, Meyer M: Direct multiplex sequencing (DMPS)-a novel method for targeted high-throughput sequencing of ancient and highly degraded DNA. Genome Res. 2009, 10: 1843-1848.View ArticleGoogle Scholar
- Hamady M, Walker JJ, Harris CK, Gold NJ, Knight R: Error-correcting barcoded primers for pyrosequencing hundreds of samples in multiplex. Nat Methods. 2008, 5: 235-237. 10.1038/nmeth.1184.PubMedPubMed CentralView ArticleGoogle Scholar
- Külheim C, Yeoh SH, Maintz J, Foley WJ, Moran GF: Comparative SNP diversity among fourEucalyptusspecies for genes from secondary metabolite biosynthetic pathways. BMC Genomics. 2009, 10: 452-10.1186/1471-2164-10-452.PubMedPubMed CentralView ArticleGoogle Scholar
- Neigenfind J, Gyetvai G, Basekow R, Diehl S, Achenbach U, Gebhardt C, Selbig J: Haplotype inference from unphased SNP data in heterozygous polyploids based on SAT. BMC Genomics. 2008, 9: 356-10.1186/1471-2164-9-356.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang IL, Chang CY: Mathematical properties and bounds on haplotyping populations by pure parsimony. Math Biosci. 2011, 231: 120-125. 10.1016/j.mbs.2011.02.008.PubMedView ArticleGoogle Scholar
- Kirkpatrick BB: Haplotypes versus genotypes on pedigrees. Algorithms Mol Biol. 2011, 6: 10-10.1186/1748-7188-6-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Barker MS, Dlugosch KM, Reddy ACC, Amyotte SN, Rieseberg LH: SCARF: maximizing next-generation EST assemblies for evolutionary and population genomic analyses. Bioinformatics. 2009, 25: 535-536. 10.1093/bioinformatics/btp011.PubMedView ArticleGoogle Scholar
- Spieß N, Oufir M, Matušíková I, Stierschneider M, Kopecky D, Homolka A, Burg K, Fluch S, Hausman JF, Wilhelm E: Ecophysiological and transcriptomic responses of oak (Quercus robur) to long-term drought exposure and rewatering. Environ Exp Bot. 2012, 77: 117-126.View ArticleGoogle Scholar
- Neale DB, Kremer A: Forest tree genomics: growing resources and applications. Nat Rev Genet. 2011, 12: 111-122.PubMedView ArticleGoogle Scholar
- Lepoittevin C, Harvengt L, Plomion C, Garnier-Géré P: Association mapping for growth, straightness and wood chemistry traits in thePinus pinasterAquitaine breeding population. Tree Genet Genomes. 2011, 8: 113-126.View ArticleGoogle Scholar
- González-Martínez SC, Krutovsky KV, Neale DB: Forest-tree population genomics and adaptive evolution. New Phytol. 2006, 170: 227-238. 10.1111/j.1469-8137.2006.01686.x.PubMedView ArticleGoogle Scholar
- Petit RJ, Csaikl UM, Bordács S, Burg K, Coart E, Cottrell J, Van Dam B, Deans JD, Dumolin-Lapègue S, Fineschi S, et al: Chloroplast DNA variation in European white oaks. Phylogeography and patterns of diversity based on data from over 2600 populations. Forest Ecol Manag. 2002, 156: 5-26. 10.1016/S0378-1127(01)00645-4.View ArticleGoogle Scholar
- Scotti-Saintagne C, Mariette S, Porth I, Goicoechea PG, Barreneche T, Bodénès C, Burg K, Kremer A: Genome scanning for interspecific differentiation between two closely related oak species [Quercus roburL. andQ. petraea(Matt.) Liebl.]. Genetics. 2004, 168: 1615-1626. 10.1534/genetics.104.026849.PubMedPubMed CentralView ArticleGoogle Scholar
- Derory J, Scotti-Saintagne C, Bertocchi E, Le Dantec L, Graignic N, Jauffres A, Casasoli M, Chancerel E, Bodénès C, Alberto F, Kremer A: Contrasting relations between diversity of candidate genes and variation of bud burst in natural and segregating populations of European oaks. Heredity. 2010, 105: 401-411. 10.1038/hdy.2009.170.PubMedView ArticleGoogle Scholar
- Quang ND, Ikeda S, Harada K: Nucleotide variation inQuercus crispulaBlume. Heredity. 2008, 101: 166-174. 10.1038/hdy.2008.42.PubMedView ArticleGoogle Scholar
- Anderson LL, Hu FS, Paige KN: Phylogeographic history of white spruce during the last glacial maximum: uncovering cryptic refugia. J Hered. 2011, 102: 207-216. 10.1093/jhered/esq110.PubMedView ArticleGoogle Scholar
- Ching A, Caldwell KS, Jung M, Dolan M, Smith OS, Tingey S, Morgante M, Rafalski AJ: SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines. BMC Genet. 2002, 3: 19-PubMedPubMed CentralView ArticleGoogle Scholar
- Chu Y, Su X, Huang Q, Zhang X: Patterns of DNA sequence variation at candidate gene loci in black poplar (Populus nigraL.) as revealed by single nucleotide polymorphism. Genetica. 2009, 137: 141-150. 10.1007/s10709-009-9371-1.PubMedView ArticleGoogle Scholar
- Joron M, Frezal L, Jones RT, Chamberlain NL, Lee SF, Haag CR, Whibley A, Becuwe M, Baxter SW, Ferguson L, Wilkinson PA, Salazar C, Davidson C, Clark R, Quail MA, Beasley H, Glithero R, Lloyd C, Sims S, Jones MC, Rogers J, Jiggins CD, Ffrench-Constant RH: Chromosomal rearrangements maintain a polymorphic supergene controlling butterfly mimicry. Nature. 2011, 477: 203-206. 10.1038/nature10341.PubMedPubMed CentralView ArticleGoogle Scholar
- Whittemore AT, Schaal B: Interspecific gene flow in sympatric oaks. Proc Nat Acad Sci USA. 1991, 88: 2540-2544. 10.1073/pnas.88.6.2540.PubMedPubMed CentralView ArticleGoogle Scholar
- Valbuena-Carabana M, González-Martínez SC, Sork VL, Collada C, Soto A, Goicoechea PG, Gil L: Gene flow and hybridization in a mixed oak forest (Quercus pyrenaicaWilld. andQuercus petraea(Matts.) Liebl.) in central Spain. Heredity. 2005, 95: 457-465. 10.1038/sj.hdy.6800752.PubMedView ArticleGoogle Scholar
- Thapa G, Dey M, Sahoo L, Panda SK: An insight into the drought stress induced alterations in plants. Biol Plantarum. 2011, 55: 603-613. 10.1007/s10535-011-0158-8.View ArticleGoogle Scholar
- Lorenz WW, Alba R, Yu Y, Bordeaux JM, Simões M, Dean JFD: Microarray analysis and scale-free gene networks identify candidate regulators in drought-stressed roots of loblolly pine (P. taedaL.). BMC Genomics. 2011, 12: 264-10.1186/1471-2164-12-264.PubMedPubMed CentralView ArticleGoogle Scholar
- González-Martínez SC, Ersoz E, Brown GR, Wheeler NC, Neale DB: DNA sequence variation and selection of tag single-nucleotide polymorphisms at candidate genes for drought-stress response inPinus taedaL. Genetics. 2006, 172: 1915-1926.PubMedPubMed CentralView ArticleGoogle Scholar
- Eckert AJ, Wegrzyn JL, Pande B, Jermstad KD, Lee JM, Liechty JD, Tearse BR, Krutovsky KK, Neale DB: Multilocus patterns of nucleotide diversity and divergence reveal positive selection at candidate genes related to cold hardiness in coastal Douglas fir (Pseudotsuga menziesiivar.menziesii). Genetics. 2009, 183: 289-298. 10.1534/genetics.109.103895.PubMedPubMed CentralView ArticleGoogle Scholar
- Csaikl UM, Burg K, Fineschi S, König AO, Matyas C, Petit RJ: Chloroplast DNA variation of white oaks in the alpine region. Forest Ecol Manag. 2002, 156: 131-1454. 10.1016/S0378-1127(01)00639-9.View ArticleGoogle Scholar
- Schüler S, Weißenbacher L: Herkunftsversuche mit Stiel- und Traubeneiche aus Österreich und angrenzenden Ländern = Provenance trials with seed sources of pendunculate and sessile oak originating from Austria and neighboring countries. BFW-Dokumentation, Wien. 2010, 13: 40-Google Scholar
- Schuelke M: An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000, 18: 233-234. 10.1038/72708.PubMedView ArticleGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Humana Press, Totowa, 365-386.Google Scholar
- Ihaka R, Gentleman R: R: a language for data analysis and graphics. J Comput Graph Stat. 1996, 5: 399-414.Google Scholar
- Corpet F: Multiple sequence alignment with hierarchical clustering. Nucl Acids Res. 1988, 16: 10881-10890. 10.1093/nar/16.22.10881.PubMedPubMed CentralView ArticleGoogle Scholar
- Stajich J, Block D, Boulez K, Brenner S, Chervitz S, Dagdigian C, Fuellen G, Gilbert JG, Korf I, Lapp H, Lehväslaiho H, Matsalla C, Mungall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E: The Bioperl toolkit: perl modules for the life sciences. Genome Res. 2002, 12: 1611-1618. 10.1101/gr.361602.PubMedPubMed CentralView ArticleGoogle Scholar
- Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Heredity. 1995, 86: 248-249.Google Scholar
- Chapuis MP, Estoup A: Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007, 24: 621-631.PubMedView ArticleGoogle Scholar
- Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar