Genotyping results from skin swabs
Over nine different sessions, we processed 470 skin swab samples of wild type and two different genetically altered lines (Table 1). Six different animal technicians obtained the swab samples across all sessions. The overall success rate of obtaining a definitive genotyping result was 86.3%.
During session 1, all the fish were swabbed and fin clipped to verify the correct genotype results we obtained. There were no errors. The total length of the fish swabbed in this session ranged from 21.53 mm to 34.28 mm.
Session 4 was performed in a small quarantine room with limited area flow. We believe this impacted the success rate based on similar error rates with fin clips collected in the same room at later dates. To rectify this issue, we added the 60-min drying time to the rest of the sessions. Removing this data point raises the average success rate of obtaining a definitive genotyping result to 90.9% (Table 2).
From session 5, a 60-min drying time was added to this procedure. The average success rate of obtaining a definitive genotyping result from sessions 5 to 9 is 90.8% (Table 3).
During sessions 8 and 9, the same technician collected all the swab samples. The average success rate across these two sessions is 92.1% (Table 4).
Genotyping results determination
The reactions were run on a thermocycler, and the data were analyzed using the thermocycler software. The cycle thresholds from that software are then used to calculate a delta Ct using both the target and control Cts:
ΔCt = Ct (gene of interest) - Ct (housekeeping gene) = RCN
The results for each target assay are then compared with RCN ranges previously established for the presence or absence of the alleles based on zebrafish tissue samples. These calculations and results review occur in a proprietary software system developed by the automated genotyping provider.
The graphs in Figure 1A–C show the curves produced by the on-board capture and analysis software of the PCR machines, depicting the visual amplification curve for the PCR assay, which measures changes in fluorescence over the PCR cycles.

Figure 1: The visual amplification curve for the PCR assay run for comparing skin swabs versus fin clips. (A) Comparison cycle of threshold (Ct) of skin swabs to fin clips for internal control assay. Each line represents the amplification curve of a sample (each sample was run in duplicate). (B) Comparison cycle of threshold (Ct) of skin swabs to fin clips for Trp53 mutant allele assay. Each line represents the amplification curve of a sample (each sample was run in duplicate). Blue lines indicate the mutant allele and green lines the wildtype allele. (C) Comparison cycle of threshold (Ct) of skin swabs to fin clips for Zfp36l1b mutant allele assay. Each line represents the amplification curve of a sample (each sample was run in duplicate). Blue lines indicate the mutant allele and green lines the wildtype allele. Fin clips exhibit earlier Ct and a tighter range of values. Please click here to view a larger version of this figure.
| Batch # | Total Samples | Successful | UD2 | UD4 | % Success | Notes |
| 1 | 15 | 15 | 0 | 0 | 100.0 | Results compared to fin clip from same animal. Body lengths ranged from 21.53 mm to 34.28 mm. |
| 2 | 54 | 44 | 9 | 1 | 81.5 | |
| 3 | 77 | 71 | 6 | 0 | 92.2 | |
| 4 | 71 | 35 | 36 | 0 | 49.3 | Higher error rate due to location sampling took place. |
| 5 | 50 | 39 | 11 | 0 | 78.0 | |
| 6 | 71 | 65 | 6 | 0 | 91.5 | |
| 7 | 29 | 29 | 0 | 0 | 100.0 | |
| 8 | 82 | 73 | 9 | 0 | 89.0 | Technician G swabbed all fish. |
| 9 | 21 | 20 | 0 | 1 | 95.2 | Technician G swabbed all fish. |
| 470 | | | | | |
| | | | Average | 86.3 | |
Table 1: Summaries of results, along with calculated averages. Summary of nine swabbing sessions held and how many errors (UD2 and UD4) occurred. Success was determined by whether a result could be reported. UD2 samples did not pass the Transnetyx internal control requirements for sample integrity. These samples will have been tested four times, but there was not sufficient DNA to determine a result. UD4 samples did not return a result, possibly due to contamination of the sample.
| Batch # | Total Samples | Successful | UD2 | UD4 | % Success | Notes |
| 1 | 15 | 15 | 0 | 0 | 100.0 | Results compared to fin clip from same animal. Body lengths ranged from 21.53 mm to 34.28 mm. |
| 2 | 54 | 44 | 9 | 1 | 81.5 | |
| 3 | 77 | 71 | 6 | 0 | 92.2 | |
| 5 | 50 | 39 | 11 | 0 | 78.0 | |
| 6 | 71 | 65 | 6 | 0 | 91.5 | |
| 7 | 29 | 29 | 0 | 0 | 100.0 | |
| 8 | 82 | 73 | 9 | 0 | 89.0 | Technician G swabbed all fish. |
| 9 | 21 | 20 | 0 | 1 | 95.2 | Technician G swabbed all fish. |
| 399 | | | Average | 90.9 | |
Table 2: Results from 9 sessions except session 4. Results from session 4 were removed because the sample collection took place in a different room with less circulation, which is hypothesized to have impacted the error rate.
| Batch # | Total Samples | Successful | UD2 | UD4 | % Success | |
| 5 | 50 | 39 | 11 | 0 | 78.0 | |
| 6 | 71 | 65 | 6 | 0 | 91.5 | |
| 7 | 29 | 29 | 0 | 0 | 100.0 | |
| 8 | 82 | 73 | 9 | 0 | 89.0 | Technician G swabbed all fish. |
| 9 | 21 | 20 | 0 | 1 | 95.2 | Technician G swabbed all fish. |
| 253 | | | Average | 90.8 | |
Table 3: Sessions 5–9. A 60-min drying period was added to sessions 5–9. The improved average success rate calculated here demonstrates the effectiveness of a longer drying period.
| Batch # | Total Samples | Successful | UD2 | UD4 | % Success | |
| 8 | 82 | 73 | 9 | 0 | 89.0 | Technician G swabbed all fish. |
| 9 | 21 | 20 | 0 | 1 | 95.2 | Technician G swabbed all fish. |
| 103 | | | Average | 92.1 | |
Table 4: Sessions 8 and 9. The same technician swabbed all the fish in sessions 8 and 9, eliminating the variable of multiple technicians.
Supplementary Table 1: List of housekeeper and reporter genes used in this study. Please click here to download this file.