Evolutionary trends of reproductive phenotype in Cycadales: an analysis of morphological evolution in Ceratozamia.

BACKGROUND AND AIMS: The size and shape of reproductive structures is especially relevant in evolution because these characters are directly related to the capacity of pollination and seed dispersal, a process that plays a basic role in evolutionary patterns. The evolutionary trajectories of reproductive phenotypes in gymnosperms have received special attention in terms of pollination and innovations related to the emergence of the Spermatophytes. However, variability of reproductive structures, evolutionary trends and the role of environment in the evolution of cycad species have not been well documented and explored. This study considered this topic under an explicitly phylogenetic and evolutionary approach that included a broad sampling of reproductive structures in the genus Ceratozamia. METHODS: We sampled 1400 individuals of 36 Ceratozamia species to explore evolutionary pattern and identify and evaluate factors that potentially drove their evolution. We analyzed characters for both pollen and ovulate strobili within a phylogenetic framework using different methods and characters (i. e., molecular and both quantitative and qualitative morphological) to infer phylogenetic relationships. Using this phylogenetic framework, evolutionary models of trait evolution for strobilar size were evaluated. In addition, quantitative morphological variation and its relation to environmental variables across species were analyzed. KEY RESULTS: We found contrasting phylogenetic signals between characters of pollen and ovulate strobili. These structures exhibited high morphological disparity in several characters related to size. Results of analyses of evolutionary trajectories suggested a stabilizing selection model. In regards to phenotype-environment, the analysis produced mixed results and differences for groups in the vegetation type where the species occur; however, a positive relationship with climatic variables was found. CONCLUSIONS: The integrated approach synthesized reproductive phenotypic variation with current phylogenetic hypotheses and provided explicit statements of character evolution. The characters of volume for ovulate strobili were the most informative, which could provide a reference for further study of the evolutionary complexity in Ceratozamia. Finally, heterogeneous environments, which are under changing weather conditions, promote variability of reproductive structures.

Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae).

BACKGROUND AND AIMS: Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico. METHODS: Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages. KEY RESULTS: The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages. CONCLUSIONS: The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

Stochastic processes and changes in evolutionary rate are associated with diversification in a lineage of tropical hard pines (Pinus).

The study of the patterns of polymorphism and molecular evolution among closely related species is key to understanding the evolutionary forces involved in the diversification of lineages. This point is a big challenge in species with slow evolutionary rates, long life cycles, and ancient, shared polymorphisms such as conifers. Under the premise of divergence in a stepwise migration process, we expect clinal geographical patterns of purifying selection efficiency, and genetic structure related to latitude or longitude. If migration is accompanied by changes in the environment, we could further expect a role of positive selection in driving species divergence. Here, we infer patterns of polymorphism, efficiency of purifying selection, and molecular evolution using a dataset of 161 nuclear genes (∼71 Kb) in a lineage of hard pines from North America, the Caribbean, Mexico, and Central America presumed to have migrated from North America toward lower latitudes with tropical conditions. Under the premise of differences in selective pressures, we also look for possible signals of positive selection. To test our hypothesis, first we estimated different indices to infer patterns of polymorphism and efficiency of purifying selection (Ka, Ks, Ka/Ks, dN, dS, dN/dS, and dxy) and compared these metrics across five clades. Also, we investigated possible clinal patterns in these indices and morphological traits (needle length and cone length). Then we inferred genetic structure and environmental differences among species to test for possible signals of positive selection using phylogenetic methods in specific clades. We found differences among clades using Ka, Ks, and Ka/Ks with a relaxation of purifying selection, especially in the Elliotti and Patula clades. We also found environmental differences related to geographic distance, and among clades suggesting differences in selective pressures. The indices Ks, dxy, and needle length had relationships with geography but not ovulate cone length. Finally, we found that most analyzed genes are under purifying selection, but there was an exception of faster evolutionary rate in some pine species, suggesting the possible action of positive selection in divergence. Our study indicated that stochastic processes have played a key role in the diversification of the group, with a possible input of positive selection in pines from Mexico and Central America.

