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.

Phylogeny and evolutionary history of Pinaceae updated by transcriptomic analysis.

Pinaceae comprises 11 genera, and represents the largest family of conifers with an extensive wild distribution in the Northern Hemisphere. Intergeneric relationships of Pinaceae have been investigated using many morphological characters and molecular markers, but phylogenetic positions of four genera, including Cathaya, Cedrus, Nothotsuga and Pseudolarix, remain controversial or have not been completely resolved. To completely resolve the intergeneric relationships of Pinaceae, we conducted a comparative transcriptomic study of 14 species representing all Pinaceae genera. Multiple data sets, containing up to 6,369,681 sites across 4676 loci, were analyzed using concatenation and coalescent methods. Our study generated a robust topology, which divides Pinaceae into two clades, one (pinoid) including Cathaya, Larix, Picea, Pinus, and Pseudotsuga, and the other (abietoid) including Abies, Cedrus, Keteleeria, Nothotsuga, Pseudolarix, and Tsuga. Cathaya and Pinus form a clade sister to Picea; Cedrus is sister to the remaining abietoid genera, and the two genera Nothotsuga and Tsuga form a clade sister to Pseudolarix. The discordant positions of Cathaya, Cedrus and Pseudolarix in different gene trees could be explained by ancient radiation and/or molecular homoplastic evolution. The hybrid origin hypothesis of Nothotsuga is not supported. Based on molecular dating, extant Pinaceae genera diverged since about 206 Mya, earlier than the break-up of Pangea, and the divergence among the pinoid genera occurred earlier than the split among the abietoid genera. Moreover, our study indicates that two radiation events occurred in the evolution of Pinaceae genera, and some important morphological characters evolved multiple times based on ancestral state reconstruction.

From America to Eurasia: a multigenomes history of the genus Abies.

The origin of conifer genera, the main components of mountain temperate and boreal forests, was deemed to arise in the Mesozoic, although paleontological records and molecular data point to a recent diversification, presumably related to Neogene cooling. The geographical area(s) where the modern lines of conifers emerged remains uncertain, as is the sequence of events leading to their present distribution. To gain further insights into the biogeography of firs (Abies), we conducted phylogenetic analyses of chloroplast, mitochondrial and nuclear markers. The species tree, generated from ten single-copy nuclear genes, yielded probably the best phylogenetic hypothesis available for Abies. The tree obtained from five regions of chloroplast DNA largely corresponded to the nuclear species tree. Ancestral area reconstructions based on fossil calibrated chloroplast DNA and nuclear DNA trees pointed to repeated intercontinental migrations. The mitochondrial DNA haplotype tree, however, disagreed with nuclear and chloroplast DNA trees. It consisted of two clusters: one included mainly American haplotypes, while the other was composed of only Eurasian haplotypes. Presumably, this conflict is due to inter-continental migrations and introgressive hybridization, accompanied by the capture of the mitotypes from aboriginal species by the invading firs. Given that several species inhabiting Northeastern Asia carry American mitotypes and mutations typical for the American cluster, whereas no Asian mitotypes were detected within the American species, we hypothesize that Abies migrated from America to Eurasia, but not in the opposite direction. The direction and age of intercontinental migrations in firs are congruent with other conifers, such as spruces and pines of subsection Strobus, suggesting that these events had the same cause.

Identifying Genetic Signatures of Natural Selection Using Pooled Population Sequencing in Picea abies.

The joint inference of selection and past demography remain a costly and demanding task. We used next generation sequencing of two pools of 48 Norway spruce mother trees, one corresponding to the Fennoscandian domain, and the other to the Alpine domain, to assess nucleotide polymorphism at 88 nuclear genes. These genes are candidate genes for phenological traits, and most belong to the photoperiod pathway. Estimates of population genetic summary statistics from the pooled data are similar to previous estimates, suggesting that pooled sequencing is reliable. The nonsynonymous SNPs tended to have both lower frequency differences and lower FST values between the two domains than silent ones. These results suggest the presence of purifying selection. The divergence between the two domains based on synonymous changes was around 5 million yr, a time similar to a recent phylogenetic estimate of 6 million yr, but much larger than earlier estimates based on isozymes. Two approaches, one of them novel and that considers both FST and difference in allele frequencies between the two domains, were used to identify SNPs potentially under diversifying selection. SNPs from around 20 genes were detected, including genes previously identified as main target for selection, such as PaPRR3 and PaGI.

