Introduction

Synthetic biology doesn’t stop at the level of molecular systems. To expand the scope of this project we have chosen for a multi-level approach, in which we are working on biobricks, proteins, a pathway and also our host. In order to achieve the latter two strategies have been conceived of. In the first we have chosen to harness the power of directed evolution, a powerful tool that not often used in this competition. An explanation for this might lie in the fact that this approach is only semi-rational at best, however we like to argue that this does not make it any less of a powerful mechanism, and neither within the field of synthetic biology. Nonetheless, in order to explore new territories we have chosen for a second, fully rational approach in which we analyze the transcriptome of two distinct phenotypes; the mycelial and the single cell.

Rationale

The secretion capacity of Aspergillus niger is the feature mainly contributing to its status as excellent industrial workhorse. However, when we investigate this process in more detail, we find that only the hyphal tips of the mycelium are actively secreting. Since the vegetative mycelium poses a burden, generating single cells.

From a process-oriented perspective the potential of host engineering is also interesting. A single cellular phenotype will result in a higher surface to volume ratio, thus effectively increasing the exchange area, while at the same time it annihilates pore clogging. A unicellular strain could be cultured in a similar fashion as a yeast, such as Saccharomyces cerevisiae . The fact that this research intersects with a fundamental topic, the evolution of multicellularity, makes it even more interesting.

Figure 1) Schematic representation bioreactor

Strategy 1: Directed evolution

The power of experimental evolution with regard to complex adaptations has been demonstrated in recent research, allowing acquisition of multicellular Saccharomyces cerevisiae. The other way around one can envision evolution of a unicellular mutant of the filamentous fungus Aspergillus niger. Obtaining such a mutant could help greatly in identifying genes related to multicellularity in complex fungi, a field that is largely unexplored.

A scientific paper on multicellular yeast by directed evolution

Ratcliff, W. C., R. F. Denison, et al. (2012). "Experimental evolution of multicellularity." Proc Natl Acad Sci U S A 109(5): 1595-1600.

Ratcliff et al. subjected the unicellular yeast Saccharomyces cerevisiae to an environment in which they expected multicellularity to be adaptive. Rapid evolution of clustering genotypes that display a novel multicellular life history characterized by reproduction via multicellular propagules were observed. Simple among-cell division of labor rapidly evolved. Early multicellular strains were composed of physiologically similar cells, but these subsequently evolved higher rates of programmed cell death (apoptosis), an adaptation that increases propagule production. This shows that key aspects of multicellular complexity can readily evolve from unicellular eukaryotes. This made us ponder and let to the idea of generating a single cellular phenotype in a reverse approach.

Aim

Obtain a single cell phenotypic Aspergillus niger strain by directed evolution

Approach

In this approach we select for cells with a reduced mycelial cohesiveness by using filters with different pore-sizes. The procedure is iterative; we grow the cells, vortex them, filter them and then grow the cells that were able to get trough the filter for the next round.

Figure 2) The iterative procedure in this directed evolution approach.

Research Methods

Mutagenised spores of A. niger N593 from which evolved strains were cultivated are stored at -80°C according to Appendix C for possible genomic comparison. Spores from A. niger N593 spores were mutagenized by exposure to a Philips TUV 30W lamp for 10 to 60 seconds. In order to apply selective pressure for the desired phenotypic trait, filters with different pore sizes are used. For filter steps in the evolution experiments BD Falcon® Cell Strainers of 40, 70 and 100 µm pore size and SEFAR NITEX® 3-150/50 nylon filter gauze of 150 µm pore size were used. Stainless steel mesh ‘cups’ with pore size 20 µm from Anping Yuansheng Mesh Cooperation were also used.

Figure 3) Materials used as filters

Strategy 2: Comparative transcriptomics

RNA sequencing is a next generation sequencing technology that is rapidly replacing the conventional DNA microarrays. Because, unlike microarrays, RNA sequencing does not rely on probes or primers, there is a smaller bias. It allows for more than analyzing differential gene expression, as it also allows for discovery of novel RNAs and analysis of isoforms of genes. Another advantage over DNA microarrays of RNA sequencing is that data can be reanalyzed once more information on the transcriptome becomes available. Combined with the existence of dimorphisms this allows for an interesting opportunity when it comes to the investigation of the genes involved in this phenotypic distinction.

