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.

Unveiling a classical mutant in the context of the GH3 β-glucosidase family in Neurospora crassa

Classical fungal mutant strains obtained by mutagenesis have helped to elucidate fundamental metabolic pathways in the past. In the filamentous fungus Neurospora crassa, the gluc-1 strain was isolated long ago and characterized by its low level of β-glucosidase activity, which is essential for the degradation of cellulose, the most abundant biopolymer on Earth and the main polymeric component of the plant cell wall. Based on genomic resequencing, we hypothesized that the causative mutation resides in the β-glucosidase gene gh3-3 (bgl6, NCU08755). In this work, growth patterns, enzymatic activities and sugar utilization rates were analyzed in several mutant and overexpression strains related to gluc-1 and gh3-3. In addition, different mutants affected in the degradation and transport of cellobiose were analyzed. While overexpression of gh3-3 led to the recovery of β-glucosidase activity in the gluc-1 mutant, as well as normal utilization of cellobiose, the full gene deletion strain Δgh3-3 was found to behave differently than gluc-1 with lower secreted β-glucosidase activity, indicating a dominant role of the amino acid substitution in the point mutated gh3-3 gene of gluc-1. Our results furthermore confirm that GH3-3 is the major extracellular β-glucosidase in N. crassa and demonstrate that the two cellodextrin transporters CDT-1 and CDT-2 are essential for growth on cellobiose when the three main N. crassa β-glucosidases are absent. Overall, these findings provide valuable insight into the mechanisms of cellulose utilization in filamentous fungi, being an essential step in the efficient production of biorefinable sugars from agricultural and forestry plant biomass.

Tuning a high performing multiplexed-CRISPRi Pseudomonas putida strain to further enhance indigoidine production

In this study, a 14-gene edited Pseudomonas putida KT2440 strain for heterologous indigoidine production was examined using three distinct omic datasets. Transcriptomic data indicated that CRISPR/dCpf1-interference (CRISPRi) mediated multiplex repression caused global gene expression changes, implying potential undesirable changes in metabolic flux. ¹³C-metabolic flux analysis (¹³C-MFA) revealed that the core P. putida flux network after CRISPRi repression was conserved, with moderate reduction of TCA cycle and pyruvate shunt activity along with glyoxylate shunt activation during glucose catabolism. Metabolomic results identified a change in intracellular TCA metabolites and extracellular metabolite secretion profiles (sugars and succinate overflow) in the engineered strains. These omic analyses guided further strain engineering, with a random mutagenesis screen first identifying an optimal ribosome binding site (RBS) for Cpf1 that enabled stronger product-substrate pairing (1.6-fold increase). Then, deletion strains were constructed with excision of the PHA operon (ΔphaAZC-IID) resulting in a 2.2-fold increase in indigoidine titer over the optimized Cpf1-RBS construct at the end of the growth phase (∼6 h). The maximum indigoidine titer (at 72 h) in the ΔphaAZC-IID strain had a 1.5-fold and 1.8-fold increase compared to the optimized Cpf1-RBS construct and the original strain, respectively. Overall, this study demonstrated that integration of omic data types is essential for understanding responses to complex metabolic engineering designs and directly quantified the effect of such modifications on central metabolism.

Online Hydrophilic Interaction Chromatography (HILIC) Enhanced Top-Down Mass Spectrometry Characterization of the SARS-CoV-2 Spike Receptor-Binding Domain

SARS-CoV-2 cellular infection is mediated by the heavily glycosylated spike protein. Recombinant versions of the spike protein and the receptor-binding domain (RBD) are necessary for seropositivity assays and can potentially serve as vaccines against viral infection. RBD plays key roles in the spike protein's structure and function, and thus, comprehensive characterization of recombinant RBD is critically important for biopharmaceutical applications. Liquid chromatography coupled to mass spectrometry has been widely used to characterize post-translational modifications in proteins, including glycosylation. Most studies of RBDs were performed at the proteolytic peptide (bottom-up proteomics) or released glycan level because of the technical challenges in resolving highly heterogeneous glycans at the intact protein level. Herein, we evaluated several online separation techniques: (1) C2 reverse-phase liquid chromatography (RPLC), (2) capillary zone electrophoresis (CZE), and (3) acrylamide-based monolithic hydrophilic interaction chromatography (HILIC) to separate intact recombinant RBDs with varying combinations of glycosylations (glycoforms) for top-down mass spectrometry (MS). Within the conditions we explored, the HILIC method was superior to RPLC and CZE at separating RBD glycoforms, which differ significantly in neutral glycan groups. In addition, our top-down analysis readily captured unexpected modifications (e.g., cysteinylation and N-terminal sequence variation) and low abundance, heavily glycosylated proteoforms that may be missed by using glycopeptide data alone. The HILIC top-down MS platform holds great potential in resolving heterogeneous glycoproteins for facile comparison of biosimilars in quality control applications.