Changes in demography and geographic distribution in the weeping pinyon pine (Pinus pinceana) during the Pleistocene.

Climate changes, together with geographical barriers imposed by the Sierra Madre Oriental and the Chihuahuan Desert, have shaped the genetic diversity and spatial distribution of different species in northern Mexico. Pinus pinceana Gordon & Glend. tolerates extremely arid conditions. Northern Mexico became more arid during the Quaternary, modifying ecological communities. Here, we try to identify the processes underlying the demographic history of P. pinceana and characterize its genetic diversity using 3100 SNPs from genotyping by sequencing 90 adult individuals from 10 natural populations covering the species' entire geographic distribution. We inferred its population history and contrasted possible demographic scenarios of divergence that modeled the genetic diversity present in this restricted pinyon pine; in support, the past distribution was reconstructed using climate from the Last Glacial Maximum (LGM, 22 kya). We inferred that P. pinceana diverged into two lineages ~2.49 Ma (95% CI 3.28-1.62), colonizing two regions: the Sierra Madre Oriental (SMO) and the Chihuahuan Desert (ChD). Our results of population genomic analyses reveal the presence of heterozygous SNPs in all populations. In addition, low migration rates across regions are probably related to glacial-interglacial cycles, followed by the gradual aridification of the Chihuahuan Desert during the Holocene.

Coalescent-based species delimitation in North American pinyon pines using low-copy nuclear genes and plastomes.

PREMISE: Accurate species delimitation is essential for evolutionary biology, conservation, and biodiversity management. We studied species delimitation in North American pinyon pines, Pinus subsection Cembroides, a natural group with high levels of incomplete lineage sorting. METHODS: We used coalescent-based methods and multivariate analyses of low-copy number nuclear genes and nearly complete high-copy number plastomes generated with the Hyb-Seq method. The three coalescent-based species delimitation methods evaluated were the Generalized Mixed Yule Coalescent (GMYC), Poisson Tree Process (PTP), and Trinomial Distribution of Triplets (Tr2). We also measured admixture in populations with possible introgression. RESULTS: Our results show inconsistencies among GMYC, PTP, and Tr2. The single-locus based GMYC analysis of plastid DNA recovered a higher number of species (up to 24 entities, including singleton lineages and clusters) than PTP and the multi-locus coalescent approach. The PTP analysis identified 10 species whereas Tr2 recovered 13, which agreed closely with taxonomic treatments. CONCLUSIONS: We found that PTP and GMYC identified species with low levels of ILS and high morphological divergence (P. maximartinezii, P. pinceana, and P. rzedowskii). However, GMYC method oversplit species by identification of more divergent samples as singletons. Moreover, both PTP and GMYC were incapable of identifying some species that are readily identified morphologically. We suggest that the divergence times between lineages within North American pinyon pines are so disparate that GMYC results are unreliable. Results of the Tr2 method coincided well with previous delimitations based on morphology, DNA, geography, and secondary chemistry.

Deep sequencing across multiple host species tests pine-endophyte specificity.

PREMISE: Foliar fungal endophytes vary in their distributions across landscapes or plant host taxa, indicative of specialized ecologies and host specific adaptations. Accounts of specialization, however, depend on the taxonomic breadth and geographic range of the host plants included in each study. A broad region-scale study or deep sampling of diverse potential host species still remains relatively rare but is becoming increasingly possible with high-throughput sequencing. METHODS: Amplicon sequencing was used to rapidly identify the fungal endophytic community among six pine (Pinus, Pinaceae) species co-occurring across northeastern United States and to test for site and host specialization. We focused on the endophytic genus Lophodermium (Rhytismataceae), whose species members are thought to specialize on different pine species, to test if amplicon sequencing could rapidly verify previously implied or discover new patterns of host specificity. RESULTS: While amplicon sequencing could analyze more samples at greater depths and recover greater numbers of unique Lophodermium taxa than when endophyte communities were surveyed with traditional culturing methods, patterns of specialization were not better supported. This may be because amplicon sequencing can indiscriminately capture non-host specific organisms found incidentally from plant tissues or because we have overestimated host-specificity in the past with biased culturing techniques. CONCLUSIONS: Amplicon sequencing can quickly identify patterns of host specificity by allowing large-scale surveys but has limitations in quantifying the level of intimacy of these relationships.

