Genetic, morphological, geographical and ecological approaches reveal phylogenetic relationships in complex groups, an example of recently diverged pinyon pine species (Subsection Cembroides)

Elucidating phylogenetic relationships and species boundaries within complex taxonomic groups is challenging for intrinsic and extrinsic (i.e., technical) reasons. Mexican pinyon pines are a complex group whose phylogenetic relationships and species boundaries have been widely studied but poorly resolved, partly due to intrinsic ecological and evolutionary features such as low morphological and genetic differentiation caused by recent divergence, hybridization and introgression. Extrinsic factors such as limited sampling and difficulty in selecting informative molecular markers have also impeded progress. Some of the Mexican pinyon pines are of conservation concern but others may remain unprotected because the species boundaries have not been established. In this study we combined approaches to resolve the phylogenetic relationships in this complex group and to establish species boundaries in four recently diverged taxa: P. discolor, P. johannis, P. culminicola and P. cembroides. We performed phylogenetic analyses using the chloroplast markers matK and psbA-trnH as well as complete and partial chloroplast genomes of species of Subsection Cembroides. Additionally, we performed a phylogeographic analysis combining genetic data (18 chloroplast markers), morphological data and geographical data to define species boundaries in four recently diverged taxa. Ecological divergence was supported by differences in climate among localities for distinct genetic lineages. Whereas the phylogenetic analysis inferred with matK and psbA-trnH was unable to resolve the relationships in this complex group, we obtained a resolved phylogeny with the use of the chloroplast genomes. The resolved phylogeny was concordant with a haplotype network obtained using chloroplast markers. In species with potential for recent divergence, hybridization or introgression, nonhierarchical network-based approaches are probably more appropriate to protect against misclassification due to incomplete lineage sorting. The boundaries among genetic lineages were delimited by the inclusion of morphological, geographical and ecological data in the haplotype network. These multiple lines of evidence help to assign species boundaries in this complex group. P. johannis, P. discolor, P. culminicola and P. cembroides are different species based on their genetic, morphological and ecological niche differences. We suggest a reevaluation of the conservation status of these species considering the information generated in this study.

The Pinus taeda genome is characterized by diverse and highly diverged repetitive sequences

BACKGROUND

In today's age of genomic discovery, no attempt has been made to comprehensively sequence a gymnosperm genome. The largest genus in the coniferous family Pinaceae is Pinus, whose 110-120 species have extremely large genomes (c. 20-40 Gb, 2N = 24). The size and complexity of these genomes have prompted much speculation as to the feasibility of completing a conifer genome sequence. Conifer genomes are reputed to be highly repetitive, but there is little information available on the nature and identity of repetitive units in gymnosperms. The pines have extensive genetic resources, with approximately 329000 ESTs from eleven species and genetic maps in eight species, including a dense genetic map of the twelve linkage groups in Pinus taeda.

RESULTS

We present here the Sanger sequence and annotation of ten P. taeda BAC clones and Genome Analyzer II whole genome shotgun (WGS) sequences representing 7.5% of the genome. Computational annotation of ten BACs predicts three putative protein-coding genes and at least fifteen likely pseudogenes in nearly one megabase of sequence. We found three conifer-specific LTR retroelements in the BACs, and tentatively identified at least 15 others based on evidence from the distantly related angiosperms. Alignment of WGS sequences to the BACs indicates that 80% of BAC sequences have similar copies (> or = 75% nucleotide identity) elsewhere in the genome, but only 23% have identical copies (99% identity). The three most common repetitive elements in the genome were identified and, when combined, represent less than 5% of the genome.

CONCLUSIONS

This study indicates that the majority of repeats in the P. taeda genome are 'novel' and will therefore require additional BAC or genomic sequencing for accurate characterization. The pine genome contains a very large number of diverged and probably defunct repetitive elements. This study also provides new evidence that sequencing a pine genome using a WGS approach is a feasible goal.

Comparative mapping among subsectionAustrales(genusPinus, family Pinaceae)

