Overview

Western honey bees (Apis mellifera) suffer under the influence of the Varroa mite, traditional control methods bear large disadvantages for beekeepers, as they accumulate in the honey and beeswax and spoil the final product. We wanted to create a new method on how to combat the Varroa mite using the egg yolk protein Vitellogenin.

Since Vitellogenin's involvement in transgenerational immune priming was only recently discovered [1], there is little information on where and how the PAMPs are bound to and transported by Vitellogenin. Hence, the heist to create a Vitellogenin variant optimized for our specific needs was a challenging task. Where in the protein should we start searching if we don't know what to look for? Only after heated discussions and extensive research did we finally develop a comprehensive plan aimed at optimizing Vitellogenin to facilitate binding to pathogen-associated molecular patterns (PAMPs) from fungi, thereby enhancing immunity to entomopathogenic fungi, allowing them to be used to combat Varroa mite infestation.

Starting from the wildtype Vitellogenin, we searched for optimized Vitellogenin in two ways: For the wet lab, we used directed evolution to introduce mutations into the protein and sort the proteins by their binding affinity through fluorescence activated cell sorting (FACS). As a dry lab approach, we did multiple predictions and simulations with individual Vitellogenin variants from the Western honey bee and other wild bee species. If we found a variant of Vitellogenin that exhibits the desired properties, the next step would be to express it in Komagataella phaffii and then test the binding affinity towards fungal PAMPs in a specifically designed binding assay.


Directed Evolution

Designing a protein whose mechanism of function is poorly understood presents numerous challenges, particularly when dealing with a complex protein like Vitellogenin. Its many different domains and functional regions [1] make it almost impossible to identify the residues responsible for PAMP (pathogen associated molecular pattern) affinity, without complex and lengthy trials. This made the prospect of success using exclusively rational design approaches limited.

That is why we decided to also pursue a different approach in the form of directed evolution. For this, a library of sequence variants is generated via random mutagenesis and is screened with a high throughput screening method to isolate mutants with the desired properties. The major advantage of this approach is that no prior knowledge of the proteins mechanisms is needed [2].

Directed evolution has been employed for a wide range of applications including optimization of wholecell biotechnological processes like fermentation [3–6], plant breeding [7], protein optimization [2] and biosensors [8], as well as new opportunities like the prediction of antibiotic resistances [9, 10] and many more [11, 12]. Its relevance is further underlined by the award of the Nobel prize in chemistry to Frances H. Arnold in 2018. She was honoured for developing the first procedures for enzyme directed evolution.

Many commercial proteins have been optimized via directed evolution. Among some of the modified parameters are protein stability and activity at different pH and temperature ranges [13–15], substrate turnover optimization [16, 17] and substrate specificity [18], but basically all properties of a protein can be influenced through a fitting directed evolution approach.

As previously discussed, the parameter we needed to improve was the zymosan-affinity of A. mellifera's’ Vitellogenin. We hypothesized that improving zymosan binding would enable the queen bee to bind and transport more zymosan into the oocytes [19] thus, mediate a higher protection of all offsprings against entomopathogenic fungi, which can be used for treating Varroa mite infestation.

With this aim, we needed to decide on a suitable workflow for generating and screening a library. We quickly discovered that surface display is one of the most widely used and effective methods for affinity engineering. It uses the capability of cells to produce and transport proteins onto their cell surface where they are easily accessible for screening, circumventing laborious purification. Single cell flow cytometry and fluorescent activated cell sorting (FACS) downstream make it possible to screen large libraries and isolate mutants that display proteins with the desired traits [20–24].

Multiple organisms can display mutated libraries such as phages [25, 26], bacterial [27], insect [28], mammalian [29, 30] and yeast cells [31–33]. Out of these options, we chose yeast as the system to use, as it offered the most complete platform with the least disadvantages. They are easier to cultivate and employ in experiments compared to insect or mammalian cells, while still being a eukaryotic expression system, meaning they can incorporate posttranslational modifications and guarantee correct protein conformation. Exactly these properties are not offered by phages and bacteria, making yeast the ideal organism for the screening of Vitellogenin mutant libraries.

Among yeast species, Saccharomyces cerevisiae is the most common choice for the display of protein libraries [31–33]. In general, a yeast surface display construct is composed of several elements (Fig. 1). A cell wall protein (CWP) interacts with a polypeptide on the cell wall, which in turn is typically linked to Glycosylphosphatidylinositol (GPI) [34, 35], thereby anchoring the construct onto the cell surface. The protein of interest (POI) is fused to the CWP via a flexible linker. It allows both the CWP and POI to fold and move independently from one another, increasing the viability of the displayed protein [36–38]. Lastly, an epitope tag is included in the construct. The epitope-tag allows the quantification of the fusion protein expression and the subsequent normalization of the measured binding affinity [24, 39].

There are several CWPs available as anchors. We decided to use the widely employed agglutinin mating protein Aga2p, which forms two covalent disulfide bonds with Agα1p on the cell surface. Aga2p also serves as a secretion factor for the construct, reducing the size and complexity of the expressed protein, as there is no need for a separate factor [31, 34, 40].

We decided to use a (G4S)3 linker as it offers a good compromise between flexibility and length. The polar amino acids also form hydrogen bonds with the surrounding water molecules, reducing the interaction of the linker with the protein domains [41].

The epitopes we used for quantification were the HA epitope and the Myc tag, which are very prevalent in many antibody engineering studies. These often have similar aims to our project in increasing binding affinity [24, 42, 43]. We used the two epitopes as a form of redundancy with the HA at the N-terminus and the cMyc C-terminal. If steric interference would have made the cMyc tag inaccessible for antibody binding, we could use the HA tag as a substitute.

The last decision we had to make for our construct was which parts of the protein to display. Vitellogenin is a large protein with a size of 202 kDA. Considering there was no literature regarding displayed proteins of this size and that full Vitellogenin secretion from yeasts was previously shown to be unsuccessful [44], we decided to display the different Vitellogenin domains.

This figure shows a schematic view of the used surface display construct. It shows a yeast cell with Agα1p on the cell surface. It is connected to Aga2p through to disulfide briges. The Aga2p is connected to HA which is connected to a Vitellogenin domain. The domain is binding a fluorescent labled PAMP. At the c-terminus is a myc tag, binding a flourescent marker. . Next is the HA epitope and the Vitellogenin domain of interest. Attached to the Vitellogenin is the FTSC labeled PAMP zymosan. The cMyc tag is shown at the c terminus with antibodies attached to it, indicated by the red point.
Figure 1: Schematic view of the used surface display construct. Aga2p is anchored to the cell wall via two disulfide bridges formed with Agα1p. The linker (not shown) connects to the rest of the construct. Next is the HA epitope and the Vitellogenin domain of interest. Attached to the Vitellogenin is the FTSC labeled PAMP zymosan. The cMyc tag is shown at the c‑terminus with antibodies attached to it, indicated by the red point.

Analyzing the functional domains [1] (Fig. 2), we discovered, a long sequence between the DUF1943 domain and the van Willebrand factor (vWf), which had no proposed function. As the vWf is predicted to be primarily involved in protein stabilization [1] we deemed this section of the protein non‑relevant to our engineering effort. Additionally, the N-terminal β-barrel domain is a putative transcription factor [45] linked to the protein via the cleavable poly-serin linker [46]. Considering its role, it will not necessarily be transported into the oocytes and likely serves no function in PAMP binding. Therefore, it was also excluded from the surface display constructs alongside the poly-serin linker. These decisions left us with only two domains of interest. The alpha-helical domain which was already suggested to bind PAMPs as a sole domain but with reduced affinity [47] and the adjacent DUF1943 domain which is known to contribute to PAMP binding in fish and coral species [48, 49]. We decided that reviewing the binding affinity of both domains separately, and combined could yield us valuable insight into the mechanism of PAMP binding, which was previously not described for A. mellifera's Vitellogenin. We also decided to try a display of a truncated Vitellogenin lacking the β-barrel domain, poly-serin linker, and ending directly after the vWf. This resulted in four constructs to be assembled and displayed. One displaying just the alpha-helical domain (45.5  kDA), one containing the DUF1943 domain (34.27 kDA), a third with the alpha-helical & DUF1943 region (97.75 kDA) and a lastly a truncated Vitellogenin version far larger than the other domains at 144.99 kDA.

The figure shows the simplyfied amino acid sequence of Vitellogenin and its domains. 1: Vitellogenin sequence. 2: β-barrel + polyserin linker, acts as a transcription factor linked via a cleavable polyserine linker. 3: Alphahelical, Heavily positive charged domain thought to be involved in PAMP binding. 4: DUF1943, contributes to PAMP-binding eith te alpha-helical domain. 5: vWF, suggested to be involved in proetin stabilitation
Figure 2: Apis mellifera Vitellogenin domains. The proposed functions of each domain are listed below. The numbers annotate the amino acids corresponding to the domains.

