Development of a gene doping detection method to detect overexpressed human follistatin using an adenovirus vector in mice

Background Gene doping is the misuse of genome editing and gene therapy technologies for the purpose of manipulating specific genes or gene functions in order to improve athletic performance. However, a non-invasive detection method for gene doping using recombinant adenoviral (rAdV) vectors containing human follistatin (hFST) genes (rAdV) has not yet been developed. Therefore, the aim of this study was to develop a method to detect gene doping using rAdV. Methods First, we generated rAdV and evaluated the overexpression of the hFST gene, FST protein, and muscle protein synthesis signaling using cell lines. Next, rAdV was injected intravenously or intramuscularly into mice, and whole blood was collected, and hFST and cytomegalovirus promoter (CMVp) gene fragments were detected using TaqMan-quantitative polymerase chain reaction (qPCR). Finally, to confirm the specificity of the primers and the TaqMan probes, samples from each experiment were pooled, amplified using TaqMan-qPCR, and sequenced using the Sanger sequencing. Results The expression of hFST and FST proteins and muscle protein synthesis signaling significantly increased in C2C12 cells. In long-term, transgene fragments could be detected until 4 days after intravenous injection and 3 days after intramuscular injection. Finally, the Sanger sequencing confirmed that the primers and TaqMan probe specifically amplified the gene sequence of interest. Conclusions These results indicate the possibility of detecting gene doping using rAdV using TaqMan-qPCR in blood samples. This study may contribute to the development of detection methods for gene doping using rAdV.


INTRODUCTION
Gene therapy using gene editing technology is emerging as a new treatment strategy for various genetic and acquired diseases (Wang et al., 2019). Gene therapy drugs, One gene that can be exploited for gene doping is the human follistatin gene (hFST). In humans, FST is produced in the liver and secreted into the bloodstream by stimuli such as exercise (Hansen & Plomgaard, 2016). The secreted FST is delivered to skeletal muscles throughout the body where it binds to activin A and myostatin, inhibiting their action and suppressing TGF-β signaling (DePaolo, 1997;Patal, 1998). As a result, it alleviates Smad3-dependent suppression of Akt/mTOR/p70S6K signaling and induces skeletal muscle hypertrophy (Bodine et al., 2001). FST is on the WADA prohibited list (The World Anti-Doping Agency (WADA), 2020), and there are concerns that it may be misused for gene doping.
To detect gene doping, it is necessary to detect fragments of the transgene or vector. In a previous study, quantitative polymerase chain reaction (qPCR) was shown to be the easiest and the most specific method for detecting gene fragments (Sugasawa et al., 2019). Therefore, in this study, we developed a qPCR detection method. qPCR mainly uses SYBR-Green or TaqMan polymerases for amplification. A previous study reported that TaqMan-qPCR has higher sensitivity and specificity than SYBR-Green-qPCR (Aoki et al., 2020). Therefore, we developed a detection method using TaqMan-qPCR.
Therefore, the purpose of this study was to develop a method using TaqMan-qPCR to detect gene doping of human FST gene using the rAdV vector.

Cloning of recombinant adenovirus vector containing the human follistatin gene
The following plasmids were used in this study: Gen EZ TM ORF clone hFST in pcDNA3.1 (+) (GenScript, Piscataway, NJ, USA), pENTR4 (Thermo Fisher Scientific, Waltham, MA, USA), and pAd/CMV/V5-DEST (Thermo Fisher Scientific, Waltham, MA, USA). HEK 293A cells (Thermo Fisher Scientific) were used to clone and amplify rAdV vectors. Gen EZ TM ORF clone FST in pcDNA3.1(+) was used as a template to amplify the hFST gene with 5′-EcoRI and 3′-BamHI restriction sites using PCR. PCR product was inserted into the pENTR4 plasmid using T4 ligase (Promega, 92 Madison, WI, USA). The sequences of inserted hFST genes in the pENTR4 plasmids were read using sanger sequencing and confirmed to be correct sequences. The following experiments were performed in accordance with a previous study (Sugasawa et al., 2019).