Histoplasma capsulatum Isolated from Tadarida brasiliensis Bats Captured in Mexico Form a Sister Group to North American Class 2 Clade.

Histoplasma capsulatum is a dimorphic fungus associated with respiratory and systemic infections in mammalian hosts that have inhaled infective mycelial propagules. A phylogenetic reconstruction of this pathogen, using partial sequences of arf, H-anti, ole1, and tub1 protein-coding genes, proposed that H. capsulatum has at least 11 phylogenetic species, highlighting a clade (BAC1) comprising three H. capsulatum isolates from infected bats captured in Mexico. Here, relationships for each individual locus and the concatenated coding regions of these genes were inferred using parsimony, maximum likelihood, and Bayesian inference methods. Coalescent-based analyses, a concatenated sequence-types (CSTs) network, and nucleotide diversities were also evaluated. The results suggest that six H. capsulatum isolates from the migratory bat Tadarida brasiliensis together with one isolate from a Mormoops megalophylla bat support a NAm 3 clade, replacing the formerly reported BAC1 clade. In addition, three H. capsulatum isolates from T. brasiliensis were classified as lineages. The concatenated sequence analyses and the CSTs network validate these findings, suggesting that NAm 3 is related to the North American class 2 clade and that both clades could share a recent common ancestor. Our results provide original information on the geographic distribution, genetic diversity, and host specificity of H. capsulatum.

Phylogenomic and ecological analyses reveal the spatiotemporal evolution of global pines.

How coniferous forests evolved in the Northern Hemisphere remains largely unknown. Unlike most groups of organisms that generally follow a latitudinal diversity gradient, most conifer species in the Northern Hemisphere are distributed in mountainous areas at middle latitudes. It is of great interest to know whether the midlatitude region has been an evolutionary cradle or museum for conifers and how evolutionary and ecological factors have driven their spatiotemporal evolution. Here, we investigated the macroevolution of Pinus, the largest conifer genus and characteristic of northern temperate coniferous forests, based on nearly complete species sampling. Using 1,662 genes from transcriptome sequences, we reconstructed a robust species phylogeny and reestimated divergence times of global pines. We found that ∼90% of extant pine species originated in the Miocene in sharp contrast to the ancient origin of Pinus, indicating a Neogene rediversification. Surprisingly, species at middle latitudes are much older than those at other latitudes. This finding, coupled with net diversification rate analysis, indicates that the midlatitude region has provided an evolutionary museum for global pines. Analyses of 31 environmental variables, together with a comparison of evolutionary rates of niche and phenotypic traits with a net diversification rate, found that topography played a primary role in pine diversification, and the aridity index was decisive for the niche rate shift. Moreover, fire has forced diversification and adaptive evolution of Pinus Our study highlights the importance of integrating phylogenomic and ecological approaches to address evolution of biological groups at the global scale.

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.

Population structure, diversifying selection, and local adaptation in Pinus patula.

PREMISE: Climate change is predicted to affect natural and plantation forests. The responses of conifers to overcome changing environments will depend on their adaptation to local conditions; however, intraspecific adaptive genetic variation is unknown for most gymnosperms. Studying genetic diversity associated with phenotypic variability along environmental gradients will enhance our understanding of adaptation and may reveal genetic pools important for conservation and management. METHODS: We used target enrichment and genome skimming to obtain single nucleotide polymorphisms (SNPs) from 61 individuals of Pinus patula, a pine tree native to Mexico widely used in plantation forestry. We investigated the adaptive genetic variation of two varieties with morphological and distributional differences potentially related to genetic and adaptive divergence. RESULTS: Population structure and haplotype network analyses revealed that genetic diversity between P. patula var. patula and P. patula var. longipedunculata was structured, even within populations of P. patula var. longipedunculata. We observed high genetic diversity, low inbreeding rate, and rapid linkage disequilibrium (LD) decay in the varieties. Based on outlier tests, loci showing signatures of natural selection were detected in geographically distant P. patula var. longipedunculata populations. For both varieties, we found significant correlations between climate-related environmental variation and SNP diversity at loci involved in abiotic stress, cell transport, defense, and cell wall biogenesis, pointing to local adaptation. CONCLUSIONS: Overall, significant intraspecific adaptive genetic variation in P. patula was detected, highlighting the presence of different genetic pools and signs of local adaptation that should be considered in forestry and conservation.

