Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Efficient genome editing using tRNA promoter-driven CRISPR/Cas9 gRNA in Aspergillus niger

Abstract

As a powerful tool for fast and precise genome editing, the CRISPR/Cas9 system has been applied in filamentous fungi to improve the efficiency of genome alteration. However, the method of delivering guide RNA (gRNA) remains a bottleneck in performing CRISPR mutagenesis in Aspergillus species. Here we report a gRNA transcription driven by endogenous tRNA promoters which include a tRNA gene plus 100 base pairs of upstream sequence. Co-transformation of a cas9-expressing plasmid with a linear DNA coding for gRNA demonstrated that 36 of the 37 tRNA promoters tested were able to generate the intended mutation in A. niger. When gRNA and cas9 were expressed in a single extra-chromosomal plasmid, the efficiency of gene mutation was as high as 97%. Co-transformation with DNA template for homologous recombination, the CRISPR/Cas9 system resulted ~42% efficiency of gene replacement in a strain with a functioning non-homologous end joining machinery (kusA+), and an efficiency of >90% gene replacement in a kusA- background. Our results demonstrate that tRNA promoter-mediated gRNA expressions are reliable and efficient in genome editing in A. niger.

Introduction

The filamentous fungus Aspergillus niger is an important cell factory used in industry for the production of enzymes and organic acids. Owing to its genetic tractability, A. niger is widely used for research in fungal physiology, biochemistry and biotechnology. The availability of genome sequence of this organism [1,2] has facilitated numerous studies in gene function, gene regulation, primary and secondary metabolism [3,4]. The utility of A. niger as an industrial cell factory and as a model organism for research can be further improved by the development of a highly efficiency genome-editing system. https://www.ncbi.nlm.nih.gov/pubmed/17259976

Multiple nuclease-based gene targeting methods have been developed to improve the efficiency and precision of generating genetic modifications. These include Zinc-finger nuclease (ZFN) [5], Transcription Activator-Like Effector Nuclease (TALEN) [6], I-SceI meganuclease [7], and Clustered Regularly Interspaced Short Palindromic Repeat / associated nuclease 9 (CRISPR/Cas9) [8]. With the ability of nucleases to make DNA double strand break (DSB) in the host organism, nuclease-based gene targeting methods have been successfully used for gene disruption, knock-in mutation as well as improving heterologous protein production by integrating foreign genes to defined genomic loci [5,9,10]. The type II CRISPR/Cas9 system from Streptococcus pyogenes has revolutionized genome editing by providing a simple, flexible and efficient strategy to precisely edit the genomes of mammals [11], plant [12], fly [13], fungi [14] and bacteria [15]. Two components are essential in type II CRISPR/Cas9 system: a functional Cas9 nuclease and a chimeric guide RNA (gRNA) consisted of two regions, a CRISPR RNA (crRNA) harboring 20-nucleotide target-recognizing sequence at the 5’-end and a trans-activating crRNA (tracrRNA) for Cas9 binding. The crRNA is user-defined to match the target genomic locus. It guides the gRNA to form a RNA/DNA hybrid at the target genomic locus and recruits the Cas9 nuclease to generate DNA DSB [16]. Consequently, the CRISPR/Cas9 approach supports the editing of numerous genomic loci.

An efficient promoter that can facilitate gRNA transcription in vivo is a bottleneck of adoption of CRISPR/Cas9 system in Aspergilli. So far, in the reported CRISPR/Cas9 platform for A. niger, RNA polymerase II (Pol II) promoter and terminator were employed to drive the transcription of the gRNA, which required the inclusion of two self-splicing ribozyme sequences flanking the gRNA sequence to free the mature gRNA transcript [14]. The endogenous U6 promoter, an RNA polymerase III (Pol III) promoter, has been used for CRISPR/Cas9 genome editing in A. oryzae with an efficiency lower than 20% in generating gene mutations [17]. Another Pol III promoter, U3 promoter of A. fumigatus (AfU3), was successfully used to transcribe gRNA in A. nidulans. However, endogenous U6 and AfU6 promoters failed in generating functional gRNA in A. nidulans, which indicates that U6 promoters have limited function as gRNA promoters in fungi [18,19]. Therefore, it is important to find reliable promoters to drive gRNA transcription in vivo that can result in efficient genome editing of A. niger to expand its utility as a model and industrial organism.

In eukaryote cells, the transfer ribonucleic acid (tRNA) genes are transcribed by Pol III [20]. The transcription of fungal tRNA and other eukaryotes genes is constitutive and is independent of carbon source and cultivation conditions [21,22]. The tRNA promoter includes a region centered approximately 30 base pairs (bp) upstream of the tRNA gene start site and two intragenic sequences, Box A and Box B [20,23,24]. Their short sequences [25] and their self-splicing capacity [26,27] make tRNA promoters suitable for transcribing the short and non-coding gRNA sequence. The self-splicing ability of tRNA has been used to separate multiplexed gRNAs in plants [28] and recently in filamentous fungi [18]. To the best of our knowledge, using tRNA promoter to independently drive gRNA transcription has not been described in filamentous fungi. In this study, we describe the use of tRNA promoters and terminators to drive the gRNA transcription in vivo and show that this method is highly efficient in genome editing in A. niger.

Results

The use of tRNA promoters and terminators to drive gRNA transcription in A. niger