CELL EXPERIMENTS Cell culture
HuH7 cells (RIKEN Bio Resource Center, Tsukuba, Japan) were cultured in Dulbecco's modified Eagle's medium (DMEM) (Thermo Fisher Scientific, Waltham, MA, USA) containing 10% fetal bovine serum (Wako, Osaka, Japan) and 1% penicillin/streptomycin (Wako), in an incubator maintained at 37 C with 5% CO 2 gas. HuH7 cell line established from a highly differentiated human hepatocellular carcinoma. HuH7 cells were used to mimic the high expression of FST in the liver. Cells were seeded in 12-well plates at a density of 1.0 × 10 5 cells and six-well plates at a density of 2.5 × 10 5 cells. After the cells reached 90% confluence, the medium was changed to serum-free DMEM, and the cells were infected by rAdV.
C2C12 cells (RIKEN Bio Resource Center, Tsukuba, Japan) were cultured under the same conditions. C2C12 is an immortalized mouse myoblast cell line. C2C12 cells were used to mimic the high expression of FST in the local muscle. After the cells reached 70-80% confluence, the medium was replaced with DMEM containing 2% horse serum and 1% penicillin/streptomycin (differentiation medium). Differentiation into myotubes was induced by culturing in the differentiation medium (changed every 2 days) for 7 days. After incubation, the medium was changed to serum-free DMEM, and the cells were infected by rAdV.

rAdV infection
HuH7 and C2C12 seeded in 12-well plates (for total RNA/DNA extraction) were infected with 1.0 × 10 10 vp, and those seeded in six-well plates (for total protein extraction) were infected with 4.0 × 10 10 vp of rAdV<hFST>. A group with the same amount of phosphate-buffered saline (PBS) and a group infected with the same amount of rAdV<mCherry> were considered as controls. Total RNA/DNA extraction and total protein extraction were performed 4 days after infection in HuH7 and 3 days in C2C12 cells.
Total RNA/DNA extraction Total RNA was extracted from culture cells. Sepasol RNA I Super G (Nacalai Tesque, Kyoto, Japan) was used for the extraction, as per the manufacturer's protocol. After extraction, the concentration of RNA in each sample was estimated using a microspectrophotometer (NanoDrop 1,000; Thermo Fisher Scientific, Waltham, MA, USA), and 500 ng of total RNA from each sample was reverse transcribed using 5 × PrimeScript RT Master Mix (Takara Bio, Shiga, Japan) to cDNA. The cDNA was diluted with a four-fold amount of purified water.
Total DNA was extracted from cell culture samples. A phenol/chloroform/isoamyl alcohol solution (Nacalai Tesque, Kyoto, Japan) was used for extraction, as per the manufacturer's protocol. After extraction, total DNA concentration of each sample was adjusted to 10 ng/µL.
The synthesized cDNA was used to measure transgene expression, and total DNA was used to detect the transgene fragments using SYBR-Green-qPCR.

Immunoblotting
Total proteins were extracted from the culture cells using lysis buffer (1% NP-40, 0.1% SDS, 20 mM Tris-HCl [pH 8.0], five mM ethylenediamine tetraacetic acid [EDTA], 150 mM NaCl, and proteinase inhibitor (Nacalai Tesque, Kyoto, Japan)). Lysates were centrifuged at 12,000×g for 15 min at 4 C. Total protein concentration for each sample was measured using the BCA protein assay kit (Takara Bio, Shiga, Japan), and five µg/lane (HuH7) or 10 µg/lane (C2C12) of total protein was used for gradient gel electrophoresis. For western blotting, the blots were incubated with primary antibodies. Table S1 lists the antibodies used in this study. Horseradish peroxidase-conjugated antirabbit IgG and anti-mouse IgG were used as the secondary antibodies. Signals were detected using a chemiluminescent reagent (EzWestLumi One; ATTO, Tokyo, Japan). Blots were scanned using a Light-Capture Cooled CCD Camera System (Image Quant LAS-4,000; GE Healthcare, Chicago, IL, USA).

Immunofluorescence
Immunofluorescence was performed 3 days after rAdV<hFST> infection of C2C12 cells. After removing the supernatant, cells were fixed with 4% paraformaldehyde (Wako) for 15 min at 20 C. Subsequently, the cells were permeabilized using 0.1% TritonX-100 (Nacalai Tesque) in PBS and blocked with 5% goat serum (Sigma-Aldrich, St. Louis, MO, USA) in PBS for 1 h at 20 C. The primary antibody against Follistatin (Proteintech Group, Chicago, IL, USA, 60060-6-lg, 1:300) was diluted in 5% goat serum in PBS and incubated with the samples for 1 h at 20 C. Wells were washed with PBS three times and incubated with the secondary antibody (Alexa Fluor 488 anti-mouse IgG; Jackson ImmunoResearch Laboratories, Inc., West Grove, PA, USA, 1:400) for 1 h at 20 C. Finally, the wells were washed with PBS three times, stained with Dapi-Fluoromount-G Ò (SouthernBiotech, Birmingham, AL, USA), and subjected to fluorescence microscopy.

