Genetic resources of common ash (Fraxinus excelsior L.) in Poland

BACKGROUND

Knowledge of genetic structure and the factors that shape it has an impact on forest management practices. European ash (Fraxinus excelsior L.) has declined dramatically throughout its range as a result of a disease caused by the fungus Hymenoscyphus fraxineus. Despite the need for conservation and restoration of the species, genetic data required to guide these efforts at the country level are scarce. Thereofore, we studied the chloroplast and nuclear genetic diversity of 26 natural common ash populations (1269 trees) in Poland.

RESULTS

Chloroplast polymorphisms grouped the populations into two geographically structured phylogenetic lineages ascribed to different glacial refugia (the Balkans and the Eastern Alps). However, the populations demonstrated high genetic diversity (mean AR = 12.35; mean Ho = 0.769; mean He = 0.542) but low differentiation based on nuclear microsatellites (FST = 0.045). Significant spatial genetic structure, consistent with models of isolation by distance, was detected in 14 out of 23 populations. Estimated effective population size was moderate-to-high, with a harmonic mean of 57.5 individuals per population.

CONCLUSIONS

Genetic diversity was not homogeneously distributed among populations within phylogenetic gene pools, indicating that ash populations are not equal as potential sources of reproductive material. Genetic differences among populations could be related to their histories, including founder effects or gene flow between evolutionary lineages (admixture). Our results suggest that ash stands across Poland could be treated as two main management units (seed zones). Therefore, despite the homogenizing effect of pollen gene flow known for this species, the genetic structure should be taken into account in the management of the genetic resources of the common ash. Although ash dieback poses an additional challenge for the management of genetic resources, efforts should be directed towards protecting populations with high genetic diversity within defined phylogenetic units, as they may be an important source of adaptive variation for future stands.

High quality genomes produced from single MinION flow cells clarify polyploid and demographic histories of critically endangered Fraxinus (ash) species

With populations of threatened and endangered species declining worldwide, efforts are being made to generate high quality genomic records of these species before they are lost forever. Here, we demonstrate that data from single Oxford Nanopore Technologies (ONT) MinION flow cells can, even in the absence of highly accurate short DNA-read polishing, produce high quality de novo plant genome assemblies adequate for downstream analyses, such as synteny and ploidy evaluations, paleodemographic analyses, and phylogenomics. This study focuses on three North American ash tree species in the genus Fraxinus (Oleaceae) that were recently added to the International Union for Conservation of Nature (IUCN) Red List as critically endangered. Our results support a hexaploidy event at the base of the Oleaceae as well as a subsequent whole genome duplication shared by Syringa, Osmanthus, Olea, and Fraxinus. Finally, we demonstrate the use of ONT long-read sequencing data to reveal patterns in demographic history.

The chromosome-level genome assembly of Fraxinus americana provides insights into the evolution of Oleaceae plants

White ash (Fraxinus americana linn.) originates from the southeastern United States. It is a tall and fast-growing tree species with strong salt-alkali resistance and cold tolerance, making it an important reforestation species and widely planted worldwide. Here, we completed the chromosome-level reference genome assembly of F. americana based on Illumina, PacBio, and Hi-C reads, with a genome size of 878.98 Mb, an N50 of 3.27 Mb, and a heterozygosity rate of 0.3 %. Based on de novo prediction, transcriptome prediction, and homology-based protein prediction, we obtained 39,538 genes. Approximately 843.21 Mb of the assembly genome was composed of 37,928 annotated protein-coding genes, with a gene function annotation rate of 95.93 %. 99.94 % of the overlap clusters (877.44 Mb) were anchored to 23 chromosomes. Synteny analysis of F. americana and other Oleaceae plants showed that F. americana underwent frequent chromosome rearrangements. The amplification of the Ale transposons effectively promoted the genome size of F. americana. Compared with other Oleaceae plants, the Glutathione S-transferase (GST) gene family in the F. americana genome has undergone significant expansion, which may help F. americana cope with adverse natural environments. Furthermore, we found that key enzyme-coding gene families related to lignin biosynthesis were expanded and highly expressed in F. americana leaves. These key genes drive lignin synthesis and benefit F. americana in fast-growing, as well as resisting biotic and abiotic stress. Overall, the F. americana genome assembly provides insights into the evolution of Oleaceae plants and provides abundant resources for breeding and germplasm conservation of white ash.

A glimmer of hope - ash genotypes with increased resistance to ash dieback pathogen show cross-resistance to emerald ash borer

Plants rely on cross-resistance traits to defend against multiple, phylogenetically distinct enemies. These traits are often the result of long co-evolutionary histories. Biological invasions can force naïve plants to cope with novel, coincident pests, and pathogens. For example, European ash (Fraxinus excelsior) is substantially threatened by the emerald ash borer (EAB), Agrilus planipennis, a wood-boring beetle, and the ash dieback (ADB) pathogen, Hymenoscyphus fraxineus. Yet, plant cross-resistance traits against novel enemies are poorly explored and it is unknown whether naïve ash trees can defend against novel enemy complexes via cross-resistance mechanisms. To gain mechanistic insights, we quantified EAB performance on grafted replicates of ash genotypes varying in ADB resistance and characterized ash phloem chemistry with targeted and untargeted metabolomics. Emerald ash borer performed better on ADB-susceptible than on ADB-resistant genotypes. Moreover, changes in EAB performance aligned with differences in phloem chemical profiles between ADB-susceptible and ADB-resistant genotypes. We show that intraspecific variation in phloem chemistry in European ash can confer increased cross-resistance to invasive antagonists from different taxonomic kingdoms. Our study suggests that promotion of ADB-resistant ash genotypes may simultaneously help to control the ADB disease and reduce EAB-caused ash losses, which may be critical for the long-term stability of this keystone tree species.

Transcriptome profiling of Fraxinus excelsior genotypes infested by emerald ash borer

European ash, Fraxinus excelsior is facing the double threat of ongoing devastation by the invasive fungal pathogen, Hymenoscyphus fraxineus and the imminent arrival of the non-native emerald ash borer (EAB), Agrilus planipennis. The spread of EAB which is currently moving westwards from European Russia and Ukraine into central Europe, poses an additional substantial threat to European ash, F. excelsior. While the molecular basis for resistance or variation in resistance among European ash genotypes is heavily investigated, comparatively little is known about the molecular ash traits involved in resistance against EAB. In this study we have gathered transcriptomic data from EAB inoculated genotypes of F. excelsior that have previously shown different levels of susceptibility to EAB. Resultant datasets show differential gene expression in susceptible and resistant genotypes in response to EAB infestation. This data will provide important information on the molecular basis of resistance to the EAB and allow the development of management plans to combat a pending threat of a culturally and ecologically important European tree species.

Traces of Hymenoscyphus fraxineus in Northeastern Europe Extend Further Back in History than Expected

Herbaria are a promising but still poorly applied information source for retrospective microbiological studies. In order to find any evidence of the virulent European origin of ash dieback agent Hymenoscyphus fraxineus and other fungal pathogens, we analyzed 109 leaf samples from three different Estonian botanical herbaria, sampled during 171 years from 20 ash species and cultivars, using a PacBio third-generation sequencing of the fungal internal transcribed spacer ITS1-5.8S-ITS2 ribosomal DNA region. We identified a large amount of saprotrophic fungi naturally colonizing ash leaves. Hymenoscyphus fraxineus colonized a Fraxinus chinensis subsp. rhynchophylla specimen and a F. chinensis specimen collected from Tallinn Botanic Garden in July 1978 and July 1992, respectively. The samples originated from trees grown in this garden from seeds collected from Shamora, Far-East Russia, in 1961 and from a Beijing botanical garden in eastern China in 1985, respectively. Repeated subsequent DNA extraction, real-time quantitative PCR, and Sanger and Illumina sequencing confirmed our findings of these apparently oldest cases of the ash dieback agent in Europe. These results show that H. fraxineus evidently was present in Estonia 19 years earlier than our previous data from fungal herbaria documented and 14 years before the first visible damage of ash trees was registered in Poland. Because we found no evidence of the saprotrophic H. albidus from earlier mycological and botanical herbarium specimens, the presence of H. albidus in Estonia remains questionable.

Genome-wide analysis and molecular dissection of the SPL gene family in Fraxinus mandshurica

BACKGROUND

SQUAMOSA promoter binding protein-like (SPL) is a unique family of transcription factors in plants, which is engaged in regulating plant growth and development, physiological and biochemical processes. Fraxinus mandshurica is an excellent timber species with a wide range of uses in northeastern China and enjoys a high reputation in the international market. SPL family analysis has been reported in some plants while SPL family analysis of Fraxinus mandshurica has not been reported.

