INTRODUCTION
Shochu is a traditional Japanese distilled spirit (
1). The black koji fungus
Aspergillus luchuensis and its albino mutant, white koji fungus
Aspergillus luchuensis mut.
kawachii (
A. kawachii), are primarily used for the production of shochu.
A. luchuensis and
A. kawachii produce enzymes that decompose the starch contained in ingredients such as rice, barley, buckwheat, and sweet potato (
2). In addition, they also produce a large amount of citric acid during the fermentation process, which prevents microbial contamination.
A. luchuensis is phylogenetically related to
Aspergillus niger, which has been used for industrial citric acid fermentation (
3–5). Studies investigating citric acid production have been performed for
A. niger with respect to various aspects (
6–8). In previous studies, nonacidifying mutant strains of
A. niger were analyzed and
laeA was found to play a significant role in the production of citric acid and secondary metabolites (
9). In addition,
laeA disruption also caused a significant decrease in citric acid production by
Aspergillus carbonarius, a species closely related to
A. niger (
10). LaeA was initially identified as a regulator of secondary metabolism in
Aspergillus spp. (
11). Subsequently, LaeA has been primarily studied as a regulator of secondary metabolic and development processes in filamentous fungi (
12,
13). A transcriptomic study also supported the idea that
laeA overexpression and disruption caused the production of secondary metabolites to dramatically change in
A. niger (
14). However, why LaeA is required for citric acid production in
A. niger and
A. carbonarius remains to be determined.
In this study, we characterized three
laeA homologous genes, namely,
laeA,
laeA2, and
laeA3, to determine the regulatory mechanism underlying citric acid production in
A. kawachii. Study of gene disruption indicated that
laeA significantly reduced citric acid production; therefore, we further analyzed LaeA-dependent gene expression by cap analysis gene expression (CAGE) and found that
laeA is required for expression of a putative citrate exporter-encoding
cexA gene, which plays a crucial role in citric acid production (
15,
16). Further analysis via complementation test and chromatin immunoprecipitation coupled with quantitative PCR (ChIP-qPCR) indicated that LaeA is required for citric acid production through epigenetic regulation of
cexA in
A. kawachii.
DISCUSSION
In this study, we characterized three LaeA-like methyltransferases, LaeA, LaeA2, and LaeA3, in
A. kawachii to reveal their role in citric acid production. LaeA2 and LaeA3 are not conserved in all
Aspergillus spp. (
Fig. 1B), implying that these two LaeA-like methyltransferases are paralogs. For example, LaeA-like methyltransferases, similar to
A. kawachii LaeA2, were only conserved in the
Aspergillus section
Nigri (
A. carbonarius ITEM 5010,
Aspergillus acidus CBS 106.47, and
Aspergillus tubingensis CBS 134.48), with the exception of
A. niger CBS 513.88. In addition, LaeA-like methyltransferases similar to
A. kawachii LaeA3 were identified only in
A. fumigatus A1163.
Among the three LaeA-like methyltransferase genes, only
laeA disruption caused significant reduction in citric acid production by the
A. kawachii strain (
Fig. 2D). This result is consistent with the results of previous studies showing that LaeA is required for citric acid accumulation in culture supernatants of
A. niger and
A. carbonarius (
9,
10).
As LaeA is known to control secondary metabolism and development in filamentous fungi (
11–13), the Δ
laeA strain showed deficient production of yellow mycelium pigment (
Fig. 2A) and formation of conidia (
Fig. 4B). The Δ
laeA strain exhibited deficient formation of conidia in the presence of sodium chloride but not in the presence of potassium chloride (
Fig. 4B), thereby indicating that LaeA is required for the formation of conidia, particularly in the presence of sodium-specific stress. Sodium is known to demonstrate toxic effects above a certain Na
+/K
+ ratio and has to be excreted from the cytosol to prevent the toxic effects caused by sodium transporters (
27–29). However, there is no significant change in the expression of genes encoding putative sodium transporters Nha1 (AKAW_09133 and AKAW_03679) and Nhx1 (AKAW_07982) in the CAGE data set (Data Set S1). In addition, CAGE could not link the transcription start sites with predicted coding sequences of Ena1 (AKAW_08368 and AKAW_08719), which are also the homologs of sodium transporters.
CAGE identified 1,248 genes differentially expressed between the control and Δ
laeA strains (Data Set S1). These genes included approximately equal numbers of upregulated and downregulated genes (590 and 658, respectively). LaeA should regulate these genes not only directly but also indirectly because they include at least 9 transcription factor gene homologs. For example, the upregulation of
lreB (AKAW_00032) (
30),
tamA (AKAW_02564) (
31), and
dbaA (AKAW_01762) (
32) and the downregulation of
amyR (AKAW_09854) (
33),
zipA (AKAW_05495) (
34),
oefC (AKAW_03045) (
35),
farB (AKAW_07498) (
36),
cpcA (AKAW_00007) (
37), and
devR (AKAW_08419) (
38) were identified in this process.
CAGE identified considerable reduction in the gene expression of the citrate exporter-encoding
cexA gene by
laeA disruption (
Table 2). A further study showed that
cexA overexpression remedied deficient citric acid production by Δ
laeA and Δ
cexA strains at similar levels (
Fig. 6A). These results indicated that the deficient gene expression of
cexA led to the deficient citric acid production by the Δ
laeA strain. For the evaluation of citric acid production,
A. kawachii strains were cultivated in CAP medium after precultivation in M medium. CAP medium was used for enhancing citric acid production concomitant with reduced hyphal growth (
24). Fungal pellet morphology is known to affect the production of organic acids, proteins, and secondary metabolites (
39), and elevated citric acid production occurs with shorter hyphae or hyperbranched phenotypes in
A. niger (
40,
41). The production of mycelial biomass by the Δ
laeA strain did not considerably change during precultivation in M medium compared with that by the control strain (
Fig. 6C); however, it considerably increased after cultivation in CAP medium (
Fig. 2C,
Fig. 6C), thereby indicating that the hyphal growth of the Δ
laeA strain was enhanced compared with that of the control strain under the condition of citric acid production. On the other hand, the production of mycelial biomass by the Δ
cexA plus OE
cexA and Δ
laeA plus OE
cexA strains reduced during precultivation in M medium compared with that by the control strain (
Fig. 6C). In addition, we observed that there were no considerable changes in the pellet sizes of the control, Δ
laeA, Δ
cexA, Δ
cexA plus OE
cexA, and Δ
laeA plus OE
cexA strains after cultivation in CAP medium; however, the surface roughness of the pellets of the Δ
cexA plus OE
cexA and Δ
laeA plus OE
cexA strains was smoother than those of the control strains (data not shown). Unlike in the control strain, these growth phenotypes and morphological changes due to
laeA deletion and
cexA overexpression might be related to citric acid production (
Fig. 6A).
Subcultivation of the control strain from M medium to CAP medium led to elevated transcript levels of
cexA, whereas the Δ
laeA strain lost its transcriptional response, thereby indicating that LaeA is involved in the transcriptional regulation of
cexA (
Fig. 6D). Although details concerning the molecular function of LaeA remain unclear (
11–13), LaeA is known to counteract the heterochromatinization of the promoter region of secondary metabolite gene clusters by histone H3K9 methylation via heterochromatin protein 1 (HepA) and a H3K9 methyltransferase (ClrD) in
A. nidulans (
42). Our findings indicated that the euchromatin structure of the
cexA promoter nearly disappeared with
laeA disruption (
Fig. 7B), and therefore, the heterochromatin level of the
cexA promoter might be enriched by HepA (AKAW_02119) and ClrD (AKAW_07568) orthologs in
A. kawachii (
Fig. 7C). Gene expression levels of AKAW_02119 and AKAW_07568 were not significantly altered in the Δ
laeA strain (Data Set S1), but loss of LaeA might affect histone modification balance and euchromatin/heterochromatin ratios. The molecular mechanism of LaeA-dependent histone modification should be confirmed through additional experiments because the methyl-accepting substrate of LaeA still remains elusive (
43). In addition, whether gene expression of
cexA requires a specific DNA-binding transcriptional factor remains unclear and should be further studied.
In conclusion, LaeA plays a significant role in citric acid production in A. kawachii by controlling cexA expression via histone modification at the cexA promoter region. Because A. kawachii is widely used in the production of shochu and elsewhere in the fermentation industry, our findings are expected to enhance the understanding of citric acid production mechanism(s) and facilitate optimization of strategies for controlling A. kawachii activity.
MATERIALS AND METHODS
Strains and growth conditions.
The
Aspergillus kawachii strains used in this study are listed in
Table 3, and strain SO2 (
44) was used as the parental strain.
For construction and characterization, the strains were grown in minimal medium (1% [wt/vol] glucose, 0.6% [wt/vol] NaNO3, 0.052% [wt/vol] KCl, 0.052% [wt/vol] MgSO4·7H2O, and 0.152% [wt/vol] KH2PO4, plus Hutner’s trace elements [pH 6.5]). The medium was adjusted to the required pH using HCl or NaOH. For the cultivation of sC− and argB− strains, 0.02% (wt/vol) methionine and/or 0.211% (wt/vol) arginine was added to M medium, respectively.
To evaluate acidification occurring on agar medium, strains were grown in YPD with methyl red (2% [wt/vol] glucose, 1% [wt/vol] yeast extract, 2% [wt/vol] peptone, and 2% methyl red solution) prepared as follows: 100 mg of methyl red (Nakalai Tesque, Kyoto, Japan) was dissolved in 100 ml of ethanol and titrated by 0.1% (wt/vol) NaOH solution until observation of an obvious color change from red to yellow; the solution was then sterile filtered (0.2-μm pore size; Toyo Roshi Kaisha, Japan).
To investigate citric acid production, A. kawachii strains were also grown in CAP medium (10% [wt/vol] glucose, 0.3% [wt/vol] (NH4)2SO4, 0.001% [wt/vol] KH2PO4, 0.05% [wt/vol] MgSO4·7H2O, 0.000005% [wt/vol] FeSO4·7H2O, 0.00025% [wt/vol] ZnSO4·5H2O, 0.00006% [wt/vol] CuSO4·5H2O [pH 4.0]). CAP medium was adjusted to the required pH with HCl.
Construction of putative methyltransferase gene disruption strain.
laeA,
laeA2, and
laeA3 were disrupted in
A. kawachii SO2 (
44) by insertion of
argB. A gene replacement cassette encompassing the homology arm at the 5′ end of the putative methyltransferase genes, an
argB selection marker, and the homology arm at the 3′ end of the putative methyltransferase genes was constructed by recombinant PCR using the primer pairs AKlaeX-FC/AKlaeX-R1, argB-F2/argB-R2, and AKlaeX-F3/AKlaeX-RC, respectively (where “X” indicates A, A2, or A3) (
Table 4). For amplification of the
argB gene, the pDC1 plasmid was used as the template DNA (
45). For amplification of other DNA fragments,
A. kawachii IFO 4308 wild-type genomic DNA was used as the template DNA. The resultant DNA fragment amplified with primers AKlaeX-F1 and AKlaeX-R3 was used to transform
A. kawachii SO2, yielding Δ
sC Δ
laeA, Δ
sC Δ
laeA2, and Δ
sC Δ
laeA3 strains. M agar plates lacking arginine were used for the selection of transformants. Introduction of
argB into each methyltransferase gene locus was confirmed by PCR using the primer pair AKlaeX-FC and AKlaeX-RC (Fig. S1).
The SO2 strain was transformed using the
argB cassette to employ the same auxotrophic genetic background strains for comparative study. This
argB gene cassette was generated with PCR using
A. kawachii genomic DNA and pDC1 as the template DNA and was used to transform the SO2 strain, yielding the Δ
sC strain (
Table 3). Transformants were selected on M agar medium lacking arginine.
Complementation of the laeA disruption strain.
For complementation analysis of
laeA disruption using wild-type
laeA, a gene replacement cassette encompassing a homology arm at the 5′ end of
laeA, wild-type
laeA, an
sC selection marker, and a homology arm at the
argB locus was constructed with recombinant PCR using the primer pairs AKlaeA-FC/AKlaeAcomp-R1 and AKlaeAcomp-F2/argB-R2 (
Table 4). For amplification of DNA fragments,
A. kawachii IFO 4308 wild-type genomic DNA and a plasmid carrying tandemly connected
sC and
argB were used as the template DNA. The resultant DNA fragment amplified with primers AKlaeA-F1/argB-R2 was used to transform
laeA disruption, yielding the Δ
laeA plus
laeA strain. Transformants were selected on M agar medium lacking methionine. Introduction of
laeA and
sC into the target locus was confirmed by PCR using primers AKlaeA-FC and argB-R2 (Fig. S2).
The Δ
sC Δ
laeA strain was transformed using the
sC cassette for use of the same auxotrophic genetic background strains for the comparative study. The
sC cassette was generated by PCR using
A. kawachii genomic DNA as the template and primers sC-comp-F and sC-comp-R (
Table 4) and was used to transform the Δ
sC Δ
laeA strain and yield the Δ
laeA strain (
Table 3). Transformants were selected on M agar medium lacking methionine.
Construction of putative citrate exporter gene disruption strain.
cexA was disrupted in
A. kawachii SO2 (
44) by insertion of
argB. A gene replacement cassette encompassing the homology arm at the 5′ end of the
cexA,
argB selection marker, and homology arm at the 3′ end of
cexA was constructed using recombinant PCR with the primer pairs AKcexA-FC/AKcexA-R1, AKcexA-F2/AKcexA-R2, and AKcexA-F3/AKcexA-RC, respectively (
Table 4). For amplification of
argB, plasmid pDC1 was used as the template DNA (
45). For amplification of the other DNA fragment,
A. kawachii IFO 4308 wild-type genomic DNA was used as a template. The resultant DNA fragment was amplified with primers AKcexA-F1 and AKcexA-R3 and was used to transform
A. kawachii SO2 and yield the Δ
sCΔ
cexA strain. M agar plates lacking arginine were used for selection of transformants. Introduction of
argB into the
cexA locus was confirmed by PCR using the primer pair AKcexA-FC and AKcexA-RC (Fig. S3). After confirmation of gene disruption, the Δ
sC Δ
cexA strain was transformed with an
sC cassette to use the same auxotrophic genetic background strains for comparative study. This cassette was synthesized by PCR using primers sC-comp-F and sC-comp-R and
A. kawachii genomic DNA as the template (
Table 4). Transformants were selected on M agar medium lacking methionine.
Construction of the putative citrate exporter overexpression strain.
Plasmid pGS-PgpdA (
46), which carries
A. kawachii sC (
28), was used for overexpression of
cexA. The
cexA gene was amplified by PCR with primers pGSG-cexA-inf-F/pGSG-cexA-inf-R using
A. kawahcii genomic DNA as the template (
Table 4). The amplicon was inserted into the SalI site of pGS-PgpdA, thereby yielding pGS-PgpdA-cexA, which was used to transform the Δ
sC Δ
cexA and Δ
sC Δ
laeA strains, yielding the Δ
cexA plus OE
cexA and Δ
laeA plus OE
cexA strains, respectively (
Table 3). Transformants were selected on M agar medium lacking methionine. Single-copy integration of pGS-PgpdA-cexA was confirmed by real-time RT-PCR (data not shown). In addition, the transformants were transferred to YPD with methyl red; then the appearance of a red border around the colonies was confirmed.
Measurement of extracellular and intracellular organic acids and extracellular glucose.
Levels of extracellular and intracellular organic acids were measured as described previously (
24). Briefly, 2 × 10
7 conidial cells of
A. kawachii were inoculated into 100 ml of M medium, precultured with shaking (180 rpm) at 30°C for 36 h, and then transferred into 50 ml of CAP medium and further cultured with shaking (163 rpm) at 30°C for 12, 24, 48, or 72 h. Culture supernatant was harvested as the extracellular fraction. Mycelia were used for preparation of the intracellular fraction using a hot-water extraction method (
47) with modifications. The freeze-dried and wet mycelial weights were measured as extracellular and intracellular fractions, respectively. Wet mycelia were ground into a powder using mortar and pestle in the presence of liquid nitrogen and then dissolved in 10 ml of hot water (80°C) per 1 g of mycelial powder, vortexed, and centrifuged at 18,800 ×
g at 4°C for 30 min. The supernatant was taken as the intracellular fraction.
To measure organic acid level, extracellular and intracellular fractions were filtered through a polytetrafluoroethylene (PTFE) filter (0.2-μm pore size; Toyo Roshi Kaisha) and analyzed with HPLC on a Prominence HPLC system (Shimadzu, Kyoto, Japan) equipped with a CDD-10AVP conductivity detector (Shimadzu). The organic acids were separated using two tandem Shimadzu Shim-pack SCR-102H columns (internal diameter, 8 by 300 mm; Shimadzu) at 50°C using 4 mM p-toluenesulfonic acid monohydrate as the mobile phase at a flow rate of 0.8 ml/min. The flow rate of the postcolumn reaction solution (4 mM p-toluenesulfonic acid monohydrate, 16 mM bis-Tris, and 80 μM EDTA) was 0.8 ml/min.
To measure the glucose level, the filtrated extracellular fraction was analyzed with the Prominence HPLC system equipped with a RID-10A refractive index detector (Shimadzu). Glucose was separated using a COSMOSIL Sugar-D column (internal diameter, 4.6 by 250 mm; Nacalai Tesque, Kyoto, Japan) at 40°C using acetonitrile and water (3:1) as the mobile phase at a flow rate of 1.0 ml/min.
CAGE analysis.
Total RNA was extracted from mycelia. Briefly, 2 × 107 conidial cells of the A. kawachii strains were inoculated into 100 ml of M medium, precultured with shaking (180 rpm) at 30°C for 36 h, and then transferred to 50 ml of CAP medium and further cultured with shaking (163 rpm) at 30°C for 12 h. The mycelia were ground to a powder as described above. Then, RNA was extracted using RNAiso Plus reagent (TaKaRa Bio, Shiga, Japan). RNA samples were treated with the SV total RNA isolation system (Promega, Madison, WI) according to the manufacturer’s protocol.
Library preparation, sequencing, and data analysis for CAGE (
48) were performed by Kabushiki Kaisha DNAFORM (Kanagawa, Japan). All CAGE experiments were performed three times with RNA samples obtained from independently prepared mycelia. First-strand cDNAs were transcribed to the 5′ end of capped RNAs and attached to CAGE barcode tags, and these tags were sequenced using the NextSeq 500 system (Illumina, San Diego, CA) and mapped to the
A. kawachii IFO 4308 genome (
49) using BWA software (v0.5.9) after discarding ribosomal or non-A/C/G/T-base-containing RNAs. For tag clustering, CAGE tag 5′ coordinates were input for RECLU clustering (
50). The criteria for linking transcriptional start sites and predicted coding sequences were within 600 bp upstream or downstream of the predicted start codon. Triplicate data were analyzed and the expression ratio was also calculated as the log (base 2) ratio through the RECLU pipeline.
Transcriptional analysis.
For RNA extraction from mycelia, conidia (2 × 10
7 cells) of the
A. kawachii control, Δ
laeA, and Δ
laeA plus OE
cexA strains were inoculated into 100 ml of M medium and cultured with shaking (180 rpm) for 36 h at 30°C. After incubation, mycelia were collected and divided into two equal portions, transferred individually to CAP medium, and cultured with shaking (163 rpm) for 12 h at 30°C. Mycelia were ground as described above and RNA was extracted using RNAiso Plus (TaKaRa Bio), and then cDNA was synthesized from total RNA using a PrimeScript Perfect real-time reagent kit (TaKaRa Bio) according to manufacturer’s protocols. Real-time RT-PCR was performed using a Thermal Cycler Dice real-time system MRQ (TaKaRa Bio) with SYBR Premix
Ex Taq II (Tli RNaseH Plus) (TaKaRa Bio). The following primer sets were used: AKcexA-RT-F and AKcexA-RT-R for
cexA and AKactA-RT-F and AKactA-RT-R for
actA (
Table 4) (
51).
ChIP-qPCR.
A. kawachii conidia were cultured as described above. ChIP was performed as previously described (
52) using normal anti-mouse IgG as a negative control (Cosmo Bio, Tokyo, Japan), as well as anti-histone H3 (Medical and Biological Laboratories, Nagoya, Japan), anti-H3K4me3 (Medical and Biological Laboratories), and anti-H3K9me3 (Medical and Biological Laboratories) antibodies. Two micrograms of antibody was used with 200 mg of total protein in each ChIP experiment. DNA quantification was performed with real-time qPCR using SYBR Premix
Ex Taq II (Tli RNaseH Plus) (TaKaRa Bio) and the primer set cexA-ChIP-F and cexA-ChIP-R (
Table 4). Positions of the primers relative to the ATG site of
cexA were +2 to +26 for cexA-ChIP-F and +238 to +257 for cexA-ChIP-R. Relative amounts of DNA (i.e., percent input DNA) were calculated by dividing immunoprecipitated DNA by input DNA.