The complete mitochondrial genome of the vascular wilt fungus Verticillium dahliae: a novel gene order for Verticillium and a diagnostic tool for species identification

The complete mitogenome of the entomopathogenic fungus Metarhizium pinghaense 15R

In this study, the complete mitogenome of the entomopathogenic fungus Metarhizium pinghaense 15 R, which is highly virulent to aphids and was isolated from Korean soil, was assembled and annotated for three ATP synthase subunits (atp6, atp8, and atp9), three cytochrome oxidase subunits (cox1, cox2, and cox3), apocytochrome b (cob), seven subunits of NADH dehydrogenase (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), two ribosomal RNAs (rnl and rns), and 19 tRNA genes. Five genes were carrying a total of eight introns, and they may encode ribosomal protein S3, LAGLIDADG and GIY-YIG endonucleases. Phylogenetic analysis based on the mitochondrial nucleotide sequence confirmed that the M. pinghaense 15 R is a member of the Clavicipitaceae, and is closely related to the species M. anisopliae, M. robertsii, and M. brunneum. The mtDNA base sequence of the M. pinghaense 15 R strain reported in this study is thought to be useful for biological resource genetic data.

Functional annotation and comparative analysis of four Botrytis cinerea mitogenomes reported from Punjab, Pakistan

Botrytis cinerea is one of the top phytopathogenic fungus which ubiquitously cause grey mold on a variety of horticultural plants. The mechanism of respiration in the fungus occurs within the mitochondria. Mitogenomes serve as a key molecular marker for the investigation of fungal evolutionary patterns. This study aimed at the complete assembly, characterization, and comparative relationship of four mitogenomes of Botrytis cinerea strains including Kst5C, Kst14A, Kst32B, Kst33A, respectively. High throughput sequencing of four mitogenomes allowed the full assembly and annotation of these sequences. The total genome length of these 4 isolates Kst5C Kst14A, Kst32B, Kst33A was 69,986 bp, 77,303 bp, 76,204 bp and 55, 226 bp respectively. The distribution of features represented 2 ribosomal RNA genes,14 respiration encoding proteins, 1 mitochondrial ribosomal protein-encoding gene, along with varying numbers of transfer RNA genes, protein-coding genes, mobile intronic regions and homing endonuclease genes including LAGLIDADG and GIY-YIG domains were found in all four mitogenomes. The comparative analyses performed also decipher significant results for four mitogenomes among fungal isolates included in the study. This is the first report on the detailed annotation of mitogenomes as a proof for investigation of variation patterns present with in the B. cinerea causing grey mold on strawberries in Pakistan. This study will also contribute to the rapid evolutionary analysis and population patterns present among Botrytis cinerea.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the pathogenic fungus Ilyonectria destructans

Ilyonectria destructans is a pathogenic fungus causing root rot and other symptoms on trees and many crops. This paper analyses the mitochondrial genome of I. destructans and compares it with other published Nectriaceae mitogenomes. The I. destructans mitogenome appears as a circular DNA molecule of 42,895 bp and an overall GC content of 28.23%. It contains 28 protein-coding genes (15 core protein genes and 13 free-standing ORFs), two rRNAs and 27 tRNAs. The gene content and order were found to be conserved in the mitogenome of I. destructans and other Nectriaceae, although the genome size varies because of the variation in the number and length of intergenic regions and introns. For most core protein-coding genes in Nectriaceae species, Ka/Ks < 1 indicates purifying selection. Among some Nectriaceae representatives, only the rps3 gene was found under positive selection. Phylogenetic analyses based on nucleotide sequences of 15 protein-coding genes divided 45 Hypocreales species into six major clades matching the families Bionectriaceae, Cordycipitaceae, Clavicipitaceae, Ophiocordycipitaceae, Hypocreaceae and Nectriaceae. I. destructans appeared as a sister species to unidentified Ilyonectia sp., closely related to C. ilicicola, N. cinnabarina and a clad of ten Fusarium species and G. moniliformis. The complete mitogenome of I. destructans reported in the current paper will facilitate the study of epidemiology, biology, genetic diversity of the species and the evolution of family Nectriace and the Hypocreales order.

The First Mitochondrial Genome of Ciborinia camelliae and Its Position in the Sclerotiniaceae Family

Ciborinia camelliae is the causal agent of camellia flower blight (CFB). It is a hemibiotrophic pathogen, inoperculate Discomycete of the family Sclerotiniaceae. It shows host and organ specificity infecting only flowers of species belonging to the genus Camellia, causing serious damage to the ornamental component of the plant. In this work, the first mitochondrial genome of Ciborinia camellia is reported. The mitogenome was obtained by combining Illumina short read and Nanopore long read technology. To resolve repetitive elements, specific primers were designed and used for Sanger sequencing. The manually curated mitochondrial DNA (mtDNA) of the Italian strain DSM 112729 is a circular sequence of 114,660 bp, with 29.6% of GC content. It contains two ribosomal RNA genes, 33 transfer RNAs, one RNase P gene, and 62 protein-coding genes. The latter include one gene coding for a ribosomal protein (rps3) and the 14 typical proteins involved in the oxidative metabolism. Moreover, a partial mtDNA assembled from a contig list was obtained from the deposited genome assembly of a New Zealand strain of C. camelliae. The present study contributes to understanding the mitogenome arrangement and the evolution of this phytopathogenic fungus in comparison to other Sclerotiniaceae species and confirms the usefulness of mitochondrial analysis to define phylogenetic positioning of this newly sequenced species.

Comparative Mitogenomics of Fungal Species in Stachybotryaceae Provides Evolutionary Insights into Hypocreales

