The complete mitochondrial genome of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae: gene order and trn gene clusters reveal a common evolutionary course for all Sordariomycetes, while intergenic regions show variation

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.

Trichoderma koningiopsis (Hypocreaceae) has the smallest mitogenome of the genus Trichoderma

Introduction

Fungal mitogenomes exhibit remarkable variation in conformation, size, gene content, arrangement and expression, including their intergenic spacers and introns.

Methods

The complete mitochondrial genome sequence of the mycoparasitic fungus Trichoderma koningiopsis was determined using the Illumina next-generation sequencing technology. We used data from our recent Illumina NGS-based project of T. koningiopsis genome sequencing to study its mitochondrial genome. The mitogenome was assembled, annotated, and compared with other fungal mitogenomes.

Results

T. koningiopsis strain POS7 mitogenome is a circular molecule of 27,560 bp long with a GC content of 27.80%. It harbors the whole complement of the 14 conserved mitochondrial protein-coding genes (PCG) such as atp6, atp8, atp9, cox1, cox2, cox3, cob, nad1, nad2, nad3, nad4, nad4L, nad5, and nad6, also found in the same gene order to other Hypocreales. The mitogenome also contains 26 transfer RNA genes (tRNAs), 5 of them with more than one copy. Other genes also present in the assembled mitochondrial genome are a small rRNA subunit and a large rRNA subunit containing ribosomal protein S3 gene. Despite the small genome size, two introns were detected in the T. koningiopsis POS7 mitogenome, one of them in cox3 gene and the other in rnl gene, accounting 7.34% of this mitogenome with a total size of 2,024 bp. A phylogenetic analysis was done using the 14 PCGs genes of T. koningiopsis strain POS7 mitogenome to compare them with those from other fungi of the Subphyla Pezizomycotina and Saccharomycotina. T. koningiopsis strain POS7 was clustered together with other representatives of Trichoderma lineage, within the Hypocreales group, which is also supported by previous phylogenetic studies based on nuclear markers.

Discussion

The mitochondrial genome of T. koningiopsis POS7 will allow further investigations into the taxonomy, phylogenetics, conservation genetics, and evolutionary biology of this important genus as well as other closely related species.

Mitochondrial Transcription of Entomopathogenic Fungi Reveals Evolutionary Aspects of Mitogenomes

Entomopathogenic fungi and more specifically genera Beauveria and Metarhizium have been exploited for the biological control of pests. Genome analyses are important to understand better their mode of action and thus, improve their efficacy against their hosts. Until now, the sequences of their mitochondrial genomes were studied, but not at the level of transcription. Except of yeasts and Neurospora crassa, whose mt gene transcription is well described, in all other Ascomycota, i.e., Pezizomycotina, related information is extremely scarce. In this work, mt transcription and key enzymes of this function were studied. RT-PCR experiments and Northern hybridizations reveal the transcriptional map of the mt genomes of B. bassiana and M. brunneum species. The mt genes are transcribed in six main transcripts and undergo post-transcriptional modifications to create single gene transcripts. Promoters were determined in both mt genomes with a comparative in silico analysis, including all known information from other fungal mt genomes. The promoter consensus sequence is 5'-ATAGTTATTAT-3' which is in accordance with the definition of the polycistronic transcripts determined with the experiments described above. Moreover, 5'-RACE experiments in the case of premature polycistronic transcript nad1-nad4-atp8-atp6 revealed the 5' end of the RNA transcript immediately after the in silico determined promoter, as also found in other fungal species. Since several conserved elements were retrieved from these analyses compared to the already known data from yeasts and N. crassa, the phylogenetic analyses of mt RNA polymerase (Rpo41) and its transcriptional factor (Mtf1) were performed in order to define their evolution. As expected, it was found that fungal Rpo41 originate from the respective polymerase of T7/T3 phages, while the ancestor of Mtf1 is of alpha-proteobacterial origin. Therefore, this study presents insights about the fidelity of the mt single-subunit phage-like RNA polymerase during transcription, since the correct identification of mt promoters from Rpo41 requires an ortholog to bacterial sigma factor, i.e., Mtf1. Thus, a previously proposed hypothesis of a phage infected alpha-proteobacterium as the endosymbiotic progenitor of mitochondrion is confirmed in this study and further upgraded by the co-evolution of the bacterial (Mtf1) and viral (Rpo41) originated components in one functional unit.

Comparative analysis of Malassezia furfur mitogenomes and the development of a mitochondria-based typing approach

