Materials and reagents

Biological materials

Trypanosoma brucei 427 strain bloodstream forms or conditional null for the gene PIP5Pase derived from the 427 strain (Cestari and Stuart, 2015)

Reagents

M-MuLV reverse transcriptase (New England Biolabs Ltd, catalog number: M02535)
M-MuLV reverse transcriptase 10× buffer (New England Biolabs Ltd, catalog number: B02535)
MgCl₂ (Invitrogen, catalog number: R0971)
NucleoMag NGS Clean-up and Size Selection beads (Takara, catalog number: 744970.5)
NEBNext end repair module (New England Biolabs Ltd, catalog number: E6050S)
Native Barcoding Expansion 1-12 (Oxford Nanopore Technologies, catalog number: EXP-PBC001)
10 M Sodium hydroxide solution (NaOH) (Sigma, catalog number: 1310-73-2)
NEBNext Quick Ligation Module (New England Biolabs Ltd, catalog number: E6056S)
NEBNext FFPE DNA Repair Mix 24 reactions (New England Biolabs Ltd, catalog number: M6630S)
Taq DNA polymerase with ThermoPol buffer (New England Biolabs Ltd, catalog number: M0267S)
Deoxyribonucleotide triphosphate (dNTPs) mixture 10 mM (Biobasic, catalog number: DD0056)
Agarose (Bioshop Canada Inc., catalog number: AGA002.250)
Iscove's modified Dulbecco's medium (IMDM) powder (Life Technologies, catalog number: 12200069)
Sodium bicarbonate (Fisher Scientific, catalog number: S233-3)
Hypoxanthine (Millipore, catalog number: 4010CBC-25GM)
Sodium pyruvate (Wisent Inc., catalog number: 600-110-EL)
Bathocuproin sulfonate disodium salt hydrate (Thermo Scientific, catalog number: B22550.MD)
L-cysteine hydrochloride monohydrate (Thermo Scientific, catalog number: J14035-22)
2-mercaptoethanol (Sigma, catalog number: M722)
Penicillin streptomycin (MP Biomedicals, catalog number: 1670049)
Heat-inactivated fetal bovine serum (FBS) (Life Technologies, catalog number:12484-028)
Neomycin (G418) (Sigma, catalog number: N6386)
Tetracycline hydrochloride (Fisher Scientific, catalog number: BP912-100)
Sodium phosphate dibasic anhydrous (Na₂HPO₄) (Fisher Scientific, catalog number: S374-500)
Sodium phosphate monobasic anhydrous (NaH₂PO₄) (Sigma, catalog number: S-0751)
Sodium chloride (NaCl) (Fisher Scientific, catalog number: S-271-3)
Glucose (Sigma, catalog number: G5400-250G)
Sodium citrate (dihydrate) (Fisher Scientific, catalog number: S279-500)
Bovine serum albumin (BSA) (Biobasic, catalog number: A500023-0100)
Glycerol (Fisher Bioreagents, catalog number: BP229-1)
Tris base (Fisher Scientific, catalog number: BP152-5)
Glacial acetic acid (Fisher Chemical, catalog number: A38-212)
Ethylenediaminetetraacetic acid (EDTA) (Fisher Scientific, catalog number: BP120-500)
Trypan Blue solution, 0.4% (Amresco, catalog number: K940-100ML)
Ethyl alcohol anhydrous solution (Commercial Alcohols, catalog number: P016EAAN)
96-well plate bacterial total RNA mini-prep super kit (Biobasic, catalog number: BS585-5)
Oxford Nanopore Ligation Sequencing Kit (Oxford Nanopore Technologies, catalog number: SQK-LSK109)
AMPure XP beads (Beckman Coulter Inc., catalog number: A63880)
Phleomycin (Bioshop Canada Inc., catalog number: PEO422.10)

Solutions

HMI-9 cell culture medium (see Recipes)
Drug selectable markers (see Recipes)
1. Neomycin (G-418) stock solution (20 mg/mL)
2. Tetracycline stock solution (5 mg/mL)
3. Phleomycin stock solution (1 mg/mL)
Phosphate-buffered saline-glucose (PBS-G) (see Recipes)
Bloodstream stabilate freezing solution (see Recipes)
50× Tris acetate EDTA solution (TAE) (see Recipes)
80% Ethanol (see Recipes)

Recipes

HMI-9 cell culture medium (1,000 mL)

Reagent	Final concentration	Quantity or Volume
IMDM	n/a	17.7 g
Sodium bicarbonate	0.035 M	3 g
Hypoxanthine	10 mM	10 mL
Sodium pyruvate	100 mM	10 mL
Bathocuproine sulfonate disodium salt hydrate	49.95 µM	28.2 mg
L-cysteine hydrochloride monohydrate	1.04 mM	182 mg
2-mercaptoethanol	14.3 M	14 µL
Nanopure MilliQ water	n/a	879.9 mL
Total volume	n/a	900 mL

Dissolve the IMDM and salts in 500 mL of water and stir to mix thoroughly. Add 10 mL of hypoxanthine, 10 mL of sodium pyruvate, and 14 µL of 2-mercaptoethanol. Adjust water volume to 900 mL. Filter the medium using a 0.22 µm filter system of 250 mL. Add 100 units of sterile penicillin and 100 µg/mL sterile streptomycin (optional). Complete the media by adding 10% heat-inactivated FBS. Store at 4 °C. The shelf life of the medium is four months.

Drug-selectable markers
1. Neomycin (G-418) stock solution (20 mg/mL)
  Reagent Final concentration Quantity or Volume
  Neomycin (G-418) n/a 1 g
  Nanopure MilliQ water n/a 50 mL
  Total n/a 50 mL
2. Tetracycline stock solution (5 mg/mL)
  Reagent Final concentration Quantity or Volume
  Tetracycline hydrochloride n/a 0.5 g
  Nanopure MilliQ water n/a 10 mL
  Total n/a 10 mL
3. Phleomycin stock solution (1 mg/mL)
  Reagent Final concentration Quantity or Volume
  Phleomycin n/a 0.05 g
  Nanopure MilliQ water n/a 50 mL
  Total n/a 50 mL
