Genomic vulnerability to rapid climate warming in a tree species with a long generation time.

The ongoing increase in global temperature affects biodiversity, especially in mountain regions where climate change is exacerbated. As sessile, long-lived organisms, trees are especially challenged in terms of adapting to rapid climate change. Here, we show that low rates of allele frequency shifts in Swiss stone pine (Pinus cembra) occurring near the treeline result in high genomic vulnerability to future climate warming, presumably due to the species' long generation time. Using exome sequencing data from adult and juvenile cohorts in the Swiss Alps, we found an average rate of allele frequency shift of 1.23 × 10-2 /generation (i.e. 40 years) at presumably neutral loci, with similar rates for putatively adaptive loci associated with temperature (0.96 × 10-2 /generation) and precipitation (0.91 × 10-2 /generation). These recent shifts were corroborated by forward-in-time simulations at neutral and adaptive loci. Additionally, in juvenile trees at the colonisation front we detected alleles putatively beneficial under a future warmer and drier climate. Notably, the observed past rate of allele frequency shift in temperature-associated loci was decidedly lower than the estimated average rate of 6.29 × 10-2 /generation needed to match a moderate future climate scenario (RCP4.5). Our findings suggest that species with long generation times may have difficulty keeping up with the rapid climate change occurring in high mountain areas and thus are prone to local extinction in their current main elevation range.

Genomic signatures of convergent adaptation to Alpine environments in three Brassicaceae species.

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole-genome pooled population sequencing to study genome-wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over-representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species-specific.

Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine.

It is generally accepted that the spatial distribution of neutral genetic diversity within a species' native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome-wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, FST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (He ) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.

Accumulation of Mutational Load at the Edges of a Species Range.

Why species have geographically restricted distributions is an unresolved question in ecology and evolutionary biology. Here, we test a new explanation that mutation accumulation due to small population size or a history of range expansion can contribute to restricting distributions by reducing population growth rate at the edge. We examined genomic diversity and mutational load across the entire geographic range of the North American plant Arabidopsis lyrata, including old, isolated populations predominantly at the southern edge and regions of postglacial range expansion at the northern and southern edges. Genomic diversity in intergenic regions declined toward distribution edges and signatures of mutational load in exon regions increased. Genomic signatures of mutational load were highly linked to phenotypically expressed load, measured as reduced performance of individual plants and lower estimated rate of population growth. The geographic pattern of load and the connection between load and population growth demonstrate that mutation accumulation reduces fitness at the edge and helps restrict species' distributions.

Improved recovery of ancient DNA from subfossil wood - application to the world's oldest Late Glacial pine forest.

Ancient DNA from historical and subfossil wood has a great potential to provide new insights into the history of tree populations. However, its extraction and analysis have not become routine, mainly because contamination of the wood with modern plant material can complicate the verification of genetic information. Here, we used sapwood tissue from 22 subfossil pines that were growing c. 13 000 yr bp in Zurich, Switzerland. We developed and evaluated protocols to eliminate surface contamination, and we tested ancient DNA authenticity based on plastid DNA metabarcoding and the assessment of post-mortem DNA damage. A novel approach using laser irradiation coupled with bleaching and surface removal was most efficient in eliminating contaminating DNA. DNA metabarcoding confirmed which ancient DNA samples repeatedly amplified pine DNA and were free of exogenous plant taxa. Pine DNA sequences of these samples showed a high degree of cytosine to thymine mismatches, typical of post-mortem damage. Stringent decontamination of wood surfaces combined with DNA metabarcoding and assessment of post-mortem DNA damage allowed us to authenticate ancient DNA retrieved from the oldest Late Glacial pine forest. These techniques can be applied to any subfossil wood and are likely to improve the accessibility of relict wood for genome-scale ancient DNA studies.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization.

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined FST and DXY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

Cysteine-rich angiogenic inducer 61 (Cyr61): a novel soluble biomarker of acute myocardial injury improves risk stratification after acute coronary syndromes.