We manually curated the genome of the A. niger NRRL3 (http://gbrowse.fungalgenomics.ca/cgi-bin/gb2/gbrowse/Aspni_nrrl3_public) and identified 284 tRNA genes (chromosomal locations and sequences were listed in S1 Table). Six of the 284 tRNA genes are pseudogenes and three encode selenocystein tRNA. There are 8–24 tRNA genes for 17 of the 20 main amino acids. The three amino acids with fewer than eight corresponding tRNA genes are: tryptophan with one, asparagine with three, and cysteine with four tRNA genes (S1 Table). To determine the ability of tRNA promoters to drive the functional transcription of gRNA, we chose two random tRNA genes from every amino acid except that of asparagine and tryptophan. Three valine tRNAs were chosen due to the lack of intron in one and the difference in intron lengths in the two others. In total, 37 tRNA genes and their flanking sequences to use as promoters and terminators were selected (Table 1). Each gRNA cassette contains the following elements: 1) tRNA promoter which comprises the entire tRNA gene including introns, plus 100 nt upstream sequence; 2) the gRNA sequence that comprises a 20 nt crRNA guiding sequence that targets the albA gene (gene ID: NRRL3_00462) and a 80 nt tracrRNA sequence for Cas9 protein binding [29]; and 3) the 7–42 nt of tRNA terminator that spans from the end of tRNA gene to the poly-T stretch. The gene albA encodes a polyketide synthase involved in both dihydroxynaphthalene melanin and naphtho-gamma-pyrone synthesis that are required for spore pigmentation [30]. The disruption of albA prevents biosynthesis of black pigment, and results in white mutant conidia [14,30]. This phenotypic change provided a visual scoring of mutants directly on the transformation plate. For a quick, initial test of the ability of tRNA promoters and terminators to regulate gRNA transcription, we co-transformed the 37 linear gRNA cassettes independently into A. niger with cas9-bearing plasmid ANEp8-Cas9 (Fig 1A). The ANEp8-Cas9 plasmid was derived from the plasmid ANEp8 [31] that contains AMA1 sequence allowing extra-chromosomal replication [32] and the pyrG selection marker.

thumbnail
Fig 1. The schematic representation of tRNA-driven CRISPR/Cas9 system adapted for A. niger.

(A) The ANEp8-Cas9 plasmid used to express cas9 gene. (B) The construct of gRNA cassette and the ANEp8-Cas9-gRNA plasmid bearing both gRNA cassette and cas9 gene. The tRNA promoter is composed of a tRNA gene and 100 bp upstream sequences. Arrow indicates the start site of transcription. Abbreviations: bla, beta-lactamase gene conferring ampicillin resistance; pkiA Prom, promoter of the pyruvate kinase gene; glaA Term, terminator of the glucoamylase gene; NLS: nuclear localization signal.

https://doi.org/10.1371/journal.pone.0202868.g001

thumbnail
Table 1. Functional testing of 37 tRNA promoters and terminators in CRISPR/cas9 genome editing.

https://doi.org/10.1371/journal.pone.0202868.t001

In the control plate where the cells were transformed with ANEp8-Cas9 alone, only black colonies appeared (S1A Fig). Except in the case of tRNAGly5, all co-transformation yielded a mixture of black and white mutant colonies. These results suggest that almost all endogenous tRNA promoters and terminators are able to transcribe functional gRNA in the CRISPR/Cas9 system in A. niger. However, the co-transformation method resulted low inactivation rate of albA gene, ranging between 2 to 13% (Table 1). The relatively low rates of mutation using this co-transformation method is expected because the frequency of maintaining extra-chromosomal plasmid harboring the selectable marker is much higher than that of maintaining the integrated exogenous DNA [32]. In other words, only an uncontrolled subpopulation of transformants in these experiments is predicted to contain the linear gRNA cassette. Consequently, this initial experiment only reveals that most endogenous tRNA promoters and terminators are able to drive gRNA transcription and support CRISPR/Cas9 mutagenesis in A. niger, it does not provide quantitative information on the efficiency of different tRNA promoters in generating functional gRNA.

tRNA promoters are highly efficient in driving gRNA transcription

To compare the efficiency of different tRNA promoters in CRISPR/Cas9 system, the gRNA cassette was cloned into ANEp8-Cas9 to generate plasmid ANEp8-Cas9-gRNA (Fig 1B). We used extra-chromosomal plasmid for three reasons: 1) it supports high frequency of transformation [33]; 2) the transcription level of cas9 and gRNA are expected to be similar among transformants as they are not influenced by the site of integration; and 3) the plasmid can be removed easily by counter selection, leaving no exogenous DNA.

From the 36 positive tRNA-driven gRNAs (targeting albA) tested in the above co-transformation experiment, we chose eleven gRNA cassettes to construct plasmids containing both cas9 and gRNA. The eleven cassettes include tRNAs carrying hydrophilic (tRNAArg21, tRNACys1 and tRNAGln2), polar aliphatic (tRNAAsp5), aliphatic hydrophobic (tRNAAla5, tRNALeu6, tRNAPro1 tRNAVal1; 4; 15) and aromatic hydrophobic (tRNATyr3) amino acids. Among them three valine tRNA promoters were included because their tRNA genes contain either no intron or different intron lengths (Table 1).

The eleven ANEp8-Cas9-gRNA plasmids were transformed individually into A. niger strain N593 (Table 2). The results showed that ten out of eleven cassettes gave rise to transformants of which 82–96% were mutants (Table 3 and S1A Fig). The tRNAGln2 cassette resulted in significantly lower mutation rate of 15±6%. Additionally, the three tRNAVal promoters displayed a similarly high level of mutational efficiency, suggesting that tRNA introns do not influence its promoter capability in gRNA transcription.

thumbnail
Table 3. Gene disruption rate mediated by using eleven different tRNA promoters in ANEp8-Cas9-gRNA plasmid.

https://doi.org/10.1371/journal.pone.0202868.t003

To test the reliability of tRNA promoters in driving functional gRNA transcription in the CRISPR/Cas9 system, we targeted another gene olvA (ID: NRRL3_01039) that is also involved in the conidial pigmentation pathway. The loss of OlvA function results in olive-colored conidial mutants [30] (S1 Fig). We used five different tRNA promoters to drive gRNA to target olvA, four of which were randomly selected from the ten promoters that yielded >80% mutational frequency. We included the tRNAGln2 promoter to assess whether the low frequency of obtaining mutants is promoter dependent. The mutation rate of olvA was very similar to that observed towards albA (Table 3). A mutation rate of 13±2% was achieved by using tRNAGln2, while tRNAArg21, tRNACys1, tRNALeu6, and tRNAPro1 were able to induce mutation rates of between 88 to 97%. These results indicate that 1) the efficiency of gene editing by the CRISPR/Cas9 system is promoter-dependent; and 2) most endogenous tRNAs are very efficient and reliable in driving gRNA transcription in the CRISPR/Cas9 system.

Both deletions and insertions are generated by tRNA promoter-driven CRISPR/Cas9 mutagenesis