Incorporating fossils into the Pinaceae tree of life.

PREMISE OF THE STUDY: Pinaceae have a rich but enigmatic early fossil record, much of which is represented by permineralized seed cones. Our incomplete knowledge of morphology and anatomy in living and extinct species poses an important barrier to understanding their phylogenetic relationships and timing of diversification. METHODS: We expanded a morphology matrix to 46 fossil and 31 extant Pinaceae species, mainly adding characters from stem and leaf anatomy and seed cones. Using parsimony and Bayesian inference, we compared phylogenetic relationships for extant taxa with and without fossils from the morphology matrix combined with an alignment of plastid gene sequences. KEY RESULTS: Combined analysis of morphological and molecular characters resulted in a phylogeny of extant Pinaceae that was robust at all nodes except those relating to the interrelationships of Pinus, Picea, and Cathaya and the position of Cedrus. Simultaneous analysis of all fossil and extant species did not result in changes in the relationships among the extant species but did greatly reduce branch support. We found that the placement of most fossils was sensitive to the method of phylogenetic reconstruction when analyzing them singly with the extant species. CONCLUSIONS: A robust phylogenetic hypothesis for the main lineages of Pinaceae is emerging. Most Early Cretaceous fossils are stem or crown lineages of Pinus, but close relationships also were found between fossils and several other extant genera. The phylogenetic position of fossils broadly supports the existence of extant genera in the Lower Cretaceous.

Multi-locus phylogenetics, lineage sorting, and reticulation in Pinus subsection Australes.

PREMISE OF THE STUDY: Both incomplete lineage sorting and reticulation have been proposed as causes of phylogenetic incongruence. Disentangling these factors may be most difficult in long-lived, wind-pollinated plants with large population sizes and weak reproductive barriers. METHODS: We used solution hybridization for targeted enrichment and massive parallel sequencing to characterize low-copy-number nuclear genes and high-copy-number plastomes (Hyb-Seq) in 74 individuals of Pinus subsection Australes, a group of ~30 New World pine species of exceptional ecological and economic importance. We inferred relationships using methods that account for both incomplete lineage sorting and reticulation. KEY RESULTS: Concatenation- and coalescent-based trees inferred from nuclear genes mainly agreed with one another, but they contradicted the plastid DNA tree in recovering the Attenuatae (the California closed-cone pines) and Oocarpae (the egg-cone pines of Mexico and Central America) as monophyletic and the Australes sensu stricto (the southern yellow pines) as paraphyletic to the Oocarpae. The plastid tree featured some relationships that were discordant with morphological and geographic evidence and species limits. Incorporating gene flow into the coalescent analyses better fit the data, but evidence supporting the hypothesis that hybridization explains the non-monophyly of the Attenuatae in the plastid tree was equivocal. CONCLUSIONS: Our analyses document cytonuclear discordance in Pinus subsection Australes. We attribute this discordance to ancient and recent introgression and present a phylogenetic hypothesis in which mostly hierarchical relationships are overlain by gene flow.

Complete plastomes of three endemic Mexican pine species (Pinus subsection Australes).