AIMS: We aimed to identify a novel biomarker involved in the early events leading to an acute coronary syndrome (ACS) and evaluate its role in diagnosis and risk stratification. METHODS AND RESULTS: Biomarker identification was based on gene expression profiling. In coronary thrombi of ACS patients, cysteine-rich angiogenic inducer 61 (Cyr61, CCN1) gene transcripts were highly up-regulated compared with peripheral mononuclear cells. In a murine ischaemia-reperfusion model (I/R), myocardial Cyr61 expression was markedly increased compared with the controls. Cyr61 levels were determined in human serum using an enzyme-linked immunosorbent assay. Cohorts of ACS (n = 2168) referred for coronary angiography, stable coronary artery disease (CAD) (n = 53), and hypertrophic obstructive cardiomyopathy (HOCM) patients (n = 15) served to identify and evaluate the diagnostic and prognostic performance of the biomarker. Cyr61 was markedly elevated in ST-elevation myocardial infarction patients compared with non-ST-elevation myocardial infarction/unstable angina or stable CAD patients, irrespective of whether coronary thrombi were present. Cyr61 was rapidly released after occlusion of a septal branch in HOCM patients undergoing transcoronary ablation of septal hypertrophy. Cyr61 improved risk stratification for all-cause mortality when added to the reference GRACE risk score at 30 days (C-statistic 0.88 to 0.89, P = 0.001) and 1 year (C-statistic 0.77 to 0.80, P < 0.001) comparable to high-sensitivity troponin T (30 days: 0.88 to 0.89, P < 0.001; 1 year: 0.77 to 0.79, P < 0.001). Similar results were obtained for the composite endpoint of all-cause mortality or myocardial infarction. Conversely, in a population-based case-control cohort (n = 362), Cyr61 was not associated with adverse outcome. CONCLUSION: Cyr61 is a novel early biomarker reflecting myocardial injury that improves risk stratification in ACS patients.

High-density genetic mapping identifies the genetic basis of a natural colony morphology mutant in the root rot pathogen Armillaria ostoyae.

Filamentous fungi exhibit a broad spectrum of heritable growth patterns and morphological variations reflecting the adaptation of the different species to distinct ecological niches. But also within species, isolates show considerable variation in growth rates and other morphological characteristics. The genetic basis of this intraspecific variation in mycelial growth and morphology is currently poorly understood. By chance, a growth mutant in the root rot pathogen Armillaria ostoyae was discovered. The mutant phenotype was characterized by extremely compact and slow growth, as well as shorter aerial hyphae and hyphal compartments in comparison to the wildtype phenotype. Genetic analysis revealed that the abnormal phenotype is caused by a recessive mutation, which segregates asa single locus in sexual crosses. In order to identify the genetic basis of the mutant phenotype, we performed a quantitative trait locus (QTL) analysis. A mapping population of 198 haploid progeny was genotyped at 11,700 genome-wide single nucleotide polymorphisms (SNPs) making use of double digest restriction site associated DNA sequencing (ddRADseq). In accordance with the genetic analysis, a single significant QTL was identified for the abnormal growth phenotype. The QTL confidence interval spans a narrow, gene dense region of 87kb in the A. ostoyae genome which contains 37 genes. Overall, our study reports the first high-density genetic map for an Armillaria species and shows its successful application in forward genetics by resolving the genetic basis of a mutant phenotype with a severe defect in hyphal growth.

Maximum-Likelihood Tree Estimation Using Codon Substitution Models with Multiple Partitions.

Many protein sequences have distinct domains that evolve with different rates, different selective pressures, or may differ in codon bias. Instead of modeling these differences by more and more complex models of molecular evolution, we present a multipartition approach that allows maximum-likelihood phylogeny inference using different codon models at predefined partitions in the data. Partition models can, but do not have to, share free parameters in the estimation process. We test this approach with simulated data as well as in a phylogenetic study of the origin of the leucin-rich repeat regions in the type III effector proteins of the pythopathogenic bacteria Ralstonia solanacearum. Our study does not only show that a simple two-partition model resolves the phylogeny better than a one-partition model but also gives more evidence supporting the hypothesis of lateral gene transfer events between the bacterial pathogens and its eukaryotic hosts.

TRAL: tandem repeat annotation library.

MOTIVATION: Currently, more than 40 sequence tandem repeat detectors are published, providing heterogeneous, partly complementary, partly conflicting results. RESULTS: We present TRAL, a tandem repeat annotation library that allows running and parsing of various detection outputs, clustering of redundant or overlapping annotations, several statistical frameworks for filtering false positive annotations, and importantly a tandem repeat annotation and refinement module based on circular profile hidden Markov models (cpHMMs). Using TRAL, we evaluated the performance of a multi-step tandem repeat annotation workflow on 547 085 sequences in UniProtKB/Swiss-Prot. The researcher can use these results to predict run-times for specific datasets, and to choose annotation complexity accordingly. AVAILABILITY AND IMPLEMENTATION: TRAL is an open-source Python 3 library and is available, together with documentation and tutorials via http://www.vital-it.ch/software/tral. CONTACT: elke.schaper@isb-sib.ch.