[Phylogeny of firs (genus Abies, Pinaceae) based on multilocus nuclear markers (AFLP)].

To study the phylogenetic relationships, evolutionary history, and molecular systematics of firs (genus Abies), the phylogenetic reconstruction, based on nuclear multilocus markers­amplified fragment length polymorphism (AFLP)­was conducted. Using seven combinations of selective primers, 84 samples of 39 taxa were genotyped for 553 polymorphic AFLP loci. A comparison with our earlier chloroplast and mitochondrial phylogenies of the genus (in 2014) shows that the nuclear phylogeny generally is more congruent to the chloroplast tree. Most of the clades resolved by the chloroplast phylogeny were supported also in the AFLP tree. Employing the nuclear DNA-based tree, we revealed the presence of new groups and the differences in the topology of several clades. AFLP confirmed the monophyly of Asian species of section Balsamea and their sister position in relation to the American group of species of this section. As shown by the tree of chloroplast DNA, Asian species of section Balsamea do not form a monophyletic group, but belong to the clade comprising the majority of Asian species. Phylogenetically mitochondrial DNA data to a large extent are not congruent to the nuclear and chloroplast DNA trees, and are more in line with geographical distribution of species. Conflicts between nuclear and cytoplasmic phylogeny were analyzed. Taking them into account, we consider the hypothesis of a hybrid origin of particular groups of firs, including ancient hybridization in section Balsamea. A comparison of molecular data with traditional taxonomy of the genus is discussed.

The complete plastid genome sequence of Abies koreana (Pinaceae: Abietoideae).

The nucleotide sequence of the chloroplast genome from Abies koreana is the first to have complete genome sequence from genus Abies of family Pinaceae. The circular double-stranded DNA, which consists of 121,373 base pairs (bp), contains a pair of very short inverted repeat regions (IRa and IRb) of 264 bp each, which are separated by a small and large single-copy regions (SSC and LSC) of 54,197 and 66,648 bp, respectively. The genome contents of 114 genes (68 peptide-encoding genes, 35 tRNA genes, four rRNA genes, six open reading frames and one pseudogene) are similar to the chloroplast DNA of other species of Abietoideae. Loss of ndh genes was also identified in the genome of A. koreana like other genomes in the family Pinaceae. Thirteen genes contain one (11 genes) or two (rps12 and ycf3 genes) introns. In phylogenetic analysis, the tree confirms that Abies, Keteleeria and Cedrus are strongly supported as monophyletic. Other inverted repeat sequences located in 42-kb inversion points (1186 bp) include trnS-psaM-ycf12- ψtrnG genes.

Integrating phylogeography and paleoecology to investigate the origin and dynamics of hybrid zones: insights from two widespread North American firs.

Secondary contact between closely related taxa routinely occurs during postglacial migrations. After initial contact, the location of hybrid zones may shift geographically or remain spatially stable over time in response to various selective pressures or neutral processes. Studying the extent and direction of introgression using markers having contrasted levels of gene flow can help unravel the historical dynamics of hybrid zones. Thanks to their contrasted maternal and paternal inheritance, resulting in different levels of gene flow for mitochondrial and chloroplast DNA (mtDNA and cpDNA), the Pinaceae stand out as a relevant biological model for this purpose. The objective of the study was to assess whether the hybrid zone between Abies balsamea and Abies lasiocarpa (two largely distributed Pinaceae) has moved or remained stable over time by analysing the distribution of cytoplasmic DNA variation as well as published palaeobotanical data. Interspecific gene flow was higher for cpDNA than mtDNA markers; hence, the geographic distribution of mitotypes was more congruent with species distributions than chlorotypes. This genetic signature was contrary to expectations under a moving hybrid zone scenario, as well as empirical observations in other conifers. Genetic evidence for this rare instance of stable hybrid zone was corroborated by the colonization chronology derived from published fossil data, indicating that the two fir species initially came into contact in the area corresponding to the current sympatric zone 11 kyr ago. While an explanatory analysis suggested the putative influence of various environmental factors on the relative abundance of cytoplasmic genome combinations, further research appears necessary to assess the role of both demographic history and selective factors in driving the dynamics of hybrid zones.

Phylogenetic relationships, possible ancient hybridization, and biogeographic history of Abies (Pinaceae) based on data from nuclear, plastid, and mitochondrial genomes.