RESULTS

We used phylogeny, conserved motifs, gene structure, secondary structure prediction, miR156 binding sites, promoter cis elements and GO annotation to systematically analyze the FmSPLs family. This was followed by expression analysis by subcellular localization, expression patterns at various tissue sites, abiotic stress and hormone induction. Because FmSPL2 is highly expressed in flowers it was selected to describe the SPL gene family of Fraxinus mandshurica by ectopic expression. Among them, 10 FmSPL genes that were highly expressed at different loci were selected for expression analysis under abiotic stress (NaCl and Cold) and hormone induction (IAA and ABA). These 10 FmSPL genes showed corresponding trends in response to both abiotic stress and hormone induction. We showed that overexpression of FmSPL2 in transgenic Nicotiana tabacum L. resulted in taller plants, shorter root length, increased root number, rounded leaves, and earlier flowering time.

CONCLUSIONS

We identified 36 SPL genes, which were classified into seven subfamilies based on sequence analysis. FmSPL2 was selected for subsequent heterologous expression by analysis of expression patterns in various tissues and under abiotic stress and hormone induction, and significant phenotypic changes were observed in the transgenic Nicotiana tabacum L. These results provide insight into the evolutionary origin and biological significance of plant SPL. The aim of this study was to lay the foundation for the genetic improvement of Fraxinus mandshurica and the subsequent functional analysis of FmSPL2.

Heterologous expression of a Fraxinus velutina SnRK2 gene in Arabidopsis increases salt tolerance by modifying root development and ion homeostasis

FvSnRK2182 is involved in regulating the growth and stress response. SnRK2 family members are positive regulators of downstream signals in the abscisic acid (ABA) signaling pathway, playing key roles in the plant responses to abiotic stresses. Fraxinus velutina Torr. is a candidate phytoremediator of saline-alkali areas, and is a valuable research subject because of its adaptability in saline soil. We identified a SnRK2 gene in F. velutina (named FvSnRK2182), which was significantly upregulated under salt stress. A bioinformatics analysis showed that FvSnRK2182 has a Ser/Thr kinase domain typical of the SnRK2 subfamily. Compared with wild-type (WT) Arabidopsis, its heterologous expression in Arabidopsis resulted in higher auxin content during seed germination and seedling growth, leading to longer primary roots and more lateral roots. The transgenic lines were better able to tolerate treatments with NaCl (100 mM) and/or ABA (0.2 and 0.5 µM), producing a greater biomass than the WT plants. Under NaCl treatment, the shoots of the transgenic lines had lower Na⁺ contents and higher K⁺ contents than the WT plants, and the genes encoding the ion transport-related proteins SOS1, HKT1, NHX1, and AKT1 were significantly upregulated. In addition, the expression of the genes functioning downstream of SnRK2 in the ABA signaling pathway (Rboh, AREB4, ABF2, and ABF3) were significantly upregulated in transgenic lines under NaCl stress. These results showed that expressing FvSnRK2182 in Arabidopsis significantly increased their resistance to ABA and salt stress by regulating root development and maintaining ion homeostasis, which suggests that FvSnRK2182 may be involved in regulating the growth and stress response of F. velutina.

First Report of Diplodia fraxini and Diplodia subglobosa Causing Canker and Dieback of Fraxinus excelsior in Slovenia

In recent decades the vitality and productivity of European ash trees in Slovenia have been reduced by the onset of canker and dieback disease symptoms on young and old trees, identified primarily as ash dieback caused by Hymenoscyphus fraxineus. Given the limited information available about the etiology of this emerging disease, a study was carried out to isolate, identify, and characterize the fungal species involved in the observed ash symptoms. Field surveys were conducted in five forest sites where 50 symptomatic branch samples were collected. All samples were inspected and used for fungal isolation. Based on morphology, colony appearance, and DNA sequence data of the internal transcribed spacer region, 125 fungal colonies belonging to five species were isolated and identified. Only a few symptomatic ash samples yielded colonies of H. fraxineus, whereas Botryosphaeriaceae species were isolated with a high frequency, with Diplodia fraxini as the dominant species. A pathogenicity test proved that all isolated species were pathogenic on European ash, causing bark lesions and wood discoloration. All Botryosphaeriaceae species isolated in this study are reported for the first time on European ash in Slovenia.

Canditate metabolites for ash dieback tolerance in Fraxinus excelsior

Ash dieback, a forest epidemic caused by the invasive fungus Hymenoscyphus fraxineus, threatens ash trees throughout Europe. Within Fraxinus excelsior populations, a small proportion of genotypes show a low susceptibility to the pathogen. We compared the metabolomes from a cohort of low-susceptibility ash genotypes with a cohort of high-susceptibility ash genotypes. This revealed two significantly different chemotypes. A total of 64 candidate metabolites associated with reduced or increased susceptibility in the chemical families secoiridoids, coumarins, flavonoids, phenylethanoids, and lignans. Increased levels of two coumarins, fraxetin and esculetin, were strongly associated with reduced susceptibility to ash dieback. Both coumarins inhibited the growth of H. fraxineus in vitro when supplied at physiological concentrations, thereby validating their role as markers for low susceptibility to ash dieback. Similarly, fungal growth inhibition was observed when the methanolic bark extract of low-susceptibility ash genotypes was supplied. Our findings indicate the presence of constitutive chemical defense barriers against ash dieback in ash.

Transcriptional responses in developing lesions of European common ash (Fraxinus excelsior) reveal genes responding to infection by Hymenoscyphus fraxineus

BACKGROUND

With the expanding ash dieback epidemic that has spread across the European continent, an improved functional understanding of the disease development in afflicted hosts is needed. The study investigated whether differences in necrosis extension between common ash (Fraxinus excelsior) trees with different levels of susceptibility to the fungus Hymenoscyphus fraxineus are associated with, and can be explained by, the differences in gene expression patterns. We inoculated seemingly healthy branches of each of two resistant and susceptible ash genotypes with H. fraxineus grown in a common garden.

RESULTS

Ten months after the inoculation, the length of necrosis on the resistant genotypes were shorter than on the susceptible genotypes. RNA sequencing of bark samples collected at the border of necrotic lesions and from healthy tissues distal to the lesion revealed relatively limited differences in gene expression patterns between susceptible and resistant genotypes. At the necrosis front, only 138 transcripts were differentially expressed between the genotype categories while 1082 were differentially expressed in distal, non-symptomatic tissues. Among these differentially expressed genes, several genes in the mevalonate (MVA) and iridoid pathways were found to be co-regulated, possibly indicating increased fluxes through these pathways in response to H. fraxineus. Comparison of transcriptional responses of symptomatic and non-symptomatic ash in a controlled greenhouse experiment revealed a relatively small set of genes that were differentially and concordantly expressed in both studies. This gene-set included the rate-limiting enzyme in the MVA pathway and a number of transcription factors. Furthermore, several of the concordantly expressed candidate genes show significant similarity to genes encoding players in the abscisic acid- or Jasmonate-signalling pathways.

CONCLUSIONS

A set of candidate genes, concordantly expressed between field and greenhouse experiments, was identified. The candidates are associated with hormone signalling and specialized metabolite biosynthesis pathways indicating the involvement of these pathways in the response of the host to infection by H. fraxineus.

Convergent molecular evolution among ash species resistant to the emerald ash borer

Recent studies show that molecular convergence plays an unexpectedly common role in the evolution of convergent phenotypes. We exploited this phenomenon to find candidate loci underlying resistance to the emerald ash borer (EAB, Agrilus planipennis), the United States' most costly invasive forest insect to date, within the pan-genome of ash trees (the genus Fraxinus). We show that EAB-resistant taxa occur within three independent phylogenetic lineages. In genomes from these resistant lineages, we detect 53 genes with evidence of convergent amino acid evolution. Gene-tree reconstruction indicates that, for 48 of these candidates, the convergent amino acids are more likely to have arisen via independent evolution than by another process such as hybridization or incomplete lineage sorting. Seven of the candidate genes have putative roles connected to the phenylpropanoid biosynthesis pathway and 17 relate to herbivore recognition, defence signalling or programmed cell death. Evidence for loss-of-function mutations among these candidates is more frequent in susceptible species than in resistant ones. Our results on evolutionary relationships, variability in resistance, and candidate genes for defence response within the ash genus could inform breeding for EAB resistance, facilitating ecological restoration in areas invaded by this beetle.

Molecular cloning and functional analysis of 4-Coumarate:CoA ligase 4(4CL-like 1)from Fraxinus mandshurica and its role in abiotic stress tolerance and cell wall synthesis