Signatures of local adaptation in candidate genes of oaks (Quercus spp.) with respect to present and future climatic conditions.

Testing how populations are locally adapted and predicting their response to their future environment is of key importance in view of climate change. Landscape genomics is a powerful approach to investigate genes and environmental factors involved in local adaptation. In a pooled amplicon sequencing approach of 94 genes in 71 populations, we tested whether >3500 single nucleotide polymorphisms (SNPs) in the three most common oak species in Switzerland (Quercus petraea, Q. pubescens, Q. robur) show an association with abiotic factors related to local topography, historical climate and soil characteristics. In the analysis including all species, the most frequently associated environmental factors were those best describing the habitats of the species. In the species-specific analyses, the most important environmental factors and associated SNPs greatly differed among species. However, we identified one SNP and seven genes that were associated with the same environmental factor across all species. We finally used regressions of allele frequencies of the most strongly associated SNPs along environmental gradients to predict the risk of nonadaptedness (RONA), which represents the average change in allele frequency at climate-associated loci theoretically required to match future climatic conditions. RONA is considerable for some populations and species (up to 48% in single populations) and strongly differs among species. Given the long generation time of oaks, some of the required allele frequency changes might not be realistic to achieve based on standing genetic variation. Hence, future adaptedness requires gene flow or planting of individuals carrying beneficial alleles from habitats currently matching future climatic conditions.

Circulating FABP4 is a prognostic biomarker in patients with acute coronary syndrome but not in asymptomatic individuals.

OBJECTIVE: Blood-borne biomarkers reflecting atherosclerotic plaque burden have great potential to improve clinical management of atherosclerotic coronary artery disease and acute coronary syndrome (ACS). APPROACH AND RESULTS: Using data integration from gene expression profiling of coronary thrombi versus peripheral blood mononuclear cells and proteomic analysis of atherosclerotic plaque-derived secretomes versus healthy tissue secretomes, we identified fatty acid-binding protein 4 (FABP4) as a biomarker candidate for coronary artery disease. Its diagnostic and prognostic performance was validated in 3 different clinical settings: (1) in a cross-sectional cohort of patients with stable coronary artery disease, ACS, and healthy individuals (n=820), (2) in a nested case-control cohort of patients with ACS with 30-day follow-up (n=200), and (3) in a population-based nested case-control cohort of asymptomatic individuals with 5-year follow-up (n=414). Circulating FABP4 was marginally higher in patients with ST-segment-elevation myocardial infarction (24.9 ng/mL) compared with controls (23.4 ng/mL; P=0.01). However, elevated FABP4 was associated with adverse secondary cerebrovascular or cardiovascular events during 30-day follow-up after index ACS, independent of age, sex, renal function, and body mass index (odds ratio, 1.7; 95% confidence interval, 1.1-2.5; P=0.02). Circulating FABP4 predicted adverse events with similar prognostic performance as the GRACE in-hospital risk score or N-terminal pro-brain natriuretic peptide. Finally, no significant difference between baseline FABP4 was found in asymptomatic individuals with or without coronary events during 5-year follow-up. CONCLUSIONS: Circulating FABP4 may prove useful as a prognostic biomarker in risk stratification of patients with ACS.

Genetic evidence for reproductive isolation among sympatric Epichloë endophytes as inferred from newly developed microsatellite markers.

