Evolutionary relics dominate the small number of secondary metabolism genes in the hemibiotrophic fungus Dothistroma septosporum

Fungal secondary metabolites have important functions for the fungi that produce them, such as roles in virulence and competition. The hemibiotrophic pine needle pathogen Dothistroma septosporum has one of the lowest complements of secondary metabolite (SM) backbone genes of plant pathogenic fungi, indicating that this fungus produces a limited range of SMs. Amongst these SMs is dothistromin, a well-characterised polyketide toxin and virulence factor that is required for expansion of disease lesions in Dothistroma needle blight disease. Dothistromin genes are dispersed across six loci on one chromosome, rather than being clustered as for most SM genes. We explored other D. septosporum SM genes to determine if they are associated with gene clusters, and to predict what their likely products and functions might be. Of nine functional SM backbone genes in the D. septosporum genome, only four were expressed under a range of in planta and in culture conditions, one of which was the dothistromin PKS backbone gene. Of the other three expressed genes, gene knockout studies suggested that DsPks1 and DsPks2 are not required for virulence and attempts to determine a functional squalestatin-like SM product for DsPks2 were not successful. However preliminary evidence suggested that DsNps3, the only SM backbone gene to be most highly expressed in the early stage of disease, appears to be a virulence factor. Thus, despite the small number of SM backbone genes in D. septosporum, most of them appear to be poorly expressed or dispensable for virulence in planta. This work contributes to a growing body of evidence that many fungal secondary metabolite gene clusters might be non-functional and may be evolutionary relics.     

Beyond the genomes of Fulvia fulva (syn. Cladosporium fulvum) and Dothistroma septosporum: New insights into how these fungal pathogens interact with their host plants

Fulvia fulva and Dothistroma septosporum are closely related apoplastic pathogens with similar lifestyles but different hosts: F. fulva is a pathogen of tomato, whilst D. septosporum is a pathogen of pine trees. In 2012, the first genome sequences of these pathogens were published, with F. fulva and D. septosporum having highly fragmented and near-complete assemblies, respectively. Since then, significant advances have been made in unravelling their genome architectures. For instance, the genome of F. fulva has now been assembled into 14 chromosomes, 13 of which have synteny with the 14 chromosomes of D. septosporum, suggesting these pathogens are even more closely related than originally thought. Considerable advances have also been made in the identification and functional characterization of virulence factors (e.g., effector proteins and secondary metabolites) from these pathogens, thereby providing new insights into how they promote host colonization or activate plant defence responses. For example, it has now been established that effector proteins from both F. fulva and D. septosporum interact with cell-surface immune receptors and co-receptors to activate the plant immune system. Progress has also been made in understanding how F. fulva and D. septosporum have evolved with their host plants, whilst intensive research into pandemics of Dothistroma needle blight in the Northern Hemisphere has shed light on the origins, migration, and genetic diversity of the global D. septosporum population. In this review, we specifically summarize advances made in our understanding of the F. fulva–tomato and D. septosporum–pine pathosystems over the last 10 years.     

Population structure and diversity of an invasive pine needle pathogen reflects anthropogenic activity

Dothistroma septosporum is a haploid fungal pathogen that causes a serious needle blight disease of pines, particularly as an invasive alien species on Pinus radiata in the Southern Hemisphere. During the course of the last two decades, the pathogen has also incited unexpected epidemics on native and non‐native pine hosts in the Northern Hemisphere. Although the biology and ecology of the pathogen has been well documented, there is a distinct lack of knowledge regarding its movement or genetic diversity in many of the countries where it is found. In this study we determined the global population diversity and structure of 458 isolates of D. septosporum from 14 countries on six continents using microsatellite markers. Populations of the pathogen in the Northern Hemisphere, where pines are native, displayed high genetic diversities and included both mating types. Most of the populations from Europe showed evidence for random mating, little population differentiation and gene flow between countries. Populations in North America (USA) and Asia (Bhutan) were genetically distinct but migration between these continents and Europe was evident. In the Southern Hemisphere, the population structure and diversity of D. septosporum reflected the anthropogenic history of the introduction and establishment of plantation forestry, particularly with Pinus radiata. Three introductory lineages in the Southern Hemisphere were observed. Countries in Africa, that have had the longest history of pine introductions, displayed the greatest diversity in the pathogen population, indicating multiple introductions. More recent introductions have occurred separately in South America and Australasia where the pathogen population is currently reproducing clonally due to the presence of only one mating type.     

