Genome-scale model development and genomic sequencing of the oleaginous clade Lipomyces

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

Simple and Effective Squash-PCR for Rapid Genotyping of Industrial Microalgae

Microalgae are recognized for their versatility in providing renewable energy, biopharmaceuticals, and nutraceuticals, attributed to their sustainable, renewable, and cost-effective nature. Genetic engineering has proven highly effective in enhancing microalgae production. PCR-based genotyping is the primary method for screening genetically transformed microalgae cells. Recently, we developed a novel PCR method, namely Squash-PCR, and employed it for the molecular analysis of industrially important fungi and yeasts. In this study, we successfully implemented the Squash-PCR technique in 12 industrially significant algae species. This approach offers a quick and reliable means of obtaining DNA templates directly from squashed algal cells, eliminating the need for time-consuming and labor-intensive cultivation and genomic DNA extraction steps. Our results demonstrate the effectiveness of Squash-PCR in detecting and characterizing target genes of interest in 12 different algae species. Overall, this study establishes the Squash-PCR method as a valuable tool for molecular studies in algae, enabling researchers to rapidly screen and manipulate genetic traits in diverse algal species.

Advances in genome-scale metabolic models of industrially important fungi

Many fungal species have been used industrially for production of biofuels and bioproducts. Developing strains with better performance in biomanufacturing contexts requires a systematic understanding of cellular metabolism. Genome-scale metabolic models (GEMs) offer a comprehensive view of interconnected pathways and a mathematical framework for downstream analysis. Recently, GEMs have been developed or updated for several industrially important fungi. Some of them incorporate enzyme constraints, enabling improved predictions of cell states and proteome allocation. Here, we provide an overview of these newly developed GEMs and computational methods that facilitate construction of enzyme-constrained GEMs and utilize flux predictions from GEMs. Furthermore, we highlight the pivotal roles of these GEMs in iterative design-build-test-learn cycles, ultimately advancing the field of fungal biomanufacturing.

Identification of a specific exporter that enables high production of aconitic acid in Aspergillus pseudoterreus

Aconitic acid is an unsaturated tricarboxylic acid that is attractive for its potential use in manufacturing biodegradable and biocompatible polymers, plasticizers, and surfactants. Previously Aspergillus pseudoterreus was engineered as a platform to produce aconitic acid by deleting the cadA (cis-aconitic acid decarboxylase) gene in the itaconic acid biosynthetic pathway. In this study, the aconitic acid transporter gene (aexA) was identified using comparative global discovery proteomics analysis between the wild-type and cadA deletion strains. The protein AexA belongs to the Major Facilitator Superfamily (MFS). Deletion of aexA almost abolished aconitic acid secretion, while its overexpression led to a significant increase in aconitic acid production. Transportation of aconitic acid across the plasma membrane is a key limiting step in its production. In vitro, proteoliposome transport assay further validated AexA's function and substrate specificity. This research provides new approaches to efficiently pinpoint and characterize exporters of fungal organic acids and accelerate metabolic engineering to improve secretion capability and lower the cost of bioproduction.

Rapid and robust squashed spore/colony PCR of industrially important fungi

BACKGROUND

Fungi have been utilized for centuries in medical, agricultural, and industrial applications. Development of systems biology techniques has enabled the design and metabolic engineering of these fungi to produce novel fuels, chemicals, and enzymes from renewable feedstocks. Many genetic tools have been developed for manipulating the genome and creating mutants rapidly. However, screening and confirmation of transformants remain an inefficient step within the design, build, test, and learn cycle in many industrial fungi because extracting fungal genomic DNA is laborious, time-consuming, and involves toxic chemicals.

RESULTS

In this study we developed a rapid and robust technique called "Squash-PCR" to break open the spores and release fungal genomic DNA as a template for PCR. The efficacy of Squash-PCR was investigated in eleven different filamentous fungal strains. Clean PCR products with high yields were achieved in all tested fungi. Spore age and type of DNA polymerase did not affect the efficiency of Squash-PCR. However, spore concentration was found to be the crucial factor for Squash-PCR in Aspergillus niger, with the dilution of starting material often resulting in higher PCR product yield. We then further evaluated the applicability of the squashing procedure for nine different yeast strains. We found that Squash-PCR can be used to improve the quality and yield of colony PCR in comparison to direct colony PCR in the tested yeast strains.