Stachybotrys chartarum is one of the world’s ten most feared fungi within the family Stachybotryaceae, although to date, not a single mitogenome has been documented for Stachybotryaceae. Herein, six mitogenomes of four different species in Stachybotryaceae are newly reported. The S. chartarum mitogenome was 30.7 kb in length and contained two introns (one each in rnl and cox1). A comparison of the mitogenomes of three different individuals of S. chartarum showed few nucleotide variations and conservation of gene content/order and intron insertion. A comparison of the mitogenomes of four different Stachybotryaceae species (Memnoniella echinata, Myrothecium inundatum, S. chartarum, and S. chlorohalonata), however, revealed variations in intron insertion, gene order/content, and nad2/nad3 joining pattern. Further investigations on all Hypocreales species with available mitogenomes showed greater variabilities in gene order (six patterns) and nad2/nad3 joining pattern (five patterns) although a dominant pattern always existed in each case. Ancestral state estimation showed that in each case the dominant pattern was always more ancestral than those rare patterns. Phylogenetic analyses based on mitochondrion-encoded genes supported the placement of Stachybotryaceae in Hypocreales. The crown age of Stachybotryaceae was estimated to be approximately the Early Cretaceous (141–142 Mya). This study greatly promotes our understanding of the evolution of fungal species in Hypocreales.

The Gene Rearrangement, Loss, Transfer, and Deep Intronic Variation in Mitochondrial Genomes of Conidiobolus

The genus Conidiobolus s.s. was newly delimited from Conidiobolus s.l. In order to gain insight into its mitochondrial genetic background, this study sequenced six mitochondrial genomes of the genus Conidiobolus s.s. These mitogenomes were all composed of circular DNA molecules, ranging from 29,253 to 48,417 bp in size and from 26.61 to 27.90% in GC content. The order and direction for 14 core protein-coding genes (PCGs) were identical, except for the atp8 gene lost in Conidiobolus chlamydosporus, Conidiobolus polyspermus, and Conidiobolus polytocus, and rearranged in the other Conidiobolus s.s. species. Besides, the atp8 gene split the cox1 gene in Conidiobolus taihushanensis. Phylogenomic analysis based on the 14 core PCGs confirmed that all Conidiobolus s.s. species formed a monophyly in the Entomophthoromycotina lineage. The number and length of introns were the main factors contributing to mitogenomic size, and deep variations and potential transfer were detected in introns. In addition, gene transfer occurred between the mitochondrial and nuclear genomes. This study promoted the understanding of the evolution and phylogeny of the Conidiobolus s.s. genus.

Comparison of mitochondrial genomes provides insights into intron dynamics and evolution in Botryosphaeria dothidea and B. kuwatsukai

Botryosphaeria dothidea is one of the most common fungal pathogens on a large number of hosts worldwide. Botryosphaeria dothidea and B. kuwatsukai are also the main causal agents of apple ring rot. In this study, we sequenced, assembled and annotated the circular mitogenomes of 12 diverse B. dothidea isolates (105.7-114.8 kb) infecting various plants including apple, and five diverse B. kuwatsukai isolates (118.0-124.6 kb) from apple. B. dothidea mitogenomes harboured a set of 29-31 introns and 48-52 ORFs. In contrast, B. kuwatsukai mitogenomes harboured more introns (32-34) and ORFs (51-54). The variation in mitogenome sizes was associated mainly with different numbers of introns and insertions of mobile genetic elements. Interestingly, B. dothidea and B. kuwatsukai displayed distinct intron distribution patterns, with three intron loci showing presence/absence dynamics in each species. Large numbers of introns (57% in B. dothidea and 49% in B. kuwatsukai) were most likely obtained through horizontal transfer from non-Dothideomycetes. The mitochondrial gene phylogeny supported the differentiation of the two species. Overall, this study sheds light into the mitochondrial evolution of the plant pathogens B. dothidea and B. kuwatsukai, and intron distribution patterns could be useful markers for studies on population diversity.

Fungal Mitogenomes: Relevant Features to Planning Plant Disease Management

Mitochondrial genomes (mt-genomes) are characterized by a distinct codon usage and their autonomous replication. Mt-genomes encode highly conserved genes (mt-genes), like proteins involved in electron transport and oxidative phosphorylation but they also carry highly variable regions that are in part responsible for their high plasticity. The degree of conservation of their genes is such that they allow the establishment of phylogenetic relationships even across distantly related species. Here, we describe the mechanisms that generate changes along mt-genomes, which play key roles at enlarging the ability of fungi to adapt to changing environments. Within mt-genomes of fungal pathogens, there are dispensable as well as indispensable genes for survival, virulence and/or pathogenicity. We also describe the different complexes or mechanisms targeted by fungicides, thus addressing a relevant issue regarding disease management. Despite the controversial origin and evolution of fungal mt-genomes, the intrinsic mechanisms and molecular biology involved in their evolution will help to understand, at the molecular level, the strategies for fungal disease management.

Comparative Analyses of Mitochondrial Genomes Provide Evolutionary Insights Into Nematode-Trapping Fungi

Predatory fungi in Orbiliaceae (Ascomycota) have evolved a diversity of trapping devices that enable them to trap and kill nematodes, other small animals, and protozoans. These trapping devices include adhesive hyphae, adhesive knobs, adhesive networks, constricting rings, and non-constricting rings. Their diversity and practical importance have attracted significant attention from biologists, making them excellent model organisms for studying adaptative evolution and as biological control agents against parasitic nematodes. The putative origins and evolutionary relationships among these carnivorous fungi have been investigated using nuclear protein-encoding genes, but their patterns of mitogenome relationships and divergences remain unknown. Here we analyze and compare the mitogenomes of 12 fungal strains belonging to eight species, including six species representing all four types of nematode trapping devices and two from related but non-predatory fungi. All 12 analyzed mitogenomes were of circular DNA molecules, with lengths ranging from 146,101 bp to 280,699 bp. Gene synteny analysis revealed several gene rearrangements and intron transfers among the mitogenomes. In addition, the number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among these mitogenomes. The increased number and total size of introns were the main contributors to the length differences among the mitogenomes. Phylogenetic analyses of the protein-coding genes indicated that mitochondrial and nuclear genomes evolved at different rates, and signals of positive selection were found in several genes involved in energy metabolism. Our study provides novel insights into the evolution of nematode-trapping fungi and shall facilitate further investigations of this ecologically and agriculturally important group of fungi.

