AmpliconDuo: A Split-Sample Filtering Protocol for High-Throughput Amplicon Sequencing of Microbial Communities.

High throughput sequencing (HTSeq) of small ribosomal subunit amplicons has the potential for a comprehensive characterization of microbial community compositions, down to rare species. However, the error-prone nature of the multi-step experimental process requires that the resulting raw sequences are subjected to quality control procedures. These procedures often involve an abundance cutoff for rare sequences or clustering of sequences, both of which limit genetic resolution. Here we propose a simple experimental protocol that retains the high genetic resolution granted by HTSeq methods while effectively removing many low abundance sequences that are likely due to PCR and sequencing errors. According to this protocol, we split samples and submit both halves to independent PCR and sequencing runs. The resulting sequence data is graphically and quantitatively characterized by the discordance between the two experimental branches, allowing for a quick identification of problematic samples. Further, we discard sequences that are not found in both branches ("AmpliconDuo filter"). We show that the majority of sequences removed in this way, mostly low abundance but also some higher abundance sequences, show features expected from random modifications of true sequences as introduced by PCR and sequencing errors. On the other hand, the filter retains many low abundance sequences observed in both branches and thus provides a more reliable census of the rare biosphere. We find that the AmpliconDuo filter increases biological resolution as it increases apparent community similarity between biologically similar communities, while it does not affect apparent community similarities between biologically dissimilar communities. The filter does not distort overall apparent community compositions. Finally, we quantitatively explain the effect of the AmpliconDuo filter by a simple mathematical model.

Generation of a vancomycin-intermediate Staphylococcus aureus (VISA) strain by two amino acid exchanges in VraS.

OBJECTIVES: Staphylococcus aureus is a notorious bacterial pathogen and antibiotic-resistant isolates complicate current treatment strategies. We characterized S. aureus VC40, a laboratory mutant that shows full resistance to glycopeptides (vancomycin and teicoplanin MICs ≥32 mg/L) and daptomycin (MIC = 4 mg/L), to gain deeper insights into the underlying resistance mechanisms. METHODS: Genomics and transcriptomics were performed to characterize changes that might contribute to development of resistance. The mutations in vraS were reconstituted into a closely related parental background. In addition, antimicrobial susceptibility testing, growth analyses, transmission electron microscopy, lysostaphin-induced lysis and autolysis assays were performed to characterize the phenotype of resistant strains. RESULTS: Genome sequencing of strain VC40 revealed 79 mutations in 75 gene loci including genes encoding the histidine kinases VraS and WalK that control cell envelope-related processes. Transcriptomics indicated the increased expression of their respective regulons. Although not reaching the measured MIC for VC40, reconstitution of the L114S and D242G exchanges in VraS(VC40) into the susceptible parental background (S. aureus NCTC 8325) resulted in increased resistance to glycopeptides and daptomycin. The expression of VraS(VC40) led to increased transcription of the cell wall stress stimulon, a thickened cell wall, a decreased growth rate, reduced autolytic activity and increased resistance to lysostaphin-induced lysis in the generated mutant. CONCLUSIONS: We show that a double mutation of a single gene locus, namely vraS, is sufficient to convert the vancomycin-susceptible strain S. aureus NCTC 8325 into a vancomycin-intermediate S. aureus.

Size unlimited markerless deletions by a transconjugative plasmid-system in Bacillus licheniformis.

Conjugative shuttle vectors of the pKVM series, based on an IncP transfer origin and the pMAD vector with a temperature sensitive replication were constructed to establish a markerless gene deletion protocol for Bacilli without natural competence such as the exoenzyme producer Bacillus licheniformis. The pKVM plasmids can be conjugated to strains of B. licheniformis and B. subtilis. For chromosomal gene deletion, regions flanking the target gene are fused and cloned in a pKVM vector prior to conjugative transfer from Escherichia coli to B. licheniformis. Appropriate markers on the vector backbone allow for the identification of the integration at the target locus and thereafter the vector excision, both events taking place via homologous recombination. The functionality of the deletion system was demonstrated with B. licheniformis by a markerless 939 bp in-frame deletion of the yqfD gene and the deletion of a 31 kbp genomic segment carrying a PBSX-like prophage.