To determine the nature of CRISPR/Cas9-mediated mutations, we PCR-amplified the target genomic regions of the phenotypic mutants and sequenced the PCR products. Thirty transformants were analyzed, representing 15 albA and 15 olvA mutants mediated by five different tRNA promoters. These fragments covered the regions of albA and olvA that were expected to be cleaved by Cas9. Table 4 shows that deletions were the predominant mutation with 22 out of 30 of Cas9-induced mutations appeared as short deletions in the range of 1 to 27 bp at the Cas9 cleavage site. Additionally, five mutants contained long deletions of 612–1051 bp. Insertions were observed in three mutants. Sequence analysis indicated that the inserted sequences were derived from ANEp8-Cas9-gRNA plasmid. These observations suggest that double strand break mediated by Cas9 induces error-prone repair that resulted in mutation of the targeted gene in A. niger. The results also suggest that the extra-chromosomal plasmid can function as pseudo-donor template in the insertion of sequences.

thumbnail
Table 4. Mutations generated by tRNA promoter-driven CRISPR/Cas9 mutagenesis.

https://doi.org/10.1371/journal.pone.0202868.t004

CRISPR/Cas9 mediates gene replacement in A. niger kusA+ and kusA- backgrounds with different efficiencies

To evaluate how CRISPR/Cas9-mediated DSB can improve gene targeted replacement in A. niger, we constructed a ANEp8-Cas9-gRNAglaA plasmid targeting the glaA locus, which is a highly expressed gene encoding glucoamylase, and a linear gene replacement cassette as donor template (S2 Fig). The linear donor template comprised the coding region of adaR (gene ID: NRRL3_09545) flanked by 600 bp of glaA promoter and terminator sequences. The adaR gene encodes a transcription regulator involved in the biosynthesis of TAN-1612/BMS-192548, a pigmented polyketide secondary metabolite [36]. The integration and overexpression of adaR yields orange colonies because of the overproduction of TAN-1612/BMS-192548 [36]. Orange colonies can be directly counted on transformation plates. As a control, we transformed donor DNA with the plasmid ANEp8-Cas9 without the gRNA cassette. In addition, we created the kusA- mutant strain according to Meyer et al. [35] and described in S3 and S4 Figs to evaluate the role of non-homologous end joining (NHEJ) DNA repair in mediating CRISPR/Cas9 mutagenesis. Mutants deficient in KusA are defective in NHEJ and cannot repair DNA DSBs [37,38]. Previous studies showed improved homologous recombination in NHEJ-deficient A. niger strains [3840]. However, little is known about the efficiency of homologous recombination when DSB is mediated by CRISPR/Cas9.

In order to easily distinguish gene integration mutants from abortive or young transformants, nine orange colonies from each transformation experiment were chosen for PCR screening using primers outside the glaA locus (S2 Fig). Table 5 shows that in the kusA+ background, without the gRNA no gene replacement occurred. In the presence of the gRNA in the kusA+ background, about 75% of the transformants displayed orange colonies and 5 of 9 orange colonies tested were the results of homologous recombination, suggesting a frequency of gene replacement of about 42%. In the kusA- background, the absence of gRNA yielded 21 orange colonies out of 75 transformants with 6 of 9 orange colonies tested positive for homologous recombination, suggesting a frequency of gene replacement of about 20%. This finding supports previous observation that the turning off of NHEJ repair that improves homologous recombination in A. niger [39,40]. In the kusA- background and in the presence of gRNA, orange colonies comprised about 93% of the transformants with all 9 orange colonies tested positive for homologous recombination. This would imply a gene replacement rate of ~93%. In addition, in the kusA- strain, no transformant was found in the experiment in which only the ANEp8-Cas9-gRNAglaA plasmid was transformed because fungal cells could not repair the DSB in the absence of a repair donor template (Table 5). Based on this reasoning, we would expect 100%, instead of 93%, gene replacement frequency in the situation where we have functional gRNA and a gene replacement DNA template in the kusA- background. The few non-orange transformants could be the result of less than 100% efficiency of the gRNA and/or some of the young transformants scored did not have sufficient time to manifest the mutant phenotype. Since kusA- strains cannot repair DSB, constructing mutants in the kusA- genetic background using CRISPR/Cas9 and rescue templates has the advantage that off-target mutations would be minimized as ectopic mutations would result in unrepaired chromosomal breaks and cell lethality.

thumbnail
Table 5. CRISPR/Cas9-mediated targeting homologous integration in A. niger.

https://doi.org/10.1371/journal.pone.0202868.t005

Discussion

Besides an active Cas9 protein, a functional gRNA is an essential component in CRISPR/Cas9 system to engineer the host genome. For delivering gRNA in vivo, the formation of functional gRNA requires an efficient promoter to drive gRNA transcription and a splicing mechanism to form mature gRNA transcript. In this study, we describe a strategy using Pol III tRNA promoters to drive gRNA transcription. This study shows that the self-splicing ability of tRNAs is highly efficient in generating functional gRNAs.

Since 2015, the CRISPR/Cas9 system in genome editing has been established in several genera of filamentous fungi. As summarized in Table 6, the promoter selection has shifted from Pol II promoter (e.g. pgdA and trpC promoters) to Pol III promoters for the in vivo transcription of gRNA. This is because Pol III promoter is more suitable for in vivo transcription of short non-coding transcripts such as gRNA [20,41]. In addition, the short length of Pol III promoters and terminators makes it easier to construct gRNA cassette. Among Pol III promoters used in filamentous fungi CRISPR/Cas9 platforms, the species-specific U6 promoter is the most widely used [17,4247], and the SNR52 promoter from Saccharomyces cerevisiae has been tested in Neurospora crassa [48] and A. fumigatus [49] (Table 6). In this study, the tRNA promoters contain the entire sequence of tRNA gene plus 100 nt upstream sequence with a total sequence varying from 171 to 272 nt (S2 Table), which is similar as the lengths of U6 (200 to 250 nt) and S. cerevisiae SNR52 (269 nt) promoters [29,50]. The main challenge of current Pol III gRNA promoters such as U3 and U6 is that the same type of promoter results in highly variable gene editing efficiency amongst different species. Unlike U6 promoter that is an upstream element outside of the transcription start site [50], the tRNA gene is transcribed together with gRNA sequence [20,51]. In the cell the tRNA transcript is posttranscriptionally cleaved at the 3‘-end by RNase Z [27,28]. Accordingly, tRNA promoter regulated system is able to release mature functional gRNA transcript as demonstrated in this study and elsewhere [52].