Our choice to use the Aga2p anchor gave us a clear direction for which yeast strain to use. The widely employed S. cerevisiae EBY100 strain [24, 32, 39] contains the cell wall-bound Agα1p component of the anchor complex under the control of the galactose-inducible GAL1 promotor, making it ideal for our concept. The vector pCTCon2 is well suited for use in this strain, as it allows for the yeast surface display construct to also be expressed under the GAL1 promotor and auxotrophic selection via tryptophan dropout medium (Addgene #41843).

After constructing the display vectors, we needed to generate a highly diverse, non-biased library. We chose error-prone PCR with the nucleotide analogous dPTP and 8-oxo-GTP for that. It requires no prior knowledge of domain function like saturation mutagenesis [50] and is easier to execute in a short time frame compared to DNA shuffling [51]. The here utilized nucleotide analogues can mediate both transitions and transversions, while introducing less GC-bias into the mutated sequence when compared to manganese ion-based error‑prone PCRs [52, 53].

We generated the library by amplifying the plasmid backbone, excluding the insert, with a high-fidelity polymerase and the inserts in error-prone reactions with varying nucleotide analogue concentrations to cover a wider range of mutation rates. Generating libraries with the correct mutation rates is critical when considering, that 30-50% of mutations are deleterious to the protein function [54, 55]. This makes highly mutated libraries unusable for directed evolution, while too little amino acid exchanges result in only a small and less diverse variant pool. We decided to use 0.2; 2 and 20 µM of nucleotide analogues for our library generation to cover a wide range of mutation ranges.

80 bp long overhangs complementary to the plasmid backbone were introduced during the error-prone PCR to enable efficient homologous recombination of the library into the backbone inside the yeast cells [56]. The fragments were purified and concentrated via lyophilisation, pooled with the backbone, and used for yeast transformation via electroporation. The transformed yeasts were first grown in SD+CAA selective medium before being induced in SG+CAA medium inducing the presentation of the generated library (Fig. 3).

The figure shows an overview of the directed evolution approach. A library was created through error-prone PCR. Yeast cells were transformed with mutant library and labled with fluorescent zymosan. Screening of the cells using fluorescent -activated cell sorting.
Figure 3: Schematic overview of our directed evolution approach: initially, we created a Vitellogenin domain mutant library through error prone PCR using pTCON2. Following this, we introduced the library into the Saccharomyces cerevisiae EBY100 strain and induced the expression and display of the Vitellogenin domain mutants. To identify and isolate mutants with the highest affinity for the PAMP, we employed FTSC‑labeled zymosan as the target molecule.

The library was then prepared for FACS by initially labelling the cMyc tag using mouse anti cMyc (ThermoFisher # 910MYCL) and goat antimouse AlexaFluor®633 (ThermoFisher # A-21126) antibodies. We then exposed the cells to FTSC labeled zymosan in varying concentrations. Protein variants with improved binding kinetics are more likely to retain the bound PAMP on the cell surface, leading to positive selection in the FACS. With this we were able to identify Vitellogenin variants with a higher affinity for zymosan, which can improve the transgenerational immune priming of honey bees against enteropathogenic fungi. Enabling their usage against Varroa mites.

Zymosan Labeling

To enable the screening process through FACS, it was crucial to fragment and label the zymosan with a fluorescent dye. As mentioned above, zymosan is a large, water-insoluble polysaccharide (approx. 300 kDa), which is why we initially subjected it to acid hydrolysis. We chose mild conditions to keep the zymosan fragments as large as possible so that they could still be recognized by Vitellogenin as PAMPs while remaining water-soluble. This was accomplished by subjecting them to a solution of 0.5 M HCl at 70 °C for a duration of 2 hours. Following acid hydrolysis, the soluble zymosan fragments were ultracentrifuged through a 3 kDa membrane to exclude very small sugar oligos, which were unlikely to be recognized as PAMPs (Fig. 4).

The figure shows the workflow of the zymosan labeling. Zymosan were labeled with FTSC. Zymosan fragments are centrifuged to seperate zymosan fragments from not usable sugar oligos.
Figure 4: Workflow Description: Acid hydrolysis and ultrafiltration for isolation of water soluble zymosan fragments. This involved incubating the zymosan in a 0.5 M HCl solution at 70 °C for 2 hours. After acid hydrolysis, the resulting soluble zymosan fragments underwent ultracentrifugation through a 3 kDa membrane to exclude small sugar oligos, which are less likely to be recognized as PAMPs.

Subsequently, the larger zymosan fragments were labeled with the fluorescent dye fluorescein-5-thiosemicarbazide (FTSC). This dye reacts with the reducing ends of sugars, causing them to form hydrazones, which are subsequently reduced by sodium cyanoborohydride (NaBH3CN) to the fluorescent zymosan derivates. For ideal labeling conditions, a molar ratio of 8:1 (FTSC to reducing ends) and 70 µL of 0.4 M NaBH3CN per mg of zymosan, dissolved in a mixture of 7:3 DMSO:HOAc, can be used (Fig. 5).

The figure schows the reaction of the zymosan labeling with 5-FTSC (fluorescein-5-thiosemicarbazid).
Figure 5: Reaction of the FTSC labeling: FTSC chemically reacts with the reducing end of sugars, leading to the formation of hydrazones. These hydrazones are then reduced into fluorescent zymosan derivatives using sodium cyanoborohydride (NaBH3CN). For optimal labeling conditions, a molar ratio of 8:1 (FTSC to reducing ends) and the addition of 70 µL of 0.4 M NaBH3CN per milligram of zymosan, dissolved in a mixture of 7:3 DMSO to HOAc was used.

After that, the FTSC-labeled zymosan was again passed through a 3 kDa membrane to remove the unreacted chemicals of the reaction. The FTSC-labeled zymosan was then precipitated in ethanol, a solvent in which the FTSC dye is soluble, to remove residual FTSC. To validate the successful acid hydrolysis and FTSC labeling, various methods were employed. Firstly, the acid hydrolysis was verified through Matrix-Assisted-Laser-Desorption-Ionization with Time-Of-Flight analysis of released ions for mass spectrometry (MALDI TOF MS). Further confirmation was obtained by observing fluorescence of the precipitated zymosan under UV light as well as unlabeled zymosan as a negative control. The protocols for this we created ourselves, as there was no such protocol available. You can find it in contribution.

Results

Acid hydrolysis was verified by MALDI TOF MS. The results are depicted in Figure 6. Glucose has a molecular weight of 180 Da, and when forming the glycosidic bond, water with a molecular weight of 18 Da is eliminated. Considering that zymosan is a polymer comprised of beta-1,3- and beta-1,6-linked glucose monomers, the peaks in the mass spectrum exhibiting a consistent gap of 162 m/z prove that the hydrolysis was successful. Furthermore, the data also reveals a decrease in the quantity of zymosan fragments with increasing molecular size. As evident from the data, the hydrolysis indeed proved effective.

This figure shows the mass spectrometry histogram plot of zymosan fragments after acid hydrolysis. The mass spectrum displays consistent 162 m/z gaps in the peaks, highlighting the zymosan structure, composed of beta-1,3 and beta-1,6 glucose monomers.
Figure 6: Mass Spectrometry histogram plot of zymosan fragments after acid hydrolysis. Remarkably, the mass spectrum displays consistent 162 m/z gaps in the peaks, highlighting the zymosan structure, composed of beta-1,3 and beta-1,6 glucose monomers.

Following the successful fragmentation of the zymosan polymer, we proceeded to label the resulting fragments using the fluorescent dye FTSC. The labeled polymer demonstrated pronounced fluorescence at 492 nm, as depicted in figure 7. Following the washing steps, the flowthrough also exhibited fluorescence. This fluorescence decreased with each subsequent wash step, indicating a successful removal of the residual FTSC. Conversely, the negative control, consisting of unlabeled zymosan fragments, did not exhibit any fluorescence signal.

This figure shows the fluorescence at 492 nm at A: labeled and purified zymosan fragments, B: flowthrough of washing step 1, C: flowthrough of washing step 2, D: flowtrough of washing step 5, E: flowthrough of washing step 4, F) unlabeled zymosan fragments as negative control.
Figure 7: Fluorescence at 492 nm. From left to right: Labeled and purified zymosan fragments flowthrough of washing step 1 flowthrough of washing step 2 flowthrough of washing step 3 flowthrough of washing step 4 unlabeled zymosan fragments as negative control

Generation of the Base Construct

First, we constructed the base surface display plasmid via Gibbson Assembly. It contained the functional elements of the surface display, but not the Vitellogenin domains. We ran colony PCRs, amplifying the complete insert, to verify the presence of positive clones. The expected fragment size was 374 bp, which we could observe for several colonies. We send several colonies that showed the right-sized band to sequencing where the correct assembly of the construct was confirmed (Fig. 8).

Constructs

This figure shows a gel electrophoresis of a colony PCR of the pCTCon2 vector containing the yeast surface display functional units is shown. Several colonies show bands of the expected size of 374 bp. The used marker was a 1 kb plus ladder. This figure shows the sequence of the colony, confirming the correct insertion of the fragment.
Figure 8: In the picture on top, the gel electrophoresis of a colony PCR of the pCTCon2 vector containing the yeast surface display functional units is shown. Several colonies show bands of the expected size of 374 bp. The used marker was a 1 kb plus ladder. The bottom picture shows the sequence of the colony, confirming the correct insertion of the fragment.

Next, we needed to integrate the four different Vitellogenin fragments into the construct. We amplified the domains out of a fully assembled and codon optimized Vitellogenin sequence, ordered from IDT. After transformation we ran colony PCRs again, before sending the constructs to be sequenced to verify their identity.

Vg_α‑helical

This figure shows the schematic of Apis mellifera Vitellogenin. The α helical domain chosen for display is highlighted.
Figure 9: Schematic of Apis mellifera Vitellogenin. The α‑helical domain chosen for display is highlighted.

Two colonies of Vg_α‑helical showed the expected band at ~1400 bp. Both were sent to sequencing, where one was confirmed positive (Fig. 10).

This upper figure shows a gel electrophoresis of a colony PCR of the pCTCon2_Vg_α‑helical on the left. The first band shows the expected size of ~1400 bp while the right band runs to high. The right side shows the sequence result of the colony, confirming it to contain the Vg_α‑helical insert.  The lower figure shows the sequence result of the colony, confirming it to contain the Vg_α‑helical insert.
Figure 10: At the top the gel electrophoresis of a colony PCR of the pCTCon2_Vg_α‑helical is shown. In the left pocket, a 1 kb plus ladder was run. By comparison, the middle band shows the expected size of ~1400 bp while the most right band runs to high. The bottom figure shows the sequence result of the colony, confirming it to contain the Vg_α‑helical insert.

Vg_DUF1943


This figure shows a schematic of Apis mellifera Vitellogenin. The DUF1943 domain chosen for display is highlighted.
Figure 11: Schematic of Apis mellifera Vitellogenin. The DUF1943 domain chosen for display is highlighted.

Two colonies showed the expected band at ~1000 bp. Both were sent for sequencing, and both were confirmed to contain the correct insert (Fig. 12 only one shown).

The upper figure shows a gel electrophoresis of three colony PCRs of the pCTCon2_Vg_DUF1943 in comparison to a 1 kb plus ladder. The two colonies on the left show the expected band size of ~1000 bp. The right side of the figure shows the sequence result of one colony, confirming it to contain the Vg_α-DUF1943 insert The lower figure shows the sequence result of one colony, confirming it to contain the Vg_α-DUF1943 insert
Figure 12: At the top the gel electrophoresis of three colony PCRs of the pCTCon2_Vg_DUF1943 in comparison to a 1 kb plus ladder is shown. The two colonies on the left show the expected band size of ~1000 bp. The bottom figure shows the sequence result of one colony, confirming it to contain the Vg_α-DUF1943 insert.

Vg_α‑helical+DUF1943


This figure shows a schematic of Apis mellifera Vitellogenin. The α‑helical+DUF1943 domain chosen for display is highlighted.
Figure 13: Schematic of Apis mellifera Vitellogenin. The α‑helical+DUF1943 domain chosen for display is highlighted.

Again, two colonies showed the expected band at ~2200 bp. Both were sent to sequencing where they were confirmed to contain the correct construct (Fig. 14 only one sequence shown).

The upper figure shows a gel electrophoresis of five colonies PCR of the pCTCon2_Vg_α‑helical+DUF1943 in comparison to a 1 kb plus ladder. Two colonies show the expected band size of ~2200 bp. The right side of the figure shows the sequence result of one colony, confirming it to contain the Vg_α‑helical+DUF1943 insert. The lower figure shows the sequence result of one colony, confirming it to contain the Vg_α‑helical+DUF1943 insert.
Figure 14: At the top the gel electrophoresis of five colonies PCR of the pCTCon2_Vg_α‑helical+DUF1943 in comparison to a 1 kb plus ladder is shown. Two colonies (third lane, sixth lane) show the expected band size of ~2200 bp. The bottom figure shows the sequence result of one colony, confirming it to contain the Vg_α‑helical+DUF1943 insert.

Vg_truncated


This figure shows a schematic of Apis mellifera Vitellogenin. The truncated domain chosen for display is highlighted.
Figure 15: Schematic of Apis mellifera Vitellogenin. The truncated domain chosen for display is highlighted.

One colony showed a promising fragment at ~3900 bp. It was sent to sequencing where it was confirmed to contain the correct fragment (Fig. 16).

The upper figure showas a gel electrophoresis of a colony PCR of the pCTCon2_Vg_truncated in comparison to a 1 kb plus ladder. The band shows the expected size of ~3900 bp. The right side of the figure shows the sequence result of the colony, confirming it to contain the Vg_truncated insert. The lower figure shows the sequence result of the colony, confirming it to contain the Vg_truncated insert.
Figure 16: At the top the gel electrophoresis of a colony PCR of the pCTCon2_Vg_truncated in comparison to a 1 kb plus ladder is shown. The band shows the expected size of ~3900 bp. The bottom figure shows the sequence result of the colony, confirming it to contain the Vg_truncated insert.

Verification of Protein Display

The first step of our engineering was to verify the correct display of our constructs. We were particularly concerned about the display of our larger constructs, Vg_alpha+DUF and Vg_truncated.

We used eight different controls (Fig. 17) to compare our constructs with a cell control, single antibody controls with the first and secondary antibodies, and a double antibody control. Each control was also run after incubation with zymosan to check for background yeast interaction with the PAMP.

This figure shows a list of controls used to verify the correct display of Vitellogenin domains. 1: Yeast cell 2: Yeast cell + labeled zymosan 3: yeast cell+ primary antibody 4: yeast cell + primary antibody + labeled zymosan 5: yeast cell + secondary antibody 6: yeast cell + secondary antibody + labeled zymosan 7: yeast cell + primary antibody + secondary antibody 8: yeast cell + primary antibody + secondary antibody + labeled zymosan
Figure 17: Controls used to verify the correct display of the Vitellogenin domains.

Before analyzing flow cytometry data, it is important to isolate the so-called “events” relevant to the experiment. An event is any signal measured by a flow cytometer without any differentiation of its source. This means that events originating from debris, dead cells, or other sources need to be excluded before continuing with the analysis. There are several different measurements used to determine which accumulation of events, or “population,” shows the desired cells. Front scatter (FSC) measures the light passing through the event, while side scatter (SSC) measures the light being refracted sideways. The signal of each can be subdivided into a height (e. g., FSC-H), width (e.g., FSC-W), and area (e.g., FSC-A value). Height measures the signal intensity, while width measures the time the event passes through the beam. The integral of the measured event generates the area value. Different combinations of these can be used to investigate the properties of the measured events.

FSC-A vs SSC-A is used to isolate the cells from debris or other impurities. FSC-A indicates cell size, while SSC-A shows the internal complexity, also called granularity. As seen in Fig. 18, we observed a well-defined cell population with no debris and other populations that would indicate dead cells or other cell types. An area, or ‘gate,’ was set marking the cells in it for further analysis. We then checked if our sample contained cell aggregations. Cells can either stick together, or can by chance, be measured as a single event. As the presence of such multicell events would interfere with data interpretation they need to be sorted out. We did this via FSC-A vs FSC-H plots. Cell aggregations would appear as a lower population because of their different FSC-A composition. They have a larger FSC-W value due to their size with a lower FSC-H, which they can be sorted by. No subpopulations were observed in FSC‑A/FSC‑H, indicating that only single cells were measured with no cell aggregations present.

This figure shows a graphic for the cell control for both antibodies.
A
This figure shows a graphic for the cell control for both antibodies.
B
Figure 18: A) Front scatter area (FSC-A) vs side scatter area (SSC-A) plot of the cell control + double antibodies + 0.2 mg of zymosan is shown. One distinct population is visible in this plot, indicating that no other cell morphologies or dead cells are present, as well as no other debris. B) Front scatter area (FSC-A) vs front scatter height (FSC-H) plot. The presence of only one distinct population indicates that no cell aggregations are present in the sample, as these would appear lower on the plot.