An assessment of animal welfare impacts in wild Norway rat (Rattus norvegicus) management

Norway rats (Rattus norvegicus) are considered one of the most significant vertebrate pests globally, because of their impacts on human and animal health. There are legal and moral obligations to minimise the impacts of wildlife management on animal welfare, yet there are few data on the relative welfare impacts of rat trapping and baiting methods used in the UK with which to inform management decisions. Two stakeholder workshops were facilitated to assess the relative welfare impacts of six lethal rat management methods using a welfare assessment model. Fifteen stakeholders including experts in wildlife management, rodent management, rodent biology, animal welfare science, and veterinary science and medicine, participated. The greatest welfare impacts were associated with three baiting methods, anticoagulants, cholecalciferol and non-toxic cellulose baits (severe to extreme impact for days), and with capture on a glue trap (extreme for hours) with concussive killing (mild to moderate for seconds to minutes); these methods should be considered last resorts from a welfare perspective. Lower impacts were associated with cage trapping (moderate to severe for hours) with concussive killing (moderate for minutes). The impact of snap trapping was highly variable (no impact to extreme for seconds to minutes). Snap traps should be regulated and tested to identify those that cause rapid unconsciousness; such traps might represent the most welfare-friendly option assessed for killing rats. Our results can be used to integrate consideration of rat welfare alongside other factors, including cost, efficacy, safety, non-target animal welfare and public acceptability when selecting management methods. We also highlight ways of reducing welfare impacts and areas where more data are needed.

Identifying the gluc-1 and gluc-2 mutations in Neurospora crassa by genome resequencing

Genome resequencing is an efficient strategy for associating mutant phenotypes with physical genomic loci (Baker 2009). A pilot study of this approach demonstrated that the Neurospora crassa genetic map was critical in narrowing the possible candidate mutations in a strain to a small number in a limited, defined region of the genome (McCluskey et al. 2011). In this study, we utilize a resequencing strategy to identify the mutations underlying the gluc-1 and gluc-2 genes in N. crassa.

Evidence of the Involvement of a Cyclase Gene in the Biosynthesis of Ochratoxin A in Aspergillus carbonarius

Ochratoxin A (OTA) is a well-known mycotoxin with wide distribution in food and feed. Fungal genome sequencing has great utility for identifying secondary metabolites gene clusters for known and novel compounds. A comparative analysis of the OTA-biosynthetic cluster in A. steynii, A. westerdijkiae, A. niger, A. carbonarius, and P. nordicum has revealed a high synteny in OTA cluster organization in five structural genes (otaA, otaB, ota, otaR1, and otaD). Moreover, a recent detailed comparative genome analysis of Aspergilli OTA producers led to the identification of a cyclase gene, otaY, located in the OTA cluster between the otaA and otaB genes, encoding for a predicted protein with high similarity to SnoaLs domain. These proteins have been shown to catalyze ring closure steps in the biosynthesis of polyketide antibiotics produced in Streptomyces. In the present study, we demonstrated an upregulation of the cyclase gene in A. carbonarius under OTA permissive conditions, consistent with the expression trends of the other OTA cluster genes and their role in OTA biosynthesis by complete gene deletion. Our results pointed out the involvement of a cyclase gene in OTA biosynthetic pathway for the first time. They represent a step forward in the understanding of the molecular basis of OTA biosynthesis in A. carbonarius.

Major involvement of two laccase genes in conidial pigment biosynthesis in Aspergillus oryzae

Wild-type strains of Aspergillus oryzae develop yellow, yellow-green, green, or brown conidia. Previous reports suggested that the conidiation initiates with the biosynthesis of a yellow pigment YWA1 from acetyl-CoA by a polyketide synthase encoded by wA (AO090102000545). This is followed by the conversion to other pigment by a laccase encoded by yA (AO090011000755). Based on orthologous pathways in other Aspergilli, it is reasonable to hypothesize that in addition to yA, AO090102000546 encoding laccase and AO090005000332 encoding Ayg1-like hydrolase play a role in A. oryzae conidial pigment biosynthesis. However, the involvement of these two genes in conidial pigmentation remains unclear. In this study, we tested this hypothesis by assessing the conidial colors of both disruption and overexpression mutants to verify whether AO090102000546 and AO090005000332 were associated with the conidial pigmentation. Observation of single, double, and triple disruptants of these three genes suggested that conidial pigments were synthesized by two laccase genes, AO090011000755 and AO090102000546, whereas Ayg1-like hydrolase gene AO090005000332 was proven to have no obvious association with the synthesis. This was corroborated by observing the phenotype of each overexpression mutant. Interestingly, AO090005000332 overexpression mutant produced smoky yellow-green conidia, different from the wild-type strain. Thus, the AO090005000332-encoded protein is likely to maintain the enzymatic activity. However, the expression level was observed to be one-third of that of AO090102000546 and one-seventh of that of AO090011000755. Consequently, apparent lack of obvious contribution of AO090005000332 to conidial pigmentation could be attributed to its low expression level. Expression analysis indicated similar profiles in several wild-type strains. KEY POINTS: • Conidial pigment biosynthesis after YWA1 mainly involves two laccases in A. oryzae. • Ayg1-like hydrolase in A. oryzae is not obviously involved in conidial pigmentation. • Conidial color is deemed dependent on expression level of two laccases and hydrolase.

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.

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.

