Demographic modelling helps track the rapid and recent divergence of a conifer species pair from Central Mexico.

Secondary contact of recently diverged species may have several outcomes, ranging from rampant hybridization to reinforced reproductive isolation. In plants, selfing tolerance and disjunct reproductive phenology may lead to reproductive isolation at contact zones. However, they may also evolve under both allopatric or parapatric frameworks and originate from adaptive and/or neutral forces. Inferring the historical demography of diverging taxa is thus a crucial step to identify factors that may have led to putative reproductive isolation. We explored various competing demographypotheses to account for the rapid divergence of a fir species complex (Abies flinckii-A. religiosa) distributed in "sky-islands" across central Mexico (i.e., along the Trans-Mexican Volcanic Belt; TMVB). Despite co-occurring in two independent sympatric regions (west and centre), these taxa rarely interbreed because of disjunct reproductive phenologies. We genotyped 1147 single nucleotide polymorphisms, generated by GBS (genotyping by sequencing), across 23 populations, and compared multiple scenarios based on the geological history of the TMVB. The best-fitting model revealed one of the most rapid and complete speciation cases for a conifer species-pair, dating back to ~1.2 million years ago. Coupled with the lack of support for stepwise colonization, our coalescent inferences point to an early cessation of interspecific gene flow under parapatric speciation; ancestral gene flow during divergence was asymmetrical (mostly from western firs into A. religiosa) and exclusive to the most ancient (i.e., central) contact zone. Factors promoting rapid reproductive isolation should be explored in other slowly evolving species complexes as they may account for the large tropical and subtropical diversity.

Non-adaptive evolutionary processes governed the diversification of a temperate conifer lineage after its migration into the tropics.

Constructing phylogenetic relationships among closely related species is a recurrent challenge in evolutionary biology, particularly for long-lived taxa with large effective population sizes and uncomplete reproductive isolation, like conifers. Conifers further have slow evolutionary rates, which raises the question of whether adaptive or non/adaptive processes were predominantly involved when they rapidly diversified after migrating from temperate regions into the tropical mountains. Indeed, fine-scale phylogenetic relationships within several conifer genus remain under debate. Here, we studied the phylogenetic relationships of endemic firs (Abies, Pinaceae) discontinuously distributed in the montane forests from the Southwestern United States to Guatemala, and addressed several hypotheses related to adaptive and non-adaptive radiations. We derived over 80 K SNPs from genotyping by sequencing (GBS) for 45 individuals of nine Mesoamerican species to perform phylogenetic analyses. Both Maximum Likelihood and quartets-inference phylogenies resulted in a well-resolved topology, showing a single fir lineage divided in four subgroups that coincided with the main mountain ranges of Mesoamerica; thus having important taxonomic implications. Such subdivision fitted a North-South isolation by distance framework, in which non-adaptive allopatric processes seemed the rule. Interestingly, several reticulations were observed within subgroups, especially in the central-south region, which may explain past difficulties for generating infrageneric phylogenies. Further evidence for non-adaptive processes was obtained from analyses of 21 candidate-gene regions, which exhibited diminishing values of πa/πs and Ka/Ks with latitude, thus indicating reduced efficiency of purifying selection towards the Equator. Our study indicates that non-adaptive allopatric processes may be key generators of species diversity and endemism in the tropics.

Isolation and characterization of novel microsatellites for Abies koreana and A. nephrolepis (Pinaceae).

Abies koreana is an endemic and rare species from Korea and is classified as endangered by the International Union for Conservation of Nature. Although the genetic diversity assessment for current population of A. koreana needs to be performed urgently, no microsatellite markers have been developed for this species. In the present study, we developed 22 novel polymorphic microsatellite loci and the characteristics of these loci were determined in A. koreana as well as in Abies nephrolepis, the most closely related species, and these loci were compared with previously reported microsatellite markers developed for the Abies genus. Genomic sequence (161 Mbp; 325,776 reads) was obtained from one individual of A. koreana using Roche 454 GS-FLX Titanium sequencing and 19,258 repeat motifs were identified from it. A total of 288 primer pairs with high copy numbers of di-repeat motifs were evaluated for amplification in A. koreana and A. nephrolepis. A total of 71 primer pairs successfully amplified fragments, of which 22 showed polymorphisms in A. koreana and A. nephrolepis. The average expected diversity was 0.767 and 0.717 in A. koreana and A. nephrolepis, respectively; these heterozygosity levels were moderate compared to the previously reported microsatellite loci from Abies species. This is the first set of microsatellite markers developed for A. koreana as well as A. nephrolepis and further population genetic studies of both species and genetic delimitation can be carried out for the species conservation and management.

Less pollen-mediated gene flow for more signatures of glacial lineages: congruent evidence from balsam fir cpDNA and mtDNA for multiple refugia in eastern and central North America.

The phylogeographic structure and postglacial history of balsam fir (Abies balsamea), a transcontinental North American boreal conifer, was inferred using mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) markers. Genetic structure among 107 populations (mtDNA data) and 75 populations (cpDNA data) was analyzed using Bayesian and genetic distance approaches. Population differentiation was high for mtDNA (dispersed by seeds only), but also for cpDNA (dispersed by seeds and pollen), indicating that pollen gene flow is more restricted in balsam fir than in other boreal conifers. Low cpDNA gene flow in balsam fir may relate to low pollen production due to the inherent biology of the species and populations being decimated by recurrent spruce budworm epidemics, and/or to low dispersal of pollen grains due to their peculiar structural properties. Accordingly, a phylogeographic structure was detected using both mtDNA and cpDNA markers and population structure analyses supported the existence of at least five genetically distinct glacial lineages in central and eastern North America. Four of these would originate from glacial refugia located south of the Laurentide ice sheet, while the last one would have persisted in the northern Labrador region. As expected due to reduced pollen-mediated gene flow, congruence between the geographic distribution of mtDNA and cpDNA lineages was higher than in other North American conifers. However, concordance was not complete, reflecting that restricted but nonetheless detectable cpDNA gene flow among glacial lineages occurred during the Holocene. As a result, new cpDNA and mtDNA genome combinations indicative of cytoplasmic genome capture were observed.