Dissolve the drugs in water. Filter the dissolved drug solutions using a 0.22 µm syringe filter. Prepare 1 mL aliquots into 1.5 mL Eppendorf tubes in the hood and store the aliquots at -20 °C freezer. Aliquots are stable for at least one year when stored at -20 °C.

Phosphate-buffered saline-glucose (PBS-G) (1,000 mL)

Reagent	Final concentration	Quantity or Volume
Sodium phosphate dibasic anhydrous	0.01 M	1.42 g
Sodium phosphate monobasic anhydrous	0.01 M	1.20 g
Sodium chloride	0.145 M	8.5 g
D-Glucose	0.006 M	1.081 g
Nanopure MilliQ water	n/a	1,000 mL
Total	n/a	1,000 mL

Dissolve salts in 900 mL of water. Add D-Glucose and stir thoroughly to dissolve it. Adjust the pH to 7.0 and filter sterilize using a 0.22 µm filter. Store at 4 °C. This buffer is stable for three months at 4 °C.

Bloodstream stabilate freezing solution (500 mL)
Reagent Final concentration Quantity or Volume
D-Glucose 0.052 M 9.3 g
Sodium chloride 0.036 M 2.1 g
Sodium citrate 0.003 M 0.75 g
BSA n/a 0.5 g
Glycerol n/a 75.0 g
Nanopure MilliQ water n/a 500 mL
Total n/a 500 mL
Dissolve the salts in 300 mL of water. Add D-glucose and stir thoroughly to dissolve it. Weigh the volume of glycerol equivalent to 75 g in a separate container and add it to the solution. Adjust the volume of the solution to 500 mL and stir thoroughly. Filter using a 0.22 µm filter. Store it at 4 °C. The solution is stable for six months.

50× Tris Acetate EDTA solution (TAE) (500 mL)
Reagent Final concentration Quantity or Volume
Tris base 1 M 121 g
Glacial acetic acid 1 M 28.6 mL
EDTA 0.05 M 50 mL
Nanopure MilliQ water n/a 421.4 mL
Total volume n/a 500 mL
Prepare a stock solution of 0.5 M EDTA in a separate tube by dissolving 18.6 g of EDTA disodium salt in 80 mL of water. Adjust the pH to 8.0 with 10 M NaOH solution. Adjust the volume to 100 mL. Dissolve the Tris base in 300 mL of water and add 28.6 mL of 17.4 M glacial acetic acid and 50 mL of 0.5 M EDTA. Adjust the volume of the solution to 500 mL with water. Autoclave the solution. Store it at room temperature (RT). The solution is stable for at least six months.

80% Ethanol
Reagent Final concentration Quantity or Volume
Ethyl alcohol anhydrous solution (100% v/v) 80 % 80 mL
Nanopure MilliQ water n/a 20 mL
Total volume n/a 100 mL
Prepare 80% v/v ethanol by transferring 80 mL of 100% Ethyl alcohol solution to a 100 mL graduated cylinder. Then, add 20 mL of Nanopure MilliQ water. Store at 4 °C. The solution is stable for one month.

Laboratory supplies

Tissue culture flasks
1. 25 cm² TC-treated T-flask with filter cap (Biobasic, catalog number: SP81136)
2. 75 cm² TC-treated T-flask with filter cap (Biobasic, catalog number: SP81186)
Cell culture plates

96-well cell culture plates (WUXI NEST Biotechnology Co., catalog number: 101722BL01)
24-well TC plates, treated (Biobasic, catalog number: SP41135)

Sterile pipette tips

1,000 µL pre-sterile barrier tips (Neptune, catalog number: BT100.96)
200 µL pre-sterile barrier tips (Neptune, catalog number: BT200)
10 µL pre-sterile barrier tips (Neptune, catalog number: BT10)

96-well 2 mL deep plate, natural, edge filled, 10 plates/bag (Biobasic, catalog number: BR581-96NS)
96-well deep collection plates (Biobasic, catalog number: 107-E627LA2221)
Plate seals (Biobasic, catalog number: BS585-5)
96-well PCR plate (Life Technologies, catalog number: 4346907)
Foil sealing film, non-sterile (Celltreat Scientific Products, catalog number: 501535152)
Multichannel pipette (Eppendorf, catalog number: 4056991)
Heating block (Eppendorf, catalog number: 535028642)
Filter system 250 mL 0.2 µm PES (Fisher Scientific, catalog number: FB12566502)
2.0 mL microtubes (UltiDent Scientific, catalog number: 48-C200-CS)
1.5 mL Eppendorf tubes (Eppendorf, catalog number: 0030123611)
100 mL graduated cylinder (Grainger, catalog number: GGS5PTJ5)
100 mL reagent bottle (UltiDent Scientific, catalog number: 170-14170100)

Equipment

Ultrafocused sonicator (Covaris, M220, catalog number: 006168)
Magnetic rack for 1.5 mL tubes (Promega, catalog number: PR-Z5342)
T100 thermal cycler (Bio-Rad, catalog number: 1861096)
Avanti J-E centrifuge (Beckman Coulter, catalog number: JSE02M13)
Allegra^TM 25R centrifuge, TJ-25 rotor (Beckman Coulter, catalog number: AJC024001)
Nanodrop (Nanodrop Spectrophotometer, catalog number: ND-1000)
Microcentrifuge (Eppendorf, catalog number: EP5401000137)
Water system ultrapure (Nanopure MilliQ water) (Millipore Synergy, catalog number: F1CA45528 A)
Electrophoresis unit: Thermo Scientific Power Supply 400 mA 300 V (Fisher, catalog number: S65533Q)
Biological safety cabinet (NuAire Biological Safety Cabinet Class II Type AIB3, catalog number: NU-425-300)
CellDrop^TM automated cell counter (DeNovix Cell Drop FL Fluorescence Cell Counter)
HEPA CO₂ incubator (Thermo Electron Corporation; Forma Series II, Water Jacketed, catalog number: 308606-30529)
Microscope (Nikon, model: Eclipse TS100, catalog number: 302115)
-80 °C freezer (New Brunswick Ultra-Low Temperature Freezer U101 Inova, catalog number: U9420-0000)
MinION Mk1C (Oxford Nanopore Technologies, catalog number: MIN-101C)
Flongle Flow Cell (R9.4.1) (Oxford Nanopore Technologies, catalog number: FLO-FLG001)
Rotator (ManSci Inc., catalog number: A706514)
Milli-Q IQ 7000 purification system (Millipore Sigma, catalog number: ZIQ7000T0C)