Malassezia furfur is a yeast species belonging to Malasseziomycetes, Ustilaginomycotina and Basidiomycota that is found on healthy warm-blooded animal skin, but also involved in various skin disorders like seborrheic dermatitis/dandruff and pityriasis versicolor. Moreover, Malassezia are associated with bloodstream infections, Crohn's disease and pancreatic carcinoma. Recent advances in Malassezia genomics and genetics have focused on the nuclear genome. In this work, we present the M. furfur mitochondrial (mt) genetic heterogenicity with full analysis of 14 novel and six available M. furfur mt genomes. The mitogenome analysis reveals a mt gene content typical for fungi, including identification of variable mt regions suitable for intra-species discrimination. Three of them, namely the trnK–atp6 and cox3–nad3 intergenic regions and intron 2 of the cob gene, were selected for primer design to identify strain differences. Malassezia furfur strains belonging to known genetic variable clusters, based on AFLP and nuclear loci, were assessed for their mt variation using PCR amplification and sequencing. The results suggest that these mt regions are excellent molecular markers for the typing of M. furfur strains and may provide added value to nuclear regions when assessing evolutionary relationships at the intraspecies level.

Comparative Mitogenomic Analysis and the Evolution of Rhizoctonia solani Anastomosis Groups

Mitochondria are the major energy source for cell functions. However, for the plant fungal pathogens, mitogenome variations and their roles during the host infection processes remain largely unknown. Rhizoctonia solani, an important soil-borne pathogen, forms different anastomosis groups (AGs) and adapts to a broad range of hosts in nature. Here, we reported three complete mitogenomes of AG1-IA RSIA1, AG1-IB RSIB1, and AG1-IC, and performed a comparative analysis with nine published Rhizoctonia mitogenomes (AG1-IA XN, AG1-IB 7/3/14, AG3, AG4, and five Rhizoctonia sp. mitogenomes). These mitogenomes encoded 15 typical proteins (cox1-3, cob, atp6, atp8-9, nad1-6, nad4L, and rps3) and several LAGLIDADG/GIY-YIG endonucleases with sizes ranging from 109,017 bp (Rhizoctonia sp. SM) to 235,849 bp (AG3). We found that their large sizes were mainly contributed by repeat sequences and genes encoding endonucleases. We identified the complete sequence of the rps3 gene in 10 Rhizoctonia mitogenomes, which contained 14 positively selected sites. Moreover, we inferred a robust maximum-likelihood phylogeny of 32 Basidiomycota mitogenomes, representing that seven R. solani and other five Rhizoctonia sp. lineages formed two parallel branches in Agaricomycotina. The comparative analysis showed that mitogenomes of Basidiomycota pathogens had high GC content and mitogenomes of R. solani had high repeat content. Compared to other strains, the AG1-IC strain had low substitution rates, which may affect its mitochondrial phylogenetic placement in the R. solani clade. Additionally, with the published RNA-seq data, we investigated gene expression patterns from different AGs during host infection stages. The expressed genes from AG1-IA (host: rice) and AG3 (host: potato) mainly formed four groups by k-mean partitioning analysis. However, conserved genes represented varied expression patterns, and only the patterns of rps3-nad2 and nad1-m3g18/mag28 (an LAGLIDADG endonuclease) were conserved in AG1-IA and AG3 as shown by the correlation coefficient analysis, suggesting regulation of gene repertoires adapting to infect varied hosts. The results of variations in mitogenome characteristics and the gene substitution rates and expression patterns may provide insights into the evolution of R. solani mitogenomes.

Mitochondrial evolution in the entomopathogenic fungal genus Beauveria

Species in the fungal genus Beauveria are pathogens of invertebrates and have been commonly used as the active agent in biopesticides. After many decades with few species described, recent molecular approaches to classification have led to over 25 species now delimited. Little attention has been given to the mitochondrial genomes of Beauveria but better understanding may led to insights into the nature of species and evolution in this important genus. In this study, we sequenced the mitochondrial genomes of four new strains belonging to Beauveria bassiana, Beauveria caledonica and Beauveria malawiensis, and compared them to existing mitochondrial sequences of related fungi. The mitochondrial genomes of Beauveria ranged widely from 28,806 to 44,135 base pairs, with intron insertions accounting for most size variation and up to 39% (B. malawiensis) of the mitochondrial length due to introns in genes. Gene order of the common mitochondrial genes did not vary among the Beauveria sequences, but variation was observed in the number of transfer ribonucleic acid genes. Although phylogenetic analysis using whole mitochondrial genomes showed, unsurprisingly, that B. bassiana isolates were the most closely related to each other, mitochondrial codon usage suggested that some B. bassiana isolates were more similar to B. malawiensis and B. caledonica than the other B. bassiana isolates analyzed.

Mitogenome of Tolypocladium guangdongense