Re-routing of Sugar Catabolism Provides a Better Insight Into Fungal Flexibility in Using Plant Biomass-Derived Monomers as Substrates

The filamentous ascomycete Aspergillus niger has received increasing interest as a cell factory, being able to efficiently degrade plant cell wall polysaccharides as well as having an extensive metabolism to convert the released monosaccharides into value added compounds. The pentoses D-xylose and L-arabinose are the most abundant monosaccharides in plant biomass after the hexose D-glucose, being major constituents of xylan, pectin and xyloglucan. In this study, the influence of selected pentose catabolic pathway (PCP) deletion strains on growth on plant biomass and re-routing of sugar catabolism was addressed to gain a better understanding of the flexibility of this fungus in using plant biomass-derived monomers. The transcriptome, metabolome and proteome response of three PCP mutant strains, ΔlarAΔxyrAΔxyrB, ΔladAΔxdhAΔsdhA and ΔxkiA, grown on wheat bran (WB) and sugar beet pulp (SBP), was evaluated. Our results showed that despite the absolute impact of these PCP mutations on pure pentose sugars, they are not as critical for growth of A. niger on more complex biomass substrates, such as WB and SBP. However, significant phenotypic variation was observed between the two biomass substrates, but also between the different PCP mutants. This shows that the high sugar heterogeneity of these substrates in combination with the high complexity and adaptability of the fungal sugar metabolism allow for activation of alternative strategies to support growth.

Multi-Omics Driven Metabolic Network Reconstruction and Analysis of Lignocellulosic Carbon Utilization in Rhodosporidium toruloides

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides. High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

Strain Construction for Intracellular Metabolic Pathway Localization in Y. lipolytica

Pathway localization by fluorophore or epitope tagging can be accomplished through a multi-staged DNA construct and confirmation process, to generate a series of successfully tagged protein targets. Prerequisite conditions for this process in Y. lipolytica are auxotrophic selection (leu2 or ura3), impaired non-homologous end joining by deletion or impairment of ku70, and plasmids or gene pieces for epitope-selection cassette construction. The general approach for gene tagging can work for C- or N-terminal tags. Gene overexpression from an episomal plasmid can be accomplished through transcript amplification and cloning. C-terminal tagging allows expression of a gene-GFP fusion to be regulated from the endogenous promoter. The epitope-selection cassette also includes a constitutive or highly expressed promoter driving the auxotrophic or other selectable marker gene such as one conferring antifungal or antibiotic resistance. Strains for pathway localization utilize overlap PCR, PEG-based transformation, and a fast DNA preparation for rapid colony screening. Successful transformants can be used for pathway localization and condition-specific response analysis.

Comparative Genomic Analysis of Ochratoxin A Biosynthetic Cluster in Producing Fungi: New Evidence of a Cyclase Gene Involvement

The widespread use of Next-Generation Sequencing has opened a new era in the study of biological systems by significantly increasing the catalog of fungal genomes sequences and identifying gene clusters for known secondary metabolites as well as novel cryptic ones. However, most of these clusters still need to be examined in detail to completely understand the pathway steps and the regulation of the biosynthesis of metabolites. Genome sequencing approach led to the identification of the biosynthetic genes cluster of ochratoxin A (OTA) in a number of producing fungal species. Ochratoxin A is a potent pentaketide nephrotoxin produced by Aspergillus and Penicillium species and found as widely contaminant in food, beverages and feed. The increasing availability of several new genome sequences of OTA producer species in JGI Mycocosm and/or GenBank databanks led us to analyze and update the gene cluster structure in 19 Aspergillus and 2 Penicillium OTA producing species, resulting in a well conserved organization of OTA core genes among the species. Furthermore, our comparative genome analyses evidenced the presence of an additional gene, previously undescribed, located between the polyketide and non-ribosomal synthase genes in the cluster of all the species analyzed. The presence of a SnoaL cyclase domain in the sequence of this gene supports its putative role in the polyketide cyclization reaction during the initial steps of the OTA biosynthesis pathway. The phylogenetic analysis showed a clustering of OTA SnoaL domains in accordance with the phylogeny of OTA producing species at species and section levels. The characterization of this new OTA gene, its putative role and its expression evidence in three important representative producing species, are reported here for the first time.

High-Throughput Large-Scale Targeted Proteomics Assays for Quantifying Pathway Proteins in Pseudomonas putida KT2440