BACKGROUND

Four-Coumarate:CoA ligase gene (4CL) plays multiple important roles in plant growth and development by catalyzing the formation of CoA ester. 4CL belongs to the plant phenylpropane derivative, which is related to the synthesis of flavonoids and lignin and is a key enzyme in the biosynthetic pathway.

RESULTS

In this study, 12 4CL genes of Fraxinus mandschurica were identified and named Fm4CL1-Fm4CL12, respectively. The analysis of the expression pattern of Fm4CL genes indicate that Fm4CL-like 1 gene may play a role in the lignin synthesis pathway. Our study indicate that overexpression of Fm4CL-like 1 increases the lignin content of transgenic tobacco by 39.5% compared to WT, and the S/G ratio of transgenic tobacco increased by 19.7% compared with WT. The xylem cell layer of transgenic line is increased by 40% compared to WT, the xylem cell wall thickness increased by 21.6% compared to the WT. Under mannitol-simulated drought stress, the root length of transgenic tobacco is 64% longer than WT, and the seed germination rate of the transgenic lines is 47% higher than that of WT. In addition, the H2O2 content in the transgenic tobacco was 22% lower than that of WT, while the POD and SOD content was higher than WT by 30 and 24% respectively, which showed Fm4CL-like 1 affect the accumulation of reactive oxygen species (ROS). The MDA content and relative conductivity was 25 and 15% lower than WT, respectively. The water loss rate is 16.7% lower than that of WT. The relative expression levels of stress-related genes NtHAK, NtAPX, NtCAT, NtABF2, and NtZFP were higher than those of WT under stress treatment. The stomatal apertures of OE (Overexpression) were 30% smaller than those of WT, and the photosynthetic rate of OE was 48% higher than that of WT. These results showed that the overexpression line exhibited stronger adaptability to osmotic stress than WT.

CONCLUSIONS

Our results indicate that Fm4CL-like 1 is involved in secondary cell wall development and lignin synthesis. Fm4CL-like 1 play an important role in osmotic stress by affecting cell wall and stomatal development.

Comparative transcriptome analysis revealing the potential mechanism of seed germination stimulated by exogenous gibberellin in Fraxinus hupehensis

BACKGROUND

Fraxinus hupehensis is an endangered tree species that is endemic to in China; the species has very high commercial value because of its intricate shape and potential to improve and protect the environment. Its seeds show very low germination rates in natural conditions. Preliminary experiments indicated that gibberellin (GA₃) effectively stimulated the seed germination of F. hupehensis. However, little is known about the physiological and molecular mechanisms underlying the effect of GA₃ on F. hupehensis seed germination.

RESULTS

We compared dormant seeds (CK group) and germinated seeds after treatment with water (W group) and GA₃ (G group) in terms of seed vigor and several other physiological indicators related to germination, hormone content, and transcriptomics. Results showed that GA₃ treatment increases seed vigor, energy requirements, and trans-Zetain (ZT) and GA₃ contents but decreases sugar and abscisic acid (ABA) contents. A total of 116,932 unigenes were obtained from F. hupehensis transcriptome. RNA-seq analysis identified 31,856, 33,188 and 2056 differentially expressed genes (DEGs) between the W and CK groups, the G and CK groups, and the G and W groups, respectively. Up-regulation of eight selected DEGs of the glycolytic pathway accelerated the oxidative decomposition of sugar to release energy for germination. Up-regulated genes involved in ZT (two genes) and GA₃ (one gene) biosynthesis, ABA degradation pathway (one gene), and ABA signal transduction (two genes) may contribute to seed germination. Two down-regulated genes associated with GA₃ signal transduction were also observed in the G group. GA₃-regulated genes may alter hormone levels to facilitate germination. Candidate transcription factors played important roles in GA₃-promoted F. hupehensis seed germination, and Quantitative Real-time PCR (qRT-PCR) analysis verified the expression patterns of these genes.

CONCLUSION

Exogenous GA₃ increased the germination rate, vigor, and water absorption rate of F. hupehensis seeds. Our results provide novel insights into the transcriptional regulation mechanism of effect of exogenous GA₃ on F. hupehensis seed germination. The transcriptome data generated in this study may be used for further molecular research on this unique species.

The first genetic linkage map for Fraxinus pennsylvanica and syntenic relationships with four related species

The genetic linkage map for green ash (Fraxinus pennsylvanica) contains 1201 DNA markers in 23 linkage groups spanning 2008.87cM. The green ash map shows stronger synteny with coffee than tomato. Green ash (Fraxinus pennsylvanica) is an outcrossing, diploid (2n = 46) hardwood tree species, native to North America. Native ash species in North America are being threatened by the rapid spread of the emerald ash borer (EAB, Agrilus planipennis), an invasive pest from Asia. Green ash, the most widely distributed ash species, is severely affected by EAB infestation, yet few genomic resources for genetic studies and improvement of green ash are available. In this study, a total of 5712 high quality single nucleotide polymorphisms (SNPs) were discovered using a minimum allele frequency of 1% across the entire genome through genotyping-by-sequencing. We also screened hundreds of genomic- and EST-based microsatellite markers (SSRs) from previous de novo assemblies (Staton et al., PLoS ONE 10:e0145031, 2015; Lane et al., BMC Genom 17:702, 2016). A first genetic linkage map of green ash was constructed from 90 individuals in a full-sib family, combining 2719 SNP and 84 SSR segregating markers among the parental maps. The consensus SNP and SSR map contains a total of 1201 markers in 23 linkage groups spanning 2008.87 cM, at an average inter-marker distance of 1.67 cM with a minimum logarithm of odds of 6 and maximum recombination fraction of 0.40. Comparisons of the organization the green ash map with the genomes of asterid species coffee and tomato, and genomes of the rosid species poplar and peach, showed areas of conserved gene order, with overall synteny strongest with coffee.

Decoupling photo- and thermoperiod by projected climate change perturbs bud development, dormancy establishment and vernalization in the model tree Populus

BACKGROUND

The performance and survival of deciduous trees depends on their innate ability to anticipate seasonal change. A key event is the timely production of short photoperiod-induced terminal and axillary buds that are dormant and freezing-tolerant. Some observations suggest that low temperature contributes to terminal bud initiation and dormancy. This is puzzling because low temperatures in the chilling range universally release dormancy. It also raises the broader question if the projected climate instabilities, as well as the northward migration of trees, will affect winter preparations and survival of trees.

RESULTS

To gauge the response capacity of trees, we exposed juvenile hybrid aspens to a 10-h short photoperiod in combination with different day/night temperature regimes: high (24/24 °C), moderate (18/18 °C), moderate-low (18/12 °C) and low (12/12 °C), and analysed bud development, dormancy establishment, and marker gene expression. We found that low temperature during the bud formation period (pre-dormancy) upregulated dormancy-release genes of the gibberellin (GA) pathway, including the key GA biosynthesis genes GA20oxidase and GA3oxidase, the GA-receptor gene GID1, as well as GA-inducible enzymes of the 1,3-β-glucanase family that degrade callose at plasmodesmal Dormancy Sphincter Complexes. Simultaneously, this pre-dormancy low temperature perturbed the expression of flowering pathway genes, including CO, FT, CENL1, AGL14, LFY and AP1. In brief, pre-dormancy low temperature compromised bud development, dormancy establishment, and potentially vernalization. On the other hand, a high pre-dormancy temperature prevented dormancy establishment and resulted in flushing.

CONCLUSIONS

The results show that pre-dormancy low temperature represents a form of chilling that antagonizes dormancy establishment. Combined with available field data, this indicates that natural Populus ecotypes have evolved to avoid the adverse effects of high and low temperatures by initiating and completing dormant buds within an approximate temperature-window of 24-12 °C. Global warming and erratic temperature patterns outside this range can therefore endanger the successful propagation of deciduous perennials.

Demographically idiosyncratic responses to climate change and rapid Pleistocene diversification of the walnut genus Juglans (Juglandaceae) revealed by whole-genome sequences