Tolypocladium guangdongense is a high-value edible fungus with various medicinal and food safety properties. However, its evolutionary and genetic information is still limited. Mitochondrial genomes are potential models for molecular evolution and phylogenetic studies. In this study, we sequenced the complete mitogenome of T. guangdongense, demonstrating circular sequence of 46,102 bp, containing 14 standard protein-coding genes (PCGs), 2 ribosomal RNA subunit genes, and 28 tRNA genes. Phylogenetic analysis based on mitochondrial genes indicated that T. guangdongense was clustered into the Tolypocladium genus with high support value, based on the core PCG dataset. In addition, rps3 is also a suitable marker in the phylogenetic analysis in Hypocreales. Gene rearrangement analysis indicated that the gene order of PCGs was highly consistent in Hypocreales, and tRNA rearrangement events occurred in most species of Hypocreales; however, the rearrangement rates were not taxonomically correlated. Divergence time estimation based on the old fossil record and previous reports revealed that T. guangdongense originated approximately in the middle Cenozoic (42 Mya, 95% highest posterior density interval: 43-116) with the Tolypocladium genus differentiation. Our results provided more mitogenomic information of T. guangdongense and shed new insights into evolution of the Tolypocladium genus. KEY POINTS: • The general and unique features of T. guangdongense mitogenome are firstly reported. • Phylogenetic analysis further verified the taxonomic status of T. guangdongense. • Divergence time estimation provides more evolutionary information of T. guangdongense.

The mitochondrial genome contribution to the phylogeny and identification of Metarhizium species and strains

The genus Metarhizium is composed of entomopathogenic fungal biological control agents (BCAs) used for invertebrate pest control. The phylogenetic relationships of species within this genus are still under scrutiny as several cryptic species can be found. In this work, the mitochondrial (mt) genome of Metarhizium brunneum ARSEF 4556 was fully sequenced and a comparative genome analysis was conducted with 7 other available mt genomes, belonging to 5 Metarhizium species: M. anisopliae, M. brunneum, M. robertsii, M. guizhouense and M. majus. Results showed that Metarhizium demonstrates greater conserved stability than other fungal mt genomes. Furthermore, this analysis located 7 diverse regions in both intergenic domains and gene fragments which were ideal for species/strain discrimination. The sequencing of these regions revealed several SNPs among 38 strains tested, 11 of which were uncharacterized. Single gene phylogenies presented variable results which may be used further for intra-species discrimination. Phylogenetic trees based on the concatenation of mt domains and the nuclear ITS1-5.8S-ITS2 region showed discrimination of the species studied and allowed the identification of uncharacterized strains. These were mostly placed within species M. anisopliae and M. brunneum. Five strains clustered together in a clade related to M. brunneum, suggesting that they comprise a cryptic species.

Detecting Introgression Between Members of the Fusarium fujikuroi and F. oxysporum Species Complexes by Comparative Mitogenomics

The Fusarium fujikuroi species complex (FFSC) and F. oxysporum species complex (FOSC) are two related groups of plant pathogens causing a wide diversity of diseases in agricultural crops world wide. The aims of this study are (1) to clarify the phylogeny of the FFSC, (2) to identify potential deviation from tree-like evolution, (3) to explore the value of using mitogenomes for these kinds of analyses, and (4) to better understand mitogenome evolution. In total, we have sequenced 24 species from the FFSC and a representative set of recently analyzed FOSC strains was chosen, while F. redolens was used as outgroup for the two species complexes. A species tree was constructed based on the concatenated alignment of seven nuclear genes and the mitogenome, which was contrasted to individual gene trees to identify potential conflicts. These comparisons indicated conflicts especially within the previously described African clade of the FFSC. Furthermore, the analysis of the mitogenomes revealed the presence of a variant of the large variable (LV) region in FFSC which was previously only reported for FOSC. The distribution of this variant and the results of sequence comparisons indicate horizontal genetic transfer between members of the two species complexes, most probably through introgression. In addition, a duplication of atp9 was found inside an intron of cob, which suggests that even highly conserved mitochondrial genes can have paralogs. Paralogization in turn may lead to inaccurate single gene phylogenies. In conclusion, mitochondrial genomes provide a robust basis for phylogeny. Comparative phylogenetic analysis indicated that gene flow among and between members of FFSC and FOSC has played an important role in the evolutionary history of these two groups. Since mitogenomes show greater levels of conservation and synteny than nuclear regions, they are more likely to be compatible for recombination than nuclear regions. Therefore, mitogenomes can be used as indicators to detect interspecies gene flow.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Mt-rps3 is an ancient gene which provides insight into the evolution of fungal mitochondrial genomes