Mitochondrial genome in Hypsizygus marmoreus and its evolution in Dikarya

BACKGROUND

Hypsizygus marmoreus, a high value commercialized edible mushroom is widely cultivated in East Asia, and has become one of the most popular edible mushrooms because of its rich nutritional and medicinal value. Mitochondria are vital organelles, and play various essential roles in eukaryotic cells.

RESULTS

In this study, we provide the Hypsizygus marmoreus mitochondrial (mt) genome assembly: the circular sequence is 102,752 bp in size and contains 15 putative protein-coding genes, 2 ribosomal RNAs subunits and 28 tRNAs. We compared the mt genomes of the 27 fungal species in the Pezizomycotina and Basidiomycotina subphyla, with the results revealing that H. marmoreus is a sister to Tricholoma matsutake and the phylogenetic distribution of this fungus based on the mt genome. Phylogenetic analysis shows that Ascomycetes mitochondria started to diverge earlier than that of Basidiomycetes and supported the robustness of the hyper metric tree. The fungal sequences are highly polymorphic and gene order varies significantly in the dikarya data set, suggesting a correlation between the gene order and divergence time in the fungi mt genome. To detect the mt genome variations in H. marmoreus, we analyzed the mtDNA sequences of 48 strains. The phylogeny and variation sited type statistics of H. marmoreus provide clear-cut evidence for the existence of four well-defined cultivations isolated lineages, suggesting female ancestor origin of H. marmoreus. Furthermore, variations on two loci were further identified to be molecular markers for distinguishing the subgroup containing 32 strains of other strains. Fifteen conserved protein-coding genes of mtDNAs were analyzed, with fourteen revealed to be under purifying selection in the examined fungal species, suggesting the rapid evolution was caused by positive selection of this gene.

CONCLUSIONS

Our studies have provided new reference mt genomes and comparisons between species and intraspecies with other strains, and provided future perspectives for assessing diversity and origin of H. marmoreus.

Detection of Verticillium species in Swedish soils using real-time PCR

Verticillium species are soilborne plant pathogens, responsible for big yield losses worldwide. Here, we report improved procedures to generate DNA from Verticillium species imbedded in farm soils. Using new genomic sequence information, primers for V. dahliae, V. albo-atrum, V. tricorpus, and V. longisporum were designed. In a survey of 429 samples from intensively farmed soil of two Swedish regions, only V. dahliae and V. longisporum were identified. A bias towards V. longisporum (40%) was seen in the south, whereas V. dahliae was more frequent in the western region (19%). Analyses of soil and leaf samples from 20 sugar beet fields, where foliar wilting had been observed, revealed V. dahliae DNA in all leaf and soil samples and V. longisporum in 18 soil samples, illustrating host choice and longevity of the V. longisporum microsclerotia. This study demonstrates the applicability of new molecular diagnostic tools that are important for growers of variable crops.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

Intron Derived Size Polymorphism in the Mitochondrial Genomes of Closely Related Chrysoporthe Species

In this study, the complete mitochondrial (mt) genomes of Chrysoporthe austroafricana (190,834 bp), C. cubensis (89,084 bp) and C. deuterocubensis (124,412 bp) were determined. Additionally, the mitochondrial genome of another member of the Cryphonectriaceae, namely Cryphonectria parasitica (158,902 bp), was retrieved and annotated for comparative purposes. These genomes showed high levels of synteny, especially in regions including genes involved in oxidative phosphorylation and electron transfer, unique open reading frames (uORFs), ribosomal RNAs (rRNAs) and transfer RNAs (tRNAs), as well as intron positions. Comparative analyses revealed signatures of duplication events, intron number and length variation, and varying intronic ORFs which highlighted the genetic diversity of mt genomes among the Cryphonectriaceae. These mt genomes showed remarkable size polymorphism. The size polymorphism in the mt genomes of these closely related Chrysoporthe species was attributed to the varying number and length of introns, coding sequences and to a lesser extent, intergenic sequences. Compared to publicly available fungal mt genomes, the C. austroafricana mt genome is the second largest in the Ascomycetes thus far.

Complete mitochondrial genome of the Verticillium-wilt causing plant pathogen Verticillium nonalfalfae

Verticillium nonalfalfae is a fungal plant pathogen that causes wilt disease by colonizing the vascular tissues of host plants. The disease induced by hop isolates of V. nonalfalfae manifests in two different forms, ranging from mild symptoms to complete plant dieback, caused by mild and lethal pathotypes, respectively. Pathogenicity variations between the causal strains have been attributed to differences in genomic sequences and perhaps also to differences in their mitochondrial genomes. We used data from our recent Illumina NGS-based project of genome sequencing V. nonalfalfae to study the mitochondrial genomes of its different strains. The aim of the research was to prepare a V. nonalfalfae reference mitochondrial genome and to determine its phylogenetic placement in the fungal kingdom. The resulting 26,139 bp circular DNA molecule contains a full complement of the 14 "standard" fungal mitochondrial protein-coding genes of the electron transport chain and ATP synthase subunits, together with a small rRNA subunit, a large rRNA subunit, which contains ribosomal protein S3 encoded within a type IA-intron and 26 tRNAs. Phylogenetic analysis of this mitochondrial genome placed it in the Verticillium spp. lineage in the Glomerellales group, which is also supported by previous phylogenetic studies based on nuclear markers. The clustering with the closely related Verticillium dahliae mitochondrial genome showed a very conserved synteny and a high sequence similarity. Two distinguishing mitochondrial genome features were also found-a potential long non-coding RNA (orf414) contained only in the Verticillium spp. of the fungal kingdom, and a specific fragment length polymorphism observed only in V. dahliae and V. nubilum of all the Verticillium spp., thus showing potential as a species specific biomarker.

Characterization and phylogenetic analysis of the mitochondrial genome of Shiraia bambusicola reveals special features in the order of pleosporales