CONCLUSION

The developed technique will enhance the efficiency of screening transformants and accelerate genetic engineering in filamentous fungi and yeast.

Engineering Rhodosporidium toruloides for production of 3-hydroxypropionic acid from lignocellulosic hydrolysate

Microbial production of valuable bioproducts is a promising route towards green and sustainable manufacturing. The oleaginous yeast, Rhodosporidium toruloides, has emerged as an attractive host for the production of biofuels and bioproducts from lignocellulosic hydrolysates. 3-hydroxypropionic acid (3HP) is an attractive platform molecule that can be used to produce a wide range of commodity chemicals. This study focuses on establishing and optimizing the production of 3HP in R. toruloides. As R. toruloides naturally has a high metabolic flux towards malonyl-CoA, we exploited this pathway to produce 3HP. Upon finding the yeast capable of catabolizing 3HP, we then implemented functional genomics and metabolomic analysis to identify the catabolic pathways. Deletion of a putative malonate semialdehyde dehydrogenase gene encoding an oxidative 3HP pathway was found to significantly reduce 3HP degradation. We further explored monocarboxylate transporters to promote 3HP transport and identified a novel 3HP transporter in Aspergillus pseudoterreus by RNA-seq and proteomics. Combining these engineering efforts with media optimization in a fed-batch fermentation resulted in 45.4 g/L 3HP production. This represents one of the highest 3HP titers reported in yeast from lignocellulosic feedstocks. This work establishes R. toruloides as a host for 3HP production from lignocellulosic hydrolysate at high titers, and paves the way for further strain and process optimization towards enabling industrial production of 3HP in the future.

PeakDecoder enables machine learning-based metabolite annotation and accurate profiling in multidimensional mass spectrometry measurements

Multidimensional measurements using state-of-the-art separations and mass spectrometry provide advantages in untargeted metabolomics analyses for studying biological and environmental bio-chemical processes. However, the lack of rapid analytical methods and robust algorithms for these heterogeneous data has limited its application. Here, we develop and evaluate a sensitive and high-throughput analytical and computational workflow to enable accurate metabolite profiling. Our workflow combines liquid chromatography, ion mobility spectrometry and data-independent acquisition mass spectrometry with PeakDecoder, a machine learning-based algorithm that learns to distinguish true co-elution and co-mobility from raw data and calculates metabolite identification error rates. We apply PeakDecoder for metabolite profiling of various engineered strains of Aspergillus pseudoterreus, Aspergillus niger, Pseudomonas putida and Rhodosporidium toruloides. Results, validated manually and against selected reaction monitoring and gas-chromatography platforms, show that 2683 features could be confidently annotated and quantified across 116 microbial sample runs using a library built from 64 standards.

Metabolic engineering to improve production of 3-hydroxypropionic acid from corn-stover hydrolysate in Aspergillus species

BACKGROUND

Fuels and chemicals derived from non-fossil sources are needed to lessen human impacts on the environment while providing a healthy and growing economy. 3-hydroxypropionic acid (3-HP) is an important chemical building block that can be used for many products. Biosynthesis of 3-HP is possible; however, low production is typically observed in those natural systems. Biosynthetic pathways have been designed to produce 3-HP from a variety of feedstocks in different microorganisms.

RESULTS