Software and datasets

Minimap2, version 2.24 (Li, 2018) (https://github.com/lh3/minimap2)
Samtools, version 1.17 (Danecek et al., 2021) (https://github.com/samtools/samtools)
DeepTools, version 2.0 (Ramírez et al., 2016) (https://deeptools.readthedocs.io/en/develop/index.html)
Subread, version 2.0.3 (Liao et al., 2014) (http://subread.sourceforge.net/featureCounts.html)
Rcgrep, version 0.1 (https://github.com/dib-lab/rcgrep)
Integrative Genomics Viewer (IGV), version 2.16.2 (Robinson et al., 2011) (https://software.broadinstitute.org/software/igv/)
VSG-BarSeq (this work), version 1.0.1 (https://github.com/cestari-lab/VSG-Bar-seq)

Procedure

Parasite treatment and cloning
We recommend determining cell treatment conditions before starting this protocol. The conditions used in this protocol were optimized for T. brucei bloodstream forms of the 427 strain or conditional null (CN) cells derived from the single-marker 427 strain (Cestari and Stuart, 2015). The treatment described here is the knockdown for 24 h of the T. brucei gene encoding PIP5Pase, which results in high rates of VSG switching (Touray et al., 2023). The PIP5Pase CN cell line is grown in G418 and phleomycin to maintain the selection of the tetracycline-inducible system. Tetracycline is added to induce expression of the PIP5Pase gene under the control of a procyclin promoter and a tetracycline operator.
1. Grow 5 mL of T. brucei cells seeded at 1.0 × 10⁴ cells/mL in HMI-9 medium supplemented with 2 µg/mL G418, 2.5 µg/mL phleomycin, and 500 ng/mL tetracycline in a 37 °C incubator with 5% CO₂ for 24 h or until it reaches mid-log growth (~1.0 × 10⁶cells/mL). Throughout this protocol, cell growth will be as described above unless otherwise stated. Cells’ doubling time should be approximately 5.5–6 h and viability approximately 90%–95%. Avoid overgrowing cell culture (>1.5 × 10⁶cells/mL) because it will affect cell viability.
2. Transfer the 5 mL cell culture to a 15 mL Falcon tube and centrifuge at 3,500× g for 5 min at RT. Discard the supernatant.
3. Resuspend the pellet in 10 mL of PBS-G pre-warmed at 37 °C and then centrifuge the cells as in step A2. Discard the supernatant.
4. Repeat step A3 three times to ensure complete removal of tetracycline.
5. Split the cells into two 5 mL cell culture flasks (treated and non-treated groups), seeding each culture at 1.0 × 10⁴ cells/mL in HMI-9 medium with 2 µg/mL G418 and 2.5 µg/mL phleomycin.
6. Add 500 ng/mL tetracycline to the non-treated flask (Tet +, control) and no tetracycline to the treatment flask (Tet -, knockdown) and grow the cells for 24 h.
7. Add 500 ng/mL tetracycline to the treatment flask (Tet -, knockdown). No additional tetracycline is required for the non-treated flask (Tet +, control). Quantify the cell concentration and viability of both the treatment and control groups by mixing 10 µL of cell culture and 10 µL of 0.4% Trypan blue staining. Add 10 µL to the CellDrop^TM cell counter to obtain viability and cell concentration. Cell concentration should be approximately 1.0 × 10⁵ cells/mL, and viability should be >90%.
8. Add 9 mL of HMI-9 medium supplemented with 500 ng/mL tetracycline to a 50 mL Falcon tube. Repeat the procedure to have three flasks for the treatment group and three for the non-treatment group.
9. From a starting cell concentration of 1.0 × 10⁵ cells/mL (from A7 above), gently mix the cells by flicking the flasks five times and transfer 1 mL of the culture to the Falcon tube 1 containing 9 mL of HMI-9 media (1 in 10 dilutions) to obtain a cell concentration of 1 × 10⁴ cells/mL.
10. From the Falcon tube 1 (1 × 10⁴ cells/mL) culture, transfer 1 mL to Falcon tube 2 to obtain a cell concentration of 1 × 10³ cells/mL (1 in 10 dilutions). Repeat the same for Falcon tube 3 to obtain a cell concentration of 1 × 10² cells/mL. Perform the procedure for control and treated groups.
11. Transfer 7 mL of the culture from Falcon tube 3 to a 75 cm² cell culture flask containing 63 mL of HMI-9 medium supplemented with 2 µg/mL G418, 2.5 µg/mL phleomycin, and 500 ng/mL tetracycline, to obtain a cell concentration of 10 cells/mL (1 in 10 dilutions).
12. Transfer 60 mL of the diluted cell culture from step A11 to a 75 cm² cell culture flask and add 140 mL of HMI-9 medium supplemented with 2 µg/mL G418, 2.5 µg/mL phleomycin, and 500 ng/mL tetracycline (3 in 10 dilutions) to obtain a final cell concentration of 3 cells/mL.
13. Aliquot the diluted parasite culture (3 cells/mL) onto ten 96-well cell culture plates (10 plates per treatment group) using a multichannel pipette. Transfer 200 µL to each well so that the probability of obtaining one single cell per well (200 µL) is approximately 30%, i.e., one cell per well for a third of the wells of a 96-well cell culture plate. Ensure the cells are well mixed by gently swirling the parasite cultures while pipetting.
14. Grow the cells in 96-well cell culture plates in an incubator for 5–7 days.
15. Check the 96-well cell culture plates under a microscope to identify parasite clones. Approximately 30% of the wells should contain parasites (see Note 1). We recommend checking each well for parasite clones starting from day 5 up until day 7 post seeding.
16. Transfer 200 µL of parasite clones onto 24-well TC plates and add 1.8 mL of fresh HMI-9 with 2 µg/mL G418, 2.5 µg/mL phleomycin, and 500 ng/mL tetracycline.
17. Grow the clones in the 24-well TC plates for 24 h to increase the number of cells for RNA extraction.
18. (Optional) Aliquot 600 µL of each clonal population onto new 96-well deep collection plates and add 600 µL of freezing solution. Freeze the parasites at -80 °C for short-term storage (2–3 weeks) or liquid nitrogen for long-term storage.
19. The remaining 1.4 mL of the clonal parasite cultures are used for RNA extraction.

RNA extraction
1. Transfer 1.4 mL of each clonal culture (approximately 1.4 × 10⁶cells) into 96-well deep collection plates, centrifuge at 3,500× g for 5 min at RT, and pour off the supernatant (see Note 2).
2. Add 350 µL of buffer Rlysis-BG (provided in the 96-well plate bacterial total RNA mini-prep super kit) to each well and resuspend the cells by pipetting up and down five times.
3. Thoroughly seal the plate with a sealing film to prevent cross-contamination of the samples and immediately mix the samples by inverting the plate three times.
4. Briefly spin the plate at 1,000× g for 30 s to collect the solution to the bottom of the wells.
5. Add 175 µL of absolute ethanol to each well, tightly seal the plate with a new sealing film, and mix thoroughly by inverting five times.
6. Place the EZ-10 96-well plate (filtration column plate provided with the 96-well plate bacterial total RNA mini-prep super kit) on top of a new 96-well deep collection plate and transfer the lysate from step B5 into the columns on the EZ-10 96-well plate.
7. Centrifuge the plate at 5,000× g for 2 min at RT and discard the flowthrough.
8. Place the EZ-10 96-well plate back on the deep cell collection plate and add 500 µL of universal GT solution (provided in the kit) to each column.
9. Centrifuge the plate at 5,000× g for 1 min at RT and discard the flowthrough.
10. Place the EZ-10 96-well plate back on the deep well collection plate and add 500 µL of universal NT solution (provided in the kit) to each column.
11. Centrifuge the plate at 5,000× g for 1 min at RT and discard the flowthrough.
12. Place the EZ-10 96-well plate back on the deep well collection plate and centrifuge the column at 5,000× g for 2 min at RT to ensure complete removal of the residual ethanol.
13. Place the EZ-10 96-well plate into a new deep-well storage plate (provided in the kit), add 30 µL of RNase-free water (supplied in the kit), and then incubate the plate at RT for 5 min.
14. Centrifuge the plate at 5,000× g for 1 min at RT to elute the RNA solution.
15. Quantify the recovered RNA by measuring 1 µL of the RNA solution at 260 nm from approximately 10 random wells using a NanoDrop. This will help estimate the isolated RNA concentrations.
16. Tightly seal the plate with an adhesive cover and keep it on ice (or at -80 °C) until the RNA samples are ready for cDNA synthesis.

VSG-enriched cDNA synthesis and barcoding
This step requires a combination of primers (Figure 1) to barcode each clonal cell population cDNAs with a unique eight-nucleotide sequence for their identification during sequencing analysis. It will also provide an adapter sequence for DNA sequencing library preparation. The forward Ad-SL (5′-TTTCTGTTGGTGCTGATATTGCacagtttctgtactatattg-3′) primer is universal and includes an Oxford nanopore adaptor sequence (capital letters) followed by a sequence (small letters) that hybridizes to mRNA splice leader sequence, a 39-nt sequence added to the 5′ of all trypanosomes’ mRNAs. The reverse Ad-3endVSG primer (5′-TACTTGCCTGTCGCTCTATCTTCXXXXXXXXgtgttaaaatatatc-3′) contains an Oxford nanopore adaptor sequence (capital letters) followed by eight-nucleotide long variable sequence (barcode) unique to each clone and a sequence pairing with the conserved 3′-end of VSG mRNAs, which encodes the C-terminus of VSG proteins (Mugnier et al., 2015). See Supplementary information for the complete primer list. We recommend preparing a working primer solution containing a mix of both primers at 10 µM in a 96-well plate.

Figure 1. Schematic representation of clonal cell barcode and nanopore sequencing protocol. (A) Diagram of the clonal cell barcoding, nanopore sequencing workflow showing the parasite treatment, and cloning (Step 1), RNA isolation from the individual clones, cDNA synthesis, and clonal barcoding (Step 2), and Oxford nanopore technology (ONT) library preparation, sequencing, and data analysis using VSG-BarSeq script (Step 3). (B) Scheme of the forward Adaptor Splice Leader (Ad-SL) primer and the reverse Adaptor-barcode-3′-end VSG primer (Ad-3endVSG) used for cDNA synthesis and clonal cell barcoding. (C) Diagram describing the cDNA synthesis, clonal barcoding, ONT library barcoding, and PCR amplification. (D) Sequences of the Ad-SL and Ad-3endVSG annealing primers used for cDNA synthesis and clonal cell barcoding (see the annealing regions in C). The sequences in black correspond to ~20 bp nanopore barcode adapter sequences, and the sequences in pink and green correspond to the forward Ad-SL and reverse Ad3endVSG primers, respectively. The Ad-SL pairs to the mRNA splice leader sequence, while the reverse Ad3endVSG primer pairs to the conserved 3′-end of VSG mRNAs. The 8-nucleotide long variable sequence (clonal barcode) unique to each clone is shown in orange and depicted here as X. A complete list of primers is available in Supplementary Information. PBS-G: phosphate-buffered saline (PBS)-glucose; ONT: Oxford Nanopore Technology; PIP5Pase CN: phosphatidylinositol 5-phosphatase (PIP5Pase) conditional null.
1. Take 4 μL of primer mix from the working primer solution plate and transfer it to each well of a new 96-well PCR plate using a multichannel pipette.
2. Add 1 μL of 10 mM dNTPs mix onto each well of the same 96-well PCR plate.
3. Thaw the RNA samples (from step B16) on ice if frozen at -80 °C and add 5 μL of the samples, keeping the same orientation in the 96-well PCR plate as the original RNA plate. Mix the solutions by gently pipetting up and down five times and spin down the plate briefly at 1,000× g for 30 s at 4 °C.
4. Incubate the plate at 65 °C for 5 min in a thermocycler and then transfer the plate immediately to ice.
5. Prepare cDNA synthesis master mix by adding 2 μL of MuLV 10× Buffer, 5 U of M-MuLV Reverse Transcriptase, and 8 U of RNase inhibitor (supplied in the MuLV reverse transcriptase kit) and adjust the reaction volume to 10 μL per reaction using nuclease-free water.
6. Aliquot 10 μL of the cDNA synthesis master mix onto each well of the 96-well plate containing the RNA, dNTPs, and primer mix.
7. Thoroughly seal the plate using an adhesive plate sealing film, spin down at 1,000× g for 30 s, and incubate at 42 °C for 2 h and then 65 °C for 20 min in a thermocycler.
8. Store the synthesized cDNA samples at -80 °C or proceed to library preparation.

Oxford nanopore library preparation and DNA sequencing
1. Combine the synthesized cDNA samples from each well (from step C8) into one 1.5 mL microcentrifuge tube and mix gently five times. Avoid forming bubbles.
2. Prepare 10 barcoding PCR reactions containing 5 μL of 10× ThermoPol buffer, 400 μM dNTPs mix, 500 nM barcode primer mix (Native Barcoding Expansion 1-12), 2 U of Taq DNA Polymerase, and 750 nM MgCl₂. Add 3 μL of the pooled cDNA and adjust the final reaction volume to 50 μL with nuclease-free water.
3. Perform the PCR reaction in a thermocycler at 95 °C for 10 min, then 22 cycles at 95 °C for 1 min, 62 °C for 1 min, and 68 °C for 3:30 min, and a final extension at 68 °C for 10 min (see Notes 3 and 4).
4. Pool together all PCR amplicons in one 1.5 mL Eppendorf tube to obtain a final volume of 500 μL. Add 350 μL (0.65× beads to sample ratio) of NucleoMag NGS Clean-up and Size Selection beads to clean up the DNA.
5. Incubate beads/samples in a rotator (60 rpm) for 10 min at RT.
6. Spin down the tube briefly at 1,000× g for 30 s at RT, place it on a magnetic rack, and incubate for 1 min at RT. Pipette off and discard the supernatant without disturbing the beads.
7. Add 350 µL of freshly prepared 80% ethanol to the beads gently. Avoid disturbing the beads or taking the tube out of the magnetic rack. The volume of the 80% ethanol to add per wash should be equal to the sum of the total PCR amplicon volume and added NucleoMag NGS beads volume (e.g., 500 µL of PCR reaction + 350 µL of NucleoMag NGS magnetic beads = 850 µL of the mix, then 850 µL of 80% ethanol).
8. Incubate the mix for 1 min on the magnetic rack at RT. Collect and discard the 80% ethanol. Repeat steps D7 and D8.
9. Centrifuge the tube at 1,000× g for 30 s at RT. Place the tube back on the magnetic rack and remove any residual 80% ethanol from the tube.
10. Remove the tube from the magnetic rack and let the beads air dry for 10 min at RT.
11. Resuspend beads in 61 μL of nuclease-free water, mix gently by flicking the tube or gently pipetting up and down five times, and incubate for 5 min at RT.
12. Place the tube back on the magnetic rack and incubate for 1 min at RT.
13. Collect and transfer the clear supernatant containing the eluted DNA into a new labeled tube. Discard the tube containing the beads.
14. Measure and assess the concentration and purity of the eluted DNA by quantifying 1 µL of the DNA sample using NanoDrop. The DNA yield ranges from approximately 30 to 60 ng/μL with 260/280 and 260/230 absorbance ratios of approximately 1.8 and 2.0, respectively.
15. Prepare 1 μg of DNA in 47 μL of nuclease-free water in a 0.2 mL thin-walled PCR tube and then add the following: 3.5 μL of NEBNext FFPE DNA repair buffer and 2 μL of NEBNext FFPE DNA repair mix (from NEBNext FFPE DNA Repair Mix), 3.5 μL Ultra II End-prep reaction buffer and 3 µL Ultra II End-prep enzyme mix (from NEBNext end repair module), and 1 μL of DNA CS (provided in the Oxford Nanopore Ligation Sequencing kit).
16. Mix the solutions by pipetting up and down gently five times and incubate for 60 min at 20 °C in a thermocycler.
17. Add 39 μL of NucleoMag NGS Clean-up and Size Selection beads to the end-repaired DNAs (0.65× beads to samples ratio) to clean up and follow steps D5–D13. Elute the DNA in 61 μL of nuclease-free water.
18. Quantify 1 µL of the eluted DNA using NanoDrop. Approximately 90% of the DNA is expected to be recovered.
19. Transfer the remaining 60 µL of end-repaired DNA volume from step D18 (0.5–1 μg) to a 0.2 mL thin-walled PCR tube and then add the following components: 25 μL of ligation buffer (LNB), 10 μL of NEBNext Quick T4 DNA ligase (from NEBNext Quick Ligation Module), and 5 μL of Adapter mix F (AMX-F) (LNB and AMX-F are components of the Oxford Nanopore Ligation sequencing kit). Mix by gently pipetting the reaction up and down five times, then briefly spin it down.
20. Incubate the reaction for 2 h at 20 °C in a thermocycler.
21. Transfer the ligation reaction mix from step D20 to a new 1.5 mL Eppendorf tube and add 40 μL of well resuspended AMPure XP beads to the reaction (0.4× beads to sample ratio). Mix by gently pipetting up and down five times.
22. Incubate the mix on a rotator (60 rpm) for 10 min at RT.
23. Briefly centrifuge the sample at 1,000× g for 30 s and place the tube on a magnetic stand for 1 min at RT or until the supernatant is clear and colorless.
24. Remove and discard the supernatant without disturbing the beads.
25. Resuspend the beads in 250 μL of Short Fragment Buffer (SFB) (provided in the Oxford Nanopore Ligation Sequencing kit) and mix gently by flicking the tube five times to wash the beads.
26. Briefly centrifuge the tube at 1,000× g for 30 s. Transfer the tube to a magnetic stand to pellet the beads. Remove and discard the supernatant without disturbing the beads.
27. Repeat steps D25 and D26.
28. Briefly centrifuge the tube at 1,000× g for 30 s. Place it on the magnetic stand and remove any residual supernatant. Air dry the beads at RT and ensure that they do not dry to the point of cracking. Ten minutes is sufficient for the beads to dry. Extended drying of the beads for 20–30 min might cause them to crack.
29. Remove the tube from the magnetic stand and resuspend the beads in 12 μL of Oxford nanopore Elution Buffer (EB) (provided in the Oxford Nanopore Ligation Sequencing Kit) by gently pipetting up and down.
30. Briefly centrifuge the tube at 1,000× g for 30 s and incubate for 10 min at RT.
31. Place the tube back onto the magnetic stand for 1 min and then remove and transfer the eluate containing the DNA library into a new 1.5 mL Eppendorf tube. Dispose of the beads.
32. Quantify 1 µL of the eluted library DNA sample using a NanoDrop. Expect a DNA library yield of approximately 10–40 ng/μL (i.e., 200–1,500 fmol) of DNA.
33. Load 3–20 fmol of the prepared library into the Flongle Flow Cell (R9.4.1) and sequence the DNA according to the manufacturer's instructions.

Data analysis

Computational analysis of the sequenced library
The data analysis described here was performed using a Linux operating system (Ubuntu). The computational resources required will vary depending on the data available for analysis. The Oxford nanopore sequencer will generate fast5 files, which are basecalled to fastq files using the Guppy tool integrated in the MinKNOW software (https://nanoporetech.com/). The fastq files are the input dataset in the analysis shown here. We developed a computational pipeline for the sequencing analysis and detection of VSG switching (Figure 2). The pipeline is run via the vsg-barseq.sh script. A fastq file with 100,906 reads was analyzed using 10 threads and 4 GB of memory and completed in 5 min. The pipeline takes the DNA sequences generated by the nanopore sequencing in fastq format and splits them into subfiles according to the cell barcodes (eight mers) used for DNA-seq library preparation. The output fastq files are named by their barcodes, e.g., CCATGCAT.fastq. A split summary file (.txt) is generated and shows the number of reads per barcode file, the total count of reads analyzed, and the total amount of reads containing barcodes (see Note 5). Sequences from each file are then mapped to the organism genome (here, T. brucei 427 strain) using minimap2, outputting .sam files. The alignments are then filtered using Samtools to remove supplementary alignments and to keep alignments with a mapping quality score (mapQ) ≥ 20. The resultant .sam files are sorted and indexed with Samtools resulting in sorted.bam and bam.bai files, respectively.
The alignments are counted with featureCounts (package Subread). The top mapped reads, which correspond to the expressed VSG, are selected and compiled in a single tab-delimited file (topmapped.txt) while keeping the original output from featureCounts, which serves for analysis of other genes identified during alignment. The analysis of other identified genes is part of the quality control process. It typically shows genes with low counts reflecting low background noise from library preparation resulting primarily from sequences derived from splice-leader cDNA synthesis (Figure 3E).
The vsg-barseq.sh is open code and available on GitHub (https://github.com/cestari-lab/VSG-Bar-seq). The script takes six arguments:
1) Directory of all files (output folders will be created in this directory)
2) Directory of fastq files
3) Directory of barcode.txt files (barcode.txt file has one barcode sequence per row)
4) Directory of genome file in fasta format
5) Directory of gene transfer format file, i.e., gtf format
6) Number of threads to be passed to minimap2, samtools, and featureCounts (we recommend eight or more)

Example of barcode.txt file. Keep one barcode sequence per row.
CCGTTAGG
CCAACTAG
CAGGGCAG
CGCAGAAG

If the barcoded.txt file is generated using Windows operational system, we recommend converting from dos to Unix format before executing the analysis. The tools minimap2, samtools, subread, and rcgrep are required to run the script. After installing and/or loading the required tools, run the script as indicated below:

sh vsg-barseq.sh path/to/directory path/to/fastq path/to/barcodes path/to/genome path/to/gtf nthreads

See the example below using the structure folder exemplified (directories in bold).

sh vsg-barseq.sh \
c/user/vsg-barseq \
myfastq \
barcodes \
genomefiles/genome.fasta \
genomefiles/genome.gtf \
8

Note that T. brucei genome (.fasta) and features (.gff, general feature format) can be downloaded from TritrypDB (https://tritrypdb.org/). TriTrypDB does not provide .gtf files, but .gff files can be used to generate .gtf. We recommend using the gffread tool (https://github.com/gpertea/gffread), also available via Galaxy tools (https://usegalaxy.org/). The command below indicates how to generate a gtf file from a gff file after installing gffread tools.

gffread annotation_file.gff -T -o annotation_file.gtf

T. brucei barcodes and test fastq data files are also available for download at https://github.com/cestari-lab/VSG-Bar-seq to test the script.

Directory structure of input files required to run the script:

c/user/vsg-barseq/
----myfastq/mydata.fastq
----barcodes/barcode.txt
----genomefiles/genome.fasta, genome.gtf

After running the script, the following output directories and files are created:

c/user/vsg-barseq/
----result_split_fastq/CCGTTAGG.fastq, CCAACTAG.fastq, CAGGGCAG.fastq, CGCAGAAG.fastq, split_summary.txt
----result_mapcount/
----sam/CCGTTAGG.sam, CCAACTAG.sam, CAGGGCAG.sam, CGCAGAAG.sam,
----bam/CCGTTAGG.bam, CCAACTAG.bam, CAGGGCAG.bam, CGCAGAAG.bam
----sorted_bam/CCGTTAGG_sorted.bam, CCGTTAGG_sorted.bam.bai, etc.
----counts/CCGTTAGG.txt, CCGTTAGG.txt.summary, etc.
----topmap/topmapped.txt

The result file topmapped.txt will have information on gene id (first row, in bold below), chromosome (or contig), nucleotide position (start, end), strand (+/-), gene length in nucleotides, number of counted alignments for the gene (seventh row, in bold below), and path to the original file containing all alignment counts for the corresponding clone; the eight mers barcode identifies the cell expressing the VSG. Each row represents a cell clone, and the gene id indicates the VSG expressed by the cells in the clonal population.

Tb427_000016000 BES1_Tb427v10 75570 77000 + 1431 461 /path/to/counts/CCGTTAGG.txt
Tb427_000016000 BES1_Tb427v10 75570 77000 + 1431 363 /path/to/counts/CCAACTAG.txt
Tb427_000008000 BES12_Tb427v10 45278 46651 + 1374 310 /path/to/counts/CAGGGCAG.txt
Tb427_000284800 Chr2_5A_Tb427v10 197226 198701 + 1476 840 /path/to/counts/CGCAGAAG.txt

Figure 2. Flowchart of computational analysis using VSG-BarSeq. Multiplexed reads from clonal VSG-seq are split into clone-specific reads based on eight mers barcode. Reads are aligned to the genome using minimap2 and filtered with Samtools to remove supplementary and secondary alignments and keep alignments with mapQ ≥ 10. Then, it counts the alignments per gene and reports the top alignment per file, corresponding to the expressed variant surface glycoproteins (VSG) gene.

Anticipated results
The analysis of parasite VSG switching starts with treating the cultures to induce the switching of the VSG gene. The treatment conditions and outcomes might differ for different cell lines, so they may need to be optimized. Here, the treatment was the temporal knockdown (24 h) of the PIP5Pase gene, which results in high rates of VSG switching (Touray et al., 2023). VSG switching occurs by alternating transcription between ESs or recombining VSG genes within ESs (Figure 3A). The parasite culture was diluted to obtain approximately 30% of clones from a 96-well plate. We recommend optimizing the dilution of the cells after the treatment and checking cell viability before the VSG-BarSeq experiment. After cDNA synthesis and library amplification with ONT-barcoding primers, we recommend analysis of the product in agarose gel. We often obtain a PCR fragment smear ranging from 200 to 2,000 bp (Figure 3B).

Figure 3. Analysis of variant surface glycoproteins (VSG) switching after PIP5Pase knockdown using VSG-BarSeq. (A) Diagram of bloodstream-form expression sites (ESs) (BES) and sub-telomeric regions containing VSG genes used for recombination. The pink arrowhead represents a BES promoter. VSG genes are transcribed from BESs only. (B) Thirty-five cycles of quality-control PCR amplification of pooled clonal VSG barcoded cDNAs. PCR was amplified with ONT barcode primers. NC: negative control, M: 5 kb DNA ladder. (C) VSG genes expressed by clones of T. brucei without PIP5Pase knockdown (Tet +) and after temporary knockdown for 24 h followed by PIP5Pase re-expression and cloning for 5–7 (Tet -/+). Annotated names and corresponding chromosomes or BES identify VSG genes. The number of clones analyzed is indicated in parentheses. (D) Read coverage plots of example clones from Tet + or Tet -/+ treatment groups. The diagram on the top right summarizes the experiment and results—all clones derived from an original clone expressing VSG2 (BES 1). Graphs show read coverage of expressed VSGs on the same scale. (E) Signal-to-noise ratios of two different VSG-BarSeq libraries differing in sequencing depth. Sequencing depth is shown above plot bars as the total reads by the mean read length. Mb: megabases.

After library preparation, we usually obtain a library concentration ranging between 10 and 40 ng/μL with 260/280 and 260/230 ratios of ~1.8 and 2.0, respectively. Library concentrations below 5 ng/µL with considerable deviations from the above 260/280 and 260/230 ratios usually result in inferior quality and low throughput sequencing. Analysis of the sequencing using the vsg-barseq.sh script typically results in 70%–99% genome mapping. Figure 3C and 3D show the results of VSG switching in T. brucei after the temporary knockdown of PIP5Pase. Analysis of 117 clones from the control (no knockdown) showed that none of the cells switched VSGs (Figure 3C). However, after 24 h temporary knockdown, there were 93 switchers out of 94 clones analyzed. There was a preference for the cells to switch to VSG8 in the BES12, suggesting transcriptional switching. Moreover, switching to VSGs from sub-telomeric regions (Chr2_5A and Chr9_3A) was also detected, indicating switching by recombination (Figure 3C and D). The analysis does not require a significant amount of RNAs or sequencing throughput since the experimental setup relies on selecting VSG mRNAs, which are highly abundant (Cestari and Stuart, 2015), and on the analysis of clones rather than heterogeneous cell populations. Analysis of the expressed VSG (signal) compared to other genes (noise) showed a high signal-to-noise ratio, and increasing sequencing depth improved the signal without significantly increasing noise levels (Figure 3E). Although DNA sequencing can be costly, the amount of total RNA for cDNA synthesis per clone in this protocol is minimal (5–20 ng), and sequencing depth required for analysis can be obtained from Oxford nanopore flongle flow cells, which typically produces 500–1,000 megabases of DNA sequencing (i.e., ~100,000 reads per group), thus reducing experimental costs. The results show the method's utility in identifying VSG switching events using multiplexed clonal cell barcodes. We anticipate that, with minimal modifications of the primers used for the target sequences, the approach can be applied to other gene families or other cell types, including var genes in Plasmodium sp., variant surface proteins in Giardia sp., as well as mucins and mucin-associated proteins in Trypanosoma cruzi. Analogously, the approach can be extended to other organisms or cell lines, including mammalian cells, e.g., T-cell or B-cell repertoire analysis.

Validation of protocol

This protocol has been used and validated in the following article:

Touray et al. (2023). A PI(3,4,5)P3-dependent allosteric switch controls antigenic variation in trypanosomes. eLife 12: RP89331 (Figure 1F, and Figure 1—figure supplement 1, panel A, B, C, D).

General notes and troubleshooting

General notes

If more than 40% of the wells in a 96-well plate are positive, the parasites might not be clonal. We recommend optimizing the dilution of the parasites to obtain approximately a third of the wells from a 96-well plate containing growing parasites.
We recommend performing a quality control PCR prior to cDNA amplification. Use the same conditions indicated in step D3 but for 35 cycles. Then, analyze 30 µL the amplicons on 1% agarose/TAE electrophoresis gel with 5 µL of Ecostain at 85 Volts for 45 min and visualize using a gel imager. A smear migrating at 200–2,000 bp should be expected in the agarose gel (Figure 3B).
We recommend the volume of pooled cDNA sample added to the PCR reaction to be less than one-tenth of the total PCR reaction volume.
Reverse transcriptase enzymes are known to inhibit PCR, particularly at low template concentrations (Chandler et al., 1998). Therefore, adding more cDNA to the PCR reactions usually results in non or very little amplification.
The split summary file helps to identify the splitting of reads into multiple barcode fastq files. If the number of reads with barcodes is larger than the total number of reads, it is indicative that some reads are included in more than one barcode fastq file. Primers with longer barcodes than eight mers could be used.

Acknowledgments

This research was funded by the Canadian Institutes of Health Research (grant CIHR PJT-175222 to IC), the Natural Sciences and Engineering Research Council of Canada (grant RGPIN-2019-05271 to IC), the Fonds de Recherche du Québec - Nature et technologie (grant 2021-NC-288072 to IC), and the Canada Foundation for Innovation (grant JELF 258389 to IC). AOT is a recipient of the Islamic Development Bank Scholarship (grant 600042744). This research was enabled in part by computational resources provided by Calcul Quebec (https://www.calculquebec.ca/en/) and the Digital Research Alliance of Canada (alliancecan.ca). This protocol is derived from the original research paper by Touray et al. (2023).

Competing interests

The authors declare that they have no conflicts of interest with the contents of this article.

References

Aitcheson, N., Talbot, S., Shapiro, J., Hughes, K., Adkin, C., Butt, T., Sheader, K. and Rudenko, G. (2005). VSG switching in Trypanosoma brucei: antigenic variation analysed using RNAi in the absence of immune selection. Mol. Microbiol. 57(6): 1608–1622.
Boothroyd, C. E., Dreesen, O., Leonova, T., Ly, K. I., Figueiredo, L. M., Cross, G. A. M. and Papavasiliou, F. N. (2009). A yeast-endonuclease-generated DNA break induces antigenic switching in Trypanosoma brucei. Nature 459(7244): 278–281.
Cestari, I., McLeland-Wieser, H. and Stuart, K. (2019). Nuclear Phosphatidylinositol 5-Phosphatase Is Essential for Allelic Exclusion of Variant Surface Glycoprotein Genes in Trypanosomes. Mol. Cell. Biol. 39(3): e00395–18.
Cestari, I. and Stuart, K. (2015). Inositol phosphate pathway controls transcription of telomeric expression sites in trypanosomes. Proc. Natl. Acad. Sci. U.S.A. 112(21): E2803–E2812.
Cestari, I. and Stuart, K. (2018). Transcriptional Regulation of Telomeric Expression Sites and Antigenic Variation in Trypanosomes. Curr. Genomics 19(2): 119–132.
Chandler, D. P., Wagnon, C. A. and Bolton, H. (1998). Reverse Transcriptase (RT) Inhibition of PCR at Low Concentrations of Template and Its Implications for Quantitative RT-PCR. Appl. Environ. Microbiol. 64(2): 669–677.
Danecek, P., Bonfield, J. K., Liddle, J., Marshall, J., Ohan, V., Pollard, M. O., Whitwham, A., Keane, T., McCarthy, S. A., Davies, R. M., et al. (2021). Twelve years of SAMtools and BCFtools. GigaScience 10(2): giab008.
Faria, J., Glover, L., Hutchinson, S., Boehm, C., Field, M. C. and Horn, D. (2019). Monoallelic expression and epigenetic inheritance sustained by a Trypanosoma brucei variant surface glycoprotein exclusion complex. Nat. Commun. 10(1): e3023.
Glover, L., Alsford, S. and Horn, D. (2013). DNA Break Site at Fragile Subtelomeres Determines Probability and Mechanism of Antigenic Variation in African Trypanosomes. PLoS Pathog. 9(3): e1003260.
Jehi, S. E., Li, X., Sandhu, R., Ye, F., Benmerzouga, I., Zhang, M., Zhao, Y. and Li, B. (2014). Suppression of subtelomeric VSG switching by Trypanosoma brucei TRF requires its TTAGGG repeat-binding activity. Nucleic Acids Res. 42(20): 12899–12911.
Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34(18): 3094–3100.
Liao, Y., Smyth, G. K. and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30(7): 923–930.
López-Escobar, L., Hänisch, B., Halliday, C., Ishii, M., Akiyoshi, B., Dean, S., Sunter, J. D., Wheeler, R. J. and Gull, K. (2022). Stage-specific transcription activator ESB1 regulates monoallelic antigen expression in Trypanosoma brucei. Nat. Microbiol. 7(8): 1280–1290.
Mugnier, M. R., Cross, G. A. M. and Papavasiliou, F. N. (2015). The in vivo dynamics of antigenic variation in Trypanosoma brucei. Science 347(6229): 1470–1473.
Ramírez, F., Ryan, D. P., Grüning, B., Bhardwaj, V., Kilpert, F., Richter, A. S., Heyne, S., Dündar, F. and Manke, T. (2016). deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 44: W160–W165.
Robinson, J. T., Thorvaldsdóttir, H., Winckler, W., Guttman, M., Lander, E. S., Getz, G. and Mesirov, J. P. (2011). Integrative genomics viewer. Nat. Biotechnol. 29(1): 24–26.
Rudenko, G. (1998). Selection for activation of a new variant surface glycoprotein gene expression site in Trypanosoma brucei can result in deletion of the old one. Mol. Biochem. Parasitol. 95(1): 97–109.
Touray, A. O., Rajesh, R., Isebe, T., Sternlieb, T., Loock, M., Kutova, O. and Cestari, I. (2023). PI(3,4,5)P3-dependent allosteric switch controls antigenic variation in trypanosomes. eLife 12: RP89331.
Ulbert, S., Chaves, I. and Borst, P. (2002). Expression site activation in Trypanosoma brucei with three marked variant surface glycoprotein gene expression sites. Mol. Biochem. Parasitol. 120(2): 225–235.

Supplementary information

The following supplementary information can be downloaded here:

Table S1. List of primer sequences indicated in the VSG-BarSeq protocol.