Whether species demography and diversification are driven primarily by extrinsic environmental changes such as climatic oscillations in the Quaternary or by intrinsic biological interactions like coevolution between antagonists is a matter of active debate. In fact, their relative importance can be assessed by tracking past population fluctuations over considerable time periods. We applied the pairwise sequentially Markovian coalescent approach on the genomes of 11 temperate Juglans species to estimate trajectories of changes in effective population size (N_e ) and used a Bayesian-coalescent based approach that simultaneously considers multiple genomes (G-PhoCS) to estimate divergence times between lineages. N_e curves of all study species converged 1.0 million yr ago, probably reflecting the time when the walnut genus last shared a common ancestor. This estimate was confirmed by the G-PhoCS estimates of divergence times. But all species did not react similarly to the dramatic climatic oscillations following early Pleistocene cooling, so the timing and amplitude of changes in N_e differed among species and even among conspecific lineages. The population histories of temperate walnut species were not driven by extrinsic environmental changes alone, and a key role was probably played by species-specific factors such as coevolutionary interactions with specialized pathogens.

Gene flow of common ash (Fraxinus excelsior L.) in a fragmented landscape

Gene flow dynamics of common ash (Fraxinus excelsior L.) is affected by several human activities in Central Europe, including habitat fragmentation, agroforestry expansion, controlled and uncontrolled transfer of reproductive material, and a recently introduced emerging infectious disease, ash dieback, caused by Hymenoscyphus fraxineus. Habitat fragmentation may alter genetic connectivity and effective population size, leading to loss of genetic diversity and increased inbreeding in ash populations. Gene flow from cultivated trees in landscapes close to their native counterparts may also influence the adaptability of future generations. The devastating effects of ash dieback have already been observed in both natural and managed populations in continental Europe. However, potential long-term effects of genetic bottlenecks depend on gene flow across fragmented landscapes. For this reason, we studied the genetic connectivity of ash trees in an isolated forest patch of a fragmented landscape in Rösenbeck, Germany. We applied two approaches to parentage analysis to estimate gene flow patterns at the study site. We specifically investigated the presence of background pollination at the landscape level and the degree of genetic isolation between native and cultivated trees. Local meteorological data was utilized to understand the effect of wind on the pollen and seed dispersal patterns. Gender information of the adult trees was considered for calculating the dispersal distances. We found that the majority of the studied seeds (55-64%) and seedlings (75-98%) in the forest patch were fathered and mothered by the trees within the same patch. However, we determined a considerable amount of pollen flow (26-45%) from outside of the study site, representing background pollination at the landscape level. Limited pollen flow was observed from neighbouring cultivated trees (2%). Both pollen and seeds were dispersed in all directions in accordance with the local wind directions. Whereas there was no positive correlation between pollen dispersal distance and wind speed, the correlation between seed dispersal distance and wind speed was significant (0.71, p < 0.001), indicating that strong wind favours long-distance dispersal of ash seeds. Finally, we discussed the implications of establishing gene conservation stands and the use of enrichment planting in the face of ash dieback.

Rising Out of the Ashes: Additive Genetic Variation for Crown and Collar Resistance to Hymenoscyphus fraxineus in Fraxinus excelsior

Since the early 1990s, ash dieback due to the invasive ascomycete Hymenoscyphus fraxineus is threatening Fraxinus excelsior in most of its natural range. Previous studies reported significant levels of genetic variability in susceptibility in F. excelsior either in field or inoculation experiments. The present study was based on a field experiment planted in 1995, 15 years before onset of the disease. Crown and collar status were monitored on 777 trees from 23 open-pollinated progenies originating from three French provenances. Health status was modeled using a Bayesian approach where spatiotemporal effects were explicitly taken into account. Moderate narrow-sense heritability was found for crown dieback (h² = 0.42). This study is first to show that resistance at the collar level is also heritable (h² = 0.49 for collar lesions prevalence and h² = 0.42 for their severity) and that there is significant genetic correlation (r = 0.40) between the severities of crown and collar symptoms. There was no evidence for differences between provenances. Family effects were detected, but computing individual breeding values showed that most of the genetic variation lies within families. In agreement with previous reports, early flushing correlates with healthier crown. Implications of these results in disease management and breeding are discussed.

Improvement and transcriptome analysis of root architecture by overexpression of Fraxinus pennsylvanica DREB2A transcription factor in Robinia pseudoacacia L. 'Idaho'

Transcription factors play a key role to enable plants to cope with abiotic stresses. DREB2 regulates the expression of several stress-inducible genes and constitutes major hubs in the water stress signalling webs. We cloned and characterized a novel gene encoding the FpDREB2A transcription factor from Fraxinus pennsylvanica, and a yeast activity assay confirmed its DRE binding and transcription activation. Overexpression of FpDREB2A in R. pseudoacacia showed enhanced resistance to drought stress. The transgenic plant survival rate was significantly higher than that of WT in soil drying and re-watering treatments. Transgenic lines showed a dramatic change in root architecture, and horizontal and vertical roots were found in transgenic plants compared to WT. The vertical roots penetrated in the field soil to more than 60 cm deep, while horizontal roots expanded within the top 20-30 cm of the soil. A physiological test demonstrated that chlorophyll contents were more gradually reduced and that soluble sugars and proline levels elevated more sharply but malondialdehyde level stayed the same (P < 0.05). Plant hormone levels of abscisic acid and IAA were higher than that of WT, while gibberellins and zeatin riboside were found to be lower. The root transcriptomes were sequenced and annotated into 2011 differential expression genes (DEGs). The DEGs were categorized in 149 pathways and were found to be involved in plant hormone signalling, transcription factors, stimulus responses, phenylalanine, carbohydrate and other metabolic pathways. The modified pathways in plant hormone signalling are thought to be the main cause of greater horizontal and vertical root development, in particular.

Evidence for the long-term maintenance of a rare self-incompatibility system in Oleaceae

A rare homomorphic diallelic self-incompatibility (DSI) system discovered in Phillyrea angustifolia (family Oleaceae, subtribe Oleinae) can promote the transition from hermaphroditism to androdioecy. If widespread and stable in Oleaceae, DSI may explain the exceptionally high rate of androdioecious species reported in this plant family. Here, we set out to determine whether DSI occurs in another Oleaceae lineage. We tested for DSI in subtribe Fraxininae, a lineage that diverged from subtribe Oleinae c. 40 million yr ago. We explored the compatibility relationships in Fraxinus ornus using 81 hermaphrodites and 25 males from one natural stand and two naturalized populations using intra- and interspecific stigma tests performed on F. ornus and P. angustifolia testers. We uncovered a DSI system with hermaphrodites belonging to one of two self-incompatibility (SI) groups and males compatible with both groups, making for a truly androdioecious reproductive system. The two human-founded populations contained only one of the two SI groups. Our results provide evidence for the evolutionary persistence of DSI. We discuss how its stability over time may have affected transitions to other sexual systems, such as dioecy.

Molecular markers for tolerance of European ash (Fraxinus excelsior) to dieback disease identified using Associative Transcriptomics

Tree disease epidemics are a global problem, impacting food security, biodiversity and national economies. The potential for conservation and breeding in trees is hampered by complex genomes and long lifecycles, with most species lacking genomic resources. The European Ash tree Fraxinus excelsior is being devastated by the fungal pathogen Hymenoscyphus fraxineus, which causes ash dieback disease. Taking this system as an example and utilizing Associative Transcriptomics for the first time in a plant pathology study, we discovered gene sequence and gene expression variants across a genetic diversity panel scored for disease symptoms and identified markers strongly associated with canopy damage in infected trees. Using these markers we predicted phenotypes in a test panel of unrelated trees, successfully identifying individuals with a low level of susceptibility to the disease. Co-expression analysis suggested that pre-priming of defence responses may underlie reduced susceptibility to ash dieback.

Characterization, localization, and seasonal changes of the sucrose transporter FeSUT1 in the phloem of Fraxinus excelsior

Trees are generally assumed to be symplastic phloem loaders. A typical feature for most wooden species is an open minor vein structure with symplastic connections between mesophyll cells and phloem cells, which allow sucrose to move cell-to-cell through the plasmodesmata into the phloem. Fraxinus excelsior (Oleaceae) also translocates raffinose family oligosaccharides in addition to sucrose. Sucrose concentration was recently shown to be higher in the phloem sap than in the mesophyll cells. This suggests the involvement of apoplastic steps and the activity of sucrose transporters in addition to symplastic phloem-loading processes. In this study, the sucrose transporter FeSUT1 from F. excelsior was analysed. Heterologous expression in baker's yeast showed that FeSUT1 mediates the uptake of sucrose. Immunohistochemical analyses revealed that FeSUT1 was exclusively located in phloem cells of minor veins and in the transport phloem of F. excelsior. Further characterization identified these cells as sieve elements and possibly ordinary companion cells but not as intermediary cells. The localization and expression pattern point towards functions of FeSUT1 in phloem loading of sucrose as well as in sucrose retrieval. FeSUT1 is most likely responsible for the observed sucrose gradient between mesophyll and phloem. The elevated expression level of FeSUT1 indicated an increased apoplastic carbon export activity from the leaves during spring and late autumn. It is hypothesized that the importance of apoplastic loading is high under low-sucrose conditions and that the availability of two different phloem-loading mechanisms confers advantages for temperate woody species like F. excelsior.