Shiraia bambusicola P. Henn. is a pathogenic fungus of bamboo, and its fruiting bodies are regarded as folk medicine. We determined and analyzed its complete mitochondrial DNA sequence (circular DNA molecule of 39,030 bp, G + C content of 25.19%). It contains the typical genes encoding proteins involved in electron transport and coupled oxidative phosphorylation (nad1-6 and nad4L, cob and cox1-3), one ATP synthase subunit (atp6), 4 hypothetical proteins, and two genes for large and small rRNAs (rnl and rns). There is a set of 32 tRNA genes comprising all 20 amino acids, and these genes are evenly distributed on the two strands. Phylogenetic analyses based on concatenated mitochondrial proteins indicated that S. bambusicola clustered with members of the order Pleosporales, which is in agreement with previous results. The gene arrangements of Dothideomycetes species contained three regions of gene orders partitioned in their mitochondrial genomes, including block 1 (nad6-atp6), block 2 (nad1-cox3) and block 3 (genes around rns). S. bambusicola displayed unique special features that differed from the other Pleosporales species, especially in the coding regions around rns (trnR-trnY). Moreover, a comparison of gene orders in mitochondrial genomes from Pezizomycotina revealed that although all encoded regions are located on the same strand in most Pezizomycotina mtDNAs, genes from Dothideomycetes species had different orientations, as well as diverse positions and colocalization of genes (such as cox3, cox1-cox2 and nad2-nad3); these distinctions were regarded as class-specific features. Interestingly, two incomplete copies of the atp6 gene were found on different strands of the mitogenomic DNA, a finding that has not been observed in the other analyzed fungal species. In our study, mitochondrial genomes from Dothideomycetes species were comprehensively analyzed for the first time, including many species that have not appeared in previous reports.

Analysis of the complete mitochondrial genome of Pochonia chlamydosporia suggests a close relationship to the invertebrate-pathogenic fungi in Hypocreales

Background

The fungus Pochonia chlamydosporia parasitizes nematode eggs and has become one of the most promising biological control agents (BCAs) for plant-parasitic nematodes, which are major agricultural pests that cause tremendous economic losses worldwide. The complete mitochondrial (mt) genome is expected to open new avenues for understanding the phylogenetic relationships and evolution of the invertebrate-pathogenic fungi in Hypocreales.

Results

The complete mitogenome sequence of P. chlamydosporia is 25,615 bp in size, containing the 14 typical protein-coding genes, two ribosomal RNA genes, an intronic ORF coding for a putative ribosomal protein (rps3) and a set of 23 transfer RNA genes (trn) which recognize codons for all amino acids. Sequence similarity studies and syntenic gene analyses show that 87.02% and 58.72% of P. chlamydosporia mitogenome sequences match 90.50% of Metarhizium anisopliae sequences and 61.33% of Lecanicillium muscarium sequences with 92.38% and 86.04% identities, respectively. A phylogenetic tree inferred from 14 mt proteins in Pezizomycotina fungi supports that P. chlamydosporia is most closely related to the entomopathogenic fungus M. anisopliae. The invertebrate-pathogenic fungi in Hypocreales cluster together and clearly separate from a cluster comprising plant-pathogenic fungi (Fusarium spp.) and Hypocrea jecorina. A comparison of mitogenome sizes shows that the length of the intergenic regions or the intronic regions is the major size contributor in most of mitogenomes in Sordariomycetes. Evolutionary analysis shows that rps3 is under positive selection, leading to the display of unique evolutionary characteristics in Hypocreales. Moreover, the variability of trn distribution has a clear impact on gene order in mitogenomes. Gene rearrangement analysis shows that operation of transposition drives the rearrangement events in Pezizomycotina, and most events involve in trn position changes, but no rearrangement was found in Clavicipitaceae.

Conclusions

We present the complete annotated mitogenome sequence of P. chlamydosporia. Based on evolutionary and phylogenetic analyses, we have determined the relationships between the invertebrate-pathogenic fungi in Hypocreales. The invertebrate-pathogenic fungi in Hypocreales referred to in this paper form a monophyletic group sharing a most recent common ancestor. Our rps3 and trn gene order results also establish a foundation for further exploration of the evolutionary trajectory of the fungi in Hypocreales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0341-8) contains supplementary material, which is available to authorized users.

The 203 kbp mitochondrial genome of the phytopathogenic fungus Sclerotinia borealis reveals multiple invasions of introns and genomic duplications

Here we report the complete sequence of the mitochondrial (mt) genome of the necrotrophic phytopathogenic fungus Sclerotinia borealis, a member of the order Helotiales of Ascomycetes. The 203,051 bp long mtDNA of S. borealis represents one of the largest sequenced fungal mt genomes. The large size is mostly determined by the presence of mobile genetic elements, which include 61 introns. Introns contain a total of 125,394 bp, are scattered throughout the genome, and are found in 12 protein-coding genes and in the ribosomal RNA genes. Most introns contain complete or truncated ORFs that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Integrations of mobile elements are also evidenced by the presence of two regions similar to fragments of inverton-like plasmids. Although duplications of some short genome regions, resulting in the appearance of truncated extra copies of genes, did occur, we found no evidences of extensive accumulation of repeat sequences accounting for mitochondrial genome size expansion in some other fungi. Comparisons of mtDNA of S. borealis with other members of the order Helotiales reveal considerable gene order conservation and a dynamic pattern of intron acquisition and loss during evolution. Our data are consistent with the hypothesis that horizontal DNA transfer has played a significant role in the evolution and size expansion of the S. borealis mt genome.

A novel and rapid loop-mediated isothermal amplification assay for the specific detection of Verticillium dahliae

AIMS

In this study, a loop-mediated isothermal amplification (LAMP) assay has been developed and evaluated for the rapid and sensitive detection of Verticillium dahliae Kleb., the causal agent of vascular wilts in many economically important crops.

METHODS AND RESULTS