Abies, the second largest genus of Pinaceae, consists of approximately 48 species occurring in the north temperate region. Previous molecular phylogenetic studies improved our understanding of relationships within the genus, but were limited by relying on only DNA sequence data from single genome and low taxonomic sampling. Here we use DNA data from three genomes (sequences of internal transcribed spacer of nrITS, three chloroplast DNA intergenic spacers, and two mitochondrial intergenic spacers) from 42 species to elucidate species relationships and construct the biogeographic history of Abies. We further estimated the divergence times of intercontinental disjunction using a relaxed molecular clock calibrated with three macro-fossils. Our phylogenetic analyses recovered six robust clades largely consistent with previous classifications of sections. A sister relationship between the eastern Asian and Europe-Mediterranean clades was highly supported. The monophyly of section Balsamea, disjunct in Far East and western North America, is supported by the nrITS data but not by the cpDNA data. Discordance on placement of section Balsamea between the paternally inherited cpDNA and maternally inherited mtDNA trees was also observed. The data suggested that ancient hybridization was likely involved in the origin of sect. Balsamea. Results from biogeographic analyses and divergence time estimation suggested an origin and early diversification of Abies in an area of high latitude around the Pacific during the Eocene. The present disjunction in eastern Asia and Europe-Mediterranean area of Abies was likely the result of southward migration and isolation by the Turgai Strait in the Late Eocene. An 'out-of-America' migration, for the origin of an eastern Asian and western North American disjunct species pairs in section Amabilis was supported. The results suggested a western North American origin of the section with subsequent dispersal across the Bering Land Bridge (BLB) to Japan during the Middle Miocene.

Phylogeography of SW Mediterranean firs: different European origins for the North African Abies species.

The current distribution of Western Mediterranean Abies species is a result of complex geodynamic processes and climatic oscillations that occurred in the past. Abies sect. Piceaster offers a good study model to explore how geo-climatic oscillations might have influenced its expansion and diversification on both sides of the W Mediterranean basin. We investigated the genetic variation within and among nine populations from five Abies species by molecular markers with high and low mutation rates and contrasting inheritance (AFLP and cpSSR). Analyses revealed the opening of the Strait of Gibraltar as an effective barrier against gene flow between the Southern Iberian (A. pinsapo) and North African (A. marocana and A. tazaotana) firs. The A. pinsapo populations in Spain and likewise those of the A. marocana - A. tazaotana population complex were not differentiated, and no evidence was found to distinguish A. tazaotana at the species level. Diversification of Abies across North Africa could occur by way of at least two vicariant events from Europe, in the west, giving rise to the A. marocana - A. tazaotana complex, and in the east, giving A. numidica. Secondary contacts among species from Abies sect. Piceaster (A. pinsapo and A. numidica), and with A. alba (Abies sect. Abies) are also indicated. However, there is a closer relationship between the Algerian fir (A. numidica) and the North Mediterranean widespread A. alba, than with the Moroccan firs (A. marocana and A. tazaotana) or the Southern Iberian (A. pinsapo). We also discuss the distribution range of these taxa in its paleogeological and paleoclimatic context, and propose that part of the modern geography of the South-Western Mediterranean firs might be traced back to the Tertiary.

Climate isn't everything: competitive interactions and variation by life stage will also affect range shifts in a warming world.

PREMISE OF THE STUDY: The extent to which climate controls species' range limits is a classic biological question that is particularly relevant given anthropogenic climate change. While climate is known to play a role in species distributions, biotic interactions such as competition also affect range limits. Furthermore, climatic and biotic controls of ranges may vary in strength across life stages, implying complex range shift dynamics with climate change. METHODS: We quantified climatic and competitive influences on growth of juvenile and adult trees of three conifer species on Mt. Rainier, Washington, United States. We collected annual growth data of these trees, which we compared to the competitive environment and annual climate (100 years of data) experienced by each individual. KEY RESULTS: We found that the relationships between growth and climate and between growth and competition differed by life stage and location. Growth was sensitive to heavy snowpack and cold temperatures at high elevation upper limits (treeline), but growth was poorly explained by climate in low elevation closed-canopy forests. Competitive effects on growth were more important for saplings than adults, but did not become more important at either upper or lower range limits. CONCLUSIONS: In all, our results suggest that range shifts under climate change will differ at leading vs. trailing edges. At treeline, warmer temperatures will lead to increased growth and likely to range expansion. However, climate change will have less dramatic effects in low elevation closed-canopy forest communities, where growth is less strongly limited by climate, especially at young life stages.

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.