Relationship and genetic structure among autoecious and heteroecious populations of Cronartium pini in northern Fennoscandia

Epidemics of Scots pine blister rust, caused by Cronartium pini, have become an increasing problem in northern Finland and Sweden. The biology of the rust fungus is complex, with two different life cycle forms that cannot be morphologically distinguished, and it is unclear to what extent the two forms contribute to the epidemics. Genetic structure of fourteen populations of C. pini were investigated in Fennoscandia. Distinction between the two life cycle forms, a heteroecious and an autoecious one, was made by determining zygosity using microsatellite markers, and AFLP markers were developed to analyse population genetic relationships. The results showed that the two life cycle forms are clearly differentiated and occur in separate populations. Within the life cycle forms, geographic differentiation was evident, probably due to restricted gene flow as well as connection with different alternating hosts. The host-alternating form dominated in the epidemic regions in northern Fennoscandia. Implications for silvicultural practices are discussed.     

Dothistroma pini, a Forest Pathogen, Contains Homologs of Aflatoxin Biosynthetic Pathway Genes

Homologs of aflatoxin biosynthetic genes have been identified in the pine needle pathogen Dothistroma pini. D. pini produces dothistromin, a difuranoanthraquinone toxin with structural similarity to the aflatoxin precursor versicolorin B. Previous studies with purified dothistromin suggest a possible role for this toxin in pathogenicity. By using an aflatoxin gene as a hybridization probe, a genomic D. pini clone was identified that contained four dot genes with similarity to genes in aflatoxin and sterigmatocystin gene clusters with predicted activities of a ketoreductase (dotA), oxidase (dotB), major facilitator superfamily transporter (dotC), and thioesterase (dotD). A D. pini dotA mutant was made by targeted gene replacement and shown to be severely impaired in dothistromin production, confirming that dotA is involved in dothistromin biosynthesis. Accumulation of versicolorin A (a precursor of aflatoxin) by the dotA mutant confirms that the dotA gene product is involved in an aflatoxin-like biosynthetic pathway. Since toxin genes have been found to be clustered in fungi in every case analyzed so far, it is speculated that the four dot genes may comprise part of a dothistromin biosynthetic gene cluster. A fifth gene, ddhA, is not a homolog of aflatoxin genes and could be at one end of the dothistromin cluster. These genes will allow comparative biochemical and genetic studies of the aflatoxin and dothistromin biosynthetic pathways and may also lead to new ways to control Dothistroma needle blight.     

The composition of the free fatty acids from Dendrolimus pini exuviae

The pine moth Dendrolimus pini effectively resists many insecticides, but it can be controlled by the use of bioinsecticides such as entomopathogenic fungi. In the use of microbial agents for the biocontrol of D. pini, it is important to identify the cuticular lipids of this pest if we are to understand the factors responsible for the preferential adhesion or selective repulsion of entomopathogenic fungi that are potentially useful in biocontrol. In this work the qualitative and quantitative analyses of free fatty acids in two exuviae extracts (petroleum ether and dichloromethane) and two developmental stages (larval-larval and larval-pupal molts) were studied. The free fatty acid composition of the epicuticular lipids from exuviae of D. pini was characterized chemically using gas chromatography (GC) and gas chromatography-electron impact mass spectrometry (GC-MS). Structural analyses of the dichloromethane extracts from larval-larval exuviae (LLE) and larval-pupal exuviae (LPE) revealed that the carbon numbers for the major acid moieties ranged from C_8:0 to C_34:0. Only C_23:0 was not identified in the LPE extract. The relative contents of fatty acids in the extracts varied from trace amounts to 34%. The fatty acids extracted by dichloromethane were essentially the same as those in the petroleum ether extract. We also identified dehydroabietic acid in the exuviae of D. pini. The respective quantities of dehydroabietic acid obtained from D. pini LLE and LPE were 1763 ± 103 μg/g exuviae and 11521 ± 1198 μg/g of exuviae.     