The nuclear ribosomal protein S3 (Rps3) is implicated in the assembly of the ribosomal small subunit. Fungi and plants present a gene copy in their mitochondrial (mt) genomes. An analysis of 303 complete fungal mt genomes showed that, when rps3 is found, it is either a free-standing gene or an anchored gene within the omega intron of the rnl gene. Early divergent fungi, Basidiomycota and all yeasts but the CTG group belong to the first case, and Pezizomycotina to the second. Its position, size and genetic code employed are conserved within species of the same Order. Size variability is attributed to different number of repeats. These repeats consist of AT-rich sequences. MtRps3 proteins lack the KH domain, necessary for binding to rRNA, in their N-terminal region. Their C-terminal region is conserved in all Domains of life. Phylogenetic analysis showed that nuclear and mtRps3 proteins are descendants of archaeal and a-proteobacterial homologues, respectively. Thus, fungal mt-rps3 gene is an ancient gene which evolved within the endosymbiotic model and presents different evolutionary routes: (a) coming from a-proteobacteria, it was relocated to another region of the mt genome, (b) via its insertion to the omega intron, it was transferred to the nucleus and/or got lost, and (c) it was re-routed to the mt genome again. Today, Basidiomycota and Saccharomycetales seem to follow the first evolutionary route and almost all Pezizomycotina support the second scenario with their exceptions being the result of the third scenario, i.e., the gene's re-entry to the mt genome.

Comparative analysis of the mitochondrial genome of the fungus Colletotrichum lindemuthianum, the causal agent of anthracnose in common beans

Fungi of the genus Colletotrichum are economically important and are used as models in plant-pathogen interaction studies. In this study, the complete mitochondrial genomes of two Colletotrichum lindemuthianum isolates were sequenced and compared with the mitochondrial genomes of seven species of Colletotrichum. The mitochondrial genome of C. lindemuthianum is a typical circular molecule 37,446 bp (isolate 89 A₂ 2-3) and 37,440 bp (isolate 83.501) in length. The difference of six nucleotides between the two genomes is the result of a deletion in the ribosomal protein S3 (rps3) gene in the 83.501 isolate. In addition, substitution of adenine for guanine within the rps3 gene in the mitochondrial genome of the 83.501 isolate was observed. Compared to the previously sequenced C. lindemuthianum mitochondrial genome, an exon no annotated in the cytochrome c oxidase I (cox1) gene and a non-conserved open reading frame (ncORF) were observed. The size of the mitochondrial genomes of the seven species of Colletotrichum was highly variable, being attributed mainly to the ncORF, ranging from one to 10 and also from introns ranging from one to 11 and which encode a total of up to nine homing endonucleases. This paper reports for the first time by means of transcriptome that then ncORFs are transcribed in Colletotrichum spp. Phylogeny data revealed that core mitochondrial genes could be used as an alternative in phylogenetic relationship studies in Colletotrichum spp. This work contributes to the genetic and biological knowledge of Colletotrichum spp., which is of great economic and scientific importance.

The complete mitochondrial genome of Engyodontium album and comparative analyses with Ascomycota mitogenomes

Engyodontium album is a widespread pathogen that causes different kinds of dermatoses and respiratory tract diseases in humans and animals. In spite of its perniciousness, the basic genetic and molecular background of this species remains poorly understood. In this study, the mitochondrial genome sequence of E. album was determined using a high-throughput sequencing platform. The circular mitogenome was found to be 28,081 nucleotides in length and comprised of 17 protein-coding genes, 24 tRNA genes, and 2 rRNA genes. The nucleotide composition of the genome was A+T-biased (74.13%). Group-II introns were found in the nad1, nad5, and cob genes. The most frequently used codon of protein-coding genes was UAU. Isoleucine was identified as the most common amino acid, while proline was the least common amino acid in protein-coding genes. The gene-arrangement order is nearly the same when compared with other Ascomycota mitogenomes. Phylogenetic relationships based on the shared protein-coding genes revealed that E. album is closely related to the Cordycipitaceae family, with a high-confidence support value (100%). The availability of the mitogenome of E. album will shed light on the molecular systematic and genetic differentiation of this species.

Genomic analyses reveal low mitochondrial and high nuclear diversity in the cyclosporin-producing fungus Tolypocladium inflatum

Mitochondrial DNA is generally regarded to evolve faster than nuclear DNA in animals, whereas if this is also true in fungi remains unclear. Herein, we annotate the first complete mitochondrial genome (mitogenome) of the cyclosporin-producing fungus Tolypocladium inflatum and report the genome-wide sequence variations among five isolates originating from distantly separated localities. We found that T. inflatum has among the most compact of fungal mitogenomes; its 25 kb DNA molecule encodes all standard fungal mitochondrial genes and harbors only one intron. Transcriptional analyses validated the expression of most conserved genes. We found several uncommon repetitive elements and evidence of gene transfer from the mitochondrion to the nucleus. Phylogenetic analyses confirmed the placement of T. inflatum in the fungal order Hypocreales although there was uncertainty on its family-level affiliation. Comparative genomic analyses among the five isolates identified an overall lower level of intraspecific variation in mitogenomes than in nuclear genomes; however, both the nuclear and mitochondrial genomes revealed similar isolate relationships, not correlating with geographic sources of these isolates. Our study shed new insights into the evolution of the medicinally important ascomycete T. inflatum.