We then investigated which constructs presented the display construct. We used the two fluorescent markers fluorescein-5-thiosemicarbazide (FTSC) attached to zymosan and AlexaFluor®633 as our secondary antibody respectively to measure the presence of the PAMP in proportion to the displayed protein. A cell control with both antibodies and 0.2 mg of zymosan served as the baseline, which the displays were compared to. A gate named “high display” was set (Fig. 19) covering the whole PAMP fluorescence range, with the left border of the antibody fluorescence range at the edge of the main population. This gate included all cells that showed interaction with the antibodies and excluded those that did not. The event count in this gate was subsequently analyzed and the constructs were compared to the cell control. A higher percentage of cells in this gate indicates higher binding of the secondary antibody, which is expected if our construct is displayed.

As seen in Fig. 19 the displays of the Vg_α‑helical, Vg_α‑helical+DUF1943, and Vg_truncated domains were successful. The cell control shows a small population of cells interacting with the antibody, leading to 8.75% of the events being in the set gate. We tried normalizing the measurements with the control, but the low signal-to-noise ratio of our samples led all measured signals to be quashed to zero, making them unusable. But even without normalization we could observe a clear increase in event count across the three constructs. As expected, the size of the domains had a great impact on display efficiency with the Vg_α‑helical domain showing the strongest display, at 11.75% more events inside the chosen gate, while the larger Vg_α‑helical+DUF1943 was displayed to a slightly lower degree with an 10.05% increase in events. To our delight, we were also able to observe the display of the 150 kDA truncated sequence at a sufficient level, at 4.65% over the cell control. Only the separate display of the Vg_DUF1943 domain could not be demonstrated as the FACS runs were unusable due to technical issues of the cytometer.

This figure shows graphics for different Vitellogenin domains.
A
This figure shows graphics for different Vitellogenin domains.
B
This figure shows graphics for different Vitellogenin domains.
C
This figure shows graphics for different Vitellogenin domains.
D
Figure 19: The plots shown compare AlexaFluor®633 fluorescence of the secondary antibody (DyLight633 channel) to the FTSC fluorescence of Zymosan (AlexaFluor®488 channel) of the A) double antibody cell control + 0.2 mg/mL zymosan, B) Vg_α‑helical + 0.2 mg/mL zymosan, C) Vg_ α‑helical+DUF1943 + 0.2 mg/mL zymosan and D) Vg_truncated + 0.2 mg/mL displays. Higher fluorescence of AlexaFluor®633 corresponds to higher protein display while higher FTSC fluorescence indicates higher zymosan binding. The gate ´high display´ was set at the right edge of the cell controls main population. Vg_α‑helical (B) shows the strongest secondary antibody fluorescence with 20.5% in the ´high display´ gate. Vg_α‑helical+DUF1943 (C) and Vg_truncated (D) show lower display efficiencies at 18.8% and 13.4% respectively. Only the cells in the high display gate were used for further affinity analysis.

Vg_DUF1943 Display Analysis

Before analysis we viewed every sample in SSC-A against time plot to see if any irregularities were present. Almost all runs showed perfect readings but the VG_DUF1943 run had severe dropouts and shifts in readings (Fig. 20). We tried to negate these problems via a more extensive gating protocol. We first used a gate called ´FACS failure´ to remove all events from the measurement that fell into the timeframe of the dropout. Unfortunately, when plotting these cells using FSC-A against SSC-A we still saw a very irregular event pattern. One diffuse population was visible at a similar position to the normal cells, with a second population at SSC-A approaching zero. We used a gate to isolate the events at the usual position for cells and proceeded to analyze these. When gating them with the same ´high display´ and ´high binding´ gates as other samples, we realized that they possessed a display value about half of the cell control and only showed baseline PAMP binding affinity. We were unable to run another measurement due to time constraints and were therefore unable to make a verdict on the display of the Vg_DUF1943 domain.

