Development of a rabies virus-based retrograde tracer with high trans-monosynaptic efficiency by reshuffling glycoprotein
Molecular Brain volume 14, Article number: 109 (2021)
Rabies virus (RV) is the most widely used vector for mapping neural circuits. Previous studies have shown that the RV glycoprotein can be a target to improve the retrograde transsynaptic tracing efficiency. However, the current versions still label only a small portion of all presynaptic neurons. Here, we reshuffled the oG sequence, a chimeric glycoprotein, with positive codon pair bias score (CPBS) based on bioinformatic analysis of mouse codon pair bias, generating ooG, a further optimized glycoprotein. Our experimental data reveal that the ooG has a higher expression level than the oG in vivo, which significantly increases the tracing efficiency by up to 12.6 and 62.1-fold compared to oG and B19G, respectively. The new tool can be used for labeling neural circuits Therefore, the approach reported here provides a convenient, efficient and universal strategy to improve protein expression for various application scenarios such as trans-synaptic tracing efficiency, cell engineering, and vaccine and oncolytic virus designs.
Mapping neural circuits is a prerequisite to elucidate the mechanisms of brain functions. Several neurotropic viruses can infect neurons and spread across synapses between neurons, playing important roles in depicting the neurocircuit . Among these virus-based tools, rabies virus (RV), encoding five proteins (N, P, M, G and L), is the most popular trans-monosynaptic viral tool to map the direct input networks of specific types of neurons in specific brain regions [2,3,4,5,6]. To realize cell type specific and trans-monosynaptic tracing, (i) the RV-G gene is deleted from genome (RV-delG); (ii) a chimeric EnvA gene is created by swapping the cytoplasmic part of G and wild type EnvA (a glycoprotein from ASLV-A; (iii) a cell line expressing the chimeric EnvA is prepared to package the genome of RV-delG to produce the pseudo-typed virus (EnvA-RV-delG); (iv) two adeno-associated virus (AAV) vectors (AAV-G, and AAV-TVA-GFP), cre-dependently delivering G-gene, GFP, and TVA (a receptor of EnvA), are constructed separately; (v) EnvA-RV-delG can enter specific neurons infected by AAV-TVA-GFP; (vi) in the neurons co-infected by EnvA-RV-delG/AAV-G/AAV-TVA-GFP, the G protein trans-complements the package of RVG-RV-delG which is capable of trans-monosynaptically retrograde spreading to input neurons [2, 3]. However, current RV versions can only label a fraction of presynaptic neurons . When AAV-G infects neurons, the G expression cassette cannot be replicated because of the AAV features (or the template number for G-gene is fixed), while EnvA-RV-delG enters the same neuron, the G-deleted RV genome can be reproduced and expressed according to the RV life cycle.
Therefore, a significant mismatch between the amounts of G protein and other RV proteins might exist, leading to few active RV particles and thus low tracing efficiency, at the same time, generating excessive unused viral proteins and thus high cytotoxicity. Conventional methods for improving the protein expression level are to use stronger promoters, optimize the gene sequence, or provide more copies of the gene [8,9,10,11,12,13]. Indeed, previous reports have provided evidence that elevating the expression level of the G protein improves the tracing efficiency [14, 15]. The protein synthesis machinery has a preference for how to better translate proteins in different species. Codon pair bias (CPB) is a common phenomenon in various species for adjacent codon pairs, which appears with different frequencies in different species . The favorite or unfavorite of a codon pair can be defined by its codon pair bias score (CPBS), defined as the ratio of observed frequency to the expected frequency. Every codon pair has a CPBS. A positive CPBS means that a codon pair is overrepresented, which may be preferred by the organism, whereas a negative CPBS may be unfavorable. Therefore, recoding a gene with positive CPBS may adjust protein expression. Many viruses have been attenuated by using the negative CPBS which results in decreased protein production [17,18,19], while only limited number of studies explored the effect of positive CPBS on viral properties. RV as a retrograde trans-monosynaptic tracer is most commonly used to label neural circuits in mice, however, mice are not the reservoir host. This cross-species shift may have a profound impact on the expression of proteins encoded by RV. The amount of G is a limiting factor for trans-synaptic efficiency of RV. Here, we propose that better conditions to increase the G protein level could be achieved by using positive CPBS in mice. Indeed, this convenient method can elevate the G expression level and improve RV trans-monosynaptic efficiency.
The relationship of canine, mouse and human in CPB
RV has wide infection range, such as canine, mouse and human. Among these animals, canine is one of the RV reservoir hosts, which providing the opportunity that RV has adopted the CPB of canine during the long-term co-evolution. While, RV-based tools are usually used in mouse model in neuroscience research. Therefore, the codon pair score (CPS) of canine, mouse and human were calculated using previously reported method . Briefly, a total of 45,094, 68,272 and 110,788 genes for canine, mouse and human were analyzed to calculate codon pair bias score for each of the 3721 codon pair (61 × 61 codons). The CPS values of these three species are very similar (Fig. 1a-c), however, the relationships of the three species are different. The correlation coefficient value of codon pair preference between human and canine (Spearman rho = 0.9859) is closer than that between mouse and canine (Spearman rho = 0.9459), and also closer than that between human and mouse (Spearman rho = 0.9588). Next, we calculated the CPBS values for each gene of mouse and canine respectively and plotted the CPBS value against its gene length. We found that the majority of genes in mouse and canine have positive CPBS, averaged at 0.0651 and 0.0704, respectively, which are similar to the human 0.0703 (Additional file 1: Fig. S1). Therefore, these data provide a basic for oG reshuffling (Fig. 2a).
ooG has higher expression level than oG
The oG, the best version of G in current RV-based tools , was selected as a model for codon pair optimization based on codons pairs preference data from the mouse genome to improve its expression level. Based on the CPS obtained from the mouse, we attempted to reshuffle the oG sequence to improve oG’s expression level, and named the new oG sequence as ooG (Additional file 1: Fig. S2). The median CPBS of oG and ooG is -0.0738, 0.3114, respectively (Fig. 2b, c). The number of CpG and UpA are 73 and 45 for oG, 38 and 34 for ooG (Fig. 2c).
Next, the genes of the ooG and oG were commercially synthesized and cloned into the self-complementary AAV (scAAV) vector, respectively , in which the ooG or oG was expressed by the syn promoter (Fig. 3a). 150 μl of scAAV-hsynP-ooG-bGHpA and scAAV-hsynP-oG-bGHpAwere injected into the ventral hippocampus (vHPC) region of C57BL/6 mouse, respectively. After three weeks, the brain tissue was collected and analyzed (Fig. 3b). Firstly, protein level was determined using antibody against G, which is specific to the ooG and the oG (Fig. 3c). We found that the amount of ooG protein is higher than oG’s (Fig. 3d). Furthermore, we found that the ooG mRNA copies are higher than oG’s (Fig. 3e). Collectively, the ooG as a novel gene has higher protein expression level than oG.
ooG can enhance the trans-synaptic efficiency
Based on the data, we provide a hypothesis that the ooG has a potential ability of increasing RV trans-synaptic efficiency. To directly compare the transsynaptic efficiency of the ooG, oG, and B19G, the long-distance input data in known circuit were collected and analyzed by calculating the convergence index. ssAAV-EF1α-DIO-EGFP-F2A-TVA-WPRE-bGHpA was co-injected with one of the three AAVs: scAAV-hsynP-ooG-bGHpA, scAAV-hsynP-oG-bGHpAand scAAV-hsynP-B19G-bGHpA into the Vhpc region of Thy-1 cre mouse, respectively. After three weeks, the EnvA-pseudorabies (Enva-RV844) were injected into the same region. The brains were collected and processed seven days later.
We observed the main long-distance input regions are the paraventricular thalamic nucleus (PVA) and the medial septal nucleus (MS), which are known and consistent with the previous reports (Fig. 4a) . Then the convergence index was calculated to directly compare the retrograde trans-synaptic efficiency (14) of these three G proteins by using the number of input neurons (mRuby3 positive) divided by the number of starter cells (EGFP and mRuby3 positive). In PVA region, the convergence index is the 2.547 ± 0.05, 1.132 ± 0.05 and 0.041 ± 0.01 for ooG, oG and B19G, respectively (Fig. 4b). In MS region, the convergence index is the 1.017 ± 0.04, 0.320 ± 0.02 and 0.061 ± 0.01 for the ooG, oG and B19G, respectively (Fig. 4c). Comparison with the data from the ooG, oG and B19G, we found that the convergence index of ooG is 2.3-fold and 62.1-fold higher than the oG and B19G, and the convergence index of oG is 27.6-fold higher than the B19G in PVA region which is similar with previous report . The convergence index of the ooG is 3.2-fold and 16.7-fold higher than the oG and B19G, and the convergence index of the oG is 5.2-fold higher than the B19G in MS region. In addition, the convergence index is the 0.113 ± 0.006 and 0.009 ± 0.002 for the ooG and oG in lateral preoptic area (LPO), and the convergence index of the ooG is 12.6-fold higher than the oG (Fig. 5).
In addition, we compared the effects of the ooG and oG on labeling specificity and cellular toxicity. Firstly, the whole-brain connections to the vHPC were analyzed to investigate the labeling specificity, which revealed that the same brain regions were labeled in the ooG and oG groups (Additional file 1: Fig. S3a). These results indicated that the labeling specificity of the ooG is consistent with the oG. Secondly, the cellular toxicity of the ooG and the oG were determined by detecting the caspase-3 and cell morphology. These results showed that the signals of caspase-3 is not obvious in the ooG and oG groups (Additional file 1: Fig. S3b), and the cell morphologies in ooG and oG groups are similar (Additional file 1: Fig. S3c).
To label the complete pre-synaptic and post-synaptic partners of a specific neuron is the ultimate goal for the neural circuit tracers. RV as the most important and popular tool has been used widely. However, the current systems can only label a fraction of presynaptic neurons . Previous studies showed that the amount of the G protein is one of the key steps for adjusting the trans-synaptic efficiency . Therefore, many efforts have been taken for elevating the expression level of G protein. The direct method is using the stronger promoter to enhance the protein amount. Various promoters have different ability of improving gene expression in various cell types [22, 23]. Miyamichi et al. select CAG as a stronger promoter to drive G expression . Screening different G proteins from various rabies virus strains is another method for improving RV trans-synaptic tracing efficiency. Kim et al. compared various G (N2C strain and Pasteur strain) and engineered a chimeric G (oG) using codon-optimized strategy . In addition, providing more copies of the G gene might be a simple method for increasing gene expression. While the package capacity of AAV is limited. Therefore, screening and optimizing G gene is an appropriate method for increasing the gene expression to improve RV trans-synaptic tracing efficiency.
During the long-term evolution, virus has adapted its host’s environment for reproducing. All mammals are susceptible to rabies virus infection, only a few species are its reservoirs. Therefore, virus gene codons and gene codon pairs are adaptive to its reservoirs. Canine is one of RV reservoirs, while mouse is not. Therefore, the replicating and packaging efficiency of RV is not optimal in mouse. Mouse is the most important animal model for neuroscience research. Reshuffling virus genes according to the mouse genetic context might provide the proper condition for virus life cycle. CPB is a stable characteristic of a species. Based on the CPB strategy, attenuated viruses have been developed using under represented codon pair to deoptimize gene in RNA viruses and DNA viruses [17,18,19, 24,25,26]. In these cases, deoptimized gene has rare codon pairs and more CpG and UpA, which might cause the lower gene expression. Several mechanism of virus attenuation were appeared, (i) under represented codon pair results in decreasing efficiency of translation, such as translation elongation rate, translation initiation complex dissociation, and protein processing [17, 19, 24, 26, 27], (ii) more CpG and UpA result in inducing an innate immune response and reducing mRNA stability [28,29,30,31]. However, the exact mechanism of virus attenuation is still disputed and needs to be determined.
Unlike these previous studies, we provided a hypothesis that optimize gene using overrepresented codon pair to improve protein expression according to the mouse genome characteristic. Indeed, the ooG protein expression amount is higher than the oG (Fig. 3d, e). Therefore, ooG can significantly increase the RV trans-synaptic efficiency (Fig. 4a–c). The characteristics of the oG and ooG in mouse are that the CPBS is − 0.0738 and 0.3114, the number of CpG is 73 and 38, and the number of UpA is 45 and 34, respectively (Fig. 2c). Based on the results and possible mechanism for virus attenuated, we proposed that the over represented codon pair ooG provides a suitable situation for elevating translation efficiency, reducing an innate immune response and helping mRNA stability. However, the real mechanism should be determined in future. Collectively, whatever the mechanism for CPB, it has become increasingly clear that CPB has a profound impact on protein expression.
In the present study, using RV as a model target and trans-synaptic tracing efficiency as a readout, we verified our hypothesis that nucleic acid sequence optimization based on CPB analysis is a convenient and efficient strategy to improve protein expression level. The ooG will be a useful tool for mapping neural circuits. Importantly, the approach reported here might provide a convenient, efficient and universal strategy to improve protein expression for various application scenarios such as other tracers, cell engineering, vaccine and oncolytic virus designs.
Materials and methods
All procedures were approved by the Animal Care and Use Committees at the Shenzhen Institute of Advanced Technology or Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences. Eight-week-old male C57BL/6 mice were used for testing the protein expression and mRNA level. Eight-week-old Thy1-Cre mice were used for analyzing the retrograde trans-synaptic efficiency of rabies virus with the trans-complementary with various version G proteins, which expresses Cre in the nervous system of transgenic mice.
Rearrangement of the oG
Firstly, the codon pair bias of mouse, canine and human were calculated using the previous method  based on the gene CDS annotation of each species from the Ensembl database release 98, which were based on the genome assembly GRCm38, CanFam3.1 and GRCh38.p13, respectively. Secondly, the oG gene was optimized based on the mouse CPS according to the general rule of maintaining the mean free energy of the RNA folding within a suitable range and to avoid dramatic changes in its secondary structure.
The algorithm to calculate CPS and CPBS for coding sequence was developed after the description of Kunec et al.  and Coleman et al. . For a specific codon pair, the CPS is defined as the natural logarithm of the result of the observed codon pairs number divided by the expected number of codon pairs in all coding sequences of a species. The observed number is the actual number of the codon pair in the coding sequences. And the expected number of a codon pair is a theoretical value based on the proportion of amino acid and codon of the codon pair.
The N(codon1) and N(codon2) denote the number of occurrences of two codons in a codon pair, respectively. The N(aa1) and N(aa2) are the number of corresponding amino acids. And N(aa1aa2) is the number of the amino acid pair for the codon pair above. A positive CPS suggested the codon pair is over represented in the species, whereas a negative value is under represented.
Then, the CPBS for a coding sequence is the result of the arithmetic mean of all codon pair CPS values.
Preparation AAV and RV
scAAV-hsynP-ooG-bGHpA, scAAV-hsynP-oG-bGHpA, and scAAV-hsynP-B19G-bGHpA were constructed by inserting the expression cassette into the self-complementary AAV (scAAV) core vector . ssAAV-EF1α-ooG-WPRE-bGHpA and ssAAV-EF1α-DIO-EGFP-F2A-TVA-WPRE-bGHpA were constructed by inserting the expression cassette into the single-strand AAV2 core vector (addgene 20298). All plasmids were confirmed by sequencing. Then these recombinant AAVs were packaged with the AAV9 capsid, respectively.
Previous report shows that the mRuby3 is a new version red fluorescent protein with high lighter and stabilization . Therefore, the mRuby3 gene was selected and inserted into pSADdeltaG-F3 vector . Then the plasmid was confirmed by sequencing. The EnvA-RV-mRuby3 (EnvA-RV844) was prepared using the previous method .
Stereotaxic microinjection in mice brains
Briefly, mice were anesthetized with chloral hydrate before injection. scAAV-hsynP-ooG-bGHpA (6.9 × 1011 vg/ml), scAAV-hsynP-oG-bGHpA (1.1 × 1012 vg/ml) and scAAV-hsynP-B19G-bGHpA (2 × 1012 vg/ml) were respectively mixed (1:1, total volume 80 nl) with ssAAV-EF1α-DIO-EGFP-F2A-TVA-WPRE-bGHpA (6.2 × 1011 vg/ml) and co-injected into the vHPC region of mice. Then, the EnvA-RV-mRuby3 (EnvA-RV844, 2 × 108 FFU/ml) was injected into the same region after three weeks. During the process of virus injection, syringes were left in place for 10 min following injections to minimize diffusion. All surgical mice were back to the animal facility for extra seven days for fluorescent protein expression and RV trans-synaptic spreading.
Brain section and immunohistochemistry
Brain sections were collected referred in previous studies . Briefly, mice were anaesthetized and transcardially perfused with 0.9% saline followed by 4% paraformaldehyde solution. The brains were removed and post-fixed over-night in 4% paraformaldehyde before being sectioned into 50 µm slices, and slices were stained with DAPI. For staining the caspase-3, fixed slices were immunostained with the caspase-3 antibody (1:500, Cell Signaling, #9661) and amplified with the goat anti-rabbit secondary antibody (1:500, Jackson, #611-605-215), and slices were stained with DAPI. Slices were imaged using the Leica TCS SP8 confocal microscope or the Olympus VS 120 slide scanning system.
Preparation of antibody against the ooG
BALB/c mice were immunized directly by intramuscular injection of 200 μl of AAV mixture with phosphate buffered saline, which contains 20 μl of ssAAV-EF1α-ooG-WPRE-bGHpA (1.4 × 1012 vg/ml). Two weeks later, the second immunization was performed by intramuscular injection of the same amount AAV and the sera were collected after two weeks. Then, the specificity of the antibody was confirmed by using western blot to detect the expression of the ooG and the oG in the 293 T cells, which were separately transfected with 2 μg of plasmids containing the expression cassettes of the ooG (pAAV-Ef1a-ooG-WPRE-BGHpA) and the oG (pAAV-Ef1a-oG-WPRE-BGHpA).
150 nl of scAAV-hsynP-ooG-bGHpA (8.7 × 1011 vg/ml) and scAAV-hsynP-oG-bGHpA (1.3 × 1012 vg/ml) were injected into the vHPC region of C57BL/6 mouse, respectively. After three weeks, the tissue was collected and treated, one part for RNA extraction, and another for western blot. The tissue lysate or 293 T cells were separately analyzed on a 12% SDS-PAGE, and then electro-transferred to PVDF Immobilon-P membranes (Millipore), blocked with 5% skim milk in TBST, and then treated with antibody against ooG at 1:500 dilution for 1 h at room temperature. After washing 3 times in TBST, the secondary anti-mouse antibody conjugated to Horseradish peroxidase at 1:5000 dilution was applied to the blots for 1 h at room temperature. To quantify the expression amount of GAPDH, the protein was analyzed by using fisrt antibody (1:2000) and second antibody (1:5000). Signal was detected with ECL western blotting reagent.
Quantitative real-time PCR
Total RNA was extracted from tissue using TRIzol Reagent. DNA was digested using DNaseI for 30 min at 37 °C. 1 μg of RNA was reverse transcribed into cDNA using random hexamers. Specific primers (for ooG, forward primer, GGCCTACAACTGGAAGATGGC, for oG, forward primer, AGCTTACAACTGGAAGATGGC, and they have the same reverse primer, TAGAAGACACCGCTACTCCT, for GAPDH, forward primer, GGTGAAGGTCGGTGTGAACG, reverse primer, CTCGCTCCTGGAAGATGGTG) were used to quantify by qPCR using the SYBR Green Master Mix on the real-time PCR detection system (Bio-Rad). The copy numbers were determined based on standard curve generated for each gene using known concentration plasmids pSyn-ooG, pSyn-oG, and pcDNA-GAPDH.
Availability of data and materials
The data that used in this study are available from the corresponding author upon reasonable request.
Nassi JJ, et al. Neuroanatomy goes viral! Front Neuroanat. 2015;9:80.
Wickersham IR, et al. Monosynaptic restriction of transsynaptic tracing from single, genetically targeted neurons. Neuron. 2007;53(5):639–47.
Wickersham IR, et al. Retrograde neuronal tracing with a deletion-mutant rabies virus. Nat Methods. 2007;4(1):47–9.
Reardon TR, et al. Rabies virus CVS-N2c(DeltaG) strain enhances retrograde synaptic transfer and neuronal viability. Neuron. 2016;89(4):711–24.
Chatterjee S, et al. Nontoxic, double-deletion-mutant rabies viral vectors for retrograde targeting of projection neurons. Nat Neurosci. 2018;21(4):638–46.
Ciabatti E, et al. Life-long genetic and functional access to neural circuits using self-inactivating rabies virus. Cell. 2017;170(2):382-392 e14.
Callaway EM, Luo L. Monosynaptic circuit tracing with glycoprotein-deleted rabies viruses. J Neurosci. 2015;35(24):8979–85.
Hong S, et al. Functional analysis of various promoters in lentiviral vectors at different stages of in vitro differentiation of mouse embryonic stem cells. Mol Ther. 2007;15(9):1630–9.
Schlabach MR, et al. Synthetic design of strong promoters. Proc Natl Acad Sci U S A. 2010;107(6):2538–43.
Francis DM, Page R. Strategies to optimize protein expression in E. coli. Curr Protoc Protein Sci. 2010;24:1–29 (chapter 5, unit 5).
Zhou Z, et al. Codon usage is an important determinant of gene expression levels largely through its effects on transcription. Proc Natl Acad Sci U S A. 2016;113(41):E6117–25.
Srividya D, Shyam Mohan AH, Rao SN. Expression and purification of codon-optimized cre recombinase in E. coli. Protein Expr Purif. 2020;167:105546.
Jia F, et al. Optimization of the fluorescent protein expression level based on pseudorabies virus bartha strain for neural circuit tracing. Front Neuroanat. 2019;13:63.
Miyamichi K, et al. Dissecting local circuits: parvalbumin interneurons underlie broad feedback control of olfactory bulb output. Neuron. 2013;80(5):1232–45.
Kim EJ, et al. Improved monosynaptic neural circuit tracing using engineered rabies virus glycoproteins. Cell Rep. 2016;15(4):692–9.
Kunec D, Osterrieder N. Codon pair bias is a direct consequence of dinucleotide bias. Cell Rep. 2016;14(1):55–67.
Coleman JR, et al. Virus attenuation by genome-scale changes in codon pair bias. Science. 2008;320(5884):1784–7.
Le Nouen C, et al. Attenuation of human respiratory syncytial virus by genome-scale codon-pair deoptimization. Proc Natl Acad Sci U S A. 2014;111(36):13169–74.
Shen SH, et al. Large-scale recoding of an arbovirus genome to rebalance its insect versus mammalian preference. Proc Natl Acad Sci U S A. 2015;112(15):4749–54.
Ling C, et al. Selective in vivo targeting of human liver tumors by optimized AAV3 vectors in a murine xenograft model. Hum Gene Therapy. 2014;25(12):1023–34.
Paxinos G, Franklin KBJ. The mouse brain in stereotaxic coordinates. 2nd ed. San Diego, Calif. London: Academic. 2001; xxv, 264 p
Qin JY, et al. Systematic comparison of constitutive promoters and the doxycycline-inducible promoter. PLoS ONE. 2010;5(5): e10611.
Norrman K, et al. Quantitative comparison of constitutive promoters in human ES cells. PLoS ONE. 2010;5(8): e12413.
Wang B, et al. Recoding of the vesicular stomatitis virus L gene by computer-aided design provides a live, attenuated vaccine candidate. MBio. 2015. https://doi.org/10.1128/mBio.00237-15.
Gao L, et al. HP-PRRSV is attenuated by de-optimization of codon pair bias in its RNA-dependent RNA polymerase nsp9 gene. Virology. 2015;485:135–44.
Eschke K, et al. Attenuation of a very virulent Marek’s disease herpesvirus (MDV) by codon pair bias deoptimization. PLoS Pathog. 2018;14(1): e1006857.
Meng J, et al. Refining the balance of attenuation and immunogenicity of respiratory syncytial virus by targeted codon deoptimization of virulence genes. MBio. 2014;5(5): e01704-14.
Gaunt E, et al. Elevation of CpG frequencies in influenza A genome attenuates pathogenicity but enhances host response to infection. Elife. 2016;5: e12735.
Burns CC, et al. Genetic inactivation of poliovirus infectivity by increasing the frequencies of CpG and UpA dinucleotides within and across synonymous capsid region codons. J Virol. 2009;83(19):9957–69.
Al-Saif M, Khabar KS. UU/UA dinucleotide frequency reduction in coding regions results in increased mRNA stability and protein expression. Mol Ther. 2012;20(5):954–9.
Tulloch F, et al. RNA virus attenuation by codon pair deoptimisation is an artefact of increases in CpG/UpA dinucleotide frequencies. Elife. 2014;3: e04531.
Bajar BT, et al. Improving brightness and photostability of green and red fluorescent proteins for live cell imaging and FRET reporting. Sci Rep. 2016;6:20889.
Osakada F, et al. New rabies virus variants for monitoring and manipulating activity and gene expression in defined neural circuits. Neuron. 2011;71(4):617–31.
We are grateful to Dr. Xiao Xiao (University of North Carolina at Chapel Hill) for his kind gift of pdsAAV2–CBAp–EGFP plasmids, Dr. Robert Konrad Naumann for the proofreading.
This work was supported by the Key-Area Research and Development Program of Guangdong Province (2018B030331001), the National Natural Science Foundation of China (32071038, 31830035, 91732304, 91632303, 81661148053, and 31771156), the GuangDong Basic and Applied Basic Research Foundation (2021A1515011235), Shenzhen Key Laboratory of Viral Vectors for Biomedicine (ZDSYS20200811142401005), the SIAT Innovation Program for Excellent Young Researchers (E1G023), the National Basic Research Program (973 Program) of China (2015CB755600), the Strategic Priority Research Program (B) (XDB32030200).
Ethics approval and consent to participate
All procedures were approved by the Animal Care and Use Committees at the Shenzhen Institute of Advanced Technology or Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Jia, F., Li, L., Liu, H. et al. Development of a rabies virus-based retrograde tracer with high trans-monosynaptic efficiency by reshuffling glycoprotein. Mol Brain 14, 109 (2021). https://doi.org/10.1186/s13041-021-00821-7