van der Lee TA, Waalwijk C, van Diepeningen AD. Mitochondrial genomes reveal recombination in the presumed asexual Fusarium oxysporum species complex

BACKGROUND

The Fusarium oxysporum species complex (FOSC) contains several phylogenetic lineages. Phylogenetic studies identified two to three major clades within the FOSC. The mitochondrial sequences are highly informative phylogenetic markers, but have been mostly neglected due to technical difficulties.

RESULTS

A total of 61 complete mitogenomes of FOSC strains were de novo assembled and annotated. Length variations and intron patterns support the separation of three phylogenetic species. The variable region of the mitogenome that is typical for the genus Fusarium shows two new variants in the FOSC. The variant typical for Fusarium is found in members of all three clades, while variant 2 is found in clades 2 and 3 and variant 3 only in clade 2. The extended set of loci analyzed using a new implementation of the genealogical concordance species recognition method support the identification of three phylogenetic species within the FOSC. Comparative analysis of the mitogenomes in the FOSC revealed ongoing mitochondrial recombination within, but not between phylogenetic species.

CONCLUSIONS

The recombination indicates the presence of a parasexual cycle in F. oxysporum. The obstacles hindering the usage of the mitogenomes are resolved by using next generation sequencing and selective genome assemblers, such as GRAbB. Complete mitogenome sequences offer a stable basis and reference point for phylogenetic and population genetic studies.

SMRT Sequencing Revealed Mitogenome Characteristics and Mitogenome-Wide DNA Modification Pattern in Ophiocordyceps sinensis

Single molecule, real-time (SMRT) sequencing was used to characterize mitochondrial (mt) genome of Ophiocordyceps sinensis and to analyze the mt genome-wide pattern of epigenetic DNA modification. The complete mt genome of O. sinensis, with a size of 157,539 bp, is the fourth largest Ascomycota mt genome sequenced to date. It contained 14 conserved protein-coding genes (PCGs), 1 intronic protein rps3, 27 tRNAs and 2 rRNA subunits, which are common characteristics of the known mt genomes in Hypocreales. A phylogenetic tree inferred from 14 PCGs in Pezizomycotina fungi supports O. sinensis as most closely related to Hirsutella rhossiliensis in Ophiocordycipitaceae. A total of 36 sequence sites in rps3 were under positive selection, with dN/dS >1 in the 20 compared fungi. Among them, 16 sites were statistically significant. In addition, the mt genome-wide base modification pattern of O. sinensis was determined in this study, especially DNA methylation. The methylations were located in coding and uncoding regions of mt PCGs in O. sinensis, and might be closely related to the expression of PCGs or the binding affinity of transcription factor A to mtDNA. Consequently, these methylations may affect the enzymatic activity of oxidative phosphorylation and then the mt respiratory rate; or they may influence mt biogenesis. Therefore, methylations in the mitogenome of O. sinensis might be a genetic feature to adapt to the cold and low PO2 environment at high altitude, where O. sinensis is endemic. This is the first report on epigenetic modifications in a fungal mt genome.

Genomic characterization of two novel pathogenic avipoxviruses isolated from pacific shearwaters (Ardenna spp.)

Background

Over the past 20 years, many marine seabird populations have been gradually declining and the factors driving this ongoing deterioration are not always well understood. Avipoxvirus infections have been found in a wide range of bird species worldwide, however, very little is known about the disease ecology of avian poxviruses in seabirds. Here we present two novel avipoxviruses from pacific shearwaters (Ardenna spp), one from a Flesh-footed Shearwater (A. carneipes) (SWPV-1) and the other from a Wedge-tailed Shearwater (A. pacificus) (SWPV-2).

Results

Epidermal pox lesions, liver, and blood samples were examined from A. carneipes and A. pacificus of breeding colonies in eastern Australia. After histopathological confirmation of the disease, PCR screening was conducted for avipoxvirus, circovirus, reticuloendotheliosis virus, and fungal agents. Two samples that were PCR positive for poxvirus were further assessed by next generation sequencing, which yielded complete Shearwaterpox virus (SWPV) genomes from A. pacificus and A. carneipes, both showing the highest degree of similarity with Canarypox virus (98% and 67%, respectively). The novel SWPV-1 complete genome from A. carneipes is missing 43 genes compared to CNPV and contains 4 predicted genes which are not found in any other poxvirus, whilst, SWPV-2 complete genome was deemed to be missing 18 genes compared to CNPV and a further 15 genes significantly fragmented as to probably cause them to be non-functional.

Conclusion