LAMP primers were designed based on a previously described RAPD marker, and the LAMP assay was applied for direct detection of V. dahliae grown on medium and from soil samples without DNA purification steps (direct-LAMP). Thirty-two agricultural soil samples from various olive orchards were collected, and the presence of pathogen was detected by LAMP, direct-LAMP and nested-PCR methods. The LAMP methodology could successfully detect V. dahliae with high specificity, and cross-reaction was not observed with different pathogenic and nonpathogenic fungi and bacteria. The LAMP assay was capable of detecting a minimum of 500 and 50 fg of purified target DNA per reaction of V. dahliae ND and D pathotypes, respectively. In contrast, nested-PCR could only detect 5 pg reaction(-1) for both pathotypes. In artificially infested soil samples, the LAMP method detected 5 microsclerotia per gram of soil. Conversely, nested-PCR assay detected 50 microsclerotia g(-1) soil. The detection ratios of LAMP and direct-LAMP protocols were better (26 and 24 positive samples out of 32 agricultural soils analysed, respectively) than that obtained for nested-PCR method (22 positive results). Moreover, direct-LAMP yielded positive detection of V. dahliae in agricultural soil samples within 60-80 min.

CONCLUSIONS

The newly developed LAMP method was proved to be an effective, simple and rapid method to detect V. dahliae without the need for either expensive equipment or DNA purification.

SIGNIFICANCE AND IMPACT OF STUDY

This technique can be considered as an excellent standard alternative to plating and nested-PCR assays for the early, sensitive and low-cost detection of V. dahliae as well as other soilborne pathogens in the field.

Characterization and phylogenetic analysis of the mitochondrial genome of Glarea lozoyensis indicates high diversity within the order Helotiales

BACKGROUND

Glarea lozoyensis is a filamentous fungus used for the industrial production of non-ribosomal peptide pneumocandin B0. In the scope of a whole genome sequencing the complete mitochondrial genome of the fungus has been assembled and annotated. It is the first one of the large polyphyletic Helotiaceae family. A phylogenetic analysis was performed based on conserved proteins of the oxidative phosphorylation system in mitochondrial genomes.

RESULTS

The total size of the mitochondrial genome is 45,038 bp. It contains the expected 14 genes coding for proteins related to oxidative phosphorylation,two rRNA genes, six hypothetical proteins, three intronic genes of which two are homing endonucleases and a ribosomal protein rps3. Additionally there is a set of 33 tRNA genes. All genes are located on the same strand. Phylogenetic analyses based on concatenated mitochondrial protein sequences confirmed that G. lozoyensis belongs to the order of Helotiales and that it is most closely related to Phialocephala subalpina. However, a comparison with the three other mitochondrial genomes known from Helotialean species revealed remarkable differences in size, gene content and sequence. Moreover, it was found that the gene order found in P. subalpina and Sclerotinia sclerotiorum is not conserved in G. lozoyensis.

CONCLUSION

The arrangement of genes and other differences found between the mitochondrial genome of G. lozoyensis and those of other Helotiales indicates a broad genetic diversity within this large order. Further mitochondrial genomes are required in order to determine whether there is a continuous transition between the different forms of mitochondrial genomes or G. lozoyensis belongs to a distinct subgroup within Helotiales.

Evidence for inter-specific recombination among the mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex

Background

The availability of mitochondrial genomes has allowed for the resolution of numerous questions regarding the evolutionary history of fungi and other eukaryotes. In the Gibberella fujikuroi species complex, the exact relationships among the so-called “African”, “Asian” and “American” Clades remain largely unresolved, irrespective of the markers employed. In this study, we considered the feasibility of using mitochondrial genes to infer the phylogenetic relationships among Fusarium species in this complex. The mitochondrial genomes of representatives of the three Clades (Fusarium circinatum, F. verticillioides and F. fujikuroi) were characterized and we determined whether or not the mitochondrial genomes of these fungi have value in resolving the higher level evolutionary relationships in the complex.

Results

Overall, the mitochondrial genomes of the three species displayed a high degree of synteny, with all the genes (protein coding genes, unique ORFs, ribosomal RNA and tRNA genes) in identical order and orientation, as well as introns that share similar positions within genes. The intergenic regions and introns generally contributed significantly to the size differences and diversity observed among these genomes. Phylogenetic analysis of the concatenated protein-coding dataset separated members of the Gibberella fujikuroi complex from other Fusarium species and suggested that F. fujikuroi (“Asian” Clade) is basal in the complex. However, individual mitochondrial gene trees were largely incongruent with one another and with the concatenated gene tree, because six distinct phylogenetic trees were recovered from the various single gene datasets.

Conclusion

The mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex are remarkably similar to those of the previously characterized Fusarium species and Sordariomycetes. Despite apparently representing a single replicative unit, all of the genes encoded on the mitochondrial genomes of these fungi do not share the same evolutionary history. This incongruence could be due to biased selection on some genes or recombination among mitochondrial genomes. The results thus suggest that the use of individual mitochondrial genes for phylogenetic inference could mask the true relationships between species in this complex.

Evaluation of mitochondrial genes as DNA barcode for Basidiomycota

Our study evaluated in silico the potential of 14 mitochondrial genes encoding the subunits of the respiratory chain complexes, including cytochrome c oxidase I (CO1), as Basidiomycota DNA barcode. Fifteen complete and partial mitochondrial genomes were recovered and characterized in this study. Mitochondrial genes showed high values of molecular divergence, indicating a potential for the resolution of lower-level relationships. However, numerous introns occurred in CO1 as well as in six other genes, potentially interfering with polymerase chain reaction amplification. Considering these results and given the minimal length of 600-bp that is optimal for a fungal barcode, the genes encoding for the ATPase subunit 6, the cytochrome oxidase subunit 3 and the NADH dehydrogenase subunit 6 have the most promising characteristics for DNA barcoding among the mitochondrial genes studied. However, biological validation on two fungal data sets indicated that no single mitochondrial gene gave a better taxonomic resolution than the ITS, the region already widely used in fungal taxonomy.