From pattern to process: species and functional diversity in fungal endophytes of Abies beshanzuensis.

The biodiversity-functional relationship in fungal ecology was recently developed and debated, but has rarely been addressed in endophytes. In this study, an integrative culture system was designed to capture a rich fungal consortium from the conifer Abies beshanzuensis. Results indicate an impressive diversity of fungal lineages (a total of 84 taxa classified in Dikarya) and a relatively high proportion of hitherto unknown species (27.4%). The laccase gene was used as a functional marker due to its involvement in lignocellulose degradation. Remarkable diversity of laccase genes was found across a wide range of taxa, with at least 35 and 19 distinct sequences in ascomycetes and basidiomycetes respectively, were revealed. Many groups displayed variable ability to decompose needles. Furthermore, many ascomycetes, including three volatile-producing Muscodor species (Xylariaceae), showed the ability to inhibit pathogens. Notably, most laccase-producing species showed little or no antibiosis and vice versa. Clavicipitalean and ustilaginomycetous fungi, specifically toxic to insects, were inferred from taxonomic information. Intra-specific physiological variation in Pezicula sporulosa, a second dominant species, was clearly high. We conclude that a suite of defensive characteristics in endophytes contributes to improving host fitness under various stresses and that a diversity of laccase genes confers an ecological advantage in competition for nutrients. Intra-specific diversity may be of great ecological significance for ecotypic adaptation. These findings suggest a fair degree of functional complementarity rather than redundancy among endemic symbionts of natural plant populations.

Species delimitation and biogeography of two fir species (Abies) in central China: cytoplasmic DNA variation.

It remains unclear how speciation history might contribute to species-specific variation and affect species delimitation. We examined concordance between cytoplasmic genetic variation and morphological taxonomy in two fir species, Abies chensiensis and A. fargesii, with overlapping distributions in central China. Range-wide genetic variation was investigated using mitochondrial (mt) and plastid (pt) DNA sequences, which contrast in their rates of gene flow. Four mtDNA haplotypes were recovered and showed no obvious species' bias in terms of relative frequency. In contrast, a high level of ptDNA variation was recorded in both species with 3 common ptDNA haplotypes shared between them and 21 rare ptDNA haplotypes specific to one or other species. We argue that the lack of concordance between morphological and molecular variation between the two fir species most likely reflects extensive ancestral polymorphism sharing for both forms of cytoplasmic DNA variation. It is feasible that a relatively fast mutation rate for ptDNA contributed to the production of many species-specific ptDNA haplotypes, which remained rare due to insufficient time passing for their spread and fixation in either species, despite high levels of intraspecific ptDNA gene flow. Our phylogeographic analyses further suggest that polymorphisms in both organelle genomes most likely originated during and following glacial intervals preceding the last glacial maximum, when species distributions became fragmented into several refugia and then expanded in range across central China.

[Spatial patterns of dominant tree species in sub-alpine Betula-Abies forest in West Sichuan of China].

Based on the investigation in a 4 hm2 Betula-Abies forest plot in sub-alpine area in West Sichuan of China, and by using point pattern analysis method in terms of O-ring statistics, the spatial patterns of dominant species Betula albo-sinensis and Abies faxoniana in different age classes in study area were analyzed, and the intra- and inter-species associations between these age classes were studied. B. albo-sinensis had a unimodal distribution of its DBH frequency, indicating a declining population, while A. faxoniana had a reverse J-shaped pattern, showing an increasing population. All the big trees of B. albo-sinensis and A. faxoniana were spatially in random at all scales, while the medium age and small trees were spatially clumped at small scales and tended to be randomly or evenly distributed with increasing spatial scale. The maximum aggregation degree decreased with increasing age class. Spatial association mainly occurred at small scales. A. faxoniana generally showed positive intra-specific association, while B. albo-sinensis generally showed negative intra-specific association. For the two populations, big and small trees had no significant spatial association, but middle age trees had negative spatial association. Negative inter-specific associations of the two populations were commonly found in different age classes. The larger the difference of age class, the stronger the negative inter-specific association.

Cytogenetic and molecular characterization of the Abies alba genome and its relationship with other members of the Pinaceae.

