Adeno-associated virus site-specifically integrates into a muscle-specific DNA region
Abstract
The nonpathogenic human virus adeno-associated virus type 2 (AAV) has evolved the potentially unique strategy to establish latency by site-specifically integrating its genome into human chromosome 19 (19q13.3-qter) at a locus designated AAVS1. This nonhomologous, site-specific recombination of viral DNA with the human genome provides a basis for developing targeted gene therapy vectors. To assess whether the region surrounding AAVS1 might have contributed to the selection of the specific integration site, we have investigated this locus. Here, we show that AAVS1 is closely linked to the slow skeletal troponin T gene, TNNT1, which has been mapped previously to 19q13.4. In support of this idea, we demonstrate that site-specific AAV DNA integration can result in the formation of TNNT1-AAV junctions. The question now arises whether muscle represents a natural target tissue for latent AAV infection. This possibility is of additional interest in view of recent observations that muscle tissue is particularly well suited for AAV-mediated gene transfer. The question also occurs whether latent infection by AAV can lead to phenotypic changes of the multinucleated muscle fiber cells.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
The human parvovirus adeno-associated virus type 2 (AAV) has adopted the strategy of alternate life styles to ensure a successful relationship with its host (1). This nonpathogenic virus productively replicates only in the presence of helper virus co- or superinfection (2). In the absence of such helper conditions, AAV has evolved a strategy potentially unique among eukaryotic viruses. On infection of the non-coinfected host cell, AAV establishes latency by integrating its genome site-specifically into a region on the q arm of human chromosome 19 (19q13.3-qter) (3, 4), designated AAVS1.
The 4.7-kb linear single-stranded genome contains two ORFs and is flanked by inverted terminal repeats (ITRs) (5). The right ORF encodes three capsid proteins, and the left ORF encodes the four nonstructural proteins (Rep proteins) that are involved in regulating all aspects of the viral life style (2).
The 145-nt ITRs are the only viral sequences required in cis for DNA replication, packaging of the viral genome into the capsid, and site-specific DNA integration. Within the ITRs, a Rep binding site (RBS) allows for specific recruitment of the large Rep proteins (i.e., Rep68 and Rep78) to the origin of replication (6–8). A Rep-specific endonuclease site (terminal resolution site, TRS) is separated from the RBS by a 13-nt spacer (9–11). Together, RBS and TRS can act as a minimal origin for Rep-mediated DNA replication (12).
To dissect the requirements for site specificity in targeted AAV DNA integration, an 8.2-kb fragment of AAVS1 has been isolated from embryonic lung fibroblast (WI38) DNA, and a 4-kb sequence was determined (13). It was shown by using an episomal model system that site specificity is determined by cellular sequences (14) and that a 33-nt sequence is necessary and sufficient for this targeted nonhomologous recombination event to occur (15). This 33-nt chromosomal sequence is similar to the minimal viral origin of DNA replication, consisting of an RBS and a TRS (16), leading to the concept that Rep-mediated DNA replication is involved in the integration mechanism. Complementing this idea was the observation that the viral Rep proteins are required for site-specific integration (17, 18). Biochemical assays further showed that the Rep proteins can specifically interact with the viral and cellular RBS and TRS motifs (8, 19, 20) to mediate replication and potentially targeted integration of AAV into AAVS1 (16, 21). Through isolation and analysis of several viral–cellular junctions, it became apparent that, although all of the junctions contained AAVS1 sequences, the immediate transitions from viral to cellular sequences were scattered over a range of approximately 1,000 nucleotides downstream of the TRS-RBS motif within AAVS1 (4, 22, 23). These observations are in agreement with the hypothesis that limited cellular DNA replication is involved in the initial steps of the mechanism underlying AAV site-specific integration. This model predicts that Rep-mediated replication is using the RBS and TRS motifs present on AAVS1 as a minimal origin that is potentially unique within the human genome (16).
Several questions remain with respect to the requirements for the chromosomal target sequence. It has not been determined whether active transcription in the vicinity of AAVS1 promotes integration by providing access to the target sequence for the proteins involved in the recombination mechanism. A second question is whether AAV integration could lead to the disruption of an active cellular gene. In most of the previous studies, DNA fragments containing only one viral–cellular junction had been isolated, leaving open the position and characteristics of additional junctions.
In this report, we present our analysis of the region surrounding AAVS1. We show that the target sequence for AAV site-specific integration is closely linked to muscle-specific genes TNNT1 and TNNI3. In support of this finding, we demonstrate that site-specific AAV DNA integration can result in the formation of TNNT1-AAV junctions.
Materials and Methods
Cell Lines.
Libraries.
Similar strategies were used to construct the HA16C5 (latently infected) and a human lymphocyte (uninfected) genomic library. The genomic DNA was partially digested with Sau3A and cloned into the BamHI site of the λDASHII vector (Stratagene). The HA16C5 and lymphocyte libraries were screened with an AAV and an AAVS1 probe, respectively. From each library, several phages positive for either AAVS1 (for the uninfected lymphocyte library) or AAV hybridization (for the HA16C5 library) were isolated, digested with NotI, and subcloned into the pBluescript vector. The plasmids were grown at 30°C in SURE bacteria (Stratagene).
Plasmid pHA16–3 was isolated from the HA16C5 library and pS1 from the lymphocyte library. The pS1 plasmid contains a 14-kb insert starting at position 2083 in AAVS1.
Plasmids, Probes, and Cosmids.
pBluescript was used as a cloning vector for the following subcloning steps (exceptions are indicated).
The AAVS1pRVK plasmid (nucleotides 1–3536) was obtained by EcoRI/KpnI digestion of the pRIA plasmid (13). The pS1 plasmid was digested with PstI to generate the AAVS1pPs plasmid (nucleotides 1438 to −2953 relative to AAVS1). AAVS1pPs was digested with DraI/PvuII to generate plasmid AAVS1pPD (nucleotides −1723 to −1970 relative to AAVS1).
The expressed sequence tag (EST) clone AA176454, which is homologous to the TNNT1 mRNA M1 clone (27), was purchased from the American Type Culture Collection (catalog no. 769879) and used to generate the TNNT1 cDNA probes. This clone was digested with XhoI-BamHI to generate the 5′TNNT1 plasmid (nucleotides 56–663) (27) and the 3′TNNT1 plasmid (nucleotides 663–896) (27). The TNNT1ex11–13 insert (nucleotides 562–754) (27) was generated by PCR on the AA176454 clone with the primers RML139 (nucleotides 562–580) (27) and RML140 (nucleotides 733–754) (27), and the PCR product was subsequently cloned into the PCR2.1 vector (Invitrogen). The cosmid c23932 was digested with StuI-BglII to generate the TNNI3ex7 plasmid (nucleotides 4653–4877). All of the inserts of these clones were sequenced.
The cosmids used for this study were generously provided by the Human Genome Center at Lawrence Livermore National Laboratory (LLNL).
The TNNT1- and AAVS1-containing EcoRI fragments of cosmid 31925 were subcloned into pBluescript, and terminal sequences were determined for each fragment by using T3 and T7 primers.
Southern Analysis.
Human genomic DNA was extracted from peripheral blood by using the QIAmp blood DNA kit (Qiagen, Chatsworth, CA). Genomic DNA from cell lines was isolated by standard procedures. Genomic DNA (10 μg) was digested, separated by gel electrophoresis on a 0.8% agarose gel, transferred onto a nylon membrane, and probed with a radiolabeled [α-32P]dCTP probe. Probes are indicated in the figure legends. The membranes were hybridized overnight at 65°C in a 0.75× nylon buffer (1× = 14% SDS/80 mM Na2HPO4/14 mM Na2H2EDTA⋅2H2O) and washed for 20 min each at 65°C in 0.5×, 0.1×, and 0.05× nylon buffer.
GenBank Accession Numbers.
AAV2, AF043303; AAVS1, S51329; cosmid 31855, AC005782; TNNT1 M1 mRNA clone, M19308; TNNT1 gene exons 1–11, AJ011712; TNNT1 gene exons 12–14, AJ011713; TNNI3 gene, X90780.1.
Results and Discussion
To determine whether AAV integration can lead to the disruption of a functional gene closely linked to AAVS1, a genomic library of HeLa cells latently infected with AAV (HA16C5) was screened with an AAV probe, and junction sequences were analyzed. All of the plasmids positive for an AAV-specific probe contained AAVS1 junction sequences (data not shown). One plasmid contained a fragment with a second junction between AAV and cellular DNA in addition to the AAV-AAVS1 junction. Sequence analysis revealed that this cellular fragment was derived from the TNNT1 gene, suggesting the possibility that AAVS1 and TNNT1 are linked. This possibility is in agreement with the previously mapped position of TNNT1 at 19q13.4 (28), which is included within the region reported for AAVS1 (3, 4).
Characterization of the Region Surrounding AAVS1.
Although the presence of AAV-AAVS1 and AAV-TNNT1 junctions within the same fragment was indicative for colocalization of AAVS1 with TNNT1, the possibilities could not be excluded that both sequences were brought to close vicinity either by cellular rearrangements inherent to the transformed HeLa cell line or by rearrangements induced by AAV integration. To test the hypothesis that AAVS1 and TNNT1 are linked, several AAVS1 probes and TNNT1 cDNA probes (Fig. 1) were used for Southern analyses of nine cosmids that had been preselected for the presence of AAVS1. These cosmids were derived from a library from flow-sorted chromosomes of a human–hamster hybrid, which carries chromosome 19 as its only human chromosome (constructed by the Human Genome Center at LLNL). In addition, the plasmid pS1, which was derived from a different, independently generated genomic library was used in these analyses. The cosmids and pS1 DNA were digested with EcoRI, and the resulting fragments were hybridized to the probes indicated in Fig. 1B. The resulting hybridization patterns were compared with those derived from genomic DNA isolated from both human peripheral blood and HEK293 cells. The analysis of genomic DNA was added to exclude rearrangements of cosmid sequences. Cosmid 31855 (insert size 35.2 kb) (Fig. 1A) extended from AAVS1 into the centromeric region. During the course of our studies, the complete sequence of cosmid 31855 was determined by the LLNL. Sequence comparison of AAVS1 with cosmid 31855 clearly linked AAVS1 to the EcoRI map graphically displayed in Fig. 1A. Note that, to date, within the 35.2 kb contained within this cosmid, no known cellular gene is found. Of eight cosmids extending into the telomeric direction of AAVS1, two cosmids hybridized to TNNT1- and AAVS1-specific probes (cosmids 31925, insert size: 37.5 kb; and 28847, insert size: 37 kb). However, cosmid 28847 had undergone rearrangements as determined by restriction patterns inconsistent with the respective EcoRI map (Fig. 1A, BC788605) provided by the LLNL and was therefore omitted (data not shown). EcoRI fragments of the inserts from cosmid 31925 and plasmid pS1 were subcloned, and terminal sequences were determined. These data allowed positioning of the inserts with respect to the previously published AAVS1 sequence (data not shown) (13). Consistent hybridization patterns were observed with cosmid 31925 and genomic DNA with the AAVS1 and TNNT1 probes (Fig. 1C). Note that the hybridization pattern of cosmid 31925 and plasmid pS1 with the AAVS1 probe is different from genomic DNA. As seen in Fig. 1A, cosmid 31925 contains only 4.5 kb of the 8-kb EcoRI fragment of AAVS1 in addition to the full TNNT1 gene and pS1 contains ca. 2 kb of the AAVS1 EcoRI fragment. In summary, our data support a close linkage between AAVS1 and TNNT1.
Figure 1
The LLNL EcoRI map (Fig. 1A; BC788605) was used to organize the fragments derived from our restriction and Southern analyses of the cosmids, pS1, and genomic DNA, resulting in an estimated distance between AAVS1 and TNNT1 of 14.7 kb. To confirm this distance between TNNT1 and AAVS1, cosmid 31925, plasmid pS1, and genomic DNA were digested with HindIII and hybridized to TNNT1 and to AAVS1 probes (Fig. 1D). As evident from the hybridization patterns shown in Fig. 1D, the distances deduced from the HindIII digests are consistent with those from the EcoRI digests, and we therefore conclude that AAVS1 and TNNT1 are separated by 14.7 kb.
TNNT1 Rearrangements in Latently Infected Cells.
The proximity of AAVS1 and TNNT1 led us to examine several cell lines for evidence that site-specific integration of AAV DNA into AAVS1 could lead to the disruption of the TNNT1 gene. Southern blots were performed on genomic DNA isolated from several latently infected cell lines, the Detroit 6 (D6) line 7374 (26), the KB lines M19 (24), M21 (24), and the HeLa line HA16C5 (25). These cell lines have previously been characterized as containing at least one copy of AAV DNA integrated into AAVS1 (refs. 23 and 29, and unpublished observations). As shown in Fig. 2, the 8-kb and 6-kb bands were detected after EcoRI and HindIII digestion, respectively, followed by hybridization to the pRVK probe. These bands correspond to the undisrupted AAVS1 locus, apparently present in all of the latently infected cells as well as in the control DNA isolated from peripheral blood (4, 13, 23, 29) (Fig. 2 A and B). As previously shown, only the infected cells showed supernumerous bands with an AAVS1 probe (4, 23, 29); some of these bands cohybridized with an AAV probe, indicating a possible AAV-AAVS1 junction fragment (Fig. 2 A and B). Hybridization of the same DNA to a 3′TNNT1 probe resulted in the detection of the undisrupted 7-kb (EcoRI, Fig. 2A) and 23-kb (HindIII, Fig. 2B) fragments (27, 30). After an EcoRI and a HindIII digestion, the M19 and HA16C5 cells showed additional bands, indicating a possible disruption of the TNNT1 gene (Fig. 2 A and B); some of these bands cohybridize with an AAV probe, suggesting linkage between AAV DNA and the TNNT1 gene. Moreover, in the HA16C5 DNA, one band cohybridized with the AAV, AAVS1, and 3′TNNT1 probes (Fig. 2A). These data indicate that integration of AAV into the AAVS1 locus can lead to rearrangements of the TNNT1 gene.
Figure 2
Recently, Barton et al. (30) found that the TNNT1 gene is closely linked to the cardiac troponin I gene (TNNI3), which had previously been mapped to 19q13.4 (31). The authors concluded that the distance between TNNT1 and TNNI3 is only 2.6 kb. Based on these findings, we performed Southern analysis on genomic DNA of the four latently infected cell lines. However, when the TNNI3 probe indicated in Materials and Methods was used, no evidence for disruption was found (data not shown).
Sequence Analysis of an AAV-TNNT1 Junction Fragment.
The junction sequences isolated from the HA16 library (and cloned to give plasmid pHA16-3) provide direct support for the hypothesis that AAV integration was the cause for the observed TNNT1 rearrangements (Fig. 3 A and B). As indicated above, this fragment contains two junctions, one between the right AAV ITR and AAVS1 (near the RBS) and a second between another right AAV ITR and the last intron (intron 13) of the TNNT1 gene. This sequence analysis demonstrates that integration of AAV into AAVS1 can indeed disrupt the TNNT1 gene. The junction fragment showed characteristics similar to those previously described for AAV-AAVS1 junctions (4, 13, 32): the ITR was partially deleted, the immediate junction consisted of nucleotides common to both donor and acceptor DNA, and the AAV RBS was found at the viral–cellular junctions.
Figure 3
To rule out the possibility that the junction analyzed (and presented in Fig. 3) might be the result of recombination in Escherichia coli, PCR analysis was performed on genomic DNA from HA16 cells. PCR using a TNNT1 (derived from exon 14) and an AAV-specific primer (derived from the cap gene) resulted in one band of the expected size, the identity of which was confirmed by restriction analysis (data not shown).
In summary, we conclude that AAVS1 is located 14.7 kb and 34.1 kb centromeric to TNNT1 and TNNI3 (Fig. 4), respectively, at position 19q13.4. Southern analysis of latently infected cells and sequence analysis of a viral-cellular junction showed that the site-specific integration of AAV DNA at AAVS1 led to rearrangements of the TNNT1 gene in two of the four cell lines studied.
Figure 4
Troponin and tropomyosin constitute the Ca2+-sensitive complex regulating the contraction of muscle fibers. The troponin complex consists of three subunits, troponin C (Ca2+ binding), troponin I (inhibitory subunit), and troponin T (tropomyosin binding subunit). There are three genes encoding troponin I isoforms, TNNI1–3, and three genes encoding troponin T isoforms, TNNT1–3. The three respective isoforms represent the slow (TNNT1, TNNI1), fast (TNNT2, TNNI2), and cardiac troponins (TNNT3, TNNI3). The expression of the troponin genes is tightly regulated in a tissue-dependent manner, as well as during development. In adults, TNNT1 is exclusively expressed in skeletal muscle, whereas TNNI3 is exclusively expressed in the heart. It is of note that TNNT1 expression was detected by Northern blot analysis in two latently infected cell lines, Detroit 6 (7374) and HA16 (data not shown). It remains to be determined whether induction of TNNT1 expression occurred as a result of AAV DNA integration or whether the parental, uninfected cells expressed TNNT1 before AAV infection.
The site-specific integration of AAV DNA has attracted considerable attention with respect to its potential for a targeted gene transfer that would minimize the risk for insertional mutagenesis. Our report highlights that the TNNT1 gene can be disrupted as a result of wild-type AAV DNA integration, raising several issues. Studies of AAV integration have predominantly used immortalized or transformed cells that have lost chromosome stability. The extent of AAV-induced rearrangements, including the disruption of TNNT1 observed in our study, may be related to loss of these, as yet unknown, control mechanisms. Furthermore, the troponin genes are telomeric with respect to AAVS1. Most of the previously characterized junctions (4, 23, 29, 32), however, were found centromeric (“downstream”) of the AAVS1 TRS (see also Fig. 3A). A model for the mechanism of AAV integration must incorporate the potential for “upstream” rearrangements.
It is noteworthy that AAV has evolved to integrate close to, and potentially disrupt, the highly muscle-specific gene, TNNT1. Based on the idea that, in the adult, TNNT1 is exclusively expressed in skeletal muscle tissues (33), the question arises whether muscle represents a natural target tissue for latent AAV infection. In support of this possibility, (E. P. Hoffman and coworkers, personal communication) find AAV DNA in 17% of normal human muscle biopsies tested.
It is possible that integration of AAV into or close to TNNT1 might change the level of expression of either the TNNT1 or the closely linked TNNI3 genes. It can now also be asked whether altered mRNA species and potentially proteins are generated in human muscle cells as a result of AAV integration. However, skeletal muscle fiber cells are multinucleated, so the disruption of a subset of TNNT1 alleles might not lead to phenotypic changes in these cells. Hence, AAV site-specific integration into this locus might represent yet another strategy by the virus to optimize the virus-host relationship.
Abbreviations
- AAV
- adeno-associated virus type 2
- ITR
- inverted terminal repeat
- RBS
- Rep binding site
- TRS
- terminal resolution site
- LLNL
- Lawrence Livermore National Laboratory
Notes
Article published online before print: Proc. Natl. Acad. Sci. USA, 10.1073/pnas.080079397.
Article and publication date are at www.pnas.org/cgi/doi/10.1073/pnas.080079397
Acknowledgments
We thank Agnes Begue for her help in constructing the genomic library HA16C5 and for providing the lymphocyte library. Furthermore, we thank Dr. Linda Ashworth for helpful discussions and for providing the cosmids, and Dr. Jörg R. Schlehofer for providing the parental HA16 cells from which HA16C5 was subcloned. We thank Drs. Mary Klotman, Nuria Morral, Nick Muzyzcka, and Peter Ward for critical reading of the manuscript. This work was supported in part by Grant DK55609 (to R.M.L.) from the National Institutes of Health.
References
1
K I Berns, R M Linden BioEssays 17, 237–245 (1995).
2
K I Berns Fields Virology, eds B N Fields, D M Knipe, P M Howley (Lippincott–Raven, Philadelphia) 2, 2173–2197 (1996).
3
R M Kotin, J C Menninger, D C Ward, K I Berns Genomics 10, 831–834 (1991).
4
R J Samulski, X Zhu, X Xiao, J D Brook, D E Housman, N Epstein, L A Hunter EMBO J 10, 3941–3950 (1991).
5
A Srivastava, E W Lusby, K I Berns J Virol 45, 555–564 (1983).
6
D M McCarty, D J Pereira, I Zolotukhin, X Zhou, J H Ryan, N Muzyczka J Virol 68, 4988–4997 (1994).
7
J H Ryan, S Zolotukhin, N Muzyczka J Virol 70, 1542–1553 (1996).
8
J A Chiorini, L Yang, B Safer, R M Kotin J Virol 69, 7334–7338 (1995).
9
D S Im, N Muzyczka Cell 61, 447–457 (1990).
10
D M McCarty, J H Ryan, S Zolotukhin, X Zhou, N Muzyczka J Virol 68, 4998–5006 (1994).
11
J R Brister, N Muzyczka J Virol 73, 9325–9336 (1999).
12
D H Smith, P Ward, R M Linden J Virol 73, 2930–2937 (1999).
13
R M Kotin, R M Linden, K I Berns EMBO J 11, 5071–5078 (1992).
14
C Giraud, E Winocour, K I Berns Proc Natl Acad Sci USA 91, 10039–10043 (1994).
15
R M Linden, E Winocour, K I Berns Proc Natl Acad Sci USA 93, 7966–7972 (1996).
16
R M Linden, P Ward, C Giraud, E Winocour, K I Berns Proc Natl Acad Sci USA 93, 11288–11294 (1996).
17
C Balague, M Kalla, W W Zhang J Virol 71, 3299–3306 (1997).
18
R T Surosky, M Urabe, S G Godwin, S A McQuiston, G J Kurtzman, K Ozawa, G Natsoulis J Virol 71, 7951–7959 (1997).
19
M D Weitzman, S R Kyostio, R M Kotin, R A Owens Proc Natl Acad Sci USA 91, 5808–5812 (1994).
20
J A Chiorini, M D Weitzman, R A Owens, E Urcelay, B Safer, R M Kotin J Virol 68, 797–804 (1994).
21
R M Kotin Hum Gene Ther 5, 793–801 (1994).
22
X Zhu Ph.D. thesis (Univ. of Pittsburgh, Pittsburgh, 1993).
23
R M Kotin, K I Berns Virology 170, 460–467 (1989).
24
C A Laughlin, C B Cardellichio, H C Coon J Virol 60, 515–524 (1986).
25
C Walz, J R Schlehofer J Virol 66, 2990–3002 (1992).
26
K I Berns, T C Pinkerton, G F Thomas, M D Hoggan Virology 68, 556–560 (1975).
27
R Gahlmann, A B Troutt, R P Wade, P Gunning, L Kedes J Biol Chem 262, 16122–16126 (1987).
28
F Samson, P J de Jong, B J Trask, P Koza-Taylor, M C Speer, T Potter, A D Roses, J R Gilbert Genomics 13, 1374–1375 (1992).
29
R M Kotin, M Siniscalco, R J Samulski, X D Zhu, L Hunter, C A Laughlin, S McLaughlin, N Muzyczka, M Rocchi, K I Berns Proc Natl Acad Sci USA 87, 2211–2215 (1990).
30
P J Barton, M E Cullen, P J Townsend, N J Brand, A J Mullen, D A Norman, P K Bhavsar, M H Yacoub Genomics 57, 102–109 (1999).
31
J Mogensen, T A Kruse, A D Borglum Cytogenet Cell Genet 79, 272–273 (1997).
32
C Giraud, E Winocour, K I Berns J Virol 69, 6917–6924 (1995).
33
M A Sabry, G K Dhoot J Muscle Res Cell Motil 12, 262–270 (1991).
34
P K Bhavsar, N J Brand, M H Yacoub, P J R Barton Genomics 35, 11–23 (1996).
Information & Authors
Information
Published in
Classifications
Copyright
Copyright © The National Academy of Sciences.
Submission history
Received: November 16, 1999
Accepted: February 24, 2000
Published online: April 11, 2000
Published in issue: April 25, 2000
Acknowledgments
We thank Agnes Begue for her help in constructing the genomic library HA16C5 and for providing the lymphocyte library. Furthermore, we thank Dr. Linda Ashworth for helpful discussions and for providing the cosmids, and Dr. Jörg R. Schlehofer for providing the parental HA16 cells from which HA16C5 was subcloned. We thank Drs. Mary Klotman, Nuria Morral, Nick Muzyzcka, and Peter Ward for critical reading of the manuscript. This work was supported in part by Grant DK55609 (to R.M.L.) from the National Institutes of Health.
Authors
Metrics & Citations
Metrics
Citation statements
Altmetrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFGet Access
Login options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to access the full text.