Open Access

Genome sequences of Knoxdaviesia capensis and K. proteae (Fungi: Ascomycota) from Protea trees in South Africa

  • Janneke Aylward1Email author,
  • Emma T. Steenkamp2,
  • Léanne L. Dreyer1,
  • Francois Roets3,
  • Brenda D. Wingfield4 and
  • Michael J. Wingfield2
Standards in Genomic Sciences201611:22

https://doi.org/10.1186/s40793-016-0139-9

Received: 13 November 2015

Accepted: 18 February 2016

Published: 29 February 2016

Abstract

Two closely related ophiostomatoid fungi, Knoxdaviesia capensis and K. proteae, inhabit the fruiting structures of certain Protea species indigenous to southern Africa. Although K. capensis occurs in several Protea hosts, K. proteae is confined to P. repens. In this study, the genomes of K. capensis CBS139037 and K. proteae CBS140089 are determined. The genome of K. capensis consists of 35,537,816 bp assembled into 29 scaffolds and 7940 predicted protein-coding genes of which 6192 (77.98 %) could be functionally classified. K. proteae has a similar genome size of 35,489,142 bp that is comprised of 133 scaffolds. A total of 8173 protein-coding genes were predicted for K. proteae and 6093 (74.55 %) of these have functional annotations. The GC-content of both genomes is 52.8 %.

Keywords

Knoxdaviesia Gondwanamycetaceae MicroascalesOphiostomatoid fungi Protea

Introduction

Two lineages of the polyphyletic assemblage known as ophiostomatoid fungi [1] are associated with the fruiting structures (infructescences) of serotinous Protea L. plants [2]. Protea species are a key component of the fynbos vegetation in the Core Cape Subregion (CCR) of South Africa [3] and the genus is predominantly encountered in South Africa [4, 5]. The Protea-associated ophiostomatoid fungi are, therefore, believed to be endemic to this region, similar to their hosts. This association of ophiostomatoid fungi with a keystone plant genus in a biodiversity hotspot is intriguing [6], as many ophiostomatoid fungi are notorious pathogens of trees [710], yet the Protea ophiostomatoid species are not associated with disease symptoms [11].

Ophiostomatoid fungi are characterized by the flask-shaped morphology of their sexual fruiting structures and their association with arthropods [1, 12]. The Protea-associated members of this assemblage are primarily dispersed by mites that come into contact with fungal spores in the Protea infructescences [13, 14]. These mites have limited dispersal ability, but use beetles and possibly larger vertebrates (such as birds) as vehicles for long-distance dispersal [15, 16].

The three Knoxdaviesia M.J. Wingf., P.S. van Wyk & Marasas species associated with Protea have intriguing host ranges. K. capensis M.J. Wingf. & P.S. van Wyk occurs in at least eight different Protea hosts, whereas K. proteae M.J. Wingf., P.S. van Wyk & Marasas and K. wingfieldii (Roets & Dreyer) Z.W. de Beer & M.J. Wingf. are confined to single host species, respectively P. repens L. and P. caffra Meisn.[1720]. An investigation of the population biology of K. proteae , revealed that this fungus has a high level of intra-specific genetic diversity and that it is extensively dispersed within the CCR of South Africa [16, 21]. However, other than host range and dispersal mechanisms, little is known about the biology and ecology of Knoxdaviesia in general [11]. Here we present the description of the first drafts of the genome sequences of the two CCR species, K. capensis and K. proteae , as well as their respective annotations.

Organism information

Classification and features

The one lineage of Protea-associated ophiostomatoid fungi resides in the Ophiostomataceae (Ophiostomatales, Ascomycota), while the second resides in the Gondwanamycetaceae (Microascales, Ascomycota) [11, 22]. The latter group includes three closely related Protea-associated species in the genus Knoxdaviesia (Fig. 1). This genus was initially described to accommodate the asexual state of the first species in the genus, K. proteae [23]. Under the dual nomenclature system of fungi, the sexual state of this fungus was described in the same paper as Ceratocystiopsis proteae M.J. Wingf., P.S. van Wyk & Marasas [23]. A new genus, Gondwanamyces G.J. Marais & M.J. Wingf., was later described to accommodate the sexual state of this species and that of another species, Ophiostoma capense M.J. Wingf. & P.S. van Wyk [24]. The asexual states of both remained to be treated as species of Knoxdaviesia . Since the abolishment of the dual nomenclature system of fungi, the oldest genus name takes preference, irrespective of morph [25, 26]. The name Knoxdaviesia , therefore, has priority and all species previously treated in Gondwanamyces were transferred to Knoxdaviesia [27].
Fig. 1