Global analyses of Ceratocystis cacaofunesta mitochondria: from genome to proteome

Background

The ascomycete fungus Ceratocystis cacaofunesta is the causal agent of wilt disease in cacao, which results in significant economic losses in the affected producing areas. Despite the economic importance of the Ceratocystis complex of species, no genomic data are available for any of its members. Given that mitochondria play important roles in fungal virulence and the susceptibility/resistance of fungi to fungicides, we performed the first functional analysis of this organelle in Ceratocystis using integrated “omics” approaches.

Results

The C. cacaofunesta mitochondrial genome (mtDNA) consists of a single, 103,147-bp circular molecule, making this the second largest mtDNA among the Sordariomycetes. Bioinformatics analysis revealed the presence of 15 conserved genes and 37 intronic open reading frames in C. cacaofunesta mtDNA. Here, we predicted the mitochondrial proteome (mtProt) of C. cacaofunesta, which is comprised of 1,124 polypeptides - 52 proteins that are mitochondrially encoded and 1,072 that are nuclearly encoded. Transcriptome analysis revealed 33 probable novel genes. Comparisons among the Gene Ontology results of the predicted mtProt of C. cacaofunesta, Neurospora crassa and Saccharomyces cerevisiae revealed no significant differences. Moreover, C. cacaofunesta mitochondria were isolated, and the mtProt was subjected to mass spectrometric analysis. The experimental proteome validated 27% of the predicted mtProt. Our results confirmed the existence of 110 hypothetical proteins and 7 novel proteins of which 83 and 1, respectively, had putative mitochondrial localization.

Conclusions

The present study provides the first partial genomic analysis of a species of the Ceratocystis genus and the first predicted mitochondrial protein inventory of a phytopathogenic fungus. In addition to the known mitochondrial role in pathogenicity, our results demonstrated that the global function analysis of this organelle is similar in pathogenic and non-pathogenic fungi, suggesting that its relevance in the lifestyle of these organisms should be based on a small number of specific proteins and/or with respect to differential gene regulation. In this regard, particular interest should be directed towards mitochondrial proteins with unknown function and the novel protein that might be specific to this species. Further functional characterization of these proteins could enhance our understanding of the role of mitochondria in phytopathogenicity.

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

Background

Mitochondrial (mt) markers are successfully applied in evolutionary biology and systematics because mt genomes often evolve faster than the nuclear genomes. In addition, they allow robust phylogenetic analysis based on conserved proteins of the oxidative phosphorylation system. In the present study we sequenced and annotated the complete mt genome of P. subalpina, a member of the Phialocephala fortinii s.l. – Acephala applanata species complex (PAC). PAC belongs to the Helotiales, which is one of the most diverse groups of ascomycetes including more than 2,000 species. The gene order was compared to deduce the mt genome evolution in the Pezizomycotina. Genetic variation in coding and intergenic regions of the mtDNA was studied for PAC to assess the usefulness of mt DNA for species diagnosis.

Results

The mt genome of P. subalpina is 43,742 bp long and codes for 14 mt genes associated with the oxidative phosphorylation. In addition, a GIY-YIG endonuclease, the ribosomal protein S3 (Rps3) and a putative N-acetyl-transferase were recognized. A complete set of tRNA genes as well as the large and small rRNA genes but no introns were found. All protein-coding genes were confirmed by EST sequences. The gene order in P. subalpina deviated from the gene order in Sclerotinia sclerotiorum, the only other helotialean species with a fully sequenced and annotated mt genome. Gene order analysis within Pezizomycotina suggests that the evolution of gene orders is mostly driven by transpositions. Furthermore, sequence diversity in coding and non-coding mtDNA regions in seven additional PAC species was pronounced and allowed for unequivocal species diagnosis in PAC.

Conclusions

The combination of non-interrupted ORFs and EST sequences resulted in a high quality annotation of the mt genome of P. subalpina, which can be used as a reference for the annotation of other mt genomes in the Helotiales. In addition, our analyses show that mtDNA loci will be the marker of choice for future analysis of PAC communities.

Complete mitochondrial genomes of two cockroaches, Blattella germanica and Periplaneta americana, and the phylogenetic position of termites

The mitochondrial genomes are one of the most information-rich markers in phylogenetics. The relationships within superorder Dictyoptera have been debated in the literature. However, the closely related termites (Isoptera) are retained as unranked taxon within the order Blattaria (cockroaches). In this work, we sequenced the complete mitogenomes of two cockroaches, reconstructed the molecular phylogeny and attempted to infer the phylogenetic position of termites in Blattaria more reliably. The complete mtDNA nucleotide sequences of the peridomestic American cockroach (Periplaneta americana L.) and the domestic German cockroach (Blattella germanica L.) are 15,025 and 15,584 bp in size, respectively. The genome shares the gene order and orientation with previously known Blattaria mitogenomes. Most tRNAs could be folded into the typical cloverleaf secondary structure, but the tRNA-Ser (AGN) of P. americana appears to be missing the dihydrouridine arm. Using nucleotide and amino acid sequences as phylogenetic markers, we proposed that termites should be treated as a superfamily (Termitoidea) of cockroaches. We suggested that Polyphagoidea was the sister group of Termitoidea in Blattaria and supported that the suborder Caelifera is more closely related to the Phasmatodea than to the suborder Ensifera of Orthoptera.

Mitochondrial haplotype analysis as a tool for differentiating isolates of Verticillium dahliae