thumbnail
Table 6. Summary of CRISPR/Cas9-mediated gene mutation and replacement efficiencies in filamentous fungi.

https://doi.org/10.1371/journal.pone.0202868.t006

By introducing upstream sequences, the tRNA promoters in our constructs include all external and internal binding sites of Pol III transcription factors needed for transcription to proceed normally. Except tRNAGln2, the tRNA promoters tested in this study lead to high rate of gene mutation (82–97%). In endogenous U6 promoter-driven platforms, the mutation rate is highly variable among different fungal species and ranges from 10 to 100% (Table 6). In a recent study, the endogenous U6 promoter and the heterologous AfU6 promoter failed to transcribe gRNA in CRISPR/Cas9 system in A. nidulans [18]. In addition, gRNAs transcribed by S. cerevisiae SNR52 promoter resulted in mutation rates of between 25 to 53% in A. fumigatus [49] (Table 6). These studies show that the type of promoters used to transcribe gRNA has a high impact in CRISPR/Cas9 activity. Overall results show that expressing Cas9 and the gRNA in single plasmid improve gene deletion and gene replacement. In addition no difference on growth was observed when transformed A. niger with ANEp8-Cas9 or ANEp8-Cas9-gRNA in the present study, except the color phenotype (S1A and S1B Fig). This study provides a broad list of highly efficient tRNA promoters for gRNA transcription.

In A. niger and other eukaryotes, DSB created in the genome is repaired by NHEJ or by homologous recombination [55]. The NHEJ is the most dominant pathway in A. niger and the frequency of homologous recombination increases with the deletion of the kusA gene, which is involved in the NHEJ repair mechanism [56]. In this study, we tested CRISPR/Cas9 mediated gene replacement by homologous recombination in both kusA- and kusA+ strains. By targeting the glaA locus, we demonstrated that the CRISPR/Cas9 method can increase the frequency of gene replacement in A. niger (Table 5). In gene replacement mediated by CRISPR/Cas9 activity, our results are comparable to the results obtained in filamentous fungi such as A. nidulans, Trichoderma reesei and Myceliophthora thermophila, and show higher frequency comparing to findings in A. aculeatus, A. fumigatus, Penicillium chrysogenum and A. carbonarius (Table 6). In summary, the tRNA promoter-driven gRNA transcription system is very efficient to generate DNA double-strand break for gene mutation mediated by NHEJ and for gene replacement by homologous recombination in kusA+ and kusA- strains. This powerful tool will facilitate further improvement of the industrial fungus A. niger.

Material and methods

CRISPR design

tRNA gene prediction and selection.