In this study, the 3-HP β-alanine pathway consisting of aspartate decarboxylase, β-alanine-pyruvate aminotransferase, and 3-hydroxypropionate dehydrogenase from selected microorganisms were codon optimized for Aspergillus species and placed under the control of constitutive promoters. The pathway was introduced into Aspergillus pseudoterreus and subsequently into Aspergillus niger, and 3-HP production was assessed in both hosts. A. niger produced higher initial 3-HP yields and fewer co-product contaminants and was selected as a suitable host for further engineering. Proteomic and metabolomic analysis of both Aspergillus species during 3-HP production identified genetic targets for improvement of flux toward 3-HP including pyruvate carboxylase, aspartate aminotransferase, malonate semialdehyde dehydrogenase, succinate semialdehyde dehydrogenase, oxaloacetate hydrolase, and a 3-HP transporter. Overexpression of pyruvate carboxylase improved yield in shake-flasks from 0.09 to 0.12 C-mol 3-HP C-mol^-1 glucose in the base strain expressing 12 copies of the β-alanine pathway. Deletion or overexpression of individual target genes in the pyruvate carboxylase overexpression strain improved yield to 0.22 C-mol 3-HP C-mol^-1 glucose after deletion of the major malonate semialdehyde dehydrogenase. Further incorporation of additional β-alanine pathway genes and optimization of culture conditions (sugars, temperature, nitrogen, phosphate, trace elements) for 3-HP production from deacetylated and mechanically refined corn stover hydrolysate improved yield to 0.48 C-mol 3-HP C-mol^-1 sugars and resulted in a final titer of 36.0 g/L 3-HP.

CONCLUSIONS

The results of this study establish A. niger as a host for 3-HP production from a lignocellulosic feedstock in acidic conditions and demonstrates that 3-HP titer and yield can be improved by a broad metabolic engineering strategy involving identification and modification of genes participated in the synthesis of 3-HP and its precursors, degradation of intermediates, and transport of 3-HP across the plasma membrane.

Bayesian Inference for Integrating Yarrowia lipolytica Multiomics Datasets with Metabolic Modeling

Optimizing the metabolism of microbial cell factories for yields and titers is a critical step for economically viable production of bioproducts and biofuels. In this process, tuning the expression of individual enzymes to obtain the desired pathway flux is a challenging step, in which data from separate multiomics techniques must be integrated with existing biological knowledge to determine where changes should be made. Following a design-build-test-learn strategy, building on recent advances in Bayesian metabolic control analysis, we identify key enzymes in the oleaginous yeast Yarrowia lipolytica that correlate with the production of itaconate by integrating a metabolic model with multiomics measurements. To this extent, we quantify the uncertainty for a variety of key parameters, known as flux control coefficients (FCCs), needed to improve the bioproduction of target metabolites and statistically obtain key correlations between the measured enzymes and boundary flux. Based on the top five significant FCCs and five correlated enzymes, our results show phosphoglycerate mutase, acetyl-CoA synthetase (ACSm), carbonic anhydrase (HCO3E), pyrophosphatase (PPAm), and homoserine dehydrogenase (HSDxi) enzymes in rate-limiting reactions that can lead to increased itaconic acid production.

Genetically Engineered Oleaginous Yeast Lipomyces starkeyi for Sesquiterpene α-Zingiberene Production

Oleaginous yeast, such as Lipomyces starkeyi, are logical organisms for production of higher energy density molecules like lipids and terpenes. We demonstrate that transgenic L. starkeyi strains expressing an α-zingiberene synthase gene from lemon basil or Hall's panicgrass can produce up to 17 mg/L α-zingiberene in yeast extract peptone dextrose (YPD) medium containing 4% glucose. The transgenic strain was further examined in 8% glucose media with C/N ratios of 20 or 100, and YPD. YPD medium resulted in 59 mg/L α-zingiberene accumulation. Overexpression of selected genes from the mevalonate pathway achieved 145% improvement in α-zingiberene synthesis. Optimization of the growth medium for α-zingiberene production led to 15% higher titer than YPD medium. The final transgenic strain produced 700 mg/L α-zingiberene in fed-batch bioreactor culture. This study opens a new synthetic route to produce α-zingiberene or other terpenoids in L. starkeyi and establishes this yeast as a platform for jet fuel biosynthesis.

Integration of Proteomics and Metabolomics Into the Design, Build, Test, Learn Cycle to Improve 3-Hydroxypropionic Acid Production in Aspergillus pseudoterreus

