Discovery of a gene cluster for the biosynthesis of novel cyclic peptide compound, KK-1, in Curvularia clavata.

KK-1, a cyclic depsipeptide with 10 residues produced by a filamentous fungus Curvularia clavata BAUA-2787, is a promising pesticide active compound with high activity against many plant pathogens, especially Botrytis cinerea. As a first step toward the future mass production of KK-1 through synthetic biological approaches, we aimed to identify the genes responsible for the KK-1 biosynthesis. To achieve this, we conducted whole genome sequencing and transcriptome analysis of C. clavata BAUA-2787 to predict the KK-1 biosynthetic gene cluster. We then generated the overexpression and deletion mutants for each cluster gene using our originally developed transformation system for this fungus, and analyzed the KK-1 production and the cluster gene expression levels to confirm their involvement in KK-1 biosynthesis. As a result of these, a region of approximately 71 kb was found, containing 10 open reading frames, which were co-induced during KK-1 production, as a biosynthetic gene cluster. These include kk1B, which encodes nonribosomal peptide synthetase with a domain structure that is consistent with the structural features of KK-1, and kk1F, which encodes a transcription factor. The overexpression of kk1F increased the expression of the entire cluster genes and, consequently, improved KK-1 production, whereas its deletion decreased the expression of the entire cluster genes and almost eliminated KK-1 production, demonstrating that the protein encoded by kk1F regulates the expressions of the other nine cluster genes cooperatively as the pathway-specific transcription factor. Furthermore, the deletion of each cluster gene caused a reduction in KK-1 productivity, indicating that each gene is involved in KK-1 production. The genes kk1A, kk1D, kk1H, and kk1I, which showed a significant decrease in KK-1 productivity due to deletion, were presumed to be directly involved in KK-1 structure formation, including the biosynthesis of the constituent residues. kk1C, kk1E, kk1G, and kk1J, which maintained a certain level of KK-1 productivity despite deletion, were possibly involved in promoting or assisting KK-1 production, such as extracellular transportation and the removal of aberrant units incorporated into the peptide chain.

Comparative analysis of the relationship between translation efficiency and sequence features of endogenous proteins in multiple organisms.

The translation efficiency of protein genes is known to be affected by sequence features. Previous studies have found that various sequence features based on codon usage and mRNA secondary structure contribute to translation efficiency. However, most studies have focused on a specific organism, usually a model organism such as Escherichia coli or Saccharomyces cerevisiae. Here, we investigate whether the relationship between translation efficiency and sequence features is conserved among multiple organisms using publicly available ribosome profiling data and RNA-Seq data. We analyze nine organisms from various taxa: Staphylococcus aureus, five species of Streptomyces, two strains of E. coli, and S. cerevisiae. We reveal that the relationship between translation efficiency and sequence features differs across organisms, partly reflecting their taxonomy. The codon adaptation index shows high correlation in all analyzed organisms. Our study provides an insight into the diversity and commonality of sequence determinants of protein expression in these organisms.

L-Alanine Prototrophic Suppressors Emerge from L-Alanine Auxotroph through Stress-Induced Mutagenesis in Escherichia coli.

In Escherichia coli, L-alanine is synthesized by three isozymes: YfbQ, YfdZ, and AvtA. When an E. coli L-alanine auxotrophic isogenic mutant lacking the three isozymes was grown on L-alanine-deficient minimal agar medium, L-alanine prototrophic mutants emerged considerably more frequently than by spontaneous mutation; the emergence frequency increased over time, and, in an L-alanine-supplemented minimal medium, correlated inversely with L-alanine concentration, indicating that the mutants were derived through stress-induced mutagenesis. Whole-genome analysis of 40 independent L-alanine prototrophic mutants identified 16 and 18 clones harboring point mutation(s) in pyruvate dehydrogenase complex and phosphotransacetylase-acetate kinase pathway, which respectively produce acetyl coenzyme A and acetate from pyruvate. When two point mutations identified in L-alanine prototrophic mutants, in pta (D656A) and aceE (G147D), were individually introduced into the original L-alanine auxotroph, the isogenic mutants exhibited almost identical growth recovery as the respective cognate mutants. Each original- and isogenic-clone pair carrying the pta or aceE mutation showed extremely low phosphotransacetylase or pyruvate dehydrogenase activity, respectively. Lastly, extracellularly-added pyruvate, which dose-dependently supported L-alanine auxotroph growth, relieved the L-alanine starvation stress, preventing the emergence of L-alanine prototrophic mutants. Thus, L-alanine starvation-provoked stress-induced mutagenesis in the L-alanine auxotroph could lead to intracellular pyruvate increase, which eventually induces L-alanine prototrophy.