Agrobacterium-mediated genetic transformation and plant regeneration of the hardwood tree species Fraxinus profunda

This transformation and regeneration protocol provides an integral framework for the genetic improvement of Fraxinus profunda (pumpkin ash) for future development of plants resistant to the emerald ash borer. Using mature hypocotyls as the initial explants, an Agrobacterium tumefaciens-mediated genetic transformation system was successfully developed for pumpkin ash (Fraxinus profunda). This transformation protocol is an invaluable tool to combat the highly aggressive, non-native emerald ash borer (EAB), which has the potential to eliminate native Fraxinus spp. from the natural landscape. Hypocotyls were successfully transformed with Agrobacterium strain EHA105 harboring the pq35GR vector, containing an enhanced green fluorescent protein (EGFP) as well as a fusion gene between neomycin phosphotransferase (nptII) and gusA. Hypocotyls were cultured for 7 days on Murashige and Skoog (MS) medium with 22.2 μM 6-benzyladenine (BA), 4.5 μM thidiazuron (TDZ), 50 mg L(-1) adenine hemisulfate (AS), and 10 % coconut water (CW) prior to transformation. Hypocotyls were transformed using 90 s sonication plus 10 min vacuum infiltration after Agrobacterium was exposed to 100 μM acetosyringone for 1 h. Adventitious shoots were regenerated on MS medium with 22.2 μM BA, 4.5 μM TDZ, 50 mg L(-1) AS, 10 % CW, 400 mg L(-1) timentin, and 20 mg L(-1) kanamycin. Timentin at 400 and 20 mg L(-1) kanamycin were most effective at controlling Agrobacterium growth and selecting for transformed cells, respectively. The presence of nptII, GUS (β-glucuronidase), and EGFP in transformed plants was confirmed using polymerase chain reaction (PCR), while the expression of EGFP was also confirmed through fluorescent microscopy and reverse transcription-PCR. This transformation protocol provides an integral foundation for future genetic modifications of F. profunda to provide resistance to EAB.

Genotypes of Fraxinus excelsior with different susceptibility to the ash dieback pathogen Hymenoscyphus pseudoalbidus and their response to the phytotoxin viridiol - a metabolomic and microscopic study

Eight European ash (Fraxinus excelsior) genotypes with different known susceptibility to Hymenoscyphus pseudoalbidus were tested against the phytotoxin viridiol and their response described at the microscopic and metabolomic level. All ash genotypes were sensitive to the toxin and necrosis was detectable after 24h. Among the three viridiol concentrations used in the experiment, the lowest concentration (14.7μM) yielded markedly lower mean damage scores compared to those resulting from seedlings tested at higher dosages. The highest damage scores were associated with the susceptible ash clones S-101, S-106 and S-125, but also with resistant clone R-104. Three resistant clones (R-131, R-121, and R-118) had lower mean damage scores compared to susceptible clones. Wilting of leaves was more common 48h after treatment and more pronounced on seedlings with high damage scores. The resulting lesions generally lacked browning of tissue and displayed only surface disruption of cells in direct contact with the toxin. A delay in symptom development was evident on all five resistant clones tested with the two higher concentrations of viridiol. LC-HRMS and MS/MS analyses of ash seedling extracts suggest several secoiridoid compounds as well as compounds related to abscisic acid (ABA) to be produced in response to viridiol. ABA-cysteine and xanthoxin were found at significantly higher concentrations in susceptible clones compared to resistant clones after treatment with viridiol, suggesting a primary role of ABA in response to stress. The results observed in this study suggest that genetic resistance to H. pseudoalbidus among ash genotypes may be explained, in part, by the varied response to phytotoxins produced by the fungus.

The phylogeography of Eurasian Fraxinus species reveals ancient transcontinental reticulation

To investigate the biogeographical history of ashes species of the Eurasian section Fraxinus and to test the hypothesis of ancient reticulations, we sequenced nuclear DNA (nETS and nITS, 1075 bp) for 533 samples and scored AFLP for 63 samples of Eurasian ashes within the section Fraxinus. The nITS phylogeny retrieved the classical view of the evolution of the section, whereas nETS phylogeny indicated an unexpected separation of F. angustifolia in two paraphyletic groups, respectively found in southeastern Europe and in the other parts of the Mediterranean basin. In the nETS phylogeny, the former group was closely related to F. excelsior, whereas the later was closely related to F. mandshurica, a species which is restricted nowadays to northeastern Asia. This topological incongruence between the two loci indicated the occurrence of an ancient reticulation between European and Asian ash species. Several other ancient reticulation events between the two European species and the other species of the section were supported by the posterior predictive checking method. Some of these reticulation events would have occurred during the Miocene, when climatic variations may have lead these species to expand their distribution range and come into contact. The recurrent reticulations observed among Eurasian ash species indicate that they should be considered as conspecific taxa, with subspecific status for some groups. Altogether, the results of the present study provide a rare documented evidence for the occurrence of multiple ancient reticulations within a group of temperate tree taxa with modern disjunct distributions in Eurasia. These ancient reticulation events indicate that the speciation process is slow in ashes, necessitating long periods of geographical isolation. The implications for speciation processes in temperate trees with similar life history and reproductive biology are discussed.

Apoplastic and symplastic phloem loading in Quercus robur and Fraxinus excelsior

Whereas most of the research on phloem loading is performed on herbaceous plants, less is known about phloem loading strategies in trees. In this study, the phloem loading mechanisms of Quercus robur and Fraxinus excelsior were analysed. The following features were examined: the minor vein structure, the sugar concentrations in phloem sap by the laser-aphid-stylet technique, the distribution of photoassimilates in the mesophyll cells by non-aqueous fractionation, gradients of sugar concentrations and osmotic pressure, and the expression of sucrose transporters. The minor vein configurations of Q. robur and F. excelsior belong to the open type. Quercus robur contained companion cells in the minor veins whereas F. excelsior showed intermediary cells in addition to ordinary companion cells. The main carbon transport form in Q. robur was sucrose (~1M). In F. excelsior high amounts of raffinose and stachyose were also transported. However, in both tree species, the osmolality of phloem sap was higher than the osmolality of the mesophyll cells. The concentration gradients between phloem sap and the cytoplasm of mesophyll cells for sucrose were 16-fold and 14-fold for Q. robur and F. excelsior, respectively. Independent of the type of translocated sugars, sucrose transporter cDNAs were cloned from both species. The results indicate that phloem loading of sucrose and other metabolites must involve active loading steps in both tree species. Quercus robur seems to be an apoplastic phloem loader while F. excelsior shows indications of being a symplastic or mixed symplastic-apoplastic phloem loader.

Emerald ash borer invasion of North America: history, biology, ecology, impacts, and management

Since its accidental introduction from Asia, emerald ash borer (EAB), Agrilus planipennis Fairmaire (Coleoptera: Buprestidae), has killed millions of ash trees in North America. As it continues to spread, it could functionally extirpate ash with devastating economic and ecological impacts. Little was known about EAB when it was first discovered in North America in 2002, but substantial advances in understanding of EAB biology, ecology, and management have occurred since. Ash species indigenous to China are generally resistant to EAB and may eventually provide resistance genes for introgression into North American species. EAB is characterized by stratified dispersal resulting from natural and human-assisted spread, and substantial effort has been devoted to the development of survey methods. Early eradication efforts were abandoned largely because of the difficulty of detecting and delineating infestations. Current management is focused on biological control, insecticide protection of high-value trees, and integrated efforts to slow ash mortality.

Environmental heterogeneity explains the genetic structure of Continental and Mediterranean populations of Fraxinus angustifolia Vahl

Tree species with wide distributions often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in a wind-pollinated Mediterranean tree species, Fraxinus angustifolia Vahl, within a recognised glacial refugium in Croatia. We sampled 11 populations from environmentally divergent habitats within the Continental and Mediterranean biogeographical regions. We combined genetic data analyses based on nuclear microsatellite loci, multivariate statistics on environmental data and ecological niche modelling (ENM). We identified a geographic structure with a high genetic diversity and low differentiation in the Continental region, which contrasted with the significantly lower genetic diversity and higher population divergence in the Mediterranean region. The positive and significant correlation between environmental and genetic distances after controlling for geographic distance suggests an important influence of ecological divergence of the sites in shaping genetic variation. The ENM provided support for niche differentiation between the populations from the Continental and Mediterranean regions, suggesting that contemporary populations may represent two divergent ecotypes. Ecotype differentiation was also supported by multivariate environmental and genetic distance analyses. Our results suggest that despite extensive gene flow in continental areas, long-term stability of heterogeneous environments have likely promoted genetic divergence of ashes in this region and can explain the present-day genetic variation patterns of these ancient populations.