Biological engineering of microorganisms to produce value-added chemicals is a promising route to sustainable manufacturing. However, overproduction of metabolic intermediates at high titer, rate, and yield from inexpensive substrates is challenging in non-model systems where limited information is available regarding metabolic flux and its control in production conditions. Integrated multi-omic analyses of engineered strains offers an in-depth look at metabolites and proteins directly involved in growth and production of target and non-target bioproducts. Here we applied multi-omic analyses to overproduction of the polymer precursor 3-hydroxypropionic acid (3HP) in the filamentous fungus Aspergillus pseudoterreus. A synthetic pathway consisting of aspartate decarboxylase, beta-alanine pyruvate transaminase, and 3HP dehydrogenase was designed and built for A. pseudoterreus. Strains with single- and multi-copy integration events were isolated and multi-omics analysis consisting of intracellular and extracellular metabolomics and targeted and global proteomics was used to interrogate the strains in shake-flask and bioreactor conditions. Production of a variety of co-products (organic acids and glycerol) and oxidative degradation of 3HP were identified as metabolic pathways competing with 3HP production. Intracellular accumulation of nitrogen as 2,4-diaminobutanoate was identified as an off-target nitrogen sink that may also limit flux through the engineered 3HP pathway. Elimination of the high-expression oxidative 3HP degradation pathway by deletion of a putative malonate semialdehyde dehydrogenase improved the yield of 3HP by 3.4 × after 10 days in shake-flask culture. This is the first report of 3HP production in a filamentous fungus amenable to industrial scale biomanufacturing of organic acids at high titer and low pH.

Transcriptomic analysis of the oleaginous yeast Lipomyces starkeyi during lipid accumulation on enzymatically treated corn stover hydrolysate

BACKGROUND

Efficient and economically viable production of biofuels from lignocellulosic biomass is dependent on mechanical and chemical pretreatment and enzymatic hydrolysis of plant material. These processing steps yield simple sugars as well as plant-derived and process-added organic acids, sugar-derived dehydration products, aldehydes, phenolics and other compounds that inhibit the growth of many microorganisms. Lipomyces starkeyi is an oleaginous yeast capable of robust growth on a variety of sugars and lipid accumulation on pretreated lignocellulosic substrates making it attractive as an industrial producer of biofuels. Here, we examined gene expression during batch growth and lipid accumulation in a 20-L bioreactor with either a blend of pure glucose and xylose or pretreated corn stover (PCS) that had been enzymatically hydrolyzed as the carbon sources.

RESULTS

We monitored sugar and ammonium utilization as well as biomass accumulation and found that growth of L. starkeyi is inhibited with PCS hydrolysate as the carbon source. Both acetic acid and furfural are present at concentrations toxic to L. starkeyi in PCS hydrolysate. We quantified gene expression at seven time-points for each carbon source during batch growth and found that gene expression is similar at physiologically equivalent points. Analysis of promoter regions revealed that gene expression during the transition to lipid accumulation is regulated by carbon and nitrogen catabolite repression, regardless of carbon source and is associated with decreased expression of the translation machinery and suppression of the cell cycle. We identified 73 differentially expressed genes during growth phase in the bioreactor that may be involved in detoxification of corn stover hydrolysate.

CONCLUSIONS

Growth of L. starkeyi is inhibited by compounds present in PCS hydrolysate. Here, we monitored key metabolites to establish physiologically equivalent comparisons during a batch bioreactor run comparing PCS hydrolysate and purified sugars. L. starkeyi's response to PCS hydrolysate is primarily at the beginning of the run during growth phase when inhibitory compounds are presumably at their highest concentration and inducing the general detoxification response by L. starkeyi. Differentially expressed genes identified herein during growth phase will aid in the improvement of industrial strains capable of robust growth on substrates containing various growth inhibitory compounds.

Multimodal microfluidic platform for controlled culture and analysis of unicellular organisms

