- Research
- Open access
- Published:
New insights into the DNA extraction and PCR amplification of minute ascomycetes in the genus Laboulbenia (Pezizomycotina, Laboulbeniales)
IMA Fungus volume 15, Article number: 14 (2024)
Abstract
Molecular studies of fungi within the order Laboulbeniales (Ascomycota, Pezizomycotina) have been hampered for years because of their minute size, inability to grow in axenic culture, and lack of reliable and cost-efficient DNA extraction protocols. In particular, the genus Laboulbenia is notorious for low success with DNA extraction and polymerase chain reaction (PCR) amplification. This is attributed to the presence of melanin, a molecule known to inhibit PCR, in the cells. We evaluated the efficacy of a standard single cell-based DNA extraction protocol by halving the recommended amount of reagents to reduce the cost per extraction and adding bovine serum albumin (BSA) during the multiple displacement amplification step to reverse the effect of melanin. A total of 196 extractions were made, 111 of which were successful. We found that halving the reagents used in the single cell-based extraction kit did not significantly affect the probability of successful DNA extraction. Using the halved protocol reduces cost and resource consumption. Moreover, there was no significant difference in the probability of successfully extracting DNA based on whether BSA was added or not, suggesting that the amount of melanin present in cells of the thallus has no major inhibitory effect on PCR. We generated 277 sequences from five loci, but amplification of the internal transcribed spacer region, the mitochondrial small subunit rDNA, and protein-coding genes remains challenging. The probability of successfully extracting DNA from Laboulbeniales was also impacted by specimen storage methods, with material preserved in > 95% ethanol yielding higher success rates compared to material stored in 70% ethanol and dried material. We emphasize the importance of proper preservation of material and propose the design of Laboulbeniales-specific primers to overcome the problems of primer mismatches and contaminants. Our new insights apply not only to the genus Laboulbenia; Laboulbeniales generally are understudied, and the vast majority of species remain unsequenced. New and approachable molecular developments will benefit the study of Laboulbeniales, helping to elucidate the true diversity and evolutionary relationships of these peculiar microfungi.
Introduction
Traditionally, fungal species identification and delimitation relied on morphological characteristics, but phenotypic plasticity within species and cryptic species posed challenges (Bridge et al. 2005; Cao et al. 2021; Maharachchikumbura et al. 2021). Some groups used alternative traits like enzymatic activity or mating compatibility (Perkins and Raju 1986; Pincus et al. 2007), but these methods are not applicable to unculturable species, thus other unambiguous characters should be used. Early molecular methods based on GC-content and DNA hybridization had limited use, except in yeast studies (Bridge et al. 2005). The introduction of PCR enabled the discrimination of closely related taxa based on molecular data, revolutionizing molecular systematics in mycology (White et al. 1990; Bridge et al. 1998, 2005; Cao et al. 2021; Maharachchikumbura et al. 2021). The internal transcribed spacer of the nuclear ribosomal DNA (ITS) was proposed as the universal fungal DNA barcode marker (Schoch et al. 2012). While the ITS is nowadays often used to delimit and identify fungi, the interspecific variation in this region is too low in some groups, necessitating the use of additional markers (Stielow et al. 2015).
A particularly understudied group of fungi is the order Laboulbeniales (Ascomycota, Pezizomycotina). These microfungi have an obligate association with arthropod hosts. Instead of producing hyphae, Laboulbeniales develop a 3-D structure called a thallus from a two-celled ascospore, which attaches externally to an arthropod’s integument (Blackwell et al. 2020). Molecular studies of Laboulbeniales have been challenging due to the minute size of these fungi (200–300 μm on average), their melanized cells, and their inability to grow in axenic culture (Haelewaters et al. 2015; Sundberg et al. 2018a). Initial attempts using a dry ice protocol by Weir and Blackwell (2001) only had a 25% success rate. Updated versions of the same protocol were used in a few studies (Goldmann and Weir 2012, 2018; Goldmann et al. 2013). However, these protocols are time-consuming, require more than one thallus per extraction, and result in only limited success of extraction and PCR amplification (Sundberg et al. 2018a). Haelewaters et al. (2015) evaluated four DNA extraction protocols and different pre-treatments with mixed success. They had difficulties extracting DNA of Laboulbenia species, with success rates between 0 and 20%. DNA extraction and amplification of Laboulbenia species has been suggested to be particularly difficult, because many species in the genus have melanized cells, which is known to interfere with DNA polymerase during PCR (Eckhart et al. 2000; Gibson 2012; Haelewaters et al. 2015; Sundberg et al. 2018a). Several PCR inhibitor removal kits are available to remove molecules such as melanin and phenols (Hu et al. 2015; Vicente et al. 2019). Similar results can be reached by adding bovine serum albumin (BSA) to the PCR mixture (Giambernardi et al. 1998), but this has not yet been tested for Laboulbenia species nor Laboulbeniales in general.
Sundberg et al. (2018a) developed a new DNA extraction protocol, in which a manual press system was used to disrupt the tough cell walls of Laboulbeniales. Using one thallus per DNA extraction, they successfully obtained 156 sequences: 20 nrSSU, 56 ITS, 59 nrLSU, and 21 mtSSU. Their protocol does not involve any other treatment of the thalli, which is a major advantage compared to the other described methods. However, drawbacks for their protocol include the need for custom-made components and the fact that DNA extractions are fully depleted during PCR amplification and thus cannot be stored.
Haelewaters et al. (2018a) published another method to extract DNA of Laboulbeniales. Using the REPLI-g Single Cell Kit (Qiagen), a sufficient amount of DNA can be obtained from a single thallus. This protocol is different from the previous methods, as it involves a whole-genome amplification (WGA) step. Due to the WGA, there is a higher risk to amplify contaminants. Yet, it has been successfully used in many other studies by Haelewaters and colleagues (Haelewaters et al. 2018b, 2019a, 2019b, 2022; Walker et al. 2018; Haelewaters and Pfister 2019; Haelewaters and De Kesel 2020; Liu et al. 2020; Van Caenegem et al. 2023a, 2023b). Haelewaters et al. (2019b) used a modified protocol, in which they halved the use of every component, to save products and reduce costs per extraction. There are doubts about the effectiveness of this modified protocol (D. Haelewaters and P. Mironova, pers. comm.), but no formal tests have been performed to evaluate the significance of these doubts.
Currently, the nuclear small and large subunit ribosomal RNA (nrSSU and nrLSU), the internal transcribed spacer region (ITS), the minichromosome maintenance complex component 7 protein-coding gene (MCM7), the translation elongation factor 1α protein-coding gene (TEF1), and the mitochondrial small subunit rRNA (mtSSU) have been sequenced for several species of Laboulbeniales (Goldmann and Weir 2012; Goldmann et al. 2013; Haelewaters et al. 2015, 2018a, 2019b, 2022; Sundberg et al. 2018a, 2018b; Liu et al. 2020; Van Caenegem et al. 2023a, 2023b). General fungal primers designed by White et al. (1990) have mainly been used to amplify regions of ribosomal RNA (nrSSU, ITS, nrLSU). Haelewaters et al. (2015) developed a Laboulbeniomycetes-specific nrSSU primer pair. Additionally, recent research reported low amplification of the ITS region using general fungal primers, which resulted in the design of a Hesperomyces-specific ITS primer pair (ITShespL and ITShespR) and the Laboulbeniomycetes-specific LabITS1 forward primer (Haelewaters et al. 2018a, 2019b).
Only 10 nrSSU sequences of 9 species, 12 ITS sequences of 4 species, and 34 nrLSU sequences of 12 species of Laboulbenia are available in GenBank. Given that almost 700 species are described (Haelewaters et al. 2024), there is a huge discrepancy between described and sequenced species of Laboulbenia. The main goal of this study was to generate DNA sequences of Laboulbenia species for future molecular studies, given their paucity in public sequence databases. We had the opportunity to test questions regarding DNA extraction protocols, primer pairs, and PCR protocols. During the quest to consistently extract and sequence DNA of Laboulbenia, we (1) researched how the preservation methods correlate with DNA extraction success, (2) explored the boundaries of the REPLI-g Single Cell Kit by halving the amount of reagents per extraction, (3) tested whether the addition of BSA results in more successful DNA extractions, and (4) identified the usefulness of different primer pairs and PCR protocols for multiple loci to successfully generate high-quality DNA sequences of Laboulbenia species.
Methods
Collection and processing of beetles
Beetles (Coleoptera) were collected using different entomological methods (light traps, pitfall traps, mouth-operated aspirator, and hand collection) by the authors and by entomologists who sent specimens for study of their Laboulbeniales. Specimens were collected in 70% to 99% ethanol. All specimens were transferred to 99% ethanol upon arrival in the lab at Ghent University. Beetles were screened for the presence of Laboulbeniales using an RZB-PL 65.500 stereoscope (Novex, Arnhem, The Netherlands). Infected specimens were identified by the authors (using Muilwijk et al. 2015) or their collectors. Other host specimens were sent to Dr. Menno Schilthuizen (Taxon Expeditions, Leiden, The Netherlands) for identification and subsequent vouchering. Specimens are stored in the collection of Taxon Expeditions (TXEX, Leiden, The Netherlands) or the entomology collections of the Royal Belgian Institute of Natural Sciences (KBIN, Brussels, Belgium).
Morphological study of Laboulbeniales
Laboulbeniales microfungi were mounted in permanent slides as described by Liu et al. (2020). A 1:1 mixture of Hoyer’s medium and glycerin was used, as pure Hoyer’s medium dries too quickly. A small droplet of water was placed on a microscope slide, on which a 22 × 22 mm coverslip was put. The purpose of this was to ensure that the 22 × 22 mm coverslip was somewhat stuck to the microscope slide during further manipulation and thus could not move unexpectedly or fall. A droplet of the Hoyer’s/glycerin mixture was placed off-center on the coverslip. A hypodermic needle was used to remove Laboulbeniales thalli from the host and place them in the droplet. Thalli were arranged in one row or multiple rows in the middle of the coverslip. A smaller 18 × 18 mm coverslip with a drop of Amann’s medium was flipped upside down (drop of Amann’s medium facing down) and positioned over the thalli by gently lowering it with a dissecting pin. Next, the corners of the 18 × 18 mm coverslip were sealed to the larger coverslip with nail polish. Solakryl BMX (Ento Sphinx, Pardubice, Czech Republic) was applied to the microscope slide, and the coverslip assembly with the smaller coverslip facing downwards was slowly lowered and gently placed sideways on the microscope slide. Our permanent slides are each composed of a 22 × 22 mm coverslip on top of an 18 × 18 mm coverslip, with the thalli in between those two coverslips, and the microscope slide at the bottom.
Mounted thalli were viewed at 100–400 × magnification under an Olympus BH-2 microscope (Olympus, Center Valley, PA). Thalli were identified based on relevant literature (Thaxter 1896; Majewski 1994; De Kesel et al. 2020; Haelewaters and De Kesel 2020; Santamaria and Pedersen 2021) and supplementary papers (Santamaria et al. 1991). Permanent slides of Laboulbeniales are deposited in the Herbarium Universitatis Gandavensis (GENT).
DNA extraction, PCR amplification, and sequencing
DNA extractions were done using the REPLI-g Single Cell Kit (Qiagen, Stanford, CA). The initial steps are described in Fig. 1. Molecular work was done at the Centre for Molecular Phylogeny and Evolution (CeMoFE) at the Ledeganck Campus of Ghent University. All steps were performed wearing disposable latex gloves. To avoid contamination, hypodermic needles for micromanipulation of thalli were thoroughly cleaned with 70% ethanol and bleach before and after every prepared extraction. Thalli of Laboulbeniales were removed from their host under a dissecting microscope, using a needle inserted onto a glass syringe for holdfast. The tip of the needle was submerged in glycerin to prevent thalli from flying away during transfer. The removed thalli were placed in a droplet of glycerin on a microscope slide. Appendages were often cut off to avoid downstream contamination, as they may harbor debris including fungal propagules. The thalli were either cut into multiple smaller pieces using the sharp tip of the hypodermic needle (sensu Haelewaters et al. 2018a), crushed by pressing the tip of the hypodermic needle onto the thalli (sensu Weir and Blackwell 2001; Sundberg et al. 2018a), or a combination of both techniques was used. These pieces were then placed in a 0.2 ml PCR tube with 4 µl of phosphate-buffered saline (PBS). Next, 3 µl of prepared D2 buffer was pipetted against the inner wall of the PCR tube to prevent accidental removal of thallus fragments, and the tube was shortly centrifuged. The tube was then incubated at 65 °C for 30 min. After incubation, the tube was centrifuged for 20 min at 4000 RPM, and 3 µl of STOP solution was added. Again, to prevent accidental removal of thallus fragments, the STOP solution was pipetted against the inner wall of the PCR tube, followed by a brief centrifugation step. From here on, two different protocols were followed:
-
1. The normal REPLI-g protocol: in the tube, 29 µl Reaction Buffer; 9 µl ddH2O; and 2 µl REPLI-g sc DNA Polymerase was added, as indicated in the manufacturer’s instructions (Qiagen).
-
2. The halved REPLI-g protocol: in the tube, 14.5 µl Reaction Buffer; 4.5 µl ddH2O; and 1 µl REPLI-g sc DNA Polymerase was added.
We also made 32 extractions that we supplied with 21.75 µl Reaction Buffer; 6.75 µl ddH2O; and 1.5 µl REPLI-g sc DNA Polymerase (a total of 30 µl). Initially, we wanted to test whether the probability to successfully extract DNA differed significantly between these three protocols (normal protocol, halved protocol, and 30 µl protocol). After the first tests, we already observed that this probability did not differ significantly, and we decided to simplify and streamline our experimental setup by focusing on the two extremes (comparing the normal and the halve protocol).
To test the effect of BSA during the whole genome amplification step of the REPLI-g Single Cell Kit, we arbitrarily added 5 µl BSA (20 mg/ml, stock concentration).
After the addition of all reagents, the samples were incubated at 30 °C for 8 h. During this incubation step, whole-genome amplification (WGA) took place: the whole genome DNA in the tubes was amplified using Multiple Displacement Amplification (MDA) (see Discussion). After the WGA, the polymerase was inactivated at 65 °C for 3 min. DNA extractions were stored at -20 °C. DNA quantification was done using the Qubit 2.0 fluorescence spectrometer (Thermo Fisher Scientific, Waltham, MA) and measurements of the A260/A280 and A260/230 absorbance ratios were taken using a NanoDrop 2000 (Thermo Fisher Scientific).
The nrSSU, ITS, and nrLSU were amplified. Additionally, we attempted to amplify MCM7, TEF1, and mtSSU. All primer pairs used are given in Table 1. PCR reactions (25 µl total) consisted of 13.3 µl of RedExtract Taq polymerase (Sigma-Aldrich), 2.5 µl of each 10 µM primer, 5.45 µl of ddH2O, and 1 µl of DNA extract. Before pipetting 1 µl of DNA extract, the PCR tube was vortexed briefly. PCR conditions are listed in Table 2. PCR products were stored at -20 °C. We also attempted to amplify additional nrLSU, ITS, and TEF1 sequences of older, preserved extractions, which were made during former studies (Haelewaters 2018; De Weggheleire 2019; Haelewaters et al. 2019a).
Following PCR, gel electrophoresis was performed by loading the PCR products on a Tris–acetate-EDTA (TAE) 1% agarose gel at 130 V for 30 min. The gels were placed in an ethidium bromide solution for 15 min to visualize the PCR products. Purification of PCR products was done using 1.5 µl of Exo-FAP (0.5 µl Exonuclease I, 1 µl FAST Alkaline Phosphatase) (Thermo Fisher Scientific) per 10 µl of PCR product, at 37 °C for 15 min, followed by deactivation at 85 °C for 15 min. Purified PCR products were sequenced using an automated ABI 3730XL capillary sequencer at Macrogen (Amsterdam, The Netherlands). Sequencing primers were the same as the primers used to amplify the region of interest. Forward and reverse sequence reads were assembled and edited in Sequencher version 5.4.6 (Gene Codes Corporation, Ann Arbor, MI). Newly generated sequences were submitted to NCBI GenBank.
To assess the identity of the newly generated sequences, we constructed four alignments (nrSSU, ITS nrLSU, and TEF1) of these sequences, supplemented with a broad range of Laboulbeniomycetes sequences available on NCBI Genbank. We aligned nrSSU, nrLSU, and TEF1 sequences by locus with the G-INS-i strategy and ITS sequences with the E-INS-i strategy using the online version 7 of MAFFT (Kuraku et al. 2013; Katoh et al. 2019). Models for nucleotide substitution were selected for each partition with ModelFinder (Kalyaanamoorthy et al. 2017) according to the corrected Akaike Information Criterion (AICc). Maximum likelihood (ML) was inferred using IQ-TREE (Nguyen et al. 2015) under partitioned models (Chernomor et al. 2016). Ultrafast bootstrapping was performed with 1000 replicates (Hoang et al. 2018). Alignments and resulting phylogenetic trees are available on GitHub: https://github.com/dannyhaelewaters/teamlaboul/tree/main/molecular_laboulbenia_paper.
Statistical analyses
To test the difference in success/fail ratio between methods of preservation and whether the probability to successfully extract DNA differs between the halved and the normal REPLI-g protocol, we used generalized linear mixed models as implemented in the lme4 package in R (Bates et al. 2015; R Core Team 2021). We tested the assumptions for these models using the DHARMa package (Hartig 2022). Three methods of preservations were defined: ‘doubtful’ (specimens that were collected and preserved in 70% ethanol or dried and pinned for a prolonged time); ‘uncertain’ (specimens for which the preservation method was unknown); and ‘good’ (specimens that were collected and preserved in > 95% ethanol). The ‘uncertain’ group was created to include specimens from which the preservation method was unknown, to prevent assigning a specimen to a wrong group (to ‘doubtful’ or ‘good’). To compare the means of DNA concentration between the halved and the normal REPLI-g protocol and whether the addition of BSA increases the probability to successfully extract DNA, we used linear mixed models as implemented in the lme4 package in R (Bates et al. 2015; R Core Team 2021). Graphical representation of data was made using the ggplot2 package (Wickham 2016).
We included the species of Laboulbenia as a random effect, because success of DNA extraction might be correlated with species, as each species has a different degree of melanization. It is important to note that most species were only represented once or only a few times, so there was no balanced design. We assigned an extraction as ‘successful’ when there was at least a clear single band of the nrSSU amplicon on the gel or when the ITS or nrLSU sequence matched with Laboulbenia species. We used the emmeans package to obtain the Estimated Marginal Means (EMM) for each group and to compare the means between groups (Lenth et al. 2024). R code, output of the analyses, and additional exploratory figures can be found on GitHub: https://github.com/dannyhaelewaters/teamlaboul/tree/main/molecular_laboulbenia_paper.
Results
Comparison of preservation methods
A total of 196 extractions were made, of which 111 were successful. The probability of successfully extracting DNA of thalli from ‘doubtful’ specimens was significantly lower (π = 0.167, 95% confidence interval (CI): 0.0669–0.360) compared to ‘good’ specimens (π = 0.824, 95% CI: 0.6581–0.919) (Fig. 2A). There was no obvious pattern visible between DNA extraction success and the time (both in months and years) between collection and DNA extraction (Fig. 2B).
Comparison of extraction protocols and DNA concentration
There was no significant difference in probability of ‘success’ between protocols, except between the ‘No BSA and halved’ and ‘BSA and halved’ groups (Fig. 2C). The EMM of the probabilities and their CIs are given in Table 3. There was a significant difference in DNA concentration between the halved and normal REPLI-g protocols (p < 0.001) (Fig. 2D). The DNA concentration was significantly higher with the normal REPLI-g protocol (29 µl Reaction Buffer, 9 µl ddH2O, and 2 µl REPLI-g sc DNA Polymerase added, totaling 40 µl) (EMM = 63.4 µg/ml, 95% CI: 56.1–70.7) compared to the halved REPLI-g protocol (only 20 µl mixture added) (EMM = 37.2 µg/ml, 95% CI: 30.9–43.5). The A260/A280 and A260/A230 ratios of 45 DNA extractions (33 undiluted and 12 1/10 diluted) were measured. Only one of the undiluted DNA extractions had an A260/A280 ratio within the accepted range of ‘pure’ DNA (1.7–2.0). For the diluted DNA extractions, 10 extractions had an A260/A280 ratio of around 1.8, and 4 of these were also within the range of the ideal A260/A230 ratio (1.9–2.2).
Evaluation of primer pairs, PCR protocols, and sequence success
Of the 111 successful extractions, 16 were contaminated and we could only generate sequences of the nrSSU region with the Laboulbeniales-specific NLS1/NSL2 primers. The nrSSU sequences of these contaminated extractions were of good quality. The contaminants were identified using ITS or nrLSU sequences and are shown in Table 4. From the successful extractions, we generated 104 nrSSU, 64 ITS, 92 nrLSU, 1 MCM7, and 16 TEF1 sequences (Tables 5 and 6). We also generated 1 nrSSU, 1 ITS, 4 nrLSU, and 5 TEF1 sequences of preserved extractions from former studies (Haelewaters 2018; De Weggheleire 2019; Haelewaters et al. 2019a). All generated sequences were submitted to NCBI Genbank and their accession numbers are presented in Table 6.
The ‘CombSIL’ PCR protocol (acronym for Combination of nrSSU, ITS, nrLSU) is a combination of the PCR protocols for nrSSU, ITS, and nrLSU. These separate protocols all use the same annealing temperatures, and similar timings for each step. By combining them into one protocol, we can use the same PCR machine to amplify these different regions at the same time. When using the old protocol to amplify TEF1 (Table 2: TEF1), if bands were visible, there were often multiple bands. When increasing the annealing temperature, clear single bands were observed on the gel after staining (Table 2: TEF1 New).
Discussion
Comparison of preservation methods
The preservation method had a significant effect on the probability to successfully extract DNA of Laboulbeniales. Thalli stored in a doubtful way (in 70% ethanol or dried and pinned for a prolonged time) had a significantly lower probability of successful DNA extraction than thalli collected and stored in > 95% ethanol. This was already reported in a few studies on Laboulbeniales (Weir and Blackwell 2001; Haelewaters et al. 2015, 2019a). Anecdotally, we observed lower extraction success for thalli that were stored for more than 4–6 months (between collection and DNA extraction) in 70% ethanol. A more detailed and standardized experiment should be conducted to study the effect of different preservations methods on DNA extractions success of Laboulbeniales, including short-term versus longer-term preservation on 70% ethanol, 96% ethanol (expensive molecular grade and cheap denatured bio-ethanol), RNAlater, CTAB, isopropanol, and on -20 °C.
The REPLI-g Kit utilizes a Multiple Displacement Amplification (MDA) to amplify DNA during the WGA step, which involves random hexamer primers and phi29 polymerase (Long et al. 2020). MDA makes use of primers that randomly link to multiple sites of the DNA template and thus no target-specific primers are needed. Disadvantages of MDA include incomplete coverage and over-representation of certain fragments in the resulting DNA extract by chance (e.g., multi-copy regions). If DNA is fragmented, amplicons may be short of even absent, leading to incomplete amplification of the whole-genomic DNA. Subsequent PCR of a specific region is likely to fail due to fragmentation in the primer annealing sites or in the target amplicon. DNA fragmentation is expected in thalli when hosts were stored dried or in 70% ethanol for a prolonged time (> 4 months), as observed in other organisms (Bruns et al. 1990; Kates et al. 2021). For successful molecular work on Laboulbeniales, it is essential to perform the first steps of DNA extraction protocols immediately after collection, or to directly transfer of host specimens to > 95% ethanol, according to our current knowledge.
While this result was expected, it is crucial to emphasize the value of collections made by collaborators and entomologists as a critical resource for Laboulbeniales research (Haelewaters et al. 2015, 2021). Unfortunately, these collections are often inadequately preserved in 70% ethanol or dried and pinned for extended periods, as it is standard procedure in entomological research. For instance, in a 2022 collection of carabid beetles, only 12 out of 41 DNA extractions were successful, likely due to preservation in 70% ethanol for 6–12 months, even though the specimens were collected and processed within the past year (Fig. 2B). For other organisms it is known that their DNA will degrade in 70% ethanol after three months (Flournoy et al. 1996) to one year (Nagy 2010).
For the purpose of molecular work involving Sanger sequencing, we encourage collectors and collaborators to store their collections in > 95% ethanol, with a single ethanol refreshment (1 to 4 weeks after collection) to maintain the required concentration, as ethanol can extract water from host specimens and Laboulbeniales (Nagy 2010; Marquina et al. 2021). Taking proactive measures and collaborating with potential partners can maximize the utility of future collections for various entomological and mycological research purposes. Laboulbeniologists should seize these opportunities to investigate diverse aspects of these understudied insect-associated fungi, encompassing alpha taxonomy, ecology, evolution, and molecular research.
However, preliminary data suggest that preservation on > 95% ethanol and long-term preservation of DNA extracts (post-MDA) at -20 °C may not not useful for other applications such as whole genome sequencing (D. Haelewaters and C.A. Quandt, unpubl. data). Additionally, preservation in ethanol destroys the biofilms present on the thallus surfaces of Laboulbeniales, impeding the study of these unknown communities (Lubbers et al. 2022).
Comparison of extraction protocols and DNA concentrations
We had a high probability to successfully extract DNA of Laboulbenia species (0.824 < π < 0.873) if material was preserved correctly, compared to the 0% and 20% reported previously (Haelewaters et al. 2015). There was no difference in the probability of obtaining a successful DNA extraction between the different protocols (halved and normal) and with or without the addition of BSA. There was one exception: this probability is significantly lower for the ‘No BSA and halved’ group compared to the and ‘BSA and halved’ group, which is likely the result of the small sample size of the former group and an artifact of the number of samples used that were stored in a doubtful way. First of all, this means that we can use the halved protocol to make DNA extractions. This way, twice the number of extractions can be made with the same kit. A single DNA extraction using the REPLI-g Single Cell Kit costs between 24 and 30 euros. When halving the amount of resources per extraction, the costs per extraction are also halved. Yet, the extracts using the halved protocol are still more expensive than other widely used extraction protocols (Lickfeldt et al. 2002; Romanelli et al. 2014; Lahuf et al. 2019). The search for time- and cost-efficient and reliable DNA extraction methods for Laboulbeniales continues.
Secondly, the addition of BSA did not result in a significantly higher probability to obtain a successful extraction. There are multiple questions and considerations that arise with this result. It is harder to extract DNA of fungi that contain high concentrations of melanin(-like) molecules. These molecules confer rigidity and protection to the cells (Butler and Day 1998). In Laboulbeniales, melanin has been suggested to be the reason for the low DNA extraction successes (Haelewaters et al. 2015; Sundberg et al. 2018a). Thalli of Laboulbenia contain variable amounts of melanin(-like) molecules in their cell walls, but their concentration is unknown. Giambernardi et al. (1998) found that the addition of more than 0.5 µg melanin to a 25-µl assay results in the inhibition of Taq polymerase. Assuming that a thallus of an average species of Laboulbenia is a cylinder with a length of 350 µm, a diameter of 60 µm, and a density between 0.1 and 1.3 g/cm3 (Bakken and Olsen 1983), results in an estimated thallus weight between 0.099 and 1.29 µg. This would mean that the amount of melanin needed to inhibit PCR (0.5 µg/25 µl) is either higher than the lowest estimated weight of a thallus (0.099 µg) or more than a third of the highest estimated weight (1.29 µg), which seems unlikely. In other words, we think that the melanin content in the thallus cells is insufficient to hinder WGA or PCR. Consequently, the addition of BSA would have little impact on the probability of successful DNA extraction. It is possible, though, that other molecules such as phenols or proteins, potentially in combination with melanin(-like) molecules, could impede WGA or PCR.
The difference in DNA concentration between the halved and normal REPLI-g protocol was expected, as reagents used to amplify the genome are halved, only half of the amount of DNA can be amplified. There is still some variation in the DNA concentration within each group, which is probably because the amount of the REPLI-g mix pipetted into the PCR tubes slightly differed between samples due to small pipetting errors. The number of thalli used did not have a significant effect on the obtained DNA concentrations. This means that using one thallus is sufficient to make a successful DNA extraction using the REPLI-g Single Cell Kit. However, some species are very small, and using multiple thalli ensures that at least a few pieces of thalli will end up in the PBS solution when transferring them into the PCR tube.
We measured the A260/A280 and A260/A230 ratios to assess the purity of the extracted and amplified DNA. The measurements of all except one undiluted DNA extractions were not within the ideal absorbance ratio range for ‘pure’ DNA of 1.7–2.0 (for A260/A280) and 1.9–2.2 (for A260/A230). This is not surprising, as BSA was added to all these extractions, and other contaminants (other proteins, melanin, and phenols) are potentially present. Most diluted DNA extractions were, based on these absorbance ratios, more purified. As the DNA is diluted, the amount of contaminants is also diluted. There is an interest to sequence whole genomes of Laboulbeniomycetes to study the evolution of fungal genomes, population genetics, speciation patterns, parasite–host interactions, and the loss of hyphae for which good quality and ‘pure’ DNA extractions are needed. Haelewaters et al. (2020) sequenced the first Laboulbeniomycetes genome, of Herpomyces periplanetae. To further purify DNA extractions of Laboulbeniales, PCR inhibitor removal kits like the OneStep PCR Inhibitor Removal Kit (Zymo Research, Irvine, CA) (Hu et al. 2015; Vicente et al. 2019; Lubbers 2021) and SPRI bead cleaning (Beckman Coulter, Brea, CA) (B. Young and W. Van Caenegem, unpubl. data) can be used. This approach seems promising for future applications.
Evaluation of primer pairs, PCR protocols, and sequence success
Generating sequences of Laboulbenia species has generally been regarded as difficult (Haelewaters et al. 2015; Sundberg et al. 2018a). Here we generated the highest number of sequences of Laboulbenia in a single study, using the REPLI-g Single Cell Kit.
The NSL1/NSL2 primer pair works very well for species in the genus Laboulbenia. This pair specifically amplifies Laboulbeniomycetes DNA, thus it can also be used even when DNA extractions are contaminated (Haelewaters et al. 2015). While the nrSSU marker is generally not useful for species delimitation as it is a very conserved region (but see Paloi et al. 2022 for a discussion on group I introns), it can be used to distinguish higher taxa (genera, family, order).
The general fungal primers for ITS do not sufficiently amplify the whole ITS region of Laboulbenia species. Although we did have some success using ITS1f/ITS4 and ITS5/ITS4, there is a considerable difference with the amplification success of nrSSU and nrLSU region. At the end of the trial, we used the ITS3/ITS4 pair to amplify the ITS2 region. Surprisingly, this generated positive results for most taxa. It is likely that there is a primer mismatch at the primer annealing sites of ITS1f and ITS5, and, for some species, of ITS4 (e.g., Laboulbenia cristata, L. galeritae, and L. gyrinicola), as previously suggested (Liu et al. 2020). Bellemain et al. (2010) reported a bias for Ascomycota in eDNA studies when using ITS2, ITS3, and ITS4, while the forward primers ITS1, ITS1f, and ITS5 show a bias towards Basidiomycota. This might explain why we successfully amplified the ITS2 region, but not the ITS1 nor the whole ITS region for most species.
For the nrLSU region, the primer pair NL1/NL4 works well. The amplicon is around 300 base pairs shorter than the one amplified by LR0R/LR5. However, LR0R/LR5 did not work for most taxa of Laboulbenia and if it worked, multiple faint bands were often observed on the gel after visualization. Similar observations were made with extractions of Gloeandromyces and Hesperomyces (W. Van Caenegem and D. Haelewaters, unpubl. data). Multiple bands on the gel might indicate that the annealing temperature used (50 °C; Table 2) was suboptimal for this primer pair (Rychlik et al. 1990). Increasing the annealing temperature might overcome this problem, but there is also likely a primer mismatch in most species of Laboulbenia for one or both primers.
Although NL1/NL4 works very well for Laboulbenia species, it also amplifies the DNA of contaminants. Species of two genera (Malassezia and Priceomyces) were observed multiple times. Malassezia species (Ustilaginomycotina, Malasseziales) are basidiomycetous yeasts that live on the skin of humans and other vertebrates (Theelen et al. 2018). These species are probably present in our DNA extractions due to the many handlings needed during the extraction process (Fig. 1). Priceomyces vitoshaensis (Saccharomycotina, Saccharomycetales) is an ascomycetous yeast described from the carabid beetle Pterostichus melas (Crous et al. 2016). There is not much known about the ecology of this species, but other species of Priceomyces are found in eDNA studies of soil, on beetles, in beetle guts, and in insect frass (Kurtzman 2011; Groenewald et al. 2018; Kudo et al. 2019). The micromanipulation needed to remove thalli from the host, may increase contamination.
We could not successfully amplify the three other markers (MCM7, TEF1, and mtSSU). Some sequences of MCM7 and TEF1 were generated for Laboulbenia species, but we decided to not explore these markers and their primers further to save time, money, and resources. Recently, MCM7 was introduced as a secondary marker in Hesperomyces (Haelewaters et al. 2022; Van Caenegem et al. 2023b) and TEF1 is easily amplified for Gloeandromyces species (Liu et al. 2020; Van Caenegem et al. 2023a). Amplifying these regions is more challenging for species of Laboulbenia compared to our experiences with Gloeandromyces, Hesperomyces, and Nycteromyces (W. Van Caenegem and D. Haelewaters, unpubl. data). The first amplification trial of mtSSU yielded no sequences and we therefore made no further use of the MS1/MS2 primer pair. Sundberg et al. (2018b, 2018a) generated mtSSU sequences of Laboulbeniales and reported this region to be “the easiest to amplify and as well as sequence.” However, they used the primer pair mrSSU1/mrSSU3R and they mostly sequenced Coreomyces species, which in part may explain the difference in amplification success.
To overcome these challenges of contaminants, primer mismatches, and low amplification success of several loci, we propose to design Laboulbeniales- and more specifically Laboulbenia-specific primers, especially for the ITS and nrLSU regions, similar to the Laboulbeniomycetes-specific nrSSU primers (NSL1/NSL2) (Haelewaters et al. 2015) and the Hesperomyces-specific ITS primers (ITSHespL/ITSHespR) (Haelewaters et al. 2018a). These developments will not only help to overcome these aforementioned challenges within the genus Laboulbenia, but they will also contribute to the study of Laboulbeniales in general, as most species remain unsequenced, resulting in the underestimation of the true diversity of these understudied microfungi as well as a poor understanding of evolutionary relationships due to under-sampling.
Conclusions
Using the REPLI-g Single Cell Kit, we made 111 successful DNA extractions from Laboulbenia species. The preservation method had a major effect on the success of DNA extraction; specimens stored for extended periods in 70% ethanol or dried and pinned were found unsuitable for molecular work. For future research, we encourage entomologists, collaborators, and collectors to preserve infected host specimens in > 95% ethanol. Our findings revealed no significant differences in DNA extraction protocols, indicating that utilizing half the recommended amount is sufficient to successfully extract DNA—saving costs. Further, the addition of BSA did not significantly impact the probability of obtaining successful DNA extractions, suggesting that the melanin content in Laboulbenia species is negligible. We generated 104 nrSSU, 64 ITS, 92 nrLSU, 1 MCM7, and 16 TEF1 sequences. The limited success in amplifying protein-coding genes can likely be attributed to primer mismatches. Furthermore, a disparity in amplification success between ITS and the nrSSU and nrLSU regions was observed, likely stemming from primer mismatches. It is imperative to develop Laboulbeniales-specific ITS and nrLSU primers to tackle contaminations and improve amplification efficiency. These new insights do not only apply to the genus Laboulbenia; the vast majority Laboulbeniales species remain unsequenced. The study of Laboulbeniales will benefit from the molecular developments reported here. And while we keep critically evaluating and improving our methods, we hope that other research groups will be inspired to start molecular work with these peculiar microfungi.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the GitHub repository, https://github.com/dannyhaelewaters/teamlaboul/tree/main/molecular_laboulbenia_paper. Newly generated sequences were submitted to the National Center for Biotechnology Information (NCBI) GenBank database (https://www.ncbi.nlm.nih.gov/genbank/), under the accession numbers indicated in Table 6.
Abbreviations
- BSA:
-
Bovine serum albumin
- CI:
-
Confidence interval
- eDNA:
-
Environmental DNA
- EMM:
-
Estimated Marginal Means
- ITS:
-
Internal transcribed spacer
- MDA:
-
Multiple Displacement Amplification
- mtSSU:
-
Mitochondrial small subunit
- NCBI:
-
National Center for Biotechnology Information
- nrLSU:
-
Nuclear large subunit ribosomal RNA
- nrSSU:
-
Nuclear small subunit ribosomal RNA
- PBS:
-
Phosphate-buffered saline
- RPM:
-
Revolutions per minute
- sc:
-
Single cell
- WGA:
-
Whole-genome amplification
References
Bakken LR, Olsen RA (1983) Buoyant densities and dry-matter contents of microorganisms: conversion of a measured biovolume into biomass. Appl Environ Microbiol 45(4):1188–1195. https://doi.org/10.1128/aem.45.4.1188-1195.1983
Bates D, Maechler M, Bolker B, Walker S (2015) Fitting linear mixed-effects models using lme4. J Stat Softw 67(1):1–48. https://doi.org/10.18637/jss.v067.i01
Bellemain E, Carlsen T, Brochmann C, Coissac E, Taberlet P, Kauserud H (2010) ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases. BMC Microbiol 10:189. https://doi.org/10.1186/1471-2180-10-189
Blackwell M, Haelewaters D, Pfister DH (2020) Laboulbeniomycetes: Evolution, natural history, and Thaxter’s final word. Mycologia 112(6):1048–1059. https://doi.org/10.1080/00275514.2020.1718442
Bridge PD, Arora DK, Reddy CA, Elander RP (1998) Applications of PCR in mycology. CAB International, Wallingford
Bridge PD, Spooner BM, Roberts PJ (2005) The impact of molecular data in fungal systematics. Adv Bot Res 42:33–67. https://doi.org/10.1016/S0065-2296(05)42002-9
Bruns TD, Fogel R, Taylor JW (1990) Amplification and sequencing of DNA from fungal herbarium specimens. Mycologia 82(2):175–184. https://doi.org/10.2307/3759846
Butler MJ, Day AW (1998) Fungal melanins: a review. Can J Microbiol 44(12):1115–1136. https://doi.org/10.1139/w98-119
Cao B, Haelewaters D, Schoutteten N, Begerow D, Boekhout T, Giachini AJ, Gorjón SP, Gunde-Cimerman N, Hyde KD, Kemler M, Li GJ, Liu DM, Liu XZ, Nuytinck J, Papp V, Savchenko A, Savchenko K, Tedersoo L, Theelen B, Thines M, Tomšovský M, Toome-Heller M, Urón JP, Verbeken A, Vizzini A, Yurkov AM, Zamora JC, Zhao RL (2021) Delimiting species in Basidiomycota: a review. Fungal Divers 109:181–237. https://doi.org/10.1007/s13225-021-00479-5
Chernomor O, von Haeseler A, Minh BQ (2016) Terrace aware data structure for phylogenomic inference from supermatrices. Syst Biol 65(6):997–1008. https://doi.org/10.1093/sysbio/syw037
Crous PW, Wingfield MJ, Richardson DM, Leroux JJ, Strasberg D, Edwards J, Roets F, Hubka V, Taylor PWJ, Heykoop M, Martín MP, Moreno G, Sutton DA, Wiederhold NP, Barnes CW, Carlavilla JR, Gené J, Giraldo A, Guarnaccia V, Guarro J, Hernández-Restrepo M, Kolařík M, Manjón JL, Pascoe IG, Popov ES, Sandoval-Denis M, Woudenberg JHC, Acharya K, Alexandrova AV, Alvarado P, Barbosa RN, Baseia IG, Blanchette RA, Boekhout T, Burgess TI, Cano-Lira JF, Čmoková A, Dimitrov RA, Dyakov MYu, Dueñas M, Dutta AK, Esteve-Raventós F, Fedosova AG, Fournier J, Gamboa P, Gouliamova DE, Grebenc T, Groenewald M, Hanse B, Hardy GEStJ, Held BW, Jurjević Ž, Kaewgrajang T, Latha KPD, Lombard L, Luangsa-ard JJ, Lysková P, Mallátová N, Manimohan P, Miller AN, Mirabolfathy M, Morozova OV, Obodai M, Oliveira NT, Ordóñez ME, Otto EC, Paloi S, Peterson SW, Phosri C, Roux J, Salazar WA, Sánchez A, Sarria GA, Shin HD, Silva BDB, Silva GA, Smith MTh, Souza-Motta CM, Stchigel AM, Stoilova-Disheva MM, Sulzbacher MA, Telleria MT, Toapanta C, Traba JM, Valenzuela-Lopez N, Watling R, Groenewald JZ, (2016) Fungal Planet description sheets: 400–468. Persoonia 36:316–458. https://doi.org/10.3767/003158516X692185
De Kesel A, Gerstmans C, Haelewaters D (2020) Catalogue of the Laboulbeniomycetes of Belgium. Sterbeeckia 36:3–143
De Weggheleire S (2019) Biodiversity of the genus Laboulbenia: are there cryptic species? MSc thesis. University of Antwerp, Antwerp, pp 1–48
Eckhart L, Bach J, Ban J, Tschachler E (2000) Melanin binds reversibly to thermostable DNA polymerase and inhibits its activity. Biochem Biophys Res Commun 271(3):726–730. https://doi.org/10.1006/bbrc.2000.2716
Flournoy LE, Adam RP, Pandy RN (1996) Interim and archival preservation of plant specimens in alcohols for DNA studies. Biotechniques 20(4):657–660. https://doi.org/10.2144/19962004657
Gardes M, Bruns TD (1993) ITS primers with enhanced specificity for basidiomycetes - application to the identification of mycorrhizae and rusts. Mol Ecol 2(2):113–118. https://doi.org/10.1111/j.1365-294X.1993.tb00005.x
Giambernardi TA, Rodeck U, Klebe RJ (1998) Bovine serum albumin reverses inhibition of RT-PCR by melanin. Biotechniques 25(4):564–566. https://doi.org/10.2144/98254bm03
Gibson JL (2012) Ultrastructural observations of the Laboulbeniales. Master thesis. State University of New York, Syracuse, pp 1–70
Goldmann L, Weir A (2012) Position specificity in Chitonomyces (Ascomycota, Laboulbeniomycetes) on Laccophilus (Coleoptera, Dytiscidae): a molecular approach resolves a century-old debate. Mycologia 104(5):1143–1158. https://doi.org/10.3852/11-358
Goldmann L, Weir A (2018) Molecular phylogeny of the Laboulbeniomycetes (Ascomycota). Fungal Biol 122(2–3):87–100. https://doi.org/10.1016/j.funbio.2017.11.004
Goldmann L, Weir A, Rossi W (2013) Molecular analysis reveals two new dimorphic species of Hesperomyces (Ascomycota, Laboulbeniomycetes) parasitic on the ladybird Coleomegilla maculata (Coleoptera, Coccinellidae). Fungal Biol 117(11–12):807–813. https://doi.org/10.1016/j.funbio.2013.10.004
Groenewald M, Lombard L, de Vries M, Lopez AG, Smith M, Crous PW (2018) Diversity of yeast species from Dutch garden soil and the description of six novel Ascomycetes. FEMS Yeast Res 18(7):foy076. https://doi.org/10.1093/femsyr/foy076
Haelewaters D (2018) Studies of the Laboulbeniomycetes: Diversity, evolution, and patterns of speciation. Ph.D. dissertation. Harvard University, Cambridge, Massachusetts, pp 1–287
Haelewaters D, De Kesel A (2020) Checklist of thallus-forming Laboulbeniomycetes from Belgium and the Netherlands, including Hesperomyces halyziae and Laboulbenia quarantenae spp. nov. MycoKeys 71:23–86. https://doi.org/10.3897/mycokeys.71.53421
Haelewaters D, De Kesel A, Gorczak M, Bao K, Gort G, Zhao SY, Pfister DH (2019a) Laboulbeniales (Ascomycota) of the Boston Harbor Islands II (and other localities): species parasitizing Carabidae, and the Laboulbenia flagellata species complex. Northeast Nat 25(Sp9):110–149. https://doi.org/10.1656/045.025.s906
Haelewaters D, De Kesel A, Pfister DH (2018a) Integrative taxonomy reveals hidden species within a common fungal parasite of ladybirds. Sci Rep 8(1):15966. https://doi.org/10.1038/s41598-018-34319-5
Haelewaters D, Gorczak M, Kaishian P, De Kesel A, Blackwell M (2021) Laboulbeniomycetes, enigmatic fungi with a turbulent taxonomic history. In: Encyclopedia of mycology, vol 1. Elsevier, Oxford, pp 263–283. https://doi.org/10.1016/B978-0-12-819990-9.00049-4
Haelewaters D, Gorczak M, Pfliegler WP, Tartally A, Tischer M, Wrzosek M, Pfister DH (2015) Bringing Laboulbeniales into the 21st century: enhanced techniques for extraction and PCR amplification of DNA from minute ectoparasitic fungi. IMA Fungus 6(2):363–372. https://doi.org/10.5598/imafungus.2015.06.02.08
Haelewaters D, Matthews TJ, Wayman JP, Cazabonne J, Heyman F, Quandt CA, Martin TE (2024) Laboulbeniomycetes as a case study for biodiversity shortfalls in poorly studied groups. J Biogeogr 51(1):29–39. https://doi.org/10.1111/jbi.14725
Haelewaters D, Okrasińska A, Gorczak M, Pfister DH (2020) Draft genome sequence of the globally distributed cockroach-infecting fungus Herpomyces periplanetae strain D. Haelew. 1187d. Microbiol Resour Announc 9(6):e01458–19. https://doi.org/10.1128/MRA.01458-19
Haelewaters D, Page RA, Pfister DH (2018b) Laboulbeniales hyperparasites (Fungi, Ascomycota) of bat flies: independent origins and host associations. Ecol Evol 8(16):8396–8418. https://doi.org/10.1002/ece3.4359
Haelewaters D, Pfister DH (2019) Morphological species of Gloeandromyces (Ascomycota, Laboulbeniales) evaluated using single-locus species delimitation methods. Fungal Syst Evol 3:19–34. https://doi.org/10.3114/fuse.2019.03.03
Haelewaters D, Pfliegler WP, Gorczak M, Pfister DH (2019b) Birth of an order: Comprehensive molecular phylogenetic study excludes Herpomyces (Fungi, Laboulbeniomycetes) from Laboulbeniales. Mol Phylogenet Evol 133:286–301. https://doi.org/10.1016/j.ympev.2019.01.007
Haelewaters D, Van Caenegem W, De Kesel A (2022) Hesperomyces harmoniae, a new name for a common ectoparasitic fungus on the invasive alien ladybird Harmonia axyridis. Sydowia 75:53–74. https://doi.org/10.12905/0380.sydowia75-2023-0053
Hartig F (2022) DHARMa: Residual Diagnostics for Hierarchical (Multi-Level / Mixed) Regression Models. https://cran.r-project.org/package=DHARMa
Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS (2018) UFBoot2: Improving the ultrafast bootstrap approximation. Mol Biol Evol 35(2):518–522. https://doi.org/10.1093/molbev/msx281
Hopple JS (1994) Phylogenetic investigations in the genus Coprinus based on morphological and molecular characters. Ph.D. dissertation. Duke University, Durham, North Carolina, pp 1–402
Hu Q, Liu Y, Yi S, Huang D (2015) A comparison of four methods for PCR inhibitor removal. Forensic Sci Int Genet 16:94–97. https://doi.org/10.1016/j.fsigen.2014.12.001
Kalyaanamoorthy S, Minh B, Wong T, von Haeseler A, Jermiin LS (2017) ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat Methods 14(6):587–589. https://doi.org/10.1038/nmeth.4285
Kates HR, Doby JR, Siniscalchi CM et al (2021) The effects of herbarium specimen characteristics on short-read NGS sequencing success in nearly 8000 specimens: Old, degraded samples have lower DNA yields but consistent sequencing success. Front Plant Sci 12:669064. https://doi.org/10.3389/fpls.2021.669064
Katoh K, Rozewicki J, Yamada KD (2019) MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform 20(4):1160–1166. https://doi.org/10.1093/bib/bbx108
Kudo R, Masuya H, Endoh R, Kikuchi T, Ikeda H (2019) Gut bacterial and fungal communities in ground-dwelling beetles are associated with host food habit and habitat. ISME J 13(3):676–685. https://doi.org/10.1038/s41396-018-0298-3
Kuraku S, Zmasek CM, Nishimura O, Katoh K (2013) aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity. Nucleic Acids Res 41:22–28. https://doi.org/10.1093/nar/gkt389
Kurtzman CP (2011) Chapter 59 - Priceomyces M. Suzuki and Kurtzman 2010. In: Kurtzman CP, Fell JW, Boekhout T (eds) The yeasts, 5th edn. Elsevier, London, pp 719–724
Kurtzman CP, Robnett CJ (1997) Identification of clinically important ascomycetous yeasts based on nucleotide divergence in the 5’ end of the large-subunit (26S) ribosomal DNA gene. J Clin Microbiol 35:1216–1223. https://doi.org/10.1128/jcm.35.5.1216-1223.1997
Lahuf AA, Jaafar OH, Hameed ZL (2019) A simple, rapid, safe and low-cost method to extract DNA from phytopathogenic fungi. Asian J Agric Biol 7(2):197–203
Landvik S, Egger KN, Schumacher T (1997) Towards a subordinal classification of the Pezizales (Ascomycota): phylogenetic analyses of SSU rDNA sequences. Nord J Bot 17:403–418. https://doi.org/10.1111/j.1756-1051.1997.tb00337.x
Lenth RV, Bolker B, Buerkner P, Giné-Vázquez I, Herve M, Jung M, Love J, Miguez F, Riebl H, Singmann H (2024) emmeans: Estimated Marginal Means, aka Least-Squares Means. https://cran.r-project.org/package=emmeans
Lickfeldt D, Hofmann NE, Jones JD, Hamblin AM, Voigt TB (2002) Comparing three DNA extraction procedures for cost, efficiency, and DNA yield. HortScience 37(5):822–825. https://doi.org/10.21273/HORTSCI.37.5.822
Liu J, Haelewaters D, Pfliegler WP, Page RA, Dick CW, Aime MC (2020) A new species of Gloeandromyces from Ecuador and Panama revealed by morphology and phylogenetic reconstruction, with a discussion of secondary barcodes in Laboulbeniomycetes taxonomy. Mycologia 112(6):1192–1202. https://doi.org/10.1080/00275514.2020.1781496
Long N, Qiao Y, Xu Z, Tu J, Lu Z (2020) Recent advances and application in whole-genome multiple displacement amplification. Quant Biol 8(4):279–294. https://doi.org/10.1007/s40484-020-0217-2
Lubbers M (2021) Peculiar parasites: on the ecology and genetics of Laboulbeniales. Master thesis. Leiden University, Leiden, pp 1–28
Lubbers M, Lamers G, De Kesel A, Nedvěd O, Schilthuizen M, Haelewaters D (2022) Bacterial biofilms on thalli of Laboulbeniales: a community uncovered. Sydowia 74:335–342. https://doi.org/10.12905/0380.sydowia74-2022-0335
Maharachchikumbura SSN, Chen Y, Ariyawansa HA, Hyde KD, Haelewaters D, Perera RH, Samarakoon MC, Wanasinghe DN, Bustamante DE, Liu JK, Lawrence DP, Cheewangkoon R, Stadler M (2021) Integrative approaches for species delimitation in Ascomycota. Fungal Divers 109:155–179. https://doi.org/10.1007/s13225-021-00486-6
Majewski T (1994) The Laboulbeniales of Poland. Polish Botanical Studies 7:3–466
Marquina D, Buczek M, Ronquist F, Łukasik P (2021) The effect of ethanol concentration on the morphological and molecular preservation of insects for biodiversity studies. PeerJ 9:e10799. https://doi.org/10.7717/peerj.10799
Miadlikowska J, Lutzoni F (2000) Phylogenetic revision of the genus Peltigera (lichen-forming Ascomycota) based on morphological, chemical, and large subunit nuclear ribosomal DNA data. Int J Plant Sci 161(6):925–958. https://doi.org/10.1086/317568
Muilwijk J, Felix R, Dekoninck W, Bleich O (2015) De loopkevers van Nederland en België (Carabidae). Entomol Tabellen 9:1–215
Nagy ZT (2010) A hands-on overview of tissue preservation methods for molecular genetic analyses. Org Divers Evol 10:91–105. https://doi.org/10.1007/s13127-010-0012-4
Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ (2015) IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32(2):268–274. https://doi.org/10.1093/molbev/msu300
Paloi S, Luangsa-ard JJ, Mhuantong W, Stadlet M, Kobmoo N (2022) Intragenomic variation in nuclear ribosomal markers and its implication in species delimitation, identification and barcoding in fungi. Fungal Biol Rev 42:1–33. https://doi.org/10.1016/j.fbr.2022.04.002
Perkins DD, Raju NB (1986) Neurospora discreta, a new heterothallic species defined by its crossing behavior. Exp Mycol 10(4):323–338. https://doi.org/10.1016/0147-5975(86)90019-8
Pincus DH, Orenga S, Chatellier S (2007) Yeast identification—past, present, and future methods. Med Mycol 45(2):97–121. https://doi.org/10.1080/13693780601059936
R Core Team (2021) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. http://www.R-project.org
Romanelli AM, Fu J, Herrera ML, Wickes BL (2014) A universal DNA extraction and PCR amplification method for fungal rDNA sequence-based identification. Mycoses 57(10):612–622. https://doi.org/10.1111/myc.12208
Rychlik W, Spencer WJ, Rhoads RE (1990) Optimization of the annealing temperature for DNA amplification in vitro. Nucleic Acids Res 18(21):6409–6412. https://doi.org/10.1093/nar/18.21.6409
Santamaria S, Balazuc J, Tavares II (1991) Distribution of the European Laboulbeniales (Fungi, Ascomycotina): an annotated list of species. Treb Inst Bot Barcelona 14:1–123
Santamaria S, Pedersen J (2021) Laboulbeniomycetes (Fungi, Ascomycota) of Denmark. Eur J Taxon 781:1–425. https://doi.org/10.5852/ejt.2021.781.1583
Schmitt I, Crespo A, Divakar PK, Fankhauser JD, Herman-Sackett E, Kalb K, Nelsen MP, Nelson NA, Rivas-Plata E, Shimp AD, Widhelm T, Lumbsch HT (2009) New primers for promising single-copy genes in fungal phylogenetics and systematics. Persoonia 23:35–40. https://doi.org/10.3767/003158509X470602
Schoch C, Seifert K, Huhndorf SM, Robert V, Spouge JL, Levesque CA, Chen W, Fungal Barcoding Consortium (2012) Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc Natl Acad Sci U S A 109(16):6241–6246. https://doi.org/10.1073/pnas.1117018109
Stielow B, Lévesque CA, Seifert K, Meyer W, Irinyi L, Smits D, Renfurm R, Verkley GJM, Groenewald M, Chaduli D, Lomascolo A, Welti S, Lesage-Meessen L, Favel A, Al-Hatmi AMS, Damm U, Yilmaz N, Houbraken J, Lombard L, Quaedvlieg W, Binder M, Vaas LAI, Vu D, Yurkov A, Begerow D, Roehl O, Guerreiro M, Fonseca A, Samerpitak K, van Diepeningen AD, Dolatabadi S, Moreno LF, Casaregola S, Mallet S, Jacques N, Roscini L, Egidi E, Bizet C, Garcia-Hermoso D, Martín MP, Deng S, Groenewald JZ, Boekhout T, de Beer ZW, Barnes I, Duong TA, Wingfield MJ, de Hoog GS, Crous PW, Lewis CT, Hambleton S, Moussa TAA, Al-Zahrani HS, Almaghrabi OA, Louis-Seize G, Assabgui R, McCormick W, Omer G, Dukik K, Cardinali G, Eberhardt U, de Vries M, Robert V (2015) One fungus, which genes? Development and assessment of universal primers for potential secondary fungal DNA barcodes. Persoonia 35:242–263. https://doi.org/10.3767/003158515X689135
Sundberg H, Ekman S, Kruys Å (2018a) A crush on small fungi: An efficient and quick method for obtaining DNA from minute ascomycetes. Methods Ecol Evol 9(1):148–158. https://doi.org/10.1111/2041-210X.12850
Sundberg H, Kruys Å, Bergsten J, Ekman S (2018b) Position specificity in the genus Coreomyces (Laboulbeniomycetes, Ascomycota). Fungal Syst Evol 1:217–228. https://doi.org/10.3114/fuse.2018.01.09
Thaxter R (1896) Contributions towards a monograph of the Laboulbeniaceae: Part I. Mem Am Acad Arts Sci N S 12:189–429
Theelen B, Cafarchia C, Gaitanis G, Bassukas ID, Boekhout T, Dawson TL Jr (2018) Malassezia ecology, pathophysiology, and treatment. Med Mycol 56:S10–S25. https://doi.org/10.1093/mmy/myx134
Van Caenegem W, Blondelle A, Dumolein I, Santamaria B, Dick CW, Hiller T, Liu J, Quandt CA, Villarreal Saucedo RV, Verbeken A, Haelewaters D (2023a) Five new species of Gloeandromyces (Fungi, Laboulbeniales) from tropical American bat flies (Diptera, Streblidae), revealed by morphology and phylogenetic reconstruction. Mycologia 115(5):714–737. https://doi.org/10.1080/00275514.2023.2230114
Van Caenegem W, Ceryngier P, Romanowski J, Pfister DH, Haelewaters D (2023b) Hesperomyces (Fungi, Ascomycota) associated with Hyperaspis ladybirds (Coleoptera, Coccinellidae): rethinking host specificity. Front Fungal Biol 3:1040102. https://doi.org/10.3389/ffunb.2022.1040102
Vicente ALSA, Bianchini RA, Laus AC, Macedo G, Reis RM, Vasquez VL (2019) Comparison of protocols for removal of melanin from genomic DNA to optimize PCR amplification of DNA purified from highly pigmented lesions. Histol Histopathol 34(9):1089–1096. https://doi.org/10.14670/hh-18-112
Vilgalys R, Hester M (1990) Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species. J Bacteriol 172(8):4238–4246. https://doi.org/10.1128/jb.172.8.4238-4246.1990
Walker MJ, Dorrestein A, Camacho JJ, Hiller T, Meckler LA, Silas KA, Haelewaters D (2018) A tripartite survey of hyperparasitic fungi associated with ectoparasitic flies on bats (Mammalia: Chiroptera) in a neotropical cloud forest in Panama. Parasite 25:19. https://doi.org/10.1051/parasite/2018017
Weir A, Blackwell M (2001) Extraction and PCR amplification of DNA from minute ectoparasitic fungi. Mycologia 93(4):802–806. https://doi.org/10.1080/00275514.2001.12063212
White TJ, Bruns TD, Lee SB, Taylor JW (1990) Analysis of phylogenetic relationships by amplification and direct sequencing of ribosomal RNA genes. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ (eds) PCR protocols: a guide to methods and applications. Academic Press, San Diego, pp 315–322
Wickham H (2016) ggplot2: elegant graphics for data analysis. Springer, New York, pp 1–213. https://doi.org/10.1007/978-0-387-98141-3
Acknowledgements
The authors are thankful to the members of #TeamLaboul for constructive feedback on previous versions of this manuscript and to two anonymous reviewers for constructive comments and suggestions. Benjamin Young and C. Alisha Quandt (University of Colorado Boulder) are thanked for sharing unpublished data. We would like to thank the many collaborators who collected and identified most infected host specimens for this study: Yusseff P. Aguirre, Aimée Blondelle, Barend de Boer, André De Kesel, Wouter Dekoninck, Albert Dees, Carl W. Dick, Colin Easton, Ron Felix, Cyrille Gerstmans, Felix Heyman, Helen Law, Lieke Moereels, Remco Ploeg, Alisha C. Quandt, Rob Rossel, Menno Schilthuizen, Henk Spijkers, Kobejoren Van Vooren, Paul S. van Wielink, Jan-Willem van Zuijlen, Ally M. Willis, and James D. Young.
Adherence to national and international regulations
Not applicable.
Funding
Danny Haelewaters is supported by the Research Foundation – Flanders (Senior Postdoctoral Fellowship 1206024N).
Author information
Authors and Affiliations
Contributions
Conceptualization, Methodology, and Writing—Review and Editing: W.V.C. and D.H. Formal Analysis, Investigation, Writing—Original Draft, and Visualization: W.V.C. Resources, Supervision, and Funding Acquisition: D.H. Both authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
Danny Haelewaters is an Associate Editor of IMA Fungus. The authors declare no other competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Van Caenegem, W., Haelewaters, D. New insights into the DNA extraction and PCR amplification of minute ascomycetes in the genus Laboulbenia (Pezizomycotina, Laboulbeniales). IMA Fungus 15, 14 (2024). https://doi.org/10.1186/s43008-024-00146-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s43008-024-00146-9