The tRNA gene and promoter sequences for this study were retrieved from the genome sequence of A. niger NRRL3 (data not published, available at http://genome.fungalgenomics.ca). We used four tRNA gene electronic tools tRNAscan-SE version 1.3 [57], SPLITSX [58], ARAGORN [59] and tRNAfinder [60] to predict all potential cytoplasmic tRNA genes in the nuclear genome of A. niger NRRL3. A total of 451 tRNA gene models were predicted by the four tools. We manually reviewed the predicted models to ensure that they conform to the length and structure of tRNAs and that the terminator contains a stretch of consecutive thymidines (poly-T) [61,62]. Where multiple non-identical models are detected in the same genomic region, we chose the consensus model predicted by more than one program. Since no fungal tRNA gene has been found to have intron exceeding 100 nucleotides [63], we removed 14 tRNA models predicted by ARAGORN to have introns >100 nucleotides. The final set contains 284 tRNA gene models (S1 Table).

Identification of guide sequences.

We used Geneious R9.1 software [64] to search 20 bp guide sequences against the A. niger NRRL3 genome. S3 Table shows the guide sequences used in this study, and none of them possesses off-target potential.

Construction of gRNA cassette.

The gRNA cassettes were assembled by fusion PCR of two DNA fragments, and the primers used were listed in S4 Table. The first fragment was the tRNA promoter that was amplified from A. niger genomic DNA and the 20 bp guide sequence was added at the 3’-end through PCR using as overlapping sequence. The second fragment was amplified from the synthesized 80 bp tracRNA template [29] with the 20 bp guide and tRNA terminator sequences (S3 Table) added at 5‘and 3‘ends respectively via PCR. In the fusion PCR, 0.5μL of each unpurified PCR fragments were used as template in a 100 μL reaction and performed with Phusion DNA polymerase (Thermo Scientific).

Plasmid construction

Escherichia coli DH5ɑ was used as the host strain for maintenance and production of plasmid DNA.

Construction of plasmid ANEp8-Cas9.

The coding sequence of S. pyogenes cas9 (Uniprot Q99ZW2) was codon optimized for translation in A. niger. A nucleotide sequence encoding a nuclear localization signal (NLS) peptide SRADPKKKRKV was fused to the carboxy terminus of cas9. The DNA sequence of Cas9-NLS was synthesized by IDT Inc. (IA, USA), and amplified by primers Cas9-F and Cas9-R (S4 and S5 Tables) to add PacI and FseI restriction sites at the 5’ and 3’ ends, respectively. Restriction enzyme digests by PacI and FseI were used to clone cas9-NLS downstream of the pkiA promoter in the expression vector ANEp8 [31]. The resulting cas9-containing plasmid was designated as ANEp8-Cas9 (S5 Table).

Construction of plasmid ANEp8-Cas9-gRNA.

A 38-bp ligation independent cloning (LIC) site centered with SwaI restriction site (S5 Fig) was introduced into ANEp8-Cas9 plasmid via PCR amplification. The resulting plasmid was used as host vector to harbor the gRNA cassette by using the LIC method [65]. The gRNA cassette used for plasmid construction was amplified with a pair of end primers to link with LIC sequence sites at the both sides as illustrated in S5 Fig, to generate complementary single-strand overhangs between ANEp8-Cas9 vector and gRNA cassette insert. The linearized ANEp8-Cas9 and the gRNA cassette DNA (ending with LIC tails) were treated by T4 DNA polymerase in the presence of dGTP and dCTP, respectively. The 20 μL reaction mixture contained 0.2 pmol of DNA, 0.8 μL of 100 mM dithiothreitol, 2 μL of 25 mM dGTP or dCTP, and 3 U of T4 DNA polymerase in NEB buffer 2.1. The reaction was carried out at 22°C for 30 minutes followed by enzyme inactivation by heating at 75°C for 20 minutes. The insert and vector were mixed in a 3:1 molar ratio. To achieve annealing, the mixture was first heated at 60°C for 5 minutes and then gradually decreased to 4°C (reduce 0.1°C per second). The annealed products were transformed into E.coli DH5α competent cells to generate plasmid ANEp8-Cas9-gRNA. All plasmids are available upon request.

Construction of adaR cassette for gene replacement at the glaA locus

The linear adaR (gene ID: NRRL3_09545) cassette with glaA flanking regions was constructed as shown in S2 Fig and S5 Table. Using genomic DNA of the A. niger strain NRRL2270 as template, 600 bp of 5’and 3’ regions flanking the coding region of the glaA gene as well as the adaR gene were amplified by primers with complementary ends (S4 Table). Based on their terminal overlaps, the three fragments were joined through fusion PCR amplification, resulting in an adaR cassette that was used for gene replacement at the glaA locus. A total of 5 μg of purified linear adaR cassettes were used in transformation. The PCR screening for targeted integration was carried out by using primers Fw_glaA-ext and Rv_glaA-ext (S2 Fig and S4 Table).

CRISPR/Cas9 genome editing in A. niger

Strain, media and growth conditions.

Table 2 lists the A. niger strains used in this study. The uridine auxotrophic A. niger strain N593 was used for gene editing and to test tRNA promoter efficiency by CRISPR/Cas9. The strains NRRL2270ΔpyrG and NRRL2270ΔpyrGΔkusA were used for testing gene replacement by CRISPR/Cas9. All A. niger strains were grown on minimal medium [66] or on complete medium (minimal medium supplemented with 5 g/L of yeast extract and 1 g/L of casamino acids). As required, a final concentration of 1.3 mg/mL of 5-fluoroorotic acid (5-FOA) and/or 10 mM of uridine were added. The liquid media were supplemented with 15% D-maltose. The transformants were cultured at 30°C for 4 days.

Transformation, mutant purification and screening.

Transformation and genomic DNA isolation of A. niger strains were performed as described previously [3]. All transformations were carried out using 1 μg of plasmid DNA together with 5 μg of linear cassette DNA when needed. Triplicate transformations were performed for each tRNA promoter tested, and duplicate transformations were carried out for the gene replacement of glaA. Mutant clones were visually identified according to conidia color. The CRISPR efficiency was calculated by dividing the number of mutant colonies by the total number of transformants. Three mutant colonies were randomly selected from the transformation plate and streak-purified twice on new minimal medium plates without uridine. Mutation patterns were analyzed by PCR amplification of albA (1.8 kb) and olvA (2.2 kb) sequences surrounding the CRISPR/Cas9 cut site using primers listed in S4 Table. Amplified DNA bands were visualized on agarose gel. Detailed information on the mutations was determined by DNA sequencing of the PCR products.

Supporting information

S1 Fig. An example of transformation plates for gene mutations and growth of purified transformants.

A) From left to right are the transformation plates of A. niger N593 cells transformed with ANEp8-Cas9 plasmid (only expressed cas9), and ANEp8-Cas9-gRNA plasmid bearing tRNAPro1-driven gRNA to disrupt albA (98% efficiency) and olvA (95% efficiency) respectively. B) Growth phenotype of purified colonies of A. niger transformed with ANEp8-Cas9 and ANEp8-Cas9-gRNAolvA.

https://doi.org/10.1371/journal.pone.0202868.s001

(PDF)

S2 Fig. Gene replacement of glaA gene.

(A) Schematic illustration of targeted replacement by homologous recombination. (B) PCR screening for CRISPR/Cas9-mediated homologous integration among the nine-selected orange transformants from kusA+ (left) and kusA- (right) strains. The PCR was performed using their genomic DNA as template with primers outside of glaA locus. Replacement of glaA by adaR gave a 4.2 kb PCR fragment. Random integration of adaR cassette would leave the glaA locus intact, resulting in a 3.8 kb PCR fragment.

https://doi.org/10.1371/journal.pone.0202868.s002

(PDF)

S3 Fig. Construction of kusA deletion cassette and strategy to generate NRRL2270ΔkusAΔpyrG strain.

(A) Fusions PCR to construct the kusA deletion cassette by including 500 bps repeat of the 3’ UTR in 5’ of the pyrG selection marker to facilitate the loop out of the selection marker. (B) Loop out strategy of the pyrG selection marker to generate ΔkusAΔpyrG strain for subsequent transformation with the same selection marker: (1) Homologous integration of the kusA deletion cassette to remove the kusA gene in A. niger; (2) Generation of ΔkusAΔpyrG genotype in presence of 5-FOA plus uridine to facilitate pyrG loop out (3).

https://doi.org/10.1371/journal.pone.0202868.s003

(PDF)

S4 Fig. PCR screening of transformants of ΔkusApyrG+ and ΔkusAΔpyrG genotypes respectively and the growth comparison with parental strain.