A novel GFP-based approach for screening biocontrol microorganisms in vitro against Dothistroma septosporum

Dothistroma septosporum is the causal agent of Dothistroma needle blight of pine trees. A novel green fluorescent protein (GFP)-based screening method was developed to assess the potential of microorganisms for biocontrol of Dothistroma. The screen utilizes GFP expression as an indicator of metabolic activity in the pathogen and hygromycin resistance selection to determine if the interaction is fungistatic or fungicidal. Results suggested that six of eight Trichoderma isolates tested have the potential to control Dothistroma in vitro, via a fungicidal action. Because D. septosporum produces a broad-spectrum toxin, dothistromin, the inhibition of Trichoderma spp. by D. septosporum was determined by growth rate measurements compared to controls. Inhibition of the Trichoderma spp. ranged from no inhibition to 30% inhibition and was influenced by the assay medium used. The GFP screening method was also assessed to determine if it was suitable for screening bacteria as potential biocontrol candidates. Although a method involving indirect-contact had to be used, two of four Bacillus strains showed antagonistic activity against D. septosporum in vitro, via a fungistatic interaction. The four bacterial strains inhibited D. septosporum growth by 14.0 to 39.8%. This GFP-based method represents a novel approach to screening fungi and bacteria for antagonistic activity.     

Diaporthe amygdali,a species complex or a complex species?

Delimitation of species boundaries within the fungal genus Diaporthe has been challenging, but the analyses of combined multilocus DNA sequences has become an important tool to infer phylogenetic relationships and to circumscribe species. However, analyses of congruence between individual gene genealogies and the application of the genealogical concordance principle have been somehow overlooked. We noted that a group of species including D. amygdali, D. garethjonesii, D. sterilis, D. kadsurae, D. ternstroemia, D. ovoicicola, D. fusicola, D. chongqingensis and D. mediterranea, commonly known as D. amygdali complex, occupy a monophyletic clade in Diaporthe phylogenies but the limits of all species within the complex are not entirely clear. To assess the boundaries of species within this complex we employed the Genealogical Concordance Phylogenetic Species Recognition principle (GCPSR) and coalescence-based models: General Mixed Yule-Coalescent (GMYC) and Poisson Tree Processes (PTP). The incongruence detected between individual gene phylogenies, as well as the results of coalescent methods do not support the recognition of lineages within the complex as distinct species. Moreover, results support the absence of reproductive isolation and barriers to gene flow in this complex, thus providing further evidence that the D. amygdali species complex constitutes a single species. This study highlights the relevance of the application of the GCPSR principle, showing that concatenation analysis of multilocus DNA sequences, although being a powerful tool, might lead to an erroneous definition of species limits. Additionally, it further shows that coalescent methods are useful tools to assist in a more robust delimitation of species boundaries in the genus Diaporthe.     

Habitat indicators for cavity-nesters: The polypore Phellinus pini in pine forests

Foresters and arborists have long used fruit-bodies of heart-rot fungi as signs of advanced live-tree decay, but such usage has not been elaborated for forest conservation. I analysed relationships between a heart-rot fungus, tree-cavity supply, and cavity-nesting bird assemblage in wet hemiboreal Scots pine (Pinus sylvestris) forests in Estonia. The focal species, Phellinus pini, is the main heartwood decayer of live pines; it typically forms fruit-bodies at the stage of advanced decay on old trees. I found that the pine wetlands had few tree-cavities (mostly in snags) and cavity-nesting birds (woodpeckers being almost absent) despite abundant snag supply. Only one fruit-body of P. pini was found on cavity-tree but stand-scale abundance of the fruit-bodies correlated well with cavity-nester densities. Multifactor models indicated that cavity formation, not tree death, was the limiting process for secondary cavity-nesters, and P. pini could indeed be used (in combination with old-pine abundance) for assessing their habitat quality. This fungus could also serve as an educational flagship species to bridge conservation biology and forest pathology, and its fruit-bodies can signal trees to be retained at harvesting in pine forests. The conclusion is that there is hope for developing practical indicators to manage for the hidden decay processes that govern tree-cavity development.     