Targeted proteomics is a mass spectrometry-based protein quantification technique with high sensitivity, accuracy, and reproducibility. As a key component in the multi-omics toolbox of systems biology, targeted liquid chromatography-selected reaction monitoring (LC-SRM) measurements are critical for enzyme and pathway identification and design in metabolic engineering. To fulfill the increasing need for analyzing large sample sets with faster turnaround time in systems biology, high-throughput LC-SRM is greatly needed. Even though nanoflow LC-SRM has better sensitivity, it lacks the speed offered by microflow LC-SRM. Recent advancements in mass spectrometry instrumentation significantly enhance the scan speed and sensitivity of LC-SRM, thereby creating opportunities for applying the high speed of microflow LC-SRM without losing peptide multiplexing power or sacrificing sensitivity. Here, we studied the performance of microflow LC-SRM relative to nanoflow LC-SRM by monitoring 339 peptides representing 132 enzymes in Pseudomonas putida KT2440 grown on various carbon sources. The results from the two LC-SRM platforms are highly correlated. In addition, the response curve study of 248 peptides demonstrates that microflow LC-SRM has comparable sensitivity for the majority of detected peptides and better mass spectrometry signal and chromatography stability than nanoflow LC-SRM.

Integrative quantitative-phase and airy light-sheet imaging

Light-sheet microscopy enables considerable speed and phototoxicity gains, while quantitative-phase imaging confers label-free recognition of cells and organelles, and quantifies their number-density that, thermodynamically, is more representative of metabolism than size. Here, we report the fusion of these two imaging modalities onto a standard inverted microscope that retains compatibility with microfluidics and open-source software for image acquisition and processing. An accelerating Airy-beam light-sheet critically enabled imaging areas that were greater by more than one order of magnitude than a Gaussian beam illumination and matched exactly those of quantitative-phase imaging. Using this integrative imaging system, we performed a demonstrative multivariate investigation of live-cells in microfluidics that unmasked that cellular noise can affect the compartmental localization of metabolic reactions. We detail the design, assembly, and performance of the integrative imaging system, and discuss potential applications in biotechnology and evolutionary biology.

Colonies of the fungus Aspergillus niger are highly differentiated to adapt to local carbon source variation

Saprobic fungi, such as Aspergillus niger, grow as colonies consisting of a network of branching and fusing hyphae that are often considered to be relatively uniform entities in which nutrients can freely move through the hyphae. In nature, different parts of a colony are often exposed to different nutrients. We have investigated, using a multi-omics approach, adaptation of A. niger colonies to spatially separated and compositionally different plant biomass substrates. This demonstrated a high level of intra-colony differentiation, which closely matched the locally available substrate. The part of the colony exposed to pectin-rich sugar beet pulp and to xylan-rich wheat bran showed high pectinolytic and high xylanolytic transcript and protein levels respectively. This study therefore exemplifies the high ability of fungal colonies to differentiate and adapt to local conditions, ensuring efficient use of the available nutrients, rather than maintaining a uniform physiology throughout the colony.

101 Dothideomycetes genomes: A test case for predicting lifestyles and emergence of pathogens

Dothideomycetes is the largest class of kingdom Fungi and comprises an incredible diversity of lifestyles, many of which have evolved multiple times. Plant pathogens represent a major ecological niche of the class Dothideomycetes and they are known to infect most major food crops and feedstocks for biomass and biofuel production. Studying the ecology and evolution of Dothideomycetes has significant implications for our fundamental understanding of fungal evolution, their adaptation to stress and host specificity, and practical implications with regard to the effects of climate change and on the food, feed, and livestock elements of the agro-economy. In this study, we present the first large-scale, whole-genome comparison of 101 Dothideomycetes introducing 55 newly sequenced species. The availability of whole-genome data produced a high-confidence phylogeny leading to reclassification of 25 organisms, provided a clearer picture of the relationships among the various families, and indicated that pathogenicity evolved multiple times within this class. We also identified gene family expansions and contractions across the Dothideomycetes phylogeny linked to ecological niches providing insights into genome evolution and adaptation across this group. Using machine-learning methods we classified fungi into lifestyle classes with >95 % accuracy and identified a small number of gene families that positively correlated with these distinctions. This can become a valuable tool for genome-based prediction of species lifestyle, especially for rarely seen and poorly studied species.

A new approach to Cas9-based genome editing in Aspergillus niger that is precise, efficient and selectable

Aspergillus niger and other filamentous fungi are widely used in industry, but efficient genetic engineering of these hosts remains nascent. For example, while molecular genetic tools have been developed, including CRISPR/Cas9, facile genome engineering of A. niger remains challenging. To address these challenges, we have developed a simple Cas9-based gene targeting method that provides selectable, iterative, and ultimately marker-free generation of genomic deletions and insertions. This method leverages locus-specific “pop-out” recombination to suppress off-target integrations. We demonstrated the effectiveness of this method by targeting the phenotypic marker albA and validated it by targeting the glaA and mstC loci. After two selection steps, we observed 100% gene editing efficiency across all three loci. This method greatly reduces the effort required to engineer the A. niger genome and overcomes low Cas9 transformations efficiency by eliminating the need for extensive screening. This method represents a significant addition to the A. niger genome engineering toolbox and could be adapted for use in other organisms. It is expected that this method will impact several areas of industrial biotechnology, such as the development of new strains for the secretion of heterologous enzymes and the discovery and optimization of metabolic pathways.