The ability to monitor mitochondrial background in Verticillium dahliae may provide an additional tool for population studies and monitoring clonal populations. Published mitochondrial genome sequences of V. dahliae (DQ351941) were used to design primers for amplification of spacer regions for assessment of mitochondrial haplotype differences among isolates. Five regions were examined (5,229 bp, or 19% of the total genome size) for 30 isolates representing a range in vegetative compatibility group (VCG), host, and geographic origin. Observed differences among isolates were due to single nucleotide polymorphisms, different numbers of bases in specific homopolymeric regions, and copies of subrepeated sequences. When the differences observed for each locus were totaled there were 28 total groupings; when the results of each locus for individual isolates were combined there were 15 mitochondrial haplotypes. Some of the observed groupings correlated with VCG. For example, five VCG-1A and VCG-1B isolates from California, Spain, and Greece had identical haplotypes; however, this was not observed for VCG-2 or VCG-4 isolates. While some VCG-2 isolates also were identical and fell into a single haplotype, five haplotypes were found for this group (five other haplotypes were observed for other isolates that had not been characterized to VCG but grouped with VCG-2 isolates in the phylogenetic analysis). Likewise, five VCG-4 isolates fell into four mitochondrial haplotypes, one of which was identical to the largest VCG-2 grouping. A heterokaryon self-incompatible isolate that was reported in the literature to cluster with VCG-2 isolates by amplified fragment length polymorphism analysis was identical with VCG-1 isolates for four of the five loci, but was intermediate between VCG-1 and VCG-2 in the haplotype analysis. Phylogenetic analysis with these regions revealed the mitochondrial background of VCG-1 and VCG2-B to be monophyletic but VCG-2A and VCG-4 could not be separated. The results obtained indicate that there is variation in mitochondrial haplotypes and this type of analysis may be a useful for characterization of isolates. While data from all five regions was used for the haplotype separation in this study, depending on the VCG or the level of variability observed within a population it is possible to use fewer loci.

Phylogenetic and biogeographic implications inferred by mitochondrial intergenic region analyses and ITS1-5.8S-ITS2 of the entomopathogenic fungi Beauveria bassiana and B. brongniartii

BACKGROUND

The entomopathogenic fungi of the genus Beauveria are cosmopolitan with a variety of different insect hosts. The two most important species, B. bassiana and B. brongniartii, have already been used as biological control agents of pests in agriculture and as models for the study of insect host - pathogen interactions. Mitochondrial (mt) genomes, due to their properties to evolve faster than the nuclear DNA, to contain introns and mobile elements and to exhibit extended polymorphisms, are ideal tools to examine genetic diversity within fungal populations and genetically identify a species or a particular isolate. Moreover, mt intergenic region can provide valuable phylogenetic information to study the biogeography of the fungus.

RESULTS

The complete mt genomes of B. bassiana (32,263 bp) and B. brongniartii (33,920 bp) were fully analysed. Apart from a typical gene content and organization, the Beauveria mt genomes contained several introns and had longer intergenic regions when compared with their close relatives. The phylogenetic diversity of a population of 84 Beauveria strains -mainly B. bassiana (n = 76) - isolated from temperate, sub-tropical and tropical habitats was examined by analyzing the nucleotide sequences of two mt intergenic regions (atp6-rns and nad3-atp9) and the nuclear ITS1-5.8S-ITS2 domain. Mt sequences allowed better differentiation of strains than the ITS region. Based on mt and the concatenated dataset of all genes, the B. bassiana strains were placed into two main clades: (a) the B. bassiana s. l. and (b) the "pseudobassiana". The combination of molecular phylogeny with criteria of geographic and climatic origin showed for the first time in entomopathogenic fungi, that the B. bassiana s. l. can be subdivided into seven clusters with common climate characteristics.

CONCLUSIONS

This study indicates that mt genomes and in particular intergenic regions provide molecular phylogeny tools that combined with criteria of geographic and climatic origin can subdivide the B. bassiana s.l. entomopathogenic fungi into seven clusters with common climate characteristics.

Streptomyces lividans inhibits the proliferation of the fungus Verticillium dahliae on seeds and roots of Arabidopsis thaliana

Verticillium wilt, a vascular disease in more than 200 dicotyledonous plants, is due to the ascomycete fungus Verticillium dahliae. As documented by video-microscopy, the soil bacterium Streptomyces lividans strongly reduces the germination of V. dahliae conidia, and the subsequent growth of hyphae. Quantification by the use of DNA-intercalating dyes and Calcofluor-staining revealed that during prolonged co-cultivation, bacterial hyphae proliferate to a dense network, provoke a poor development of V. dahliae vegetative hyphae and lead to an enormous reduction of conidia and microsclerotia. Upon individual application to seeds of the model plant Arabidopsis thaliana, either the bacterial spores or the fungal conidia germinate at or within the mucilage, including its volcano-shaped structures. The extension of hyphae from each individual strain correlates with the reduction of the pectin-containing mucilage-layer. Proliferating hyphae then spread to roots of the emerging seedlings. Plants, which arise in the presence of V. dahliae within agar or soil, have damaged root cells, an atrophied stem and root, as well as poorly developed leaves with chlorosis symptoms. In contrast, S. lividans hyphae settle in bunches preferentially at the outer layer near tips and alongside roots. Resulting plants have a healthy appearance including an intact root system. Arabidopsis thaliana seeds, which are co-inoculated with V. dahliae and S. lividans, have preferentially proliferating bacterial hyphae within the mucilage, and at roots of the outgrowing seedlings. As a result, plants have considerably reduced disease symptoms. As spores of the beneficial S. lividans strain are obtainable in large quantity, its application is highly attractive.

The Cryphonectria parasitica mitochondrial rns gene: plasmid-like elements, introns and homing endonucleases

The mt-rns gene of Cryphonectria parasitica is 9872bp long and includes two group I and two group II introns. An analysis of intronic protein-encoding sequences revealed that LAGLIDADG ORFs, which usually are associated with group I introns, were transferred at least twice into group II introns. A plasmid-like mitochondrial element (plME) that appears in high amounts in previously mutagen-induced mit1 and mit2 hypovirulent mutants of the Ep155 standard virulent strain of C. parasitica was found to be derived from a short region of the mt-rns gene, including the exon 1 and most of the first intron. The plME is a 4.2-kb circular, multimeric DNA and an autonomously-replicating mtDNA fragment. Although sexual transmission experiments indicate that the plME does not directly cause hypovirulence, its emergence is one manifestation of the many complex molecular and genetic events that appear to underlie this phenotype.