Oviposition preference and larval survival of Diprion pini on Scots pine clones in relation to foliage characteristics

1 When laying eggs, adult female Diprion pini discriminated between clones of Scots pine, Pinus sylvestris. 2 Larval survival was affected by Scots pine clone but was not correlated with oviposition preference. 3 Diprion pini laid most eggs on tougher needles of Scots pine. 4 There was no evidence to suggest that D. pini is responding to the monoterpene composition of the foliage.     

Forest defoliator outbreaks under climate change: effects on the frequency and severity of outbreaks of five pine insect pests

To identify general patterns in the effects of climate change on the outbreak dynamics of forest‐defoliating insect species, we examined a 212‐year record (1800–2011) of outbreaks of five pine‐defoliating species (Bupalus piniarius, Panolis flammea, Lymantria monacha, Dendrolimus pini, and Diprion pini) in Bavaria, Germany for the evidence of climate‐driven changes in the severity, cyclicity, and frequency of outbreaks. We also accounted for historical changes in forestry practices and examined effects of past insecticide use to suppress outbreaks. Analysis of relationships between severity or occurrence of outbreaks and detrended measures of temperature and precipitation revealed a mixture of positive and negative relationships between temperature and outbreak activity. Two moth species (P. flammea and Dendrolimus pini) exhibited lower outbreak activity following years or decades of unusually warm temperatures, whereas a sawfly (Diprion pini), for which voltinism is influenced by temperature, displayed increased outbreak occurrence in years of high summer temperatures. We detected only one apparent effect of precipitation, which showed Dendrolimus pini outbreaks tending to follow drought. Wavelet analysis of outbreak time series suggested climate change may be associated with collapse of L. monacha and Dendrolimus pini outbreak cycles (loss of cyclicity and discontinuation of outbreaks, respectively), but high‐frequency cycles for B. piniarius and P. flammea in the late 1900s. Regional outbreak severity was generally not related to past suppression efforts (area treated with insecticides). Recent shifts in forestry practices affecting tree species composition roughly coincided with high‐frequency outbreak cycles in B. piniarius and P. flammea but are unlikely to explain the detected relationships between climate and outbreak severity or collapses of outbreak cycles. Our results highlight both individualistic responses of different pine‐defoliating species to climate changes and some patterns that are consistent across defoliator species in this and other forest systems, including collapsing of population cycles.     

Incomplete barriers to mitochondrial gene flow between pheromone races of the North American pine engraver,Ips pini (Say) (Coleoptera,Scolytidae)

The pine engraver Ips pini (Say) is known to include three pheromone races, but gene flow between these races has not been investigated. We used maternally inherited mitochondrial DNA (mtDNA) variation to infer gene flow between 22 widely distributed North American populations of I. pini for a total of 217 individuals, based on 354 bp of the cytochrome oxidase I gene. Gene flow was estimated cladistically as migrants per generation (Nm) and as haplotype variation between populations (N_st). Three distinct mtDNA haplotype lineages, generally corresponding to eastern (I), Rocky Mountain (II) and western (III) regions of North America, were resolved with a total of 34 distinct I. pini haplotypes. The distributions of these lineages were largely congruent with the geographical ranges of the ''New York'', ''California'' and ''Idaho–Montana'' pheromone races. Only individuals with lineage I mtDNA were observed among eastern populations, whereas individuals with lineage II or III mtDNA predominated among western populations. Gene flow (Nm and N_st) was generally moderate between all populations. However, the presence of lineage I mtDNA on the eastern side of western North America and the absence of lineage II and III mtDNA in eastern North America suggest directional gene flow from east to west. These results indicate that female-controlled assortative mating among pheromone races may disrupt gene flow between conspecifics, reflecting incomplete pre-mating barriers.     

Characterization of the ascomycetes Therrya fuckelii and T. pini fruiting on Scots pine branches in Nordic countries