Whole-Genome Sequence of Monascus purpureus GB-01, an Industrial Strain for Food Colorant Production.

We report the draft genome sequence of Monascus purpureus GB-01, an industrial strain used as a food colorant. De novo assembly of long reads resulted in 121 chromosomal contigs and 1 mitochondrial contig, and sequencing errors were corrected by paired-end short reads. This genome sequence will provide useful information for azaphilone pigments and mycotoxin citrinin biosynthesis.

Production of recombinant extracellular cholesterol esterase using consistently active promoters in Burkholderia stabilis.

Burkholderia stabilis FERMP-21014 produces highly active cholesterol esterase in the presence of fatty acids. To develop an overexpression system for cholesterol esterase production, we carried out RNA sequencing analyses to screen strongly active promoters in FERMP-21014. Based on gene expression consistency analysis, we selected nine genes that were consistently expressed at high levels, following which we constructed expression vectors using their promoter sequences and achieved overproduction of extracellular cholesterol esterase under fatty acid-free conditions. Of the tested promoters, the promoter of BSFP_0720, which encodes the alkyl hydroperoxide reductase subunit AhpC, resulted in the highest cholesterol esterase activity (24.3 U mL-1). This activity level was 243-fold higher than that of the wild-type strain under fatty acid-free conditions. We confirmed that cholesterol esterase was secreted without excessive accumulation within the cells. The gene expression consistency analysis will be useful to screen promoters applicable to the overexpression of other industrially important enzymes.

Developing a codon optimization method for improved expression of recombinant proteins in actinobacteria.

Codon optimization by synonymous substitution is widely used for recombinant protein expression. Recent studies have investigated sequence features for codon optimization based on large-scale expression analyses. However, these studies have been limited to common host organisms such as Escherichia coli. Here, we develop a codon optimization method for Rhodococcus erythropolis, a gram-positive GC-rich actinobacterium attracting attention as an alternative host organism. We evaluate the recombinant protein expression of 204 genes in R. erythropolis with the same plasmid vector. The statistical analysis of these expression data reveals that the mRNA folding energy at 5' regions as well as the codon frequency are important sequence features for codon optimization. Intriguingly, other sequence features such as the codon repetition rate show a different tendency from the previous study on E. coli. We optimize the coding sequences of 12 genes regarding these sequence features, and confirm that 9 of them (75%) achieve increased expression levels compared with wild-type sequences. Especially, for 5 genes whose expression levels for wild-type sequences are small or not detectable, all of them are improved by optimized sequences. These results demonstrate the effectiveness of our codon optimization method in R. erythropolis, and possibly in other actinobacteria.

Biosynthesis of Novel Statins by Combining Heterologous Genes from Xylaria and Aspergillus.

For many secondary metabolites, heterologous synthesis is the definitive step to determine their required biosynthetic genes. Using a multivector expression system in Saccharomyces cerevisiae, we reconstituted not only two natural statins from two fungal species, i.e., lovastatin from Aspergillus terreus and FR901512 from Xylaria grammica, but also new statin structures by mixing their genes. Combinatorial gene exchange experiments revealed the functional promiscuity of two polyketide synthases in A. terreus, lovB, and lovF; they could synthesize FR901512 with Xylaria genes. Key structure determinants of statins are essential accessory genes that are irreplaceable across species.

Knockout of the SREBP system increases production of the polyketide FR901512 in filamentous fungal sp. No. 14919 and lovastatin in Aspergillus terreus ATCC20542.