This figure shows graphics for different Vitellogenin domains.
A
This figure shows graphics for different Vitellogenin domains.
B
This figure shows graphics for different Vitellogenin domains.
C
This figure shows graphics for different Vitellogenin domains.
D
Figure 20: Irregular FACS run of VG_DUF1943. A) A long signal dropout is visible in the SSC‑A/time plot which was eliminated via the gate ´FACS failure´. B) The following FSC‑A/SSC‑A plot showed a very irregular event pattern with one diffuse cell population and a second population at a very low SSC‑A. The cell population was used and analysed as previously using the C) ´high display´ and D) ´high binding´ gates where the cells showed antibody binding below the level of the cell control and only baseline PAMP affinity.

Measurement of Domain PAMP Affinity

After verifying the correct display of our domains of interest we investigated which of them were involved in Vitellogenin’s PAMP binding affinity. A new gate was created in the “high display” gate called ´high binding´ (Fig. 21). It was set to include the top ~50% events of the highest binding domain. This range was chosen as it reduced the effect of both positive and negative outliers. A gate set too high would increase the share of positive outliers in the overall events, while a broader gate, more focused on lower binding domains, would include almost all events of strong binding domains, making these domains more difficult to characterize.

The cell control again showed some affinity which could not be compensated. We expected the Vg_αhelical domain to have decent binding affinity, as it was previously suggested to be the main factor for PAMP recognition [57], but we could not confirm this. With 19.2% of events in the set gate it was only marginally more than the cell control. We suspected potential misfolding to be responsible for its poor binding characteristics, but in silico analysis did not support this. The Vg_alpha+DUF1943 domain showed the strongest binding, with 52.7% inside its gate, strongly suggesting a role of the Vg_DUF1943 domain in PAMP binding. This was further supported by the measurements of Vg_truncated. It had just 35.3% of events in the ´high binding´ gate, indicating that the remaining protein does not, or only to a small degree, participate in zymosan PAMP recognition.

As we had no ability to investigate the Vg_DUF1943 affinity we were unable to determine if zymosan binding is only mediated by the Vg_DUF1943 domain, or if the Vg_α‑helical contributes to it as well. Based on earlier studies describing Vg_α‑helical PAMP recognition and the marginal binding affinity observed in our measurement we suggest that both domains are cooperatively binding the polymer.

A
B
C
D
Figure 21: The plots shown compare AlexaFluor®633 fluorescence of the secondary antibody (DyLight633 channel) to the FTSC fluorescence of Zymosan (AlexaFluor®488 channel) in the ´high display´ gate of the A) double antibody cell control + 200 µg/mL zymosan, B) Vg_α‑helical + 200  µg/mL zymosan, C) Vg_α‑helical+DUF1943 + 200 µg/mL zymosan and D) Vg_truncated + 200 µg/mL. Higher fluorescence of AlexaFluor®633 corresponds to higher protein display while higher FTSC fluorescence indicates higher zymosan binding. The gate ´high binding´ was set to include half the events of the strongest binding domain. This guaranteed that the focus on outliers of both strong and weak binding domains was minimized and that lower binding domains were included in the measurement to a sufficient degree. Vg_ αhelical+DUF1943 (C) showed the strongest binding with 52.7% of events inside the gate, binding stronger than Vg_truncated (D) at 35.3% and much stronger than Vg_ αhelical (B) at 19.2% which was previously suggested for PAMP affinity.

Based on these results, we chose Vg_alpha+DUF1943 for further work. It presented the best balance between display efficiency and strong binding affinity. Before screening the library, we needed to determine the binding kinetics of Vg_alpha+DUF1943. Without the correct PAMP concentration, we could either saturate the different protein variants regardless of their affinity or present the variants with too little PAMP to generate a sufficient signal to sort the cells. We measured a range of zymosan dilutions between 0.02 µg/mL and 200 µg/mL with the display and gated them as detailed above. We then used the event ratio to determine the degree of PAMP saturation of the presented domains. When fitting a curve to our data points assuming Michaelis‑Menten-like kinetics we calculated a KD = 0.61 µg/mL. As there is no substrate turnover with Vitellogenin KD is equal to KM (Fig. 22). Molar concentrations could not be determined, as we were unsure about the ratio of different zymosan oligomers in the labeled solution. Based on this we chose a concentration of 2 µg/mL for library screening as it presented a concentration which should saturate cells to ~75%. We expected that this would improve library sorting as binding cells would be strongly labeled and easier to isolate. Lower zymosan concentrations around KD or below would be more sensitive but could compromise sorting as binding cells are harder to differentiate from low or non‑binding ones.

This figure shows a curve of the zymosan concentration in µg/mL against the Vitellogenin zymosan saturation. The curve is exponentiel until around (10|0.9), then it runs against 1.
Figure 22: Kinetics fit to the binding affinity of Vg_α‑helical+DUF1943 surface display. Different zymosan concentrations between 0.2 µg/mL and 200 µg/mL were measured and normalized using the cell control. Based on this kinetics fit, the KD was determined to be 0.61 µg/ml.

We generated the library via error‑prone PCR using 0.1, 1, and 10 µM of 8‑o‑xoGTP and dPTP each. The fragments and a backbone amplified by a high-fidelity polymerase were used for the transformation of yeast. Subsequently, the cells were prepared for sorting. We decided to use two gates for sorting. One gate isolated cells with the combined properties of a high display with high binding affinity, while the other was set to sort out cells with slightly lower display and strong binding. These were chosen to not favor cells displaying variants with improved display characteristics, as these would be enriched in a high display gate. Due to time constraints, we were limited to a single FACS run and had to sort single cells immediately. Iterative FACS runs would have enabled us to sort out the top percent of binding cells and regrow them to a higher quality library, which would be sorted again. This would enrich improved variants to a point where they could be easily isolated with a lower risk of losing improved variants from a too aggressive sorting approach.

Unfortunately, when analyzing the data of the cell sorting it became apparent that the binding affinity of all libraries was not significantly different from the level of the cell control (Fig. 23).

This figure shows a curve of zymosan concentration against the percentage of cells in the high display gate.
Figure 23: Binding affinities of the Vg_α‑helical+DUF1943 libraries which were generated using 0.2, 2 and 20 µM of nucleotide analogues. All libraries were measured using 2 µg/mL and 20 µg/mL. The cell control also included a measurement at 10 µg/mL. The high display gate was set to include about half of the strongest binding variants. As evident no library showed increased PAMP affinity over the cell control which would indicate display of functional protein.

To investigate the source of our findings, we cloned fraction of each library into a pBSK vector and send them for sequencing. It was revealed that the error-prone PCR amplified a range of unspecific and gapped vector fragments (Fig. 24). This, of course, made the library unusable as the lack of the homologous ends, only present to the left and right border of Vg_α‑helical+DUF1943, made recombination of the fragments inside the yeast impossible. The reason for the poor performance of the primers was likely in the switch from the NEB Q5® polymerase, where they were very reliable, to a Taq polymerase required for error-prone PCR and the accompanying changes in annealing temperature. The new reaction conditions apparently led to unspecific primer annealing. In the PCR, smaller fragments are generated with a higher amplification efficiency, which leads to them enriching in the reaction. Thereby, the desired sequence is largely replaced.

This figure shows the alignment of the fragments amplified via error-prone PCR
Figure 24: Sequences of fragments amplified via error-prone PCR. The fragments were cloned into pBSK via blunt-end ligation and send for sequencing. It is evident, that the primers did not bind specifically, with only one fragment aligning to the intended Vg_αhelical+DUF1943. The frequent occurrence of gaps is also a severe problem for library generation.

Outlook

While our assay proved to be a reliable and easy-to‑perform, there are some adjustments that can be made to improve it for future applications. We would advise the addition of a terminator like PGK1 or CYC1 in the display constructs. It could increase protein yield and display capability of the cells by terminating transcription directly after the protein, reducing the metabolic burden of expression [58].

Another issue was the poor separation of displaying cells from non‑displaying ones in the cell staining. This appears to be a potential issue for large or flexible display constructs [39, 59]. The first approach for improving cell staining would be to utilize the second immune epitope to avoid potential steric hinderance of antibody binding at the c‑Myc tag. If this yields no improvement changes in the staining protocol could prove promising. Increased antibody incubation time or decreasing the incubation temperature present viable approaches [39]. Wen et al. (2010) also encountered poor staining separation issues similar to ours, which they resolved by using a single antibody stain instead of the standard two-step antibody protocol. We think these adjustments would greatly enhance the precision of the analysis and cell sorting which leads to an improved identification of the displayed cells while the carryover of non‑displaying cells is greatly improved.

For library generation we advise to check the error-prone PCR on a gel to see if unspecific fragments are present. We did not find sources suggesting that nucleotide analogues are responsible for unspecific primer binding or generation of gapped PCR sequences. If our issues reoccur with adjusted annealing temperature and primers it would be advisable to change the error-prone PCR to a variant employing manganese ions and adjusted dNTP mixes to counteract GC-bias.



Expression of Helicoverpa armigera C-type Lectin 14