Deciduous trees and the application of universal DNA barcodes: a case study on the circumpolar Fraxinus

The utility of DNA barcoding for identifying representative specimens of the circumpolar tree genus Fraxinus (56 species) was investigated. We examined the genetic variability of several loci suggested in chloroplast DNA barcode protocols such as matK, rpoB, rpoC1 and trnH-psbA in a large worldwide sample of Fraxinus species. The chloroplast intergenic spacer rpl32-trnL was further assessed in search for a potentially variable and useful locus. The results of the study suggest that the proposed cpDNA loci, alone or in combination, cannot fully discriminate among species because of the generally low rates of substitution in the chloroplast genome of Fraxinus. The intergenic spacer trnH-psbA was the best performing locus, but genetic distance-based discrimination was moderately successful and only resulted in the separation of the samples at the subgenus level. Use of the BLAST approach was better than the neighbor-joining tree reconstruction method with pairwise Kimura's two-parameter rates of substitution, but allowed for the correct identification of only less than half of the species sampled. Such rates are substantially lower than the success rate required for a standardised barcoding approach. Consequently, the current cpDNA barcodes are inadequate to fully discriminate Fraxinus species. Given that a low rate of substitution is common among the plastid genomes of trees, the use of the plant cpDNA "universal" barcode may not be suitable for the safe identification of tree species below a generic or sectional level. Supplementary barcoding loci of the nuclear genome and alternative solutions are proposed and discussed.

Interspecific proteomic comparisons reveal ash phloem genes potentially involved in constitutive resistance to the emerald ash borer

The emerald ash borer (Agrilus planipennis) is an invasive wood-boring beetle that has killed millions of ash trees since its accidental introduction to North America. All North American ash species (Fraxinus spp.) that emerald ash borer has encountered so far are susceptible, while an Asian species, Manchurian ash (F. mandshurica), which shares an evolutionary history with emerald ash borer, is resistant. Phylogenetic evidence places North American black ash (F. nigra) and Manchurian ash in the same clade and section, yet black ash is highly susceptible to the emerald ash borer. This contrast provides an opportunity to compare the genetic traits of the two species and identify those with a potential role in defense/resistance. We used Difference Gel Electrophoresis (DIGE) to compare the phloem proteomes of resistant Manchurian to susceptible black, green, and white ash. Differentially expressed proteins associated with the resistant Manchurian ash when compared to the susceptible ash species were identified using nano-LC-MS/MS and putative identities assigned. Proteomic differences were strongly associated with the phylogenetic relationships among the four species. Proteins identified in Manchurian ash potentially associated with its resistance to emerald ash borer include a PR-10 protein, an aspartic protease, a phenylcoumaran benzylic ether reductase (PCBER), and a thylakoid-bound ascorbate peroxidase. Discovery of resistance-related proteins in Asian species will inform approaches in which resistance genes can be introgressed into North American ash species. The generation of resistant North American ash genotypes can be used in forest ecosystem restoration and urban plantings following the wake of the emerald ash borer invasion.

Presence of natural genetic resistance in Fraxinus excelsior (Oleraceae) to Chalara fraxinea (Ascomycota): an emerging infectious disease

Fraxinus excelsior, common ash native to Europe, is threatened by a recently identified pathogenic fungus Chalara fraxinea, which causes extensive damage on ash trees across Europe. In Denmark, most stands are severely affected leaving many trees with dead crowns. However, single trees show notably fewer symptoms. In this study, the impact of the emerging infectious disease on native Danish ash trees is assessed by estimating presence of inherent resistance in natural populations. Disease symptoms were assessed from 2007 to 2009 at two different sites with grafted ramets of 39 selected clones representing native F. excelsior trees. A strong genetic variation in susceptibility to C. fraxinea infections was observed. No genetic or geographic structure can explain the differences, but strong genetic correlations to leaf senescence were observed. The results suggest that a small fraction of trees in the Danish population of ash possess substantial resistance against the damage. Though this fraction is probably too low to avoid population collapse in most natural or managed ash forests, the observed presence of putative resistance against the emerging infectious disease in natural stands is likely to be of evolutionary importance. This provides prospects of future maintenance of the species through natural or artificial selection in favour of remaining healthy individuals.

Transcriptomic signatures of ash (Fraxinus spp.) phloem

BACKGROUND

Ash (Fraxinus spp.) is a dominant tree species throughout urban and forested landscapes of North America (NA). The rapid invasion of NA by emerald ash borer (Agrilus planipennis), a wood-boring beetle endemic to Eastern Asia, has resulted in the death of millions of ash trees and threatens billions more. Larvae feed primarily on phloem tissue, which girdles and kills the tree. While NA ash species including black (F. nigra), green (F. pennsylvannica) and white (F. americana) are highly susceptible, the Asian species Manchurian ash (F. mandshurica) is resistant to A. planipennis perhaps due to their co-evolutionary history. Little is known about the molecular genetics of ash. Hence, we undertook a functional genomics approach to identify the repertoire of genes expressed in ash phloem.

METHODOLOGY AND PRINCIPAL FINDINGS

Using 454 pyrosequencing we obtained 58,673 high quality ash sequences from pooled phloem samples of green, white, black, blue and Manchurian ash. Intriguingly, 45% of the deduced proteins were not significantly similar to any sequences in the GenBank non-redundant database. KEGG analysis of the ash sequences revealed a high occurrence of defense related genes. Expression analysis of early regulators potentially involved in plant defense (i.e. transcription factors, calcium dependent protein kinases and a lipoxygenase 3) revealed higher mRNA levels in resistant ash compared to susceptible ash species. Lastly, we predicted a total of 1,272 single nucleotide polymorphisms and 980 microsatellite loci, among which seven microsatellite loci showed polymorphism between different ash species.

CONCLUSIONS AND SIGNIFICANCE

The current transcriptomic data provide an invaluable resource for understanding the genetic make-up of ash phloem, the target tissue of A. planipennis. These data along with future functional studies could lead to the identification/characterization of defense genes involved in resistance of ash to A. planipennis, and in future ash breeding programs for marker development.

Isolation and characterisation of bacterial colonies from seeds and in vitro cultures of Fraxinus spp. from Italian sites

Culturable bacteria were isolated from seeds, embryos and contaminated in vitro cultures of ash (Fraxinus excelsior L., F. ornus L. and F. angustifolia L.) and were identified using morphological and molecular analyses. Fourteen morphologically distinct isolates were recovered from seeds of Fraxinus spp. 16S rDNA sequencing categorised these isolates into ten separate genera. Three strains isolated from contaminated in vitro cultures, Pantoea agglomerans, Staphylococcus succinus and Aerococcus viridans, were used for comparative analysis with isolates from seeds. Antibiotic sensitivity testing of the isolated contaminants, including phytotoxicity of antibiotics on in vitro cultures of ash, was also investigated. Phytotoxic effects on explants immersed in ampicillin or cultured on medium containing ampicillin were negligible, however tetracycline, either alone or in combination with other antibiotics, had phytotoxic effects. We conclude that ampicillin is a suitable antibiotic to limit the growth of contaminating bacteria during the in vitro culture of ash.

Cryopreservation of dormant buds from diverse Fraxinus species

Ash (Fraxinus) is an economically important tree genus in the landscape industry, as well as a key component of North American forests, especially in the North Central United States and adjacent regions in Canada. In recent years, the Emerald Ash Borer beetle (Agrilus planipennis) has significantly threatened the survival of native North American Fraxinus species. A dormant-bud cryopreservation technique has been developed as a method to conserve specific clones of ash. Dormant buds of three ash species were successfully cryopreserved when desiccated on their stem sections to 30 percent moisture content (w/v) and then cooled at rates of either -1 degree C per h or -5 degree C per day to either -30 or -35 degree C before immersion in liquid nitrogen vapor (LNV). Stem sections were removed from LNV, warmed, and rehydrated, and their buds grafted onto rootstocks to evaluate survival. Recovery percentages ranged from 34 to 100 percent after LNV exposure and were dependent upon accession and cooling rate. The cryopreservation methods proposed herein can complement seed-collection efforts aimed at conserving diversity, supplementing ex situ genebank and botanic-garden collections.

