Adaptive Variation and Introgression of a CONSTANS-Like Gene in North American Red Oaks

Oaks provide a model system to study maintenance of species identity by divergent selection since they maintain morphological differences and ecological adaptations despite interspecific hybridization. The genome of closely related interfertile oak species was shown to be largely homogeneous, with a few genomic areas exhibiting high interspecific differentiation possibly as result of strong divergent selection. Previously, a genic microsatellite was identified as under strong divergent selection, being nearly fixed on alternative alleles in the two interfertile North American red oak species: Quercus rubra L. and Quercus ellipsoidalis E.J. Hill. Further genotyping in two other red oak species—Quercus velutina Lam. and Quercus coccinea Munchh.—revealed a similar bias for the Q. ellipsoidalis-specific allele. To further elucidate the basis of this differentiation, we sequenced the microsatellite in individuals from all four red oak species. Sequence variability was observed in the microsatellite motif which encodes a poly-Q repeat in a COL gene involved in phenology and growth. Furthermore, in neighboring (parapatric) Q. rubra/Q. ellipsoidalis populations, introgression of the Q. ellipsoidalis-specific allele into Q. rubra occurred at a lower rate than introgression of the Q. rubra-specific allele into Q. ellipsoidalis despite symmetric interspecific gene flow, indicating potential adaptive introgression. Introgression of adaptive alleles can be an important mechanism for rapid adaptation to new environmental conditions (e.g., climate change).