Evolution of an ancient microsatellite hotspot in the conifer mitochondrial genome and comparison with other plants.

In plants, mitochondrial sequence tandem repeats (STRs) have been associated with intragenomic recombination, a process held responsible for evolutionary outcomes such as gene regulation or cytoplasmic male-sterility. However, no link has been established between the recurrent accumulation of STRs and increased mutation rates in specific regions of the plant mtDNA genome. Herein, we surveyed this possibility by comparing, in a phylogenetic context, the variation of a STR-rich mitochondrial intron (nad5-4) with eleven mtDNA genes devoid of STRs within Abies (Pinaceae) and its related genera. This intron has been accumulating repeated stretches, generated by at least three-independent insertions, before the split of the two Pinaceae subfamilies, Abietoideae and Pinoideae. The last of these insertions occurred before the divergence of Abies and produced, exclusively within this genus, a tenfold increase of both the indel and substitution rates in the STR hotspot of the intron. The regions flanking the STRs harbored mutation rates as low as those estimated in mitochondrial genes devoid of repeated stretches. Further searches in complete plant mtDNA genomes, and previous studies reporting polymorphic mtSTRs, revealed that repeated stretches are common in all sorts of plants, but their accumulation in STR hotspots appears to be taxa specific. Our study suggests a new mutagenic role for repeated sequences in the plant mtDNA.

Phylogeny, diversification rates and species boundaries of Mesoamerican firs (Abies, Pinaceae) in a genus-wide context.

The genus Abies is distributed discontinuously in the temperate and subtropical montane forests of the northern hemisphere. In Mesoamerica (Mexico and northern Central America), modern firs originated from the divergence of isolated mountain populations of migrating North American taxa. However, the number of ancestral species, migratory waves and diversification speed of these taxa is unknown. Here, variation in repetitive (Pt30204, Pt63718, and Pt71936) and non-repetitive (rbcL, rps18-rpl20 and trnL-trnF) regions of the chloroplast genome was used to reconstruct the phylogenetic relationships of the Mesoamerican Abies in a genus-wide context. These phylogenies and two fossil-calibrated scenarios were further employed to estimate divergence dates and diversification rates within the genus, and to test the hypothesis that, as in many angiosperms, conifers may exhibit accelerated speciation rates in the subtropics. All phylogenies showed five main clusters that mostly agreed with the currently recognized sections of Abies and with the geographic distribution of species. The Mesoamerican taxa formed a single group with species from southwestern North America of sections Oiamel and Grandis. However, populations of the same species were not monophyletic within this group. Divergence of this whole group dated back to the late Paleocene and the early Miocene depending on the calibration used, which translated in very low diversification rates (r(0.0)=0.026-0.054, r(0.9)=0.009-0.019 sp/Ma). Such low rates were a constant along the entire genus, including both the subtropical and temperate taxa. An extended phylogeographic analysis on the Mesoamerican clade indicated that Abies flinckii and A. concolor were the most divergent taxa, while the remaining species (A. durangensis, A. guatemalensis, A. hickelii, A. religiosa and A. vejari) formed a single group. Altogether, these results show that divergence of Mesoamerican firs coincides with a model of environmental stasis and decreased extinction rate, being probably prompted by a series of range expansions and isolation-by-distance.

Ancestry and divergence of subtropical montane forest isolates: molecular biogeography of the genus Abies (Pinaceae) in southern México and Guatemala.

The genus Abies has a complex history in southern México and Guatemala. In this region, four closely related species, Abies flinckii, A. guatemalensis, A. hickelii, and A. religiosa, are distributed in fragmented and isolated montane populations. Range-wide genetic variation was investigated across species using cytoplasmic DNA markers with contrasted inheritance. Variation at two maternally inherited mitochondrial DNA markers was low. All species shared two of the nine mitotypes detected, while the remaining seven mitochondrial DNA types were restricted to a few isolated stands. Mitochondrial genetic differentiation across taxa was high (G(ST) = 0.933), it was not related to the taxonomic identity (amova; P > 0.05) of the populations, and it was not phylogeographically structured (G(ST) approximately N(ST)). In contrast, variation at three paternally inherited chloroplast DNA microsatellites was high. Chloroplast genetic differentiation was lower (G(ST) = 0.402; R(ST) = 0.547) than for mitochondrial DNA, but it was significantly related to taxonomy (amova; P < 0.001), and exhibited a significant phylogeographical structure (G(ST) < R(ST)). Different analyses of population structure indicated that A. flinckii was the most divergent taxon, while the remaining three species formed a relatively homogeneous group. However, a small number of the populations of these three taxa, all located at the limits of their respective ranges or in the Transverse Volcanic Belt, diverged from this main cluster. These trends suggest that the Mesoamerican Abies share a recent common ancestor and that their divergence and speciation is mainly driven by genetic drift and isolation during the warm interglacial periods.