Literature DB >> 24761428

Targeted CNS Delivery Using Human MiniPromoters and Demonstrated Compatibility with Adeno-Associated Viral Vectors.

Charles N de Leeuw1, Frank M Dyka2, Sanford L Boye2, Stéphanie Laprise3, Michelle Zhou3, Alice Y Chou3, Lisa Borretta3, Simone C McInerny3, Kathleen G Banks3, Elodie Portales-Casamar3, Magdalena I Swanson3, Cletus A D'Souza4, Shannon E Boye2, Steven J M Jones5, Robert A Holt6, Daniel Goldowitz1, William W Hauswirth2, Wyeth W Wasserman1, Elizabeth M Simpson7.   

Abstract

Critical for human gene therapy is the availability of small promoter tools to drive gene expression in a highly specific and reproducible manner. We tackled this challenge by developing human DNA MiniPromoters using computational biology and phylogenetic conservation. MiniPromoters were tested in mouse as single-copy knock-ins at the Hprt locus on the X Chromosome, and evaluated for lacZ reporter expression in CNS and non-CNS tissue. Eighteen novel MiniPromoters driving expression in mouse brain were identified, two MiniPromoters for driving pan-neuronal expression, and 17 MiniPromoters for the mouse eye. Key areas of therapeutic interest were represented in this set: the cerebral cortex, embryonic hypothalamus, spinal cord, bipolar and ganglion cells of the retina, and skeletal muscle. We also demonstrated that three retinal ganglion cell MiniPromoters exhibit similar cell-type specificity when delivered via adeno-associated virus (AAV) vectors intravitreally. We conclude that our methodology and characterization has resulted in desirable expression characteristics that are intrinsic to the MiniPromoter, not dictated by copy number effects or genomic location, and results in constructs predisposed to success in AAV. These MiniPromoters are immediately applicable for pre-clinical studies towards gene therapy in humans, and are publicly available to facilitate basic and clinical research, and human gene therapy.

Entities:  

Keywords:  AAV; gene therapy; promoters

Year:  2014        PMID: 24761428      PMCID: PMC3992516          DOI: 10.1038/mtm.2013.5

Source DB:  PubMed          Journal:  Mol Ther Methods Clin Dev        ISSN: 2329-0501            Impact factor:   6.698


Introduction

Several research groups have focused on genome-wide expression analyses in mouse brain.[1-4] However, these projects are limited in their ability to provide information on the location and function of specific regulatory elements that drive the expression pattern. Recently, the VISTA enhancer project has generated data regarding specific putative regulatory elements.[5] Identification of such regulatory elements facilitates the design of compact promoters that could be used in downstream clinical applications. Gene therapies for severe brain and eye disorders hold great therapeutic promise. Adeno-associated virus (AAV) is likely to be a key delivery mechanism due to its nonpathogenic, noninsertional, and low immunogenicity characteristics.[6] However, because of its small size, the DNA payload is severely limited. To accommodate such space restrictions, compact promoters will need to be developed. Many gene therapy studies have employed ubiquitous promoters to drive expression; however, this strategy can be restricted by off-target side effects. To limit such effects, the development of region-specific or cell type–specific promoters will be crucial. In addition, physiological levels of gene expression may be most appropriate for gene-based therapeutics. The use of functional endogenous promoters, which confer physiological levels of expression, may even result in higher expression than ubiquitous promoters in particular cell types.[7] Finally, the presence of viral sequences may enhance transgene inactivation, and off-target expression may increase immunogenicity, resulting in failure of long-term therapeutic results.[8] The Pleiades Promoter Project aims to overcome these biological challenges by generating small promoters (MiniPromoters (MiniPs)) of purely human DNA content, which exhibit highly specific expression patterns, and that could be used in space-constrained viral vectors. We previously published the first set of such Pleiades (Ple) MiniPs, each of which was 4 kilobases or less in length and able to drive regional or cell type–specific expression from the mouse genome in the mouse brain.[9] With the new work presented here, we have undertaken additional promoter development, characterize more deeply some previous designs, and, most importantly, tested a subset of MiniPs in AAV. While our primary target tissue remains adult brain for gene therapy, other central nervous system (CNS) tissues use similar transcriptional programs and are important for gene therapy. Therefore, MiniPromoter characterization has been extended to the spinal cord and retina. Eighteen novel brain MiniPs are described herein, with 2 identified as delivering near pan-neuronal expression in the adult mouse brain, 13 in the spinal cord, and 6 validated for use in developmental studies. We further demonstrate 17 MiniPs with expression in the eye, mostly in the retina, but including a subset directing expression to the cornea or lens. Three MiniPs are demonstrated to retain the capacity to target expression to the ganglion cell layer when delivered to the eye in AAV vectors. These novel tools will significantly improve current methodologies of gene therapy molecular medicine through increased specificity in vector constructs for key anatomical regions of therapeutic interest.

Results

33 MiniPs characterized

All MiniPs contain either a promoter fragment (Prom) or a longer Prom (LongProm) segment usually spanning a known transcription start site from a mammalian endogenous gene. In general, using computational biology and phylogenetic conservation, we designed and constructed four MiniPs for each source gene. While some MiniPs only contained the minimal promoter element, most contained additional conserved putative regulatory elements, or regulatory regions (RRs), in a configuration 5′ of the minimal promoter. All of the sequences used in generating the MiniPs are composed of entirely human sequence generated by polymerase chain reaction (PCR) from the RP11 BAC library or directly synthesized using the human genome reference sequence (hg18, March 2006 reference genome). The 33 MiniP designs studied herein cover 75,777 base pairs (bp) of the human genome (Table 1). Most importantly, we have developed 18 novel MiniPs for expression in the brain. In addition, we have undertaken and expanded characterization of 12 previously described MiniPs from our group.[9] Finally, we characterize three previously positive MiniPs using a more sensitive reporter, lacZ.
Table 1

Summary of MiniPromoters

Ple 
number[a]Gene RRsPlasmid 
pEMS number[b]Strain 
categoryBrain 
expression[c]Eye 
expression
12AVP1482Novel MiniPromoterThis studyNegative
15C8ORF461485Novel MiniPromoterThis studyNegative[d]
16C8ORF461486Novel MiniPromoterThis studyNegative[d]
17C8ORF461487Expanded characterizationPortales-Casamar et al.[9]Negative
22CCKBR1493Novel MiniPromoterThis studyNegative[d]
23CCKBR1494Novel MiniPromoterThis studyNegative[d]
24CCKBR1495Expanded characterizationPortales-Casamar et al.[9]GCL, rare INL
25CCKBR1496Expanded characterizationPortales-Casamar et al.[9]GCL, rare INL
26CCL271497Expanded characterizationPortales-Casamar et al.[9]All layers; cornea; optic nerve
27CCL271498Novel MiniPromoterThis studyNegative[d]
28CCL271499Novel MiniPromoterThis studyAll layers; lens
29CCL271500Novel MiniPromoterThis studyNegative[d]
32CLDN51503Expanded characterizationPortales-Casamar et al.[9]GCL, INL
34CLDN51505Expanded characterizationPortales-Casamar et al.[9]GCL, INL, IPL, PR
53DCX1524Different reporterThis studyGCL, rare INL; optic nerve
55DCX1526Expanded characterizationPortales-Casamar et al.[9]GCL, rare INL; optic nerve
67FEV1538Different reporterThis studyGCL, rare INL
88GFAP1559Expanded characterizationPortales-Casamar et al.[9]Astrocytic layer; optic nerve
97GPX31568Novel MiniPromoterThis studyNegative[d]
112HCRT1583Different reporterThis studyNegative[d]
122ICMT1593Novel MiniPromoterThis studyNegative[d]
123ICMT1594Expanded characterizationPortales-Casamar et al.[9]All layers; cornea; optic nerve
131MKI671602Expanded characterizationPortales-Casamar et al.[9]Cornea
140NR2E11612Expanded characterizationPortales-Casamar et al.[9]Negative
146NTSR11617Novel MiniPromoterThis studyGCL, rare INL
155PCP21626Novel MiniPromoterThis studyINL, IPL
167POGZ1091Expanded characterizationPortales-Casamar et al.[9]All layers; cornea
170POGZ1641Novel MiniPromoterThis studyAll layers; cornea
201SLC6A51673Novel MiniPromoterThis studyNegative[d]
232TNNT11704Novel MiniPromoterThis studyNegative
235TRH1707Novel MiniPromoterThis studyNegative[d]
238TRH1710Novel MiniPromoterThis studyNegative
240UGT81712Novel MiniPromoterThis studyGCL, INL, PR; cornea; optic nerve

GCL, ganglion cell layer; INL, inner nuclear layer; IPL, inner plexiform layer; PR, photoreceptors; RRs, regulatory regions.

Ple MiniPromoters drive lacZ, except Ple167 which drives EGFP/cre.