Introduction
The transfer of adaptive alleles and traits by hybridization might be an important mechanism of rapid evolution and adaptation to changing environments, e.g., in the face of climate change, and evidence for adaptive trait transfer has been reported in both plant and animal species (reviewed in [1,2]).The availability of genomic resources and analytical methods for the identification of loci under strong divergent selection (outlier loci) [3] allows us to trace the introgression of potentially adaptive alleles in interspecific hybrid zones.Oaks generally reveal porous species boundaries, but morphological species' identity and ecological adaptations (e.g., soil moisture) are generally maintained despite recurrent interspecific gene flow [4][5][6].Thus, oaks provide a model for the identification of outlier loci under divergent selection with annotated functions and potential roles in stress tolerance and reproductive isolation between species (e.g., [7]).
The four interfertile North American red oak species-Quercus rubra L., Quercus velutina Lam., Quercus coccinea Münchh., and Quercus ellipsoidalis E.J. Hill-exhibit porous species boundaries and recent studies found strong evidence for contemporary interspecific gene flow in sympatric and parapatric stands based on genetic assignment and parentage analyses [5,[8][9][10].Most recently, Owusu et al. used both genic and nuclear microsatellite markers to better resolve the taxonomic relationships between these four species [10].The genetic assignment analysis showed a clear separation of all four species, with Q. rubra being the most differentiated.After excluding genetically intermediate individuals identified in the genetic assignment analysis, a phylogenetic tree based on population distances at nuclear microsatellite markers revealed Q. ellipsoidalis and Q. coccinea as most closely related and a clear separation of Q. rubra from the Q. velutina/Q.ellipsoidalis/Q.coccinea clade [10].This result is also supported by previous studies and a recent restriction site associated sequencing (RAD-seq) phylogeny of the genus ( [9,[11][12][13]).Genetically intermediate forms were found in contact zones between Q. ellipsoidalis, Q. rubra and Q. velutina, in sympatric stands indicating gene flow among these three species [8,10,14].High levels of interspecific gene flow were detected especially between Q. ellipsoidalis and Q. velutina.Thus about 20% of Q. velutina individuals had a recent Q. ellipsoidalis ancestor and about 30% of Q. ellipsoidalis individuals had a recent Q. velutina ancestor [8].This supports previous indications of recent and ongoing hybridization between Q. ellipsoidalis and Q. velutina [9,11].Recent studies have shown gene flow between Q. rubra and the other three species to be more restricted, yet generally symmetric [9,10,15].These four species also represent a drought tolerance gradient, with Q. rubra as the least drought tolerant followed by Q. velutina, Q. coccinea and lastly, Q. ellipsoidalis as the most drought tolerant [16,17].
Previously, we discovered among 36 genic and eight non-genic microsatellites one genic microsatellite (FIR013), located in the coding region of a CONSTANS-like (COL) gene, which was nearly fixed on alternative alleles in multiple population pairs of two of these species (Q.rubra and Q. ellipsoidalis; [7]).This genic microsatellite was originally developed for Quercus robur L. and the trinucleotide microsatellite encodes a poly-Q-repeat (poly-glutamine repeat) [18].Quercus ellipsoidalis is characterized by the lack of one repeat unit (138 base pair (bp) allele) as compared to Q. rubra (141 bp allele, [7]).The locus has a putative function as a COL gene, which is thought to be involved in flowering time and growth [7], both of which can be impacted by water availability [19].In the European white oak species Quercus petraea (Matt.)Liebl., a single nucleotide polymorphism (SNP) located in the same COL gene, was identified as significantly associated with bud burst along an altitudinal gradient [20].Latitudinal and altitudinal gradients, as well as local environmental conditions such as water availability, have been shown to impact bud burst timing in other oak species [21,22].Quercus rubra and Q. ellipsoidalis seedlings grown in a common garden exhibited differences in bud burst timing and leaf fall over two consecutive years [23] although this did not hold in natural populations of the same provenance [24].Thus, differences in flowering time could still contribute to limit gene flow between these two species [24,25] and highly divergent markers, such as FIR013, located in the first exon of a COL gene, might be involved in their adaptive divergence and/or partial reproductive isolation.COL genes have been linked to adaptive divergence in other species, including the European white oak Q. petraea and Populus tremula L. [20,26,27].For example, in P. tremula, an allele of the coding poly-E repeat (poly-glutamic acid) microsatellite in COL2B, was associated with growth cessation across a latitudinal gradient [26].Interestingly, Q. ellipsoidalis exhibited later leaf fall than Q. rubra in a common garden, indicating genetic differences between species in leaf fall [22], a trait which is positively associated with growth cessation [28].The biological function of single amino acid repeats (e.g., poly-E and poly-Q) has been studied mostly in animal species, with variation in these repeats often associated with genetic disorders [29,30].However, variation in poly-glutamine repeats, such as the poly-Q repeat found in the COL gene, have been shown to be under selection in various species including fish, birds, plants and fungi, making this particular marker an excellent candidate gene, potentially contributing to reproductive isolation and adaptive divergence between the two oak species [31].
In this study, we asked two questions: (1) is the microsatellite size variation due to poly-Q repeat variation in these four red oak species; and (2) is there introgression of potentially adaptive alleles between drought averse Q. rubra and drought tolerant Q. ellipsoidalis, Q. velutina and Q. coccinea?
To answer the first question, we have sequenced part of the COL gene including the poly-Q repeat in a total of 46 samples representing each species and allele at FIR013.In reference to the second question, we assessed genetic variation in 16 populations (Table S1, Figure 1) at 12 microsatellite markers (including FIR013) and examined relative allelic and genotypic frequencies in each genetically assigned species at FIR013.

Materials and Methods
DNA samples were obtained from 16 populations of four red oak species (Q.rubra, Q. ellipsoidalis, Q. velutina, and Q. coccinea) from the Great Lakes region (Wisconsin, Michigan, Illinois and Indiana; see Table S1 and Figure 1) [8,10,25].The samples consisted of eight population pairs of Q. rubra and Q. ellipsoidalis (parapatric populations), six mixed stands of Q. ellipsoidalis and Q. velutina in proximity to unsampled Q. rubra populations (parapatric populations) and one mixed stand of Q. ellipsoidalis, Q. velutina, and Q. rubra (sympatric population).Finally, one Q.coccinea population within the distribution range of Q. rubra and Q. velutina was also included.Leaf material was stored at −20 • C until DNA extraction using the DNeasy96 Plant Kit (Qiagen, Valencia, CA, USA) following the manufacturer's instructions.
Forests 2017, 8, 322 3 of 13 total of 46 samples representing each species and allele at FIR013.In reference to the second question, we assessed genetic variation in 16 populations (Table S1, Figure 1) at 12 microsatellite markers (including FIR013) and examined relative allelic and genotypic frequencies in each genetically assigned species at FIR013.

Materials and Methods
DNA samples were obtained from 16 populations of four red oak species (Q.rubra, Q. ellipsoidalis, Q. velutina, and Q. coccinea) from the Great Lakes region (Wisconsin, Michigan, Illinois and Indiana; see Table S1 and Figure 1) [8,10,25].The samples consisted of eight population pairs of Q. rubra and Q. ellipsoidalis (parapatric populations), six mixed stands of Q. ellipsoidalis and Q. velutina in proximity to unsampled Q. rubra populations (parapatric populations) and one mixed stand of Q. ellipsoidalis, Q. velutina, and Q. rubra (sympatric population).Finally, one Q.coccinea population within the distribution range of Q. rubra and Q. velutina was also included.Leaf material was stored at −20 °C until DNA extraction using the DNeasy96 Plant Kit (Qiagen, Valencia, CA, USA) following the manufacturer's instructions.S1; approximate species ranges are from Little (1971).
Using the primer pair for microsatellite FIR013 (Table S2), we sequenced the microsatellite and flanking regions in 46 samples (all homozygous for either the 138 bp, 141 bp or 144 bp allele at FIR013) from sampled populations of the four red oak species: Q. ellipsoidalis (11), Q. rubra (9), Q. velutina (18), and Q. coccinea (8) (see Table S3).Sequencing reactions were run on an ABI3730 DNA Analyzer using the ABI BigDye Terminator Cycle Sequencing Ready Reaction Kit v.3.1 (Applied Biosystems, Foster City, CA, USA) at the Nevada Genomics Center at the University of Nevada in Reno.Sequence data were aligned using the BioEdit Sequence Alignment Editor [32].Alignment of forward and reverse sequences was completed by using pairwise alignment with the option to allow the two ends to slide.Multiple fragment alignments were completed using the Muscle algorithm (with default settings) in MEGA (Molecular Evolutionary Genetics Analysis) 6.06 [33] and final alignments were generated after careful visual inspection and manual re-editing.Sequence data were translated into amino acid sequences in MEGA using the standard genetic code option.
All individuals from the 16 sampled populations were previously genotyped at 11 microsatellite markers (six expressed sequence tag-simple sequence repeats (EST-SSRs) and five nuclear simple sequence repeats (nSSRs); Table S2).Genotypic data on EST-SSR FIR013 were available for the eight Q. rubra and Q. ellipsoidalis populations, and we genotyped the remaining seven populations (six sympatric Q. ellipsoidalis/Q.velutina populations and the sympatric Q. rubra/Q.ellipsoidalis/ Q. velutina population; Table S2) at FIR013 [8,10,25].A single population of Q. coccinea was included for the sequencing effort and for the amplification of FIR013, where eight samples were sequenced and all 20 samples were genotyped at FIR013.
All genetic variation analyses were conducted both with and without the outlier FIR013.Genetic variation parameters were calculated in GeneAlEx 6.41 [34] including the number of alleles per locus (N a ), Nei's unbiased gene diversity (H e ), observed heterozygosity (H o ) and Wright's inbreeding coefficient (F) by population and species.Using only pure individuals from each species (based on previous genetic assignment analyses exluding outlier loci [8,10,25]), F ST and pairwise F ST between species with corresponding significances were calculated in GenePop 4.1 [35].F ST -based outlier screens were conducted in the program LOSITAN [36] between Q. rubra and Q. ellipsoidalis, Q. rubra and Q. velutina, and Q. velutina and Q. ellipsoidalis using the same settings as detailed in [7].Specifically, we applied the stepwise mutation model with 50,000 simulations at the 99% confidence level and false discovery rate of 0.10.
Relative frequencies of introgressive forms and hybrids between drought averse Q. rubra and the more drought tolerant Q. ellipsoidalis, Q. velutina, and Q. coccinea were calculated based on results from two previous genetic assignment analyses [10,25].Although slightly different marker sets were used in the genetic assignments, we have shown previously that very similar results are obtained using marker sets of 10, 16 or 44 markers [7,24].The same method and criteria for genetic assignment analysis in the program STRUCTURE [37] were used as described in [25] and did not include the marker FIR013 or any other outlier marker.Pure species, introgressive forms and hybrids were classified as having >0.90, 0.61-0.89,and 0.4-0.6 proportion of their ancestry in one genetic cluster, respectively.Introgression of the species characteristic FIR013 alleles 138 bp (Q.ellipsoidalis) and 141 bp (Q.rubra) [7] between species with different drought tolerance was determined by relative frequencies of the 141 bp allele in genetically assigned Q. ellipsoidalis, Q. velutina and Q. coccinea individuals and as the relative frequency of the 138 bp allele in genetically assigned Q. rubra individuals.

Nucleotide Sequence Analysis
Sequence variation confirmed that allele size differences were due to variation in a poly-Q repeat (Figure 2).Quercus velutina and Q. coccinea share the same common allele (138 bp) with Q. ellipsoidalis (Figure 3) which is reflected in the sequence variation between these three species and Q. rubra (Figure 2).In addition, there was a non-synonymous single nucleotide polymorphism (SNP) within the nucleotide repeat region present in all four species at position 9 (Figure 2).This SNP displayed different allelic patterns in Q. rubra as compared to the other three more drought tolerant species, with almost all Q. rubra individuals homozygous for the nucleotide "C" (88.9%), while heterozygotes (C/G) are only present in Q. velutina (44.4%),Q. ellipsoidalis (18.2%) and Q. coccinea (50.0%).The frequency of homozygotes (G/G) was similar (between 10%-20%) for all but Q. ellipsoidalis (45.5%) (Table S4).Sample sizes were small for all four species and may bias these results.
(Figure 3) which is reflected in the sequence variation between these three species and Q. rubra (Figure 2).In addition, there was a non-synonymous single nucleotide polymorphism (SNP) within the nucleotide repeat region present in all four species at position 9 (Figure 2).This SNP displayed different allelic patterns in Q. rubra as compared to the other three more drought tolerant species, with almost all Q. rubra individuals homozygous for the nucleotide "C" (88.9%), while heterozygotes (C/G) are only present in Q. velutina (44.4%),Q. ellipsoidalis (18.2%) and Q. coccinea (50.0%).The frequency of homozygotes (G/G) was similar (between 10%-20%) for all but Q. ellipsoidalis (45.5%) (Table S4).Sample sizes were small for all four species and may bias these results.

Genetic Diversity and Differentiation Analysis
Levels of genetic diversity were similar across populations and species (Table S6).Quercus rubra showed significant, but low differentiation from the other three species at the 11 non-outlier microsatellites (3% to 4.3%) (Table 1), and very high and significant differentiation at the FIR013 locus (62.5% to 67.2%, Table 1).Genetic diversity and differentiation at all 12 microsatellite markers showed

Genetic Diversity and Differentiation Analysis
Levels of genetic diversity were similar across populations and species (Table S6).Quercus rubra showed significant, but low differentiation from the other three species at the 11 non-outlier microsatellites (3% to 4.3%) (Table 1), and very high and significant differentiation at the FIR013 locus (62.5% to 67.2%, Table 1).Genetic diversity and differentiation at all 12 microsatellite markers showed similar patterns (Tables S5 and S7).Differentiation at FIR013 between the drought tolerant species was low at 0.9% (Table 1).All outlier screens between Q. rubra and other more drought tolerant species identified FIR013 as an outlier, consistent with the allelic patterns observed for all four species (Table S8, Figures S2-S4).In the only Q. coccinea population, the 141 bp (Q.rubra) allele was rare (3%) and no homozygotes for this allele were found (Table 2; Figure 3).The frequency of the 141 bp allele in Q. velutina and Q. ellipsoidalis was higher (22% and 15%, respectively) and homozygotes for the 141 bp allele were observed in both species (6% in Q. velutina, 5% in Q. ellipsoidalis) (Figure 3).
For sympatric Q. velutina/Q.ellipsoidalis populations that were in close proximity to Q. rubra populations, introgression of the 141 bp (Q.rubra) allele was considerably higher into Q.velutina than into Q.ellipsoidalis for three out of the five populations (Figure 4; Table 2: populations MI-NC, MI-OGC, WRPNP-IL-CC).The level of introgression of the 141 bp allele was similar for Q. velutina and Q. ellipsoidalis in population LS-IN-LC and higher in Q. ellipsoidalis than in Q. velutina in population HPNSP-IN-LC (Figure 4; Table 2).

Discussion
Quercus rubra is the only species nearly fixed on the 141 bp allele at outlier locus FIR013.The other three species are largely fixed on the 138 bp allele, and this difference is reflected in the sequence variation of the microsatellite itself.Previous work on these four species has also indicated that Q. rubra is the most diverged of the four species and that there is contemporary gene flow occurring between all four of these species [8,10].The lack of one glutamine residue in the poly-Q repeat is a shared character for the three more drought tolerant species, with interspecific gene flow likely responsible for incomplete differentiation.The phylogeny of section Lobatae is in progress [13] and can be used to evaluate whether the character is a synapomorphy for the three more drought tolerant species.The four species represent a drought tolerance gradient with Q. ellipsoidalis being the most drought tolerant, followed by Q. coccinea and then Q. velutina, while Q. rubra is considered the most mesophilic red oak species [16,38].Despite morphological and ecophysiological differences [8,16,38,39] between the four oak species, genetic differentiation between them is low at most nuclear genetic markers [7,8,10,25] as well as at chloroplast markers [14].This is consistent with the tendency for closely related oak species to frequently hybridize with one other [40].Thus, genomic regions that display high levels of genetic differentiation are of great interest to elucidate how oaks maintain their species identity despite interspecific gene flow.
Poly-Q repeat function and variation have been studied extensively in regards to human disease, but there is scant information about their function in plant and animal species [29][30][31].Thus, further investigation of their functional role and variation is warranted.In addition to the repeat number variation, a non-synonymous SNP within the poly-Q repeat (Figure 2) was found in all four species (C/G) and there was a bias towards both the heterozygotes (C/G) and homozygotes (G/G) in the three more drought tolerant species.The observed SNP change from C to G results in a change of the protein sequence from a glutamine (uncharged polar) to a histidine (positive polar).These sequence differences between species with different drought tolerance might be related to functional differences.Association studies in full-sib families derived from controlled intra-and interspecific crosses, and gene expression analyses could be used to associate sequence variation in the COL gene with phenotypic variation.
Genetic markers that exhibit extreme differentiation between species as signatures of divergent selection can be used as diagnostic markers to assess introgression of potentially adaptive alleles between interfertile species in parapatric and sympatric populations.FIR013 provides such a marker and our examination of the relative allelic frequencies at the FIR013 locus in parapatric Q. rubra/ Q. ellipsoidalis populations has shown that introgression of the 141 bp allele is consistently higher into the drought tolerant Q. ellipsoidalis than introgression of the 138 bp allele into the mesophilic species, Q. rubra (Table 2, Figure 4; [7]).This pattern of introgression of outlier alleles suggests different strength of selection against the 138 bp and 141 bp alleles in non-parental environments if interspecific gene flow is symmetric [41][42][43] as observed in the present study between Q. rubra and Q. ellipsoidalis.In contrast, in the only sympatric population, introgression of the 141 bp into Q.ellipsoidalis is much lower (5%) than introgression of the 138 bp into Q.rubra (19%), suggesting that competition between species might affect introgression at FIR013 outlier alleles.Also, the extent of interspecific gene flow between Q. rubra/Q.ellipsoidalis population pairs was not related to the introgression of alleles 138 bp and 141 bp, again suggesting that introgression of these alleles was also affected by environment-dependent selection.For example, the introgression of allele 138 bp into Q.rubra was 2% for N-QR and 18% for C-QR, while interspecific gene flow was very similar (Table 3).In future studies, a systematic assessment of introgression rates and soil and climate variables is needed in both sympatric and parapatric Q. rubra/Q.ellipsoidalis populations to analyze the effect of interspecific competition on introgression rates at outlier alleles and to test for an association between introgression rates at outlier alleles and environmental variables (e.g., climate, soil, habitat characteristics).
In a seedling common garden experiment, Q. ellipsoidalis had a later bud burst and leaf fall, a much slower growth rate and higher mortality than Q. rubra over two growing seasons [23].A slower growth strategy may be an adaptation to a xeric environment [44] and Q. ellipsoidalis is generally restricted to dry, sandy pine barrens [7,25] and is considered to be the most drought tolerant red oak species [16,45].In most mixed Q. ellipsoidalis/Q.velutina populations and the sympatric population (PV-DC), introgression of the 141 bp allele was higher into Q.velutina than into Q.ellipsoidalis.Quercus velutina is less water efficient than Q. ellipsoidalis and generally found in savannas [16,46,47], while Q. ellipsoidalis is generally found on very dry sandy sites, possibly creating differences in introgression patterns between the three species in their sympatric range.
Introgression of adaptive alleles results in adaptive trait transfer between species and may be crucial for the adaptation to rapidly changing environmental conditions (e.g., [2]).In addition to heterozygous genotypes (138 bp/141 bp) at outlier locus FIR013, we found rare homozygotes for the 141 bp (Q.rubra) allele in both Q. velutina and Q. ellipsoidalis as well as homozygotes for the 138 bp allele in Q. rubra (see Figure 3).Common garden studies with different drought treatments could reveal whether introgression of outlier alleles is associated with differences in drought tolerance.For example, assessments of phenotypic traits, such as growth and water use efficiency for different genotypes at FIR013 (138 bp/138 bp, 138 bp/141 bp, 141 bp/141 bp) in each species, could allow quantification of relative fitness related to variation in the poly-Q repeat of COL.
In conclusion, we have identified a genic microsatellite marker with repeat number variation resulting in two major alleles differentiating between the drought averse Q. rubra and the drought tolerant oak species Q. ellipsoidalis, Q. velutina and Q. coccinea.This candidate gene may be involved in adaptive differences between the species.Common garden experiments of seedlings with all possible genotype combinations at FIR013 in each species can elucidate the allelic and genotypic effects at FIR013 on survival and other fitness related traits and thus the relation between introgression of FIR013 alleles between species and adaptive trait introgression.

Supplementary Materials:
The following are available online at www.mdpi.com/1999-4907/8/1/3/s1, Figure S1: Relative genotypic frequencies in pure species of Q. rubra, Q. ellipsoidalis and Q. velutina in parapatric populations and one sympatric population of Q. rubra, Q. ellipsoidalis and Q. velutina (starred) at the FIR013 locus, Figure S2: Outlier loci analysis of 12 markers between Q. rubra and Q. ellipsoidalis as calculated in the program LOSITAN [7] with a 99% confidence interval; FIR013 was found to be an outlier even after the application of 10% false discovery rate (FDR), Figure S3: Outlier loci analysis of 12 markers between Q. rubra and Q. velutina as calculated in LOSITAN [7] with a 99% confidence interval; 3A05 was found not to be an outlier after the application of 10% FDR, Figure S4: Outlier loci analysis of 12 markers between Q. velutina and Q. ellipsoidalis as calculated in LOSITAN [7] with a 99% confidence interval; 3A05 was found not to be an outlier after the application of 10% FDR, Table S1: Location of collection sites for the 16 Quercus spp.populations, Table S2: Microsatellite marker characteristics, Table S3: Characteristics of samples used in the sequencing effort, Table S4: Relative frequency of a single nucleotide polymorphism in the coding region of a COL gene (Position 9, see Figure 2), Table S5: Genetic variation parameters by population and species across the 12 microsatellite markers, Table Genetic variation parameters by population and species across 11 non-outlier microsatellite markers, Table S7: Pairwise F ST value by species across all 12 microsatellite markers using GeneAlEx [6], Table S8: F ST -based outlier screening results across the 12 microsatellite markers using LOSITAN [7].

Figure 1 .
Figure 1.Map of sampled locations; areas 1-2 each include one population pair of Q. ellipsoidalis and Q. rubra, area 3 includes two population pairs of Q. ellipsoidalis and Q. rubra, area 4 includes two mixed Q. ellipsoidalis and Q. velutina populations, areas 5-8 each include one mixed Q. ellipsoidalis and Q. velutina population each, area 9 includes the single Q. coccinea population, and area 10 includes the

Figure 1 .
Figure 1.Map of sampled locations; areas 1-2 each include one population pair of Q. ellipsoidalis and Q. rubra, area 3 includes two population pairs of Q. ellipsoidalis and Q. rubra, area 4 includes two mixed Q. ellipsoidalis and Q. velutina populations, areas 5-8 each include one mixed Q. ellipsoidalis and Q. velutina population each, area 9 includes the single Q. coccinea population, and area 10 includes the single sympatric population of Q. rubra/Q.ellipsoidalis/Q.velutina; area number associated with individual populations is shown in TableS1; approximate species ranges are from Little (1971).

Figure 2 .
Figure 2. Sequence variation at the poly-Q repeat in the COL gene in four red oak species (S = G or C; Y = C or T); codons alternately highlighted in medium grey and white.

Figure 2 . 13 Figure 3 .
Figure 2. Sequence variation at the poly-Q repeat in the COL gene in four red oak species (S = G or C; Y = C or T); codons alternately highlighted in medium grey and white.Forests 2017, 8, 322 6 of 13

Figure 3 .
Figure 3. Relative genotypic frequencies in four oak species at the FIR013 locus.

Figure 4 .
Figure 4. FIR013 allele frequencies in each population separated by species.Orange represents the 138 bp allele and blue represents the 141 bp allele.Full names of populations listed in TableS1.

Figure 4 .
Figure 4. FIR013 allele frequencies in each population separated by species.Orange represents the 138 bp allele and blue represents the 141 bp allele.Full names of populations listed in TableS1.

Table 3 .
Relative frequency of hybrids and introgressive forms *.
Q. rubra with Introgression from Q. ellipsoidalis or Q. velutina Q. ellipsoidalis/Q.velutina with Introgression from Q. rubra F1 Hybrids of Q. rubra with Q. ellipsoidalis or Q.

Table 3 .
Relative frequency of hybrids and introgressive forms *.