ABSTRACT
The use of microsatellites as highly polymorphic DNA markers for the typing of isolates of Aspergillus fumigatus was investigated. Four CA repeats were selected by screening an A. fumigatus DNA library with a (CA)10 oligonucleotide. Primers flanking these CA repeats were designed to amplify each locus. One primer of each pair was labeled with a fluorophore, and the PCR products were analyzed with an automatic sequencer and the GeneScan software. For each primer set and for a given isolate, one band was detected and was assigned to an allele because A. fumigatus is haploid. With 50 clinical isolates, 50 environmental isolates, and 2 reference strains we obtained 12, 11, 10, and 23 different alleles for the four CA microsatellites, respectively (discriminatory power, 0.994). The results were identical by whatever DNA extraction technique was used. Interestingly, no clustering between environmental and clinical isolates was observed, suggesting that every isolate is potentially pathogenic. Microsatellite markers appear suitable for use in large epidemiological studies of invasive aspergillosis.
Invasive aspergillosis due toAspergillus fumigatus has become the leading cause of death in immunocompromised patients such as allogeneic bone marrow transplant recipients. Invasive aspergillosis can be contracted during hospitalization, especially by patients during their first month after bone marrow transplantation (26). Because inhalation of spores into the respiratory tract is the usual route of contamination, prevention can be achieved by using rooms with laminar air flow, but all patients at risk do not benefit from such conditions. Therefore, tracking the sources of contamination is an essential preventive measure. Molecular tools provide a means of comparing isolates from patients and from their environment to determine the source of contamination.
DNA fingerprinting techniques include the direct detection of restriction fragment length polymorphisms (RFLPs) (6, 8) and the detection of RFLPs by Southern hybridization with the M13 phage (1), ribosomal probes (22), or a retrotransposon-like element (19). The last probe has proved to be highly discriminatory for A. fumigatus isolates (7). However, RFLP analysis is a very time-consuming technique. Large amounts of DNA are needed, and DNA electrophoresis, blotting, and probing take more than 5 days for a limited number of isolates.
In contrast, the PCR-based methods are attractive because of their rapidity. Randomly amplified polymorphic DNA (RAPD) analysis has been applied to A. fumigatus. The primers either were short and arbitrarily chosen (2, 13, 15, 21) or consisted of repetitive motifs known to detect variable DNA sequences in lower eukaryotes and prokaryotes (13, 25). However, RAPD analysis exhibits a low level of reproducibility because of the low-stringency conditions used in the PCR, and these conditions lead to mismatched pairings (28). Therefore, the patterns may be complex and hardly comparable between laboratories. More recently, a high-stringency PCR technique with primers specific for RAPD products has been proposed (16). This strategy requires the sequencing of RAPD products in a first step and overcomes the shortcomings of low-stringency conditions, but several sets of primers are needed to obtain high discriminatory power. However, whatever RAPD technique is used, the nature of RAPD polymorphisms is too poorly understood to be useful without limitations in phylogenetic studies (3).
Microsatellites represent another class of genetic markers that have not yet been used to differentiate A. fumigatusstrains. Microsatellites are short tandem repeats of two to six nucleotides that are known to be highly polymorphic as well as numerous and spread equally throughout the human genome (27). Microsatellites have also been found in lower eukaryotic organisms including fungi (4, 9). The polymorphism of microsatellites can be evaluated by PCR, and precise allele sizing can be achieved with fluorescent primers and an automatic sequencer. Therefore, we investigated whether microsatellites could be found in A. fumigatus and used as DNA markers to differentiate isolates of this fungal species.
MATERIALS AND METHODS
Screening of libraries, hybridization, and sequencing.Total DNA of A. fumigatus IP 2279.94 (Pasteur Institute, Paris, France) was partially digested with either TaqI orMboI and was used to create two genomic libraries. The 150- to 1,000-bp fragments were inserted into the corresponding enzymatic site of bacteriophage M13mp18 (Appligene Oncor, Illkirch, France) before transformation of Epicurian Coli XL1-Blue MRF′ electroporation-competent cells (Stratagene, La Jolla, Calif.). Bacteriophage plaques were transferred onto Hybond-N+filters (Amersham, Les Ulis, France) and were probed with a (CA)10 oligonucleotide. The oligonucleotide was labeled with fluorescein-11-dUTP by using the ECL 3′-oligolabelling system (Amersham) and was diluted in hybridization buffer to a final concentration of 10 ng per ml. After hybridization at 37°C for 2 h, the filters were washed twice at room temperature for 5 min each time in 5× saline sodium citrate (SSC; 1× SSC is 0.15 M NaCl plus 0.015 M sodium citrate)–0.1% sodium dodecyl sulfate (SDS) and twice at 37°C for 15 min each time in 1× SSC–0.1% SDS. Incubation with anti-fluorescein-horseradish peroxidase conjugate, detection with the ECL detection system (Amersham), and exposure to Hyperfilm-ECL (Amersham) were carried out according to the manufacturer’s procedures. Clones showing hybridization signals were picked. Single-stranded DNAs from the bacteriophage were extracted, and automated sequencing was performed with the ABI dye primer kit (Perkin-Elmer, Courtaboeuf, France).
PCR amplification and analysis.For the clones containing a (CA)n microsatellite, primer pairs complementary to the flanking sequences of each microsatellite were designed. One primer of each set was labeled with a fluorescent dye, either 6-carboxyfluorescein (6-FAM) or 4,7,2′,4′,5′,7′-hexachloro-6-carboxyfluorescein (HEX) (Oligo-Express, Paris, France), for detection with an automated DNA sequencer. PCR amplification was performed in a 20-μl volume containing 1.5 mM MgCl2, 10 mM Tris-HCl (pH 9.0), and 50 mM KCl and with the forward and reverse primers at concentrations of 100 nM each, the deoxynucleoside triphosphates at concentrations of 100 μM each, 0.5 U of Taq DNA polymerase (Pharmacia Biotech, Orsay, France), 5% (vol/vol) dimethyl sulfoxide, and 50 ng of DNA template. Amplification was carried out in a Perkin-Elmer Cetus system 480 thermocycler, with denaturation for 5 min at 94°C, 30 times cycles of 30 s at 94°C, 30 s at 59°C, and 30 s at 72°C, and a final extension step at 72°C for 30 min. The PCR products were diluted 1/10 in water, and 1 μl of each was run on a 36-cm-long acrylamide-urea gel (4.25% polyacrylamide, 8.3 M urea, and 1× Tris-borate-EDTA for 2 h at 3,000 V and 51°C). TheN,N,N′,N′-tetramethyl-6-carboxyrhodamine (TAMRA)-labeled GeneScan size standard (Perkin-Elmer) was loaded into each well along with the PCR products. Signals were read with an automatic sequencer (ABI 377; Applied BioSystems), and the data were stored and analyzed with GeneScan software (version 2.0.2; Perkin-Elmer) by the local Southern sizing method.
A. fumigatus isolates.One hundred isolates (50 environmental isolates and 50 clinical isolates) were studied. The isolates were identified on the basis of culture characteristics and the morphologies of the conidiophores and conidia. The environmental isolates were collected from three Parisian hospitals over a 2-year period and were epidemiologically unrelated to the clinical isolates (7). DNA was extracted from a single colony by a phenol-chloroform technique at the Pasteur Institute (10). Additionally, 20 of the same isolates were also sent to our laboratory as a culture on a 2% malt slant, and a new DNA extraction was performed by a cetyltrimethylammonium bromide technique (18).
Of the 50 clinical isolates, 39 were obtained from different patients with invasive aspergillosis followed at the Henri Mondor Hospital over a 6-year period, 3 were recovered from patients with invasive aspergillosis hospitalized in Germany, and 8 were from patients with cystic fibrosis followed in Parisian hospitals. Two reference strains originally isolated from two patients (strains IP 2279.94 and CBS 143.89) were also studied. DNA was extracted by a cetyltrimethylammonium bromide technique. Twenty of the clinical isolates were also tested without any specific DNA extraction. About 104 to 105 spores were suspended in 50 μl of distilled water, frozen at −80°C for 10 min, thawed, and spun, and 5 μl of the supernatant was taken for PCR.
Discriminatory power and cluster analysis.The ability of microsatellite polymorphism to discriminate between isolates was assessed by using Simpson’s index of diversity:
RESULTS
Among 17,000 clones of the two A. fumigatuslibraries screened, 9 clones hybridized with the (CA)10 oligonucleotide probe and were sequenced. One of them was rejected as a TaqI cloning site and raised the possibility that recombination had occurred during the ligation process. Four other clones with less than eight CA repeats were also rejected because no polymorphism was detected upon a preliminary PCR screening with 10 isolates. Therefore, only four microsatellites were retained for further analysis, and there are henceforth referred to as microsatellites A, B, C, and D (Table 1). A search of the sequences in the GenBank and EMBL databases with the Blastn search program did not yield any sequence homology between the microsatellite flanking sequences and prokaryotic or eukaryotic sequences.
Features of the four polymorphic microsatellite sequences of A. fumigatus upon analysis of 100 isolates and 2 reference strains
PCR products consisted mainly of a single band, but additional bands shorter by 2 to 4 bp were seen for some alleles (Fig.1). These extra bands were due to the slippage of the Taq DNA polymerase and were always less intense than the longest band. Only the latter band was used for scoring and was assigned to an allele since A. fumigatus is thought to be haploid. For one isolate, we found two completely different peaks. We therefore performed a monospore culture and reamplified the DNA obtained for 10 different colonies. Each colony gave an unique peak. This isolate was actually a mixture of twoA. fumigatus isolates and was therefore excluded from the analysis. For the 102 isolates tested in the study, each PCR amplification was repeated at least twice with the same DNA template, and identical profiles were obtained. In addition, identical results were obtained regardless of the DNA extraction technique used. Besides, the amplifications were specific to A. fumigatusbecause no band was observed upon amplification of otherAspergillus species: A. flavus IP 597.69,A. niger IP 1431.83, A. terreus IP 1136.76, and A. nidulans IP 17.60.
GeneScan analysis of PCR profiles obtained with two independent isolates (isolates A560 and A422). PCRs were performed with a 6-FAM-labeled primer (loci A, B, and C) or a HEX-labeled primer (locus D), and the PCR products for each microsatellite locus (A, B, C, and D) were run in an acrylamide-urea gel. Bands produced fluorescent peaks, and their molecular sizes were automatically determined by comparison to the TAMRA-labeled GeneScan internal size standards loaded in each well (data not shown). The numbers refer to the sizes (in base pairs) of the PCR products by considering the fluorescent peak with the maximum height.
Upon the analysis of the 2 reference strains and the 100 unrelated isolates, 12, 11, 10, and 23 alleles were detected for microsatellites A, B, C, and D, respectively (Fig.2). The combination of the four markers led to 80 different allelic types, which corresponded to a discriminatory power of 0.994 (Table 2).
Allele size distributions of A. fumigatus isolates at microsatellite loci A (A), B (B), C (C), and D (D) upon analysis of 50 environmental isolates (solid bars) and 52 clinical isolates including 2 reference strains (striped bars).
Discriminatory power of microsatellite markers for the 100 A. fumigatus isolates and the 2 reference strains tested
To check that the differences in the lengths of the PCR products were due to differences in the number of CA repeats at each microsatellite locus, PCR products of strain IP 2279.94 and an isolate (A.224) with distinct alleles were inserted into pGEM-T Easy vector (Promega, Madison, Wis.). After transformation of competent XL1-Blue cells (Stratagene), two white colonies were selected from each PCR product and the inserts were sequenced. The sequences obtained showed the different numbers of CA repeats which were responsible for the different sizes of alleles: (CA)9(GA)25 and (CA)9(GA)8 for microsatellite A, (CA)2C(CA)23 and (CA)2C(CA)10 for microsatellite B, (CA)8 and (CA)13 for microsatellite C, and (CA)21 and (CA)32 for microsatellite D for strain IP 2279.94 and isolate A.224, respectively. Subsequently, we decided to systematically run reference strain IP 2279.94 on each gel to check the correct automatic sizing with the GeneScan software.
To check the stability of the microsatellite repeats over time, DNAs from four subcultures of strain CBS 143.89 obtained in 1971, 1986, 1989, and 1995 (provided by J.-P. Latgé) were tested. This strain has been regularly subcultured since 1989 at the Pasteur Institute. Identical microsatellite profiles were observed between the four subcultures tested. The same results were obtained with strain IP 2279.94, which has been regularly subcultured in our laboratory since 1994.
Of the 100 isolates and 2 reference strains studied, 64 had a unique type and 38 were grouped into 16 types (12 types of two isolates each, 3 types of 3 isolates each, and 1 type of 5 isolates). Some isolates had a common type, although they were undoubtedly independent because they were collected at completely different places. For instance, one environmental isolate from a Parisian hospital had the same type as a clinical isolate from a German patient. For other isolates with a common type, their independence was questionable because they were collected at the same place, although at different times. For instance, the same type was recovered from four patients hospitalized in the same hospital but in different wards over a 22-month period.
The isolates were analyzed by correspondence analysis according to their clinical or environmental origin. Correspondence analysis produced a multidimensional representation of the data, which were reduced in plane to be more easily visualized. The first three planes were defined by the three major axes which explained most of the variance (relative inertia, 19%). The two ellipses which contained 90% of the individuals of each subgroup (patient and environmental isolates) widely overlapped and their centers were very close (Fig.3). The same results were obtained for the other projections. Therefore, the diversity of patient isolates was of the same order as that of the environmental isolates.
Correspondence analysis performed for all alleles showing the dispersion observed for the 102 A. fumigatus isolates. The projection in the plane defined by the two most informative axes is shown. Squares, environmental isolates; circles, clinical isolates; empty symbols, unique profiles; shaded symbols, profiles shared by two or more isolates. The two ellipses contain 90% of the individuals of each subgroup (dotted line for patient isolates; continuous line for environmental isolates). The centers of the ellipses are indicated, with a black circle for the patient isolate ellipse and a black square for the environmental isolate ellipse.
DISCUSSION
For the first time, the characterization of microsatellites ofA. fumigatus is reported and the use of microsatellites as a typing system was investigated. These microsatellites are polymorphic, and reproducible typing can be obtained by analysis with an automatic sequencer. Our results are consistent with those obtained with microsatellites of other species. In the human genome, the higher the number of uninterrupted CA repeats, the more polymorphic the microsatellite (27). Indeed, no polymorphism was detected with three microsatellites with CA repeat numbers of less than 8; 10 alleles were found with microsatellite C (CA repeat range, 7 to 16), and 23 alleles were found with microsatellite D (CA repeat range, 7 to 36). The reproducibility of this technique is high because (i) the PCR amplifications are performed under high-stringency conditions, (ii) the results do not depend on the DNA isolation technique, (iii) and the sizing is automatic, with an internal size standard systematically loaded along with the amplified products. Thus, the coding of alleles for further analysis is simple. Moreover, these microsatellite markers appear stable upon subculturing. Another advantage is that there is no need for monospore culture. Indeed, a single peak is expected for a given locus becauseA. fumigatus is haploid. If several peaks are observed, the isolate tested is actually a mixture of several isolates and can be rejected.
The combined discriminatory power of the four microsatellites described is 0.994 and outscores the discriminatory power obtained with other PCR-based markers (14, 16). Therefore, these markers appear promising for large epidemiological studies. From a technical point of view, the throughput of A. fumigatusisolates could be increased if several fluorescent dyes were used to label the primers in multiplex PCRs.
However, some technical problems can be encountered with this technique. Amplification of dinucleotide repeats can lead to PCR artifacts with several peaks separated by intervals of two nucleotides. These extra bands are mainly due to the slippage of Taq DNA polymerase during amplification (11). In our experience, this trouble was overcome by performing the amplification with dimethyl sulfoxide. Thus, the longest band was always the most intense and was the only one retained for analysis. Another PCR artifact is the addition of an extra A to the 3′ end of product strands, more specifically, when the 5′ end of the reverse primer is a G, which is the preferential template for the addition of extra nucleotides by theTaq polymerase (5). The final extension at 72°C for 30 min in our PCRs was aimed at ensuring the systematic addition of an extra A to make the PCR products homogeneous. Nevertheless, a reference strain must be systematically included as an internal control in all amplification runs and analyses. Because the results expected for the reference strain are known, it is possible to detect an artifact of the PCR or a sizing error.
Another question related to the microsatellite markers is whether the variations in the sizes of the PCR products are due to variations in the number of CA repeats. It has recently been reported (20) that the addition or deletion of nucleotides in the flanking regions of the dinucleotide repeat instead of in the microsatellite sequence itself could account for length variations. For the microsatellites tested in this work, the primers for amplification were chosen close to the CA repeat in order to generate PCR products shorter than 200 bp. Thus, the probability that mutations outside the microsatellite account for the length polymorphism is low. Therefore, in contrast to RAPD techniques, the nature of the differences is known. Currently, this information is not usable because it is impossible to assess the mutation rate at a given microsatellite in the absence of sexual reproduction for A. fumigatus. Nevertheless, it is reasonable to think that knowledge obtained with a model species with a sexuality, such as Emericella nidulans, could be transposable to A. fumigatus.
By RFLP, Debeaupuis et al. (7) have shown that no particular genotype was associated with virulence. We were also unable to cluster isolates depending on their clinical or environmental origin. This finding may indicate that every A. fumigatusisolate is potentially pathogenic and that the risk for the patient is inhalation of the conidia of any A. fumigatus isolate. This hypothesis is in contrast to the hypothesis that certain isolates are more pathogenic than others (17, 23).
The microsatellites described in this work represent a new class of highly polymorphic markers for A. fumigatus. Epidemiological studies by automated procedures with large numbers of isolates can be designed. Moreover, microsatellites offer perspectives for studying the genetic relatedness among A. fumigatusisolates and for genome mapping, as has already been done for the genomes of other species.
ACKNOWLEDGMENTS
We thank J.-P. Latgé from the Pasteur Institute for providing A. fumigatus isolates and R. Calderone from Georgetown University, Washington, D.C., for critical reading of the manuscript.
This work was supported by grant BQR 9014 R1 from the Paris XII University.
FOOTNOTES
- Received 16 March 1998.
- Returned for modification 22 April 1998.
- Accepted 1 June 1998.
- Copyright © 1998 American Society for Microbiology