As an alternative approach to increase pathogen associated molecular pattern (PAMP) binding affinity, we wanted to express the C-type lectin 14 (CTL14) of Helicoverpa armigera in Escherichia coli for testing in the zymosan binding assay. C-type lectins are a broad superfamily of proteins characterized by their C-type lectin domain (CTLD). A common feature of CTLs is calcium dependent binding of ligands [56]. These proteins play important roles in immunity. For example, they recognize microbes, such as fungi and bacteria, through PAMPs through various glycan structures [56]. CTLs are key factors in the immune system of invertebrates and most importantly for insects, such as Apis mellifera [57, 58]. We found out that the CTL14 of H. armigera can bind the PAMP zymosan [59]. Zymosan is also a major cell wall component of Beauveria bassiana, one of the entomopathogenic fungi which could be used to fight against the Varroa mite [60], thus making the CTL14 a suitable protein to test in the zymosan binding assay. Furthermore, we wanted to test a fusion protein of Apis mellifera Vitellogenin and CTL14. We hypothesized that a fusion protein would show improved zymosan binding activity, due to the binding capabilities of CTL14 and Vitellogenin. We started planning our lab work and looking into possible expression vectors. We identified the pET28b as a suitable expression vector for the expression of recombinant proteins. The pET28b has the advantage that expression can be induced through the lac-operon with isopropyl‑beta‑D‑thiogalactopyranoside (IPTG). Additionally, we needed a possible linker for the A. mellifera Vitellogenin and H. armigera CTL14 to construct a fusion protein. Our research was focused on possible linkers, which would keep the activity of both proteins intact and their structures distinct. We simulated the protein conformation of the fusion protein of Vitellogenin and CTL14 with AlphaFold using different linkers. In the end, we chose to express the fusion protein with a helical linker with the amino acid sequence A(EAAAK)2A. This linker was chosen as it was simulated to keep the proteins structures distinct and the His-Tag accessible for purification. More information regarding the fusion protein modeling is available on our modeling page. For the construction of the vector, we chose the E. coli Top10 strain for its high transformation rates and for the expression of the proteins we chose E. coli BL21, as it has reduced protein degrading and the T7 RNA-polymerase of the T7 phage under the control of a lac operon, making it ideal for high protein expression using the pET28b vector. For this project, we chose E. coli strains instead of Komagataella phaffii, as we did for the expression solely of the A. mellifera Vitellogenin. E. coli is often used in the literature for the expression of CTL´s [61]. Consequently, we aimed to investigate whether the interaction between Vitellogenin and CTL14 could lead to a substantial change in binding affinity and potentially enhance the bee's immune response positively. Thus, we opted to evaluate this combination using a fusion protein.

Having decided on the expression strain, linker, and vector for the expression, we started the construction of the pET28b-vector containing the CTL14, Vitellogenin and them linked together, so three different vectors. The proteins CTL14 and Vitellogenin were cloned solely into the pET28b vector to normalize the binding of zymosan by them. We decided to do Gibson assemblies for the construction of the expression vectors as the pET28b lacks the necessary restriction sites for golden gate assemblies.

pET28b vector graphics

Figure 25: pET28b_CTL14 vector map. Incorporated in the PET28b backbone was the Helicoverpa armigera ctl14 gene. Method used for the assembly of the expression vector was Gibson Assembly. BsaI restriction sites flanking the ctl14 gene were added for iGEM compatibility. A C-terminal 10x his-tag was added to the ctl14 for purification of the recombinant protein. For the expression the Escherichia coli BL21 strain was to be used.
Figure 26: pET28b_Vg vector map. Incorporated into the pET28b backbone was the Apis mellifera vitellogenin gene. Method used for the assembly of the expression vector was Gibson Assembly. BsaI restriction sites flanking the vitellogenin were added for iGEM compatibility. A C-terminal 10x his-tag was added to the vitellogenin for purification of the recombinant protein. For the expression of the Vitellogenin Escherichia coli BL21 was to be used.
Figure 27: pET28b_Vg-linker-CTL14 vector map. Incorporated into the pET28b backbone was the linker protein of Apis mellifera Vitellogenin and Helicoverpa armigera CTL14. The Vitellogenin and CTL14 were linked through a helical linker with the sequence A(EAAAK)2A. BsaI restriction sites flanking the linker protein were added for iGEM compatibility. The method used for the construction of the expression vector was Gibson Assembly. A C-terminal 10x his-tag was added to the ctl14 for purification of the recombinant linker protein. For the expression of the linker protein Escherichia coli BL21 was to be used.

Our plan was to add BsaI restriction sites flanking the genes during the construction of the expression vectors for RCF[1000]compatibility. The protein expression of the three recombinant proteins would then have been induced through the lac-operon located on the pET28b vector.

We began by preparing the vitellogenin and ctl14 genes, cloning a 10x His-tag onto the C-terminal ends of the genes by PCR amplification. This would allow us to purify the proteins later through affinity chromatography. Gel electrophoresis was used to verify the length of the DNA fragments and is shown in Fig. 28.

all relevant information in caption
Figure 28: Agarose gel of the PCR, which added the His10 tag to the Vitellogenin and CTL14 gene. The used ladder is the 1 kb Plus ladder from New England Biolabs. Nr.1 shows the Vitellogenin PCR and Nr.2 shows the CTL14 PCR. Both show the desired length of around 5,4 kb for Vitellogenin and 1 kb for CTL14.

all relevant information in caption
Figure 29: Agarose gel of the linearization of the pET28b vector. The used ladder is the 1 kb Plus ladder from New England Biolabs. Lane 1 shows the linearized vector at around 5,5 kb.

In the next step, the previously amplified fragments were appended to the pET28b vector using PCR with optimized primers to create the necessary overhangs for Gibson Assembly (Fig. 29). Initially, primers were used to create an overhang that included specific bases from the pET28b vector, ensuring that the fusion sides for the Gibson Assembly were located in the pET28b vector and not in the genes. For this purpose, the pET28b empty vector was linearized at the desired locations by PCR, with the ends containing the necessary overhangs for Gibson Assembly.

Following the amplification, a DpnI digestion was done to remove the backbone from the reaction mixture and increase the cloning efficiency. DpnI is a restriction enzyme cutting only methylated DNA, which is only present in the plasmid isolated from E. coli and not the PCR product. Subsequently, two Gibson assemblies were performed for the linearized vectors pET28b_Vg and pET28b_CTL14, followed by a transformation into E. coli Top10. Since the pET28b vector contains an antibiotic resistance gene for chloramphenicol, selection was performed on LB agar plates supplemented with chloramphenicol. The grown colonies were checked for the presence of the desired genes using a colony PCR. No positive colonies were found during screening.

all relevant information in caption
Figure 30: Agarose gel of the pET28b vector linearized with the new primers for Vitellogenin. The used ladder is the 1 kb Plus ladder from New England Biolabs. Lane 1 shows the linearized pET28b_Vg vector after DpnI Digestion at around 5,5 kb.

In consultation with our instructors, we decided to subject the pET28b vector to extended DpnI digestion overnight. Additionally, we optimized the Gibson assembly overhang primers to not only extend the overhang in the pET28b vector but also extend the vector into the genes. These overhangs were designed to be 20 to 44 bases in length. Unfortunately, the results with the extended Gibson assembly overhang primers were unsatisfactory, prompting us to attempt gradient PCR and touchdown PCR to create optimal conditions for the primers. After many optimization steps, we were able to successfully complete the PCRs and amplify the genes with the necessary overhangs we needed for cloning them into the pET28b vector. One result is shown in Fig. 30.


At this point, we shifted our focus on the project involving a vector in which Vitellogenin is attached to CTL14 using a protein linker. A His-tag was first cloned onto the C-terminal end of the CTL14 gene through PCR. Subsequently, optimized primers were used again to clone the protein linker onto the C‑terminal end of the vitellogenin gene. Due to the linker between the Vitellogenin and CTL14, different Gibson overhangs were cloned via PCR. We used the same linearized vector from the pET28b_Vg plasmid for the Gibson Assembly. The Gibson Assembly for the pET28b_Vg_alphahelical_CTL14His10xC was then carried out and transformed into E. coli Top10. Unfortunately, no positive colonies were observed in subsequent colony PCRs. We repeated the colony PCRs using optimized primers for the incorporated genes, but still, no positive colonies could be detected.

The picture shows the four steps in which the CTL14 project workflow was
  split. The first step was gathering information and therefore discovering the
  CTL14 protein of <em>Helicoverpa armigera a</em>nd a useful expression vector.
  The second step was the designing of all the laboratory work. The biggest part
  was the simulations of suitable linkers for the fusion protein. In the third
  step the designed vectors were constructed and controled. The last step
  contains the expression of the fusion protein in <em>E. coli</em> BL21.
Figure 31: Workflow for the expression of the Apis mellifera Vitellogenin and Helicoverpa armigera CTL14 linker protein.

Despite encountering similar problems, we diligently attempted to troubleshoot and rectify the issues. However, we regret to report that we were not successful in assembling the desired vector, and further investigation was impeded by time constraints. Therefore, we reached a point where we were unable to explore additional potential sources of error due to the limited available time for our research.



Expression of Apis mellifera Vitellogenin

The figure shows a flow chart of the workflow for the expression of Vitellogenin. 1. gathering information, 2. planning the
lab work, 3. contructing the vectors, 4. expression in
Komagataella phaffii
Figure 32: Workflow for the expression of Apis mellifera Vitellogenin in Komagataella phaffii CBS7435 PEP4 KO.

The aim of our project was to express the optimized Apis mellifera Vitellogenin, which was created either through directed evolution or modeled in the dry lab via rational design choices, for further analysis and characterization. The expression system used was the yeast Komagataella phaffii, formerly Pichia pastoris.

Our research into the A. mellifera Vitellogenin began after we had decided to tackle the problem Varroa mites and other pests pose to Western honey bees. We became aware that Vitellogenin is glycosylated under native conditions, thus directing our research towards an organism capable of glycosylation [1]. As a eukaryote, Komagataella phaffii is capable of glycosylation and is often used for the expression of glycosylated proteins instead of Escherichia coli, as glycosylation does not occur in the latter [62–66]. For this reason, we chose K. phaffii as the organism for expression instead of E. coli. We acquired the K. phaffii CBS7435 PEP4 KO strain in cooperation with the Institute for Microbiology and Biotechnology from the University of Vienna. The strain exhibits reduced protein degradation, as the knocked-out pep4 gene is the precursor of a major protease.