Comparative mapping in conifers has not yet been used to test for small-scale genomic disruptions such as inversions, duplications, and deletions occurring between closely related taxa. Using comparative mapping to probe this smaller scale of inquiry may provide clues about speciation in a phylogenetically problematic taxon, the diploxylon pine subsection Australes (genus Pinus , family Pinaceae). Genetic maps were constructed for two allopatric species of Australes, P. elliottii var. elliottii and P. caribaea var. hondurensis , using microsatellites and an F1hybrid. A third map was generated directly from the meiotic products of an adult F1hybrid, eliminating the need for an F2generation. Numerous small-scale disruptions were detected in addition to synteny and collinearity, and these included (1) map shrinkage, (2) a paracentric inversion, (3) transmission ratio distortion, and (4) mild selection against a parental haplotype. Such cryptic signatures of genomic divergence between closely related interfertile species are useful in elucidating this problematic evolutionary history.

Cross-species transferability and mapping of genomic and cDNA SSRs in pines

Two unigene datasets of Pinus taeda and Pinus pinaster were screened to detect di-, tri- and tetranucleotide repeated motifs using the SSRIT script. A total of 419 simple sequence repeats (SSRs) were identified, from which only 12.8% overlapped between the two sets. The position of the SSRs within their coding sequences were predicted using FrameD. Trinucleotides appeared to be the most abundant repeated motif (63 and 51% in P. taeda and P. pinaster, respectively) and tended to be found within translated regions (76% in both species), whereas dinucleotide repeats were preferentially found within the 5'- and 3'-untranslated regions (75 and 65%, respectively). Fifty-three primer pairs amplifying a single PCR fragment in the source species (mainly P. taeda), were tested for amplification in six other pine species. The amplification rate with other pine species was high and corresponded with the phylogenetic distance between species, varying from 64.6% in P. canariensis to 94.2% in P. radiata. Genomic SSRs were found to be less transferable; 58 of the 107 primer pairs (i.e. 54%) derived from P. radiata amplified a single fragment in P. pinaster. Nine cDNA-SSRs were located to their chromosomes in two P. pinaster linkage maps. The level of polymorphism of these cDNA-SSRs was compared to that of previously and newly developed genomic-SSRs. Overall, genomic SSRs tend to perform better in terms of heterozygosity and number of alleles. This study suggests that useful SSR markers can be developed from pine ESTs.

Genetic maps for Pinus elliottii var. elliottii and P. caribaea var. hondurensis using AFLP and microsatellite markers

Genetic maps for individual Pinus elliottii var. elliottii and P. caribaea var. hondurensis trees were generated using a pseudo-testcross mapping strategy. A total of 329 amplified fragment length polymorphic (AFLP) and 12 microsatellite markers were found to segregate in a sample of 93 interspecfic F(1) progeny. The male P. caribaea var. hondurensis parent was more heterozygous than the female P. elliottii var. elliottii parent with 19% more markers segregating on the male side. Framework maps were constructed using a LOD 5 threshold for grouping and interval support threshold of LOD 2. The framework map length for the P. elliottii var. elliottii megagametophyte parent (1,170 cM Kosambi; 23 linkage groups) was notably smaller than the P. caribaea var. hondurensis pollen parent (1,658 cM Kosambi; 27 linkage groups). The difference in map lengths was assumed to be due to sex-related recombination variation, which has been previously reported for pines, as the difference in map lengths not be accounted for by the larger number of markers mapping to the P. caribaea var. hondurensis parent - 109 compared with 78 in P. elliottii var. elliottii parent. Based on estimated genome sizes for these species, the framework maps for P. elliottii var. elliottii and P. caribaea var. hondurensis covered 82% and 88% of their respective genomes. The pseudo-testcross strategy was extended to include AFLP and microsatellite markers in an intercross configuration. These comprehensive maps provided further genome coverage, 1,548 and 1,828 cM Kosambi for P. elliottii var. elliottii and P. caribaea var. hondurensis, respectively, and enabled homologous linkage groups to be identified in the two parental maps. Homologous linkage groups were identified for 11 out of 24 P. elliottii var. elliottii and 10 out of 25 P. caribaea var. hondurensis groups. A higher than expected level of segregation distortion was found for both AFLP and microsatellite markers. An explanation for this segregation distortion was not clear, but it may be at least in part due to genetic mechanisms for species isolation in this wide cross.