Modern live-cell imaging approaches permit real-time visualization of biological processes, yet limitations exist for unicellular organism isolation, culturing, and long-term imaging that preclude fully understanding how cells sense and respond to environmental perturbations and the link between single-cell variability and whole-population dynamics. Here, we present a microfluidic platform that provides fine control over the local environment with the capacity to replace media components at any experimental time point, and provides both perfused and compartmentalized cultivation conditions depending on the valve configuration. The functionality and flexibility of the platform were validated using both bacteria and yeast having different sizes, motility, and growth media. The demonstrated ability to track the growth and dynamics of both motile and non-motile prokaryotic and eukaryotic organisms emphasizes the versatility of the devices, which should enable studies in bioenergy and environmental research.

Erratum to: A molecular genetic toolbox for Yarrowia lipolytica

[This corrects the article DOI: 10.1186/s13068-016-0687-7.].

A molecular genetic toolbox for Yarrowia lipolytica

BACKGROUND

Yarrowia lipolytica is an ascomycete yeast used in biotechnological research for its abilities to secrete high concentrations of proteins and accumulate lipids. Genetic tools have been made in a variety of backgrounds with varying similarity to a comprehensively sequenced strain.

RESULTS

We have developed a set of genetic and molecular tools in order to expand capabilities of Y. lipolytica for both biological research and industrial bioengineering applications. In this work, we generated a set of isogenic auxotrophic strains with decreased non-homologous end joining for targeted DNA incorporation. Genome sequencing, assembly, and annotation of this genetic background uncovers previously unidentified genes in Y. lipolytica. To complement these strains, we constructed plasmids with Y. lipolytica-optimized superfolder GFP for targeted overexpression and fluorescent tagging. We used these tools to build the "Yarrowia lipolytica Cell Atlas," a collection of strains with endogenous fluorescently tagged organelles in the same genetic background, in order to define organelle morphology in live cells.

CONCLUSIONS

These molecular and isogenetic tools are useful for live assessment of organelle-specific protein expression, and for localization of lipid biosynthetic enzymes or other proteins in Y. lipolytica. This work provides the Yarrowia community with tools for cell biology and metabolism research in Y. lipolytica for further development of biofuels and natural products.

Ochratoxin A production by Penicillium thymicola

Ochratoxin A (OTA) is a mycotoxin produced by some Aspergillus and Penicillium species that grow on economically important agricultural crops and food products. OTA is classified as Group 2B carcinogen and is potently nephrotoxic, which is the basis for its regulation in some jurisdictions. Using high resolution mass spectroscopy, OTA and ochratoxin B (OTB) were detected in liquid culture extracts of Penicillium thymicola DAOMC 180753 isolated from Canadian cheddar cheese. The genome of this strain was sequenced, assembled and annotated to probe for putative genes involved in OTA biosynthesis. Known OTA biosynthetic genes from Penicillium verrucosum or Penicillium nordicum, two related Penicillium species that produce OTA, were not found in P. thymicola. However, a gene cluster containing a polyketide synthase (PKS) and PKS-nonribosomal peptide synthase (NRPS) hybrid encoding genes were located in the P. thymicola genome that showed a high degree of similarity to OTA biosynthetic enzymes of Aspergillus carbonarius and Aspergillus ochraceus. This is the first report of ochratoxin from P. thymicola and a new record of the species in Canada.

Regulation of amino-acid metabolism controls flux to lipid accumulation in Yarrowia lipolytica

Yarrowia lipolytica is a promising microbial cell factory for the production of lipids to be used as fuels and chemicals, but there are few studies on regulation of its metabolism. Here we performed the first integrated data analysis of Y. lipolytica grown in carbon and nitrogen limited chemostat cultures. We first reconstructed a genome-scale metabolic model and used this for integrative analysis of multilevel omics data. Metabolite profiling and lipidomics was used to quantify the cellular physiology, while regulatory changes were measured using RNAseq. Analysis of the data showed that lipid accumulation in Y. lipolytica does not involve transcriptional regulation of lipid metabolism but is associated with regulation of amino-acid biosynthesis, resulting in redirection of carbon flux during nitrogen limitation from amino acids to lipids. Lipid accumulation in Y. lipolytica at nitrogen limitation is similar to the overflow metabolism observed in many other microorganisms, e.g. ethanol production by Sacchromyces cerevisiae at nitrogen limitation.