These are the first avipoxvirus complete genome sequences that infect marine seabirds. In the comparison of SWPV-1 and −2 to existing avipoxvirus sequences, our results indicate that the SWPV complete genome from A. carneipes (SWPV-1) described here is not closely related to any other avipoxvirus genome isolated from avian or other natural host species, and that it likely should be considered a separate species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3680-z) contains supplementary material, which is available to authorized users.

The mitochondrial genome of the ethanol-metabolizing, wine cellar mold Zasmidium cellare is the smallest for a filamentous ascomycete

Fungi in the class Dothideomycetes often live in extreme environments or have unusual physiology. One of these, the wine cellar mold Zasmidium cellare, produces thick curtains of mycelia in cellars with high humidity, and its ability to metabolize volatile organic compounds is thought to improve air quality. Whether these abilities have affected its mitochondrial genome is not known. To fill this gap, the circular-mapping mitochondrial genome of Z. cellare was sequenced and, at only 23 743 bp, is the smallest reported for a filamentous fungus. Genes were encoded on both strands with a single change of direction, different from most other fungi but consistent with the Dothideomycetes. Other than its small size, the only unusual feature of the Z. cellare mitochondrial genome was two copies of a 110-bp sequence that were duplicated, inverted and separated by approximately 1 kb. This inverted-repeat sequence confused the assembly program but appears to have no functional significance. The small size of the Z. cellare mitochondrial genome was due to slightly smaller genes, lack of introns and non-essential genes, reduced intergenic spacers and very few ORFs relative to other fungi rather than a loss of essential genes. Whether this reduction facilitates its unusual biology remains unknown.

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.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Complete mitochondrial genomes of chimpanzee- and gibbon-derived Ascaris isolated from a zoological garden in southwest China

Roundworms (Ascaridida: Nematoda), one of the most common soil-transmitted helminths (STHs), can cause ascariasis in various hosts worldwide, ranging from wild to domestic animals and humans. Despite the veterinary and health importance of the Ascaridida species, little or no attention has been paid to roundworms infecting wild animals including non-human primates due to the current taxon sampling and survey bias in this order. Importantly, there has been considerable controversy over the years as to whether Ascaris species infecting non-human primates are the same as or distinct from Ascaris lumbricoides infecting humans. Herein, we first characterized the complete mitochondrial genomes of two representative Ascaris isolates derived from two non-human primates, namely, chimpanzees (Pan troglodytes) and gibbons (Hylobates hoolock), in a zoological garden of southwest China and compared them with those of A. lumbricoides and the congeneric Ascaris suum as well as other related species in the same order, and then used comparative mitogenomics, genome-wide nucleotide sequence identity analysis, and phylogeny to determine whether the parasites from chimpanzees and gibbons represent a single species and share genetic similarity with A. lumbricoides. Taken together, our results yielded strong statistical support for the hypothesis that the chimpanzee- and gibbon-derived Ascaris represent a single species that is genetically similar to A. lumbricoides, consistent with the results of previous morphological and molecular studies. Our finding should enhance public alertness to roundworms originating from chimpanzees and gibbons and the mtDNA data presented here also serves to enrich the resource of markers that can be used in molecular diagnostic, systematic, population genetic, and evolutionary biological studies of parasitic nematodes from either wild or domestic hosts.

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.

Multiple copies of cytochrome oxidase 1 in species of the fungal genus Fusarium

Using data from published mitochondrial or complete genomes, we developed and tested primers for amplification and sequencing of the barcode region of cytochrome oxidase 1 (COX1) of the fungal genus Fusarium, related genera of the order Hypocreales, and degenerate primers for fungi in the subdivision Pezizomycotina. The primers were successful for amplifying and sequencing COX1 barcodes from 13 genera of Hypocreales (Acremonium, Beauveria, Clonostachys, Emericellopsis, Fusarium, Gliocladium, Hypocrea, Lanatonectria, Lecanicillium, Metarhizium, Monocillium, Neonectria and Stilbella), 22 taxa of Fusarium, and two genera in other orders (Arthrosporium, Monilochaetes). Parologous copies of COX1 occurred in several strains of Fusarium. In some, copies of the same length were detected either by heterozygous bases in otherwise clean sequences or in different replicates of amplification and sequencing events; this may indicate multiple transcribed copies. Other strains included one or two introns. Two intron insertion sites had at least two nonhomologous intron sequences among Fusarium species. Irrespective of whether the multiple copy issue could be resolved by sequencing RNA transcripts, developing a precise COX1-based barcoding system for Fusarium may not be feasible. The overall divergence among homologous COX1 sequences obtained so far is rather low, with many species sharing identical sequences.

A phylogenetic analysis of Greek isolates of Aspergillus species based on morphology and nuclear and mitochondrial gene sequences