ANIMAL EXPERIMENTS Animal
The animal experiments conducted in this study were approved by the Animal Experiment Committee of the University of Tsukuba (Approval Number: 20-378). A total of 6-week-old male ICR mice were purchased from Central Laboratories for Experimental Animals (Tokyo, Japan) and subjected to a 1-week acclimation period in a cage (Maximum of five mice in one cage). Mice were housed in an air-conditioned, pathogen-free animal room with a 12/12 h light/dark cycle. The mice were allowed to consume normal solid food and water ad libitum. At the beginning of the experiment, the mice weighed 34.1~37.4 g and were 7 weeks old. If weight loss was less than 20% during the experiment, or if there was obvious injury or illness, euthanasia was to be performed. There were no mice that corresponded to the above in this experiment. The sample size for each experiment was determined with previous study (Sugasawa et al., 2019).

Long-term detection of transgene fragments for intravenous injection
Approximately, 50 µL of whole blood was collected from the tail vein of 10 untreated mice aged 7 weeks into a 1.5-mL tube containing 150 µL of EDTA-disodium salt (EDTA-2Na)/PBS mixture, which was used as a pre-injection sample. A total of 4 days later, 2.0 × 10 11 vp of rAdV<hFST> was injected through the orbital venous plexus. Thereafter, approximately 50 µL of whole blood was collected from the tail vein at 3, 6, and 12 h, and 1, 2, 3, 4, 5, 6, and 7 days after injection.
Total DNA was extracted using a phenol/chloroform/isoamyl alcohol solution (Nacalai tesque, Kyoto, Japan). After extraction, the concentration of DNA each sample was adjusted to 10 ng/µL, and DNA was subjected to TaqMan qPCR to detect hFST and CMVp gene fragments.

Detection of transgene fragments in each specimen for intravenous injection
Fourteen untreated mice aged 7 weeks were randomly divided into two groups: group 1 (G1) (n = 6) and group 2 (G2) (n = 8). G1 was injected with 100 µL of PBS, and G2 was injected with 100 µL of 2.0 × 10 12 vp/mL rAdV<hFST> through the orbital venous plexus. A total of 6 h after injection, the mice were administered general anesthesia by isoflurane inhalation, followed by whole blood collection with EDTA-2Na. The mice were euthanized by cervical dislocation. Whole blood was centrifuged at 5,000 rpm for 10 min at 4 C, and plasma and blood cell fractions were collected separately. Each sample was stored at −20 C until the next analysis.
Total DNA was extracted from each sample using NucleoSpinÒ cfDNA XS (Macherey-Nagel, Düren, Germany) for plasma fractions and NucleoSpinÒ Blood (Macherey-Nagel, Düren, Germany) for blood cell fractions, as per the manufacturer's protocol. After extraction, the concentration of DNA each sample was adjusted to 10 ng/µL. The adjusted DNA was used to detect hFST and CMVp gene fragments using TaqMan qPCR.

Determination of specificity of the primers and TaqMan probe
The DNA samples from each specimen (whole blood, plasma, and blood cells) mentioned in "Long-term detection of transgene fragments for intravenous injection" and "Detection of transgene fragments in each specimen for intravenous injection" were pooled in a 1.5-mL tube. The pooled DNA was subjected to TaqMan qPCR to detect the transgene fragment using the same primers and TaqMan probes as in "Long-term detection of transgene fragments for intravenous injection" and "Detection of transgene fragments in each specimen for intravenous injection". After the TaqMan qPCR, the amplified product was collected in a 1.5-mL tube and purified with NucleoSpinÒ Gel and PCR Clean-up (Macherey-Nagel). The purified and amplified products were sequence by an external vendor (FASMAC, Kanagawa, Japan) using the Sanger Sequence Analysis Service. The obtained data were analyzed using CLC Sequence Viewer ver. 8.0 (QIAGEN, Hilden, Germany) and BioEdit ver. 7.2.5 (developer: Tom Hall).