The Genomes of Three Uneven Siblings: Footprints of the Lifestyles of Three Trichoderma Species

The genus Trichoderma contains fungi with high relevance for humans, with applications in enzyme production for plant cell wall degradation and use in biocontrol. Here, we provide a broad, comprehensive overview of the genomic content of these species for "hot topic" research aspects, including CAZymes, transport, transcription factors, and development, along with a detailed analysis and annotation of less-studied topics, such as signal transduction, genome integrity, chromatin, photobiology, or lipid, sulfur, and nitrogen metabolism in T. reesei, T. atroviride, and T. virens, and we open up new perspectives to those topics discussed previously. In total, we covered more than 2,000 of the predicted 9,000 to 11,000 genes of each Trichoderma species discussed, which is >20% of the respective gene content. Additionally, we considered available transcriptome data for the annotated genes. Highlights of our analyses include overall carbohydrate cleavage preferences due to the different genomic contents and regulation of the respective genes. We found light regulation of many sulfur metabolic genes. Additionally, a new Golgi 1,2-mannosidase likely involved in N-linked glycosylation was detected, as were indications for the ability of Trichoderma spp. to generate hybrid galactose-containing N-linked glycans. The genomic inventory of effector proteins revealed numerous compounds unique to Trichoderma, and these warrant further investigation. We found interesting expansions in the Trichoderma genus in several signaling pathways, such as G-protein-coupled receptors, RAS GTPases, and casein kinases. A particularly interesting feature absolutely unique to T. atroviride is the duplication of the alternative sulfur amino acid synthesis pathway.

Multi-omics analysis reveals regulators of the response to nitrogen limitation in Yarrowia lipolytica

Background

Yarrowia lipolytica is an oleaginous ascomycete yeast that stores lipids in response to limitation of nitrogen. While the enzymatic pathways responsible for neutral lipid accumulation in Y. lipolytica are well characterized, regulation of these pathways has received little attention. We therefore sought to characterize the response to nitrogen limitation at system-wide levels, including the proteome, phosphoproteome and metabolome, to better understand how this organism regulates and controls lipid metabolism and to identify targets that may be manipulated to improve lipid yield.

Results

We found that ribosome structural genes are down-regulated under nitrogen limitation, during which nitrogen containing compounds (alanine, putrescine, spermidine and urea) are depleted and sugar alcohols and TCA cycle intermediates accumulate (citrate, fumarate and malate). We identified 1219 novel phosphorylation sites in Y. lipolytica, 133 of which change in their abundance during nitrogen limitation. Regulatory proteins, including kinases and DNA binding proteins, are particularly enriched for phosphorylation. Within lipid synthesis pathways, we found that ATP-citrate lyase, acetyl-CoA carboxylase and lecithin cholesterol acyl transferase are phosphorylated during nitrogen limitation while many of the proteins involved in β-oxidation are down-regulated, suggesting that storage lipid accumulation may be regulated by phosphorylation of key enzymes. Further, we identified short DNA elements that associate specific transcription factor families with up- and down-regulated genes.

Conclusions

Integration of metabolome, proteome and phosphoproteome data identifies lipid accumulation in response to nitrogen limitation as a two-fold result of increased production of acetyl-CoA from excess citrate and decreased capacity for β-oxidation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2471-2) contains supplementary material, which is available to authorized users.

Comprehensive Metabolomic, Lipidomic and Microscopic Profiling of Yarrowia lipolytica during Lipid Accumulation Identifies Targets for Increased Lipogenesis

Yarrowia lipolytica is an oleaginous ascomycete yeast that accumulates large amounts of lipids and has potential as a biofuel producing organism. Despite a growing scientific literature focused on lipid production by Y. lipolytica, there remain significant knowledge gaps regarding the key biological processes involved. We applied a combination of metabolomic and lipidomic profiling approaches as well as microscopic techniques to identify and characterize the key pathways involved in de novo lipid accumulation from glucose in batch cultured, wild-type Y. lipolytica. We found that lipids accumulated rapidly and peaked at 48 hours during the five day experiment, concurrent with a shift in amino acid metabolism. We also report that exhaustion of extracellular sugars coincided with thickening of the cell wall, suggesting that genes involved in cell wall biogenesis may be a useful target for improving the efficiency of lipid producing yeast strains.