Mitochondrial gene sequences alone or combined with ITS region sequences provide firm molecular criteria for the classification of Lecanicillium species

The recent revision of Verticillium sect. Prostrata led to the introduction of the genus Lecanicillium, which comprises the majority of the entomopathogenic strains. Sixty-five strains previously classified as Verticillium lecanii or Verticillium sp. from different geographical regions and hosts were examined and their phylogenetic relationships were determined using sequences from three mitochondrial (mt) genes [the small rRNA subunit (rns), the NADH dehydrogenase subunits 1 (nad1) and 3 (nad3)] and the ITS region. In general, single gene phylogenetic trees differentiated and placed the strains examined in well-supported (by BS analysis) groups of L. lecanii, L. longisporum, L. muscarium, and L. nodulosum, although in some cases a few uncertainties still remained. nad1 was the most informative single gene in phylogenetic analyses and was also found to contain group I introns with putative open reading frames (ORFs) encoding for GIY-YIG endonucleases. The combined use of mt gene sequences resolved taxonomic uncertainties arisen from ITS analysis and, alone or in combination with ITS sequences, helped in placing uncharacterised Verticillium lecanii and Verticillium sp. firmly into Lecanicillium species. Combined gene data from all the mt genes and all the mt genes and the ITS region together, were very similar. Furthermore, a relaxed correlation with host specificity -- at least for Homoptera -- was indicated for the rns and the combined mt gene sequences. Thus, the usefulness of mt gene sequences as a convenient molecular tool in phylogenetic studies of entomopathogenic fungi was demonstrated.

Molecular characterization of isolates of Beauveria bassiana obtained from overwintering and summer populations of Sunn Pest (Eurygaster integriceps)

AIM

To examine whether isolates of the entomopathogenic fungus Beauveria bassiana are more closely associated to their summer hosts compared with overwintering hosts, with recently developed molecular tools based on mitochondrial regions.

METHODS AND RESULTS

Primers for the traditional ITS1-5.8S-ITS2 region and two mitochondrial intergenic regions, namely, nad3-atp9 and atp6-rns, were used. All amplified products were sequenced, aligned and Neighbour-Joining (NJ), parsimony and Bayesian phylogenetic inference analyses were performed. The isolates examined were grouped with very good support into three distinct groups, two of them showed geographical correlation, but no clear association to their host.

CONCLUSIONS

The mitochondrial intergenic regions used were more informative than the nuclear ITS1-5.8S-ITS2 sequences. The sequence variability observed, that allowed the phylogenetic placement of the isolates into distinct groups, depended on the geographical origin of the isolates and can be exploited for designing group-specific and isolate-specific primers for their genetic fingerprinting. No clear associations with summer Sunn Pest populations were observed.

SIGNIFICANCE AND IMPACT OF THE STUDY

Studies on the genetic variability of biocontrol agents like B. bassiana are indispensable for the development of molecular tools for their future monitoring.

The sporadic occurrence of a group I intron-like element in the mtDNA rnl gene of Ophiostoma novo-ulmi subsp. americana

The presence of group I intron-like elements within the U7 region of the mtDNA large ribosomal subunit RNA gene (rnl) was investigated in strains of Ophiostoma novo-ulmi subsp. americana from Canada, Europe and Eurasia, and in selected strains of O. ips, O. minus, O. piceae, O. ulmi, and O. himal-ulmi. This insertion is of interest as it has been linked previously to the generation of plasmid-like mtDNA elements in diseased strains of O. novo-ulmi. Among 197 O. novo-ulmi subsp. americana strains tested, 61 contained a 1.6kb insertion within the rnl-U7 region and DNA sequence analysis suggests the presence of a group I intron (IA1 type) that encodes a potential double motif LAGLIDADG homing endonuclease-like gene (HEG). Phylogenetic analysis of rnl-U7 intron encoded HEG-like elements supports the view that double motif HEGs originated from a duplication event of a single-motif HEG followed by a fusion event that combined the two copies into one open reading frame (ORF). The data also show that rnl-U7 intron encoded ORFs belong to a clade that includes ORFs inserted into different types of group I introns, e.g. IB, ID, IC3, IA1, present within a variety of different mtDNA genes, such as the small ribosomal subunit RNA gene (rns), apo-cytochrome b gene (cob), NADH dehydrogenase subunit 5 (nad5), cytochrome oxidase subunit 1 gene (coxI), and ATPase subunit 9 gene (atp9). We also compared the occurrence of the rnl-U7 intron in our collection of 227 strains with the presence of the rnl-U11 group I intron and concluded that the U7 intron appears to be an optional element and the U11 intron is probably essential among the strains tested.

Complete mitochondrial genome sequence of the wine yeast Candida zemplinina: intraspecies distribution of a novel group-IIB1 intron with eubacterial affiliations

The mtDNA of the ascomycetous wine yeast Candida zemplinina is a circularly mapping genome of 23,114 bp. It contains 35 genes coding for the seven basic subunits of oxidative phosporylation found in yeasts (the genes encoding for NADH oxidoreductase subunits are absent), the ribosomal protein Var1, two rRNAs and 25 tRNA genes. Although protein phylogenetic analysis showed a divergent mitochondrial genome, several traits appeared preserved. The conserved gene blocks between the mtDNAs of C. zemplinina and Candida glabrata were maintained and changes in gene order and putative promoters were due to restricted genome reshuffling. New heterogeneous hairpin elements were identified scattered throughout cox1 introns. The large subunit rRNA gene harboured the first group-IIB1 intron containing a putative active reverse transcriptase (RT) in mitochondrial genomes of fungi. Phylogenetic analysis of the RT protein confirmed its closer relationship to eubacterial intronic RTs, while being only distantly related to all other fungal mitochondrial group-II introns and RTs. The findings point towards an early migration event of a eubacterial group-II intron to the mitochondrial genome of C. zemplinina.