Reproductive isolation is central to the maintenance of species, and especially in sympatry, effective barriers to prevent interspecific crosses are expected. Host specificity is thought to constitute an effective mechanism for the formation of barriers in different genera of Fungi, but evidence for endophytes is so far lacking. Sexual Epichloë species (Ascomycota, Clavicipitaceae) represent an ideal study system to investigate the mechanisms underlying speciation as mediated by host specificity because they include species complexes with several host-specific taxa. Here, we studied genetic differentiation of three host-specific Epichloë species using microsatellite markers that were newly in silico identified on the genome of Epichloë poae. Among these, 15 were experimentally tested and applied to study an extensive sampling of isolates representing Epichloë typhina infecting Dactylis glomerata and Epichloë clarkii infecting Holcus lanatus from a site with sympatric populations in Switzerland, as well as a reduced sampling of E. poae infecting Poa nemoralis to create a three-taxon dataset. Both principal coordinate analysis and Bayesian clustering algorithm showed three genetically distinct groups representing the three host-specific species. High pairwise F ST values among the three species, as well as sequencing data of the tefA gene revealing diagnostic single nucleotide polymorphisms (SNPs), further support the hypothesis of genetic discontinuities among the taxa. These results provide genotypic evidence of the maintenance of reproductive isolation of the species in a context of sympatry. In silico testing of 885 discovered microsatellites on the genome of Epichloë festucae extend their applicability to a wider taxonomic range of Epichloë.

Improving phylogenetic inference with a semiempirical amino acid substitution model.

Amino acid substitution matrices describe the rates by which amino acids are replaced during evolution. In contrast to nucleotide or codon models, amino acid substitution matrices are in general parameterless and empirically estimated, probably because there is no obvious parametrization for amino acid substitutions. Principal component analysis has previously been used to improve codon substitution models by empirically finding the most relevant parameters. Here, we apply the same method to amino acid substitution matrices, leading to a semiempirical substitution model that can adjust the transition rates to the protein sequences under investigation. Our new model almost invariably achieves the best likelihood values in large-scale comparisons with established amino acid substitution models (JTT, WAG, and LG). In particular for longer alignments, these likelihood gains are considerably larger than what could be expected from simply having more parameters. The application of our model differs from that of mixture models (such as UL2 or UL3), as we optimize one rate matrix per alignment, whereas mixture models apply the variation per alignments site. This makes our model computationally more efficient, while the performance is comparable to that of UL3. Applied to the phylogenetic problem of the origin of placental mammals, our new model and the UL3 mixed model are the only ones of the tested models that cluster Afrotheria and Xenarthra into a clade called Atlantogenata, which would be in correspondence with recent findings using more sophisticated phylogenetic methods.

CodonPhyML: fast maximum likelihood phylogeny estimation under codon substitution models.

Markov models of codon substitution naturally incorporate the structure of the genetic code and the selection intensity at the protein level, providing a more realistic representation of protein-coding sequences compared with nucleotide or amino acid models. Thus, for protein-coding genes, phylogenetic inference is expected to be more accurate under codon models. So far, phylogeny reconstruction under codon models has been elusive due to computational difficulties of dealing with high dimension matrices. Here, we present a fast maximum likelihood (ML) package for phylogenetic inference, CodonPhyML offering hundreds of different codon models, the largest variety to date, for phylogeny inference by ML. CodonPhyML is tested on simulated and real data and is shown to offer excellent speed and convergence properties. In addition, CodonPhyML includes most recent fast methods for estimating phylogenetic branch supports and provides an integral framework for models selection, including amino acid and DNA models.

A new semiempirical codon substitution model based on principal component analysis of mammalian sequences.

Codon substitution models have traditionally been parametric Markov models, but recently, empirical and semiempirical models also have been proposed. Parametric codon models are typically based on 61×61 rate matrices that are derived from a small number of parameters. These parameters are rooted in experience and theoretical considerations and generally show good performance but are still relatively arbitrary. We have previously used principal component analysis (PCA) on data obtained from mammalian sequence alignments to empirically identify the most relevant parameters for codon substitution models, thereby confirming some commonly used parameters but also suggesting new ones. Here, we present a new semiempirical codon substitution model that is directly based on those PCA results. The substitution rate matrix is constructed from linear combinations of the first few (the most important) principal components with the coefficients being free model parameters. Thus, the model is not only based on empirical rates but also uses the empirically determined most relevant parameters for a codon model to adjust to the particularities of individual data sets. In comparisons against established parametric and semiempirical models, the new model consistently achieves the highest likelihood values when applied to sequences of vertebrates, which include the taxonomic class where the model was trained on.

Comparative genomics approach to detecting split-coding regions in a low-coverage genome: lessons from the chimaera Callorhinchus milii (Holocephali, Chondrichthyes).