Maximum Likelihood tree illustrating the phylogenetic position of K. capensis and K. proteae in the Gondwanamycetaceae (grey block). The Protea-associated species are shaded red and the two isolates for which genome sequences were determined are indicated with a box. The sequences of the Internal Transcribed Spacer (ITS) region (available from GenBank®, accession numbers in brackets following isolate numbers) were aligned in MAFFT 7 [55]. The phylogeny was calculated in MEGA6 [56] using the Tamura-Nei substitution model [57], 1000 bootstrap replicates and Ceratocystis fimbriata (Ceratocystidaceae) as an outgroup

In a study determining the genome sequence of any fungus, it is advisable to use a living isolate connected to the type specimen. However, the ex-type isolate of K. proteae (CMW738 = CBS486.88) is more than 20 years old and does not display the characteristic morphological features of the fungus in culture anymore. No living ex-type isolate exists for K. capensis . We thus collected fresh isolates of both species for this study in order to eliminate possible mutations or degradation that may have occurred though continual artificial propagation in culture media. The new isolates (Figs. 1 & 2) were collected from the same localities and hosts as the holotype specimens: K. capensis (CMW40890 = CBS139037) from the infructescences of P. longifolia Andrews in Hermanus, and K. proteae (CMW40880 = CBS140089) from P. repens infructescences in Stellenbosch, both locations in the Western Cape Province of South Africa. General features of these isolates are outlined in Table 1.
Fig. 2

Sexual sporing structures of the two Knoxdaviesia species sequenced in this study. K. capensis (a) and K. proteae (b) were sampled from Protea longifolia and P. repens flowers, respectively. Scale bars = 1 mm

Table 1

Classification and general features of K. capensis and K. proteae [29]

MIGS ID

Property

K. capensis Term

K. proteae Term

Evidence codea

 

Classification

Domain Fungi

Domain Fungi

TAS [19, 23]

  

Phylum Ascomycota

Phylum Ascomycota

TAS [19, 23]

  

Class Sordariomycetes

Class Sordariomycetes

TAS [19, 23]

  

Order Microascales

Order Microascales

TAS [2]

  

Family Gondwanamycetaceae

Family Gondwanamycetaceae

TAS [22]

  

Genus Knoxdaviesia

Genus Knoxdaviesia

TAS [27]

  

Species K. capensis

Species K. proteae

TAS [27]

  

Strain: CMW40890 = CBS139037

Strain: CMW40880 = CBS140089

 
 

Cell shape

septate, smooth-walled hyphae

septate, smooth-walled hyphae

TAS [19, 23]

 

Motility

Non-motile

Non-motile

NAS

 

Sporulation

Unsheathed allantoid ascospores

Falcate ascospores

TAS [19, 23]

 

Temperature range

15–30 °C

15–30 °C

TAS [19, 23]

 

Optimum temperature

25 °C

25 °C

TAS [19, 23]

 

pH range; Optimum

Unknown

Unknown

 
 

Carbon source

Unknown

Unknown

 

MIGS-6

Habitat

Seed cones (infructescences) of Protea spp.

Seed cones (infructescences) of Protea repens L.

TAS [19, 23]

MIGS-6.3

Salinity

Unknown

Unknown

 

MIGS-22

Oxygen requirement

Aerobic; requirement/tolerance unknown

Aerobic; requirement/tolerance unknown

 

MIGS-15

Biotic relationship

Plant-associated

Plant-associated

TAS [24]

MIGS-14

Pathogenicity

None known

None known

 

MIGS-4

Geographic location