Investigation of inter- and intraspecies variation through genome sequencing of Aspergillus section Nigri

Aspergillus section Nigri comprises filamentous fungi relevant to biomedicine, bioenergy, health, and biotechnology. To learn more about what genetically sets these species apart, as well as about potential applications in biotechnology and biomedicine, we sequenced 23 genomes de novo, forming a full genome compendium for the section (26 species), as well as 6 Aspergillus niger isolates. This allowed us to quantify both inter- and intraspecies genomic variation. We further predicted 17,903 carbohydrate-active enzymes and 2,717 secondary metabolite gene clusters, which we condensed into 455 distinct families corresponding to compound classes, 49% of which are only found in single species. We performed metabolomics and genetic engineering to correlate genotypes to phenotypes, as demonstrated for the metabolite aurasperone, and by heterologous transfer of citrate production to Aspergillus nidulans. Experimental and computational analyses showed that both secondary metabolism and regulation are key factors that are significant in the delineation of Aspergillus species.

Protein hyperproduction in fungi by design

The secretion of enzymes used by fungi to digest their environment has been exploited by humans for centuries for food and beverage production. More than a century after the first biotechnology patent, we know that the enzyme cocktails secreted by these amazing organisms have tremendous use across a number of industrial processes. Secreting the maximum titer of enzymes is critical to the economic feasibility of these processes. Traditional mutagenesis and screening approaches have generated the vast majority of strains used by industry for the production of enzymes. Until the emergence of economical next generation DNA sequencing platforms, the majority of the genes mutated in these screens remained uncharacterized at the sequence level. In addition, mutagenesis comes with a cost to an organism’s fitness, making tractable rational strain design approaches an attractive alternative. As an alternative to traditional mutagenesis and screening, controlled manipulation of multiple genes involved in processes that impact the ability of a fungus to sense its environment, regulate transcription of enzyme-encoding genes, and efficiently secrete these proteins will allow for rational design of improved fungal protein production strains.

Prediction of Metabolite Concentrations, Rate Constants and Post-Translational Regulation Using Maximum Entropy-Based Simulations with Application to Central Metabolism of Neurospora crassa

We report the application of a recently proposed approach for modeling biological systems using a maximum entropy production rate principle in lieu of having in vivo rate constants. The method is applied in four steps: (1) a new ordinary differential equation (ODE) based optimization approach based on Marcelin’s 1910 mass action equation is used to obtain the maximum entropy distribution; (2) the predicted metabolite concentrations are compared to those generally expected from experiments using a loss function from which post-translational regulation of enzymes is inferred; (3) the system is re-optimized with the inferred regulation from which rate constants are determined from the metabolite concentrations and reaction fluxes; and finally (4) a full ODE-based, mass action simulation with rate parameters and allosteric regulation is obtained. From the last step, the power characteristics and resistance of each reaction can be determined. The method is applied to the central metabolism of Neurospora crassa and the flow of material through the three competing pathways of upper glycolysis, the non-oxidative pentose phosphate pathway, and the oxidative pentose phosphate pathway are evaluated as a function of the NADP/NADPH ratio. It is predicted that regulation of phosphofructokinase (PFK) and flow through the pentose phosphate pathway are essential for preventing an extreme level of fructose 1,6-bisphophate accumulation. Such an extreme level of fructose 1,6-bisphophate would otherwise result in a glassy cytoplasm with limited diffusion, dramatically decreasing the entropy and energy production rate and, consequently, biological competitiveness.

Blocking hexose entry into glycolysis activates alternative metabolic conversion of these sugars and upregulates pentose metabolism in Aspergillus nidulans

Background

Plant biomass is the most abundant carbon source for many fungal species. In the biobased industry fungi, are used to produce lignocellulolytic enzymes to degrade agricultural waste biomass. Here we evaluated if it would be possible to create an Aspergillus nidulans strain that releases, but does not metabolize hexoses from plant biomass. For this purpose, metabolic mutants were generated that were impaired in glycolysis, by using hexokinase (hxkA) and glucokinase (glkA) negative strains. To prevent repression of enzyme production due to the hexose accumulation, strains were generated that combined these mutations with a deletion in creA, the repressor involved in regulating preferential use of different carbon catabolic pathways.

Results