In the production of useful microbial secondary metabolites, the breeding of strains is generally performed by random mutagenesis. However, because random mutagenesis introduces many mutations into genomic DNA, the causative mutations leading to increased productivity are mostly unknown. Therefore, although gene targeting is more efficient for breeding than random mutagenesis, it is difficult to apply. In this study, a wild-type strain and randomly mutagenized strains of fungal sp. No. 14919, a filamentous fungus producing the HMG-CoA reductase inhibitor polyketide FR901512, were subjected to point mutation analysis based on whole genome sequencing. Among the mutated genes found, mutation of the sterol regulatory element-binding protein (SREBP) cleavage-activating protein (SCAP) had a positive effect on increasing FR901512 productivity. By complementing the SCAP gene in the SCAP-mutated strain, productivity was decreased to the level of the SCAP-intact strain. Conversely, when either the SCAP or SREBP gene was deleted, the productivity was significantly increased. By genomic transcriptional analysis, the expression levels of three enzymes in the ergosterol biosynthesis pathway were shown to be decreased by SCAP mutation. These findings led to the hypothesis that raw materials of polyketides, such as acetyl-CoA and malonyl-CoA, became more available for FR901512 biosynthesis due to depression in sterol biosynthesis caused by knockout of the SREBP system. This mechanism was confirmed in Aspergillus terreus producing the polyketide lovastatin, which is structurally similar to FR901512. Thus, knockout of the SREBP system should be considered significant for increasing the productivities of polyketides, such as HMG-CoA reductase inhibitors, by filamentous fungi.

Complete Genome Sequence of Burkholderia stabilis FERMP-21014.

Cholesterol esterase (EC 3.1.1.13) was identified in a bacterium, Burkholderia stabilis strain FERMP-21014. Here, we report the complete genome sequence of B. stabilis FERMP-21014, which has been used in the commercial production of cholesterol esterase. The genome sequence information may be useful for improving production levels of cholesterol esterase.

Genome Sequence of the Fungal Strain 14919 Producing 3-Hydroxy-3-Methylglutaryl-Coenzyme A Reductase Inhibitor FR901512.

Fungal strain 14919 was originally isolated from a soil sample collected at Mt. Kiyosumi, Chiba Prefecture, Japan. It produces FR901512, a potent and strong 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) reductase inhibitor. The genome sequence of fungal strain 14919 was determined and annotated to improve the productivity of FR901512.

Further enhancement of FR901469 productivity by co-overexpression of cpcA, a cross-pathway control gene, and frbF in fungal sp. No. 11243.

FR901469 is a secondary metabolite with antifungal activity, produced by fungal sp. No. 11243. In our previous study, we constructed the frbF overexpression mutant (TFH2-2) from the wild-type strain. FR901469 productivity of TFH2-2 was 3.4 times higher than that of the wild-type strain. To further enhance FR901469 productivity in TFH2-2, we attempted to find genes from the genome that limited the productivity as bottlenecks in this study. Based on both correlation analysis of gene expression level against FR901469 productivity and genome annotation information, the cross-pathway control gene A (cpcA) was most predicted as the bottleneck. The cpcA and frbF co-overexpression mutant named TFCH3 was then constructed from TFH2-2. As a result, FR901469 productivity of TFCH3 was enhanced at 1.8 times higher than that of TFH2-2. Transcriptome analysis revealed that many genes involved in amino acid biosynthesis and encoding tRNA ligases were significantly upregulated in TFCH3, which implied increase of amino acids as the substrates of FR901469 would be a reason of further productivity enhancement.

Identification of a putative FR901469 biosynthesis gene cluster in fungal sp. No. 11243 and enhancement of the productivity by overexpressing the transcription factor gene frbF.