Aspergillus species originating from Greece were examined by morphological and molecular criteria to explore the diversity of this genus. The phylogenetic relationships of these species were determined using sequences from the ITS and IGS region of the nuclear rRNA gene complex, two nuclear genes ( β -tubulin (benA) and RNA polymerase II second largest subunit (rpb2)) and two mitochondrial genes (small rRNA subunit (rns) and cytochrome oxidase subunit I (cox1)) and, where available, related sequences from databases. The morphological characters of the anamorphs and teleomorphs, and the single gene phylogenetic trees, differentiated and placed the species examined in the well-supported sections of Aenei, Aspergillus, Bispori, Candidi, Circumdati, Clavati, Cremei, Flavi, Flavipedes, Fumigati, Nidulantes, Nigri, Restricti, Terrei, Usti, and Zonati, with few uncertainties. The combined use of the three commonly employed nuclear genes (benA, rpb2, and ITS), the IGS region, and two less often used mitochondrial gene sequences (rns and cox1) as a single unit resolved several taxonomic ambiguities. A phylogenetic tree was inferred using Neighbour-Joining, Maximum Parsimony, and Bayesian methods. The strains examined formed seven well-supported clades within the genus Aspergillus. Altogether, the concatenated nuclear and mitochondrial sequences offer additional tools for an improved understanding of phylogenetic relationships within this genus.

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 genome of compactin-producing fungus Penicillium solitum and comparative analysis of Trichocomaceae mitochondrial genomes

We determined the complete mitochondrial genome sequence of the compactin-producing fungus Penicillium solitum strain 20-01. The 28 601-base pair circular-mapping DNA molecule encodes a characteristic set of mitochondrial proteins and RNA genes and is intron-free. All 46 protein- and RNA-encoding genes are located on one strand and apparently transcribed in one direction. Comparative analysis of this mtDNA and previously sequenced but unannotated mitochondrial genomes of several medically and industrially important species of the Aspergillus/Penicillium group revealed their extensive similarity in terms of size, gene content and sequence, which is also reflected in the almost perfect conservation of mitochondrial gene order in Penicillium and Aspergillus. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed the monophyletic origin of Eurotiomycetes.

Cordyceps militaris (Hypocreales: Cordycipitaceae): transcriptional analysis and molecular characterization of cox1 and group I intron with putative LAGLIDADG endonuclease

The full-length cytochrome c oxidase subunit I gene (cox1) containing a group I intron was isolated from an important medical fungus Cordyceps militaris (Cordycipitaceae). The open reading frame (ORF) of 1,593 nucleotides encoded a predicted protein COX1 of 530 amino acids. The group I intron encoded a putative homing endonuclease (HE) with two LAGLIDADG motifs. RT-PCR and Northern analysis showed a mature transcript of spliced cox1. Both 5'exon-intron and intron-3'exon junctions were also found by RT-PCR, suggesting the possible presence of unspliced cox1 RNA in total RNA. Sequence comparison by BLASTn showed that the coding region of cox1 (CRcox1) of C. militaris had significant similarities to those of related species (such as Cordyceps bassiana and C. brongniartii), while the intron had no significant homologous sequences of Cordycipitaceae fungi in NCBI database. The phylogenetic tree based on the CRcox1 confirmed the present taxonomic status of related species, but the cox1 introns were phylogenetically distinct. Compared to C. bassiana and C. brongniartii, the cox1 intron of C. militaris exhibited specific splicing site and different intronic ORF. The analysis of the folding RNA structures of the known cox1 introns from Cordyceps species showed different base pairs and conserved regions (P1-P10) in their structures. The present results provide useful information on the studies of cox1 intron splicing and Cordyceps evolution.

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.

Recent dermatophyte divergence revealed by comparative and phylogenetic analysis of mitochondrial genomes

Background

Dermatophytes are fungi that cause superficial infections of the skin, hair, and nails. They are the most common agents of fungal infections worldwide. Dermatophytic fungi constitute three genera, Trichophyton, Epidermophyton, and Microsporum, and the evolutionary relationships between these genera are epidemiologically important. Mitochondria are considered to be of monophyletic origin and mitochondrial sequences offer many advantages for phylogenetic studies. However, only one complete dermatophyte mitochondrial genome (E. floccosum) has previously been determined.

Results

The complete mitochondrial DNA sequences of five dermatophyte species, T. rubrum (26,985 bp), T. mentagrophytes (24,297 bp), T. ajelloi (28,530 bp), M. canis (23,943 bp) and M. nanum (24,105 bp) were determined. These were compared to the E. floccosum sequence. Mitochondrial genomes of all 6 species were found to harbor the same set of genes arranged identical order indicating that these dermatophytes are closely related. Genome size differences were largely due to variable lengths of non-coding intergenic regions and the presence/absence of introns. Phylogenetic analyses based on complete mitochondrial genomes reveals that the divergence of the dermatophyte clade was later than of other groups of pathogenic fungi.