Phenotypic analysis revealed reduced growth for the hxkA1 glkA4 mutant on wheat bran. However, hexoses did not accumulate during growth of the mutants on wheat bran, suggesting that glucose metabolism is re-routed towards alternative carbon catabolic pathways. The creAΔ4 mutation in combination with preventing initial phosphorylation in glycolysis resulted in better growth than the hxkA/glkA mutant and an increased expression of pentose catabolic and pentose phosphate pathway genes. This indicates that the reduced ability to use hexoses as carbon sources created a shift towards the pentose fraction of wheat bran as a major carbon source to support growth.

Conclusion

Blocking the direct entry of hexoses to glycolysis activates alternative metabolic conversion of these sugars in A. nidulans during growth on plant biomass, but also upregulates conversion of other sugars, such as pentoses.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4609-x) contains supplementary material, which is available to authorized users.

Forward genetics screen coupled with whole-genome resequencing identifies novel gene targets for improving heterologous enzyme production in Aspergillus niger

Plant biomass, once reduced to its composite sugars, can be converted to fuel substitutes. One means of overcoming the recalcitrance of lignocellulose is pretreatment followed by enzymatic hydrolysis. However, currently available commercial enzyme cocktails are inhibited in the presence of residual pretreatment chemicals. Recent studies have identified a number of cellulolytic enzymes from bacteria that are tolerant to pretreatment chemicals such as ionic liquids. The challenge now is generation of these enzymes in copious amounts, an arena where fungal organisms such as Aspergillus niger have proven efficient. Fungal host strains still need to be engineered to increase production titers of heterologous protein over native enzymes, which has been a difficult task. Here, we developed a forward genetics screen coupled with whole-genome resequencing to identify specific lesions responsible for a protein hyper-production phenotype in A. niger. This strategy successfully identified novel targets, including a low-affinity glucose transporter, MstC, whose deletion significantly improved secretion of recombinant proteins driven by a glucoamylase promoter.

Correction to: Cloning and Expression of Heterologous Cellulases and Enzymes in Aspergillus niger

The author's family name were incorrectly published in the original version. This has been corrected to read as.

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.

Non-steady state mass action dynamics without rate constants: dynamics of coupled reactions using chemical potentials

Comprehensive and predictive simulation of coupled reaction networks has long been a goal of biology and other fields. Currently, metabolic network models that utilize enzyme mass action kinetics have predictive power but are limited in scope and application by the fact that the determination of enzyme rate constants is laborious and low throughput. We present a statistical thermodynamic formulation of the law of mass action for coupled reactions at both steady states and non-stationary states. The formulation uses chemical potentials instead of rate constants. When used to model deterministic systems, the method corresponds to a rescaling of the time dependent reactions in such a way that steady states can be reached on the same time scale but with significantly fewer computational steps. The relationships between reaction affinities, free energy changes and generalized detailed balance are central to the discussion. The significance for applications in systems biology are discussed as is the concept and assumption of maximum entropy production rate as a biological principle that links thermodynamics to natural selection.

A parts list for fungal cellulosomes revealed by comparative genomics

Cellulosomes are large, multiprotein complexes that tether plant biomass-degrading enzymes together for improved hydrolysis¹. These complexes were first described in anaerobic bacteria, where species-specific dockerin domains mediate the assembly of enzymes onto cohesin motifs interspersed within protein scaffolds¹. The versatile protein assembly mechanism conferred by the bacterial cohesin-dockerin interaction is now a standard design principle for synthetic biology^2,3. For decades, analogous structures have been reported in anaerobic fungi, which are known to assemble by sequence-divergent non-catalytic dockerin domains (NCDDs)⁴. However, the components, modular assembly mechanism and functional role of fungal cellulosomes remain unknown^5,6. Here, we describe a comprehensive set of proteins critical to fungal cellulosome assembly, including conserved scaffolding proteins unique to the Neocallimastigomycota. High-quality genomes of the anaerobic fungi Anaeromyces robustus, Neocallimastix californiae and Piromyces finnis were assembled with long-read, single-molecule technology. Genomic analysis coupled with proteomic validation revealed an average of 312 NCDD-containing proteins per fungal strain, which were overwhelmingly carbohydrate active enzymes (CAZymes), with 95 large fungal scaffoldins identified across four genera that bind to NCDDs. Fungal dockerin and scaffoldin domains have no similarity to their bacterial counterparts, yet several catalytic domains originated via horizontal gene transfer with gut bacteria. However, the biocatalytic activity of anaerobic fungal cellulosomes is expanded by the inclusion of GH3, GH6 and GH45 enzymes. These findings suggest that the fungal cellulosome is an evolutionarily chimaeric structure-an independently evolved fungal complex that co-opted useful activities from bacterial neighbours within the gut microbiome.

Erratum to: A molecular genetic toolbox for Yarrowia lipolytica

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

Diverse data supports the transition of filamentous fungal model organisms into the post-genomics era