FR901469 is an antifungal antibiotic produced by fungal sp. No. 11243. Here, we searched for FR901469 biosynthesis genes in the genome of No. 11243. Based on the molecular structure of FR901469 and endogenous functional motifs predicted in each genomic NRPS gene, a putative FR901469 biosynthesis gene cluster harboring the most plausible NRPS gene was identified. A transcription factor gene, designated frbF, was found in the cluster. To improve FR901469 productivity, we constructed a strain in which frbF was overexpressed and named it TFH2-2. FR901469 productivity of TFH2-2 was 3.4 times higher than that of the wild-type strain. Transcriptome analysis revealed that most of the genes in the putative FR901469 biosynthesis gene cluster were upregulated in TFH2-2. It also showed that the expression of genes related to ergosterol biosynthesis, ß-1,3-glucan catabolism, and chitin synthesis was inclined to exhibit significant differences in TFH2-2.

Genome Sequence of Ustilaginoidea virens IPU010, a Rice Pathogenic Fungus Causing False Smut.

Ustilaginoidea virens is a rice pathogenic fungus that causes false smut disease, a disease that seriously damages the yield and quality of the grain. Analysis of the U. virens IPU010 33.6-Mb genome sequence will aid in the understanding of the pathogenicity of the strain, particularly in regard to effector proteins and secondary metabolic genes.

Class of cyclic ribosomal peptide synthetic genes in filamentous fungi.

Ustiloxins were found recently to be the first example of cyclic peptidyl secondary metabolites that are ribosomally synthesized in filamentous fungi. In this work, two function-unknown genes (ustYa/ustYb) in the gene cluster for ustiloxins from Aspergillus flavus were found experimentally to be involved in cyclization of the peptide. Their homologous genes are observed mainly in filamentous fungi and mushrooms. They have two "HXXHC" motifs that might form active sites. Computational genome analyses showed that these genes are frequently located near candidate genes for ribosomal peptide precursors, which have signal peptides at the N-termini and repeated sequences with core peptides for the cyclic portions, in the genomes of filamentous fungi, particularly Aspergilli, as observed in the ustiloxin gene cluster. Based on the combination of the ustYa/ustYb homologous genes and the nearby ribosomal peptide precursor candidate genes, 94 ribosomal peptide precursor candidates that were identified computationally from Aspergilli genome sequences were classified into more than 40 types including a wide variety of core peptide sequences. A set of the predicted ribosomal peptide biosynthetic genes was experimentally verified to synthesize a new cyclic peptide compound, designated as asperipin-2a, which comprises the amino acid sequence in the corresponding precursor gene, distinct from the ustiloxin precursors.

Genome Sequence of Fungal Species No.11243, Which Produces the Antifungal Antibiotic FR901469.

Fungal species No.11243 was originally isolated from a decayed leaf sample collected in Kyoto, Japan. It produces FR901469, a 1,3-beta-glucan synthase inhibitor. The genome sequence of No.11243 was determined and annotated to obtain useful information for improving productivity of the effective antifungal agent FR901469.

Ustiloxins, fungal cyclic peptides, are ribosomally synthesized in Ustilaginoidea virens.

MOTIVATION: Ustiloxins A and B are toxic cyclic tetrapeptides, Tyr-Val/Ala-Ile-Gly (Y-V/A-I-G), that were originally identified from Ustilaginoidea virens, a pathogenic fungus affecting rice plants. Contrary to our report that ustiloxin B is ribosomally synthesized in Aspergillus flavus, a recent report suggested that ustiloxins are synthesized by a non-ribosomal peptide synthetase in U.virens. Thus, we analyzed the U.virens genome, to identify the responsible gene cluster. RESULTS: The biosynthetic gene cluster was identified from the genome of U.virens based on homologies to the ribosomal peptide biosynthetic gene cluster for ustiloxin B identified from A.flavus. It contains a gene encoding precursor protein having five Tyr-Val-Ile-Gly and three Tyr-Ala-Ile-Gly motifs for ustiloxins A and B, respectively, strongly indicating that ustiloxins A and B from U.virens are ribosomally synthesized. AVAILABILITY AND IMPLEMENTATION: Accession codes of the U.virens and A.flavus gene clusters in NCBI are BR001221 and BR001206, respectively. Supplementary data are available at Bioinformatics online.

Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae.