The figure shows the
<em>Saccharomyces cerevisiae</em> shuttle vector pSB3CY_His. The plasmid
contains the <em>S. cerevisiae</em> origin of replication 2µ. The plasmids
contains the a selction marker for histidine.
Figure 33: pSB3CY_His (BBa_K4188002) plasmid map. This level two shuttle vector, also called pSB3CY_His, was designed by the previous iGEM Münster 2022 team. BBa_K4188002 was used as the backbone during the construction of the expression vectors and the 2µ origin of replication and His3 marker were interchanged for usage in Komagataella phaffii.

Next, we investigated promoters and terminators suitable for expression in K. phaffii. The AOX1 promoter, a methanol inducible promoter, and its corresponding terminator are used for their high expression rates in K. phaffii [66, 67]. For our project, we decided to build upon the work of the iGEM Münster team 2022 and use their previously designed vector pSB3CY_His (BBa_K4188002) for S. cerevisiae (Fig. 33). The pSB3CY_His has a origin of replication and a his3 auxotrophy marker. We replaced both parts with parts commonly used in literature for recombinant protein expression in K. phaffii. Most often, positive transformants are selected by introducing either an auxotrophy marker or an antibiotic resistance to the plasmid. The selection of positive transformants in the project was done by the antibiotic Zeocin®. Many expression vectors for K. phaffii carry PARS for replication. Based on the widespread usage of both Zeocin® for selection and PARS as the origin of replication, these parts were integrated into our vector to generate the pSB3CYP_Zeo plasmid (BBa_K4675009) (Fig. 34)­. Apis mellifera’s vitellogenin gene was ordered from IDT in 4 fragments and assembled via Gibson Assembly (Fig. 35).

The figure shows a Gel electrophoresis of the PCR
    amplifying the PARS and zeoR in the constructed pSB3CYP_Zeo. A band at
    approximately 700 bp was expected for the amplifications. The analysis
    showcased that the parts were integrated into the vector.
Figure 34: Gel electrophoresis of the PCR amplifying the PARS and zeoR in the constructed pSB3CYP_Zeo (BBa_K4675009). A band at approximately 700 bp was expected for the amplifications. The analysis showcased that the parts were integrated into the vector.
The figure shwos a Gel electrophoresis analyzing
    the Gibson Assembly of the Apis mellifera vitellogenin. The gene was assembled
    from four fragments synthesized by IDT. A band at 5 kb was expected for the
    assembled gene. The correct band was extracted from the gel.
Figure 35: Gel electrophoresis analyzing the Gibson Assembly of the Apis mellifera vitellogenin. The gene was assembled from four fragments synthesized by IDT. A band at 5 kb was expected for the assembled gene. The band was extracted from the gel.

To begin with the construction of the pSB3CYP_VgHis10xC, we designed primers to amplify the PARS, zeoR, promoter AOX1 and terminator AOX1 from the Open Yeast Library. Additionally, PCRs were carried out amplifying the A. mellifera vitellogenin. The E. coli Turbo strain was used to assemble and amplify our newly constructed plasmids. The primers were designed to add SapI restriction sites and fusion sites to the promotor AOX1, the vitellogenin gene and the terminator AOX1 for the assembly into the vector, replacing the mRFP gene. Additionally, the primer designed for vitellogenin added a C-terminal 10x His-Tag. The zeoR was amplified with primers adding NotI and SpeI restriction sites and PARS with primers adding SpeI and BamHI restriction sites. Before PARS could be integrated into the pSB3CY_His shuttle vector (BBa_K4188002) a site-directed mutagenesis had to be performed, deleting a BamHI restriction site in the PARS origin by changing one base of the sequence. The PCRs were done using high fidelity polymerase to minimize the introduction of additional mutations into the DNA sequences. The PCRs amplifying the parts from the distribution kit did not work and optimization of the primers yielded no positive results (Fig. 36). Additionally, we tested the plasmids in a gel electrophoresis for length, to check for dimerization or multimerization. The gel electrophoresis showcased that nearly all tested plasmid multimerized during the incubation (Fig. 37). Hence, we decided to order the genes present in the distribution kit from IDT after discussing the issue with our instructors and advisors. Using the genes synthesized by IDT as a template, the fragments containing restriction sites and fusion sites were successfully amplified. The zeoR and PARS were integrated into the vector by restriction cloning, replacing the original and his3. The next step was to exchange the mRFP with the cassette for expression of Vitellogenin. The Golden Gate assemblies (RCF[1000]) using the SapI restriction enzymes did not work using our protocol, as SapI did not seem to cut in the vector. Here again we saw that dimerization occurred in E. coli Turbo. Shortly here after we switched to E. coli Top10 (Fig. 38). The PCRs for the amplification of the parts for the expression cassettes were repeated and the protocols checked by our advisors. Still no successful construction could be achieved. We then ordered the pSB3CYP_VgHis10xC vector as fragments and assembled the vector via Gibson Assembly. We divided the vector into seven fragments and added the overhangs necessary for the Gibson Assembly directly to the ordered sequences (Fig. 39). However, sequencing of the construct showcased that the stop codon in the zeoR gene had undergone a mutation, making the gene non-functional and selection of transformed K. phaffii impossible. Due to time constraints, we were unable to conduct site-directed mutagenesis in the stop codon of zeoR.

The figure shows a gel electrophoresis of the
    PCR amplification for the PARS and ZeoR genes from the plasmids of the
    distribution kit. For the amplification of the PARS a band at approximately
    200 bp was expected and for zeoR amplification a band at approximately 400 bp.
    In the zeoR amplifications one band appeared at 500 bp, which could correspond
    to the PCR product. For the PARS amplifications, no bands were detected.
Figure 36: Gel electrophoresis of the PCR amplification for the PARS and ZeoR genes from the plasmids BBF10K 000465 and BBF10K_000475, respectively. For the amplification of the PARS a band at approximately 200 bp was expected and for zeoR amplification a band at approximately 400  bp. In the zeoR amplifications one band appeared at 500 bp, which could correspond to the PCR product. For the PARS amplifications, no bands were detected. This gel electrophoresis is showcased as an overview for all amplifications of the parts used in the Open Yeast Library.
The figure shows a Gel electrophoresis checking
    the plasmids from the Open Yeast Library. In the first lane is the pAOX1
    plasmid, in the second lane is the aMF plasmid, in the third lane is the PARS
    plasmid, in the forth lane is the tAOX1 plasmid, in the fifth lane is the zeoR
    plasmid, in the sixth lane is the zeoR plasmid. Bands were visible in samples
    1;4;5 and 6 at approximately 20 kb. The tested plasmids should be between 2000
    bp and 3000 bp in length.
Figure 37: Gel electrophoresis checking the plasmids from the Open Yeast Library. 1= pAOX1 plasmid BBF10K_000181, 2= aMF plasmid BBF10K_000355, 3= PARS plasmid BBF19K_000465, 4= tAOX1 plasmid BBF10K_000418, 5= zeoR plasmid BBF10K_000475, 6= zeoR plasmid BBF10K_000475. Bands were visible in samples 1;4;5 and 6 at approximately 20 kb. The tested plasmids are between 2000 bp and 3000 bp in length.
The figure shows a Gel electrophoresis of the
      SapI test digestion using the pSB3CY_His vector. The test digestion was done
      to check for cutting of the vector by SapI. The first lane is the reaction
      mixture containing the plasmid without SapI, the second lane is the reaction
      mixture containing the plasmid and SapI,the third lane is the negative control
      without plasmid. The pSB3CY_His has a length of 6.3 kb. All samples containing
      the plasmids showed a band corresponding to the plasmid, but also a band at 20
      kb. No bands were detected corresponding to fragments made during the
      restriction digestion. SapI does not seem to cut in the vector.
Figure 38: Gel electrophoresis of the SapI test digestion using the pSB3CY_His vector (BBa_K4188002). The test digestion was done to check for cutting of the vector by SapI. 1= reaction mixture containing the plasmid without SapI, 2= reaction mixture containing the plasmid and SapI, 3= negative control without plasmid. The pSB3CY_His (BBa_K4188002) has a length of 6.3 kb. All samples containing the plasmids showed a band corresponding to the plasmid, but also a band at 20 kb. No ba  detected corresponding to fragments made during the restriction digestion. SapI does not seem to cut in the vector.
The figure shows the
      Komagataella phaffii shuttle vector pSB3CYP_VgHis10xC. The plasmid
      contains the gene for Vitellogenin and a 10x HisTag at the Vitellogenin
      C-terminus.
Figure 39: pSB3CYP_VgHis10xC plasmid map. iGEM part BBa_K4675009 was used as the backbone for the construction of the expression vector. Method for the assembly of the expression vector was Golden Gate Assembly (RCF[1000]). A C-terminal 10x his-tag was added to the Vitellogenin for purification of the recombinant protein. For the expression, Komagataella phaffii CBS7450 PEP4 KO was to be used. PARS origin of replication and Zeocin® resistance (zeoR) were used for selection of positive Komagataella phaffii colonies.