Agrobacterium-mediated transformation of Fraxinus pennsylvanica hypocotyls and plant regeneration

A genetic transformation protocol for green ash (Fraxinus pennsylvanica) hypocotyl explants was developed. Green ash hypocotyls were transformed using Agrobacterium tumefaciens strain EHA105 harboring binary vector pq35GR containing the neomycin phosphotransferase (nptII) and beta-glucuronidase (GUS) fusion gene, and an enhanced green fluorescent protein gene. Pre-cultured hypocotyl explants were transformed in the presence of 100 microM acetosyringone using 90 s sonication plus 10 min vacuum-infiltration. Kanamycin at 20 mg l(-1) was used for selecting transformed cells. Adventitious shoots regenerated on Murashige and Skoog medium supplemented with 13.3 microM 6-benzylaminopurine, 4.5 microM thidiazuron, 50 mg l(-1) adenine sulfate, and 10% coconut water. GUS- and polymerase chain reaction (PCR)-positive shoots from the cut ends of hypocotyls were produced via an intermediate callus stage. Presence of the GUS and nptII genes in GUS-positive shoots were confirmed by PCR and copy number of the nptII gene in PCR-positive shoots was determined by Southern blotting. Three transgenic plantlets were acclimatized to the greenhouse. This transformation and regeneration system using hypocotyls provides a foundation for Agrobacterium-mediated transformation of green ash. Studies are underway using a construct containing the Cry8Da protein of Bacillus thuringiensis for genetic transformation of green ash.

Nuclear DNA microsatellites reveal genetic variation but a lack of phylogeographical structure in an endangered species, Fraxinus mandshurica, across North-east China

BACKGROUND AND AIMS

The widely accepted paradigm that the modern genetic structure of plant species in the northern hemisphere has been largely determined by recolonization from refugia after the last glacial maximum fails to explain the presence of cold-tolerant species at intermediate latitudes. Another generally accepted paradigm is that mountain ridges act as important barriers causing genetic isolation of species, but this too has been challenged in recent studies. The aims of the work reported here were to determine the genetic diversity and distribution patterns of extant natural populations of an endangered cool temperate species, Faxinus mandshurica, and to examine whether these two paradigms are appropriate when applied to this species over a wide geographical scale.

METHODS

1435 adult individuals were sampled from 30 natural populations across the main and central range of the species, covering major mountain ranges across North-east China (NEC). Genetic variation was estimated based on nine polymorphic nuclear microsatellite loci. Phylogeographical analyses were employed using various approaches, including Bayesian clustering, spatial analysis of molecular variance, Monmonier's algorithm, neighbor-joining trees, principal co-ordinate analysis and isolation by distance.

KEY RESULTS

Genetic diversity within populations was relatively high, and no significant recent bottlenecks were detected in any of the populations. A significant negative correlation between intra-population genetic diversity and latitude was identified. In contrast, genetic differentiation among all the populations examined was extremely low and no clear geographic genetic structure was identified, with the exception of one distinct population.

CONCLUSIONS

The modern genetic structure in this species can be explained by extensive gene flow, an absence of mountains acting as barriers, and the presence of a wide refuge across NEC rather than multiple small refugia. Intra-population genetic variation along latitudes is probably associated with the systematically northward shifts of forest biomes in eastern China during the mid-Holocene. To determine important genetic patterns and identify resources for conservation, however, it will be necessary to examine differentially inherited genetic markers exposed to selection pressures (e.g. chloroplast DNA) and to investigate different generations.

Temporal cline in a hybrid zone population between Fraxinus excelsior L. and Fraxinus angustifolia Vahl

The two closely related ash species Fraxinus excelsior L. (common ash) and Fraxinus angustifolia Vahl (narrow-leaved ash) have a broad contact zone in France where they hybridize. However, little is known about the local structure of hybrid zone populations and the isolation mechanisms. We assessed the potential effect of floral phenology on the structure of a riparian ash hybrid zone population in central France. The distribution of flowering times was unimodal and lay between the flowering periods of the two species. Using microsatellite markers, we detected isolation by time, which has possibly originated from assortative mating. Multivariate analyses indicated that morphological variation is not distributed at random with respect to flowering times. Spatial autocorrelation analyses showed that temporal and spatial patterns were tightly linked. Interestingly, despite the fact that the population shows isolation by time, neighbourhood size and historical dispersal variance (sigma = 63 m) are similar to those detected in pure stands of F. excelsior where individuals flower rather synchronously and hermaphrodites are not the most frequent sexual type. Trees flowering at intermediate dates, which comprised the majority of the population, produced on average more flowers and fruits. We detected no significant differences in floral parasite infections relative to reproductive timing, although there was a tendency for late flowering trees to suffer from more gall attack. We discuss the impact of temporal variation in fitness traits and their possible role in the maintenance of the hybrid zone.

Assortative mating and differential male mating success in an ash hybrid zone population

BACKGROUND

The structure and evolution of hybrid zones depend mainly on the relative importance of dispersal and local adaptation, and on the strength of assortative mating. Here, we study the influence of dispersal, temporal isolation, variability in phenotypic traits and parasite attacks on the male mating success of two parental species and hybrids by real-time pollen flow analysis. We focus on a hybrid zone population between the two closely related ash species Fraxinus excelsior L. (common ash) and F. angustifolia Vahl (narrow-leaved ash), which is composed of individuals of the two species and several hybrid types. This population is structured by flowering time: the F. excelsior individuals flower later than the F. angustifolia individuals, and the hybrid types flower in-between. Hybrids are scattered throughout the population, suggesting favorable conditions for their local adaptation. We estimate jointly the best-fitting dispersal kernel, the differences in male fecundity due to variation in phenotypic traits and level of parasite attack, and the strength of assortative mating due to differences in flowering phenology. In addition, we assess the effect of accounting for genotyping error on these estimations.

RESULTS

We detected a very high pollen immigration rate and a fat-tailed dispersal kernel, counter-balanced by slight phenological assortative mating and short-distance pollen dispersal. Early intermediate flowering hybrids, which had the highest male mating success, showed optimal sex allocation and increased selfing rates. We detected asymmetry of gene flow, with early flowering trees participating more as pollen donors than late flowering trees.

CONCLUSION

This study provides striking evidence that long-distance gene flow alone is not sufficient to counter-act the effects of assortative mating and selfing. Phenological assortative mating and short-distance dispersal can create temporal and spatial structuring that appears to maintain this hybrid population. The asymmetry of gene flow, with higher fertility and increased selfing, can potentially confer a selective advantage to early flowering hybrids in the zone. In the event of climate change, hybridization may provide a means for F. angustifolia to further extend its range at the expense of F. excelsior.

Differential patterns of morphological and molecular hybridization between Fraxinus excelsior L. and Fraxinus angustifolia Vahl (Oleaceae) in eastern and western France

We examined large-scale patterns of morphology, genetic structure and ecological correlates of Fraxinus excelsior and the closely related species Fraxinus angustifolia in France, in order to determine the degree of hybridization between them. We sampled 24 populations in two putative hybrid zones (Loire and Saône), and five control populations of each species. We measured foliar characteristics of adult trees and used five nuclear microsatellites as molecular markers. Canonical discriminant analysis indicated that the two species differ in morphology, but that intermediate types are common in the Loire region but less frequent in the Saône region. Bayesian population assignment identified one F. angustifolia and two F. excelsior gene pools. Most Loire individuals clustered genetically with the F. angustifolia gene pool. In contrast, the Saône region presented individuals belonging mostly to F. excelsior pools, although the F. angustifolia type was frequent in certain populations. The lowest F(ST) values were found between the Loire and F. angustifolia controls that also exhibited no significant isolation by distance. The proportion of the F. angustifolia gene pool in each locality was negatively correlated with winter temperatures, suggesting that a cold climate may be limiting. Hybridization is probably favoured by the intermediate climatic conditions in the Loire region that allow both species to occur, but is somewhat hampered by the harsher winters in the Saône area where morphological introgression has apparently not yet occurred.

Ovule discounting in an outcrossing, cryptically dioecious tree