Hermanus, South Africa

Stellenbosch, South Africa

 

MIGS-5

Sample collection

February 2014

January 2014

 

MIGS-4.1

Latitude

-34.4093

-33.9430

 

MIGS-4.2

Longitude

19.2150

18.8802

 

MIGS-4.4

Altitude

20 m

140 m

 

aEvidence codes - IDA inferred from direct assay, TAS traceable author statement (i.e., a direct report exists in the literature), NAS non-traceable author statement (i.e., not directly observed for the living, isolated sample, but based on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from http://www.geneontology.org/GO.evidence.shtml of the Gene Ontology project [58]

Genome sequencing information

Genome project history

Considering the lack of ecological information on the genus Knoxdaviesia and the close relationship these Microascalean fungi have to important plant pathogens, two Protea-associated Knoxdaviesia species, believed to be native to the CCR in South Africa, were selected for genome sequencing. Both species were sequenced at Fasteris in Switzerland. The genome projects are listed in the Genomes OnLine Database [28] and the whole genome shotgun (WGS) project has been deposited at DDBJ/EMBL/GenBank (Table 2). Table 2 presents the project information and its association with the minimum information about a genome sequence version 2.0 compliance [29]. The full MIGS records for K. capensis and K. proteae are available in Additional file 1: Table S1 and Additional file 2: Table S2, respectively.
Table 2

Project information

MIGS ID

Property

K. capensis Term

K. proteae Term

MIGS 31

Finishing quality

High quality draft

High quality draft

MIGS-28

Libraries used

2x paired-end (PE) (350 and 550 bp) and 1x mate-pair (MP) (3 kbp)

2x paired-end (PE) (350 and 550 bp) and 1x mate-pair (MP) (3 kbp)

MIGS 29

Sequencing platforms

Illumina Hiseq 2500

Illumina Hiseq 2500

MIGS 31.2

Fold coverage

PE library 1: 91.6 x

PE library 1: 142 x

PE library 2: 80 x

PE library 2: 79.3 x

MP library: 17 x

MP library: 50.2 x

MIGS 30

Assemblers

ABySS 1.5.2; SSPACE 3.0

ABySS 1.5.2; SSPACE 3.0

MIGS 32

Gene calling method

MAKER 2.31.8

MAKER 2.31.8

 

Genbank ID

LNGK00000000

LNGL00000000

 

GenBank Date of Release

11th January 2016

11th January 2016

 

GOLD ID

Gp0093999

Gp0110284

 

BIOPROJECT

PRJNA246171

PRJNA275563

MIGS 13

Source Material Identifier

CMW40890/CBS139037

CMW40880/CBS 140089

 

Project relevance

Biodiversity, evolution

Biodiversity, evolution

Growth conditions and genomic DNA preparation

Both K. capensis and K. proteae were cultured on Malt Extract Agar (MEA; Merck, Wadeville, South Africa) overlaid with sterile cellophane sheets (Product no. Z377597, Sigma-Aldrich, Steinham, Germany). After 10 days of growth at 25 °C, mycelia was scraped from the cellophane and DNA was extracted according to Aylward et al. [30]. Approximately 5 μg DNA from each species was used to prepare the three Illumina libraries (Table 2).

RNA was extracted from the K. proteae genome isolate to use as evidence for gene prediction. After growth on MEA at 25 °C for approximately 10 days, total RNA was isolated from the mycelia with the PureLink™ RNA Mini Kit (Ambion, Austin, TX, USA). Quality control was performed on the Agilent 2100 Bioanalyzer (Agilent Technologies, USA) using the RNA 6000 Nano Assay kit (Agilent Technologies, USA). The mRNA component of the total RNA was subsequently extracted with the Dynabeads® mRNA purification kit (Ambion, Austin, TX, USA).

Genome sequencing and assembly

The genomes of K. capensis and K. proteae were sequenced with the Illumina HiSeq 2500 platform at Fasteris, Switzerland, using two paired-end and one Nextera mate-pair library (Table 2). More than 60 million paired-end and 8 million mate-pair reads were obtained for each species. These reads were trimmed in CLC Genomics Workbench 6.5 (CLC bio, Aarhus, Denmark) so that the Phred Q (quality) score of each base was at least Q20. VelvetOptimiser (Gladman & Seeman, unpublished), a Perl script used as part of the Velvet assembler [31, 32], was initially used to optimize the assembly parameters. Assembly of contigs was performed in ABySS 1.5.2 [33] using the optimal parameters suggested by VelvetOptimiser as a starting point. Several assemblies were computed using kmer-values slightly higher and lower than the kmer-value suggested by VelvetOptimiser. The assembly with the lowest number of contigs was used to build scaffolds in SSPACE 3.0 [34], discarding scaffolds smaller than 1000 bp. Automatic gap closure was performed in GapFiller 1.10 [35]. The average genome coverage of each library was estimated using the Lander-Waterman equation (total sequenced nucleotides/genome size) (Table 2), which yielded a combined average coverage for the three libraries of 188.5x ( K. capensis ) and 271.5x ( K. proteae ).

The K. capensis genome consists of 29 scaffolds ranging between 1226 and 5,637,848 bp, whereas the 133 scaffolds of K. proteae are sized between 1022 and 2,610,973 bp. A search for the 1438 fungal universal single-copy ortholog genes with BUSCO 1.1b1 [36] identified 1355 complete and 67 partial genes in K. capensis and 1366 complete and 57 partial genes in K. proteae . The two genomes are therefore estimated to be >98 % complete.

The extracted mRNA of K. proteae was sequenced using an Ion PI™ Chip on the Ion Proton™ System (Life Technologies, Carlsbad, CA) at the Central Analytical Facility (CAF), Stellenbosch University, South Africa. The >49 million raw RNA-Seq reads were mapped to the K. capensis genome in CLC Genomics Workbench and assembled with Trinity 2.0.6 [37] using the genome-guided option.

Genome annotation

Genome annotation was performed with the MAKER 2.31.8 pipeline [38, 39], using custom repeat libraries for each species constructed with RepeatScout 1.0.5 [40] and two de novo gene predictors, SNAP 2006-07-28 [41] and AUGUSTUS 3.0.3 [42]. The assembled K. proteae RNA-Seq and predicted protein and/or transcript sequences from 22 sequenced Sordariomycete species (Additional file 3: Table S3), including two Microascalean fungi, were provided as additional evidence. AUGUSTUS was trained with the assembled K. proteae RNA-Seq data and subsequently MAKER was used to annotate the largest scaffold of the K. capensis and the largest scaffold of the K. proteae assembly, independently. After manually curating all the gene predictions on these scaffolds with Apollo 1.11.8 [43], SNAP was trained with the curated gene predictions of each scaffold and the scaffolds were re-annotated. SNAP was re-trained for each species individually and subsequently both genomes were annotated. EuKaryotic Orthologous Group (KOG) classifications were assigned to the predicted proteins through the WebMGA [44] portal that performs reverse-position-specific BLAST [45] searches on the KOG database [46]. Additional functional annotations were predicted with InterProScan 5.13-52.0 [47, 48], SignalP 4.1 [49] and TMHMM 2.0 [50].

Genome properties

K. capensis and K. proteae have similar genome sizes at 35.54 and 35.49 Mbp, respectively. It was possible to assemble the K. capensis genome into 29 scaffolds larger than 1000 bp, whereas the number of scaffolds above this threshold achieved for K. proteae was 133. Both genomes had a GC content of 52.8 %.

A total of 7940 protein-coding genes were predicted for K. capensis and 8174 for K. proteae . Additionally 137 and 116 tRNA and 30 and 27 rRNA genes were predicted for each species, respectively. More than 74 % of the protein-coding genes of each species could be assigned to a putative function via the KOG and Pfam databases. The content of the two genomes are summarized in Tables 3 and 4.
Table 3

Genome statistics

Species

K. capensis

K. proteae

Attribute

Value

% of Totala

Value

% of Totala

Genome size (bp)

35,537,816

100.00

35,489,142

100.00

DNA coding (bp)

12,640,368

35.57

12,542,580

35.34

DNA G + C (bp)

18,774,628

52.83

18,745,365

52.82

DNA scaffolds

29

 

133

 

Total genes

8107

100.00

8316

100.00

Protein coding genes

7940

97.94

8173

98.28

RNA genesb

167

2.06

143

1.72

Pseudo genes

unknown

 

unknown

 

Genes in internal clusters

unknown

 

unknown

 

Genes with function prediction

6192

77.98

6093

74.55

Genes assigned to KOGs

6059

76.31

6015

73.60

Genes with Pfam domains

5455

68.70

5335

65.28

Genes with signal peptides

354

4.46

335

4.10

Genes with transmembrane helices

1510

19.02

1527

18.68

CRISPR repeats

N/A

 

N/A

 

aThe total is based on either the size of the genome in base pairs or the total number of protein-coding genes in the annotated genome

bBased on tRNA and rRNA genes only

Table 4

Number of genes associated with the 25 general KOG functional categories

Species

K. capensis

K. proteae

 

Code

Value

% of totala

Value

% of totala

Description

J

359

4.52

371

4.54

Translation, ribosomal structure and biogenesis

A

280

3.53

273

3.34

RNA processing and modification

K

475

5.98

484

5.92

Transcription

L

196

2.47

198

2.42

Replication, recombination and repair

B

109

1.37

99

1.21

Chromatin structure and dynamics

D

209

2.63

227

2.78

Cell cycle control, cell division, chromosome partitioning

Y

34

0.43

32

0.39

Nuclear structure

V

32

0.40

32

0.39

Defence mechanisms

T

505

6.36

586

5.95

Signal transduction mechanisms

M

69

0.87

76

0.93

Cell wall/membrane/envelope biogenesis

N

6

0.08

6

0.07

Cell motility

Z

279

3.51

289

3.54

Cytoskeleton

W

10

0.13

12

0.15

Extracellular structures

U

539

6.79

543

6.64

Intracellular trafficking, secretion, and vesicular transport

O

502

6.32

495

6.06

Post-translational modification, protein turnover, chaperones

C

265

3.34

256

3.13

Energy production and conversion

G

202

2.54

202

2.47

Carbohydrate transport and metabolism

E

227

2.86

228

2.79

Amino acid transport and metabolism

F

76

0.96

74

0.91

Nucleotide transport and metabolism

H

87

1.10

85

1.04

Coenzyme transport and metabolism

I

234

2.95

234

2.86

Lipid transport and metabolism

P

144

1.81

151

1.85

Inorganic ion transport and metabolism

Q

139

1.75

137

1.68

Secondary metabolites biosynthesis, transport and catabolism

R

735

9.26

694

8.49

General function prediction only

S

344

4.33

330

4.04

Function unknown

X

2

0.03

1

0.01

Multiple functions

-

1881

23.69

2159

26.41

Not in KOGs

aThe total is based on the total number of protein coding genes in the genome

Conclusions

At least six Microascalean fungi currently have publically accessible genomes [5154]. K. capensis and K. proteae , however, represent the first sequenced genomes from the Microascalean family Gondwanamycetaceae . The genomes of these two species will not only enable further understanding of the unique ecology of Protea-inhabiting fungi, but will also be valuable in taxonomic and evolutionary studies.

Abbreviations

CCR: 

core cape subregion

MEA: 

malt extract agar

KOG: 

EuKaryotic Orthologous Groups of proteins

Declarations

Acknowledgements

We are grateful to Quentin Santana and Dr. Lieschen Bahlmann for their guidance in the genome assembly and annotation procedures and to Dr. Wilhelm de Beer for the taxonomic information he contributed to this manuscript. This research was funded by the National Research Foundation (NRF) and the Department of Science and Technology/NRF Centre of Excellence in Tree Health Biotechnology. We also thank the Cape Nature Conservation Board for supplying the necessary collection permits.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Botany and Zoology, Stellenbosch University
(2)
Department of Microbiology and Plant Pathology, University of Pretoria
(3)
Department of Conservation Ecology and Entomology, Stellenbosch University
(4)
Department of Genetics, University of Pretoria

References

  1. Spatafora JW, Blackwell M. The polyphyletic origins of ophiostomatoid fungi. Mycol Res. 1994;98:1–9.View ArticleGoogle Scholar
  2. Wingfield BD, Viljoen CD, Wingfield MJ. Phylogenetic relationships of ophiostomatoid fungi associated with Protea infructescences in South Africa. Mycol Res. 1999;103:1616–20.View ArticleGoogle Scholar
  3. Manning J, Goldblatt P. Plants of the Greater Cape Floristic Region. 1: The Core Cape Flora. Strelitzia 29, vol 29. Pretoria: South African National Biodiversity Institute; 2012.Google Scholar
  4. Gibbs RG. Analysis of the size and composition of the southern African flora. Bothalia. 1984;15:613–29.Google Scholar
  5. Rebelo T. Proteas: a field guide to the Proteas of Southern Africa. Vlaeberg: Fernwood Press; 1995.Google Scholar
  6. Mittermeier RA, Myers N, Thomsen JB, Da Fonseca GAB, Olivieri S. Biodiversity hotspots and major tropical wilderness areas: approaches to setting conservation priorities. Conserv Biol. 1998;12:516–20.View ArticleGoogle Scholar
  7. Paine TD, Raffa KF, Harrington TC. Interactions among scolytid bark beetles, their associated fungi, and live host conifers. Annu Rev Entomol. 1997;42:179–206.View ArticlePubMedGoogle Scholar
  8. Brasier CM. Ophiostoma novo-ulmi sp. nov., causative agent of current Dutch elm disease pandemics. Mycopathologia. 1991;115:151–61.View ArticleGoogle Scholar
  9. Harrington TC, Wingfield MJ. The Ceratocystis species on conifers. Can J Botany. 1998;76:1446–57.Google Scholar
  10. Roux J, Wingfield MJ. Ceratocystis species: emerging pathogens of non-native plantation Eucalyptus and Acacia species. South Forests. 2009;71:115–20.View ArticleGoogle Scholar
  11. Roets F, Wingfield MJ, Crous PW, Dreyer LL. Taxonomy and ecology of ophiostomatoid fungi associated with Protea infructescences. In: Seifert KA, de Beer ZW, Wingfield MJ, editors. Ophiostomatoid fungi: expanding frontiers. Utrecht: CBS Biodiversity Series; 2013. p. 177–87.Google Scholar
  12. Malloch D, Blackwell M. Dispersal biology of the ophiostomatoid fungi. In: Wingfield MJ, Seifert KA, Webber JF, editors. Ceratocystis and Ophiostoma: taxonomy, ecology and pathology. St. Paul: APS Press; 1993. p. 195–206.Google Scholar
  13. Roets F, Wingfield MJ, Crous PW, Dreyer LL. Discovery of Fungus-Mite Mutualism in a Unique Niche. Environ Entomol. 2007;36:1226–37.View ArticlePubMedGoogle Scholar
  14. Roets F, Wingfield MJ, Wingfield BD, Dreyer LL. Mites are the most common vectors of the fungus Gondwanamyces proteae in Protea infructescences. Fungal Biol. 2011;115:343–50.View ArticlePubMedGoogle Scholar
  15. Roets F, Crous PW, Wingfield MJ, Dreyer LL. Mite-mediated hyperphoretic dispersal of Ophiostoma spp. from the Infructescences of South African Protea spp. Environ Entomol. 2009;28:143–52.View ArticleGoogle Scholar
  16. Aylward J, Dreyer LL, Steenkamp ET, Wingfield MJ, Roets F. Long-distance dispersal and recolonization of a fire-destroyed niche by a mite-associated fungus. Fungal Biol. 2015;119:245–56.View ArticlePubMedGoogle Scholar
  17. Roets F, Theron N, Wingfield MJ, Dreyer LL. Biotic and abiotic constraints that facilitate host exclusivity of Gondwanamyces and Ophiostoma on Protea. Fungal Biol. 2011;116:49–61.View ArticlePubMedGoogle Scholar
  18. Roets F, Dreyer LL, Crous PW. Seasonal trends in colonisation of Protea infructescences by Gondwanamyces and Ophiostoma spp. S Afr J Bot. 2005;71:307–11.View ArticleGoogle Scholar
  19. Wingfield MJ, Van Wyk PS. A new species of Ophiostoma from Protea infructescences in South Africa. Mycol Res. 1993;97:709–16.View ArticleGoogle Scholar
  20. Crous PW, Summerell BA, Shivas RG, Burgess TI, Decock CA, Dreyer LL, et al. Fungal Planet description sheets: 107-127. Persoonia. 2012;28:138–82.PubMed CentralView ArticlePubMedGoogle Scholar
  21. Aylward J, Dreyer LL, Steenkamp ET, Wingfield MJ, Roets F. Panmixia defines the genetic diversity of a unique arthropod-dispersed fungus specific to Protea flowers. Ecol Evol. 2014;4:3444–55.PubMed CentralView ArticlePubMedGoogle Scholar
  22. Réblová M, Gams W, Seifert KA. Monilochaetes and allied genera of the Glomerellales, and a reconsideration of families in the Microascales. Stud Mycol. 2011;68:163–91.PubMed CentralView ArticlePubMedGoogle Scholar
  23. Wingfield MJ, Wyk PSV, Marasas WFO. Ceratocystiopsis proteae sp. nov. with a new anamorph genus. Mycologia. 1988;80:23–30.Google Scholar
  24. Marais GJ, Wingfield MJ, Viljoen CD, Wingfield BD. A new ophiostomatoid genus from Protea infructescences. Mycologia. 1998;90:136–41.View ArticleGoogle Scholar
  25. Hawksworth DL. A new dawn for the naming of fungi: impacts of decisions made in Melbourne in July 2011 on the future publication and regulation of fungal names. IMA Fungus. 2011;2:155.PubMed CentralView ArticlePubMedGoogle Scholar
  26. McNeill J, Barrie F, Buck W, Demoulin V, Greuter W, Hawksworth D, et al. International Code of Nomenclature for algae, fungi, and plants (Melbourne Code) adopted by the Eighteenth International Botanical Congress Melbourne, Australia, July 2011. Regnum Vegetabile 154, Koeltz Scientific Books. 2012.Google Scholar
  27. De Beer ZW, Seifert KA, Wingfield MJ. A nomenclature for ophiostomatoid genera and species in the Ophiostomatales and Microascales. In: Seifert KA, de Beer ZW, Wingfield MJ, editors. Ophiostomatoid fungi: expanding frontiers, CBS Biodiversity Series. 2013. p. 245–322.Google Scholar
  28. Pagani I, Liolios K, Jansson J, Chen I-MA, Smirnova T, Nosrat B, et al. The Genomes OnLine Database (GOLD) v. 4: status of genomic and metagenomic projects and their associated metadata. Nucleic Acids Res. 2012;40:D571–9.PubMed CentralView ArticlePubMedGoogle Scholar
  29. Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, et al. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol. 2008;26:541–7.PubMed CentralView ArticlePubMedGoogle Scholar
  30. Aylward J, Dreyer LL, Steenkamp ET, Wingfield MJ, Roets F. Development of polymorphic microsatellite markers for the genetic characterisation of Knoxdaviesia proteae (Ascomycota: Microascales) using ISSR-PCR and pyrosequencing. Mycol Prog. 2014;13:439–44.View ArticleGoogle Scholar
  31. Zerbino DR. Using the Velvet de novo Assembler for Short‐Read Sequencing Technologies. Curr Protoc Bioinformatics. 2010; Unit 11.5:1-13.Google Scholar
  32. Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.PubMed CentralView ArticlePubMedGoogle Scholar
  33. Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol İ. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–23.PubMed CentralView ArticlePubMedGoogle Scholar
  34. Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27:578–9.View ArticlePubMedGoogle Scholar
  35. Boetzer M, Pirovano W. Toward almost closed genomes with GapFiller. Genome Biol. 2012;13:R56.PubMed CentralView ArticlePubMedGoogle Scholar
  36. Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Supplementary Online Materials: http://busco.ezlab.org/files/BUSCO-SOM.pdf. 2015.Google Scholar
  37. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.View ArticlePubMedGoogle Scholar
  38. Cantarel BL, Korf I, Robb SM, Parra G, Ross E, Moore B, et al. MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 2008;18:188–96.PubMed CentralView ArticlePubMedGoogle Scholar
  39. Holt C, Yandell M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinformatics. 2011;12:491.PubMed CentralView ArticlePubMedGoogle Scholar
  40. Price AL, Jones NC, Pevzner PA. De novo identification of repeat families in large genomes. Bioinformatics. 2005;21:i351–8.View ArticlePubMedGoogle Scholar
  41. Korf I. Gene finding in novel genomes. BMC Bioinformatics. 2004;5:59.PubMed CentralView ArticlePubMedGoogle Scholar
  42. Stanke M, Steinkamp R, Waack S, Morgenstern B. AUGUSTUS: a web server for gene finding in eukaryotes. Nucleic Acids Res. 2004;32:W309–12.PubMed CentralView ArticlePubMedGoogle Scholar
  43. Lewis SE, Searle S, Harris N, Gibson M, Lyer V, Richter J, et al. Apollo: a sequence annotation editor. Genome Biol. 2002;3:1–14.View ArticleGoogle Scholar
  44. Wu S, Zhu Z, Fu L, Niu B, Li W. WebMGA: a customizable web server for fast metagenomic sequence analysis. BMC Genomics. 2011;12:444.PubMed CentralView ArticlePubMedGoogle Scholar
  45. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
  46. Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41.PubMed CentralView ArticlePubMedGoogle Scholar
  47. Goujon M, McWilliam H, Li W, Valentin F, Squizzato S, Paern J, et al. A new bioinformatics analysis tools framework at EMBL–EBI. Nucleic Acids Res. 2010;38:W695–9.PubMed CentralView ArticlePubMedGoogle Scholar
  48. Zdobnov EM, Apweiler R. InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001;17:847–8.View ArticlePubMedGoogle Scholar
  49. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Meth. 2011;8:785–6.View ArticleGoogle Scholar
  50. Krogh A, Larsson B, Von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.View ArticlePubMedGoogle Scholar
  51. Van der Nest MA, Bihon W, De Vos L, Naidoo K, Roodt D, Rubagotti E, et al. Draft genome sequences of Diplodia sapinea, Ceratocystis manginecans, and Ceratocystis moniliformis. IMA Fungus. 2014;5:135–40.PubMed CentralView ArticlePubMedGoogle Scholar
  52. Van der Nest MA, Beirn LA, Crouch JA, Demers JE, de Beer ZW, De Vos L, et al. IMA Genome-F 3: Draft genomes of Amanita jacksonii, Ceratocystis albifundus, Fusarium circinatum, Huntiella omanensis, Leptographium procerum, Rutstroemia sydowiana, and Sclerotinia echinophila. IMA Fungus. 2014;5:473–86.PubMed CentralPubMedGoogle Scholar
  53. Vandeputte P, Ghamrawi S, Rechenmann M, Iltis A, Giraud S, Fleury M, et al. Draft genome sequence of the pathogenic fungus Scedosporium apiospermum. Genome Announc. 2014;2:e00988-14.PubMed CentralView ArticlePubMedGoogle Scholar
  54. Wilken PM, Steenkamp ET, Wingfield MJ, De Beer ZW, Wingfield BD. IMA Genome-F 1: Ceratocystis fimbriata: draft nuclear genome sequence for the plant pathogen, Ceratocystis fimbriata. IMA Fungus. 2013;4:357.PubMed CentralView ArticlePubMedGoogle Scholar
  55. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.PubMed CentralView ArticlePubMedGoogle Scholar
  56. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.PubMed CentralView ArticlePubMedGoogle Scholar
  57. Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10:512–26.PubMedGoogle Scholar
  58. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.PubMed CentralView ArticlePubMedGoogle Scholar

Copyright

© Aylward et al. 2016