We assembled the plastomes of Pinus greggii, P. jaliscana and P. oocarpa from 100 bp paired-end Illumina reads. We combined de novo (comparing Velvet and SPAdes) with reference-guided assembly and a final step of gap filling. SPAdes performed better than Velvet based on scaffold number (180 vs. 263) and mean length (1886 vs. 560 bp), and number of gaps (2 vs. 4). Annotations were automatically transferred from P. taeda NC_021440 and carefully revised by hand. Phylogenetic analysis with additional plastomes revealed very short branch lengths, supporting a rapid diversification within Australes and close relatedness among pines from Western Mexico.

Pinus ponderosa: A checkered past obscured four species.

PREMISE OF THE STUDY: Molecular genetic evidence can help delineate taxa in species complexes that lack diagnostic morphological characters. Pinus ponderosa (Pinaceae; subsection Ponderosae) is recognized as a problematic taxon: plastid phylogenies of exemplars were paraphyletic, and mitochondrial phylogeography suggested at least four subdivisions of P. ponderosa. These patterns have not been examined in the context of other Ponderosae species. We hypothesized that putative intraspecific subdivisions might each represent a separate taxon. METHODS: We genotyped six highly variable plastid simple sequence repeats in 1903 individuals from 88 populations of P. ponderosa and related Ponderosae (P. arizonica, P. engelmannii, and P. jeffreyi). We used multilocus haplotype networks and discriminant analysis of principal components to test clustering of individuals into genetically and geographically meaningful taxonomic units. KEY RESULTS: There are at least four distinct plastid clusters within P. ponderosa that roughly correspond to the geographic distribution of mitochondrial haplotypes. Some geographic regions have intermixed plastid lineages, and some mitochondrial and plastid boundaries do not coincide. Based on relative distances to other species of Ponderosae, these clusters diagnose four distinct taxa. CONCLUSIONS: Newly revealed geographic boundaries of four distinct taxa (P. benthamiana, P. brachyptera, P. scopulorum, and a narrowed concept of P. ponderosa) do not correspond completely with taxonomies. Further research is needed to understand their morphological and nuclear genetic makeup, but we suggest that resurrecting originally published species names would more appropriately reflect the taxonomy of this checkered classification than their current treatment as varieties of P. ponderosa.

Phylogenetic relationships and species delimitation in pinus section trifoliae inferrred from plastid DNA.

Recent diversification followed by secondary contact and hybridization may explain complex patterns of intra- and interspecific morphological and genetic variation in the North American hard pines (Pinus section Trifoliae), a group of approximately 49 tree species distributed in North and Central America and the Caribbean islands. We concatenated five plastid DNA markers for an average of 3.9 individuals per putative species and assessed the suitability of the five regions as DNA bar codes for species identification, species delimitation, and phylogenetic reconstruction. The ycf1 gene accounted for the greatest proportion of the alignment (46.9%), the greatest proportion of variable sites (74.9%), and the most unique sequences (75 haplotypes). Phylogenetic analysis recovered clades corresponding to subsections Australes, Contortae, and Ponderosae. Sequences for 23 of the 49 species were monophyletic and sequences for another 9 species were paraphyletic. Morphologically similar species within subsections usually grouped together, but there were exceptions consistent with incomplete lineage sorting or introgression. Bayesian relaxed molecular clock analyses indicated that all three subsections diversified relatively recently during the Miocene. The general mixed Yule-coalescent method gave a mixed model estimate of only 22 or 23 evolutionary entities for the plastid sequences, which corresponds to less than half the 49 species recognized based on morphological species assignments. Including more unique haplotypes per species may result in higher estimates, but low mutation rates, recent diversification, and large effective population sizes may limit the effectiveness of this method to detect evolutionary entities.

Plastid DNA diversity is higher in the island endemic Guadalupe cypress than in the continental Tecate cypress.

