- Short Report
- Open Access
Next-generation genome sequencing can be used to rapidly characterise sequences flanking T-DNA insertions in random insertional mutants of Leptosphaeria maculans
Fungal Biology and Biotechnologyvolume 1, Article number: 10 (2014)
Banks of mutants with random insertions of T-DNA from Agrobacterium tumefaciens are often used in forward genetics approaches to identify phenotypes of interest. Upon identification of mutants of interest, the flanking sequences of the inserted T-DNA must be identified so that the mutated gene can be characterised. However, for many fungi, this task is not trivial as widely used PCR-based methods such as thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR) are not successful.
Next-generation Illumina sequencing was used to locate T-DNA insertion sites in four mutants of Leptosphaeria maculans, a fungal plant pathogen. Sequence reads of up to 150 bp and coverage ranging from 6 to 24 times, were sufficient for identification of insertion sites in all mutants. All T-DNA border sequences were truncated to different extents. Additionally, next-generation sequencing revealed chromosomal rearrangements associated with the insertion in one of the mutants.
Next-generation sequencing is a cost-effective and rapid method of identifying sites of T-DNA insertions, and associated genomic rearrangements in Leptosphaeria maculans and potentially in other fungal species.
Forward genetics approaches, where banks of mutants with random insertions of DNA are screened for phenotype changes, are often used to study fungi that are recalcitrant to reverse genetics approaches such as gene deletion by homologous recombination -. The inserted DNA, often T-DNA from Agrobacterium tumefaciens, usually integrates at a single site . Multiple isolates with random insertions are generated and screened for the phenotype of interest. The insertion site is then identified by thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR) or plasmid rescue ,. For TAIL-PCR, sequence flanking the T-DNA insertion is amplified using a set of nested T-DNA-specific primers paired with a short degenerate primer. TAIL-PCR has been used to obtain 74 to 100% of flanking sequences in fungi such as Fusarium oxysporum, Magnaporthe oryzae, Histoplasma capsulatum, as well as the lichen fungus Umbilicaria muehlenbergii ,-. However, the rate of success in identifying insertions in other fungi such as Trichoderma reesei, and the ascomycete plant pathogens Colletotrichum higginsianum and Leptosphaeria maculans, ranges from 30 to 66% ,,. For example; only 135 out of 400 (33.7%) flanking sequences were retrieved from the right borders in an analysis of T-DNA insertional mutants of the L. maculans . In plasmid rescue, genomic DNA of the T-DNA-containing mutant is cut with a restriction enzyme and then circularized with T4 DNA ligase, thus releasing a fragment of the T-DNA, a selectable marker and flanking genomic sequences ,. This method requires the use of a T-DNA plasmid containing a bacterial origin of replication and selectable marker and its successful use depends upon the proximity of appropriate restriction enzyme sites to the T-DNA insertion. The use of next-generation sequencing, for identifying the T-DNA insertion site, can overcome these limitations. In this paper we report the use of Illumina sequencing to identify T-DNA flanking sequences in four loss-of-pathogenicity mutants identified from a bank of L. maculans insertional mutants.
L. maculans isolate IBCN18 was transformed with plasmids pKHT or pKO21 ,. Plasmid pKO21 is a modified version of pPZPtk8.10, which has a hygromycin resistance gene regulated by the A. nidulans trpC promoter and terminator sequences inserted into the multiple cloning site. Four insertional mutants, UM176, UM314, UM316, UM317, were identified with a loss-of-pathogenicity phenotype on cotyledons of Brassica napus, as described previously . Genomic DNA (4 μg) was prepared from mycelia , digested with restriction enzymes, electrophoresed and then probed with the hygromycin resistance gene on a Southern analysis blot (data not shown). This confirmed that UM176, UM316 and UM317 each carried a single T-DNA whereas UM314 had two independent insertions of T-DNA.
Identification of the corresponding insertion sites in the four mutants was attempted using TAIL-PCR. Six different degenerate primers were used in combination with border-specific primers to amplify either the left border (LB) or right border (RB) sequences ,. No amplicons containing the corresponding T-DNA sequence were produced using any combination of TAIL-PCR primers. Plasmid rescue failed to retrieve flanking sequence from UM176, in which the T-DNA was derived from pKHT. Since both TAIL-PCR and plasmid rescue failed to identify insertion sites, Illumina next-generation sequencing was then undertaken.
Genomic DNA was extracted from mycelia, then treated with RNaseA (Invitrogen). Sequencing libraries (NexteraXT; Illumina) were prepared, pooled and run on a MiSeq desktop sequencer (Illumina) . Paired-end reads (150 bp) were generated, with genome coverage for the four mutants ranging from 6.6×-24× (Table 1). Sequencing reads were filtered to remove those of low quality and/or high percentage of Ns; adapter sequences were also removed. Separate BLASTN databases were generated for the sequence reads of each of the mutants using Geneious R6.0.5 (BLAST version Blast 2.2.28+). These BLAST databases were then queried using the sequences of the region between the LB and RB of the T-DNA insert. All reads that matched the T-DNA query, with expect scores of 1e−10, were retrieved and aligned to the query sequence. Flanking sequences that extended beyond the T-DNA query sequence were then used for sequence walking by querying the BLAST database again and thus retrieving more flanking sequences. This process of sequence walking was continued until enough flanking sequence was obtained to BLAST the GenBank nr database to identify neighboring regions in the L. maculans genome. Based on previous studies, a minimum alignment length of 15 bp unambiguously identifies read locations . Thus to identify an insertion site, a read needs to contain at least 15 bp of the T-DNA sequence, as well as at least 15 bp of L. maculans genome sequence. Therefore for each 150 bp read, there are 120 consecutive binding sites that meet these criteria. The genome of L. maculans is approximately 45 Mbp . If each read can bind to 120 consecutive loci in the genome then the probability of identifying a T-DNA flanking sequence is 2.6 × E−06 (120 bp/45,000,000 bp) or 1 in 345,000 reads. The predicted number of reads that should flank the T-DNA insertion site for each mutant, based on these criteria, is very similar to the actual number of reads identified (Table 1). A T-DNA insertion site was located in each of the four mutants (Table 1).
In three of the four mutants, the T-DNA was located between genes and in all mutants, the LB and/or RB sequences were partially or completely deleted, possibly explaining why TAIL-PCR failed to identify the T-DNA insertion site. In UM316, no sequences flanking the RB were identified probably due to its proximity to a region of repetitive DNA. In UM176 a large deletion (26 kb) appeared to be associated with a translocation between two Supercontigs (17 and 5), which in the wild-type are on different chromosomes. This rearrangement resulted in deletion of at least 12 genes.
To confirm the location of insertion sites within the genomes of mutant UM176, UM316 and UM317, primers were designed to amplify from inside the T-DNA insert out into the L. maculans genome. For all three mutants, PCR amplicons were successfully obtained and sequenced, confirming the location of the insertions (data not shown).
With the decreasing costs of next-generation sequencing, this method of identification of insertion sites in T-DNA mutants is a cost and time-effective solution for both fungi and for plants, where PCR-based methods have not been successful. Furthermore, next-generation sequencing allows identification of chromosomal rearrangements associated with the insertion but not specifically at the insertion site. Such rearrangements are common in both plants  and fungi such as the basidiomycete human pathogen, Cryptococcus neoformans and the ascomycete rice pathogen Magnaporthe oryzae ,. However, when using Illumina sequencing in situations where the T-DNA has inserted into a highly repetitive region of the genome, as seen in mutant UM316, it might not be possible to identify the precise location of the T-DNA junction and to quantify the amount of fungal DNA deleted during the T-DNA insertion event. Additionally, the presence of multiple T-DNA insertions or repeats within the vector sequences, such as in the case of UM314, can confound the alignment of the border sequences. Therefore all border sequences might not be identified by Illumina sequencing. However, if there is a single T-DNA insertion, in most cases the identification of just one border is sufficient to locate the insertion site and only a small deletion at the insertion site. An advantage of Illumina sequencing, as illustrated in this study, is that as little as 20–30 bp of flanking sequence was sufficient to identify the insertion site of the T-DNA. In contrast, in a previous study, TAIL-PCR of a set of L. maculans mutants resulted in 33 flanking sequences of insufficient length to identify the insertion site, and a further 49 sequences gave ambiguous BLAST matches .
Similarly, identification of T-DNA insertion sites in mutant lines of the plant Arabidopsis thaliana has proven problematic. For example, Ji and Braam (2010) developed restriction site extension PCR, which relies on multiple PCR steps, but this method resulted in recovery of only 21 of 37 flanking sequences . Lepage et al. (2013) used Roche 454 genomic sequencing and T-DNA specific primers to identify insertion sites in 55 out of 64 mutants . However, this approach relies on the binding of T-DNA specific primers, which can be problematic if the T-DNA borders are mutated during integration. Furthermore, Roche454 sequencing is relatively expensive. In the current study, Illumina reads of up to 150 bp containing as little as 20–30 bp of L. maculans sequence were sufficient for identification of insertion sites; furthermore this approach overcomes potential issues with loss of T-DNA primer binding sites, and chromosomal rearrangements can be identified. With organisms of larger genome size, the amount of sequence data must be proportionately larger to maintain the same likelihood of locating the flanking regions. With the compact size of fungal genomes and the popularity of T-DNA mutagenesis -, Illumina sequencing is an ideal method for characterization of insertion sites.
Blaise F, Remy E, Meyer M, Zhou L, Narcy JP, Roux J, Balesdent MH, Rouxel T: A critical assessment of Agrobacterium tumefaciens-mediated transformation as a tool for pathogenicity gene discovery in the phytopathogenic fungus Leptosphaeria maculans. Fungal Genet Biol 2007, 44: 123–138. 10.1016/j.fgb.2006.07.006
Choi J, Park J, Jeon J, Chi MH, Goh J, Yoo SY, Jung K, Kim H, Park SY, Rho HS, Kim S, Kim BR, Han SS, Kang S, Lee YH: Genome-wide analysis of T-DNA integration into the chromosomes of Magnaporthe oryzae. Mol Microbiol 2007, 66: 371–382. 10.1111/j.1365-2958.2007.05918.x
Leclerque A, Wan H, Abschutz A, Chen S, Mitina GV, Zimmermann G, Schairer HU: Agrobacterium-mediated insertional mutagenesis (AIM) of the entomopathogenic fungus Beauveria bassiana. Curr Genet 2004, 45: 111–119. 10.1007/s00294-003-0468-2
Mullins ED, Chen X, Romaine P, Raina R, Geiser DM, Kang S: Agrobacterium-Mediated Transformation of Fusarium oxysporum: An Efficient Tool for Insertional Mutagenesis and Gene Transfer. Phytopathology 2001, 91: 173–180. 10.1094/PHYTO.2001.91.2.173
Rho HS, Kang S, Lee YH: Agrobacterium tumefaciens-mediated transformation of the plant pathogenic fungus, Magnaporthe grisea. Mol Cells 2001, 12: 407–411.
Walton FJ, Idnurm A, Heitman J: Novel gene functions required for melanization of the human pathogen Cryptococcus neoformans. Mol Microbiol 2005, 57: 1381–1396. 10.1111/j.1365-2958.2005.04779.x
Michielse CB, van Wijk R, Reijnen L, Cornelissen BJC, Rep M: Insight into the molecular requirements for pathogenicity of Fusarium oxysporum f. sp lycopersici through large-scale insertional mutagenesis. Genome Biol 2009, 10: R4. 10.1186/gb-2009-10-1-r4
Kemski MM, Stevens B, Rappleye CA: Spectrum of T-DNA integrations for insertional mutagenesis of Histoplasma capsulatum. Fungal Biol 2013, 117: 41–51. 10.1016/j.funbio.2012.11.004
Park S-Y, Jeong M-H, Wang H-Y, Kim JA, Yu N-H, Kim S, Cheong YH, Kang S, Lee Y-H, Hur J-S: Agrobacterium tumefaciens-Mediated Transformation of the Lichen Fungus, Umbilicaria muehlenbergii. Plos One 2013, 8: e83896. 10.1371/journal.pone.0083896
Zhong YH, Wang XL, Wang TH, Jiang Q: Agrobacterium-mediated transformation (AMT) of Trichoderma reesei as an efficient tool for random insertional mutagenesis. Appl Microbiol Biotechnol 2007, 73: 1348–1354. 10.1007/s00253-006-0603-3
Huser A, Takahara H, Schmalenbach W, O'Connell R: Discovery of Pathogenicity Genes in the Crucifer Anthracnose Fungus Colletotrichum higginsianum, Using Random Insertional Mutagenesis. Mol Plant Microbe Interact 2009, 22: 143–156. 10.1094/MPMI-22-2-0143
Gardiner DM, Howlett BJ: Negative selection using thymidine kinase increases the efficiency of recovery of transformants with targeted genes in the filamentous fungus Leptosphaeria maculans . Curr Genet 2004, 45: 249–255. 10.1007/s00294-004-0488-6
Van de Wouw AP, Pettolino FA, Howlett BJ, Elliott CE: Mutations to LmIFRD affect cell wall integrity, development and pathogenicity of the ascomycete Leptosphaeria maculans. Fungal Genet Biol 2009, 46: 695–706. 10.1016/j.fgb.2009.06.004
Sexton AC, Howlett BJ: Characterisation of a cyanide hydratase gene in the phytopathogenic fungus Leptosphaeria maculans. Mol Gen Genet 2000, 263: 463–470. 10.1007/s004380051190
McElver J, Tzafrir I, Aux G, Rogers R, Ashby C, Smith K, Thomas C, Schetter A, Zhou Q, Cushman MA, Tossberg J, Nickle T, Levin JZ, Law M, Meinke D, Patton D: Insertional mutagenesis of genes required for seed development in Arabidopsis thaliana. Genetics 2001, 159: 1751–1763.
Zander M, Patel DA, Van de Wouw A, Lai K, Lorenc MT, Campbell E, Hayward A, Edwards D, Raman H, Batley J: Identifying genetic diversity of avirulence genes in Leptosphaeria maculans using whole genome sequencing. Funct Integr Genomics 2013, 13: 295–308. 10.1007/s10142-013-0324-5
Unneberg P, Wennborg A, Larsson M: Transcript identification by analysis of short sequence tags–influence of tag length, restriction site and transcript database. Nucleic Acids Res 2003, 31: 2217–2226. 10.1093/nar/gkg313
Rouxel T, Grandaubert J, Hane JK, Hoede C, van de Wouw AP, Couloux A, Dominguez V, Anthouard V, Bally P, Bourras S, Cozijnsen A.J, Ciuffetti LM, Degrave A, Dilmaghani A, Duret L, Fudal I, Goodwin SB, Gout L, Glaser N, Linglin J, Kema GH, Lapalu N, Lawrence CB, May K, Meyer M, Ollivier B, Poulain J, Schoch CL, Simon A, Spatafora JW, et al.: Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations. Nature Commun 2011, 2: n202. 10.1038/ncomms1189
Clark KA, Krysan PJ: Chromosomal translocations are a common phenomenon in Arabidopsis thaliana T-DNA insertion lines. Plant J 2010, 64: 990–1001. 10.1111/j.1365-313X.2010.04386.x
Idnurm A, Walton FJ, Floyd A, Reedy JL, Heitman J: Identification of ENA1 as a virulence gene of the human pathogenic fungus Cryptococcus neoformans through signature-tagged insertional mutagenesis. Eukaryot Cell 2009, 8: 315–326. 10.1128/EC.00375-08
Bourras S, Meyer M, Grandaubert J, Lapalu N, Fudal I, Linglin J, Ollivier B, Blaise F, Balesdent MH, Rouxel T: Incidence of genome structure, DNA asymmetry, and cell physiology on T-DNA integration in chromosomes of the phytopathogenic fungus Leptosphaeria maculans. G3 (Bethesda) 2012, 2: 891–904. 10.1534/g3.112.002048
Ji J, Braam J: Restriction site extension PCR: a novel method for high-throughput characterization of tagged DNA fragments and genome walking. PLoS One 2010, 5: e10577. 10.1371/journal.pone.0010577
Lepage E, Zampini E, Boyle B, Brisson N: Time- and cost-efficient identification of T-DNA insertion sites through targeted genomic sequencing. PLoS One 2013, 8: e70912. 10.1371/journal.pone.0070912
We thank the Australian Grains Research and Development Corporation and the Australian Research Council (Project DP0985953) for funding.
The authors declare that they have no competing interests.
KC analysed the mutants. RL designed the bioinformatics analyses. BH conceived the study. MZ and JB constructed sequencing libraries and performed the next generation sequencing. AV conceived the study, participated in its design and coordination, and drafted the manuscript. CE prepared the mutant banks, carried out southern analyses, and drafted the manuscript. All authors read and approved the final manuscript.