Ovule discounting denotes the reduction in the number of ovules available for cross-fertilization due to the interference of inferior pollen. Traditionally, ovule discounting has been discussed solely from the perspective of compromised outcrossing opportunities as a result of selfing, but the principle is more general. Here, we extend its applicability beyond the simple contrast between selfing and outcrossing by showing that, in the cryptically dioecious tree species Fraxinus ornus, ovule discounting through frequent outcrossing with inferior fathers also constitutes a substantial cost of mating. In F. ornus, hermaphrodites produce pollen capable of siring offspring, but these offspring are less viable than those sired by males and are inferred to produce few, if any, surviving progeny. In this paper, we used microsatellite markers to analyze the mating system and paternity in a wild population of F. ornus. We found that the effective number of sires per mother was low (N(ep) = 2.93 to 4.95), and that paternity was correlated among progeny sampled from the same mother, but not among progeny sampled from neighboring mothers. Despite the existence of a local spatial genetic structure (up to 30 m), we found no evidence of biparental inbreeding. There was negligible selfing by hermaphrodites, but they sired approximately one fourth of the seeds produced by other hermaphrodites. Given that these progeny are not inferred to reach reproductive maturity, this constitutes a substantial cost of ovule discounting in the broad sense. We discuss the possible reasons for why hermaphrodites invest resources into inferior pollen.

Fat-tailed gene flow in the dioecious canopy tree species Fraxinus mandshurica var. japonica revealed by microsatellites

Pollen flow, seed dispersal and individual reproductive success can be simultaneously estimated from the genotypes of adults and offspring using stochastic models. Using four polymorphic microsatellite loci, gene flow of the wind-pollinated and wind-seed-dispersed dioecious tree species, Fraxinus mandshurica var. japonica, was quantified in a riparian forest, in northern Japan. In a 10.5-ha plot, 74 female adults, 76 male adults and 292 current-year seedlings were mapped and genotyped, together with 200 seeds. To estimate dispersal kernels of pollen and seeds, we applied normal, exponential power, Weibull, bivariate t-distribution kernels, and two-component models consisting of two normal distribution functions, one with a small and one with a large variance. A two-component pollen flow model with a small contribution (26.1%) from short-distance dispersal (sigma = 7.2 m), and the rest from long-distance flow (sigma = 209.9 m), was chosen for the best-fitting model. The average distance that integrated pollen flows inside and outside the study plot was estimated to be 196.8 m. Tree size and flowering intensity affected reproduction, and there appeared to be critical values that distinguished reproductively successful and unsuccessful adults. In contrast, the gene flow model that estimated both pollen and seed dispersal from established seedlings resulted in extensive seed dispersal, and the expected spatial genetic structures did not satisfactorily fit with the observations, even for the selected model. Our results advanced small-scale individual-based parentage analysis for quantifying fat-tailed gene flow in wind-mediated species, but also clarified its limitations and suggested future possibilities for gene flow studies.

Chloroplast DNA phylogeography of European ashes,Fraxinussp. (Oleaceae): roles of hybridization and life history traits

We investigated range-wide phylogeographic variation in three European ash species (Fraxinus sp., Oleaceae). Chloroplast DNA (cpDNA) microsatellites were typed in the thermophilous Fraxinus angustifolia and Fraxinus ornus and the observed haplotypes and the geographic distribution of diversity were compared to cpDNA data previously obtained in the more cold-tolerant Fraxinus excelsior. We found wide-ranging haplotype sharing between the phylogenetically close F. angustifolia and F. excelsior, suggesting hybridization (i) in common glacial refuges in the Iberian Peninsula, northern Italy, the eastern and/or Dinaric Alps and the Balkan Peninsula, and/or (ii) during postglacial recolonization. The data allowed us to propose additional glacial refuges for F. angustifolia in southern Italy and in Turkey, and populations from the latter region were particularly polymorphic. There was evidence for refuge areas in Italy, the Balkan Peninsula and Turkey for F. ornus, which did not share any single chloroplast haplotype with the other species. In both F. angustifolia and F. ornus, cpDNA diversity (h(S) = 0.027 and h(S) = 0.009, respectively) was lower and fixation levels (G(ST) = 0.964 and G(ST) = 0.983, respectively) higher than in sympatric F. excelsior (h(S) = 0.096, G(ST) = 0.870). These diversity patterns could be due to temperature tolerance or the demographic history.

Effective seed dispersal across a fragmented landscape

The role of seed dispersal in maintaining genetic connectivity among forest fragments has largely been ignored because gene flow by pollen is expected to predominate. By using genealogical reconstruction, we investigated gene flow after establishment of seeds in a wind-pollinated, wind-dispersed tree. Our data show that seed dispersal is the main vector of gene flow among remnants and that long-distance dispersal is common across a chronically fragmented landscape. The relative importance of seed-mediated gene flow may have been underemphasized in other fragmented systems, and diagnosing the response of forest trees to current anthropogenic disturbances requires the assessment of phenomena after establishment.

Characterization and primer development for amplification of chloroplast microsatellite regions of Fraxinus excelsior

This study reports the characterization of chloroplast DNA (cpDNA) variation of Fraxinus excelsior at five loci and the successful development of primer pairs for the amplification of three of these containing mononucleotide microsatellites. We detected high levels of haplotype variation among provenances of Fraxinus around Europe and within Ireland.

Historical and contemporary mating patterns in remnant populations of the forest tree Fraxinus excelsior L

Genetic variation at microsatellite markers was used to quantify genetic structure and mating behavior in a severely fragmented population of the wind-pollinated, wind-dispersed temperate tree Fraxinus excelsior in a deforested catchment in Scotland. Remnants maintain high levels of genetic diversity, comparable with those reported for continuous populations in southeastern Europe, and show low interpopulation differentiation (E = 0.080), indicating that historical gene exchange has not been limited (Nm = 3.48). We estimated from seeds collected from all trees producing fruits in three of five remnants that F. excelsior is predominantly outcrossing (t(m) = 0.971 +/- 0.028). Use of a neighborhood model approach to describe the relative contribution of local and long-distance pollen dispersal indicates that pollen gene flow into each of the three remnants is extensive (46-95%) and pollen dispersal has two components. The first is very localized and restricted to tens of meters around the mother trees. The second is a long-distance component with dispersal occurring over several kilometers. Effective dispersal distances, accounting for the distance and directionality to mother trees of sampled pollen donors, average 328 m and are greater than values reported for a continuous population. These results suggest that the opening of the landscape facilitates airborne pollen movement and may alleviate the expected detrimental genetic effects of fragmentation.

Chloroplast DNA variation and postglacial recolonization of common ash (Fraxinus excelsior L.) in Europe

We used chloroplast polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and chloroplast microsatellites to assess the structure of genetic variation and postglacial history across the entire natural range of the common ash (Fraxinus excelsior L.), a broad-leaved wind-pollinated and wind-dispersed European forest tree. A low level of polymorphism was observed, with only 12 haplotypes at four polymorphic microsatellites in 201 populations, and two PCR-RFLP haplotypes in a subset of 62 populations. The clear geographical pattern displayed by the five most common haplotypes was in agreement with glacial refugia for ash being located in Iberia, Italy, the eastern Alps and the Balkan Peninsula, as had been suggested from fossil pollen data. A low chloroplast DNA mutation rate, a low effective population size in glacial refugia related to ash's life history traits, as well as features of postglacial expansion were put forward to explain the low level of polymorphism. Differentiation among populations was high (GST= 0.89), reflecting poor mixing among recolonizing lineages. Therefore, the responsible factor for the highly homogeneous genetic pattern previously identified at nuclear microsatellites throughout western and central Europe (Heuertz et al. 2004) must have been efficient postglacial pollen flow. Further comparison of variation patterns at both marker systems revealed that nuclear microsatellites identified complex differentiation patterns in south-eastern Europe which remained undetected with chloroplast microsatellites. The results suggest that data from different markers should be combined in order to capture the most important genetic patterns in a species.

Nuclear microsatellites reveal contrasting patterns of genetic structure between western and southeastern European populations of the common ash (Fraxinus excelsior L.)

To determine extant patterns of population genetic structure in common ash and gain insight into postglacial recolonization processes, we applied multilocus-based Bayesian approaches to data from 36 European populations genotyped at five nuclear microsatellite loci. We identified two contrasting patterns in terms of population genetic structure: (1) a large area from the British Isles to Lithuania throughout central Europe constituted effectively a single deme, whereas (2) strong genetic differentiation occurred over short distances in Sweden and southeastern Europe. Concomitant geographical variation was observed in estimates of allelic richness and genetic diversity, which were lowest in populations from southeastern Europe, that is, in regions close to putative ice age refuges, but high in western and central Europe, that is, in more recently recolonized areas. We suggest that in southeastern Europe, restricted postglacial gene flow caused by a rapid expansion of refuge populations in a mountainous topography is responsible for the observed strong genetic structure. In contrast, admixture of previously differentiated gene pools and high gene flow at the onset of postglacial recolonization of western and central Europe would have homogenized the genetic structure and raised the levels of genetic diversity above values in the refuges.