Recent development of deep sequencing technologies has facilitated de novo genome sequencing projects, now conducted even by individual laboratories. However, this will yield more and more genome sequences that are not well assembled, and will hinder thorough annotation when no closely related reference genome is available. One of the challenging issues is the identification of protein-coding sequences split into multiple unassembled genomic segments, which can confound orthology assignment and various laboratory experiments requiring the identification of individual genes. In this study, using the genome of a cartilaginous fish, Callorhinchus milii, as test case, we performed gene prediction using a model specifically trained for this genome. We implemented an algorithm, designated ESPRIT, to identify possible linkages between multiple protein-coding portions derived from a single genomic locus split into multiple unassembled genomic segments. We developed a validation framework based on an artificially fragmented human genome, improvements between early and recent mouse genome assemblies, comparison with experimentally validated sequences from GenBank, and phylogenetic analyses. Our strategy provided insights into practical solutions for efficient annotation of only partially sequenced (low-coverage) genomes. To our knowledge, our study is the first formulation of a method to link unassembled genomic segments based on proteomes of relatively distantly related species as references.

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps.

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool-Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using FST -based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo-climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo-climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.

Mitochondrial genome evolution in species belonging to the Phialocephala fortinii s.l. - Acephala applanata species complex.

BACKGROUND: Mitochondrial (mt) markers are successfully applied in evolutionary biology and systematics because mt genomes often evolve faster than the nuclear genomes. In addition, they allow robust phylogenetic analysis based on conserved proteins of the oxidative phosphorylation system. In the present study we sequenced and annotated the complete mt genome of P. subalpina, a member of the Phialocephala fortinii s.l. - Acephala applanata species complex (PAC). PAC belongs to the Helotiales, which is one of the most diverse groups of ascomycetes including more than 2,000 species. The gene order was compared to deduce the mt genome evolution in the Pezizomycotina. Genetic variation in coding and intergenic regions of the mtDNA was studied for PAC to assess the usefulness of mt DNA for species diagnosis. RESULTS: The mt genome of P. subalpina is 43,742 bp long and codes for 14 mt genes associated with the oxidative phosphorylation. In addition, a GIY-YIG endonuclease, the ribosomal protein S3 (Rps3) and a putative N-acetyl-transferase were recognized. A complete set of tRNA genes as well as the large and small rRNA genes but no introns were found. All protein-coding genes were confirmed by EST sequences. The gene order in P. subalpina deviated from the gene order in Sclerotinia sclerotiorum, the only other helotialean species with a fully sequenced and annotated mt genome. Gene order analysis within Pezizomycotina suggests that the evolution of gene orders is mostly driven by transpositions. Furthermore, sequence diversity in coding and non-coding mtDNA regions in seven additional PAC species was pronounced and allowed for unequivocal species diagnosis in PAC. CONCLUSIONS: The combination of non-interrupted ORFs and EST sequences resulted in a high quality annotation of the mt genome of P. subalpina, which can be used as a reference for the annotation of other mt genomes in the Helotiales. In addition, our analyses show that mtDNA loci will be the marker of choice for future analysis of PAC communities.

Rust fungi on Annonaceae II: the genus Dasyspora Berk. & M.A. Curtis.

Dasyspora gregaria, the single species of the allegedly monotypic rust genus Dasyspora (Basidiomycota, Pucciniales), was investigated by light microscopy and DNA sequencing (ITS1-5.8S-ITS2 region, partial LSU and SSU of the nuclear rDNA, mt cytochrome oxidase subunit 3). Both methods indicated that D. gregaria is not a single species but can be split in 11 distinct taxa, each of which appear confined to a single Xylopia species (Annonaceae) host. Herein nine of these are described as new. Both the phylogenetic analyses and morphology show that the species are grouped into two main clades designated Dasyspora gregaria and D. winteri. The first comprises D. gregaria, the type species of the genus, which is restricted to X. cayennensis, two new species on X. aromatica, D. segregaria from northern South America and D. echinata from Brazil. The second clade is formed by D. winteri, recombined from Puccinia winteri on X. sericea, and the new species D. amazonica on X. amazonica, D. emarginatae on X. emarginata, D. frutescentis on X. frutescens, D. ferrugineae on X. frutescens var. ferruginea, D. guianensis on X. benthamii, D. mesoamericana on X. frutescens, and D. nitidae on X. nitida. Dasyspora frutescentis and D. mesoamericana were not clearly distinguishable by their morphology and host associations but differed from another in their sequences and geographic distributions. They are considered cryptic species. An identification key and the distributions are given for all recognized species. Along with molecular data we discuss the systematic position of Dasyspora in the Pucciniales.