Methylome reorganization during in vitro dedifferentiation and regeneration of Populus trichocarpa

BACKGROUND

Cytosine DNA methylation (5mC) is an epigenetic modification that is important to genome stability and regulation of gene expression. Perturbations of 5mC have been implicated as a cause of phenotypic variation among plants regenerated through in vitro culture systems. However, the pattern of change in 5mC and its functional role with respect to gene expression, are poorly understood at the genome scale. A fuller understanding of how 5mC changes during in vitro manipulation may aid the development of methods for reducing or amplifying the mutagenic and epigenetic effects of in vitro culture and plant transformation.

RESULTS

We investigated the in vitro methylome of the model tree species Populus trichocarpa in a system that mimics routine methods for regeneration and plant transformation in the genus Populus (poplar). Using methylated DNA immunoprecipitation followed by high-throughput sequencing (MeDIP-seq), we compared the methylomes of internode stem segments from micropropagated explants, dedifferentiated calli, and internodes from regenerated plants. We found that more than half (56%) of the methylated portion of the genome appeared to be differentially methylated among the three tissue types. Surprisingly, gene promoter methylation varied little among tissues, however, the percentage of body-methylated genes increased from 9% to 14% between explants and callus tissue, then decreased to 8% in regenerated internodes. Forty-five percent of differentially-methylated genes underwent transient methylation, becoming methylated in calli, and demethylated in regenerants. These genes were more frequent in chromosomal regions with higher gene density. Comparisons with an expression microarray dataset showed that genes methylated at both promoters and gene bodies had lower expression than genes that were unmethylated or only promoter-methylated in all three tissues. Four types of abundant transposable elements showed their highest levels of 5mC in regenerated internodes.

CONCLUSIONS

DNA methylation varies in a highly gene- and chromosome-differential manner during in vitro differentiation and regeneration. 5mC in redifferentiated tissues was not reset to that in original explants during the study period. Hypermethylation of gene bodies in dedifferentiated cells did not interfere with transcription, and may serve a protective role against activation of abundant transposable elements.

Bacterial resistance to antisense peptide phosphorodiamidate morpholino oligomers

Peptide phosphorodiamidate morpholino oligomers (PPMOs) are synthetic DNA mimics that bind cRNA and inhibit bacterial gene expression. The PPMO (RFF)(3)RXB-AcpP (where R is arginine, F, phenylalanine, X is 6-aminohexanoic acid, B is β-alanine, and AcpP is acyl carrier protein) is complementary to 11 bases of the essential gene acpP (which encodes acyl carrier protein). The MIC of (RFF)(3)RXB-AcpP was 2.5 μM (14 μg/ml) in Escherichia coli W3110. The rate of spontaneous resistance of E. coli to (RFF)(3)RXB-AcpP was 4 × 10(-7) mutations/cell division. A spontaneous (RFF)(3)RXB-AcpP-resistant mutant (PR200.1) was isolated. The MIC of (RFF)(3)RXB-AcpP was 40 μM (224 μg/ml) for PR200.1. The MICs of standard antibiotics for PR200.1 and W3110 were identical. The sequence of acpP was identical in PR200.1 and W3110. PR200.1 was also resistant to other PPMOs conjugated to (RFF)(3)RXB or peptides with a similar composition or pattern of cationic and nonpolar residues. Genomic sequencing of PR200.1 identified a mutation in sbmA, which encodes an active transport protein. In separate experiments, a (RFF)(3)RXB-AcpP-resistant isolate (RR3) was selected from a transposome library, and the insertion was mapped to sbmA. Genetic complementation of PR200.1 or RR3 with sbmA restored susceptibility to (RFF)(3)RXB-AcpP. Deletion of sbmA caused resistance to (RFF)(3)RXB-AcpP. We conclude that resistance to (RFF)(3)RXB-AcpP was linked to the peptide and not the phosphorodiamidate morpholino oligomer, dependent on the composition or repeating pattern of amino acids, and caused by mutations in sbmA. The data further suggest that (RFF)(3)R-XB PPMOs may be transported across the plasma membrane by SbmA.