The relative frequency of Therrya fuckelii and T. pini fruiting on dead branches of Scots pine was investigated in southern Norway by examining lightning-damaged and wind-fallen trees, randomly collected branches and Nordic herbarium collections of these ascomycetes representing the order Rhytismatales. Ascus, ascospore, and subhymenium characteristics were used as criteria for species identification, while a sequence analysis of ITS rDNA gene cluster was performed to compare the relatedness of the species to each other and to corresponding fungal sequences available at the NCBI GenBank Sequence Database. In a few cases, the two Therrya species co-occurred on the same branch, but in general, whether field or herbarium material, T. fuckelii was clearly more common than T. pini. Within the Nordic countries, both species occurred throughout the natural distribution area of Scots pine. The ITS rDNA sequence of T. pini strains was 91% similar to T. fuckelii strains, the differences locating both within the internal transcribed spacers ITS1 and ITS2 and the 5.8 S rDNA gene. More variation in the ITS1-5.8S-ITS2 sequence was observed among T. pini than T. fuckelii samples; genetic implications of this finding are discussed. Upon sequence analysis, we discovered that a T. pini sequence has been deposited in the NCBI GenBank under a false identity. We emphasize the importance of co-examining strains that originate from mature fruit bodies with fully developed morphologic features as reference samples.     

Genetic structure of loggerhead turtle (Caretta caretta) populations in Turkey

This study investigates genetic structure of loggerhead turtle (Caretta caretta) for conservation purposes using both mitochondrial and nuclear DNA markers. To evaluate overall genetic structure of C. caretta at the Turkish nesting beaches and to compare it with previously published data from other nesting sites of the Mediterranean Sea, we have studied 256 hatchlings from 18 different nesting beaches. Seven distinct haplotypes were detected, of which three have previously been reported from the Mediterranean and one from the Atlantic. The remaining three haplotypes are described for the first time for C. caretta in the present study. Distribution of these haplotypes among the nesting sites showed a significant genetic structuring, indicating that females are philopatric and that gene flow among populations is restricted. Both mtDNA and microsatellite analyses determined genetic structuring (mtDNA: γ_st = 0.214, p < 0.01; nDNA F_st = 0.0004 p < 0.05) among nesting aggregates of C. caretta throughout the study area and enabled the detection of the different haplotypes to inform conservation strategies.     

Biosynthesis of dothistromin

Dothistromin is a mycotoxin that is remarkably similar in structure to versicolorin B, a precursor of both aflatoxin and sterigmatocystin. Dothistromin-producing fungi also produce related compounds, including some aflatoxin precursors as well as alternative forms of dothistromin. Dothistromin is synthesized by pathogenic species of Dothistroma in the red bands of pine needles associated with needle blight, but is also made in culture where it is strongly secreted into the surrounding medium. Orthologs of aflatoxin and sterigmatocystin biosynthetic genes have been found that are required for the biosynthesis of dothistromin, along with others that are speculated to be involved in the same pathway on the basis of their sequence similarity to aflatoxin genes. An epoxide hydrolase gene that has no homolog in the aflatoxin or sterigmatocystin gene clusters is also clustered with the dothistromin genes, and all these genes appear to be located on a minichromosome in Dothistroma septosporum. The dothistromin genes are expressed at an early stage of growth, suggesting a role in the first stages of plant invasion by the fungus. Future studies are expected to reveal more about the role of dothistromin in needle blight and about the genomic organization and expression of dothistromin genes: these studies will provide for interesting comparisons with these aspects of aflatoxin and sterigmatocystin biosynthesis.     

Cytonuclear discordance and the species status of Myotis myotis and Myotis blythii (Chiroptera)