(A) Screening of transformants for ΔkusApyrG+ genotype with primers binding externally to the kusA locus and primers binding within the pyrG selection marker. Two PCR bands with size of 3 kb and 2.5 kb respectively are expected by using primer pairs P1_ext/P6_int and P5_int/P8_ext. Transformants 3, 6, 9 and 13 have the deletion cassette integrated in the kusA locus and have ΔkusApyrG+ genotype, when transformants 4, 5, 7, 8, 10, 11, 12 and 14 have partial integration of the deletion cassette. For transformants lacking pyrG, the primers are not able to bind and amplify the gDNA of the parental strain (PS) NRRL2270pyrG-. (B) Screening for transformants with pyrG loop out. Spores of transformant 3 (NRRL2270ΔkusApyrG+) have been plated on minimal medium containing 5-FOA and uridine for 5 days. Colonies growing on 5-FOA are selected for screening pyrG selection marker loop out phenotype by PCR. A correct loop out of the pyrG marker is characterized by a PCR giving a band size of 2.2 kb, while a strain failed in pyrG loop out will display a band of 4 kb. Transformant 3’ with the correct loop out of the pyrG selection marker is NRRL2270ΔkusAΔpyrG and has been selected for further work. (C) Growth comparison of NRRL2270ΔpyrG and NRRL2270ΔkusAΔpyrG. One million spores of strains NRRL2270ΔkusAΔpyrG and NRRL2270ΔpyrG were plated on the uridine-contained minimal or complete medium, and incubated at 30°C for 3 days.

https://doi.org/10.1371/journal.pone.0202868.s004

(PDF)

S5 Fig. The construction of ANEp8-Cas9-gRNA plasmid using LIC method.

In ANEp8-Cas9 linearized vector, the single-strand 5' overhangs for LIC cloning were achieved by 3' →5' exonuclease activity of T4 DNA polymerase in the presence of dGTP. For gRNA cassette, the LIC tails were added to the ends via PCR, and the reverse complementing overhangs were generated by the same T4 DNA polymerase treatment process with dCTP. Through base pairing, plasmid ANEp8-Cas9-gRNA was assembled by annealing the complementary sequences between gRNA insert and plasmid vector. The complementary LIC sequences in the gRNA insert and ANEp8-Cas9 vector are shown in the figure.

https://doi.org/10.1371/journal.pone.0202868.s005

(PDF)

S1 Table. The 284 predicted tRNA genes from A. niger NRRL3 genome.

https://doi.org/10.1371/journal.pone.0202868.s006

(XLSX)

S2 Table. DNA sequence of the 37 tRNA promoters and terminators used in the functional test.

https://doi.org/10.1371/journal.pone.0202868.s007

(XLSX)

S3 Table. The guide sequences used in this study.

https://doi.org/10.1371/journal.pone.0202868.s008

(DOCX)

S5 Table. DNA sequences of plasmids ANEp8-Cas9 and ANEp8-Cas9-gRNAglaA, and linear construct of adaR.

https://doi.org/10.1371/journal.pone.0202868.s010

(PDF)

Acknowledgments

We thank Sylwester Palys for reviewing the manuscript.

