Identification and characterization of a galacturonic acid transporter from Neurospora crassa and its application for Saccharomyces cerevisiae fermentation processes
© Benz et al.; licensee BioMed Central Ltd. 2014
Received: 5 November 2013
Accepted: 15 January 2014
Published: 6 February 2014
Pectin-rich agricultural wastes potentially represent favorable feedstocks for the sustainable production of alternative energy and bio-products. Their efficient utilization requires the conversion of all major constituent sugars. The current inability of the popular fermentation host Saccharomyces cerevisiae to metabolize the major pectic monosaccharide D-galacturonic acid (D-GalA) significantly hampers these efforts. While it has been reasoned that the optimization of cellular D-GalA uptake will be critical for the engineering of D-GalA utilization in yeast, no dedicated eukaryotic transport protein has been biochemically described. Here we report for the first time such a eukaryotic D-GalA transporter and characterize its functionality in S. cerevisiae.
We identified and characterized the D-GalA transporter GAT-1 out of a group of candidate genes obtained from co-expression analysis in N. crassa. The N. crassa Δgat-1 deletion strain is substantially affected in growth on pectic substrates, unable to take up D-GalA, and impaired in D-GalA-mediated signaling events. Moreover, expression of a gat-1 construct in yeast conferred the ability for strong high-affinity D-GalA accumulation rates, providing evidence for GAT-1 being a bona fide D-GalA transport protein. By recombinantly co-expressing D-galacturonate reductase or uronate dehydrogenase in yeast we furthermore demonstrated a transporter-dependent conversion of D-GalA towards more reduced (L-galactonate) or oxidized (meso-galactaric acid) downstream products, respectively, over a broad concentration range.
By utilizing the novel D-GalA transporter GAT-1 in S. cerevisiae we successfully generated a transporter-dependent uptake and catalysis system for D-GalA into two products with high potential for utilization as platform chemicals. Our data thereby provide a considerable first step towards a more complete utilization of biomass for biofuel and value-added chemicals production.
KeywordsPectin D-galacturonic acid Neurospora crassa Sugar transport Saccharomyces cerevisiae Metabolic engineering Bioconversion Meso-galactaric acid L-galactonic acid
The development of sustainable energy and chemical production methods are among the greatest societal challenges facing our generation. Biofuel and green chemistry efforts utilizing renewable lignocellulosic plant biomass feedstocks contribute substantially towards these goals. However, to maximize economic viability, these processes must implement the complete conversion of all major constituent feedstock sugars. Pectin represents a main plant cell-wall polysaccharide that, unlike cellulose and hemicellulose, is not commonly utilized, especially in Saccharomyces cerevisiae (S. cerevisiae), for biofuel and platform chemical production. Since pectin is most abundant in the primary cell walls of soft and growing tissues, agricultural waste streams from fruits and vegetables are particularly pectin-rich. Crops with a high pectin content of 20 to >40% include sugar beet pulp, citrus peels, and apple pomace [1–3], and the use of their respective waste streams for the production of biofuels has been attempted in several studies [4–8]. These residues are naturally largely devoid of lignin, already collected at the processing plant, and partly pretreated during the sugar and juice extraction procedure, warranting their high potential as biofuel feedstocks [9, 10]. The pectic cell wall fraction can be inexpensively and efficiently hydrolyzed into its component monosaccharides by enzymatic processes, as it is considerably less recalcitrant than other plant polysaccharides, such as cellulose. Despite these advantages, pectin-rich sugar beet pulp and citrus peels are currently disposed of in landfills or dried for use as a low-value cattle feed, an expensive procedure that significantly reduces profit margins. Utilizing these waste streams would augment current biofuel production efforts without creating direct or indirect land-use changes.
Pectin comprises a heterogeneous polysaccharide family with four main structural classes: homogalacturonan (HG), rhamnogalacturonan I (RG-I), and the substituted HGs rhamnogalacturonan II (RG-II) and xylogalacturonan (XG) (for review, see [11, 12]). All pectic structures have α-D-galacturonic acid (D-GalA) as a major backbone component linked at the O-1 and the O-4 positions. This monosaccharide alone can comprise up to 70% of pectin; the remaining sugars include L-arabinose (L-Ara), D-galactose (D-Gal), L-rhamnose (L-Rha), and D-xylose (D-Xyl).
Engineering an organism to utilize D-GalA promises an attractive route for converting a currently underutilized, relatively abundant sugar into a biofuel or platform chemical. The yeast S. cerevisiae is currently the most attractive production host and remains the most popular microorganism for industrial fermentation strategies to produce bioethanol. Its advantages include a high tolerance to growth inhibitors from lignocellulose hydrolysates as well as ethanol, the ability to withstand low pH conditions that eradicate many bacterial contaminants, fast fermentation kinetics, and the suitability for many rounds of recycling [13, 14]. Ample engineering efforts have already been undertaken to utilize glucose, xylose, and arabinose. Unfortunately, S. cerevisiae cannot metabolize D-GalA, because it lacks the genes encoding a catabolic pathway [15–17]. When fermenting hydrolysates from pectin-rich feedstocks, this could, therefore, lead to the accumulation of D-GalA in the broth, which was shown to be inhibitory to the fermentation of D-Gal, L-Ara, and D-Xyl . A possible approach to overcome this problem is metabolic engineering. In this case, the genes encoding the necessary enzymes for D-GalA metabolism derived from organisms capable of utilizing this sugar could be heterologously expressed in yeast. Such pathways have been described in bacteria, such as Escherichia coli and Erwinia[19–22], as well as more recently in filamentous fungi, such as Aspergillus niger, Trichoderma reesei (anamorph of Hypocrea jecorina) and Botrytis cinerea[23–28] (also see  for further review).
A recent attempt to engineer S. cerevisiae strains carrying a bacterial D-GalA catabolic pathway met with considerable challenges in expressing functional enzymes [16, 30]. Moreover, even though D-GalA was demonstrated to be able to enter S. cerevisiae cells under certain conditions (through an as-yet unidentified, low-affinity and channel-like pore at acidic pH close to its pKa of about 3.5 ), it was reasoned that the optimization of D-GalA transport will be essential for the successful engineering of D-GalA utilization in yeast [14, 16]. Although prokaryotic D-GalA transport systems are well known [32–35], these are notoriously difficult to express functionally in a eukaryotic host. However, so far no such transport protein has been described in a Eukaryote. Here we report such a eukaryotic D-GalA transporter. The corresponding gene was identified through a transcriptomics analysis of pectin degradation by the model filamentous fungus Neurospora crassa (N. crassa) , followed by a genetic analysis and the biochemical characterization of the encoded protein. Moreover, we show the transporter to be functional when heterologously expressed in S. cerevisiae and useful for the conversion of D-GalA to downstream products. Our findings are therefore an important step towards the effective utilization of pectin-rich feedstocks for the production of platform chemicals or biofuels.
Identification of NCU00988 from N. crassa as a candidate D-galacturonic acid transporter
The N. crassa deletion strain of NCU00988 (Δgat-1) is impaired in growth on pectins and unable to take up D-galacturonic acid
We next set out to test whether the growth inhibition on pectin could be directly linked to loss of D-GalA transport activity of the deletion strain. For this, sucrose pre-grown cultures of WT and Δgat-1 were first transferred for an additional 4 h to Vogel’s solution containing 0.5% pectin as the sole carbon source to induce the pectic response, before being transferred to the reaction solution containing 90 μM of either D-Gal, D-GalA, or glucuronic acid (D-GlcA). The residual concentrations of the monosaccharides in the reaction supernatant were measured by high pH-anion exchange-chromatography with pulsed amperometric detection (HPAEC-PAD) and plotted over time in Figures 2B and C. Whereas the transport speed for D-Gal was found to be WT-like in the Δgat-1 strain (Figure 2C), the transport of both D-GalA and D-GlcA was almost completely abolished (Figure 2B). These findings provide evidence that the physiological function of GAT-1 is indeed the transport of D-GalA (and the minor component D-GlcA) during fungal growth on pectin.
GAT-1 is required for D-GalA-mediated signaling processes in N. crassa
The transport of D-galacturonic acid and quinic acid are separate physiological processes
The (generic) annotation of the N. crassa genome by the Broad Institute classifies the gat-1 gene (NCU00988) as “MFS quinate transporter”, which presumably can be attributed to the fact that GAT-1 and the quinate permease (NCU06026; quinate-Y[39–41]; Broad Institute of Harvard and MIT [http://www.broadinstitute.org/]) indeed share considerable homology (39% identity, 54% similarity ; Additional file 2: Figure S2). Intriguingly, D-GalA and quinate are also similar in their structure and conformation: quinic acid (qa) being a cyclic polyol, with the carboxyl and at least one hydroxyl substituent in identical equatorial conformation and ring position when presented in six-membered ring chair form structures as shown in Additional file 3: Figure S3A. To test whether the two transporters are also functionally redundant, we performed transport assays with N. crassa WT mycelia that had been pre-grown in sucrose and were then transferred for 4 h to Vogel’s medium containing both 0.5% pectin and 100 μM qa for the parallel induction of both uptake systems. When the uptake of both substrates was measured over the course of 40 minutes, the absence of GAT-1 did not affect qa uptake rates (Additional file 3: Figure S3B), just as the absence of the quinate permease did not negatively affect the uptake of D-GalA (but rather led to a slightly accelerated D-GalA consumption; Additional file 3: Figure S3C). In contrast, the deletion strain for the quinate permease (Δ6026) was almost completely unable to take up qa, indicating that it is the only physiologically relevant transporter for this compound. Therefore, despite their homology, our data suggest a distinct, not overlapping role for GAT-1 and the quinate permease in N. crassa.
Heterologous expression of GAT-1 in S. cerevisiae confers the ability to take up D-galacturonic acid with high affinity
Engineering of S. cerevisiae for the production of L-galactonic acid and meso-galactaric acid from D-galacturonic acid
Passive uptake of D-GalA by WT yeast has been observed at high concentrations of this substrate (50 to 200 mM), which is presumably within the concentration range found in hydrolysates from pectin-rich feedstocks. This uptake phenomenon was shown to be largely dependent on the concentration of D-GalA and pH . The expression of a primary or secondary active transporter could plausibly improve D-GalA transport rates even at high concentrations and/or sustain import as extracellular concentrations decrease. To test whether heterologous expression of the high-affinity transporter GAT-1 will also have beneficial effects in high-D-GalA media conditions we measured L-GalOA and GalAA production in our bioconversion strains at 100 mM D-GalA and pH 6.0. After 1 h incubation, we observed an average 1.8-fold improvement in L-GalOA and an average 2.1-fold improvement in GalAA production with expression of GAT-1-GFP over GFP alone (Figure 5D).
In the present study, we identified and characterized the first D-GalA transporter from a eukaryotic organism out of a group of candidate genes derived from a co-expression analysis performed in N. crassa. Genetic and biochemical analyses show that the N. crassa Δgat-1 deletion strain is substantially affected in growth on pectic substrates and is unable to take up D-GalA from the medium. Moreover, expression of a gat-1 construct in yeast conferred the ability for robust intracellular D-GalA accumulation rates, providing evidence for GAT-1 being a bona fide D-GalA transport protein. In view of the fact that the Δgat-1 deletion strain is also impaired in D-GalA-mediated signaling events, this furthermore corroborates the notion of an activation system for pectinolytic genes responding to D-GalA ( and references therein) and strongly suggests that GAT-1 is an integral part of the D-GalA signaling system in N. crassa - either by taking up D-GalA as the signaling molecule (or a precurser thereof) or by acting as a receptor itself.
Several bacterial D-GalA transport systems have already been described, such as ExuT, TogT, TogMNAB and KdgM from E. coli and Erwinia, a bacterial soft-rot causing plant pathogen. KdgM (a porin of the outer membrane ), TogMNAB (an ABC superfamily member in the inner membrane ) and TogT (an MFS-type transporter of the glycoside-pentoside-hexuronide (GPH) sub-family also located in the inner membrane ) are responsible for D-GalA oligosaccharide transport, whereas ExuT was shown to exclusively transport monomeric D-GalA . Functionally, therefore, GAT-1 seems most similar to ExuT. Indeed, both are high-affinity transporters (KM of 25 μM for ExuT versus 1.2 μM for GAT-1), display uptake rates in a very similar range (Vmax of 38 nmol min-1 mg protein-1 for ExuT versus 12.8 nmol min-1 mg protein-1 for GAT-1), transport D-GlcA in addition to D-GalA and are sensitive to the presence of uncouplers . However, it is noteworthy to mention that the kinetic parameters of GAT-1 were recorded as super-folder (sf)GFP fusion protein heterologously expressed in S. cerevisiae, and that therefore the behavior of the native protein in N. crassa might differ from these values. Despite belonging to the same MFS-type superfamily, ExuT (345 aa in E. chrysanthemi) and GAT-1 (537 aa in N. crassa) do not share any recognizable protein or sequence similarity and, therefore, have to be considered functional analogs.
Orthologs of GAT-1 (as well as the quinate permease) are found throughout species of Ascomycota and Basidiomycota (see Additional file 2: Figure S2), with the notable exception of Tuber melanosporum, which only seems to have a quinate permease and no additional homolog of GAT-1, and S. cerevisiae, which does not have a clear ortholog of either. Our phylogenetic analysis shows that GAT-1 homologs can be clustered into three clades, each with representatives from both Ascomycota and Basidiomycota (Additional file 2: Figure S2). It furthermore suggests that the GAT-1-containing clade 1 and the quinate permease clade both share a common ancestor and are equally derived from clade 2 and the more distantly related clade 3. While this manuscript was in proof, Enquist-Newman et al. (2014) discovered another uronic acid transporter (DHT1 from Asteromyces cruciatus) for the uptake of the alginate monomer 4-deoxy-L-erythro-5-hexoseulose uronate (DEHU), which represents a homolog to GAT-1 and the quinic acid permeases and is another good example of the fact that these two transporter families are evolutionarily well related .
Considering the number of GAT-1 homologs per genome, Botrytis cinerea seems most diverse, having GAT-1 homologs in all three clades, whereas others, such as Aspergillus nidulans and Trichoderma reesei, have one GAT-1 ortholog (in clade 1) and one or two paralogs in clade 2, respectively. Interestingly, Aspergillus niger only has a clade-2 homolog to GAT-1. However, since the corresponding gene (An14g04280) has been previously shown to be co-regulated with genes involved in pectin degradation and D-GalA catabolism and failed to express in absence of gaaB, the hypothesis of it being a D-GalA transporter has been brought forward [25, 45]. Additionally, Zhang et al. (2013) recently reported the An14g04280 homolog in Botrytis cinerea (BcHxt15) to be specifically induced by D-GalA and the corresponding mutant to be impaired in growth on D-GalA, which led them to suggest that BcHXT15 contributes to D-GalA uptake . Even though no representative has been properly biochemically characterized yet, considering these lines of evidence, it seems possible that the clade-2 homologs are also bona fide D-GalA transport proteins.
Regarding biotechnological applicability, we demonstrated a transporter-dependent conversion of D-GalA into two different molecules, L-GalOA and GalAA, by recombinantly expressing D-galacturonate reductase (GAAA) and uronate dehydrogenase (UDH) in yeast, respectively. GalAA is an isomer of glucaric acid, which was identified as a “top value-added chemical” in a report by the Pacific Northwest National Laboratory and the National Renewable Energy Laboratory for its potential use as a building block for hydroxylated nylons [47–50]. Consequently, GalAA is also likely to be capable of producing desirable hydroxylated nylons. L-GalOA on the other hand is reported to have similar properties to the commodity compound D-gluconic acid and correspondingly holds potential as a chelator in a variety of industries (for example cosmetic, pharmaceutical, food) and can also be used as a chemical precursor to L-ascorbic acid [45, 51, 52]. Production of both of these chemicals was completely GAT-1-dependent at the low D-GalA concentrations typical for high-affinity transport (90 μM). Even at high D-GalA concentrations (100 mM/2%), at which non-specific uronic sugar uptake has been shown to be substantial , expression of the transporter increased product accumulation roughly 2-fold over that observed in control strains at pH 6.0. Thus, D-GalA utilization can be improved over a broad concentration range in engineered yeast strains expressing a heterologous transporter. Interestingly, whereas production of GalAA was dependent on both GAT-1 and UDH, production of L-GalOA did not require heterologous expression of GAAA. Therefore, a basal level of native D-GalA reductase activity is present in S. cerevisiae, though its origin remains unknown.
Utilization of D-GalA for the biosynthesis of biofuels, specifically ethanol, from pectin-rich feedstocks has been a focus of increasing attention in recent years. The isolation and demonstration of activity of this transporter in S. cerevisiae provides an important contribution towards realization of this goal. Although functional expression of the requisite pathway enzymes remains an unsolved challenge, an increasing number of alternative orthologous enzymes and pathways are currently being explored [53–55], which will potentially allow this issue to be overcome in the near future.
In the present manuscript we have for the first time unambiguously identified and biochemically characterized a D-GalA transporter from a eukaryotic origin. GAT-1 was found to be the physiologically relevant D-GalA transporter in N. crassa, involved both in the uptake of D-GalA as a metabolite as well as for pectinase signaling events. Moreover, by heterologous expression of GAT-1 in combination with two downstream enzymes, we demonstrate successful transport and catalysis of D-GalA in the fermentation host S. cerevisiae, and thereby provide a considerable first step towards a more complete utilization of biomass for biofuel and value-added chemicals production.
Strains, media, and growth conditions
N. crassa WT (FGSC #2489) and gene deletion strains used in this study were obtained from the Fungal Genetics Stock Center (FGSC; [http://www.fgsc.net]). N. crassa was grown on 1 × Vogel’s salts  with either sucrose (Fisher S3-12; Fisher Scientific, Chicago, IL, USA), or pectin (Sigma P9135; Sigma-Aldrich, St Louis, MO, USA), xylan (beechwood, Sigma X4252), orange peel powder (OPP ), or Avicel (PH 101, Fluka 11365) at the indicated concentrations (w/v) at 25°C and 200 rpm with constant light unless stated otherwise.
The yeast strain used in this study was BY4742 (MATα his3Δ1 leu2Δ0 lys2Δ0 ura3Δ0) . Growth was usually performed in yeast extract peptone dextrose (YPD) media supplemented with 100 mg/L adenine hemisulfate and transformed strains were grown in the appropriate complete minimal dropout media, supplemented with 100 mg/L adenine hemisulfate. Yeast cultures were grown in 250-mL baffled flasks at 30°C and 250 rpm for microscopy and uptake experiments.
Plasmids and cloning
For yeast expression purposes, the 2 μ plasmid pRS426, equipped with the PGK1 promoter and sfGFP as generated and described in , was used. The gat-1 cDNA (with optimized Kozak sequence) was inserted between SpeI and ClaI sites using the primers ATACTAGT AAAAATGGGTCTTTCGATAGGAAATAGG and ATATCGAT AACATAAACCTCCACATGCTTCG (restriction sites in boldface). For genomic integration, the gat-1 sequence (RefSeq: XM_958805.2) with a C-terminal GFP fusion, or GFP alone, with an upstream TDH3 promoter was flanked by respective 5′ and 3′ ura3 homology regions using golden-gate cloning. The integration cassette was amplified by PCR and the resulting fragments were used in yeast transformations. The Agrobacterium tumefacians udh open reading frame [GenBank: BK006462.1] and a codon-optimized Aspergillus niger GAAA [GenBank: ABQ53587.1] were cloned into a Cen6/leu2 vector derived from pRS316  with a TDH3 promoter. For used gene sequences, see Table S1 in Additional file 4.
Phenotypic analyses of deletion strains were usually performed in 24 deep-well plates in a volume of 3 mL over the course of 4 days. Dry weight was determined after an overnight incubation of the mycelial mass in aluminum pans in a 105°C oven. Pectin consumption over time was followed using the phenol-sulfuric acid assay (PSA) . For this, pectin-grown N. crassa cultures were initially cleared by centrifugation at 20,000 × g for 5 minutes: 2.5 μL of the supernatant was diluted into 150 μL de-ionized (DI) water and mixed with 150 μL of a 5% phenol solution. Finally, 750 μL of concentrated sulfuric acid was quickly added, and the suspension was vortexed for 5 to 10 seconds. After about 7 minutes incubation at room temperature, 200 μL of the reactions were analyzed in a plate-reader spectrophotometer (Paradigm; Beckman-Coulter, Brea, CA, USA) at 487 nm.
Protein sequences of GAT-1, NCU06026 (quinate permease), and fungal homologs from representative ascomycete and basidiomycete species were retrieved from the NCBI database. Sequences were aligned and a phylogenetic tree constructed using the program Phylogeny.fr  (alignment by MUSCLE - excluding curation by Gblocks, phylogeny calculation by PhyML, and tree rendering by TreeDyn). Gal2p from S. cerevisiae was used as the out-group.
Media shift assays for transcriptional studies
The standard N. crassa growth conditions for any experiments involving media switches were as follows: cultures were pre-grown from 9- to 10-day-old conidia for 16 h in 3 mL of 1 × Vogel’s salts plus 2% (w/v) sucrose in 24 deep-well plates. The mycelia were then washed three times in 1 × Vogel’s salts without added carbon (NoC) and transferred to 1 × Vogel’s salts plus the respective new carbon source (as indicated). The mycelial mass was harvested after an additional 4 h by quickly blotting dry on Whatman tissues and subsequent flash-freezing in liquid nitrogen to be stored at -80°C. Total RNA from frozen samples was isolated using zirconia/silica beads (0.5 mm diameter; Biospec, Bartlesville, OK, USA) and a Mini-Beadbeater-96 (Biospec) with 1 mL TRIzol reagent (Invitrogen/Life Technologies, Carlsbad, CA, USA) according to the manufacturer’s instructions. The total RNA was further digested with TURBO DNA-free (Ambion/Life Technologies) and purified using an RNeasy kit (Qiagen, Valencia, CA, USA). RNA concentration and integrity was checked by Nanodrop and agarose gel electrophoresis.
Quantitative RT-PCR was performed using the EXPRESS One-Step SYBR GreenER with Premixed ROX kit (Invitrogen/Life Technologies) and the StepOnePlus Real-Time PCR System (Applied Biosystems/Life Technologies). Reactions were performed in triplicate (each from three biological replicates) with a total reaction volume of 10 μL including 300 nM each of forward and reverse primers (Additional file 5: Table S2) and 75 ng template RNA. Data analysis was performed by the StepOne Software (Applied Biosystems) using the Relative Quantitation/Comparative CT (ΔΔCT) setting. Data were normalized to the endogenous control actin (NCU04173) with expression on sucrose or no carbon as the reference sample (as indicated).
N. crassa cultures were pre-grown from 10-day-old conidia for 16 h in 100 mL of 1 × Vogel’s salts plus 2% (w/v) sucrose using 250-mL shake-flasks. The mycelia were then washed three times in 1 × Vogel’s salts without added carbon (NoC) and transferred to 2% sucrose or 1% (w/v) pectin for induction [37, 62]. After an additional 4 h, the mycelial mass was harvested over a Whatman glass microfiber filter (GF/F) on a Buchner funnel and subsequently flash-frozen in liquid nitrogen to be stored at -80°C. RNA was then prepared as above.
Library preparation for RNA sequencing (RNA-Seq) was performed essentially as described by Benz et al. and Coradetti et al. [37, 62]. These were sequenced in SR50 mode on an Illumina HiSeq2000 at the UC Davis Genome Center, and the files analyzed using the Illumina RTA 1.12 software. Mapping of the reads was done against the current version at the time of the N. crassa OR74A genome (v10)  using Tophat v1.2.0 (http://tophat.cbcb.umd.edu/) . Transcript abundance was estimated with Cufflinks v0.9.3 in fragments per kilobase of transcript per million mapped reads (FPKMs) using upper quartile normalization and mapping against reference isoforms from the Broad Institute (http://cufflinks.cbcb.umd.edu/) [64–66]. Profiling data are available at the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) [GEO: GSE42692]. Genes with a multiple-hypothesis-adjusted P-value of <0.05 using Cuffdiff were called as significantly differentially expressed between conditions. Independent triplicate cultures were harvested and analyzed for N. crassa WT on pectin and sucrose.
Monosaccharide uptake assays
Transport assays were performed according to . Briefly, N. crassa cultures were pre-grown from 9- to 10-day-old conidia for 16 h in 3 mL of 1 × Vogel’s salts plus 2% (w/v) sucrose using 24 deep-well plates. The mycelia were then washed three times in 1 × Vogel’s salts without added carbon (NoC) and transferred to 0.5% (w/v) pectin for induction. After an additional 4 h, the mycelia were washed again as above and transferred into the uptake buffer (1 × Vogel’s salts plus 90 μM monosaccharides or quinic acid, as indicated) for pre-equilibration and to reduce the dilution of the final uptake reaction. Two mycelia of the same genotype were combined into one well to increase the biomass at this stage. In the case of the addition of uncoupler, CCCP (in ethanol at the indicated concentrations; controls with ethanol only) was added at this stage and kept throughout the assay. After 5 minutes, the mycelia were transferred into fresh uptake buffer to start the reaction. Time points of the supernatants were usually taken at 0, 5, and up to 40 minutes. The samples were cleared by centrifugation (1 minute at 20,000 × g) and 50 μL of the supernatant diluted into 450 μL of DI water. The monosaccharide concentrations were then quantified by HPAEC-PAD (Thermo Fisher Scientific Inc., Bannockburn, IL, USA): for neutral sugars, a sample size of 25 μL was injected onto a Dionex CarboPac PA20 column (3 × 30 mm guard and 3 × 150 mm analytical) and eluted at 30°C using an isocratic mobile phase of 18 mM KOH at 0.4 mL/minute over 11 minutes. For uronic acids, 25 μL was injected onto a Dionex CarboPac PA200 column (3 × 30 mm guard and 3 × 250 mm analytical) and eluted at 30°C using a 50 mM to 170 mM sodium acetate gradient (in 0.1 M NaOH) over 8 minutes at 0.4 mL/minutes.
For mass spectrometry (MS) detection of quinic acid, samples were injected onto a Rezex RFQ Fast Fruit H+ (8%) (Phenomenex, Torrance, CA, USA) column (100 × 7.8 mm) and eluted at 55°C using an isocratic mobile phase of 0.5% formic acid at a flow rate of 0.3 mL/minute. An LTQ XL linear ion trap (Thermo Scientific) was used for MS and MS/MS detection in negative ion mode. The quinic acid precursor ion of m/z 191.2 was followed (isolation width m/z 2; CID setting 35; product scan range m/z 50 to 200), and for chromatogram processing the mass transition (MS/MS) of m/z 191.2 to m/z 127.2 was selected.
Yeast bioconversion assays
Yeast strains were grown as described above and harvested at mid-log phase. Cells were washed twice in ½ × Vogel’s salts (pH 5.8) with 50 mM ethanol added as a non-fermentable carbon source and resuspended in a final volume to yield an optical density (OD)600 of 20. To each experimental condition, 0.5 mL of concentrated cells were added to 24-well deep-well blocks and 0.5 mL of 2 × solution was added to yield a final concentration of ½ × Vogel’s salts, 90 μM D-GalA (Fluka/Sigma-Aldrich 48280-F) and 50 mM ethanol at an OD600 of 10. Blocks were incubated at 30°C with 750 rpm of shaking until the appropriate time point. Samples were cleared by centrifugation and the supernatant was prepared and analyzed as described above. Cell pellets were washed 2 × in cold DI water and freeze-dried in pre-weighed tubes to calculate dry cell weight. For assaying intracellular metabolites, desiccated yeast pellets were subjected to bead-beating (1 minute) in a Mini-Beadbeater-96 (Biospec) with an added 2:2:1 (v/v/v) ratio of chloroform, methanol and water. The aqueous phase was collected after centrifugation at 14,000 rpm for 10 minutes and was subsequently concentrated into 60 μL of DI water after vacuum centrifugation. The resulting sample was used for MS and MS/MS detection as described above. The respective precursor ions and mass transitions followed were: D-GalA: m/z 193.1 to m/z 131.1, L-GalOA: m/z 195.1 to m/z 159.1; and GalAA: m/z 209.1 to m/z 85.1. For the calculation of intracellular metabolite concentrations, dry cell weight was converted to cell volume using the factor 2.2 μL/mg dry weight [31, 67].
High-D-GalA concentration yeast bioconversions were performed in a similar manner as low-concentration uptake assays, but with some modifications. Cells were grown and harvested at mid-log and washed twice in 1 × PBS, pH 6.0, with an added 50 mM EtOH. The final cell suspensions in this buffer were adjusted to an OD600 of 20 and 0.5 mL of cells were added to each condition in a 24 deep-well block in triplicate. To this, 0.5 mL of solution was added to yield 100 mM D-GalA in 1 × PBS (pH 6.0) with an OD600 of 10. Blocks were incubated at 30°C with 750 rpm of shaking for 1 h. To assay for bioconversion products, a portion of the cell suspension was flash-frozen in liquid nitrogen and added as a 2:2:1 (v/v/v) ratio of chloroform, methanol and cell suspension. Washing the cell pellets was avoided to minimize product loss . A duplicate portion of cell suspension of each condition was lyophilized to calculate dry cell weight. Aqueous phase extraction and product detection proceeded as described above.
(3H)-Galacturonic acid transport assays in yeast
Transport assays were performed using a modified oil-stop method . In short, yeast strains expressing gat-1 fused to GFP were grown to an OD (600 nm) of 4.0 to 5.0 in selective media, washed 3 × with ice cold assay buffer (30 mM MES-NaOH (pH 5.6) and 50 mM ethanol), and resuspended to an OD of 40.0. To start transport reactions, 50 μL of cells were added to 50 μL of (3H)-D-GalA layered over 100 μL of silicone oil (Sigma 85419). Reactions were stopped by spinning cells through oil for 1 minute at 17,000 g, tubes were frozen in ethanol/dry ice, and tube-bottoms containing the cell pellets were clipped off into 1 mL of 0.5 M NaOH. The pellets were solubilized overnight, 5 mL of Ultima Gold scintillation fluid added, and counts per minute (CPM) determined in a Tri-Carb 2900TR scintillation counter. (3H)-D-GalA was purchased from ViTrax, Inc. (Placentia, CA, USA), and had a specific activity of 25 Ci/mmol and a purity of >99%. Kinetic parameters were determined by measuring the linear rate of (3H)-D-GalA uptake over 3 minutes for D-GalA concentrations between 0.1 and 10 μM. Vmax and KM values were determined by fitting a “Ligand binding”, “one site saturation” function to a plot of rates versus D-GalA concentrations by non-linear regression in SigmaPlot®. Vmax values were normalized by total protein content using a Bradford assay. Kinetic parameters reported in the text are the mean ± SD from three separate experiments.
meso-galactaric acid (mucic acid)
green fluorescent protein
high pH anion-exchange chromatography with pulsed amperometric detection
liquid chromatography coupled to tandem mass-spectrometry
major facilitator superfamily
Nicotinamide adenine dinucleotide
Nicotinamide adenine dinucleotide phosphate (reduced)
phosphoglycerate kinase 1
relative transcript quantities
reverse transcriptase PCR
super-folder green fluorescent protein
yeast extract peptone dextrose.
We want to thank Ana B Ibáñez for excellent technical assistance, Jon Galazka, Yuping Lin and Ligia Acosta-Sampson for help with the set-up of the tritiated kinetics experiments, and Shu-Lun Tang for help with the ChemDraw structures. We thank Luke Latimer for helpful advice and for his comments on the manuscript. We furthermore acknowledge the FGSC and the Neurospora Genome Project for continuous support. This work was supported by grants from the Energy Biosciences Institute to CRS and JED. JPB was supported by a Feodor-Lynen fellowship from the Humboldt-Foundation (GER).
- Doran JB, Cripe J, Sutton M, Foster B: Fermentations of pectin-rich biomass with recombinant bacteria to produce fuel ethanol. Appl Biochem Biotechnol 2000, 84–86: 141-152.View ArticleGoogle Scholar
- Kennedy M, List D, Lu Y, Foo LY, Newman RH, Sims IM, Bain PJS, Hamilton B, Fenton G: Apple Pomace and Products Derived from Apple Pomace: Uses, Composition and Analysis. In Analysis of Plant Waste Materials. Edited by: Linskens HF, Jackson JF. Berlin, Heidelberg: Springer; 1999:75-119.View ArticleGoogle Scholar
- Rivas B, Torrado A, Torre P, Converti A, Dominguez JM: Submerged citric acid fermentation on orange peel autohydrolysate. J Agric Food Chem 2008, 56: 2380-2387. 10.1021/jf073388rView ArticleGoogle Scholar
- Angel Siles Lopez J, Li Q, Thompson IP: Biorefinery of waste orange peel. Crit Rev Biotechnol 2010, 30: 63-69. 10.3109/07388550903425201View ArticleGoogle Scholar
- Grohmann K, Baldwin EA, Buslig BS: Production of ethanol from enzymatically hydrolyzed orange peel by the yeast Saccharomyces cerevisiae . Appl Biochem Biotechnol 1994, 45–46: 315-327.View ArticleGoogle Scholar
- Hang YD, Lee CY, Woodams EE, Cooley HJ: Production of alcohol from apple pomace. Appl Environ Microbiol 1981, 42: 1128-1129.Google Scholar
- Rorick R, Nahar N, Pryor SW: Ethanol production from sugar beet pulp using escherichia coli KO11 and Saccharomyces cerevisiae . Biol Eng 2011, 3: 199-209.View ArticleGoogle Scholar
- Vendruscolo F, Albuquerque PM, Streit F, Esposito E, Ninow JL: Apple pomace: a versatile substrate for biotechnological applications. Crit Rev Biotechnol 2008, 28: 1-12. 10.1080/07388550801913840View ArticleGoogle Scholar
- Edwards MC, Doran-Peterson J: Pectin-rich biomass as feedstock for fuel ethanol production. Appl Microbiol Biotechnol 2012, 95: 565-575. 10.1007/s00253-012-4173-2View ArticleGoogle Scholar
- Xiao C, Anderson CT: Roles of pectin in biomass yield and processing for biofuels. Front Plant Sci 2013, 4: 67.View ArticleGoogle Scholar
- Atmodjo MA, Hao Z, Mohnen D: Evolving views of pectin biosynthesis. Annu Rev Plant Biol 2013, 64: 747-779. 10.1146/annurev-arplant-042811-105534View ArticleGoogle Scholar
- Mohnen D: Pectin structure and biosynthesis. Curr Opin Plant Biol 2008, 11: 266-277. 10.1016/j.pbi.2008.03.006View ArticleGoogle Scholar
- Amorim HV, Lopes ML, de Castro Oliveira JV, Buckeridge MS, Goldman GH: Scientific challenges of bioethanol production in Brazil. Appl Microbiol Biotechnol 2011, 91: 1267-1275. 10.1007/s00253-011-3437-6View ArticleGoogle Scholar
- van Maris AJ, Abbott DA, Bellissimi E, van den Brink J, Kuyper M, Luttik MA, Wisselink HW, Scheffers WA, van Dijken JP, Pronk JT: Alcoholic fermentation of carbon sources in biomass hydrolysates by Saccharomyces cerevisiae : current status. Antonie Van Leeuwenhoek 2006, 90: 391-418. 10.1007/s10482-006-9085-7View ArticleGoogle Scholar
- Barnett JA, Payne RW, Yarrow D (Eds): A guide to identifying and classifying yeasts. Cambridge: Cambridge University Press; 1990.Google Scholar
- Huisjes EH, Luttik MA, Almering MJ, Bisschops MM, Dang DH, Kleerebezem M, Siezen R, van Maris AJ, Pronk JT: Toward pectin fermentation by Saccharomyces cerevisiae : expression of the first two steps of a bacterial pathway for D-galacturonate metabolism. J Biotechnol 2012, 162: 303-310. 10.1016/j.jbiotec.2012.10.003View ArticleGoogle Scholar
- Nijkamp JF, van den Broek M, Datema E, de Kok S, Bosman L, Luttik MA, Daran-Lapujade P, Vongsangnak W, Nielsen J, Heijne WHM, Klaassen P, Paddon CJ, Platt D, Kotter P, van Ham RC, Reinders MJT, Pronk JT, de Ridder D, Daran JM: De novo sequencing, assembly and analysis of the genome of the laboratory strain Saccharomyces cerevisiae CEN.PK113-7D, a model for modern industrial biotechnology. Microb Cell Fact 2012, 11: 36. 10.1186/1475-2859-11-36View ArticleGoogle Scholar
- Huisjes EH, de Hulster E, van Dam JC, Pronk JT, van Maris AJ: Galacturonic acid inhibits the growth of Saccharomyces cerevisiae on galactose, xylose, and arabinose. Appl Environ Microbiol 2012, 78: 5052-5059. 10.1128/AEM.07617-11View ArticleGoogle Scholar
- Ashwell A, Wahba AJ, Hickman J: A new pathway of uronic acid metabolism. Biochim Biophys Acta 1958, 30: 186-187. 10.1016/0006-3002(58)90257-9View ArticleGoogle Scholar
- Ashwell G: Enzymes of glucuronic and galacturonic acid metabolism in bacteria. Method Enzymol 1962, 5: 190-208.View ArticleGoogle Scholar
- Hugouvieux-Cotte-Pattat N, Robert-Baudouy J: Hexuronate catabolism in Erwinia chrysanthemi . J Bacteriol 1987, 169: 1223-1231.Google Scholar
- Kilgore WW, Starr MP: Catabolism of galacturonic and glucuronic acids by Erwinia carotovora . J Biol Chem 1959, 234: 2227-2235.Google Scholar
- Hilditch S: Identification of the fungal catabolic D-galacturonate pathway. PhD thesis. University of Helsinki, Faculty of Biosciences, Department of Biological and Environmental Sciences, Division of Biochemistry; 2010.Google Scholar
- Hilditch S, Berghall S, Kalkkinen N, Penttila M, Richard P: The missing link in the fungal D-galacturonate pathway: identification of the L-threo-3-deoxy-hexulosonate aldolase. J Biol Chem 2007, 282: 26195-26201. 10.1074/jbc.M704401200View ArticleGoogle Scholar
- Martens-Uzunova E: Assessment of the Pectinolytic Network of Aspergillus niger by Functional Genomics - Insights from the Transcriptome. PhD thesis. Wageningen University, Laboratory of Microbiology; 2008.Google Scholar
- Martens-Uzunova ES, Schaap PJ: An evolutionary conserved d-galacturonic acid metabolic pathway operates across filamentous fungi capable of pectin degradation. Fungal Genet Biol 2008, 45: 1449-1457. 10.1016/j.fgb.2008.08.002View ArticleGoogle Scholar
- Zhang L: Pectin degradation by Botrytis cinerea : recognition of endo-polygalacturonases by an Arabidopsis receptor and utilization of D-galacturonic acid. PhD thesis. Wageningen University, Laboratory of Phytopathology; 2013.Google Scholar
- Zhang L, Thiewes H, van Kan JA: The D-galacturonic acid catabolic pathway in Botrytis cinerea . Fungal Genet Biol 2011, 48: 990-997. 10.1016/j.fgb.2011.06.002View ArticleGoogle Scholar
- Richard P, Hilditch S: D-galacturonic acid catabolism in microorganisms and its biotechnological relevance. Appl Microbiol Biotechnol 2009, 82: 597-604. 10.1007/s00253-009-1870-6View ArticleGoogle Scholar
- Huisjes E: Towards fermentation of galacturonic acid-containing feedstocks with Saccharomyces cerevisiae . PhD thesis. Delft University of Technology, Industrial Microbiology Section, Department of Biotechnology; 2013.Google Scholar
- Souffriau B, den Abt T, Thevelein JM: Evidence for rapid uptake of D-galacturonic acid in the yeast Saccharomyces cerevisiae by a channel-type transport system. FEBS Lett 2012, 586: 2494-2499. 10.1016/j.febslet.2012.06.012View ArticleGoogle Scholar
- Blot N, Berrier C, Hugouvieux-Cotte-Pattat N, Ghazi A, Condemine G: The oligogalacturonate-specific porin KdgM of Erwinia chrysanthemi belongs to a new porin family. J Biol Chem 2002, 277: 7936-7944. 10.1074/jbc.M109193200View ArticleGoogle Scholar
- Hugouvieux-Cotte-Pattat N, Blot N, Reverchon S: Identification of TogMNAB, an ABC transporter which mediates the uptake of pectic oligomers in Erwinia chrysanthemi 3937. Mol Microbiol 2001, 41: 1113-1123.View ArticleGoogle Scholar
- Hugouvieux-Cotte-Pattat N, Reverchon S: Two transporters, TogT and TogMNAB, are responsible for oligogalacturonide uptake in Erwinia chrysanthemi 3937. Mol Microbiol 2001, 41: 1125-1132.View ArticleGoogle Scholar
- San Francisco MJ, Keenan RW: Uptake of galacturonic acid in Erwinia chrysanthemi EC16. J Bacteriol 1993, 175: 4263-4265.Google Scholar
- Znameroski EA, Glass NL: Using a model filamentous fungus to unravel mechanisms of lignocellulose deconstruction. Biotechnol Biofuels 2013, 6: 6. 10.1186/1754-6834-6-6View ArticleGoogle Scholar
- Benz JP, Chau BH, Zheng D, Bauer S, Glass NL, Somerville CR: A comparative systems analysis of polysaccharide-elicited responses in Neurospora crassa reveals carbon source-specific cellular adaptations. Mol Microbiol 2014, 91: 275-99. 10.1111/mmi.12459View ArticleGoogle Scholar
- Delmas S, Pullan ST, Gaddipati S, Kokolski M, Malla S, Blythe MJ, Ibbett R, Campbell M, Liddell S, Aboobaker A, Tucker GA, Archer DB: Uncovering the genome-wide transcriptional responses of the filamentous fungus Aspergillus niger to lignocellulose using RNA sequencing. Plos Genet 2012, 8: e1002875. 10.1371/journal.pgen.1002875View ArticleGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, Rehman B, Elkins T, Engels R, Wang S, Nielsen CB, Butler J, Endrizzi M, Qui D, Ianakiev P, Bell-Pedersen D, Nelson MA, Werner-Washburne M, Selitrennikoff CP, Kinsey JA, Braun EL, Zelter A, Schulte U, Kothe GO, Jedd G, Mewes W, et al.: The genome sequence of the filamentous fungus Neurospora crassa . Nature 2003, 422: 859-868. 10.1038/nature01554View ArticleGoogle Scholar
- Tang X, Dong W, Griffith J, Nilsen R, Matthes A, Cheng KB, Reeves J, Schuttler HB, Case ME, Arnold J, Logan DA: Systems biology of the qa gene cluster in Neurospora crassa . Plos One 2011, 6: e20671. 10.1371/journal.pone.0020671View ArticleGoogle Scholar
- Whittington HA, Grant S, Roberts CF, Lamb H, Hawkins AR: Identification and isolation of a putative permease gene in the quinic acid utilization (QUT) gene cluster of Aspergillus nidulans . Curr Genet 1987, 12: 135-139. 10.1007/BF00434668View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215: 403-410.View ArticleGoogle Scholar
- Yoon SH, Moon TS, Iranpour P, Lanza AM, Prather KJ: Cloning and characterization of uronate dehydrogenases from Two pseudomonads and Agrobacterium tumefaciens strain C58. J Bacteriol 2009, 191: 1565-1573. 10.1128/JB.00586-08View ArticleGoogle Scholar
- Enquist-Newman M, Faust AM, Bravo DD, Santos CN, Raisner RM, Hanel A, Sarvabhowman P, Le C, Regitsky DD, Cooper SR, Peereboom L, Clark A, Martinez Y, Goldsmith J, Cho MY, Donohoue PD, Luo L, Lamberson B, Tamrakar P, Kim EJ, Villari JL, Gill A, Tripathi SA, Karamchedu P, Paredes CJ, Rajgarhia V, Kotlar HK, Bailey RB, Miller DJ, Ohler NL, et al.: Efficient ethanol production from brown macroalgae sugars by a synthetic yeast platform. Nature 2014, 505: 239-243.View ArticleGoogle Scholar
- Kuivanen J, Mojzita D, Wang Y, Hilditch S, Penttila M, Richard P, Wiebe MG: Engineering filamentous fungi for conversion of D-galacturonic acid to L-galactonic acid. Appl Environ Microbiol 2012, 78: 8676-8683. 10.1128/AEM.02171-12View ArticleGoogle Scholar
- Zhang L, Hua C, Stassen JHM, Chatterjee S, Cornelissen M, van Kan JAL: Genome-wide analysis of pectate-induced gene expression in Botrytis cinerea : identification and functional analysis of putative D-galacturonate transporters. Fungal Genet Biol 2013. doi:10.1016/j.fgb.2013.10.002Google Scholar
- Kiely DE, Chen L, Lin TH: Synthetic polyhydroxypolyamides from galactaric, xylaric, D-glucaric, and D-mannaric acids and alkylenediamine monomers - Some comparisons. J Polym Sci Pol Chem 2000, 38: 594-603. 10.1002/(SICI)1099-0518(20000201)38:3<594::AID-POLA24>3.0.CO;2-#View ArticleGoogle Scholar
- Mojzita D, Wiebe M, Hilditch S, Boer H, Penttila M, Richard P: Metabolic engineering of fungal strains for conversion of D-galacturonate to meso -galactarate. Appl Environ Microbiol 2010, 76: 169-175. 10.1128/AEM.02273-09View ArticleGoogle Scholar
- Moon TS, Yoon SH, Lanza AM, Roy-Mayhew JD, Prather KLJ: Production of glucaric acid from a synthetic pathway in recombinant Escherichia coli . Appl Environ Microb 2009, 75: 589-595. 10.1128/AEM.00973-08View ArticleGoogle Scholar
- Werpy T, Petersen G, Aden A, Bozell J, Holladay J, White J, Manheim A, Eliot D, Lasure L, Jones S: Top Value Added Chemicals From Biomass. Volume 1 - Results of Screening for Potential Candidates From Sugars and Synthesis Gas. Department of Energy Report; 2004:1-76.Google Scholar
- Csiba M, Cleophax J, Petit S, Gero SD: An expedient and practical 3-step synthesis of vitamin-C from a by-product of the sugar-industry - the L-galactono-1,4-lactone pathway. J Org Chem 1993, 58: 7281-7282. 10.1021/jo00077a062View ArticleGoogle Scholar
- Roland JF, Cayle T, Dinwoodie RC, Mehnert DW: Fermentation production of ascorbic acid from L-galactonic substrate. 1983.Google Scholar
- Agius F, Gonzalez-Lamothe R, Caballero JL, Munoz-Blanco J, Botella MA, Valpuesta V: Engineering increased vitamin C levels in plants by overexpression of a D-galacturonic acid reductase. Nat Biotechnol 2003, 21: 177-181. 10.1038/nbt777View ArticleGoogle Scholar
- Chang YF, Feingold DS: D-Glucaric acid and galactaric acid catabolism by Agrobacterium tumefaciens . J Bacteriol 1970, 102: 85.Google Scholar
- Hamada S, Seike Y, Tanimori S, Sakamoto T, Kishida M: Characterization of D-galacturonate reductase purified from the psychrophilic yeast species Cryptococcus diffluens . J Biosci Bioeng 2011, 111: 518-521. 10.1016/j.jbiosc.2010.12.019View ArticleGoogle Scholar
- Vogel HJ: A convenient growth medium for Neurospora (medium N). Microbial Genetics Bulletin 1956, 13: 42.Google Scholar
- Brachmann CB, Davies A, Cost GJ, Caputo E, Li J, Hieter P, Boeke JD: Designer deletion strains derived from Saccharomyces cerevisiae S288C: a useful set of strains and plasmids for PCR-mediated gene disruption and other applications. Yeast 1998, 14: 115-132. 10.1002/(SICI)1097-0061(19980130)14:2<115::AID-YEA204>3.0.CO;2-2View ArticleGoogle Scholar
- Galazka JM, Tian C, Beeson WT, Martinez B, Glass NL, Cate JH: Cellodextrin transport in yeast for improved biofuel production. Science 2010, 330: 84-86. 10.1126/science.1192838View ArticleGoogle Scholar
- Lee ME, Aswani A, Han AS, Tomlin CJ, Dueber JE: Expression-level optimization of a multi-enzyme pathway in the absence of a high-throughput assay. Nucleic Acids Res 2013, 41: 10668-10678. 10.1093/nar/gkt809View ArticleGoogle Scholar
- Dubois M, Gilles KA, Hamilton JK, Rebers PA, Smith F: Colorimetric method for determination of sugars and related substances. Anal Chem 1956, 28: 350-356. 10.1021/ac60111a017View ArticleGoogle Scholar
- Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, Claverie JM, Gascuel O: Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res 2008, 36: W465-W469. 10.1093/nar/gkn180View ArticleGoogle Scholar
- Coradetti ST, Craig JP, Xiong Y, Shock T, Tian C, Glass NL: Conserved and essential transcription factors for cellulase gene expression in ascomycete fungi. Proc Natl Acad Sci USA 2012, 109: 7397-7402. 10.1073/pnas.1200785109View ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 2009, 10: R25. 10.1186/gb-2009-10-3-r25View ArticleGoogle Scholar
- Roberts A, Pimentel H, Trapnell C, Pachter L: Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics 2011, 27: 2325-2329. 10.1093/bioinformatics/btr355View ArticleGoogle Scholar
- Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L: Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol 2011, 12: R22. 10.1186/gb-2011-12-3-r22View ArticleGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 2010, 28: 511-515. 10.1038/nbt.1621View ArticleGoogle Scholar
- Sant’Ana GD, Paes LD, Paiva AFV, Fietto LG, Totola AH, Tropia MJM, Silveira-Lemos D, Lucas C, Fietto JLR, Brandao RL, Castro ID: Protective effect of ions against cell death induced by acid stress in Saccharomyces. Fems Yeast Res 2009, 9: 701-712. 10.1111/j.1567-1364.2009.00523.xView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.