Genome resequencing reveals multiscale geographic structure and extensive linkage disequilibrium in the forest tree Populus trichocarpa

• Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype-genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination. • We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29,213 single-nucleotide polymorphisms. • Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r(2) dropping below 0.2 within 3-6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (N(e) ≈ 4000-6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features. • Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

Dynamic DNA cytosine methylation in the Populus trichocarpa genome: tissue-level variation and relationship to gene expression

Background

DNA cytosine methylation is an epigenetic modification that has been implicated in many biological processes. However, large-scale epigenomic studies have been applied to very few plant species, and variability in methylation among specialized tissues and its relationship to gene expression is poorly understood.

Results

We surveyed DNA methylation from seven distinct tissue types (vegetative bud, male inflorescence [catkin], female catkin, leaf, root, xylem, phloem) in the reference tree species black cottonwood (Populus trichocarpa). Using 5-methyl-cytosine DNA immunoprecipitation followed by Illumina sequencing (MeDIP-seq), we mapped a total of 129,360,151 36- or 32-mer reads to the P. trichocarpa reference genome. We validated MeDIP-seq results by bisulfite sequencing, and compared methylation and gene expression using published microarray data. Qualitative DNA methylation differences among tissues were obvious on a chromosome scale. Methylated genes had lower expression than unmethylated genes, but genes with methylation in transcribed regions ("gene body methylation") had even lower expression than genes with promoter methylation. Promoter methylation was more frequent than gene body methylation in all tissues except male catkins. Male catkins differed in demethylation of particular transposable element categories, in level of gene body methylation, and in expression range of genes with methylated transcribed regions. Tissue-specific gene expression patterns were correlated with both gene body and promoter methylation.

Conclusions

We found striking differences among tissues in methylation, which were apparent at the chromosomal scale and when genes and transposable elements were examined. In contrast to other studies in plants, gene body methylation had a more repressive effect on transcription than promoter methylation.

Library preparation and data analysis packages for rapid genome sequencing

High-throughput sequencing (HTS) has quickly become a valuable tool for comparative genetics and genomics and is now regularly carried out in laboratories that are not connected to large sequencing centers. Here we describe an updated version of our protocol for constructing single- and paired-end Illumina sequencing libraries, beginning with purified genomic DNA. The present protocol can also be used for "multiplexing," i.e. the analysis of several samples in a single flowcell lane by generating "barcoded" or "indexed" Illumina sequencing libraries in a way that is independent from Illumina-supported methods. To analyze sequencing results, we suggest several independent approaches but end users should be aware that this is a quickly evolving field and that currently many alignment (or "mapping") and counting algorithms are being developed and tested.

Genome-wide high throughput analysis of DNA methylation in eukaryotes

Cytosine methylation is the quintessential epigenetic mark. Two well-established methods, bisulfite sequencing and methyl-DNA immunoprecipitation (MeDIP) lend themselves to the genome-wide analysis of DNA methylation by high throughput sequencing. Here we provide an overview and brief review of these methods. We summarize our experience with MeDIP followed by high throughput Illumina/Solexa sequencing, exemplified by the analysis of the methylated fraction of the Neurospora crassa genome ("methylome"). We provide detailed methods for DNA isolation, processing and the generation of in vitro libraries for Illumina/Solexa sequencing. We discuss potential problems in the generation of sequencing libraries. Finally, we provide an overview of software that is appropriate for the analysis of high throughput sequencing data generated by Illumina/Solexa-type sequencing by synthesis, with a special emphasis on approaches and applications that can generate more accurate depictions of sequence reads that fall in repeated regions of a chosen reference genome.