Conclusion

This is the first systematic comparative genomic study on dermatophytes, a highly conserved and recently-diverged lineage of ascomycota fungi. The data reported here provide a basis for further exploration of interrelationships between dermatophytes and will contribute to the study of mitochondrial evolution in higher fungi.

Associated links among mtDNA glycation, oxidative stress and colony sectorization in Metarhizium anisopliae

Mycelial colonies of filamentous fungi often deteriorate when maintained on artificial media, and this can take the form of sterile sectors. We previously established that sectorization by the entomopathogenic fungus Metarhizium anisopliae correlates with intracellular accumulation of reactive oxygen species (ROS). In this study we demonstrate that: (1) H(2)O(2) increases rates of sectorization; (2) a stable strain of M. anisopliae eliminates intracellular ROS more rapidly than an unstable strain; (3) mitochondrial DNA from sectors undergoes a non-enzymatic glycation of deoxyguanosine that is not shown by genomic DNA; (4) the membrane potential of mitochondria in sector cells is decreased in comparison to wild type cells indicating loss of function; (5) DNA glycation changes the properties of DNA and (6) treating wild type mycelia with H(2)O(2) reproduced the glycation pattern shown in sectors. H(2)O(2) also reproduced the morphological changes in mitochondria and lipid droplets that occur in sector cells. Fungal sectorization thus displays aging related developmental impairments resulting from oxidative stress, suggesting a new research direction for studies on fungal colony deterioration. Mitochondrial DNA has a very high AT bias. We speculate that reducing the consequences of glycation could provide an adaptive reason for this.

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.

Dothideomycete plant interactions illuminated by genome sequencing and EST analysis of the wheat pathogen Stagonospora nodorum

Stagonospora nodorum is a major necrotrophic fungal pathogen of wheat (Triticum aestivum) and a member of the Dothideomycetes, a large fungal taxon that includes many important plant pathogens affecting all major crop plant families. Here, we report the acquisition and initial analysis of a draft genome sequence for this fungus. The assembly comprises 37,164,227 bp of nuclear DNA contained in 107 scaffolds. The circular mitochondrial genome comprises 49,761 bp encoding 46 genes, including four that are intron encoded. The nuclear genome assembly contains 26 classes of repetitive DNA, comprising 4.5% of the genome. Some of the repeats show evidence of repeat-induced point mutations consistent with a frequent sexual cycle. ESTs and gene prediction models support a minimum of 10,762 nuclear genes. Extensive orthology was found between the polyketide synthase family in S. nodorum and Cochliobolus heterostrophus, suggesting an ancient origin and conserved functions for these genes. A striking feature of the gene catalog was the large number of genes predicted to encode secreted proteins; the majority has no meaningful similarity to any other known genes. It is likely that genes for host-specific toxins, in addition to ToxA, will be found among this group. ESTs obtained from axenic mycelium grown on oleate (chosen to mimic early infection) and late-stage lesions sporulating on wheat leaves were obtained. Statistical analysis shows that transcripts encoding proteins involved in protein synthesis and in the production of extracellular proteases, cellulases, and xylanases predominate in the infection library. This suggests that the fungus is dependant on the degradation of wheat macromolecular constituents to provide the carbon skeletons and energy for the synthesis of proteins and other components destined for the developing pycnidiospores.

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 sequence (27,184 bp) of the mitochondrial (mt) genome of the phytopathogenic fungus Verticillium dahliae has been determined. It contains 14 protein-coding genes related to oxidative phosphorylation, two rRNA genes and a set of 25 tRNA genes. A single intron, that harbors an intronic ORF coding for a putative ribosomal protein (rps), is located within the large rRNA gene (rnl). Gene order comparisons of V. dahliae mtDNA and complete mt genomes of Pezizomycotina revealed four units of synteny for Sordariomycetes, namely rnl-trn ((11-12))-nad2-nad3, nad4L-nad5-cob-cox1, nad1-nad4-atp8-atp6 and rns-trn ((1-5))-cox3-trn ((1-5))-nad6-trn ((2-5)). These four units, in different combinations, merged to single continuous unit in the orders of Hypocreales and Sordariales. V. dahliae (Phyllachorales) and all members of the genus showed a unique feature which is the translocation of the nad1-nad4-atp8-atp6-rns-cox3-nad6 region in between genes nad3 and atp9 of the Hypocreales mtDNA gene order. Analysis of mt intergenic sequences of Verticillium species permitted the design of a species-specific primer allowing the discrimination of V. longisporum against V. dahliae and V. albo-atrum. By considering the protein-coding gene sequences as one unit, a phylogenetic comparison with representatives of Ascomycota complete mtDNA was performed.