Filamentous fungi have been important as model organisms since the beginning of modern biological inquiry and have benefitted from open data since the earliest genetic maps were shared. From early origins in simple Mendelian genetics of mating types, parasexual genetics of colony colour, and the foundational demonstration of the segregation of a nutritional requirement, the contribution of research systems utilising filamentous fungi has spanned the biochemical genetics era, through the molecular genetics era, and now are at the very foundation of diverse omics approaches to research and development. Fungal model organisms have come from most major taxonomic groups although Ascomycete filamentous fungi have seen the most major sustained effort. In addition to the published material about filamentous fungi, shared molecular tools have found application in every area of fungal biology. Similarly, shared data has contributed to the success of model systems. The scale of data supporting research with filamentous fungi has grown by 10 to 12 orders of magnitude. From genetic to molecular maps, expression databases, and finally genome resources, the open and collaborative nature of the research communities has assured that the rising tide of data has lifted all of the research systems together.

Comparative genomics reveals high biological diversity and specific adaptations in the industrially and medically important fungal genus Aspergillus

BACKGROUND

The fungal genus Aspergillus is of critical importance to humankind. Species include those with industrial applications, important pathogens of humans, animals and crops, a source of potent carcinogenic contaminants of food, and an important genetic model. The genome sequences of eight aspergilli have already been explored to investigate aspects of fungal biology, raising questions about evolution and specialization within this genus.

RESULTS

We have generated genome sequences for ten novel, highly diverse Aspergillus species and compared these in detail to sister and more distant genera. Comparative studies of key aspects of fungal biology, including primary and secondary metabolism, stress response, biomass degradation, and signal transduction, revealed both conservation and diversity among the species. Observed genomic differences were validated with experimental studies. This revealed several highlights, such as the potential for sex in asexual species, organic acid production genes being a key feature of black aspergilli, alternative approaches for degrading plant biomass, and indications for the genetic basis of stress response. A genome-wide phylogenetic analysis demonstrated in detail the relationship of the newly genome sequenced species with other aspergilli.

CONCLUSIONS

Many aspects of biological differences between fungal species cannot be explained by current knowledge obtained from genome sequences. The comparative genomics and experimental study, presented here, allows for the first time a genus-wide view of the biological diversity of the aspergilli and in many, but not all, cases linked genome differences to phenotype. Insights gained could be exploited for biotechnological and medical applications of fungi.

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.

Genome sequencing and transcriptome analysis of Trichoderma reesei QM9978 strain reveals a distal chromosome translocation to be responsible for loss of vib1 expression and loss of cellulase induction

BACKGROUND

The hydrolysis of biomass to simple sugars used for the production of biofuels in biorefineries requires the action of cellulolytic enzyme mixtures. During the last 50 years, the ascomycete Trichoderma reesei, the main source of industrial cellulase and hemicellulase cocktails, has been subjected to several rounds of classical mutagenesis with the aim to obtain higher production levels. During these random genetic events, strains unable to produce cellulases were generated. Here, whole genome sequencing and transcriptomic analyses of the cellulase-negative strain QM9978 were used for the identification of mutations underlying this cellulase-negative phenotype.

RESULTS

Sequence comparison of the cellulase-negative strain QM9978 to the reference strain QM6a identified a total of 43 mutations, of which 33 were located either close to or in coding regions. From those, we identified 23 single-nucleotide variants, nine InDels, and one translocation. The translocation occurred between chromosomes V and VII, is located upstream of the putative transcription factor vib1, and abolishes its expression in QM9978 as detected during the transcriptomic analyses. Ectopic expression of vib1 under the control of its native promoter as well as overexpression of vib1 under the control of a strong constitutive promoter restored cellulase expression in QM9978, thus confirming that the translocation event is the reason for the cellulase-negative phenotype. Gene deletion of vib1 in the moderate producer strain QM9414 and in the high producer strain Rut-C30 reduced cellulase expression in both cases. Overexpression of vib1 in QM9414 and Rut-C30 had no effect on cellulase production, most likely because vib1 is already expressed at an optimal level under normal conditions.

CONCLUSION

We were able to establish a link between a chromosomal translocation in QM9978 and the cellulase-negative phenotype of the strain. We identified the transcription factor vib1 as a key regulator of cellulases in T. reesei whose expression is absent in QM9978. We propose that in T. reesei, as in Neurospora crassa, vib1 is involved in cellulase induction, although the exact mechanism remains to be elucidated. The data presented here show an example of a combined genome sequencing and transcriptomic approach to explain a specific trait, in this case the QM9978 cellulase-negative phenotype, and how it helps to better understand the mechanisms during cellulase gene regulation. When focusing on mutations on the single base-pair level, changes on the chromosome level can be easily overlooked and through this work we provide an example that stresses the importance of the big picture of the genomic landscape during analysis of sequencing data.

Expansion of Signal Transduction Pathways in Fungi by Extensive Genome Duplication