A scientific paper from 1971

Anderson, J. G. and J. E. Smith (1972). "Effects of Elevated-Temperatures on Spore Swelling and Germination in Aspergillus Niger." Canadian Journal of Microbiology 18 (3): 289-297.

Anderson and Smith found that at 44C germ-tube formation was completely inhibited in Aspergillus niger, although spherical growth could occur over a prolonged period to produce large spherical cells. More generally, there are more dimorphic fungi that display such a distinctive phenotypic transition at elevated temperatures. By comparing transcriptional profiles of distinct phenotypical morphologies we can obtain insight in the genes causative to multicellularity.

Aim

Finding sets of candidate genes causative to the single cell phenotype in Aspergillus niger by transcriptome analysis.

Approach

Spores of A. niger N400 are grown at three distinct conditions. Growing A. niger in liquid medium at 45C will yield a single cell phenotype. When grown at the same temperature on a solid medium mycelium is formed. Mycelium is also formed when A. niger is grown at 30C. Thus, by varying temperature or phase-state of the medium one can obtain either a multi- or single cellular phenotype. Transcriptome analysis allows identification of the transcripts that are uniquely present or absent in the single cell phenotype. This knowledge can be used to genetically modify A. niger to obtain a single cell phenotype at a broad range of conditions, such as at room temperature.

Figure 4)

Research Methods

After cultivation of the Aspergilli at different conditions, RNA needs to be extracted according to Appendix H. This RNA will be sent out to a sequencing company and after some amount of time the data will be returned. For the analysis of the data a transcriptomics pipeline needs to be constructed. Since there is a reference genome available we don’t need to perform a de novo assembly.

Tophat

RNA-Seq generates millions of short sequence reads and therefore the mapping of RNA-Seq is an intensive computational task. This process is performed with the use of specific software. One of these software packages is Tophat. Most mapping algorithms depend on the known splice junctions. Tophat on the other hand is designed to map RNA-Seq reads without relying on these known splice junctions. Tophat detects splice sites ab initio by identifying reads that span exon junctions. The Tophat algorithm first maps the non-junction and short reads using Bowtie. Next the program creates a consensus of mapped reads. Then TopHat breaks up the initially unaligned reads into smaller pieces, Constructing a seed table and tracking of possible junctions, enabling them to be matched onto the reference genome.

Figure 5) Mapping reads onto the Aspergillus niger reference genome allows for discovery of patterns in gene expression that are unique to the single cell phenotype.

Cufflinks

Cufflinks was developed to investigate differential gene expression. This program is designed to calculate the abundances of transcripts constructing an overlap graph. This graph is used to calculate the maximum likelihood of transcripts resulting in Fragment Per Kilobase per Milion (FPKM) values. This value, which represents the quantity of mapped fragments relative to the length of the transcript, normalizes for transcript length and machine yield such that expression of different transcripts (transcripts.gtf) can be compared. Expression levels of genes can simply be determined by summing the FPKM values of their respective isoforms.

Downstream analysis

The statistical program ‘R’ will be used to visualise the final data output (see r-project.org). Extraction of gene ID’s from the fasta files into a text file allows for manual input of this data into cytoscape (see cytoscape.org). A hierarchical network can be inferred for all three Gene Ontology (GO) annotations for visualization of overrepresented terms in specific (sets of) nodes.

Results: growth and metabolic activity

One of the main questions that arises in our comparative transcriptomics approach is whether the so-called single cells that are formed at 45C are metabolically active, or are rather just crippled swollen cells. In our investigation we first monitor the growth of the single cells by a time-series analysis. Next, we assess the metabolic activity of the cells by several methods in which we attempt to examine division of cells, their nuclei and transcriptional and translation activity.

Spherical growth