Reassessment of the Listeria monocytogenes pan-genome reveals dynamic integration hotspots and mobile genetic elements as major components of the accessory genome.

BACKGROUND: Listeria monocytogenes is an important food-borne pathogen and model organism for host-pathogen interaction, thus representing an invaluable target considering research on the forces governing the evolution of such microbes. The diversity of this species has not been exhaustively explored yet, as previous efforts have focused on analyses of serotypes primarily implicated in human listeriosis. We conducted complete genome sequencing of 11 strains employing 454 GS FLX technology, thereby achieving full coverage of all serotypes including the first complete strains of serotypes 1/2b, 3c, 3b, 4c, 4d, and 4e. These were comparatively analyzed in conjunction with publicly available data and assessed for pathogenicity in the Galleria mellonella insect model. RESULTS: The species pan-genome of L. monocytogenes is highly stable but open, suggesting an ability to adapt to new niches by generating or including new genetic information. The majority of gene-scale differences represented by the accessory genome resulted from nine hyper variable hotspots, a similar number of different prophages, three transposons (Tn916, Tn554, IS3-like), and two mobilizable islands. Only a subset of strains showed CRISPR/Cas bacteriophage resistance systems of different subtypes, suggesting a supplementary function in maintenance of chromosomal stability. Multiple phylogenetic branches of the genus Listeria imply long common histories of strains of each lineage as revealed by a SNP-based core genome tree highlighting the impact of small mutations for the evolution of species L. monocytogenes. Frequent loss or truncation of genes described to be vital for virulence or pathogenicity was confirmed as a recurring pattern, especially for strains belonging to lineages III and II. New candidate genes implicated in virulence function were predicted based on functional domains and phylogenetic distribution. A comparative analysis of small regulatory RNA candidates supports observations of a differential distribution of trans-encoded RNA, hinting at a diverse range of adaptations and regulatory impact. CONCLUSIONS: This study determined commonly occurring hyper variable hotspots and mobile elements as primary effectors of quantitative gene-scale evolution of species L. monocytogenes, while gene decay and SNPs seem to represent major factors influencing long-term evolution. The discovery of common and disparately distributed genes considering lineages, serogroups, serotypes and strains of species L. monocytogenes will assist in diagnostic, phylogenetic and functional research, supported by the comparative genomic GECO-LisDB analysis server (http://bioinfo.mikrobio.med.uni-giessen.de/geco2lisdb).

Genome sequence of Staphylococcus aureus VC40, a vancomycin- and daptomycin-resistant strain, to study the genetics of development of resistance to currently applied last-resort antibiotics.

The increasing emergence of multidrug-resistant Staphylococcus aureus is a problem of global importance. Here, we report the genome of S. aureus VC40, which is resistant to the last-resort antibiotics vancomycin and daptomycin. Its genome sequence will allow insights into the mechanisms that convey full resistance to these compounds.

Complete genome sequences of the chemolithoautotrophic Oligotropha carboxidovorans strains OM4 and OM5.

We report on genome sequencing of Oligotropha carboxidovorans strain OM4 and resequencing of strain OM5. The genomes of both are composed of one chromosome and two plasmids. The presence of two plasmids in the OM5 genome is inconsistent with the previously published sequence, for which only one plasmid was described (D. Paul, S. Bridges, S. Burgess, Y. Dandass, and M. Lawrence, BMC Genomics 11:511, 2010).

Whole-genome sequence of the transformable Neisseria meningitidis serogroup A strain WUE2594.

Serogroup A meningococci are a leading cause of bacterial meningitis in children and young adults worldwide. However, the genetic basis of serogroup A strains' virulence and their epidemiological properties remain poorly understood. Therefore, we sequenced the complete genome of the transformable Neisseria meningitidis serogroup A strain WUE2594.

Genomic features and insights into the biology of Mycoplasma fermentans.

We present the complete genomic sequence of Mycoplasma fermentans, an organism suggested to be associated with the pathogenesis of rheumatoid arthritis in humans. The genome is composed of 977,524 bp and has a mean G+C content of 26.95 mol%. There are 835 predicted protein-coding sequences and a mean coding density of 87.6 %. Functions have been assigned to 58.8 % of the predicted protein-coding sequences, while 18.4 % of the proteins are conserved hypothetical proteins and 22.8 % are hypothetical proteins. In addition, there are two complete rRNA operons and 36 tRNA coding sequences. The largest gene families are the ABC transporter family (42 members), and the functionally heterogeneous group of lipoproteins (28 members), which encode the characteristic prokaryotic cysteine 'lipobox'. Protein secretion occurs through a pathway consisting of SecA, SecD, SecE, SecG, SecY and YidC. Some highly conserved eubacterial proteins, such as GroEL and GroES, are notably absent. The genes encoding DnaK-DnaJ-GrpE and Tig, forming the putative complex of chaperones, are intact, providing the only known control over protein folding. Eighteen nucleases and 17 proteases and peptidases were detected as well as three genes for the thioredoxin-thioreductase system. Overall, this study presents insights into the physiology of M. fermentans, and provides several examples of the genetic basis of systems that might function as virulence factors in this organism.

Diversity of glycosyl hydrolases from cellulose-depleting communities enriched from casts of two earthworm species.

The guts and casts of earthworms contain microbial assemblages that process large amounts of organic polymeric substrates from plant litter and soil; however, the enzymatic potential of these microbial communities remains largely unexplored. In the present work, we retrieved carbohydrate-modifying enzymes through the activity screening of metagenomic fosmid libraries from cellulose-depleting microbial communities established with the fresh casts of two earthworm species, Aporrectodea caliginosa and Lumbricus terrestris, as inocula. Eight glycosyl hydrolases (GHs) from the A. caliginosa-derived community were multidomain endo-beta-glucanases, beta-glucosidases, beta-cellobiohydrolases, beta-galactosidase, and beta-xylosidases of known GH families. In contrast, two GHs derived from the L. terrestris microbiome had no similarity to any known GHs and represented two novel families of beta-galactosidases/alpha-arabinopyranosidases. Members of these families were annotated in public databases as conserved hypothetical proteins, with one being structurally related to isomerases/dehydratases. This study provides insight into their biochemistry, domain structures, and active-site architecture. The two communities were similar in bacterial composition but significantly different with regard to their eukaryotic inhabitants. Further sequence analysis of fosmids and plasmids bearing the GH-encoding genes, along with oligonucleotide usage pattern analysis, suggested that those apparently originated from Gammaproteobacteria (pseudomonads and Cellvibrio-like organisms), Betaproteobacteria (Comamonadaceae), and Alphaproteobacteria (Rhizobiales).

Phylogenetic diversity and metabolic potential revealed in a glacier ice metagenome.

The largest part of the Earth's microbial biomass is stored in cold environments, which represent almost untapped reservoirs of novel species, processes, and genes. In this study, the first metagenomic survey of the metabolic potential and phylogenetic diversity of a microbial assemblage present in glacial ice is presented. DNA was isolated from glacial ice of the Northern Schneeferner, Germany. Pyrosequencing of this DNA yielded 1,076,539 reads (239.7 Mbp). The phylogenetic composition of the prokaryotic community was assessed by evaluation of a pyrosequencing-derived data set and sequencing of 16S rRNA genes. The Proteobacteria (mainly Betaproteobacteria), Bacteroidetes, and Actinobacteria were the predominant phylogenetic groups. In addition, isolation of psychrophilic microorganisms was performed, and 13 different bacterial isolates were recovered. Analysis of the 16S rRNA gene sequences of the isolates revealed that all were affiliated to the predominant groups. As expected for microorganisms residing in a low-nutrient environment, a high metabolic versatility with respect to degradation of organic substrates was detected by analysis of the pyrosequencing-derived data set. The presence of autotrophic microorganisms was indicated by identification of genes typical for different ways of carbon fixation. In accordance with the results of the phylogenetic studies, in which mainly aerobic and facultative aerobic bacteria were detected, genes typical for central metabolism of aerobes were found. Nevertheless, the capability of growth under anaerobic conditions was indicated by genes involved in dissimilatory nitrate/nitrite reduction. Numerous characteristics for metabolic adaptations associated with a psychrophilic lifestyle, such as formation of cryoprotectants and maintenance of membrane fluidity by the incorporation of unsaturated fatty acids, were detected. Thus, analysis of the glacial metagenome provided insights into the microbial life in frozen habitats on Earth, thereby possibly shedding light onto microbial life in analogous extraterrestrial environments.

Rhizobium sp. strain NGR234 possesses a remarkable number of secretion systems.

Rhizobium sp. strain NGR234 is a unique alphaproteobacterium (order Rhizobiales) that forms nitrogen-fixing nodules with more legumes than any other microsymbiont. We report here that the 3.93-Mbp chromosome (cNGR234) encodes most functions required for cellular growth. Few essential functions are encoded on the 2.43-Mbp megaplasmid (pNGR234b), and none are present on the second 0.54-Mbp symbiotic plasmid (pNGR234a). Among many striking features, the 6.9-Mbp genome encodes more different secretion systems than any other known rhizobia and probably most known bacteria. Altogether, 132 genes and proteins are linked to secretory processes. Secretion systems identified include general and export pathways, a twin arginine translocase secretion system, six type I transporter genes, one functional and one putative type III system, three type IV attachment systems, and two putative type IV conjugation pili. Type V and VI transporters were not identified, however. NGR234 also carries genes and regulatory networks linked to the metabolism of a wide range of aromatic and nonaromatic compounds. In this way, NGR234 can quickly adapt to changing environmental stimuli in soils, rhizospheres, and plants. Finally, NGR234 carries at least six loci linked to the quenching of quorum-sensing signals, as well as one gene (ngrI) that possibly encodes a novel type of autoinducer I molecule.

Genome sequence of Desulfobacterium autotrophicum HRM2, a marine sulfate reducer oxidizing organic carbon completely to carbon dioxide.

Sulfate-reducing bacteria (SRB) belonging to the metabolically versatile Desulfobacteriaceae are abundant in marine sediments and contribute to the global carbon cycle by complete oxidation of organic compounds. Desulfobacterium autotrophicum HRM2 is the first member of this ecophysiologically important group with a now available genome sequence. With 5.6 megabasepairs (Mbp) the genome of Db. autotrophicum HRM2 is about 2 Mbp larger than the sequenced genomes of other sulfate reducers (SRB). A high number of genome plasticity elements (> 100 transposon-related genes), several regions of GC discontinuity and a high number of repetitive elements (132 paralogous genes Mbp(-1)) point to a different genome evolution when comparing with Desulfovibrio spp. The metabolic versatility of Db. autotrophicum HRM2 is reflected in the presence of genes for the degradation of a variety of organic compounds including long-chain fatty acids and for the Wood-Ljungdahl pathway, which enables the organism to completely oxidize acetyl-CoA to CO(2) but also to grow chemolithoautotrophically. The presence of more than 250 proteins of the sensory/regulatory protein families should enable Db. autotrophicum HRM2 to efficiently adapt to changing environmental conditions. Genes encoding periplasmic or cytoplasmic hydrogenases and formate dehydrogenases have been detected as well as genes for the transmembrane TpII-c(3), Hme and Rnf complexes. Genes for subunits A, B, C and D as well as for the proposed novel subunits L and F of the heterodisulfide reductases are present. This enzyme is involved in energy conservation in methanoarchaea and it is speculated that it exhibits a similar function in the process of dissimilatory sulfate reduction in Db. autotrophicum HRM2.

The genome of Clostridium kluyveri, a strict anaerobe with unique metabolic features.

Clostridium kluyveri is unique among the clostridia; it grows anaerobically on ethanol and acetate as sole energy sources. Fermentation products are butyrate, caproate, and H2. We report here the genome sequence of C. kluyveri, which revealed new insights into the metabolic capabilities of this well studied organism. A membrane-bound energy-converting NADH:ferredoxin oxidoreductase (RnfCDGEAB) and a cytoplasmic butyryl-CoA dehydrogenase complex (Bcd/EtfAB) coupling the reduction of crotonyl-CoA to butyryl-CoA with the reduction of ferredoxin represent a new energy-conserving module in anaerobes. The genes for NAD-dependent ethanol dehydrogenase and NAD(P)-dependent acetaldehyde dehydrogenase are located next to genes for microcompartment proteins, suggesting that the two enzymes, which are isolated together in a macromolecular complex, form a carboxysome-like structure. Unique for a strict anaerobe, C. kluyveri harbors three sets of genes predicted to encode for polyketide/nonribosomal peptide synthetase hybrides and one set for a nonribosomal peptide synthetase. The latter is predicted to catalyze the synthesis of a new siderophore, which is formed under iron-deficient growth conditions.

Metagenomics of the deep Mediterranean, a warm bathypelagic habitat.

BACKGROUND: Metagenomics is emerging as a powerful method to study the function and physiology of the unexplored microbial biosphere, and is causing us to re-evaluate basic precepts of microbial ecology and evolution. Most marine metagenomic analyses have been nearly exclusively devoted to photic waters. METHODOLOGY/PRINCIPAL FINDINGS: We constructed a metagenomic fosmid library from 3,000 m-deep Mediterranean plankton, which is much warmer (approximately 14 degrees C) than waters of similar depth in open oceans (approximately 2 degrees C). We analyzed the library both by phylogenetic screening based on 16S rRNA gene amplification from clone pools and by sequencing both insert extremities of ca. 5,000 fosmids. Genome recruitment strategies showed that the majority of high scoring pairs corresponded to genomes from Rhizobiales within the Alphaproteobacteria, Cenarchaeum symbiosum, Planctomycetes, Acidobacteria, Chloroflexi and Gammaproteobacteria. We have found a community structure similar to that found in the aphotic zone of the Pacific. However, the similarities were significantly higher to the mesopelagic (500-700 m deep) in the Pacific than to the single 4000 m deep sample studied at this location. Metabolic genes were mostly related to catabolism, transport and degradation of complex organic molecules, in agreement with a prevalent heterotrophic lifestyle for deep-sea microbes. However, we observed a high percentage of genes encoding dehydrogenases and, among them, cox genes, suggesting that aerobic carbon monoxide oxidation may be important in the deep ocean as an additional energy source. CONCLUSIONS/SIGNIFICANCE: The comparison of metagenomic libraries from the deep Mediterranean and the Pacific ALOHA water column showed that bathypelagic Mediterranean communities resemble more mesopelagic communities in the Pacific, and suggests that, in the absence of light, temperature is a major stratifying factor in the oceanic water column, overriding pressure at least over 4000 m deep. Several chemolithotrophic metabolic pathways could supplement organic matter degradation in this most depleted habitat.

Comparative analysis of the complete genome sequence of the plant growth-promoting bacterium Bacillus amyloliquefaciens FZB42.

Bacillus amyloliquefaciens FZB42 is a Gram-positive, plant-associated bacterium, which stimulates plant growth and produces secondary metabolites that suppress soil-borne plant pathogens. Its 3,918-kb genome, containing an estimated 3,693 protein-coding sequences, lacks extended phage insertions, which occur ubiquitously in the closely related Bacillus subtilis 168 genome. The B. amyloliquefaciens FZB42 genome reveals an unexpected potential to produce secondary metabolites, including the polyketides bacillaene and difficidin. More than 8.5% of the genome is devoted to synthesizing antibiotics and siderophores by pathways not involving ribosomes. Besides five gene clusters, known from B. subtilis to mediate nonribosomal synthesis of secondary metabolites, we identified four giant gene clusters absent in B. subtilis 168. The pks2 gene cluster encodes the components to synthesize the macrolactin core skeleton.

Complete nucleotide sequence of the 113-kilobase linear catabolic plasmid pAL1 of Arthrobacter nitroguajacolicus Rü61a and transcriptional analysis of genes involved in quinaldine degradation.

The nucleotide sequence of the linear catabolic plasmid pAL1 from the 2-methylquinoline (quinaldine)-degrading strain Arthrobacter nitroguajacolicus Rü61a comprises 112,992 bp. A total of 103 open reading frames (ORFs) were identified on pAL1, 49 of which had no annotatable function. The ORFs were assigned to the following functional groups: (i) catabolism of quinaldine and anthranilate, (ii) conjugation, and (iii) plasmid maintenance and DNA replication and repair. The genes for conversion of quinaldine to anthranilate are organized in two operons that include ORFs presumed to code for proteins involved in assembly of the quinaldine-4-oxidase holoenzyme, namely, a MobA-like putative molybdopterin cytosine dinucleotide synthase and an XdhC-like protein that could be required for insertion of the molybdenum cofactor. Genes possibly coding for enzymes involved in anthranilate degradation via 2-aminobenzoyl coenzyme A form another operon. These operons were expressed when cells were grown on quinaldine or on aromatic compounds downstream in the catabolic pathway. Single-stranded 3' overhangs of putative replication intermediates of pAL1 were predicted to form elaborate secondary structures due to palindromic and superpalindromic terminal sequences; however, the two telomeres appear to form different structures. Sequence analysis of ORFs 101 to 103 suggested that pAL1 codes for one or two putative terminal proteins, presumed to be covalently bound to the 5' termini, and a multidomain telomere-associated protein (Tap) comprising 1,707 amino acids. Even if the putative proteins encoded by ORFs 101 to 103 share motifs with the Tap and terminal proteins involved in telomere patching of Streptomyces linear replicons, their overall sequences and domain structures differ significantly.

Genome sequence of the bioplastic-producing "Knallgas" bacterium Ralstonia eutropha H16.

The H(2)-oxidizing lithoautotrophic bacterium Ralstonia eutropha H16 is a metabolically versatile organism capable of subsisting, in the absence of organic growth substrates, on H(2) and CO(2) as its sole sources of energy and carbon. R. eutropha H16 first attracted biotechnological interest nearly 50 years ago with the realization that the organism's ability to produce and store large amounts of poly[R-(-)-3-hydroxybutyrate] and other polyesters could be harnessed to make biodegradable plastics. Here we report the complete genome sequence of the two chromosomes of R. eutropha H16. Together, chromosome 1 (4,052,032 base pairs (bp)) and chromosome 2 (2,912,490 bp) encode 6,116 putative genes. Analysis of the genome sequence offers the genetic basis for exploiting the biotechnological potential of this organism and provides insights into its remarkable metabolic versatility.

Structural and functional characterization of three polyketide synthase gene clusters in Bacillus amyloliquefaciens FZB 42.

Although bacterial polyketides are of considerable biomedical interest, the molecular biology of polyketide biosynthesis in Bacillus spp., one of the richest bacterial sources of bioactive natural products, remains largely unexplored. Here we assign for the first time complete polyketide synthase (PKS) gene clusters to Bacillus antibiotics. Three giant modular PKS systems of the trans-acyltransferase type were identified in Bacillus amyloliquefaciens FZB 42. One of them, pks1, is an ortholog of the pksX operon with a previously unknown function in the sequenced model strain Bacillus subtilis 168, while the pks2 and pks3 clusters are novel gene clusters. Cassette mutagenesis combined with advanced mass spectrometric techniques such as matrix-assisted laser desorption ionization-time of flight mass spectrometry and liquid chromatography-electrospray ionization mass spectrometry revealed that the pks1 (bae) and pks3 (dif) gene clusters encode the biosynthesis of the polyene antibiotics bacillaene and difficidin or oxydifficidin, respectively. In addition, B. subtilis OKB105 (pheA sfp(0)), a transformant of the B. subtilis 168 derivative JH642, was shown to produce bacillaene, demonstrating that the pksX gene cluster directs the synthesis of that polyketide. The GenBank accession numbers for gene clusters pks1(bae), pks2, and pks3(dif) are AJ 634060.2, AJ 6340601.2, and AJ 6340602.2, respectively.

FLO11 mediated filamentous growth of the yeast Saccharomyces cerevisiae depends on the expression of the ribosomal RPS26 genes.

The RPS26A and RPS26B isogenes of Saccharomyces cerevisiae encode two almost identical proteins of the small 40S ribosomal subunit, which differ by only two amino acid residues. Growth of an rps26BDelta mutant strain is normal, whereas an rps26ADelta strain displays a reduced growth rate and increased sensitivity towards the specific translational inhibitor paromomycin. An rps26ADelta rps26BDelta double mutant strain is inviable. RPS26A but not RPS26B is required for haploid adhesive and diploid pseudohyphal growth mediated by FLO11, which encodes an adhesion. The RPS26A and RPS26B transcripts make up about 70 and 30% of the cellular RPS26 mRNA, respectively. Overexpression of RPS26B, as well as an RPS26B open reading frame driven by the RPS26A promoter, complements the rps26ADelta deletion and restores haploid invasive growth as well as diploid pseudohyphal growth. These results suggest that the two proteins are functionally interchangeable. FLO11-lacZ activity is not present in haploid rps26ADelta yeast mutant strains, even though FLO11 mRNA levels are not reduced. This suggests that the amount of Rps26p is critical for accurate translation of the FLO11 mRNA, and therefore for the dimorphic switch of the bakera9s yeast from a single cell yeast to an adhesive filamentous growth form.

Evolution of 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase-encoding genes in the yeast Saccharomyces cerevisiae.

The shikimate pathway resulting in three aromatic amino acids is initiated in different organisms by two and three 3-deoxy-d-arabino-heptulosonate-7-phosphate synthases, respectively. Aro3p and Aro4p are the yeast enzymes feedback-inhibited by phenylalanine and tyrosine, respectively. A yeast strain deficient in the general control transcriptional regulatory system of amino acid biosynthesis is unable to live in the presence of high amounts of phenylalanine and tyrosine. Here, we show that this yeast strain can be rescued by the expression of aroH from Escherichia coli encoding the tryptophan-regulated AroH as third isoenzyme. Yeast carrying Ec AroH as the only enzyme for the initial step of the shikimate pathway can grow in the absence of tryptophan. Without aromatic amino acids, this yeast strain survives only when the yeast ARO3 promoter instead of the ARO4 promoter drives E. coli aroH. The detailed analysis of Aro3p and Aro4p revealed a triple feedback control by tyrosine/phenylalanine and tryptophan. Dissecting this control allowed engineering of Aro4p S195A as an enzyme, which is inhibited like AroH only by tryptophan. In addition, Aro4p variants were constructed that show an equally strong inhibition by tyrosine and tryptophan (Aro4p P165G Q302R) and in which the regulation by tyrosine and tryptophan was reversed (Aro4p P165G). Our data suggest that yeast possesses only two instead of three isogenes encoding 3-deoxy-D-arabino-heptulosonate-7-phosphate synthases because both isoenzymes can be fine tuned by tryptophan as additional effector and because transcriptional regulation by the general control system can be induced as backup when aromatic amino acids in the environment are imbalanced.