Plants and fungi use light and other signals to regulate development, growth, and metabolism. The fruiting bodies of the fungus Phycomyces blakesleeanus are single cells that react to environmental cues, including light, but the mechanisms are largely unknown [1]. The related fungus Mucor circinelloides is an opportunistic human pathogen that changes its mode of growth upon receipt of signals from the environment to facilitate pathogenesis [2]. Understanding how these organisms respond to environmental cues should provide insights into the mechanisms of sensory perception and signal transduction by a single eukaryotic cell, and their role in pathogenesis. We sequenced the genomes of P. blakesleeanus and M. circinelloides and show that they have been shaped by an extensive genome duplication or, most likely, a whole-genome duplication (WGD), which is rarely observed in fungi [3-6]. We show that the genome duplication has expanded gene families, including those involved in signal transduction, and that duplicated genes have specialized, as evidenced by differences in their regulation by light. The transcriptional response to light varies with the developmental stage and is still observed in a photoreceptor mutant of P. blakesleeanus. A phototropic mutant of P. blakesleeanus with a heterozygous mutation in the photoreceptor gene madA demonstrates that photosensor dosage is important for the magnitude of signal transduction. We conclude that the genome duplication provided the means to improve signal transduction for enhanced perception of environmental signals. Our results will help to understand the role of genome dynamics in the evolution of sensory perception in eukaryotes.

The mitochondrial genome of the ethanol-metabolizing, wine cellar mold Zasmidium cellare is the smallest for a filamentous ascomycete

Fungi in the class Dothideomycetes often live in extreme environments or have unusual physiology. One of these, the wine cellar mold Zasmidium cellare, produces thick curtains of mycelia in cellars with high humidity, and its ability to metabolize volatile organic compounds is thought to improve air quality. Whether these abilities have affected its mitochondrial genome is not known. To fill this gap, the circular-mapping mitochondrial genome of Z. cellare was sequenced and, at only 23 743 bp, is the smallest reported for a filamentous fungus. Genes were encoded on both strands with a single change of direction, different from most other fungi but consistent with the Dothideomycetes. Other than its small size, the only unusual feature of the Z. cellare mitochondrial genome was two copies of a 110-bp sequence that were duplicated, inverted and separated by approximately 1 kb. This inverted-repeat sequence confused the assembly program but appears to have no functional significance. The small size of the Z. cellare mitochondrial genome was due to slightly smaller genes, lack of introns and non-essential genes, reduced intergenic spacers and very few ORFs relative to other fungi rather than a loss of essential genes. Whether this reduction facilitates its unusual biology remains unknown.

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.

Compartmentalized microchannel array for high-throughput analysis of single cell polarized growth and dynamics

Interrogating polarized growth is technologically challenging due to extensive cellular branching and uncontrollable environmental conditions in conventional assays. Here we present a robust and high-performance microfluidic system that enables observations of polarized growth with enhanced temporal and spatial control over prolonged periods. The system has built-in tunability and versatility to accommodate a variety of scientific applications requiring precisely controlled environments. Using the model filamentous fungus, Neurospora crassa, our microfluidic system enabled direct visualization and analysis of cellular heterogeneity in a clonal fungal cell population, nuclear distribution and dynamics at the subhyphal level, and quantitative dynamics of gene expression with single hyphal compartment resolution in response to carbon source starvation and exchange. Although the microfluidic device is demonstrated on filamentous fungi, the technology is immediately extensible to a wide array of other biosystems that exhibit similar polarized cell growth, with applications ranging from bioenergy production to human health.

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.

Pregnancy rates with intrauterine insemination: comparing 1999 and 2010 World Health Organization semen analysis norms

Over the past 30 years, The World Health Organization has serially measured norms for human sperm. In this study, 1999 and 2010 semen analysis norms as predictors of pregnancy were compared during intrauterine insemination (IUI). A retrospective cohort study was conducted using data collected from the Stanford Fertility Center, between 2005 and 2007, with 981 couples undergoing 2231 IUI cycles. Collected semen was categorized according to total motile sperm counts (TMSC): 'normal (N.) 1999 TMSC', 'abnormal (AbN.) 1999/N. 2010 TMSC', or 'AbN. 2010 TMSC'. Sample comparison was also based on individual semen parameters: 'N. 1999 WHO', 'AbN. 1999/N. 2010 WHO', or 'AbN. 2010 WHO'. Pregnancy (defined by beta-HCG concentration) rates were calculated. Data were compared using correlation coefficients, t-tests and chi-squared tests, with and without adjusting for confounders. Pregnancy rate comparison based on TMSC ('N. 1999 TMSC', 'AbN. 1999/N. 2010 TMSC' and 'AbN. 2010 TMSC') showed a negative correlation (r = -0.41, P = 0.05). Pregnancy rate did not differ when comparisons were based on the presence of abnormal parameters, even when controlling for confounders. Therefore, TMSC based on the 1999 parameters shows best correlation with pregnancy rate for IUI; updating these norms in 2010 has little clinical implication in infertile populations.