Over an time interval of 32 hours, every four hours a series of photo's is taken to obtain insight in the development of the cells and their spherical growth. Most if not all of the cells survive this high temperature condition. Even after 96 hours most cells are still intact, however irregular structures within the cell give the appearance that the cells are not as healthy as they used to be. However, since there is transcriptomics data available of different conditions of another strain of Aspergillus niger, we have chosen to switch strains for this approach to allow for comparison of the transcriptional landscapes.

Figure) Growth of A. niger N593 at 44C. From left to right, top to bottom: 0h, 4h, 8h, 12h, 16h, 20h, 24, 28h, 32h, 96h.

Change of Strain

As described in literature, N593 formed giant cells after 24h at 44C. Since transcriptome data on N400 from multiple stages in its life cycle is available, this strain is chosen such that this research complements the current RNA landscape profile of A. niger and allows for comparison of data from different life stages. However, unlike N400, N593 repeatedly formed mycelium at 44C. To overcome this effect the temperature was increased to 45C, at which a single cell phenotype was obtained for N400.

Figure) Distinct morphologies N400 and N593 grown at 44C.

Beside from the occurrence or germination at 44C, another noticeable different is the size. It appears that the mutations that N593 has obtained in comparison to the mother strain, N400, are causative to a reduction in size as well as a reduction in the maximum temperature at which germination can occur.

Finding the right conditions

There are three conditions from which we wish to obtain samples; 30C mycelium, 45C single cells and 45C mycelium. The first of these conditions is often encountered in the environment of A. niger and thus an familiar phenotype. The 45C single cell can be obtained as shown above. This leaves us with the challenge of obtaining a multicellular phenotype at 45C. In order to obtain this phenotype we intend to grow A. niger N400 on a plate without shaking it at 250RPM.

Figure) Searching for a multicellular phenotype at 45C

Calcofluor Staining

Figure) Calcofluor staining of giant cells

DAPI Staining

Figure) DAPI staining indicated nuclei divide

expressing GFP

Figure) GFP transformed A. niger grown at 45C.

Results directed evolution

Evolution experiments can be viewed as a cycle of growth and selection, as depicted above in figure 2. Selection is the evolutionary pressure exerted on the population of cells, which in this case is based on size. The most important thing was to find the largest constant evolutionary pressure, and thus the smallest evolutionary bottleneck, without it leading to major extinction of the parallel evolving lines. To increase the efficiency of this process we have chosen to irradiate the spores first to generate a larger genotypic diversity.

UV mutagenesis

Small petri dishes were filled with 2 ml spore suspension of 108 spores/ml and exposed to UV radiation for 10, 20 and 40 seconds respectively. Generating extra mutations gives the population a mutational ‘head start’ and may speed up the generation of beneficial mutations. On the basis of the survival rate, which was ~90% after 10 seconds, these spores were chosen for evolution experiments with 10 ml liquid cultures. For the evolution experiment of 50 ml liquid cultures diluted, larger batches of spores were used to enhance homogeneous distribution of UV radiation. A spore suspension was diluted to 5*10^6 spores/ml and batches of 10 ml were exposed to UV in standard 94 mm petri dishes for 10, 20, 40 and 60 seconds. The same approach was used to estimate survival rate and in this case the spores that had been radiated for 60 seconds showed a survival rate of 70% and were chosen.

Methodology

Since there is not a standardised protocol, the protocol that has been developed itself can be considered a result. Six parallel lines of 10 ml CM liquid cultures were used to be able to note similar or divergent patterns in evolution. Multiple lines are needed also to make sure the experiment can still continue when part of the lines are lost due to for instance severe infection with bacteria. To initiate the first round of selection, six liquid cultures of 10 ml CM were inoculated with mutagenised spores to 106 spores/ml and incubated at 30°C and 250 RPM. Selection was imposed and growth was allowed to resume under the same conditions for approximately 24 hours before the next round of selection.

The iterative cycle:
1.Vortex the culture at max speed.
2.Filter the culture.
3.Spin down the filtrate (15 min @ 4800g).
4.Decant the supernatant carefully.
5.Resuspend the pellet in CM.
6.Inoculate at least half of the resuspension in 10 ml fresh CM.
7.Use the other part to inoculate CM plates for spores or store it at 4°C as backup material.