BACKGROUND: Callitropsis guadalupensis (Guadalupe cypress) is endemic to Guadalupe Island, Mexico, where it is the dominant species of the only forest. The species has suffered declining numbers following the introduction of goats to the island over 150 years ago. Callitropsis guadalupensis is closely related to Callitropsis forbesii (Tecate cypress), distributed in small isolated populations in mainland Baja California and southern California. The objective of the present study was to compare the genetic diversity of the island endemic to the continental species. METHODOLOGY/PRINCIPAL FINDINGS: We measured genetic diversity in Callitropsis guadalupensis (n =54) from Guadalupe Island and in Callitropsis forbesii (n = 100) from five populations in mainland Baja California. The plastid DNA trnS-trnG spacer and the trnL-trnF region were chosen for characterization. Thirty-four haplotypes were observed, of which six were shared between both species. One of these haplotypes was also shared with three other species, Callitropsis lusitanica, Callitropsis montana, and Callitropsis stephensonii. Haplotype diversity (h) and nucleotide diversity (π) were significantly higher for Callitropsis guadalupensis (h = 0.698, π = 0.00071) than for Callitropsis forbesii (h = 0.337, π = 0.00024). CONCLUSIONS/SIGNIFICANCE: Callitropsis guadalupensis shows no evidence of a founder effect or of a genetic bottleneck, and can be added to a growing list of insular species with higher genetic diversity than their mainland relatives.

Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology.

Organellar DNA sequences are widely used in evolutionary and population genetic studies, however, the conservative nature of chloroplast gene and genome evolution often limits phylogenetic resolution and statistical power. To gain maximal access to the historical record contained within chloroplast genomes, we have adapted multiplex sequencing-by-synthesis (MSBS) to simultaneously sequence multiple genomes using the Illumina Genome Analyzer. We PCR-amplified approximately 120 kb plastomes from eight species (seven Pinus, one Picea) in 35 reactions. Pooled products were ligated to modified adapters that included 3 bp indexing tags and samples were multiplexed at four genomes per lane. Tagged microreads were assembled by de novo and reference-guided assembly methods, using previously published Pinus plastomes as surrogate references. Assemblies for these eight genomes are estimated at 88-94% complete, with an average sequence depth of 55x to 186x. Mononucleotide repeats interrupt contig assembly with increasing repeat length, and we estimate that the limit for their assembly is 16 bp. Comparisons to 37 kb of Sanger sequence show a validated error rate of 0.056%, and conspicuous errors are evident from the assembly process. This efficient sequencing approach yields high-quality draft genomes and should have immediate applicability to genomes with comparable complexity.

Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for pinus.

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.

Phylogenetics of Lophodermium from pine.

Lophodermium comprises ascomycetous fungi that are both needle-cast pathogens and asymptomatic endophytes on a diversity of plant hosts. It is distinguished from other genera in the family Rhytismataceae by its filiform ascospores and ascocarps that open by a longitudinal slit. Nucleotide sequences of the internal transcribed spacer (ITS) region of nuclear ribosomal DNA were used to infer phylogenetic relationships within Lophodermium. Twenty-nine sequences from approximately 11 species of Lophodermium were analyzed together with eight sequences from isolates thought to represent six other genera of Rhytismataceae: Elytroderma, Lirula, Meloderma, Terriera, Tryblidiopsis and Colpoma. Two putative Meloderma desmazieresii isolates occurred within the Lophodermium clade but separate from one another, one grouped with L. indianum and the other with L. nitens. An isolate of Elytroderma deformans also occurred within the Lophodermium clade but on a solitary branch. The occurrence of these genera within the Lophodermium clade might be due to problems in generic concepts in Rhytismataceae, such as emphasis on spore morphology to delimit genera, to difficulty of isolating Rhytismataceae needle pathogens from material that also is colonized by Lophodermium or to a combination of both factors. We also evaluated the congruence of host distribution and several morphological characters on the ITS phylogeny. Lophodermium species from pine hosts formed a monophyletic sister group to Lophodermium species from more distant hosts from the southern hemisphere, but not to L. piceae from Picea. The ITS topology indicated that Lophodermium does not show strict cospeciation with pines at deeper branches, although several closely related isolates have closely related hosts. Pathogenic species occupy derived positions in the pine clade, suggesting that pathogenicity has evolved from endophytism. A new combination is proposed, Terriera minor (Tehon) P.R. Johnst.