Long term detection of transgene fragments for intramuscular injection
Approximately, 50 µL of whole blood was collected from the tail vein of six untreated mice aged 7 weeks into a 1.5-mL tube containing 150 µL of EDTA-2Na/PBS mixture, which was used as a pre-injection sample. A total of 3 days later, 2.5 × 10 12 vp/mL of rAdV<hFST> was injected at 30 µL into the tibialis anterior and at 50 µL into the gastrocnemius. Thereafter, approximately 50 µL of whole blood was collected from the tail vein at 3, 6, and 12 h, and 1, 1.5, 2, 2.5, 3, 3.5, 4, and 5 days after injection.
Total DNA was extracted using a phenol/chloroform/isoamyl alcohol solution. After extraction, the concentration of DNA in each sample was adjusted to 10 or 100 ng/µL, and DNA was subjected to TaqMan qPCR to detect hFST and CMVp gene fragments.

PRIMER DESIGN AND QPCR
The sequences of primers and TaqMan probes are listed in Table S2. Primers and TaqMan probes were synthesized by Integrated DNA Technologies (Coralville, IA, USA).
SYBRÒ FAST qPCR Master Mix (KK 4602; Kapa Biosystems, Wilmington, MA, USA) reagent was used to perform SYBR-Green qPCR. The volume of the template solution, final concentration of the primer, negative control, and experimental conditions were taken from a previous study (Aoki et al., 2020).
PrimeTime Gene Expression Master Mix (Integrated DNA Technologies, Coralville, IA, USA) reagent was used to perform TaqMan-qPCR. The volume of template solution, final concentrations of primer and probe, negative control, and experimental conditions were taken from a previous study (Aoki et al., 2020). For standard curve for absolute quantification, 100 pg/µL Gen EZ TM ORF clone FST in pcDNA3.1(+) was used. The coefficient of determination (R 2 ) of the calibration curve was set to R 2 > 0.98.
Both SYBR-Green qPCR and TaqMan-qPCR samples were run in duplicates, and QuantStudio 5 Real-Time PCR Systems (Thermo Fisher Scientific, Waltham, MA, USA) were used.

STATISTICAL ANALYSIS
GraphPad Prism version 7.04 (GraphPad, Inc., La Jolla, CA, USA) was used for statistical analysis of the data. First, the Shapiro-Wilk normality test was performed on all data to check the normality of the distributions. Then, data following a normal distribution were subjected to an unpaired t-test or one-way analysis of variance test, and the Tukey-Kramer method was used for the post hoc test. For non-normal distribution, Mann Whitney's U test and Benjamini, Krieger and Yekutiel's two-stage method were used. p-values less than 0.05 were considered statistically significant. Cell experiments data were indicated mean ± SD. Animal experiments data were indicated mean ± SEM.

Transgene fragments were detected for 4 days after intravenous injection
Compared with the pre-infection samples, hFST transgene fragments were detected until 2 days (p < 0.05) and CMVp until 4 days (p < 0.05) after injection. In addition, both gene fragments were the highest at 6 h after injection (Fig. 3).

Transgene fragments were detected in each specimen of intravenous injection
No transgene fragments were detected in any of the specimens in G1. In contrast, transgene fragments, including of hFST and CMVp, were detected in both specimens in G2. Moreover, when plasma and blood cell fractions were compared, transgene fragments were significantly more localized in the plasma fraction (p < 0.05) (Fig. 4).

Primers and TaqMan probes were specific for the targeted transgene fragment
The results of the amplification curve of TaqMan-qPCR confirmed that the hFST and CMVp gene fragments in whole blood, plasma, and blood cell samples were fully amplified. Sanger sequencing showed that the sequences of all samples matched the reference sequences by more than 90 bases (Fig. 5). This confirmed that amplification by the prepared primers and TaqMan probe was specific.