The aspergilli comprise a diverse group of filamentous fungi spanning over 200 million years of evolution. Here we report the genome sequence of the model organism Aspergillus nidulans, and a comparative study with Aspergillus fumigatus, a serious human pathogen, and Aspergillus oryzae, used in the production of sake, miso and soy sauce. Our analysis of genome structure provided a quantitative evaluation of forces driving long-term eukaryotic genome evolution. It also led to an experimentally validated model of mating-type locus evolution, suggesting the potential for sexual reproduction in A. fumigatus and A. oryzae. Our analysis of sequence conservation revealed over 5,000 non-coding regions actively conserved across all three species. Within these regions, we identified potential functional elements including a previously uncharacterized TPP riboswitch and motifs suggesting regulation in filamentous fungi by Puf family genes. We further obtained comparative and experimental evidence indicating widespread translational regulation by upstream open reading frames. These results enhance our understanding of these widely studied fungi as well as provide new insight into eukaryotic genome evolution and gene regulation.

Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus.

Aspergillus fumigatus is exceptional among microorganisms in being both a primary and opportunistic pathogen as well as a major allergen. Its conidia production is prolific, and so human respiratory tract exposure is almost constant. A. fumigatus is isolated from human habitats and vegetable compost heaps. In immunocompromised individuals, the incidence of invasive infection can be as high as 50% and the mortality rate is often about 50% (ref. 2). The interaction of A. fumigatus and other airborne fungi with the immune system is increasingly linked to severe asthma and sinusitis. Although the burden of invasive disease caused by A. fumigatus is substantial, the basic biology of the organism is mostly obscure. Here we show the complete 29.4-megabase genome sequence of the clinical isolate Af293, which consists of eight chromosomes containing 9,926 predicted genes. Microarray analysis revealed temperature-dependent expression of distinct sets of genes, as well as 700 A. fumigatus genes not present or significantly diverged in the closely related sexual species Neosartorya fischeri, many of which may have roles in the pathogenicity phenotype. The Af293 genome sequence provides an unparalleled resource for the future understanding of this remarkable fungus.

Genome sequencing and analysis of Aspergillus oryzae.

The genome of Aspergillus oryzae, a fungus important for the production of traditional fermented foods and beverages in Japan, has been sequenced. The ability to secrete large amounts of proteins and the development of a transformation system have facilitated the use of A. oryzae in modern biotechnology. Although both A. oryzae and Aspergillus flavus belong to the section Flavi of the subgenus Circumdati of Aspergillus, A. oryzae, unlike A. flavus, does not produce aflatoxin, and its long history of use in the food industry has proved its safety. Here we show that the 37-megabase (Mb) genome of A. oryzae contains 12,074 genes and is expanded by 7-9 Mb in comparison with the genomes of Aspergillus nidulans and Aspergillus fumigatus. Comparison of the three aspergilli species revealed the presence of syntenic blocks and A. oryzae-specific blocks (lacking synteny with A. nidulans and A. fumigatus) in a mosaic manner throughout the genome of A. oryzae. The blocks of A. oryzae-specific sequence are enriched for genes involved in metabolism, particularly those for the synthesis of secondary metabolites. Specific expansion of genes for secretory hydrolytic enzymes, amino acid metabolism and amino acid/sugar uptake transporters supports the idea that A. oryzae is an ideal microorganism for fermentation.

Processing sequence annotation data using the Lua programming language.

The data processing language in a graphical software tool that manages sequence annotation data from genome databases should provide flexible functions for the tasks in molecular biology research. Among currently available languages we adopted the Lua programming language. It fulfills our requirements to perform computational tasks for sequence map layouts, i.e. the handling of data containers, symbolic reference to data, and a simple programming syntax. Upon importing a foreign file, the original data are first decomposed in the Lua language while maintaining the original data schema. The converted data are parsed by the Lua interpreter and the contents are stored in our data warehouse. Then, portions of annotations are selected and arranged into our catalog format to be depicted on the sequence map. Our sequence visualization program was successfully implemented, embedding the Lua language for processing of annotation data and layout script. The program is available at http://staff.aist.go.jp/yutaka.ueno/guppy/.