Genome size, karyotype structure, heterochromatin distribution, position and number of ribosomal genes, as well as the ITS2 sequence of the internal transcribed spacer (ITS) were analysed in silver fir (Abies alba Mill.). The analysis also included characterization of the Arabidopsis-type of telomeric repeats in silver fir and in related species. The results were compared with results from other species of the Pinaceae, to evaluate phylogeny and chromosomal and molecular evolution in the Pinaceae. Integrated chromosomal data provided insights into chromosome and karyotype evolution in the Pinaceae. The evolutionary trend for GC-rich heterochromatic blocks seems to involve loss of blocks that are not associated with rDNA. Similarly, numerous large blocks of interstitial plant telomeric repeats that are typical for all analysed species of the genus Pinus were not observed in the evolutionarily younger genera, such as Abies, Picea and Larix. On the contrary, the majority of telomeric sequences in these three genera appeared confined to the chromosome ends. We confirmed the current position of Abies and Tsuga in subfamily Abietoideae and the position of Pinus in the subfamily Pinoideae based on ITS2 sequences. Pseudotsuga is placed together with Larix into the subfamily Laricoideae. We conclude that the current position of the genus Picea in the subfamily Abietoideae should be reconsidered and, possibly, the genus Picea should be reclassified as a separate subfamily, Piceoideae, as recently proposed.

Phytochemical and biological studies of Abies species.

In the present review, the literature data on the phytochemical and biological investigations on the genus of Abies are summarized with 110 references. Up to now, 277 compounds were isolated from 19 plants of Abies species. The chemical constituents are mostly terpenoids, flavonoids, and lignans, together with minor constituents of phenols, steroids, and others. The crude extracts and metabolites have been found to possess various bioactivities including insect juvenile hormone, antitumor, antimicrobial, anti-ulcerogenic, anti-inflammatory, antihypertensive, antitussive, and CNS (central nervous system) activities.

[Population dynamics of endangered plant species Abies chensiensis].

In order to know the endangered status and causes of Abies chensiensis in Qinlin Mountains, a field investigation on 18 plots was conducted on its age structure, life table and fecundity, and its population dynamics were predicted by time sequence model. The analysis on the age structure of Abies chensiensis populations showed that there were fewer young individuals, but middle-aged and old individuals were relatively rich. The population D in Abies chensiensis-Indigofera amblyantha-Carex lanceolata association showed a relatively stable development tendency, while other four populations (A, B, C and E) in Abies chensiensis-Pinus tabulaeformis-Sinarundinaria nitida-Carex lanceolata association, Abies chensiensis-Quercus aliena var. acutserrata-Litsea pungens-Carex lanceolata association, Abies chensiensis-Betula albo-sinensis-Sinarundinaria nitida-Duchesnea indica association, and Abies chensiensis-Pinus tabulaeformis-Smilax stans-Carex lanceolata association all showed an obviously declining tendency. The analysis on the life tables and survival curves showed that the survival curve of Abies chensiensis populations belonged to Deevey III, and the death peak of different populations was in the period of 60--100 years old. The number difference among populations reflected the population habitat. Time sequence prediction indicated the numbers of old individuals would be increased at the beginning, and decreased finally in 20, 40, and 80 years. It was difficult to maintain the population stability. Analysis on 10 ecological factors showed that tree coverage, soil organisms and air humidity influenced population positively, and human disturbance and sunlight influenced population negatively. In situ conservation should be taken as the most important management countermeasure, and natural regeneration should be promoted. At the same time, artificial population should be expanded.

Taxonomy of gymnospermae: multivariate analyses of leaf fatty acid composition.

The fatty acid composition of photosynthetic tissues from 137 species of gymnosperms belonging to 14 families was determined by gas chromatography. Statistical analysis clearly discriminated four groups. Ginkgoaceae, Cycadaceae, Stangeriaceae, Zamiaceae, Sciadopityaceae, Podocarpaceae, Cephalotaxaceae, Taxaceae, Ephedraceae and Welwitschiaceae are in the first group, while Cupressaceae and Araucariaceae are mainly in the second one. The third and the fourth groups composed of Pinaceae species are characterized by the genera Larix, and Abies and Cedrus, respectively. Principal component and discriminant analyses and divisive hierarchical clustering analysis of the 43 Pinaceae species were also performed. A clear-cut separation of the genera Abies, Larix, and Cedrus from the other Pinaceae was evidenced. In addition, a mass analysis of the two main chloroplastic lipids from 14 gymnosperms was performed. The results point to a great originality in gymnosperms since in several species and contrary to the angiosperms, the amount of digalactosyldiacylglycerol exceeds that of monogalactosyldiacylglycerol.