We used both highly variable mitochondrial and nuclear loci to investigate the large mouse‐eared bat species complex in the Western Palaearctic to clarify their systematic position. Although mitochondrial lineages show no species segregation and some haplotypes are shared between Myotis myotis and Myotis blythii sensu lato, Bayesian clustering methods based on multilocus genotypes indicate highly concordant nuclear and morphological species assignment. These multilocus, nuclear analyses detected only a single putative F₁ hybrid in the extensive areas of sympatry sampled, thus confirming the biological species status of M. myotis and M. blythii s.l. We propose that the strong cytonuclear discordance in these species complex results from a combination of prior spatial isolation of the two species in different glacial refugia, followed by a succession of mitochondrial introgression events that occurred during the eastward and westward expansions of M. myotis and of M. blythii, respectively. The nuclear markers further indicate the presence of a notable genetic discontinuity within M. myotis that broadly separates populations into an eastern and a western component with an overlap zone in the Balkans. This eastern and western discontinuity is also apparent in the mitochondrial lineages with the D haplogroup largely confined to samples found in Thrace and Asia Minor. None of these genetic discontinuities correspond to the distribution of the two commonly recognized M. myotis subspecies (myotis and macrocephalicus). We also show that distinct morphological subspecies within M. blythii (oxygnathus, omari, risorius and lesviacus) in Europe and the near‐East do not correlate with significant evolutionarily units, whether identified by mitochondrial or nuclear data and thus only represent local morphological variants with little taxonomic relevance.     

Description of Meloidogyne pini n. sp., a Root-Knot Nematode Parasitic on Sand Pine (Pinus clausa), with Additional Notes on the Morphology of M. megatyla

Meloidogyne pini n. sp. is described from sand pine, Pinus clausa, in Georgia. The perineal pattern of the female has a large cuticular ridge surrounding a deeply recessed perivulval area. The lateral fields are marked by transverse striae. The female stylet is 14.6 μm long, and the knobs are small, rounded, and set off from the straight and narrow shaft. The excretory pore is near the level of the base of the stylet. The labial disc of the male is large, rounded, and fused with the crescent-shaped medial lips. The head region is smooth, the styler is 20.8 μm long, and the cone is more than twice as long as the shaft. The knobs are rounded and set off from the shaft. In the second-stage juvenile, the labial disc, medial lips, and lateral lips form one smooth, continuous, ovoid head cap. Mean juvenile length is 434 μm, stylet length is 12.8 μm, and tail length is 44.4 μm. M. pini n. sp. also parasitizes loblolly and slash pine. Additional morphological details of M. megatyla are presented.     

Differentiation of Fusarium subglutinans f. sp. pini by Histone Gene Sequence Data

Fusarium subglutinans f. sp. pini (= F. circinatum) is a pathogen of pine and is one of eight mating populations (i.e., biological species) in the Gibberella fujikuroi species complex. This species complex includes F. thapsinum, F. moniliforme (= F. verticillioides), F. nygamai, and F. proliferatum, as well as F. subglutinans associated with sugarcane, maize, mango, and pineapple. Differentiating these forms of F. subglutinans usually requires pathogenicity tests, which are often time-consuming and inconclusive. Our objective was to develop a technique to differentiate isolates of F. subglutinans f. sp. pini from other isolates identified as F. subglutinans. We sequenced the histone H3 gene from a representative set of Fusarium isolates. The H3 gene sequence was conserved and contained two introns in all the isolates studied. From both the intron and the exon sequence data, we developed a PCR-restriction fragment length polymorphism technique that reliably distinguishes F. subglutinans f. sp. pini from the other biological species in the G. fujikuroi species complex.     

Functional characterization of myrcene hydroxylases from two geographically distinct Ips pini populations

Ips pini bark beetles use myrcene hydroxylases to produce the aggregation pheromone component, ipsdienol, from myrcene. The enantiomeric ratio of pheromonal ipsdienol is an important prezygotic mating isolation mechanism of I. pini and differs among geographically distinct populations. We explored the substrate and product ranges of myrcene hydroxylases (CYP9T2 and CYP9T3) from reproductively-isolated western and eastern I. pini. The two cytochromes P450 share 94% amino acid identity. CYP9T2 mRNA levels were not induced in adults exposed to myrcene-saturated atmosphere. Functional assays of recombinant enzymes showed both hydroxylated myrcene, (+)- and (?)-α-pinene, 3-carene, and R-(+)-limonene, but not α-phellandrene, (?)-β-pinene, γ-terpinene, or terpinolene, with evidence that CYP9T2 strongly preferred myrcene over other substrates. They differed in the enantiomeric ratios of ipsdienol produced from myrcene, and in the products resulting from different α-pinene enantiomers. These data provide new information regarding bark beetle pheromone evolution and factors affecting cytochrome P450 structure–function relationships.