References

  1. 1. Pel HJ, de Winde JH, Archer DB, Dyer PS, Hofmann G, Schaap PJ, et al. Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88. Nat Biotechnol. 2007;25: 221–231. pmid:17259976
  2. 2. Meyer V. Genetic engineering of filamentous fungi—Progress, obstacles and future trends. Biotechnol Adv. 2008;26: 177–185. pmid:18201856
  3. 3. Meyer V, Wanka F, van Gent J, Arentshorst M, van den Hondel CAMJJ, Ram AFJ. Fungal gene expression on demand: an inducible, tunable, and metabolism-independent expression system for Aspergillus niger. Appl Environ Microbiol. 2011;77: 2975–2983. pmid:21378046
  4. 4. Andersen MR. Elucidation of primary metabolic pathways in Aspergillus species: orphaned research in characterizing orphan genes. Brief Funct Genomics. 2014;13: 451–455. pmid:25114096
  5. 5. Bibikova M, Beumer K, Trautman JK, Carroll D. Enhancing gene targeting with designed zinc finger nucleases. Science. 2003;300: 764–764. pmid:12730594
  6. 6. Boch J, Scholze H, Schornack S, Landgraf A, Hahn S, Kay S, et al. Breaking the code of DNA binding specificity of TAL-type III effectors. Science. 2009;326: 1509–1512. pmid:19933107
  7. 7. Bellaiche Y, Mogila V, Perrimon N. I-SceI endonuclease, a new tool for studying DNA double-strand break repair mechanisms in Drosophila. Genetics. 1999;152: 1037–1044. pmid:10388822
  8. 8. Gasiunas G, Barrangou R, Horvath P, Siksnys V. Cas9–crRNA ribonucleoprotein complex mediates specific DNA cleavage for adaptive immunity in bacteria. Proc Natl Acad Sci. 2012;109: E2579–E2586. pmid:22949671
  9. 9. Ouedraogo JP, Arentshorst M, Nikolaev I, Barends S, Ram AFJ. I-SceI enzyme mediated integration (SEMI) for fast and efficient gene targeting in Trichoderma reesei. J Biotechnol. 2016;222: 25–28. pmid:26860210
  10. 10. Christian M, Cermak T, Doyle EL, Schmidt C, Zhang F, Hummel A, et al. Targeting DNA double-strand breaks with TAL effector nucleases. Genetics. 2010;186: 757–761. pmid:20660643
  11. 11. Chu VT, Weber T, Wefers B, Wurst W, Sander S, Rajewsky K, et al. Increasing the efficiency of homology-directed repair for CRISPR-Cas9-induced precise gene editing in mammalian cells. Nat Biotechnol. 2015;33: 543–548. pmid:25803306
  12. 12. Woo JW, Kim J, Kwon SI, Corvalán C, Cho SW, Kim H, et al. DNA-free genome editing in plants with preassembled CRISPR-Cas9 ribonucleoproteins. Nat Biotechnol. 2015;33: 1162–1164. pmid:26479191
  13. 13. Port F, Chen H-M, Lee T, Bullock SL. Optimized CRISPR/Cas tools for efficient germline and somatic genome engineering in Drosophila. Proc Natl Acad Sci. 2014;111: E2967–E2976. pmid:25002478
  14. 14. Nødvig CS, Nielsen JB, Kogle ME, Mortensen UH. A CRISPR-Cas9 system for genetic engineering of filamentous fungi. Yu J-H, editor. PLOS ONE. 2015;10: e0133085. pmid:26177455
  15. 15. Jiang Y, Chen B, Duan C, Sun B, Yang J, Yang S. Multigene editing in the Escherichia coli genome via the CRISPR-Cas9 system. Appl Environ Microbiol. 2015;81: 2506–2514. pmid:25636838
  16. 16. Sander JD, Joung JK. CRISPR-Cas systems for genome editing, regulation and targeting. Nat Biotechnol. 2014;32: 347–355. pmid:24584096
  17. 17. Katayama T, Tanaka Y, Okabe T, Nakamura H, Fujii W, Kitamoto K, et al. Development of a genome editing technique using the CRISPR/Cas9 system in the industrial filamentous fungus Aspergillus oryzae. Biotechnol Lett. 2016;38: 637–642. pmid:26687199
  18. 18. Nødvig CS, Hoof JB, Kogle ME, Jarczynska ZD, Lehmbeck J, Klitgaard DK, et al. Efficient oligo nucleotide mediated CRISPR-Cas9 gene editing in Aspergilli. Fungal Genet Biol. 2018;
  19. 19. Zhang T, Gao Y, Wang R, Zhao Y. Production of guide RNAs in vitro and in vivo for CRISPR using ribozymes and RNA polymerase II promoters. BIO-Protoc. 2017;7.
  20. 20. White RJ. Transcription by RNA polymerase III: more complex than we thought. Nat Rev Genet. 2011;12: 459–463. pmid:21540878
  21. 21. Hopper AK, Pai DA, Engelke DR. Cellular dynamics of tRNAs and their genes. FEBS Lett. 2010;584: 310. pmid:19931532
  22. 22. Cieśla M, Towpik J, Graczyk D, Oficjalska-Pham D, Harismendy O, Suleau A, et al. Maf1 is involved in coupling carbon metabolism to RNA polymerase III transcription. Mol Cell Biol. 2007;27: 7693–7702. pmid:17785443
  23. 23. Male G, von Appen A, Glatt S, Taylor NMI, Cristovao M, Groetsch H, et al. Architecture of TFIIIC and its role in RNA polymerase III pre-initiation complex assembly. Nat Commun. 2015;6: 7387. pmid:26060179
  24. 24. Young LS, Rivier DH, Sprague KU. Sequences far downstream from the classical tRNA promoter elements bind RNA polymerase III transcription factors. Mol Cell Biol. 1991;11: 1382–1392. pmid:1996100
  25. 25. Geiduschek EP, Tocchini-Valentini GP. Transcription by RNA polymerase III. Annu Rev Biochem. 1988;57: 873–914. pmid:3052292
  26. 26. Schiffer S, Rösch S, Marchfelder A. Assigning a function to a conserved group of proteins: the tRNA 3′‐processing enzymes. EMBO J. 2002;21: 2769–2777. pmid:12032089
  27. 27. Kruszka K, Barneche F, Guyot R, Ailhas J, Meneau I, Schiffer S, et al. Plant dicistronic tRNA–snoRNA genes: a new mode of expression of the small nucleolar RNAs processed by RNase Z. EMBO J. 2003;22: 621–632. pmid:12554662
  28. 28. Xie K, Minkenberg B, Yang Y. Boosting CRISPR/Cas9 multiplex editing capability with the endogenous tRNA-processing system. Proc Natl Acad Sci. 2015;112: 3570–3575. pmid:25733849
  29. 29. DiCarlo JE, Norville JE, Mali P, Rios X, Aach J, Church GM. Genome engineering in Saccharomyces cerevisiae using CRISPR-Cas systems. Nucleic Acids Res. 2013;41: 4336–4343. pmid:23460208
  30. 30. Jørgensen TR, Park J, Arentshorst M, van Welzen AM, Lamers G, vanKuyk PA, et al. The molecular and genetic basis of conidial pigmentation in Aspergillus niger. Fungal Genet Biol. 2011;48: 544–553. pmid:21277986
  31. 31. Storms R, Zheng Y, Li H, Sillaots S, Martinez-Perez A, Tsang A. Plasmid vectors for protein production, gene expression and molecular manipulations in Aspergillus niger. Plasmid. 2005;53: 191–204. pmid:15848224
  32. 32. Aleksenko A, Clutterbuck AJ. The plasmid replicator AMA1 in Aspergillus nidulans is an inverted duplication of a low-copy-number dispersed genomic repeat. Mol Microbiol. 1996;19: 565–574. pmid:8830247
  33. 33. Ballance DJ, Turner G. Development of a high-frequency transforming vector for Aspergillus nidulans. Gene. 1985;36: 321–331. pmid:3000883
  34. 34. Goosen T, Bloemheuvel G, Gysler C, Bie DA de, Broek HWJ van den, Swart K. Transformation of Aspergillus niger using the homologous orotidine-5′-phosphate-decarboxylase gene. Curr Genet. 1987;11: 499–503. pmid:2836081
  35. 35. Hang YD, Woodams EE. Apple pomace: a potential substrate for citric acid production by Aspergillus niger. Biotechnol Lett. 1984;6: 763–764.
  36. 36. Li Y, Chooi Y-H, Sheng Y, Valentine JS, Tang Y. Comparative characterization of fungal anthracenone and naphthacenedione biosynthetic pathways reveals an α-hydroxylation-dependent claisen-like cyclization catalyzed by a dimanganese thioesterase. J Am Chem Soc. 2011;133: 15773–15785. pmid:21866960
  37. 37. Shrivastav M, De Haro LP, Nickoloff JA. Regulation of DNA double-strand break repair pathway choice. Cell Res. 2008;18: 134–147. pmid:18157161
  38. 38. Carvalho NDSP, Arentshorst M, Jin Kwon M, Meyer V, Ram AFJ. Expanding the ku70 toolbox for filamentous fungi: establishment of complementation vectors and recipient strains for advanced gene analyses. Appl Microbiol Biotechnol. 2010;87: 1463–1473. pmid:20422182
  39. 39. Meyer V, Arentshorst M, El-Ghezal A, Drews A-C, Kooistra R, van den Hondel CAMJJ, et al. Highly efficient gene targeting in the Aspergillus niger kusA mutant. J Biotechnol. 2007;128: 770–775. pmid:17275117
  40. 40. Zhang J, Mao Z, Xue W, Li Y, Tang G, Wang A, et al. Ku80 gene is related to non-homologous end-joining and genome stability in Aspergillus niger. Curr Microbiol. 2011;62: 1342–1346. pmid:21225265
  41. 41. Jacobs JZ, Ciccaglione KM, Tournier V, Zaratiegui M. Implementation of the CRISPR-Cas9 system in fission yeast. Nat Commun. 2014;5: ncomms6344.
  42. 42. Arazoe T, Miyoshi K, Yamato T, Ogawa T, Ohsato S, Arie T, et al. Tailor-made CRISPR/Cas system for highly efficient targeted gene replacement in the rice blast fungus. Biotechnol Bioeng. 2015;112: 2543–2549. pmid:26039904
  43. 43. Zhang C, Meng X, Wei X, Lu L. Highly efficient CRISPR mutagenesis by microhomology-mediated end joining in Aspergillus fumigatus. Fungal Genet Biol. 2016;86: 47–57. pmid:26701308
  44. 44. Schuster M, Schweizer G, Reissmann S, Kahmann R. Genome editing in Ustilago maydis using the CRISPR–Cas system. Fungal Genet Biol. 2016;89: 3–9. pmid:26365384
  45. 45. Pohl C, Kiel JAKW, Driessen AJM, Bovenberg RAL, Nygård Y. CRISPR/Cas9 based genome editing of Penicillium chrysogenum. ACS Synth Biol. 2016;5: 754–764. pmid:27072635
  46. 46. Liu Q, Gao R, Li J, Lin L, Zhao J, Sun W, et al. Development of a genome-editing CRISPR/Cas9 system in thermophilic fungal Myceliophthora species and its application to hyper-cellulase production strain engineering. Biotechnol Biofuels. 2017;10: 1. pmid:28053662
  47. 47. Sugano SS, Suzuki H, Shimokita E, Chiba H, Noji S, Osakabe Y, et al. Genome editing in the mushroom-forming basidiomycete Coprinopsis cinerea, optimized by a high-throughput transformation system. Sci Rep. 2017;7: 1260. pmid:28455526
  48. 48. Matsu-ura T, Baek M, Kwon J, Hong C. Efficient gene editing in Neurospora crassa with CRISPR technology. Fungal Biol Biotechnol. 2015;2: 4. pmid:28955455
  49. 49. Fuller KK, Chen S, Loros JJ, Dunlap JC. Development of the CRISPR/Cas9 system for targeted gene disruption in Aspergillus fumigatus. Eukaryot Cell. 2015;14: 1073–1080. pmid:26318395
  50. 50. Brow DA, Guthrie C. Spliceosomal RNA U6 is remarkably conserved from yeast to mammals. Nature. 1988;334: 213–218. pmid:3041282
  51. 51. Huang Y, Maraia RJ. Comparison of the RNA polymerase III transcription machinery in Schizosaccharomyces pombe, Saccharomyces cerevisiae and human. Nucleic Acids Res. 2001;29: 2675–2690. pmid:11433012
  52. 52. Schwartz CM, Hussain MS, Blenner M, Wheeldon I. Synthetic RNA polymerase III promoters facilitate high-efficiency CRISPR/Cas9-mediated genome editing in Yarrowia lipolytica. ACS Synth Biol. 2016;5: 356–359. pmid:26714206
  53. 53. Liu R, Chen L, Jiang Y, Zhou Z, Zou G. Efficient genome editing in filamentous fungus Trichoderma reesei using the CRISPR/Cas9 system. Cell Discov. 2015;1: celldisc20157.
  54. 54. Weyda I, Yang L, Vang J, Ahring BK, Lübeck M, Lübeck PS. A comparison of Agrobacterium-mediated transformation and protoplast-mediated transformation with CRISPR-Cas9 and bipartite gene targeting substrates, as effective gene targeting tools for Aspergillus carbonarius. J Microbiol Methods. 2017;135: 26–34. pmid:28159628
  55. 55. Huertas P. DNA resection in eukaryotes: deciding how to fix the break. Nat Struct Mol Biol. 2010;17: 11–16. pmid:20051983
  56. 56. Meyer V, Fiedler M, Nitsche B, King R. The Cell Factory Aspergillus Enters the Big Data Era: Opportunities and Challenges for Optimising Product Formation. In: Krull R, Bley T, editors. Filaments in Bioprocesses. Springer International Publishing; 2015. pp. 91–132.
  57. 57. Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25: 955–964. pmid:9023104
  58. 58. Sugahara J, Yachie N, Arakawa K, Tomita M. In silico screening of archaeal tRNA-encoding genes having multiple introns with bulge-helix-bulge splicing motifs. RNA. 2007;13: 671–681. pmid:17369313
  59. 59. Laslett D, Canback B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004;32: 11–16. pmid:14704338
  60. 60. Kinouchi M, Kanaya S, Ikemura T, Kudo Y. Detection of tRNA based on the cloverleaf secondary structure. Genome Inform. 2000;11: 301–302.
  61. 61. Sharp S, DeFranco D, Silberklang M, Hosbach HA, Schmidt T, Kubli E, et al. The initiator tRNA genes of Drosophila melanogaster: evidence for a tRNA pseudogene. Nucleic Acids Res. 1981;9: 5867–5882. pmid:6273811
  62. 62. Hamada M, Sakulich AL, Koduru SB, Maraia RJ. Transcription termination by RNA polymerase III in fission yeast. A genetic and biochemically tractable model system. J Biol Chem. 2000;275: 29076–29081. pmid:10843998
  63. 63. Jühling F, Mörl M, Hartmann RK, Sprinzl M, Stadler PF, Pütz J. tRNAdb 2009: compilation of tRNA sequences and tRNA genes. Nucleic Acids Res. 2009;37: D159–D162. pmid:18957446
  64. 64. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28: 1647–1649. pmid:22543367
  65. 65. Aslanidis C, Jong PJ de. Ligation-independent cloning of PCR products (LIC-PCR). Nucleic Acids Res. 1990;18: 6069–6074. pmid:2235490
  66. 66. Bennett JW, Lasure LL. Appendix B—Growth Media. In: Bennett JW, Lasure LL, editors. More Gene Manipulations in Fungi. San Diego: Academic Press; 1991. pp. 441–458.