References

  1. [1] V. Leipart et al., ‘Structure prediction of honey bee vitellogenin: a multi-domain protein important for insect immunity.’, FEBS Open Bio, vol. 12, no. 1, pp. 51–70, Jan. 2022, doi: 10.1002/2211-5463.13316
  2. [2] J. D. Bloom and F. H. Arnold, ‘In the light of directed evolution: pathways of adaptive protein evolution.’, Proc Natl Acad Sci U S A, vol. 106 Suppl 1, no. Suppl 1, pp. 9995–10000, Jun. 2009, doi: 10.1073/pnas.0901522106
  3. [3] S.-M. Lee, T. Jellison, and H. S. Alper, ‘Directed evolution of xylose isomerase for improved xylose catabolism and fermentation in the yeast Saccharomyces cerevisiae.’, Appl Environ Microbiol, vol. 78, no. 16, pp. 5708–16, Aug. 2012, doi: 10.1128/AEM.01419-12
  4. [4] A. L. Betteridge, K. M. Sumby, J. F. Sundstrom, P. R. Grbin, and V. Jiranek, ‘Application of directed evolution to develop ethanol tolerant Oenococcus oeni for more efficient malolactic fermentation.’, Appl Microbiol Biotechnol, vol. 102, no. 2, pp. 921–932, Jan. 2018, doi: 10.1007/s00253-017-8593-x
  5. [5] A. M. Guindal, R. Gonzalez, J. Tronchoni, J. S. Roodink, and P. Morales, ‘Directed evolution of Saccharomyces cerevisiae for low volatile acidity during winemaking under aerobic conditions.’, Food Microbiol, vol. 114, p. 104282, Sep. 2023, doi: 10.1016/j.fm.2023.104282
  6. [6] J. Jiang, K. M. Sumby, J. F. Sundstrom, P. R. Grbin, and V. Jiranek, ‘Directed evolution of Oenococcus oeni strains for more efficient malolactic fermentation in a multi-stressor wine environment.’, Food Microbiol, vol. 73, pp. 150–159, Aug. 2018, doi: 10.1016/j.fm.2018.01.005
  7. [7] E. B. Cahoon and J. Shanklin, ‘Substrate-dependent mutant complementation to select fatty acid desaturase variants for metabolic engineering of plant seed oils.’, Proc Natl Acad Sci U S A, vol. 97, no. 22, pp. 12350–5, Oct. 2000, doi: 10.1073/pnas.210276297
  8. [8] S. Li et al., ‘Directed evolution of TetR for constructing sensitive and broad-spectrum tetracycline antibiotics whole-cell biosensor.’, J Hazard Mater, vol. 460, p. 132311, Oct. 2023, doi: 10.1016/j.jhazmat.2023.132311
  9. [9] C. Furusawa, T. Horinouchi, and T. Maeda, ‘Toward prediction and control of antibiotic-resistance evolution.’, Curr Opin Biotechnol, vol. 54, pp. 45–49, Dec. 2018, doi: 10.1016/j.copbio.2018.01.026
  10. [10] Á. Nyerges et al., ‘Directed evolution of multiple genomic loci allows the prediction of antibiotic resistance.’, Proc Natl Acad Sci U S A, vol. 115, no. 25, pp. E5726–E5735, Jun. 2018, doi: 10.1073/pnas.1801646115
  11. [11] E. T. Farinas, T. Bulter, and F. H. Arnold, ‘Directed enzyme evolution.’, Curr Opin Biotechnol, vol. 12, no. 6, pp. 545–51, Dec. 2001, doi: 10.1016/s0958-1669(01)00261-0
  12. [12] S. R. Nirantar, ‘Directed Evolution Methods for Enzyme Engineering.’, Molecules, vol. 26, no. 18, Sep. 2021, doi: 10.3390/molecules26185599
  13. [13] D. Sriprapundh, C. Vieille, and J. G. Zeikus, ‘Directed evolution of Thermotoga neapolitana xylose isomerase: high activity on glucose at low temperature and low pH.’, Protein Eng, vol. 16, no. 9, pp. 683–90, Sep. 2003, doi: 10.1093/protein/gzg082
  14. [14] M. W. Traxlmayr and E. V Shusta, ‘Directed Evolution of Protein Thermal Stability Using Yeast Surface Display.’, Methods Mol Biol, vol. 1575, pp. 45–65, 2017, doi: 10.1007/978-1-4939-6857-2_4
  15. [15] C.-Q. Zhong et al., ‘Improvement of low-temperature caseinolytic activity of a thermophilic subtilase by directed evolution and site-directed mutagenesis.’, Biotechnol Bioeng, vol. 104, no. 5, pp. 862–70, Dec. 2009, doi: 10.1002/bit.22473
  16. [16] B. van Loo et al., ‘High-Throughput, Lysis-Free Screening for Sulfatase Activity Using Escherichia coli Autodisplayed Arylsulfatase.’, ACS Synth Biol, vol. 8, no. 3, pp. 516–525, Mar. 2019, doi: 10.1021/acssynbio.8b00373.
  17. [17] S. Wang et al., ‘Directed evolution of xylose isomerase for improved xylose catabolism and ethanol production in the yeast Saccharomyces cerevisiae.’, Appl Environ Microbiol, vol. 80, no. 12, pp. 3813–21, Jun. 2014, doi: 10.1128/AEM.00401-14
  18. [18] T. S. Lee, A. J. Krupa, and C. A. Zhang, ‘Characterization of Nucleotide Distributions in Attenuator Stem-Loops Using a High-Throughput Assay.’, Biophys J, vol. 111, no. 2, pp. 356–368, Jul. 2016, doi: 10.1016/j.bpj.2016.06.018.
  19. [19] M. W. Traxlmayr, T. Hasenhindl, and A. Hackl, ‘Stable Variants of Human Carbonic Anhydrase II for the Application in a High-Temperature Adsorption Process.’, Biotechnol J, vol. 14, no. 10, Oct. 2019, doi: 10.1002/biot.201800166.
  20. [20] J. S. Lee et al., ‘Directed evolution of AraC for improved compatibility of arabinose- and lactose-inducible promoters.’, Appl Environ Microbiol, vol. 81, no. 18, pp. 6318–27, Sep. 2015, doi: 10.1128/AEM.01180-15.
  21. [21] M. W. Traxlmayr et al., ‘Directed Evolution of a Secretory Leader for the Improved Expression of Heterologous Proteins and Full-Length Antibodies in Yarrowia lipolytica.’, Appl Environ Microbiol, vol. 81, no. 23, pp. 7108–18, Dec. 2015, doi: 10.1128/AEM.02105-15.
  22. [22] J. B. Matulef, S. Park, and W. N. DeGrado, ‘Probing Protein Structure and Function in Real Time and in Cells with All-Maleimide Photoswitchable Mimics of Green Fluorescent Protein.’, J Am Chem Soc, vol. 142, no. 1, pp. 214–221, Jan. 2020, doi: 10.1021/jacs.9b11215.
  23. [23] N. Steiner, S. Dootz, and H. F. Noller, ‘Directed evolution of ribosomes and the prospects for synthetic biology.’, FEBS Lett, vol. 592, no. 11, pp. 2104–2112, Jun. 2018, doi: 10.1002/1873-3468.13068.
  24. [24] T. D. Newby, J. L. Hollinshead, R. G. Bhandari, and S. Dennison, ‘A Directed Evolution Approach to Select for Improved Bacterial Proteins.’, PLoS One, vol. 11, no. 4, p. e0153875, 2016, doi: 10.1371/journal.pone.0153875.
  25. [25] J. Schilling, C. M. Schöppe, and S. A. Plückthun, ‘Co-evolution of Affinity and Stability of GxxxG-Like Pairs in Transmembrane Helices of Oligomeric Proteins.’, Biochemistry, vol. 55, no. 48, pp. 6634–6648, Dec. 2016, doi: 10.1021/acs.biochem.6b00943.
  26. [26] Y. Li, H. Mo, and X. Guo, ‘Directed Evolution of an α-Fucosidase from a GH29 Family Member for Improved Catalytic Activity.’, Molecules, vol. 26, no. 6, Mar. 2021, doi: 10.3390/molecules26061444.
  27. [27] J. T. M. Kennaway, G. S. He, D. A. Khan, and D. J. Scott, ‘On the Design of an Improved Sulfotransferase Assay with Broad Substrate Specificity.’, SLAS Discov, vol. 22, no. 1, pp. 59–66, Jan. 2017, doi: 10.1177/2472555216662066.
  28. [28] T. S. Bzymek, M. F. Hsieh, and J. G. Zahler, ‘How the Human Telomerase RNA Template-DNA Hybrid and POT1-TPP1 Promote Telomerase Primer-Template Recognition.’, J Am Chem Soc, vol. 141, no. 49, pp. 19556–19564, Dec. 2019, doi: 10.1021/jacs.9b08229.
  29. [29] T. R. Abrahamson et al., ‘A Novel Amber Allele Reveals a Key Role for Bacillus subtilis Ribonuclease III in Organismal RNA Processing.’, Mol Cell Biol, vol. 39, no. 15, Jul. 2019, doi: 10.1128/MCB.00056-19.
  30. [30] A. T. Dorrity, D. C. Zhu, M. F. Rock, and D. A. Cropek, ‘High-Throughput Screening of Bacterial Transformants Expressing Biomass-Degrading Enzymes.’, SLAS Discov, vol. 22, no. 1, pp. 11–21, Jan. 2017, doi: 10.1177/2472555216674074.
  31. [31] K. L. Stanger, A. B. Knecht, and M. R. Pelton, ‘Evaluation of a New Biosensor for the Detection of Organophosphorus Compounds.’, Anal Chem, vol. 93, no. 6, pp. 3022–3031, Mar. 2021, doi: 10.1021/acs.analchem.0c04476
  32. [32] N. F. Polizzi and M. A. DeLisa, ‘Exploiting Outer Membrane Vesicles in Native-Purpose Display Reveals a General Strategy for Developing Bacterial Vaccines.’, ACS Central Sci, vol. 6, no. 8, pp. 1343–1352, Aug. 2020, doi: 10.1021/acscentsci.0c00709
  33. [33] J. L. Hollinshead, T. D. Newby, R. G. Bhandari, and S. Dennison, ‘Directed Evolution of a Genetic Switch at the RNA Domain Level.’, Biotechnol J, vol. 13, no. 6, pp. 1700062, Jun. 2018, doi: 10.1002/biot.201700062
  34. [34] Y. He et al., ‘Directed Evolution of the 3-Hydroxypropionic Acid Production Pathway by Engineering Aldehyde Dehydrogenase Using a Key Metabolite Intermediate.’, Front Bioeng Biotechnol, vol. 7, p. 372, 2019, doi: 10.3389/fbioe.2019.00372
  35. [35] L. M. Molina-Espeja, T. Gullfot, and P. Christakopoulos, ‘A novel bacterial GH99 endo-α-1,2-mannanase from Pediococcus acidilactici that catalyses the hydrolysis of β-1,4-mannosidic linkages.’, Biochim Biophys Acta Proteins Proteomics, vol. 1866, no. 12, p. 140382, Dec. 2018, doi: 10.1016/j.bbapap.2018.08.008
  36. [36] Y. G. Kim et al., ‘Directed Evolution of Recombination Hotspots in the Lactococcal 936 Type Phage P2c.’, Front Microbiol, vol. 9, p. 2013, 2018, doi: 10.3389/fmicb.2018.02013
  37. [37] S. Zhang et al., ‘Directed Evolution of Ligninolytic Oxidoreductases: Engineering a Pycnoporus cinnabarinus Laccase for Pleurotus eryngii Dye Decolorization.’, Sci Rep, vol. 6, p. 24329, Apr. 2016, doi: 10.1038/srep24329
  38. [38] C. Xue et al., ‘Directed evolution of the cyanobacterial bicarbonate transporter BicA for high-affinity CO2 uptake.’, eLife, vol. 8, Apr. 2019, doi: 10.7554/eLife.45043
  39. [39] E. M. Harstad and S. Swanson, ‘Immunization of Rabbits for the Production of Antibodies.’, J Vis Exp, no. 63, p. 374, Mar. 2012, doi: 10.3791/374
  40. [40] Y. Zhao, H. K. Huang, and Y. Zhang, ‘Engineering lipoic acid production in Saccharomyces cerevisiae by rewiring the sulfur assimilation pathway and activating lipoic acid synthase.’, Metab Eng, vol. 34, pp. 1–12, Nov. 2015, doi: 10.1016/j.ymben.2015.12.006
  41. [41] E. Zelcbuch, R. Litsanov, and M. V. Siewers, ‘Cloning- and PCR-Based Methods for Gene Synthesis and Protein Engineering.’, Methods Mol Biol, vol. 1116, pp. 287–312, 2014, doi: 10.1007/978-1-62703-764-8_18
  42. [42] C. M. Li, P. S. Tan, and C. S. Lin, ‘A novel strategy for high-level production of heterologous proteins in Escherichia coli: the case of nucleoside diphosphate kinase.’, PLoS One, vol. 8, no. 2, p. e56472, 2013, doi: 10.1371/journal.pone.0056472
  43. [43] J. A. Wilson, A. L. Kim, and S. M. Carrington, ‘Discovery and Characterization of Protein Mutations that Alter the Phospholipase Activity of the Type III Effector AvrRpm1 from Pseudomonas syringae pv. Phaseolicola.’, PLoS One, vol. 10, no. 9, p. e0137189, 2015, doi: 10.1371/journal.pone.0137189
  44. [44] S. S. Lee et al., ‘Directed Evolution of Myocbacterium neoaurum Lipase for Enhanced Hydrolysis of Long-Chain Fatty Acids.’, J Agric Food Chem, vol. 63, no. 30, pp. 6796–6802, Aug. 2015, doi: 10.1021/acs.jafc.5b02426
  45. [45] A. A. Katrahalli et al., ‘Directed evolution of the TALE N-terminal domain for recognition of all 5’ bases.’, Nucleic Acids Res, vol. 42, no. 6, pp. 3939–3947, Apr. 2014, doi: 10.1093/nar/gku097
  46. [46] A. L. Griffiths and N. M. Anderson, ‘Challenges and Opportunities for Improving the Chick Embryo Chorioallantoic Membrane (CAM) Assay as a Model for Screening Anti-Angiogenic Agents.’, Cancers (Basel), vol. 11, no. 2, Feb. 2019, doi: 10.3390/cancers11020218
  47. [47] A. J. Waldo et al., ‘Improving protein crystal quality by selective cryo-microdissection of X-ray-induced microcrystals at a synchrotron light source.’, J Synchrotron Radiat, vol. 25, no. Pt 4, pp. 893–900, Jul. 2018, doi: 10.1107/S1600577518007914
  48. [48] R. F. H. J. Cox et al., ‘Rapid Analysis and Quantification of Nucleic Acid Samples Using UV Spectroscopy.’, Methods Mol Biol, vol. 1975, pp. 13–28, 2019, doi: 10.1007/978-1-4939-9220-1_2.
  49. [49] T. V. Nguyen and T. M. Gillis, ‘Site-saturation mutagenesis of residues at the α1–β2 interface reveals a critical role for Leu65 in human hemoglobin.’, Proteins, vol. 82, no. 1, pp. 13–29, Jan. 2014, doi: 10.1002/prot.24367.
  50. [50] M. S. Khan et al., ‘Engineering of Bacillus pumilus alkaline serine protease (AprBp) through immobilization technique for improved stability and activity.’, Int J Biol Macromol, vol. 176, pp. 120–129, May 2021, doi: 10.1016/j.ijbiomac.2021.01.017.
  51. [51] A. Kapoor, R. K. Sani, and D. S. Kumar, ‘Gene cloning, expression and characterization of a novel organic solvent stable lipase from Pseudomonas plecoglossicida JCM 10502.’, Microb Cell Fact, vol. 10, p. 33, Apr. 2011, doi: 10.1186/1475-2859-10-33.
  52. [52] K. Y. Santhi et al., ‘Site-Directed Mutagenesis of Yarrowia lipolytica Lipase Lip2 to Improve Its Thermostability.’, Appl Biochem Biotechnol, vol. 189, no. 4, pp. 1290–1302, Dec. 2019, doi: 10.1007/s12010-019-03042-x.
  53. [53] J. D. Roote, A. M. Lu, and J. B. Beckingham, ‘Tandem duplications of a 3.5 kb genomic fragment from the white locus are not found in white mutants.’, Fly (Austin), vol. 2, no. 6, pp. 310–312, Nov. 2008, doi: 10.4161/fly.7226.
  54. [54] P. Tyagi et al., ‘Engineering of lipase B from Candida antarctica (CalB) for the highly enantioselective Michael addition of 1,3-dicarbonyl compounds.’, Chem Commun (Camb), vol. 55, no. 2, pp. 196–199, Jan. 2019, doi: 10.1039/c8cc08474g.
  55. [55] A. L. B. Ambra et al., ‘Rational design of an ultrasensitive quorum-sensing switch.’, ACS Synth Biol, vol. 4, no. 5, pp. 412–419, May 2015, doi: 10.1021/sb500206d.
  56. [56] G. D. Brown, J. A. Willment, and L. Whitehead, “C-type lectins in immunity and homeostasis,” Nature reviews. Immunology, vol. 18, no. 6, pp. 374–389, 2018, doi: 10.1038/s41577-018-0004-8 Titel anhand dieser DOI in Citavi-Projekt übernehmen.
  57. [57] B. Pees, W. Yang, A. Zárate-Potes, H. Schulenburg, and K. Dierking, “High Innate Immune Specificity through Diversified C-Type Lectin-Like Domain Proteins in Invertebrates,” Journal of innate immunity, vol. 8, no. 2, pp. 129–142, 2016, doi: 10.1159/000441475.
  58. [58] Y. Lu et al., “Comparative genomic analysis of C-type lectin-domain genes in seven holometabolous insect species,” Insect biochemistry and molecular biology, vol. 126, p. 103451, 2020, doi: 10.1016/j.ibmb.2020.103451.
  59. [59] Y. Cheng, Z. Lin, J.-M. Wang, L.-S. Xing, G.-H. Xiong, and Z. Zou, “CTL14, a recognition receptor induced in late stage larvae, modulates anti-fungal immunity in cotton bollworm Helicoverpa armigera,” Developmental and comparative immunology, vol. 84, pp. 142–152, 2018, doi: 10.1016/j.dci.2018.02.010.
  60. [60] M. Lamkanfi, R. K. S. Malireddi, and T.-D. Kanneganti, “Fungal zymosan and mannan activate the cryopyrin inflammasome,” The Journal of biological chemistry, vol. 284, no. 31, pp. 20574–20581, 2009, doi: 10.1074/jbc.M109.023689.
  61. [61] D. Martínez-Alarcón, A. Blanco-Labra, and T. García-Gasca, “Expression of Lectins in Heterologous Systems,” International journal of molecular sciences, vol. 19, no. 2, 2018, doi: 10.3390/ijms19020616.
  62. [62] S. K. Prabhu, Q. Yang, X. Tong, and L.-X. Wang, “Exploring a combined Escherichia coli-based glycosylation and in vitro transglycosylation approach for expression of glycosylated interferon alpha,” Bioorganic & medicinal chemistry, vol. 33, p. 116037, 2021, doi: 10.1016/j.bmc.2021.116037.
  63. [63] V. Juturu and J. C. Wu, “Heterologous Protein Expression in Pichia pastoris: Latest Research Progress and Applications,” Chembiochem : a European journal of chemical biology, vol. 19, no. 1, pp. 7–21, 2018, doi: 10.1002/cbic.201700460.
  64. [64] G. D. Barone et al., “Industrial Production of Proteins with Pichia pastoris-Komagataella phaffii,” Biomolecules, vol. 13, no. 3, 2023, doi: 10.3390/biom13030441.
  65. [65] M. Ahmad, M. Hirz, H. Pichler, and H. Schwab, “Protein expression in Pichia pastoris: recent achievements and perspectives for heterologous protein production,” Applied microbiology and biotechnology, vol. 98, no. 12, pp. 5301–5317, 2014, doi: 10.1007/s00253-014-5732-5.
  66. [66] A. Shrivastava, M. Pal, and R. K. Sharma, “Pichia as yeast cell factory for production of industrially important bio-products: Current trends, challenges, and future prospects,” Journal of Bioresources and Bioproducts, vol. 8, no. 2, pp. 108–124, 2023, doi: 10.1016/j.jobab.2023.01.007.
  67. [67] S. Macauley-Patrick, M. L. Fazenda, B. McNeil, and L. M. Harvey, “Heterologous protein production using the Pichia pastoris expression system,” Yeast, vol. 22, no. 4, pp. 249–270, 2005, doi: 10.1002/yea.1208.