Plasmid constructs available from AddGene (http://wwww.addgene.org/).

All constructs were positive in the brain.

These constructs were analyzed in chimeras and thus might be false negatives in the eye.

While the primary aim of promoter design was to develop tools for adult brain gene therapy, we recognized the growing importance of gene therapy for the eye and made that an area of focused characterization for MiniPs (Table 1). As such, we generated a unique resource of 17 MiniPs for expression in the mouse eye, including the neural retina, lens, cornea, and optic nerve.

Novel brain MiniP expression ranged from ubiquitous to highly regionalized

Novel MiniPs were analyzed for brain expression (Figure 1), and detailed expression analysis can be found in the Supplementary Data. Positive MiniPs can be grouped by their primary site of expression in the mouse brain. Four MiniPs (Ple27 (CCL27 RRs), Ple28 (CCL27 RRs), Ple122 (ICMT RRs), and Ple170 (POGZ RRs)) expressed in nearly all brain regions. Three MiniPs, Ple15 (C8ORF46 RRs), Ple16 (C8ORF46 RRs), and Ple23 (CCKBR RRs), demonstrated expression primarily limited to the cortex. Ple240 (UGT8 RRs) expression was restricted to the white matter. Several MiniPs expressed strongest in the diencephalon (Ple12 (AVP RRs), Ple29 (CCL27 RRs), Ple97 (GPX3 RRs), Ple112 (HCRT RRs), Ple146 (NTSR1 RRs), Ple155 (PCP2 RRs), and Ple232 (TNNT1 RRs)). Ple67 was expressed in regions outside of the telencephalon. Four MiniPs expressed primarily in the hindbrain—Ple22 (CCKBR RRs), Ple201 (SLC6A5 RRs), Ple235 (TRH RRs), and Ple238 (TRH RRs). This group of MiniPs represents a novel and diverse collection of brain promoters.
Figure 1

Brain expression patterns of 18 novel MiniPromoters and 2 additional MiniPromoter strains that were remade with the lacZ reporter. Adult mice were perfused, brains postfixed, sectioned sagittally or coronally at 1 mm, and stained overnight for lacZ using the substrate X-gal. At least two germline animals were analyzed for each established strain. Ple12 (AVP RRs) expressed in amygdala and thalamus (top panel), and in the ventral hippocampus and cortex layer IV (bottom panel). Ple15 (C8ORF46 RRs) expressed throughout the cortex and throughout the hippocampal pyramidal and dentate granule neurons with strong signal in the dentate gyrus. Ple16 (C8ORF46 RRs) showed a similar but weaker expression pattern to Ple15. Ple22 (CCKBR RRs) had lacZ-positive cells located in the periaqueductal gray (top panel, arrow) and in the hypothalamus (bottom panel, arrow). Ple23 displayed cortical staining throughout the brain, colliculi staining, a strong localized area at the anterior thalamus (top panel, arrow), and in the midbrain colliculi (bottom panel, arrow head), cerebellar granular and Purkinje cells (bottom panel, arrow). Ple27 and Ple28 (CCL27 RRs) showed expression throughout most brain regions, with Ple28 being much stronger. Ple29 (CCL27 RRs) expressed weakly in small localized areas of the thalamus and hypothalamus (bottom panel, arrow and arrowhead, respectively). Ple67 (FEV RRs) expressed in cortical layers II/III, cerebellar Purkinje cells, and strongly throughout most midbrain and hindbrain areas, with a definitive exclusion of striatal tissue (top panel). Raphe nuclei are clearly positive (bottom panel). Ple97 (GPX3 RRs) stained cells localized at the anterior thalamus (top panel, arrow) and ventrolateral hippocampus (bottom panel, arrow), and also in large blood vessels (top panel, arrow head). Ple112 (HCRT RRs) expressed in the anteroventral thalamus (top panel, arrow) and lateral hypothalamic area (top panel, arrow head). The lateral geniculate complex was positive (bottom panel, arrow). Ple122 (ICMT RRs) expressed strongly throughout all major brain regions, with comparatively weaker expression found in the thalamus, and hippocampal CA2, 3, 4 and dentate gyrus (top panel). Ple146 (NTSR1 RRs) had scattered staining in the anterior thalamus (top panel, arrow) and in Purkinje cells (bottom panel, arrow). Ple155 (PCP2 RRs) expressed strongly in the anterodorsal, anteroventral, and anteromedial nuclei (top panel, arrow). Additional staining was seen in the retrosplenial area (bottom panel, arrow). Ple170 (POGZ RRs) expressed strongly throughout all brain regions, with weaker staining in the striatal nuclei. The olfactory bulb mitral, glomerular, and granular layers are all stained (bottom panel). Ple201 (SLC6A5 RRs) expressed weakly in the anterior and retrosplenial areas of the cortex (top panel, arrow). Staining was present in the medulla also (bottom panel, arrows). Ple232 (TNNT1 RRs) expressed in the anterior cingulate and zona incerta (top panel, arrow and arrow head, respectively) and in the lateral hippocampus (bottom panel, arrow). Ple235 (TRH RRs) expressed in the pontine grey and periaqueductal grey (top and bottom panels respectively, arrows). Ple238 (TRH RRs) had scattered expression in many brain regions, including the cerebellar interposed and fastigial nuclei and in parts of the medulla (bottom panel). Ple240 (UGT8 RRs) showed strong fiber-tract associated staining in all parts of the brain, with cortex layer IV strongly positive and in thalamic intralaminar and anterior nuclei groups. In addition, the olfactory bulb glomeruli and corpus callosum were clearly labeled. Bs, brainstem; CA, cornus ammonis; Cb, cerebellum; CH, chimera; Ctx, cortex; Hipp, hippocampus; Hth, hypothalamus; Mb, midbrain; RRs, regulatory regions; Str, striatum; Th, thalamus.

The majority of novel MiniPs showed unique expression patterns. However, we found 33% (i.e., Ple12, Ple22, Ple23, Ple29, Ple155, and Ple232) of the 18 novel MiniPs delivered expression in the zona incerta of the brain, suggestive of a common regulatory property or function in this brain region despite otherwise diverse expression. Previously, we encountered examples of MiniPs negative with the enhanced green fluorescent protein (EGFP) reporter gene, even with antibody-based detection methods. However, subsequent testing with lacZ found them to be positive (unpublished observation and ref. 9). This data implied that the lacZ marker is more sensitive than EGFP for evaluating MiniP expression. Thus, we set out to determine whether the use of lacZ with previous positive designs would broaden the expression pattern. We evaluated the expression of three of previously published MiniPs driving the EGFP reporter[9]: Ple53 (Doublecortin (DCX) RRs; expression in neurons), Ple67 (FEV RRs; expression in serotonergic brain regions), and Ple112 (HCRT RRs; expression in hypothalamus), but now driving the more sensitive lacZ reporter (Figure 1; and DCX results). All three demonstrated an overlapping but expanded expression pattern when compared with their published EGFP counterparts, suggesting that the previous EGFP reporter was expressed below antibody-detection limits in some cells. Variability in expression was observed for three of the eight novel MiniP and new reporter germline strains (Ple12, Ple232, and Ple240; n ≥ 4 for each strain). However, only one (Ple232) demonstrated extreme variability from nearly negative in the brain to strongly positive. Thus, we conclude that the in vivo expression variability of MiniPs docked at Hprt is generally low.

Ple26 and Ple170 MiniPs drive near pan-neuronal expression

Establishing a pan-neuronal promoter for the adult brain has proven to be challenging in the field due to the highly diverse collection of neuronal subtypes represented with numerous nonoverlapping transcriptional programs.[10] Recently, Heimer-McGinn and Young[11] published a random-insertion transgenic mouse driving cre recombinase in an apparently pan-neuronal manner. However, although a great resource for lineage tracing studies, there was no promoter fragment that can be used to create other mouse strains, or be used in viral or plasmid delivery systems. We identified two MiniPs, Ple26 (CCL27 RRs), originally described in Portales-Casamar et al.,[9] and Ple170 (POGZ RRs; a novel MiniP), as pan-neuronal candidates. In this study, we expanded our characterization and assessed the extent to which these MiniPs label neurons. These two MiniPs were assessed for neuronal costaining using NeuN and off-target glial expression using Gfap (Figure 2).
Figure 2

Two MiniPromoters, Ple26 (CCL27 RRs) and Ple170 (POGZ RRs), drive near pan-neuronal expression throughout the mouse brain. For each, a representative 1-mm sagittal and coronal whole brain section is shown, as well as cryosections of NeuN (neuronal marker) and Gfap (glial marker) antibody co-labeling with X-gal in the cortex, hippocampus, and cerebellum. (a) Ple26 (CCL27 RRs) showed marked expression throughout the adult brain with slightly less labeling in the olfactory bulb. In the cortex, NeuN closely co-labeled many of the X-gal positive cells (black arrows). However, there were some NeuN cells that do not co-label with X-gal (red arrows), as well as many smaller punctate X-gal nuclei that do not co-label with NeuN (white arrows). Nearly all hippocampal neurons co-label with X-gal, and extensive co-labeling was also observed in the cerebellum, with X-gal only staining also present in the molecular layer. Also, rare X-gal/Gfap-positive cells were observed (black arrows). (b) Ple170 (POGZ RRs) demonstrated a similar expression pattern to that of Ple26. The CA2 and CA3 regions of the hippocampus showed less staining than the dentate gyrus and CA1 regions. Many NeuN co-labeled cells were visible (black arrows). However, some NeuN-positive cells in the cortex were not labeled with X-gal (red arrows), as well as some punctate X-gal nuclei were NeuN-negative (white arrows). Reduced NeuN co-labeling was seen in the hippocampus and cerebellum, particularly in the cerebellar granule cell layer. Occasionally, a cell was found to be double labeled with Gfap and X-gal (black arrows). CA, cornus ammonis; Cb, cerebellum; Ctx, cortex; Hipp, hippocampus; RRs, regulatory regions. (Scale bars = 100 μm).

Both MiniPs were positive throughout the cortex, the hippocampus, and the cerebellum (Figure 2a,b). The olfactory bulb was more weakly stained with Ple26 as compared with the rest of the brain and with Ple170. The general overlap with NeuN was more pronounced for Ple26 than for Ple170. This was particularly true for the cerebellum where comparatively weak labeling was observed for Ple170. Using hippocampal neurons and estimating the average coverage based on co-labeling in cornus ammonis 2–3 regions, we found a significant difference of 91.34% co-labeling with NeuN for Ple26 and 82.89% for Ple170 (n = 3 animals for each, and 3 sections averaged for each animal; t-test, P < 0.05). This data suggests that Ple26 would be superior to Ple170 for use in driving near pan-neuronal expression.

MiniP expression observed in spinal cord and non-CNS tissue

For 21 MiniPs, we collected adult spinal cord, heart, liver, lung, and ear-notch tissue samples to be stained for β-galactosidase (β-gal) (see Supplementary Table S1). Within this analysis, we looked at all germline strains except for Ple25 (CCKBR RRs), which was analyzed prior to the introduction of this protocol. Thirteen MiniPs stained positive in spinal cord (Figure 3), 12 in heart, 6 in liver, 8 in lung, and 14 in ear-notch samples. Nearly half (46%) of the 13 CNS-positive MiniPs were also expressed in all four additional tissues. Interestingly, only four (31%) were restricted to spinal cord, and hence specific to the CNS.
Figure 3

Thirteen MiniPromoters expressed in the mouse spinal cord. LacZ expression was detected by X-gal immunohistochemistry in whole mount thoracic spinal cord. (a) Ple17 (C8ORF46 RRs), Ple28 (CCL27 RRs), Ple32 (CLDN5 RRs), and Ple170 (POGZ RRs) showed scattered expression in the central spinal cord and in nerve bundles. Ple24 (CCKBR RRs) expression was seen primarily in the spinal cord with rare puncta in the dorsal root ganglia (data not shown). Ple26 (CCL27 RRs), Ple53 (DCX RRs), and Ple55 (DCX RRs) expressed throughout the gray matter regions and also in dorsal root ganglia; the latter particularly strong for Ple53 and Ple55. Ple34 (CLDN5 RRs) and Ple123 (ICMT RRs) expressed throughout all of the spinal cord, while Ple123 demonstrated weaker expression in nerve bundles but strong expression in dorsal root ganglia. Ple88 (GFAP RRs) was expressed only in patches along white matter regions of the central spinal cord. Ple167 (POGZ RRs) expressed in patches in the spinal cord with some nerve fibers and occasional dorsal root ganglia puncta also positive (data not shown). Ple240 (UGT8 RRs) is expressed strongly and predominantly along spinal nerve cord bundles with weaker expression in the central spinal cord tissue and in dorsal root ganglia. (b) Summary of expression pattern for each construct in the three key regions of the spinal cord. CH, chimera; DRG, dorsal root ganglion; ND, not done; RRs, regulatory regions; √, expression present; ×, expression absent.

MiniPs that were positive in the spinal cord were often expressed in the nerve fibers exiting the spinal column and/or dorsal root ganglia. Overview images for spinal cord expression are shown in Figure 3a and summarized in Figure 3b. We observed expression patterns that consisted of only the central spinal cord (n = 1, Ple88 (GFAP RRs)), spinal cord with nerve fibers (n = 4; Ple17 (C8ORF46 RRs), Ple28 (CCL27 RRs), Ple32 (CLDN5 RRs), and Ple170 (POGZ RRs)), spinal cord with dorsal root ganglia (n = 4; Ple24 (CCKBR RRs), Ple26 (CCL27 RRs), Ple53 (DCX RRs), and Ple55 (DCX RRs)), and spinal cord with nerve fibers and dorsal root ganglia positive (n = 4; Ple34 (CLDN5 RRs), Ple123 (ICMT RRs), Ple167 (POGZ RRs), and Ple240 (UGT8 RRs)).

Embryonic analyses of six MiniPs revealed hypothalamus-specific and skeletal muscle expressing constructs

Six MiniPs were analyzed at embryonic day (E) 12.5 and assessed by X-gal staining to gain insight into their developmental expression (Figure 4). For each MiniP, at least two hemizygous males and one litter-matched wild-type were stained for X-gal. We selected three MiniPs with widespread adult brain expression data: Ple26 (CCL27 RRs), Ple123 (ICMT RRs), and Ple170 (POGZ RRs). We selected three further MiniPs of specific interest: Ple131 (MKI67 RRs) for the role of KI-67 in proliferation,[12] Ple140 (NR2E1 RRs) because of its highly regionalized hypothalamic expression in adults,[9] and Ple232 (TNNT1 RRs) for the role of TNNT1 in skeletal muscle.[13]
Figure 4

Expression of six MiniPromoters during development. E12.5 embryos were harvested and stained for lacZ expression using X-gal. Whole-mount (first panels) and partially cleared whole-mount embryos (second panels) were photographed. (a) Ple26 (CCL27 RRs) stained throughout the entire embryo. (b) Ple123 (ICMT RRs) displayed variability, with the most common expression pattern being limited to the developing heart as shown. (c) Ple131 (MKI67 RRs) was expressed throughout the embryo at varying levels, most prominently in the heart. (d) Ple140 (NR2E1 RRs) showed a very strongly defined pattern in the hypothalamus of the embryo. (e) Ple170 (POGZ RRs) demonstrated staining in all parts of the embryo, but generally stronger outside of the nervous system. (f) Ple232 (TNNT1 RRs) stained strongly in skeletal muscle and the tongue of the embryo but was not visible in the brain or spinal cord. RRs, regulatory regions.

Ple26 (CCL27 RRs) (Figure 4a), a pan-neuronal candidate, directed strong expression throughout the embryo, indicative of ubiquitous expression at this developmental stage. Ple123 (ICMT RRs) and Ple131 (MKI67 RRs) were both predominantly expressed in the heart (Figure 4b,c), with Ple131 also marking scattered cells throughout all other regions of the embryo. Our previously published highly specific adult hypothalamus promoter, Ple140 (NR2E1 RRs) (Figure 4d), demonstrated exclusive expression to the developing hypothalamus, consistent with the observation in the adult.[9] This early developmental specificity suggests expression in a cell type that maintains its identity throughout development. Interestingly, Ple170 (POGZ RRs) (Figure 4e), the second pan-neuronal candidate, is expressed more strongly throughout the limbs and tail, with a blue-green haze of expression throughout the rest of the embryo, including the nervous system. The Ple232 (TNNT1 RRs) (Figure 4f) MiniP showed strong expression in the developing skeletal muscles of the embryo. At E12.5, 50% of MiniPs (Ple123, Ple131, and Ple170) showed variability in expression, in contrast with the relative lack of variability in the adult. Figure 4 presents the most frequently observed staining pattern. For Ple123 (n = 7 embryos), expression levels ranged from primarily heart only (n = 5) to ubiquitous (n = 2). Ple131 (n = 5 embryos) displayed heart only (n = 1) to weak ubiquitous (n = 3) to strong ubiquitous (n = 1) expression. Ple170 (n = 3 embryos) ranged from negative (n = 2) to ubiquitous (n = 1) expression. Interestingly, in contrast to the adult brain, we did not observe variability in the embryo for Ple232 expression.

DCX-based MiniPs demonstrate expression driven by a “Prom” and a candidate regulatory element

We compared two MiniPs, both designed using the DCX gene, to understand differences in the RRs used (Figure 5a,b). The first design, Ple53, was based on a known promoter (LongProm) used in the literature[14] and found to be expressed as anticipated in neurogenic regions of adult brain. The second design, Ple55, featured a 5′ shortened promoter region (Prom) coupled to a highly conserved putative regulatory element from within intron 4.
Figure 5

Comparison of two DCX-based MiniPromoter constructs with similar but nonidentical expression patterns identifies putative functional TFBS. (a) Diagram showing the differences between Ple53 and Ple55. Both MiniPromoters contained a minimal “Prom” sequence of 2,353 bp; however, Ple53 contained an additional but contiguous sequence and thus was named “LongProm”, whereas Ple55 contained a noncontiguous segment 3 in addition to the Prom. (b) In the UCSC genome browser, predicted TFBS for Mafk and BRACH protein families were observed in the extended LongProm sequence specific to Ple53. Predicted RXR/RAR heterodimer and SRF sites were found in the segment 3 specific to Ple55. Comparison of the expression of lacZ using the substrate X-gal (blue) demonstrated that (c, e, g) Ple53-lacZ was expressed as a subset of the more expansive (d, f, h) Ple55-lacZ expression pattern. (c, d) The comparison demonstrated that with the Ple53 construct, the majority of cells labeled in the cortex were located in the deeper cortical layers, whereas with Ple55, they were distributed throughout the cortex. There were also more cells in Ple55 that did not co-label with NeuN. (e, f) Expression analysis in E12.5 X-gal stained whole-mount and partially cleared embryos for each MiniPromoter demonstrated that Ple55 contained all of the neuronal based expression observed in Ple53 plus staining in the precartilage or premuscle of the face and limb buds. Both MiniPromoters exhibited expression at the distal edge of the limb buds, and both included staining of the developing kidneys (white arrows). A close-up of the face and front limb showed the additional staining in Ple55 (black arrows). (g, h) (upper panels) MiniPromoter lacZ expression as detected by X-gal reaction in the retina (left panel). Expression for the DCX-based designs was limited to the ganglion cell layer (GCL) for both Ple53 and Ple55, with rare cells observed in the inner nuclear layer. NeuN (marking all retinal ganglion and amacrine cells) co-labeled with X-gal (center panel), but Gfap (retinal astrocytes) did not (right panel). Rare cells did not co-label with NeuN (white arrows). (g, h) (lower panels) Immunofluorescence co-labeling with Brn3, βtubIII, and Sox2. Co-labeling demonstrated that most lacZ-positive cells stained with βtubIII and Brn3 markers (ganglion cells) for both MiniPromoters (left and center panels), and only rarely with Sox2-positive cholinergic amacrine cells (right panel, white arrows). Ctx, cortex; Hipp, hippocampus; NR, neutral red; RAR, retinoic acid receptor; RRs, regulatory regions; RXR, retinoid X receptor; SRF, serum response factor; TFBS, transcription factor binding sites; UCSC, University of California, Santa Cruz. (Scale bars = 100 μm, except for (g, bottom panel) fluorescence images 20 μm).

We first looked at the cortical layer specificity of Ple53 and Ple55 (Figure 5c,d). Ple53 expression, as detected by X-gal and counterstained with neutral red, was limited to the deepest cortical layer VI (Figure 5c, left panel). Rare positive cells in layers II, IV, and V were observed. In contrast, Ple55 was expressed in all cortical layers, including layer I (Figure 5d, left panel). Nearly all layer II, IV, and VI neurons were labeled. Both MiniPs were expressed in the known adult neurogenic regions of the dentate subgranular zone, subventricular zone, and rostral migratory stream axis (data not shown). To assess whether the positive cells in the cortex were primarily neurons, we compared the extent of NeuN and Gfap co-labeling with X-gal in the cortex of these two strains (Figure 5c,d, center and right panels). Cells were found in both strains that co-labeled (black arrows) and did not co-label (white arrows) with NeuN (Figure 5c,d). We observed some co-labeling with Gfap in the hippocampal formation, which was more extensive with Ple55 (Figure 5d, right panel). The overlapping cortical and neurogenic expression between the two designs suggests an expansion of the expression pattern with Ple55. At E12.5 (Figure 5e,f), we noted an overall conserved staining pattern for Ple53 and Ple55, consistent with the developing nervous system. Specifically, expression was observed throughout the developing brain, eye, spinal cord, and dorsal root ganglia. Labeling extends along both the dorsal and ventral ramus. Expression in the developing kidney was found in partially cleared embryos (Figure 5e,f, white arrows). Ple55 showed additional expression in the face and limbs (black arrows). Only the nervous system expression and kidney expression have been previously documented for the mouse Dcx gene.[15] In the retina (Figure 5g,h), Ple53 and Ple55 were expressed nearly exclusively in the ganglion cell layer (GCL), with rare inner nuclear layer (INL) cells positive. Both MiniPs were also expressed in the optic nerve, axonal extensions from the GCL. Studies by others in rat retina suggested horizontal cell labeling for endogenous Dcx.[16] However, only the rare INL MiniP-expressing cells (white arrows in Figure 5g,h, top middle panel) may be consistent with their localized expression in horizontal cells. Thus, we conclude that the retinal expression for Ple53 and Ple55 is ectopic and unexpected. To further the characterization of the retinal expression, we performed co-labeling studies. Strong overlap between NeuN-positive nuclei and X-gal (Figure 5g,h, top middle panel), and with rare exceptions, a lack of overlap with the astroglial marker Gfap (Figure 5g,h, top right panel), was observed. Further analysis by immunofluorescence of β-gal using co-labeling with β-Tubulin III (Figure 5g,h, bottom left panel) and Brn3 (Figure 5g,h, bottom middle panel), both markers of retinal ganglion cells, demonstrated significant overlap, while only rare cells were positive for the cholinergic amacrine marker, Sox2 (Figure 5g,h, bottom right panel, white arrows). To account for the expression differences observed between Ple53 and Ple55, we undertook a bioinformatics analysis (Figure 5a,b). The 5′ extension of the Prom segment found in the LongProm used in Ple53 and the additional element from intron 4 added together with the Prom segment in Ple55 were used to search for predicted transcription factor binding sites (TFBS). Ple53 was predicted to contain a MafK (a transcription factor binding site based on the MafK gene family) TFBS (UCSC (University of California, Santa Cruz)/ChIP (chromatin immunoprecipitation)-seq ENCODE data) and a BRACH (a transcription factor binding site based on the brachyury gene family) TFBS (TFBIND database), which are located in the 5′ nonoverlapping region of the LongProm and thus missing in Ple55. Furthermore, Ple55 contained two predicted serum response factor (analyzed using CONSITE with profiles from the JASPAR database) and two retinoid X receptor/retinoic acid receptor (JASPAR database) TFBSs within the conserved element 3, derived from the intron 4 region. Serum response factor is a transcription factor critical for muscle development,[17] and a TFBS for it is located in element 3 of Ple55. We propose that this serum response factor TFBS in Ple55 is responsible for the additional expression pattern in embryos and that the tissue identity is premuscle. Lastly, it appeared that there was greater co-labeling with Gfap in the hippocampal formation of adult Ple55 mouse brain, which would be consistent with a role for serum response factor in astrocyte and oligodendrogial development.[18] This unique predicted TFBS may account for the expression differences observed. Interestingly, in contrast to all other MiniP strains, we were unable to maintain the DCX-based MiniP strains on the C57BL/6J background past backcross N3 due to lack of allele transmission. However, we were able to restore and maintain normal allele transmission after crossing the allele to the 129S1/SvImJ strain and hypothesized an unknown genotoxic effect for these alleles in C57BL/6J mice.

“Proms” typically confer the observed specificity

It has been proposed that the minimal promoter for a given gene confers the general specificity of expression but that additional regulatory elements either alter the timing or fine-tune the regional specificity.[19,20] To test this concept, we analyzed sets of MiniPs designed from the same gene in our combined dataset (this paper and ref. 9). For the purpose of this analysis, we classified MiniPs (with respect to brain only) as: (i) related, having at least some expression similar to the endogenous gene; or (ii) unrelated, having no pattern similar to the endogenous gene. Similarly, we classified sets of MiniPs as: (i) related (i.e., all MiniPs are related to endogenous and similar in pattern to each other); (ii) unrelated (i.e., all MiniPs do not have expression in common with the endogenous gene but are similar in pattern to each other); or (iii) mixed (i.e., a set containing related and unrelated MiniPs, or alternatively, a set containing MiniPs that are all related to the endogenous expression pattern but not similar to each other). We had 16 such sets (see Supplementary Figure S1). Nine of these sets contain MiniPs that all have an identical basal Prom (or minimal promoter) element within the set. Of these nine, six sets of MiniPs (C8ORF46, CLDN5, HCRT, ICMT, MKI67, and S100β RRs) displayed related expression, while the remaining two (POGZ and RGS16 RRs) displayed unrelated expression, and one set had one related and one unrelated MiniP (NR2E1 RRs). An additional six sets included a combination of MiniPs that either contained the Prom element or a LongProm element (usually a known promoter in the literature and contiguous but larger than the Prom element) within the given set (i.e., overlapping basal promoters). Five of these six sets (CCKBR, DBH, DCX, GFAP, and PITX3 RRs) had 5′ extensions compared with the Prom, and one set (TRH RRs) contained a 3′ extension compared with the Prom. All six sets displayed related expression. Lastly, the CCL27-based set contained MiniPs with two nonoverlapping Prom1 and Prom2 designs. Prom1 was located further upstream and annotated as a putative cortical expression promoter based on 5′ cap analysis of gene expression data. A Prom2-based MiniP, located near the known transcription start site, was still related but distinctly different from those based on Prom1. Due to this difference, we classified this as a mixed set. Seventy-five percent (12/16) of the MiniP sets displayed related expression in all positive MiniPs. Thus, we conclude that the common Prom element is generally sufficient to confer specificity, as compared with the endogenous mouse gene, and added putative regulatory elements either enhance or reduce expression levels, or alternatively expand or narrow the expression pattern.

Seventeen MiniPs deliver expression to the adult mouse eye

The gene set for MiniP design was originally chosen for region-specific and cell type–specific adult brain expression; however, since the brain and eye are both generated from diencephalon, it is likely that many designs would also be positive in the eye. Thus, we undertook a detailed investigation of eyes in otherwise brain-positive MiniPs. We analyzed the expression pattern of 17 MiniPs found to be positive in the adult mouse eye (Table 1). Within this set, we observed expression in nearly all major retinal cell types. In extra-retinal tissue, we identified seven MiniPs positive in the cornea, six in the optic nerve, and two in the lens. Fifteen eye positive MiniPs are presented in Figure 6, and two additional ones are presented in Figure 5, as part of the DCX MiniP comparison.
Figure 6

Fifteen MiniPromoters expressed in the mouse eye. Mice were perfused, eyes harvested, and stained for lacZ using the substrate X-gal (blue). Eyes were cryoprotected, embedded, and sectioned. Ple24 (CCKBR RRs), Ple25 (CCKBR RRs), Ple32 (CLDN5 RRs), Ple67 (FEV RRs), and Ple146 (NTSR1 RRs) expressed in a subset of GCL cells. Rarely, an amacrine cell in the innermost part of the INL was positive in Ple24, or a rare positive INL cell in Ple32. Ple26 and Ple28 (both based on CCL27 RRs) expressed in the INL, and in photoreceptors as evidenced by the columnar staining of the ONL and photoreceptor segments. Ple26 showed extensive IPL and OPL fiber staining and Ple28 was also positive in the GCL. Ple34 (CLDN5 RRs) and Ple155 (PCP2 RRs) have staining characteristic of bipolar cells, with cell bodies in the INL and processes extending to the GCL (Ple34) or inner half of the IPL (Ple155). Ple88 (GFAP RRs) stained in the retinal astrocytic layer. Ple123 (ICMT RRs) expressed throughout the retina. Ple131 (MKI67 RRs) expressed in cells of the corneal epithelium. Ple167 and Ple170 (both based on POGZ RRs) demonstrated similar staining patterns in the GCL, INL, OPL and photoreceptor segments. Ple240 (UGT8 RRs) had fiber-like staining spanning the IPL and clear cell body labeling in a subset of INL cells. CH, chimera; INL, inner nuclear layer; IPL, inner plexiform layer; ONL, outer nuclear layer; OPL, outer plexiform layer; RRs, regulatory regions. (Scale bar = 100 μm).

We found Ple24 and Ple25 (CCKBR RRs) were expressed in the GCL and occasionally in the inner part of the INL. Ple26 (CCL27 RRs) was expressed in the cornea, optic nerve, and all retinal layers. Another design for this gene, Ple28, showed expression in all layers and also in the lens but no expression in the cornea or optic nerve. Ple32 (CLDN5 RRs) was expressed in the GCL and INL. Ple34 (CLDN5 RRs) was positive in the INL and GCL and was expressed rarely in cells of the outer nuclear layer (photoreceptors). The morphology of cell nuclei located in the INL, with processes extending to the GCL, was suggestive of bipolar ON cells. Ple67 (FEV RRs) stained cells in the GCL and rare INL cells. Ple88, based on a known promoter for GFAP,[21] was expressed in the astrocytic layer as expected. Ple123 (ICMT RRs) stained all retinal layers including the cornea and optic nerve. Ple131 (MKI67 RRs), based on proliferative maker KI-67, stained patches of the cornea only, consistent with other studies showing expression of KI-67 in basal epithelial cells in the cornea.[22] Ple146 (NTSR1 RRs) stained rare GCL and INL cells (INL-positive cell not pictured). Ple155 (PCP2 RRs) stained predominantly in the INL and inner plexiform layer and was strongly suggestive of bipolar ON cells, which is expected expression for this gene. Ple167 and Ple 170 (POGZ RRs) stained cells in all layers of the neural retina including the cornea. Ple240 (UGT8 RRs) stained positive in the retinal inner segment (photoreceptors) and occasionally in the INL with fibers extending into the astrocytic layer, suggesting a Müller glial or bipolar cell type. The cornea and optic nerve were also positive for Ple240. This characterization has led to a set that is unique compared with currently available eye promoters (see ref. 23 for a brief list).

Knock-in MiniPs that express in retinal ganglion cells predict expression of MiniP-GFP using AAV

To assess the feasibility of using the MiniPs for AAV-mediated delivery, we cloned three (Ple25, Ple53, and Ple67) MiniPs that were expressed primarily in the retinal ganglion cell layer, into the AAV2-(quad Y-F) vector. This vector, containing four tyrosine to phenylalanine mutations, has previously been shown to transduce cells of the retina uniformly, with a transduction profile that includes photoreceptors, Müller glia, and horizontal, bipolar, amacrine, and ganglion cells.[24] Two (Ple25 and Ple53) of the three MiniPs were 5′ trimmed by ~200–400 bp to fit into the viral vector due to limits on the carrying capacity of AAV constructs. Ple25 was trimmed from 3,740 bp (KI allele) to 3,312 bp (AAV), and Ple53 was trimmed from 3,520 bp (KI allele) to 3,310 bp (AAV), resulting in Ple25(Fmd1) and Ple53(Fmd1), respectively. Ple67 (2,202 bp) did not require trimming. Mice were injected intravitreally at postnatal day (P) 28, harvested 4 weeks later, and analyzed for green fluorescent protein (GFP) expression. The transduction profile of AAV2(quad Y-F) using the ubiquitous chicken beta-actin promoter demonstrated expression in all retinal layers (Figure 7a). Figure 7b,d,f show a side-by-side comparison of the knock-in mouse strain and the AAV vector–injected mice. All three AAV vectors supported expression of the reporter gene in the retinal ganglion cell layer. For Ple53 and Ple67, we also observed rare positive cells in the innermost aspect of the INL. Importantly, the 5′ trimming of ~200–400 bp for Ple25 and Ple53 to accommodate AAV space constraints did not alter their specificity, indicating that the deleted segments were nonessential for specificity. We conclude that for at least these three MiniPs, the expression pattern in the knock-in strain predicted that obtained using a virus-based delivery to the mouse retina. These data suggest that MiniPromoter specificity may be maintained in AAV vectors.
Figure 7

Three MiniPromoters that expressed in the retinal ganglion cell layer in knock-in mice maintained this restricted expression from AAV viral vectors. (a) Transduction profile of the AAV2(quad Y-F) vector with the CBA promoter driving GFP. Three MiniPromoter constructs (b, c) Ple25, (d, e) Ple53, and (f, g) Ple67, which expressed in the GCL when docked 5′ of Hprt driving β-galactosidase, were recloned into AAV2(quad Y-F) driving hGFP. (b, d, f) Knock-in mouse retinas stained for lacZ using the substrate X-gal (blue) and counterstained with neutral red (first panel) showed similar ganglion expression to intravitreal injected AAV2(quad) with the MiniPromoter driving hGFP (second panel, epifluorescence). (c, e, g) Virus-injected retinas were immunostained for Brn3a (a marker of ganglion cells). First panel shows Brn3a (red) with DAPI (blue), second panel hGFP (green) with DAPI (blue), and third panel the merge (yellow) of Brn3a (red), hGFP (green), with DAPI (gray). CBA, chicken beta-actin promoter; DAPI, 4′,6-diamidino-2-phenylindole; GCL, ganglion cell layer; GFP, green fluorescent protein; INL, inner nuclear layer; ONL, outer nuclear layer; RRs, regulatory regions. (Scale bars = 100 μm, except for fluorescence images 20 μm).

The retinal ganglion cell (RGC) layer in mice contains ~40% ganglion cells and 60% amacrine cells.[25] Nadal-Nicolás et al.[26] estimated that Brn3a+ cells label ~92.2% of RGCs, as co-labeled with the Fluorogold tracer in rat retina. Similarly, Fluorogold was estimated to label 97.8–98.4% of RGCs in rats.[27] Thus, we considered that ~90% of RGCs in the adult mouse retina should be labeled with Brn3a and could be used as a reliable marker of retinal ganglion cells. We undertook Brn3a co-labeling in virus-injected mice for Ple25, Ple53, and Ple67 to assess the extent of expression in ganglion cells, with the remainder of GFP-positive cells in the GCL assumed to be amacrine. As shown in Figure 7c,e,g, the MiniPs ranked qualitatively from most to least overlap with Brn3a: Ple53 > Ple25 > Ple67. We were also interested in how much expression was observed in amacrine (Brn3a-negative) cells in the GCL. Similarly, qualitative ranking of the same MiniPs from most to least expression in Brn3a-negative but GFP-positive cells ordered the MiniPs: Ple67 > Ple25 > Ple53. Thus, we conclude that Ple53 supports expression in the greatest proportion of retinal ganglion cells and is most specific to this cell type, while Ple67 supports expression in a subset of retinal ganglion cells in addition to numerous amacrine cells.

Discussion

Interest in the identification of small promoters for capturing the expression pattern of biologically unique genes is longstanding. The Pleiades approach, with methodology described in detail previously,[9] has the ability to generate novel promoters in a high-throughput manner. Here, we present 18 additional promoters identified by this technique, and the further characterization of 15 previous designs. Our approach focuses on physiologically relevant expression. The reasoning for this is threefold: (i) overexpression of therapeutic molecules can be deleterious[28]; (ii) overexpression of purely research-related molecules such as cre can cause cytotoxicity[19,29,30]; and (iii) introduction of artificial sequences to enhance expression can pose regulatory challenges due to their unknown specificity. Our strategy of mimicking endogenous expression levels using endogenous sequences is aimed at overcoming these challenges and thus enabling improved gene therapy designs. However, there are some limitations to the Pleiades design pipeline. First, a small number of MiniPs displayed highly variable expression when docked at Hprt. We were unable to find previous reports of variable expression for constructs knocked-in at Hprt, which is generally a neutral docking location with consistent expression of docked constructs.[31] Expression variability was also not a general feature of our previous MiniP work.[9] Thus, we hypothesize that the variability observed is construct and sequence specific and reflects the large number of MiniPs tested. For the most part, the extent of variability we saw in embryos correlated with that observed for adults. It remains to be seen whether such construct-dependent variability is genome context specific or would be similarly reproduced in other systems (e.g., virus-based delivery or knock-in at the ROSA26 locus). Additionally, the Pleiades data was limited by the inability to attain germline animals for a small subset of MiniPs. In these cases, chimeras were studied. Although each chimera may not display the entire cellular expression profile, the observed chimeric expression has been a reliable predictor of germline expression.[9] Lastly, we were limited in characterization of the specific cell types of expression, particularly for the brain where hundreds or thousands of different neuronal cell types are predicted to exist.[32] However, this data on regionalized expression of MiniP constructs remains a valuable scientific and therapeutic resource. Importantly, this data represents numerous novel and biologically significant findings: first, we identified the zona incerta as a common region of expression for MiniP designs. We hypothesize that overrepresentation of the zona incerta brain region in this set is due to the diverse range of neurochemical cells located in this region, which includes neurons varying in soma size and shape (see ref. 33 for a comprehensive review on the zona incerta). In addition, extensive networking exists between cells of the zona incerta and numerous other brain regions (e.g., the cerebral cortex, diencephalon, basal ganglia, brainstem, and spinal cord),[33] further suggesting a heterogeneous population in this region. The persistent finding of staining in this region suggests caution for therapeutic delivery where expression in the zona incerta could cause harm. The brain-positive MiniPs provide novel tools for researchers and clinicians alike. Within our toolbox of MiniPs, we have now added constructs that would drive expression in the cortex and hippocampus (Ple15 and Ple16), in oligodendroglia (Ple240), in a near pan-neuronal expression (Ple26 and Ple170), or in the developing nervous system (Ple53 and Ple55). Selected designs were also validated for their use in developmental studies. We tested six adult brain–positive MiniPs and showed expression in the brain for five of these at E12.5. The remaining MiniP, Ple232, was expressed only outside of the embryonic brain, in contrast to its expression in the adult. This finding is not entirely unexpected, since not all genes will be expressed in the brain at the time point used in this study (E12.5, mid-point of neurogenesis). Remarkably, we showed highly specific expression of Ple140 (NR2E1 RRs) in the hypothalamus at E12.5, with no evidence of its expression in other embryonic tissue. These embryonic studies allow us to rapidly detect the functioning of MiniPs outside of the CNS, like Ple232 (TNNT1 RRs) which could prove a valuable tool for researchers studying, and developing treatments for, skeletal muscle disorders. Overall, this embryonic dataset suggests that many other Pleiades MiniPs are likely to be positive during development. We were surprised by the portion of MiniPs expressed in the mouse eye. This represents a novel toolset for driving retinal gene expression, and these MiniPs offer great therapeutic potential. Bipolar cells are the next step in retinal processing after photoreceptors are activated by light. In diseases where photoreceptors are nonfunctional or have undergone cell death, one strategy is to target the bipolar cells.[34] Similarly, since retinal ganglion cells form the last step in the visual processing pathway, expression of therapeutic molecules could be a viable option here also,[34] with perhaps more generic retinal disease application. Our efforts have developed at least two MiniPs with bipolar cell expression (Ple34 and Ple155) and six with ganglion cell layer expression (Ple24, Ple25, Ple53, Ple55, Ple67, and Ple146). The finding of such a large proportion (35.3%) of eye positives specific to the GCL implies that selection of candidates for MiniP design using brain expression data may be biased toward generating RGC-positive designs. We hypothesize this may be due to the similarities between RGCs and other CNS neurons with regards to (i) the ability to fire action potentials and (ii) similarities in gene expression between RGC-specific transcripts and other CNS neurons that fire action potentials.[35-39] Interestingly, the ganglion cell-marking Brn3 transcription factors are also highly expressed in the sensory neurons of the dorsal root ganglia.[40] Of the six GCL-positive MiniPs, we were able to assess dorsal root ganglion expression in three and found all three positive (Ple24, Ple53, and Ple55). Thus, our selection strategy has not only preferentially generated MiniPs positive in the GCL of the retina but also in the dorsal root ganglion (eight out of ten MiniPs tested). Other MiniPs that expressed in the eye also demonstrate potential utility: Ple88 (GFAP RRs) could be used to aid in the treatment of retinal disorders with increased astrogliosis (defined as increased GFAP expression), and Ple131 (MKI67 RRs) may be useful for corneal repair. Curiously, we found occasional discordance between cell type expression in the brain and eye. The Ple155 (PCP2 RRs) MiniP was expressed as expected in the mouse retina (bipolar ON cells), but with seemingly unrelated expression in the brain (MiniP expression present in cortex, thalamus, and hypothalamus but absent in Purkinje cells for which PCP/L7 is well known).[41,42] This data further supports the notion that cell type specificity is regulatory element driven in specific organs or tissues.[43] Altogether, this represents a powerful set of tools for driving gene expression in the eye. One concern was whether the observed expression pattern of MiniPs would be reliably reproduced when used in viral vectors. Viral vector gene delivery is the preferred therapeutic method for the brain and eye, primarily because the number of treatment options via gene delivery is vast. We showed that a subset of three MiniPs (Ple25, Ple53, and Ple67), which expressed in the retinal ganglion cell layer in knock-in mice, exhibited a similar expression pattern in the mouse retina when delivered intravitreally and expressing from recombinant AAV2. We hypothesize that the characterization by single-copy site-specific knock-ins, and use of exclusively human DNA, may have predisposed these MiniPs to success in AAV. Using this methodology, the desirable expression characteristics are more likely to be inherent to the MiniP only and not dictated by copy-number or the genomic location. In addition, human DNA should be less prone to expression variability compared with more distantly related sequence or viral promoters. Finally, critical to the success of this work was the choice of an AAV serotype with a broad transduction profile in the retina, which included the target cell types of the MiniPs. Interestingly, with this success, the MiniPs can also be used to map the specificity of novel serotypes, providing a valuable tool for researchers working on AAV capsid biology. An interesting pattern emerged using these three MiniPs: we observed an inverse qualitative relationship between the overlap of stained retinal ganglion cells (defined here as Brn3a-positive in the adult retina) and amacrine cells (defined as Brn3a-negative) located in the GCL. This result was surprising, as we were anticipating that MiniP expression may target different subsets of ganglion or amacrine cells, with no expression in the other cell types. However, we observed overlapping expression in both cell types, and our data implies that transcriptional regulation between amacrine and RGC cells may not be distinctly separate. Indeed, in support of this idea are the results of Kunzevitzky et al.,[44] who evaluated the transcriptional profiles using microarrays of amacrine and retinal ganglion cells at E20, P5, and P11 in the rat. They found ~70% overlap in expression between the two cell types at all time points (see Figure 5 in the study of Kunzevitsky et al.). Both amacrine and ganglion cells are born between E14 and E18. Evidence for potential overlapping transcriptional regulation between these two cell types also comes from work that showed VC1.1+ progenitors can generate both VC1.1+ and VC1.1− progenitors.[45,46] The former produces amacrine and horizontal progeny, while the latter produces retinal ganglion cells and cones.[46] Two of these three Ple MiniPs were trimmed prior to use in AAV. In the future, the process of making Pleiades MiniPs even smaller will be informed by data generated by novel technologies. Such data includes that of the ENCODE project,[47] FANTOM projects,[48-50] and VISTA Enhancer project.[5] Specifically, the ENCODE transcription factor ChIP-seq and DNase I hypersensitivity data, along with FANTOM cap analysis of gene expression,[50] and p300 enhancer protein binding[51] data will allow for anchoring of, and supportive evidence for function of, RRs. In addition, the development of ChIP-exo,[52] a refinement to ChIP-seq, brings single-nucleotide resolution and in-depth analysis of regulatory motifs using massively parallel reporter assays.[53,54] In conclusion, this work confirms that MiniP design is a tractable high-throughput option for tissues and shows that it can identify novel promoters that translate across transgene delivery systems. We attribute the success of our design strategy to: (i) the use of evolutionary conserved sequence and human-only DNA limiting transgene inactivation; (ii) paring of putative regulatory elements along with their respective core basal promoters instead of unrelated or viral core promoters, the importance of which is summarized in the study by Haeussler and Joly[20]; and (iii) characterization of reproducible, single-copy knock-ins in a consistent neutral genomic environment. These MiniPs immediately enable greater specificity in gene therapy vectors, and the resource will greatly benefit both basic research and future human gene therapy.

Materials and Methods

Ethics statement

All procedures involving animals were in accordance with the Canadian Council on Animal Care and University of British Columbia Animal Care Committee (protocols #A05-1258, A05-1748, A09-0980, and A09-0981). All virus experiments were approved by the University of Florida Institutional Animal Care and Use Committee and were conducted in accordance with the Association for Research in Vision and Opthalmology Statement for the Use of Animals in Ophthalmic and Vision Research and with National Institutes of Health regulations.

Generation of knock-in mice

Candidate genes for MiniP design were identified in previous datasets[55] as previously described.[9] MiniPs were designed and single-copy knock-in mice at the Hprt locus on the X-chromosome were generated as previously described.[9,56] Animals were housed at the Centre for Molecular Medicine and Therapeutics, in a barrier pathogen-free facility. Animals were kept on a 14-hour on light cycle and provided food and water ad libitum. MiniP strains were backcrossed to either C57BL/6J (JAX Stock#000664) or 129S1/SvImJ (JAX Stock#002448). Animals used for analysis were at least two male germline mice unless otherwise noted with the chimera indicator (CH)—chimeras were analyzed for MiniPs that never entered the germline despite extensive breeding of chimeras and often multiple rounds of microinjection of different embryonic stem cell clones. Each strain and embryonic stem cell line was deposited at The Jackson Laboratory as indicated in the Supplementary Table S2.

PCR analysis of genomic DNA

Vector NTI software (Invitrogen, Carlsbad, CA) was used to design PCR assays for the different constructs. MiniP-specific PCR genotyping assays are available on the http://www.pleiades.org Web site or by contacting the corresponding author.

Immunohistochemistry of adult brains

One millimeter sectioned tissues

The mice were anesthetized with avertin and perfused transcardially with 4% paraformaldehyde (PFA) and postfixed 2–4 hours. The brains were then removed from 4% PFA and immediately sectioned. To section, brains were placed ventral side up into the adult mouse coronal or sagittal Rodent Brain Matrix (ASI Instruments, Warren, MI). Brains were sectioned lateral through medial to lateral (for sagittal) or rostral to caudal (for coronal) using single-edge razor blades (Electron Microscopy Sciences, Hatfield, PA). All slices from one brain were placed into 1 well of a 12-well plate containing 1× phosphate-buffered saline (PBS) (Invitrogen) or 0.1% sodium azide in 1× PBS until staining. Expression of β-gal (lacZ) was detected with X-gal (5-bromo-4-chloro-3-indolyl-β-d-galactopyranoside) staining as previously described.[9] Briefly, the staining was performed with 3–5 ml of X-gal staining solution (25 mg/ml X-gal, 1 mol/l MgCl2, 50 mmol/l potassium ferri-cyanide, 50 mmol/l potassium ferro-cyanide, and 1× PBS to volume) per well in a 24-well plate. The plate was wrapped in aluminum foil and incubated at 37 °C for 10–16 hours with gentle agitation on a rotator board. Subsequently, the sections were transferred into PBS, examined under a dissecting microscope, and photographed using a CoolSNAP-Procf color camera (Media Cybernetics, Rockville, MD) mounted on a Leica MZ12.5 stereomicroscope (Leica Microsystems, Buffalo Grove, IL) and Image-Pro Express v.4.5.1.3 software (Media Cybernetics).

Cryosectioned tissues

The mice were anesthetized with avertin and perfused transcardially with 4% PFA in 0.1 mol/l phosphate buffer (pH 7.4) for 15 minutes. Brains were dissected and postfixed in the same solution for 2 hours and transferred into 25% sucrose-PBS overnight. Each brain was sectioned using a cryostat, and 20-μm sagittal sections were collected. X-gal staining and imaging was performed as described for 1-mm sections. Co-staining of X-gal and NeuN or Gfap was performed sequentially. Primary antibodies used were mouse anti-GFAP (1:1,000; Millipore, Billerica, MA) and mouse anti-NeuN (1:500; Chemicon, Billerica, MA). Bound primary antibodies were detected using the VectaStain Elite ABC Kit (Vector Labs, Burlingame, CA; PK-6100), as per manufacturer’s instructions. Slides were gradually dehydrated in ethanol and xylene series and mounted with coverslips. Bright-field images were taken with an Olympus (Center Valley, PA) BX61 motorized microscope using the DP Controller software, Olympus.

Other tissues

Spinal cord, heart, liver, lung, and ear notch were collected from perfused mice and stained as described above.

Embryo harvesting and staining

For timed pregnancies, female heterozygous MiniP mice were crowded for 7–10 days and set-up in the afternoon with C57BL/6J or 129S1/SvImJ studs. Subsequently, females were checked for copulation plugs for 5 days in the morning. The day of plug was considered E0.5. At E12.5, embryos were harvested, yolk sac tissue taken for PCR analysis, and embryos fixed in 4% PFA/PBS for 2–6 hours. Embryos were stained in X-gal solution for 18 hours, postfixed overnight at 4 °C, and stored in 70% ethanol. Pictures were taken of uncleared embryos. The embryos were then partially cleared based on the method of Schatz et al.[57] In brief, embryos were incubated in step 1 solution (20% glycerol with 0.8% KOH) for 3 days at 37 °C with swirling once daily. Further clearing was done in step 2 solution (50% glycerol with 0.5% KOH) for 3 days at 37 °C. Postfixing of fragile cleared embryos was done in 4% PFA for 2 hours at room temperature. Embryos were stored until imaging in 70% ethanol, after which they were equilibrated to 100% glycerol for imaging. Imaging was performed as per 1-mm brain sections (above). All embryo expression patterns were seen in at least two animals.

Immunohistochemistry on eyes

LacZ staining was detected either directly from MiniPs driving β-gal, or by driving expression of EGFP/cre resulting in recombination of the Gt(ROSA)26Sor locus. After whole body perfusion of mice, eyes were postfixed in 4% PFA. LacZ staining was performed using the X-gal substrate for 18 hours at 37 °C. Eyes were then postfixed in 4% PFA for 2 hours prior to cryoprotection in 25% sucrose-PBS overnight. Eyes were embedded in optimal cutting temperature (OCT), along with positive and negative controls, and cryosectioned at 12 μm and mounted on SuperFrost Plus slides, Fisher Scientific. Eye sections were counterstained with neutral red for 45 seconds to 1 minute. Eye expression was confirmed in at least two animals unless otherwise noted with the chimera indicator (CH). Bright-field images were taken using Olympus BX61 motorized microscope using the DP Controller software.

Immunofluorescence on eyes

Eyes were harvested and fixed in 4% PFA for 2 hours prior to cryoprotection in 25% sucrose-PBS overnight. Eyes were embedded in OCT, and cryosectioned at 12 μm and mounted on SuperFrost Plus slides prior to immunofluorescence. Standard immunofluorescence procedures were followed.[9] Primary antibodies include: chicken anti-β-gal (Abcam Cambridge, MA; ab9361), rabbit anti-β-tubulin III (Covance, Princeton, NJ; PRB-435P), rabbit anti-Brn3 (Santa Cruz Biotechnology, Santa Cruz, CA; sc28595), and rabbit anti-Sox2 (Abcam; ab97959). Secondary antibodies used were (all from Invitrogen): chicken anti-rabbit Alexa-594, goat anti-chicken Alexa-488, goat anti-chicken Alexa-594, goat anti-mouse Alexa-488, goat anti-mouse Alexa-594, goat anti-rabbit Alexa-488, and goat anti-rabbit Alexa-594. Nuclei were stained using 1 μg/ml of Hoechst 33342 (Sigma, St Louis, MO; 14533). Slides were mounted with Prolong Gold (Invitrogen) antifade reagent and imaged using an Olympus BX61 microscope and the InVivo software (Olympus), or alternatively using confocal laser scanning microscopy with a Leica SP5 II microscope (Leica Microsystems) and LAF software (Leica Microsystems).

Adeno-associated virus

Recombinant AAV vector plasmids containing a 2,202 bp Ple67 (FEV RRs), a 3,310 bp Ple53 (DCX RRs) or 3,312 bp Ple25 (CCKBR RRs) MiniPromoter construct, an SV40 splice donor–acceptor site followed by the cDNA for GFP were constructed. To accommodate the packaging limit of recombinant AAV, Ple25 and Ple53 were shortened relative to the knock-in MiniP allele as follows: MiniPs were amplified by PCR, and the respective restriction sites were added via the oligonucleotides. The following oligonucleotides were used: Ple25 forward KpnI: GTCGGGTACCAGAATTTCAGGAGCCCATATG, Ple25 reverse XhoI: GTCGCTCGAGTGCCTTACTCAGCTCGACC, Ple53 forward KpnI: GTCGGGTACCTTATCCCATCCATTGTTGTATC, Ple53 reverse XhoI: GTCGCTCGAGGGTGGAACCTCAGAGACCTG, Ple67 forward EcoRI: GCTGGAATTCGGAGAAGATCAAGACTCAGG, and Ple67 reverse XhoI: GCTGCTCGAGGAGGCAGGCTTTGCGCTC, (restriction sites underlined). The plasmids pEMS1496, pEMS1524, and pEMS1538 were used as templates. The shorter Ple25 promoter construct is referred to as Ple25(Fmd1), and the shorter Ple53 promoter construct as Ple53(Fmd1). Identity for all promoters was confirmed by sequencing the recombinant AAV vector plasmids before packaging into AAV2-(quad Y-F), in which the VP3 capsid protein contains four surface exposed tyrosine to phenylalanine mutations (Y272F, Y444F, Y500F, and Y730F).[24] Vectors were purified and titered according to previously published methods.[58] At least 5 × 109 vector genomes (vg) per eye were injected intravitreally in 4-week-old C57BL/6 mice. GFP expression was analyzed 4 weeks postinjection. Mice were euthanized, and the eyes were fixed in 4% PFA for 2 hours at room temperature. Eyes were then treated with a stepwise sucrose gradient (from 4 to 20% sucrose-PBS) and embedded and frozen in Tissue-Tek OCT compound (Sakura Finetek, Torrance, CA). Ten-micron-thick sections of each eye were then washed in PBS for 15 minutes, mounted in VECTASHIELD mounting medium (Vector Labs) with 4′,6-diamidino-2-phenylindole, and analyzed for GFP expression by confocal microscopy.

Image processing

Images were processed using ImageJ (http://rsbweb.nih.gov/ij/), Photoshop, and Illustrator (Adobe, San Jose, CA). Brightness, contrast, and scaling adjustments were performed. For images where brains were not properly aligned (Figure 2), they were traced with magnetic lasso and repositioned on a black background. The auto-tone feature in Photoshop was used to correct for color balance and exposure settings for images of eyes (Figure 6).

Data and resource distribution

Data on MiniPs can also be found on Pleiades Promoter Project website (http://www.pleiades.org/). All plasmids and sequence information are available from Addgene (http://www.addgene.org/). Embryonic stem cells and mice have been deposited at, and are available from, the Mutant Mouse Regional Resource Centers (http://www.mmrrc.org/) and The Jackson Laboratory (http://www.jax.org/).
  57 in total

Review 1.  Neuronal cell types.

Authors:  Richard H Masland
Journal:  Curr Biol       Date:  2004-07-13       Impact factor: 10.834

Review 2.  When needles look like hay: how to find tissue-specific enhancers in model organism genomes.

Authors:  Maximilian Haeussler; Jean-Stéphane Joly
Journal:  Dev Biol       Date:  2010-12-03       Impact factor: 3.582

Review 3.  Some certainty for the "zone of uncertainty"? Exploring the function of the zona incerta.

Authors:  J Mitrofanis
Journal:  Neuroscience       Date:  2005       Impact factor: 3.590

4.  Efficient inducible Pan-neuronal cre-mediated recombination in SLICK-H transgenic mice.

Authors:  Victoria Heimer-McGinn; Paul Young
Journal:  Genesis       Date:  2011-08-23       Impact factor: 2.487

5.  Production of a mouse monoclonal antibody reactive with a human nuclear antigen associated with cell proliferation.

Authors:  J Gerdes; U Schwab; H Lemke; H Stein
Journal:  Int J Cancer       Date:  1983-01-15       Impact factor: 7.396

6.  Comprehensive genome-wide protein-DNA interactions detected at single-nucleotide resolution.

Authors:  Ho Sung Rhee; B Franklin Pugh
Journal:  Cell       Date:  2011-12-09       Impact factor: 41.582

7.  Targeted transgene expression in neuronal precursors: watching young neurons in the old brain.

Authors:  Sebastien Couillard-Despres; Beate Winner; Claudia Karl; Gudrun Lindemann; Peter Schmid; Robert Aigner; Joern Laemke; Ulrich Bogdahn; Juergen Winkler; Josef Bischofberger; Ludwig Aigner
Journal:  Eur J Neurosci       Date:  2006-09       Impact factor: 3.386

8.  Requirement for serum response factor for skeletal muscle growth and maturation revealed by tissue-specific gene deletion in mice.

Authors:  Shijie Li; Michael P Czubryt; John McAnally; Rhonda Bassel-Duby; James A Richardson; Franziska F Wiebel; Alfred Nordheim; Eric N Olson
Journal:  Proc Natl Acad Sci U S A       Date:  2005-01-12       Impact factor: 11.205

9.  The transcriptional landscape of the mammalian genome.

Authors:  P Carninci; T Kasukawa; S Katayama; J Gough; M C Frith; N Maeda; R Oyama; T Ravasi; B Lenhard; C Wells; R Kodzius; K Shimokawa; V B Bajic; S E Brenner; S Batalov; A R R Forrest; M Zavolan; M J Davis; L G Wilming; V Aidinis; J E Allen; A Ambesi-Impiombato; R Apweiler; R N Aturaliya; T L Bailey; M Bansal; L Baxter; K W Beisel; T Bersano; H Bono; A M Chalk; K P Chiu; V Choudhary; A Christoffels; D R Clutterbuck; M L Crowe; E Dalla; B P Dalrymple; B de Bono; G Della Gatta; D di Bernardo; T Down; P Engstrom; M Fagiolini; G Faulkner; C F Fletcher; T Fukushima; M Furuno; S Futaki; M Gariboldi; P Georgii-Hemming; T R Gingeras; T Gojobori; R E Green; S Gustincich; M Harbers; Y Hayashi; T K Hensch; N Hirokawa; D Hill; L Huminiecki; M Iacono; K Ikeo; A Iwama; T Ishikawa; M Jakt; A Kanapin; M Katoh; Y Kawasawa; J Kelso; H Kitamura; H Kitano; G Kollias; S P T Krishnan; A Kruger; S K Kummerfeld; I V Kurochkin; L F Lareau; D Lazarevic; L Lipovich; J Liu; S Liuni; S McWilliam; M Madan Babu; M Madera; L Marchionni; H Matsuda; S Matsuzawa; H Miki; F Mignone; S Miyake; K Morris; S Mottagui-Tabar; N Mulder; N Nakano; H Nakauchi; P Ng; R Nilsson; S Nishiguchi; S Nishikawa; F Nori; O Ohara; Y Okazaki; V Orlando; K C Pang; W J Pavan; G Pavesi; G Pesole; N Petrovsky; S Piazza; J Reed; J F Reid; B Z Ring; M Ringwald; B Rost; Y Ruan; S L Salzberg; A Sandelin; C Schneider; C Schönbach; K Sekiguchi; C A M Semple; S Seno; L Sessa; Y Sheng; Y Shibata; H Shimada; K Shimada; D Silva; B Sinclair; S Sperling; E Stupka; K Sugiura; R Sultana; Y Takenaka; K Taki; K Tammoja; S L Tan; S Tang; M S Taylor; J Tegner; S A Teichmann; H R Ueda; E van Nimwegen; R Verardo; C L Wei; K Yagi; H Yamanishi; E Zabarovsky; S Zhu; A Zimmer; W Hide; C Bult; S M Grimmond; R D Teasdale; E T Liu; V Brusic; J Quackenbush; C Wahlestedt; J S Mattick; D A Hume; C Kai; D Sasaki; Y Tomaru; S Fukuda; M Kanamori-Katayama; M Suzuki; J Aoki; T Arakawa; J Iida; K Imamura; M Itoh; T Kato; H Kawaji; N Kawagashira; T Kawashima; M Kojima; S Kondo; H Konno; K Nakano; N Ninomiya; T Nishio; M Okada; C Plessy; K Shibata; T Shiraki; S Suzuki; M Tagami; K Waki; A Watahiki; Y Okamura-Oho; H Suzuki; J Kawai; Y Hayashizaki
Journal:  Science       Date:  2005-09-02       Impact factor: 47.728

10.  The evolving doublecortin (DCX) superfamily.

Authors:  Orly Reiner; Frédéric M Coquelle; Bastian Peter; Talia Levy; Anna Kaplan; Tamar Sapir; Irit Orr; Naama Barkai; Gregor Eichele; Sven Bergmann
Journal:  BMC Genomics       Date:  2006-07-26       Impact factor: 3.969

View more
  23 in total

1.  Intravitreal delivery of a novel AAV vector targets ON bipolar cells and restores visual function in a mouse model of complete congenital stationary night blindness.

Authors:  Miranda L Scalabrino; Sanford L Boye; Kathryn M H Fransen; Jennifer M Noel; Frank M Dyka; Seok Hong Min; Qing Ruan; Charles N De Leeuw; Elizabeth M Simpson; Ronald G Gregg; Maureen A McCall; Neal S Peachey; Shannon E Boye
Journal:  Hum Mol Genet       Date:  2015-08-26       Impact factor: 6.150

Review 2.  Gene therapy for the nervous system: challenges and new strategies.

Authors:  Casey A Maguire; Servio H Ramirez; Steven F Merkel; Miguel Sena-Esteves; Xandra O Breakefield
Journal:  Neurotherapeutics       Date:  2014-10       Impact factor: 7.620

Review 3.  Strategies for targeting primate neural circuits with viral vectors.

Authors:  Yasmine El-Shamayleh; Amy M Ni; Gregory D Horwitz
Journal:  J Neurophysiol       Date:  2016-04-06       Impact factor: 2.714

4.  Prdm1 overexpression causes a photoreceptor fate-shift in nascent, but not mature, bipolar cells.

Authors:  Noah B Goodson; Ko U Park; Jason S Silver; Vince A Chiodo; William W Hauswirth; Joseph A Brzezinski
Journal:  Dev Biol       Date:  2020-06-17       Impact factor: 3.582

Review 5.  Advances in optogenetic and chemogenetic methods to study brain circuits in non-human primates.

Authors:  Adriana Galvan; Michael J Caiola; Daniel L Albaugh
Journal:  J Neural Transm (Vienna)       Date:  2017-02-25       Impact factor: 3.575

6.  AAV capsid variants with brain-wide transgene expression and decreased liver targeting after intravenous delivery in mouse and marmoset.

Authors:  David Goertsen; Nicholas C Flytzanis; Nick Goeden; Miguel R Chuapoco; Alexander Cummins; Yijing Chen; Yingying Fan; Qiangge Zhang; Jitendra Sharma; Yangyang Duan; Liping Wang; Guoping Feng; Yu Chen; Nancy Y Ip; James Pickel; Viviana Gradinaru
Journal:  Nat Neurosci       Date:  2021-12-09       Impact factor: 24.884

7.  Targeted Expression of Retinoschisin by Retinal Bipolar Cells in XLRS Promotes Resolution of Retinoschisis Cysts Sans RS1 From Photoreceptors.

Authors:  Camasamudram Vijayasarathy; Yong Zeng; Dario Marangoni; Lijin Dong; Zhuo-Hua Pan; Elizabeth M Simpson; Robert N Fariss; Paul A Sieving
Journal:  Invest Ophthalmol Vis Sci       Date:  2022-10-03       Impact factor: 4.925

8.  Intracerebroventricular Administration of AAV9-PHP.B SYN1-EmGFP Induces Widespread Transgene Expression in the Mouse and Monkey Central Nervous System.

Authors:  Adriana Galvan; Terri L Petkau; Austin M Hill; Andrea J Korecki; Ge Lu; Diane Choi; Kazi Rahman; Elizabeth M Simpson; Blair R Leavitt; Yoland Smith
Journal:  Hum Gene Ther       Date:  2021-06-01       Impact factor: 4.793

9.  Intravenous AAV9 efficiently transduces myenteric neurons in neonate and juvenile mice.

Authors:  Sara E Gombash; Christopher J Cowley; Julie A Fitzgerald; Jodie C E Hall; Christian Mueller; Fedias L Christofi; Kevin D Foust
Journal:  Front Mol Neurosci       Date:  2014-10-15       Impact factor: 5.639

10.  AAV-mediated transduction and targeting of retinal bipolar cells with improved mGluR6 promoters in rodents and primates.

Authors:  Q Lu; T H Ganjawala; E Ivanova; J G Cheng; D Troilo; Z-H Pan
Journal:  Gene Ther       Date:  2016-04-26       Impact factor: 5.250

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.