DISCUSSION
In this study, we evaluated gene and protein expression and muscle protein synthesis signaling using rAdV<hFST> and developed a detection method for gene doping by intravenous or intramuscular injection of rAdV<hFST>. In cell experiments, the hFST gene and FST protein expression in HuH7 and C2C12 cells significantly increased after injecting rAdV<hFST>, suggesting that the hFST gene and FST protein expression in the liver or muscle might significantly increase after intravenous or intramuscular injection of rAdV<hFST>. In addition, because FST is secreted extracellularly, the significantly raised FST in the liver or muscle would be secreted into the blood or near the muscle, respectively. In addition, because the phosphorylation of Akt and p70S6K significantly increased in C2C12 cells, it is expected that muscle synthesis would have been activated to induce muscle hypertrophy. The results of the cell experiments in this study suggested that intravenous or intramuscular injection of rAdV<hFST> may be used to construct a gene-doping model. We injected rAdV<hFST> in mice intravenously or intramuscularly and then sampled blood over time to see how many days after injection it took for the transgene and rAdV vector fragments to disappear. The results showed that gene fragments could be detected up to 4 days after intravenous injection (Fig. 3) and up to 3 days after intramuscular injection (Fig. 6). To date, no study has attempted to detect the hFST gene fragment, and these results are considered important findings for future research. However, further studies are needed to clarify the detectable period of rAdV<hFST> in humans. It has been reported that drug metabolism in small experimental animals such as mice is about 10 times higher than that in humans (Kato, 1981), and it is possible that the rate of metabolism of genes introduced by adenovirus vectors differs between experimental animals and humans, but the extent of the difference is not yet clear.
In addition, compared with the amount of virus used in humans (Zhang et al., 2018), this study used approximately 1/10 for intravenous injection and approximately 1/50 for intramuscular injection. Considering the body weight ratio (2,000:1) and circulating blood volume ratio (1,500:1) of humans (approximately 60 kg) and mice (approximately 30 g), a large amount of virus was administered in this study. Doping tests require highly sensitive detection methods. To determine localization of the transgene fragments in whole blood, whole blood was collected and centrifuged into plasma and blood cell fractions and analyzed to determine which fraction contained the highest number of gene fragments. The hFST and CMVp gene fragments were localized more in the plasma fraction than in the blood cell fraction (Fig. 4). In contrast, in a previous study using rAdV, gene fragments were more localized in the blood cell fraction than in the plasma fraction (Sugasawa et al., 2019). This may be due to the time between injection and blood collection. Blood samples were collected 5 days after injection in the previous study, whereas blood was sampled 6 h after injection in this study. Because the injected rAdV vector was used as a viral solution, it is presumed to have localized in the plasma immediately after injection. It infected red blood cells or was phagocytosed by leukocytes over time. This is expected to lead to the development of more sensitive detection methods by determining the exact time points at which gene fragments localize predominantly in the plasma or in the blood cell fraction and which of the blood cell components (erythrocytes, leukocytes) or the surfaces of blood cells have the most localization.
Doping tests should have a high specificity. Therefore, it is necessary to develop a method to specifically detect target gene fragments. We used the Sanger sequencing method to verify the specificity of the detection performed in this study. As a result, we were able to specifically detect transgene fragments (Fig. 5). Furthermore, no non-specific amplification or contamination was found in the control group or samples before rAdV vector injection. These results indicate that the DNA extraction method, primer and TaqMan probe design, thermal cycling conditions, and reagents for TaqMan-qPCR in this study were accurate and optimal. Therefore, the protocol constructed in this study may be directly applicable to human specimens.
A limitation of this study was that we did not construct a gene doping model by intravenous or intramuscular injection of rAdV<hFST>. In other words, we were not able to establish a model that shows specific phenotypes of intravenous or intramuscular injection of rAdV<hFST>, such as changes in muscle strength, muscle wet weight, and muscle fiber type, owning to the increased levels of hFST gene and FST protein in the liver and muscle and activation of muscle synthesis. Therefore, it is necessary to establish a gene-doping model by intravenous or intramuscular injections in the future.
Research on methods for detecting gene doping has been conducted for only about 20 years, and it is still a developing research field with many unknowns. This study is a pioneering work in this field, and the results of this study may contribute to the development of detection methods for gene doping using rAdV<hFST>.

CONCLUSIONS
In this study, we aimed to develop a method to detect multiple transgene fragments as proof of gene doping using rAdV<hFST>. Figure 7 presents a summary. This study showed that in the rAdV<hFST> transgenic model, multiple transgene fragments could be detected using TaqMan-qPCR from as little as 50 µL of whole blood, and each gene fragment was localized predominantly in the plasma. The new findings of this study may contribute to the development of detection methods for gene doping using rAdV<hFST>.