Variables in this scheme are vortexing time, filter pore size and complete or incomplete decanting. Vortexing was done holding the liquid culture tube upright, standardising shearing rate. More vortexing time would create larger shearing stress and thus a less strict selection on more easily separable viable mycelium, so less evolutionary pressure.
Filtration creates a size threshold, excluding everything with a size above it. This threshold is mainly dependent on the filter pore size. Larger pore size would mean a less strict selection. Filtrate was collected and centrifuged in plastic 50 ml Greiner tubes. After centrifugation the samples were decanted to prevent factors affecting growth to be transfered from one round to the other as much as possible. However, complete decantation results in loss of cells because the pellet is not dense enough, thus the choice was made to let the last fraction of medium remain. Shaking the samples at 250RPM is important, cause when cells get stuck to the glass they can sporulate. If this occurs the samples is rendered unsuitable for the next round of selection, since spores can easily permeate the filter which results in a loss of evolutionary pressure. The final protocol deemed successful thus consists of the variables 10s of vortexing, filtering over 150 µm filter and incomplete decanting of the medium. A detailed version can be found under Appendix F.
A protocol for an evolution experiment with 3 larger liquid cultures of 50 ml in 250 ml Erlenmeyers was based on this. Variables are the same, except for the vortexing step, which showed to be non-effective in a test run. Instead, a stirrer rod of 50 mm length was spun at 300 RPM to disrupt and disperse mycelium. In six rounds of selection, one extinction event occurred. A detailed protocol can be found under Appendix G.

For integrity of the parallel lines employed in the evolution experiments, mycelium-inoculated plates for harvest of spores should be kept of all lines and if necessary, lines should be recovered from spores of the same line.

Characterisation evolved strains

Experimental evolution lines are coded ‘x-yTz’, where ‘x’ is the line number, ‘y’ is the number of selection rounds that have been applied, ‘T’ is short for transfer, the practical term for selection round, and ‘z’ is the filter pore size used. For evolved strains resulting from the experiment employing 50 ml liquid cultures, an ‘L’ was added for ‘Larger volume’.

Morphology

When plated on CM and incubated at 30°C, evolved strains show the same morphology and growth rate as the parent strain N593. Also at 44°C no difference was found. However, when grown in liquid culture, a different macromorphology shows between the evolved and the parent strain. All evolved strains, and 12T150 in particular, have developed a more pellet-like macromorphology. 12T150 shows the presence of pellets that appear to be of 2 distinct sizes. Microscopic analysis however showed no different micromorphology.

Figure) Parent strain N593 (left) and evolved strain 12T150 (right) in liquid culture. The lower pictures show liquid cultures poured in a petri dish.

Organic acid secretion

Oxalic acid production values calculated from UV measurements are given in figure 23. No significant decrease in production was found. 12T150 even shows a higher oxalic acid production. Data points even stronger in this direction when the final yield is calculated from the production values at 72h and the biomass measured for each liquid culture individually (figure 24). 12T150 shows a significant increase in yield of 57 % over N593 (P = 0.00456; two-tailed t-test, unequal variances).

Figure) Oxalic acid production during 72 hours of growth in 10 ml production medium.

Figure) Oxalic acid production per amount of biomass after 72 hours of growth in 10 ml production medium

Reduced mycelial cohesiveness

Quantification of mycelial cohesiveness was done by analysing the amount of CFUs per amount of mycelium able to permeate a filter, or permeative ability. All filtration steps were done using a pore size of 150 µm. The general layout for this analysis would be:

1.growth of biomass in 10 ml CM (30°C, 250 RPM)
2.determination of ungerminated spores
a. filtering supernatant (1 ml)
b. plating filtrate (100 µl)
3.determination of total permeative CFUs
a. vortexing culture for 30 seconds
b. filtering the complete culture
c. plating filtrate (several volumes)

Figure) Permeative ability of parent strain N593 and evolved strains by CFU count, including a check for ungerminated spores by plating filtered supernatant without prior disruption of the culture. Counts were made up to 200 CFUs per plate, causing a negative bias in the 1000 µl plates, which suffered from that.

Results comparative transcriptomics

RNA StdSens analysis

Pipeline